Chiral Soliton Lattice in Magnets & QCD

• Chiral Soliton Lattice is a spatially periodic array of 2π topological solitons formed in a chiral order parameter, observed in both magnetic and QCD systems.
• It emerges from the competition between Dzyaloshinskii–Moriya interactions and exchange/gradient energies in magnets and from anomalous couplings in QCD, yielding tunable soliton density.
• The CSL exhibits distinct transport, dynamical, and spectroscopic responses, with phenomena such as lock-in behavior and current-driven motion offering practical insights.

A chiral soliton lattice (CSL) is a spatially periodic array of topological solitons in a chiral order parameter. In monoaxial chiral magnets, the relevant field is the spin orientation along a distinguished crystallographic axis: a zero-field helix deforms under a magnetic field perpendicular to that axis into wide field-aligned regions separated by localized $2\pi$ twists, and these twists form the lattice. In QCD and related effective theories, the order parameter is instead a neutral pion or more general chiral field, and the CSL appears as a periodic stack of domain walls induced by anomalous couplings to magnetic field, baryon chemical potential, or rotation (Okumura et al., 2017, Brauner et al., 2016, Canfora et al., 13 Oct 2025, Eto et al., 2023).

1. Definition, topology, and symmetry content

In monoaxial chiral magnets such as CrNb$_3$S$_6$, spins favor ferromagnetic alignment but are twisted into a helix along the chiral axis by the Dzyaloshinskii–Moriya (DM) interaction. At zero field the ground state is a uniform helix with fixed pitch. Under a magnetic field perpendicular to the helical axis, the twist ceases to be uniform: the system develops nearly ferromagnetic segments separated by sharp chiral twists, each carrying an approximate $2\pi$ rotation, and the periodic repetition of these twists is the CSL. With increasing field, the soliton density decreases and the lattice period grows until the system crosses over or transitions to a forced ferromagnetic state (Okumura et al., 2017, Matsumura et al., 2017).

A convenient topological measure in the magnetic setting is the winding number

$$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$

which counts the total twist in units of $2\pi$ and, in the CSL regime, effectively counts the number of chiral solitons along the chain (Okumura et al., 2017). In continuum language, this is the discrete analogue of $\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$.

In QCD, the CSL is an inhomogeneous hadronic phase in which the neutral pion field or a related chiral angle forms a periodic soliton configuration. Here the order parameter is a pseudoscalar condensate, so the phase spontaneously breaks parity and continuous translations along the modulation direction down to a discrete subgroup. The same periodic soliton configuration carries baryon number through anomalous couplings, so the CSL is simultaneously a topological and a density-carrying phase (Brauner et al., 2016, Canfora et al., 13 Oct 2025).

A common source of confusion is the meaning of “chiral.” In monoaxial magnets it refers to handedness fixed by the DM interaction and crystal chirality. In the QCD effective-theory context emphasized for the neutral pion CSL, it refers to a parity-violating pseudoscalar condensate rather than merely to the underlying $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$ symmetry (Brauner et al., 2016).

2. Effective-field-theory descriptions

The standard continuum description of a magnetic CSL is a chiral sine-Gordon-type theory for an angular field $\theta(z)$ or $\phi(z)$. In monoaxial helimagnets, the energy contains exchange, a DM term linear in the spatial derivative, and Zeeman coupling to the transverse field. In one continuum formulation used for current-driven dynamics,

$_3$0

with $_3$1 and $_3$2 in the CSL geometry (Osorio et al., 2022). The zero-field period is $_3$3, while the field-dependent CSL period satisfies

$_3$4

so the period diverges as $_3$5 and the soliton density vanishes (Osorio et al., 2022).

For QCD with two light flavors at finite baryon chemical potential $_3$6 and external magnetic field $_3$7, the low-energy neutral-pion Hamiltonian density takes the form

$_3$8

