Chiral Soliton Lattices under Magnetic Fields and Rotation: a Holographic Analysis

Abstract: We study the chiral soliton lattice (CSL) in rotating QCD matter under a background magnetic field within holographic QCD. We show that rotation can be incorporated as a background gauge field and that the CSL in rotating matter becomes a ground state in the gravity dual. We also provide a brane interpretation of the CSL under a magnetic field and rotation. Solving the bulk equations, we derive an effective Hamiltonian where the meson decay constant $\tilde{f}$ emerges as an anisotropic, field-dependent matrix. Furthermore, we discuss the thermodynamic properties of the ground state in rotating matter and study its equation of state.

• The paper shows that rotation can be modeled at leading order as an effective baryonic magnetic field, producing a rotational chiral soliton lattice above a critical combination of chemical potential and angular velocity.
• For two flavors, magnetic and vortical effects generate four possible phases—vacuum, two single-CSL states, and a dual-CSL state—controlled by the competing field combinations B₊ and B₋.
• The full holographic calculation produces an anisotropic, field-dependent meson decay-constant matrix, shifting chiral perturbation theory phase boundaries and yielding a strong-field threshold that scales with the square of the pion mass.

Overview

This paper extends the holographic analysis of chiral soliton lattices (CSLs) in dense QCD to matter subject simultaneously to a background magnetic field and rotation, using the Sakai–Sugimoto (Witten–Sakai–Sugimoto, WSS) model. The central technical device is the identification of rotation, at leading order in the local velocity $|\vec v| = |\vec\Omega \times \vec x|$, with a spatial component of a background $\mathrm{U}(1)_B$ gauge field whose effective magnetic field is $\vec B_{\rm eff} = 2\mu_B\vec\Omega$ — the relativistic Barnett effect. This allows magnetic and vortical anomaly-induced effects to be treated on a common footing within a single five-dimensional bulk computation. The paper builds directly on the authors' earlier holographic study of the neutral-pion CSL in purely magnetic backgrounds (Amano et al., 22 Jul 2025) and their systematic re-derivation of Wess–Zumino–Witten (WZW) couplings of neutral mesons under magnetic fields and rotation (Amano et al., 6 Jul 2026).

The main results are: (i) for $N_f=1$, the $\eta$ meson obeys a chiral sine-Gordon equation and forms a rotational CSL above a critical value of $\mu_B^2\Omega_3$; (ii) for $N_f=2$, the coupled $\pi^0$–$\eta$ system decomposes into two decoupled sine-Gordon sectors with distinct phase structure; (iii) a brane interpretation is given in terms of dissolved D6-brane charge and D4-brane charge sourced by five-dimensional instanton density; and (iv) solving the full bulk gauge-field equations including the tower of massive vector mesons yields an anisotropic, background-dependent matrix of meson decay constants $\tilde f_{ij}$ that shifts the CSL phase boundaries relative to chiral perturbation theory (ChPT).

Rotation as a background gauge field

The incorporation of rotation follows from applying a local Lorentz boost with velocity $\mathrm{U}(1)_B$0 to the rest-frame chemical potential configuration $\mathrm{U}(1)_B$1. To leading order in $\mathrm{U}(1)_B$2, this produces a gauge field with spatial components whose curl gives an effective magnetic field $\mathrm{U}(1)_B$3. In the holographic setup, the baryon-number gauge field is taken as $\mathrm{U}(1)_B$4, so the induced baryonic "magnetic" field is $\mathrm{U}(1)_B$5. The authors emphasize that this identification holds only at $\mathrm{U}(1)_B$6, which bounds the validity of all subsequent results.

For $\mathrm{U}(1)_B$7, the gauged WZW term reduces to a term linear in $\mathrm{U}(1)_B$8:

$\mathrm{U}(1)_B$9

which does not affect the equation of motion but tilts the energy landscape. Combined with the DBI kinetic term and the meson mass term introduced via extra D6-branes (0803.4192, 0803.3547), the $\vec B_{\rm eff} = 2\mu_B\vec\Omega$0 dynamics becomes a chiral sine-Gordon model. The periodic soliton solution is expressed through the Jacobi amplitude function, with elliptic modulus $\vec B_{\rm eff} = 2\mu_B\vec\Omega$1 interpolating between isolated kinks ($\vec B_{\rm eff} = 2\mu_B\vec\Omega$2) and the dilute limit. Minimizing the energy per period yields the existence condition

