--- title: Skyrmion Lattice Model Overview url: https://www.emergentmind.com/topics/skyrmion-lattice-model type: topic --- # Skyrmion Lattice Model Overview A skyrmion lattice model is a theoretical description of a periodic array of topological textures, most commonly a triangular skyrmion crystal in a chiral magnet, but in the current literature also including antiskyrmion, antiferromagnetic, multipolar, and orbital analogues. The model space is not singular: published formulations range from continuum Ginzburg–Landau and CP\(^1\) field theories to discrete lattice-spin Hamiltonians, effective point-particle and soft-disk reductions, Brownian and Thiele dynamics, and itinerant-electron models in which RKKY, Kondo, Rashba, or scalar-chirality terms generate the noncollinear order. Across these settings, the recurring problem is to determine when exchange, Dzyaloshinskii–Moriya interaction, scalar chirality, anisotropy, and Zeeman coupling stabilize a periodic topological state and how that state melts, deforms, or acquires collective modes [1006.3973], [1304.6580], [2403.00151], [2605.21010]. ## 1. Foundational Hamiltonians The canonical chiral-magnet starting point is a continuum energy for a unit-vector field. In two dimensions, one representative form is \[ {\cal F}[{\bf n}] \;=\;\dfrac{J}{2}\sum_{\mu=x,y}\partial_\mu{\bf n}\!\cdot\!\partial_\mu{\bf n} \;+\;D\,{\bf n}\cdot(\nabla\times{\bf n}) \;-\;{\bf B}\cdot{\bf n}, \] with ferromagnetic exchange \(J\), Dzyaloshinskii–Moriya coupling \(D\), and Zeeman field \({\bf B}=B\hat z\). In the CP\(^1\) rewrite, the spin field is represented as \({\bf n}=z^\dagger\boldsymbol\sigma z\) with \(z^\dagger z=1\), and the resulting mean-field equation can be reduced to an Abrikosov-type problem whose lowest mode generates a triangular skyrmion crystal [1006.3973]. Continuum helimagnet models use the same exchange–DM–Zeeman structure but emphasize symmetry and perturbative departures from circular symmetry. A standard form is \[ \mathcal{H}[\mathbf{m}] = J\!\int\!d^2r\,(\nabla\mathbf{m})^2 + D\!\int\!d^2r\,\mathbf{m}\!\cdot\!(\nabla\times\mathbf{m}) -\mu_0 M_s H\!\int\!d^2r\,m^3, \] for a normalized magnetization field \(\mathbf m(\mathbf r)\), with weak-\(D/J\) perturbation theory predicting systematic asymmetric deformations for both single skyrmions and hexagonal Skyrme crystals [1510.05767]. Discrete lattice-spin Hamiltonians make the microscopic ingredients explicit. In a two-dimensional square lattice of classical spins \(\mathbf s_i\), one model is \[ \mathcal{H} = -\tfrac12\sum_{ij}J_{ij}\,\mathbf{s}_i\!\cdot\!\mathbf{s}_j -H\sum_i s_{iz} -A\sum_i\bigl[(\mathbf{s}_i\times\mathbf{s}_{i+\delta_x})_x+(\mathbf{s}_i\times\mathbf{s}_{i+\delta_y})_y\bigr] -K\sum_i(s_{iz})^2, \] with nearest-neighbor ferromagnetic exchange \(J\), Bloch-type DMI \(A\), perpendicular field \(H\), and uniaxial anisotropy \(K\). For \(A/J=0.2\), \(H/J=-0.025\), and \(K=0\), energy minimization gives a triangular skyrmion lattice of period \(a_S\approx 38.5\,a\) at \(T=0\) [2403.00151]. Three-dimensional chiral magnets are modeled on the cubic lattice by nearest- and next-nearest-neighbor exchange and DM terms together with a Zeeman field, \[ H=H_{nn}+H_{nnn}-H_Z, \] where \(J'=J/16\) and \(K'=K/8\) are introduced to cancel leading lattice-discretization anisotropies. For \(K/J=\tan(2\pi/10)\simeq0.727\), the model reproduces a helical pitch of ten sites and yields the finite-temperature A-phase as a triangular lattice of skyrmion tubes [1304.6580]. ## 2. Reduced and coarse-grained descriptions A major theme in skyrmion-lattice modeling is reduction from a spin field to effective particles. In a dilute two-dimensional skyrmion crystal, well-separated skyrmions can be treated as point particles interacting through a short-range repulsive potential \[ U(d)\simeq F\exp\!\Bigl(-\frac{d}{\delta_H}\Bigr), \qquad F = 60\,J\Bigl(\frac{A^2}{J|H|}\Bigr)^2, \qquad \delta_H=a\sqrt{\frac{J}{|H|}}, \] with skyrmion centers extracted from connected islands of spins satisfying \(s_z>0\). In the square-lattice spin model and in the derived point-particle model, the two descriptions produce similar melting behavior [2403.00151]. Thin-film multilayers admit an even more aggressive coarse-graining. There each skyrmion is replaced by a “soft disk” of characteristic size \(\sigma\) interacting by a purely repulsive power law, \[ V(r)=\epsilon(\sigma/r)^n, \] where \(\sigma\) is fixed by the first peak of the experimental pair-correlation function \(g_{\rm exp}(r)\), and \(n\) is obtained by minimizing a mean-squared deviation \(\chi^2(n)\). For the high-density lattice at \(T=338\,{\rm K}\), \(H=20\,\mu{\rm T}\), the extracted parameters are \(\sigma\simeq1.08\,\mu{\rm m}\), \(n\simeq10.3\pm0.1\), and \(\epsilon\simeq k_B(338\,{\rm K})\), while at \(T=330\,{\rm K}\), \(H=40\,\mu{\rm T}\), one finds \(n\approx9.1\) [2004.09244]. Brownian-dynamics and Thiele-type descriptions focus on collective motion rather than equilibrium energetics. In confinement studies, the steady-state equation of motion is \[ -\,\alpha\,\mathbf{v}_i-G_{\rm rel}\,\hat{\mathbf e}_z\times\mathbf v_i +\boldsymbol\xi_i(t)+\sum_{j\neq i}\mathbf F_{ij}+\mathbf F_i^{\rm conf}=0, \] with damping, Magnus force, thermal noise, pair repulsion, and boundary forces. For static ordering, \(G_{\rm rel}\) is set to zero, the pair interaction is approximated by \(V(r)=1/r^8\) with cutoff \(r_{\rm cut}=1.8\), and the confinement is modeled by a truncated Lennard–Jones repulsion from polygonal or circular boundaries [2508.15758]. Rigid-lattice Thiele models also describe collective phasons. For a synthetic antiferromagnet, coupled Thiele equations show that when the interlayer coupling becomes comparable to or larger than the intralayer coupling, the long-wavelength phason branch changes from quadratic to linear, matching the continuum result that the skyrmion-lattice phason becomes a linear gapless mode after cancellation of the Berry term in the Néel-field dynamics [2209.05893]. ## 3. Order parameters and topological diagnostics The standard structural diagnostic is sixfold bond-orientational order. In particle or center-coordinate descriptions one defines \[ \Psi_i=\frac{1}{6}\sum_{j\in{\rm n.n.}(i)} e^{6i\theta_{ij}}, \qquad O_6=\frac{1}{N_S}\sum_{i=1}^{N_S}\Psi_i, \qquad V_6=\sqrt{\frac{1}{N_S}\sum_i|\Psi_i|^2}. \] For a perfect horizontal hexagon, \(|\Psi_i|=1\), while in the high-temperature liquid \(|\Psi_i|\) fluctuates around \(1/\sqrt6\), giving \(V_6\to1/\sqrt6\) and \(|O_6|\to0\). Solid and liquid phases are then distinguished by \(|O_6|\), the radial distribution \(g(r)\), and Bragg peaks in the structure factor \(S(\mathbf q)\) [2403.00151]. In Voronoi-based analyses, the local order parameter is written as \[ \psi_6(\mathbf r_j)=\frac{1}{n_j}\sum_{k=1}^{n_j} e^{6i\theta_{jk}}, \] with global metrics \(\langle|\psi_6|\rangle=(1/N)\sum_j|\psi_6(\mathbf r_j)|\) and \(\Psi_6=\bigl|\frac1N\sum_j\psi_6(\mathbf r_j)\bigr|\). A perfect monodomain hexagonal lattice sends both quantities toward unity, whereas a polycrystal can retain a large \(\langle|\psi_6|\rangle\) while suppressing \(\Psi_6\) by phase cancellation between domains [2508.15758]. Positional and topological information comes from correlation functions and charge formulas. In two-dimensional melting studies, the pair-correlation function \(g(r)\), the hexatic correlator \[ g_6(r)=\langle\psi_6(\mathbf r_i)\psi_6^*(\mathbf r_j)\rangle_{|\mathbf r_i-\mathbf r_j|=r}, \] and the translational correlator \(g_T(r)\) separate liquid, hexatic, and solid regimes by exponential versus algebraic decay. In calibrated soft-disk simulations, \(\langle|\psi_6|\rangle\approx0.69\) marks the liquid–hexatic crossover, the interval \(0.69\lesssim\langle|\psi_6|\rangle\lesssim0.85\) is associated with a hexatic phase, and the crystal is identified by algebraic decay in both \(g_6(r)\) and \(g_T(r)\) [2004.09244]. The skyrmion number is usually defined by the continuum integral \[ Q=\frac{1}{4\pi}\int d^2r\,\mathbf S\cdot(\partial_x\mathbf S\times\partial_y\mathbf S), \] or by a lattice discretization based on signed spherical-triangle areas. In three-dimensional chiral magnets, the projected plaquette charge \(q_p\) identifies skyrmion tubes in each \(xy\)-plane; in triangular-lattice RKKY systems, the total scalar chirality \(\chi_{\rm tot}\) acts as the order parameter of the triple-\(q\) SkX, while the staggered chirality \(\chi_{\rm stg}\) diagnoses a single-\(q\) spiral [1304.6580], [2112.06414]. ## 4. Thermal phase behavior and melting The thermal behavior of skyrmion lattices is model dependent rather than universal. In a two-dimensional square-lattice spin model and in its point-particle reduction, the skyrmion lattice exhibits a sharp one-step transition between solid and liquid phases on temperature and magnetic field. At fixed \(H/J=-0.025\), the transition occurs at \(T_m\simeq0.118\,J\), where \(|O_6|\) and \(V_6\) drop sharply, the Bragg peaks in \(S(\mathbf q)\) vanish, translational order is lost, and the magnetization \(m_z(T)\) shows a clear kink [2403.00151]. By contrast, thin-film multilayers mapped to repulsive soft disks exhibit the two-step melting structure associated with Kosterlitz–Thouless–Halperin–Nelson–Young theory. In the \((\rho,T^*)\) plane, the liquid–hexatic and hexatic–solid boundaries are approximately \[ \rho_{LH}(T^*)\sigma^2 \simeq 0.83 + 0.15(1-T^*),\qquad \rho_{HS}(T^*)\sigma^2 \simeq 0.96 + 0.05(1-T^*), \] with a small region of two-phase liquid–hexatic coexistence when \(n\gtrsim6\). This directly connects skyrmion lattices in multilayers to the soft-disk melting problem [2004.09244]. Three-dimensional chiral magnets add a further layer of fluctuation physics. Classical Monte Carlo on the cubic lattice shows that the skyrmion-lattice A-phase occupies a small pocket \(0.12\,J\lesssim B\lesssim0.22\,J\) and \(0.82\,J\lesssim T\lesssim0.88\,J\), with \(T_c(B=0)\approx0.92\,J\). The A-phase vanishes as \(T\to0\) and as \(T\to T_c\), and its stabilization is attributed to thermal fluctuations beyond Gaussian order, in qualitative agreement with the Brazovskii scenario for the helical–paramagnetic transition [1304.6580]. A third pattern appears in the \(SO(3)\)-symmetric triangular-lattice SkX-2 model with scalar chirality. There the finite-temperature behavior depends on charge density: the lattice-locked SkX-2 shows a first-order transition, while the free SkX-2 exhibits a sharp, size-independent peak in \(C(T)\), a crossover of spin–spin correlations from exponential to algebraic, and proliferation or binding of \(5\)–\(7\) dislocation pairs consistent with a single Kosterlitz–Thouless–Halperin–Nelson–Young-type melting to a floating solid [2605.21010]. Field gradients generate coexistence phenomena. When a static gradient \(\nabla H\parallel\hat x\) is applied to a two-dimensional skyrmion system, the local melting temperature becomes \(T_m[H(x)]\), and for \(\min T_m < T < \max T_m\) a stationary solid–liquid interface appears near \(H(x)=H_c\) such that \(T_m[H_c]=T\). In the point-particle model, \(|O_6|(x)\) and \(V_6(x)\) fall from \(\approx1\) to \(\approx1/\sqrt6\) over a width of order a few \(a_S\) [2403.00151]. ## 5. Itinerant, frustrated, and zero-field stabilization mechanisms Not all skyrmion-lattice models rely on explicit Dzyaloshinskii–Moriya interactions. In the isotropic RKKY Heisenberg model on the triangular lattice, \[ H=-\sum_{iK_c\simeq0.42\,J\). For \(K/J=0.7\), simulated annealing yields nanoscale skyrmions at vanishingly small DMI and zero magnetic field, together with field-stabilized skyrmion-lattice phases [1911.03171]. Scalar-chirality terms provide another zero-field route. On the triangular lattice, the \(SO(3)\)-symmetric Hamiltonian \[ \mathcal H=J\sum_{\langle i,j\rangle}\mathbf S_i\!\cdot\!\mathbf S_j +K\sum_{(i,j,k)\in\Delta}\mathbf S_i\!\cdot\!(\mathbf S_j\times\mathbf S_k) \] stabilizes SkX-2, a skyrmion crystal carrying two topological charges per magnetic unit cell and no magnetization at the ferromagnetic point in reciprocal space. Near the onset, a continuum expansion gives a Landau functional in the skyrmion density \(q\) with a continuous onset \(q\propto(K/J+\sqrt3/2)\) [2605.21010]. ## 6. Generalized skyrmion lattices and collective modes Antiferromagnetic skyrmion-lattice models alter both topology and dynamics. In an intrinsic triangular-lattice antiferromagnet with \(J_1<0\), \(J_3>0\), third-neighbor DMI \(\mathbf D_{ij}=D_3\,\hat z\times\hat r_{ij}\), anisotropy \(K\), and field \(h\), isolated and interlinked AFM skyrmions are stabilized on top of a row-wise AFM background. Each ferromagnetic sublattice carries an integer \(\pm1\) charge, while the total charge over the four sublattices vanishes, implying zero net Magnus force [2302.14398]. A distinct antiferromagnetic route begins from itinerant electrons with Rashba spin–orbit coupling and strong Hund’s coupling. After a second-order expansion in the insulating limit, the resulting classical spin Hamiltonian contains antiferromagnetic superexchange, a Dzyaloshinskii–Moriya term proportional to \(t\lambda/J_H\), and a bond-anisotropy term proportional to \(\lambda^2/J_H\). Monte Carlo on the triangular lattice finds two field-stabilized AF-skyrmion crystals, AF-SkX1 and AF-SkX2, which are distinguished by peaks at distinct momenta in the spin structure factor [2102.07563]. Synthetic antiferromagnets emphasize collective excitations rather than static topology alone. In the three-coupled-helix continuum theory, a single ferromagnetic skyrmion lattice has a Berry-term-induced quadratic gapless phason, but in the bilayer SAF the Berry term is cancelled and the effective Lagrangian becomes \[ \mathcal L_{\rm SAF} =\int d^2r\, \Bigl\{ \tfrac{\tilde m_\varphi}{2}(\dot\varphi_+^2+\dot\varphi_-^2) -\tilde A_{ex}\bigl[(\nabla\varphi_+)^2+(\nabla\varphi_-)^2\bigr] \Bigr\}, \] which yields a linear gapless mode. The same work shows nonzero spin pumping for arbitrary microwave frequency in the lattice case because the excitation spectrum is gapless [2209.05893]. The notion of a skyrmion lattice also extends beyond ordinary magnetic dipoles. In non-Kramers \(\Gamma_{3g}\) doublet systems on the triangular lattice, bond-dependent hoppings generate a compass-like anisotropy in the effective RKKY interaction, stabilizing zero-field multipolar skyrmion crystals. Two phases are reported: MSkX-I with \(N_{\rm sk}=\pm1\) and MSkX-II with \(N_{\rm sk}=\pm2\), and the total charge fractionalizes locally into meron and antimeron composites subject to charge conservation within the unit cell [2404.14404]. In optical-lattice boson systems, the relevant degree of freedom is orbital pseudospin rather than magnetization. Spinless bosons in degenerate \(p\)-orbital bands of a triangular lattice are described by a multi-orbital Bose–Hubbard model whose strong-coupling expansion yields an anisotropic orbital-exchange Hamiltonian. Bosonic dynamical mean-field theory and exact diagonalization show a periodic orbital texture with skyrmion number \(Q=\pm1\) per magnetic unit cell, arising from the interplay of \(p\)-orbital symmetry and geometric frustration [2212.00167]. Symmetry also permits non-skyrmionic topological lattices. For \(D_{2d}\)-type Dzyaloshinskii–Moriya interaction, a three-dimensional lattice-spin model supports an antiskyrmion lattice pocket over \(0.05\lesssim\hat B\lesssim0.30\) and \(0.6\lesssim\hat T\lesssim1.1\) at \(\hat K=\tan(2\pi/10)\approx0.7265\). The same study treats Bloch and Néel skyrmion lattices within the same formalism, while perturbative helimagnet analysis predicts intrinsic sixfold anisotropies and asymmetric lattice distortions rather than perfectly isotropic unit cells [2109.15020], [1510.05767]. ## 7. Interpretation and scope of the model concept The published usage of “skyrmion lattice model” therefore denotes a modeling framework rather than a unique Hamiltonian. At one end are microscopic spin models with exchange, DM interaction, anisotropy, and field; at the other are reduced particle or continuum descriptions that preserve the lattice symmetry, topological charge, and relevant elastic or thermodynamic observables. Between these limits lie itinerant and multipolar formulations where the periodic topological order is selected indirectly by electronic, orbital, or chirality-mediated interactions [2403.00151], [2004.09244], [2111.13976], [2404.14404]. Two points recur across otherwise different constructions. First, the experimentally visible lattice is usually triangular, but the mechanism of stabilization is not unique: it may be field-induced in a chiral magnet, fluctuation-stabilized near \(T_c\), generated by frustration and emergent many-body terms in RKKY systems, tied to \(2k_F\) in Kondo-lattice models, or realized at zero field through scalar chirality, compass anisotropy, or orbital exchange [1304.6580], [2112.06414], [2605.21010], [2212.00167]. Second, the thermal fate of the lattice is likewise model specific: one-step solid–liquid melting, two-step KTHNY-type melting through a hexatic regime, first-order locking transitions, and floating-solid behavior all occur in the current literature [2403.00151], [2004.09244], [2508.15758], [2605.21010]. A plausible implication is that “skyrmion lattice model” is best understood as a topological-crystal modeling class whose defining ingredients are a periodic noncoplanar texture, a quantized charge diagnostic, and a symmetry-respecting effective theory for ordering, melting, and collective motion. The supplied studies show that this class now encompasses chiral magnets, centrosymmetric itinerant metals, synthetic antiferromagnets, non-Kramers multipolar systems, and orbital optical lattices, with corresponding observables ranging from Bragg peaks and magnetization kinks to Hall responses, defect proliferation, and phason spin pumping [1006.3973], [2209.05893], [2404.14404], [2508.15758]. Source: https://www.emergentmind.com/topics/skyrmion-lattice-model