The last term is anomalous and linear in $_3$9; it plays the same structural role that the DM term plays in magnets by favoring a finite twist density (Brauner et al., 2016). The Euler–Lagrange equation is the sine-Gordon equation

$_6$0

with periodic solution

$_6$1

and lattice period

$_6$2

The CSL is energetically preferred when

$_6$3

and in the chiral limit the solution reduces to the linear profile

$_6$4

These results are fully analytic at leading order in chiral perturbation theory (Brauner et al., 2016).

A more microscopic hadronic derivation starts from the gauged Skyrme model coupled to Maxwell theory. Under the ansatz

$_6$5

the full $_6$6-dimensional Skyrme–Maxwell system reduces to the $_6$7-dimensional sine-Gordon equation

$_6$8

and the baryon density is purely of Callan–Witten origin,

$_6$9

so the CSL becomes a lattice of baryonic hadronic layers (Canfora et al., 13 Oct 2025).

3. Microscopic realizations in condensed matter

A direct microscopic realization of the magnetic CSL is provided by a one-dimensional Kondo lattice with classical localized spins, ferromagnetic on-site Hund coupling $2\pi$0, DM interaction $2\pi$1, hopping $2\pi$2, and transverse magnetic field $2\pi$3,

$2\pi$4

For $2\pi$5, $2\pi$6, $2\pi$7, quarter filling, and $2\pi$8, Monte Carlo simulations find a zero-field helix with period 10 sites, a field-induced CSL, and finally a forced ferromagnetic state as $2\pi$9 approaches the DM scale (Okumura et al., 2017).

The field evolution is visible in the spin structure factor

$$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$0

At $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$1, $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$2 peaks at $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$3, indicating a helix of period $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$4. With increasing field, the helical peak moves continuously toward smaller $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$5, while the $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$6 component grows, showing the transformation into a CSL with increasing ferromagnetic component; the CSL period is $$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$7 (Okumura et al., 2017).

Experiments on Yb(Ni$$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$8Cu$$W=\frac{1}{2\pi}\sum_l \langle \Delta\theta_l\rangle,$$9)$2\pi$0Al$2\pi$1 provide direct diffraction evidence for CSL formation. Resonant X-ray diffraction detects a helical ground state with propagation vector $2\pi$2, a one-to-one correspondence between crystal chirality and magnetic helicity, and, under transverse field, the emergence of second-harmonic $2\pi$3 and third-harmonic $2\pi$4 peaks. The growth of these higher harmonics is strong evidence for a CSL, since a perfect sinusoidal helix would carry only the first harmonic (Matsumura et al., 2017). In the $2\pi$5 sample, a commensurate lock-in at $2\pi$6 was observed in field, indicating coupling between the CSL and the lattice (Matsumura et al., 2017).

The role of itinerant electrons in such lock-in phenomena is explicit in a variational Kondo-lattice study, which finds that the CSL period can be locked at particular values dictated by the Fermi wave number, in contrast to spin-only models, and that the same spin-charge coupling can stabilize a spontaneous CSL even at zero magnetic field (Okumura et al., 2018). This provides a plausible microscopic explanation for the lock-in behavior seen in Yb(Ni$2\pi$7Cu$2\pi$8)$2\pi$9Al$\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$0 (Okumura et al., 2018, Matsumura et al., 2017).

At the structural level, Cr$\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$1NbS$\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$2 has been analyzed in terms of long-range exchange pathways. A dominant antiferromagnetic coupling $\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$3 along only one of two crystallographically equivalent diagonals of the trigonal-prism network builds left-handed helices along the $\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$4 axis; weaker inter-helix couplings preserve quasi-one-dimensionality. In that picture, the DM interaction is responsible for final ordering and stabilization of these chiral helices into a CSL (Volkova et al., 2014).

4. Transport, current response, and spectroscopy

The CSL has a distinctive transport signature: in the Kondo-lattice Monte Carlo study, the low-temperature resistivity tracks the soliton number. Using the coherent optical weight $\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$5 extracted from the Kubo conductivity, one finds