$\vec B_{\rm eff} = 2\mu_B\vec\Omega$3

so the rotational CSL becomes the ground state above a critical combination of chemical potential and angular velocity, consistent with earlier ChPT analyses of rotating matter (1711.02190, Nishimura et al., 2020).

Brane interpretation

A notable feature of the paper is the Ramond–Ramond charge analysis of the rotating CSL. The sine-Gordon kink induces a localized $\vec B_{\rm eff} = 2\mu_B\vec\Omega$4 flux on the D8-brane worldvolume, sourcing D6$\vec B_{\rm eff} = 2\mu_B\vec\Omega$5-brane charge. The constant baryonic magnetic field $\vec B_{\rm eff} = 2\mu_B\vec\Omega$6 sources an additional, non-quantized D6$\vec B_{\rm eff} = 2\mu_B\vec\Omega$7 charge — interpreted as dissolved D6-branes uniformly distributed in the $\vec B_{\rm eff} = 2\mu_B\vec\Omega$8 plane rather than as independent probe branes. Their intersection generates D4$\vec B_{\rm eff} = 2\mu_B\vec\Omega$9 charge via the $N_f=1$0 coupling, i.e., nonzero instanton density in the $N_f=1$1 directions, which also encodes the baryon number carried by the CSL.

The authors highlight a qualitative contrast with the purely electromagnetic case of ref. (Amano et al., 22 Jul 2025): there, the D6-brane charges cancel, leaving only dissolved D4-brane charge, whereas under rotation both D6 and D4 charges survive. This constitutes a genuinely distinct geometric realization of the CSL in the gravity dual. The stability of the kink-as-D6 configuration relies on the quark-mass deformation shifting the D8–D6 open-string tachyon mass positive, following the mechanism established for QCD domain walls (Argurio et al., 2018). The authors note explicitly that this brane interpretation is valid only within the leading-order-in-$N_f=1$2 prescription in which vorticity is encoded at the boundary.

Two-flavor phase structure

For $N_f=1$3, the electromagnetic field enters as the Cartan generator of $N_f=1$4 with charge matrix $N_f=1$5. Defining $N_f=1$6, the equations of motion decouple into two independent sine-Gordon sectors:

$N_f=1$7

Each sector experiences an effective topological field $N_f=1$8 and $N_f=1$9, so magnetic and vortical effects combine linearly but differently in each sector. The ChPT-level condition for a CSL in sector $\eta$0 is

$\eta$1

with positive $\eta$2 producing a CSL and negative $\eta$3 an anti-CSL. Four phases result: the QCD vacuum, two single-CSL phases (in which a residual U(1) locks $\eta$4), and a dual-CSL phase where both sectors condense. Because $\eta$5 and $\eta$6 respond oppositely to the electromagnetic field, rotation can compensate a strong magnetic field in one sector while leaving the other unaffected — a competition absent in the purely magnetic problem. The net D6-brane charge of the kink array vanishes for (kink, anti-kink) combinations, in which case only the $\eta$7 CSL survives; for aligned kinks the $\eta$8 CSL carries nonzero D6 charge.

Bulk analysis and the anisotropic decay constant

The most technically substantive part of the paper solves the full five-dimensional gauge-field equations on the D8-branes, retaining the tower of massive vector mesons. Although vector modes do not contribute at two-derivative order in vacuum, the WZW term mixes them with the zero modes when $\eta$9 or $\mu_B^2\Omega_3$0. Working in the $\mu_B^2\Omega_3$1 basis, where the equations of motion decouple, the massless bulk solutions are obtained analytically in closed form in terms of hyperbolic functions of $\mu_B^2\Omega_3$2.

Evaluating the on-shell action yields a four-dimensional Hamiltonian of the form

$\mu_B^2\Omega_3$3

where the decay constant is replaced by the anisotropic matrix

$\mu_B^2\Omega_3$4

with $\mu_B^2\Omega_3$5 the projector transverse to $\mu_B^2\Omega_3$6. For weak fields, $\mu_B^2\Omega_3$7; for strong fields it grows linearly, $\mu_B^2\Omega_3$8. This generalizes the scalar, field-dependent decay constant found in the purely magnetic case (Amano et al., 22 Jul 2025), and its origin is traced to the backreaction of the bulk gauge fields rather than to any low-energy input.