$\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$6

through the CSL regime, so negative magnetoresistance is directly proportional to the number of chiral solitons (Okumura et al., 2017). This identifies individual solitons as spin-scattering centers for itinerant electrons.

The dynamical response to spin-polarized current depends on the theoretical setting. In a continuum Landau–Lifshitz–Gilbert treatment of a monoaxial helimagnet under transverse field, sufficiently small current density produces steady CSL motion with velocity proportional to current, and the mobility is independent of soliton density and magnetic field; a small conical distortion accompanies the motion. At larger current density, the spin-transfer torque destabilizes the CSL and drives the system into a ferromagnetic state parallel to the field (Osorio et al., 2022). By contrast, in a classical spin chain coupled to conduction electrons via $\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$7–$\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$8 exchange and analyzed with collective coordinates plus an SU(2) gauge method, the terminal CSL velocity decreases as the CSL period becomes longer because the nonadiabatic force is proportional to the spin-structure-induced resistivity (Tokushuku et al., 2017). Taken together, these results indicate that “current-driven CSL motion” is not a single universal phenomenon but depends on whether the dominant microscopic mechanism is phenomenological spin-transfer torque or explicit electron backaction.

Microwave spectroscopy in CrNb$\frac{1}{2\pi}\int dz\,\partial_z\theta(z)$9S$\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$0 shows that the collective excitations of the CSL are highly sensitive to magnetic disorder. Three resonance modes were observed over a wide frequency range; their predominance depends on field history, and sweeping the external field through the ideal helical state at $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$1 suppresses disorder and restores macroscopic coherence of the CSL (Goncalves et al., 2019). In finite-size CSLs, a standing-spin-wave theory generalizing Kittel–Pincus dynamics predicts two classes of modes: a Pincus mode, consisting of a long-period Bloch wave with a short-period ripple inherited from the CSL, and a short-period Kittel ripple excited only when the ac field is perpendicular to the chiral axis. Their coexistence accounts for the experimentally observed double-resonance profile (Kishine et al., 2019).

5. QCD, QCD-like theories, and extensions beyond magnets

In QCD at nonzero baryon chemical potential and magnetic field, the CSL appears as a periodic neutral-pion condensate that carries baryon density and magnetization per unit cell,

$\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$2

so each lattice period contributes $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$3 and $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$4 (Brauner et al., 2016). The fluctuation problem is analytically tractable: neutral fluctuations are governed by the $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$5 Lamé equation and yield a phonon with

$\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$6

while charged pions can undergo Bose–Einstein condensation at stronger magnetic field. In the chiral limit, the instability threshold is

$\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$7

so there is a window $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$8 in which the neutral-pion CSL is stable (Brauner et al., 2016).

The gauged Skyrme–Maxwell analysis strengthens this picture by showing that the CSL is not an artifact of a truncated low-energy theory. For the ansatz $\mathrm{SU}(N_f)_L\times \mathrm{SU}(N_f)_R$9, one has $\theta(z)$0 and hence

$\theta(z)$1

As a result, generalized Skyrme corrections built from $\theta(z)$2 vanish identically on the CSL configuration, so the sine-Gordon equation, the energy density, and the topological properties of the CSL remain unchanged even after including subleading large-$\theta(z)$3 corrections. This is the precise sense in which the CSL is “universal” in the low-energy limit of QCD coupled to electromagnetism (Canfora et al., 13 Oct 2025).

QCD-like theories with real or pseudoreal color representations provide a different extension. Their low-energy theory contains neutral pion and diquark modes, and at finite baryon chemical potential and magnetic field the phase diagram exhibits vacuum, homogeneous Bose–Einstein condensate, and CSL regions. The onset of the inhomogeneous phase is already visible in the adjacent homogeneous phase through a roton-like minimum in the lowest excitation, and these theories furnish explicit counterexamples to the conjecture that positivity of the Dirac determinant forbids spontaneous breaking of translational invariance (Brauner et al., 2019).