The minimized ground-state gradients are parallel to $\mu_B^2\Omega_3$9, recovering the previous single-component results when coordinates are aligned accordingly. The minimum energy density interpolates between quadratic behavior at small fields, $N_f=2$0, and linear behavior at large fields, $N_f=2$1 — the latter independent of $N_f=2$2. A quantitative consequence is a saturation effect: magnetization peaks at $N_f=2$3 and decreases beyond, equivalently the angular momentum susceptibility asymptotes to a constant at strong fields.

Phase boundaries: comparison with ChPT

In the massive case, the bulk equations reduce to a sine-Gordon equation for the boundary fields plus a separable inhomogeneous ODE for $N_f=2$4, solved semi-analytically via Fourier decomposition along $N_f=2$5 (using the known Fourier series of $N_f=2$6) and Chebyshev decomposition along $N_f=2$7. Numerical evaluation of the on-shell Hamiltonian reproduces the same four-phase topology as ChPT, but with boundaries shifted by the bulk backreaction encoded in $N_f=2$8.

Two analytic regimes emerge for small pion mass. At weak fields, expanding to $N_f=2$9 gives the threshold

$\pi^0$0

which is about 11% higher than the ChPT boundary — reflecting the crude approximation $\pi^0$1 used for the amplitude function; multiplying the mass-term contribution by $\pi^0$2 restores agreement with ChPT at small fields. At strong fields the threshold scales instead as

$\pi^0$3

quadratic in $\pi^0$4 and falling as $\pi^0$5, with explicit $\pi^0$6 enhancement inherited from the Chern–Simons coupling. The authors stress that this $\pi^0$7 scaling is harsher than the $\pi^0$8 behavior of ChPT, so the WSS and ChPT phase boundaries differ significantly near the pion-mass scale even though they coincide at large chemical potential where the anomalous term dominates.

For large pion mass, numerics indicate the minimizing modulus approaches $\pi^0$9; approximating the kink profile by a piecewise-linear function yields the stability condition $\eta$0, implying that no CSL$\eta$1 exists unless $\eta$2. Numerics further suggest a maximum pion mass beyond which the CSL is unstable at any field strength, with the phase boundary becoming horizontal in the $\eta$3 plane.

An additional finding concerns competing fields: offsetting magnetic field against counter-rotation can drive one of $\eta$4 toward zero, enlarging the CSL regions under strong backgrounds. The authors flag two caveats here — the effect is relatively small, and because the analysis sits in the large-$\eta$5 limit, it remains unclear whether the enhancement reflects QCD or is an artifact of the approximation.

Limitations and open questions

Several limitations are acknowledged in the paper itself. The rotation-as-gauge-field prescription is valid only to leading order in the local velocity and ignores the fully geometric treatment of rotating spacetimes (Chen et al., 2020, Braga et al., 2022); a complete geometric analysis remains open. The massive bulk equations were solved semi-analytically with numerical evaluation, and a systematic treatment in the regime where $\eta$6, $\eta$7, and $\eta$8 are all comparable has not been performed. Only neutral mesons are included; charged pions, which may introduce additional phase structure, are excluded throughout. The connection between the CSL and the domain-wall Skyrmion phase under simultaneous magnetic and rotational backgrounds is not addressed. Finally, finite-temperature effects and the large-$\eta$9 dependence of the competing-field enhancement are unresolved.

Conclusion

This work unifies magnetic and vortical anomaly-induced physics in a single holographic framework by encoding rotation as an effective $\tilde f_{ij}$0 magnetic field proportional to $\tilde f_{ij}$1. It establishes the rotational CSL as a ground state of the Sakai–Sugimoto dual for $\tilde f_{ij}$2, identifies a four-phase structure for $\tilde f_{ij}$3 governed by the combinations $\tilde f_{ij}$4, provides a brane interpretation in which rotation leaves residual dissolved D6-brane charge alongside D4-brane charge, and derives an anisotropic, field-dependent meson decay constant matrix from the full bulk dynamics. The resulting phase boundaries agree with ChPT at weak fields up to controlled corrections but deviate qualitatively at strong fields, where the threshold scales as $\tilde f_{ij}$5 rather than $\tilde f_{ij}$6. The analysis is bounded by its leading-order treatment of rotation, its restriction to neutral mesons, and the large-$\tilde f_{ij}$7 limit, each of which defines a concrete question left open by the paper.