The notion of CSL also generalizes to rotation and to internal non-Abelian structure. In rotating QCD matter with two flavors, the $\theta(z)$4 meson can form an Abelian CSL, and in a large parameter region this becomes a non-Abelian CSL in which each $\theta(z)$5 soliton splits into a pair of non-Abelian sine-Gordon solitons carrying $\theta(z)$6 moduli. Around that non-Abelian CSL there are three gapless type-A Nambu–Goldstone modes—two isospinons and a phonon—and in the deconfined phase the isospinon dispersion becomes Dirac type, linear even at large momentum (Eto et al., 2023).

Two further extensions clarify how broadly the concept propagates across field theory. A supersymmetric chiral sine-Gordon model supports a CSL ground state in the presence of strong magnetic field and/or large chemical potential, or via a fermion bilinear condensate built from the gaugino and the superpartner of a baryon gauge field (Nitta et al., 2024). In the skyrmion-crystal description of dense baryonic matter, introducing a CSL of the neutral pion induces an inverse catalysis of the skyrmion–half-skyrmion topology change, deforms the single-baryon profile into a highly oscillating structure, enhances the baryon energy, and disappears again in the high-density region after the topology change (Kawaguchi et al., 2018).

6. Conceptual clarifications, limitations, and open directions

A CSL is not merely any helical or spiral state. The decisive feature is the nonuniform twist: the soliton lattice consists of extended nearly aligned regions separated by localized topological walls. The uniform helix at zero field and the forced ferromagnet at high field are the two limiting states between which the CSL interpolates (Okumura et al., 2017, Matsumura et al., 2017).

Another misconception is that the CSL period must vary smoothly with field. Spin-only continuum models do predict a continuously varying period, but itinerant electrons can qualitatively change this conclusion by locking the period to values dictated by the Fermi wave number, and experiments on Yb(Ni$\theta(z)$7Cu$\theta(z)$8)$\theta(z)$9Al$\phi(z)$0 show commensurate lock-in behavior consistent with that mechanism (Okumura et al., 2018, Matsumura et al., 2017).

Theoretical limitations remain substantial. The Monte Carlo transport study uses a strictly one-dimensional model with classical localized spins, open boundaries, a single conduction band, and no explicit electronic spin–orbit coupling; in that setting the field-driven evolution should be interpreted as a crossover rather than a sharp thermodynamic phase transition at finite temperature, and a three-dimensional model is required to address true finite-$\phi(z)$1 critical behavior (Okumura et al., 2017). Continuum current-driven theories capture field-tunable soliton transport but neglect realistic disorder, domain structure, and fully microscopic electron dynamics (Osorio et al., 2022, Tokushuku et al., 2017). Spectroscopy in finite crystals shows that boundary pinning and field history are not peripheral complications: they are intrinsic determinants of the observed mode structure (Kishine et al., 2019, Goncalves et al., 2019).

In QCD and related theories, the open questions are of a different type. The existence of the CSL at low energy is analytically controlled in several settings, but quantitative placement in realistic phase diagrams depends on medium-dependent effective couplings, charge-neutrality constraints, and competition with other inhomogeneous phases. The fact that sign-problem-free QCD-like theories admit CSL phases suggests that lattice Monte Carlo can test anomalous translation-breaking phases directly, but the concrete realization of such simulations remains an open program (Brauner et al., 2019).

Across condensed matter and hadronic physics, the unifying lesson is precise: the CSL is the ordered phase generated when a chiral derivative term—DM in magnets, anomalous Wess–Zumino or Callan–Witten coupling in QCD—competes with exchange or gradient energy, a periodic potential, and external control parameters such as magnetic field, chemical potential, or rotation. Its importance lies in that combination of topology, tunable periodicity, and measurable dynamical and transport response (Okumura et al., 2017, Brauner et al., 2016, Canfora et al., 13 Oct 2025).