---
title: Baryonic Vortex Lattice in Dense QCD
url: https://www.emergentmind.com/topics/baryonic-vortex-lattice
type: topic
---

# Baryonic Vortex Lattice in Dense QCD

Searching arXiv for the cited works and closely related papers.
A baryonic vortex lattice is a proposed phase of low-energy dense QCD in which baryon number is carried by topological vortex configurations rather than by conventional localized baryon fields. In two-flavor chiral perturbation theory with finite isospin chemical potential, finite baryon chemical potential, finite pion mass, and dynamical electromagnetism, the relevant ground state for $\mu_I>m_\pi$ contains pion condensation and supports vortices in both charged and neutral pion sectors. The recent linked-vortex formulation identifies each baryon with a topological linking between a local $\pi^\pm$ Abrikosov–Nielsen–Olesen-like vortex and a closed $\pi^0$ global vortex line attached to a domain wall, while earlier work emphasized a vortex-Skyrmion state in which neutral-pion modulation inside a charged vortex core endows the vortex with baryon number. When the baryon chemical potential exceeds a critical value, these objects become energetically favorable and can organize into a triangular Abrikosov-type lattice, yielding a candidate dense baryonic phase of hadronic matter [2509.20844] [2403.07433].

## 1. Effective-field-theory setting

The construction is formulated in two-flavor Chiral Perturbation Theory at leading order in the derivative expansion. The pion field is parametrized by an $SU(2)$ matrix
\[
U(x)=\exp\Big(i\,\tau^a\pi^a(x)/f_\pi\Big)\in SU(2),
\]
with Pauli matrices $\tau^a$. The effective Lagrangian is written as
\[
L_{\rm eff}=L_{\rm chiral}+L_{\rm EM}+L_{\rm WZW},
\]
where
\[
L_{\rm chiral}
=\frac{f_\pi^2}{4}\,\mathrm{Tr}\!\left[D_\mu U^\dagger D^\mu U\right]
+\frac{f_\pi^2 m_\pi^2}{4}\,\mathrm{Tr}\!\left[U+U^\dagger-2\right],
\]
\[
L_{\rm EM}=-\frac{1}{4e^2}F_{\mu\nu}F^{\mu\nu},
\qquad
L_{\rm WZW}=\bigl(A_\mu^B+q\,A_\mu^{\rm EM}\bigr)j_B^\mu,
\qquad q=\tfrac12.
\]
The covariant derivative is
\[
D_\mu U=\partial_\mu U-iA_\mu[Q,U],
\qquad
Q=\frac{1}{6}\mathbf1+\frac{\tau^3}{2},
\]
with the combination $A_\mu$ including the static isospin chemical potential as $A_0=\mu_I$ plus the dynamical photon. A fictitious baryon gauge field implements the baryon chemical potential through $A_0^B=\mu_B$ [2509.20844].

The topological input is the conserved baryon current. In the absence of gauge fields it is
\[
j_B^\mu
=\frac{1}{24\pi^2}\epsilon^{\mu\nu\rho\sigma}
\mathrm{Tr}\!\left[
(U^{-1}\partial_\nu U)(U^{-1}\partial_\rho U)(U^{-1}\partial_\sigma U)
\right],
\]
up to covariant gauge-field corrections. Integrating $j_B^0$ over space yields an integer Skyrmion number. In the complementary notation used for the vortex-Skyrmion formulation, $\Sigma\in SU(2)$ with $Q=\mathrm{diag}(2/3,-1/3)$, and the same low-energy description couples chiral fields, electromagnetism, and the Goldstone–Wilczek current to finite $\mu_I$ and $\mu_B$ [2403.07433].

## 2. Homogeneous pion-condensed background

The existence of the baryonic vortex lattice relies on the structure of the homogeneous ground state at finite isospin density. For $\mu_I<m_\pi$, the minimum of the chiral Lagrangian is the trivial vacuum,
\[
U=1 \qquad (\pi^a=0).
\]
Once $\mu_I>m_\pi$, charged pions condense in the direction $\tau^1\pm i\tau^2$, and one may write the uniform condensate as
\[
U_{\rm cond}=\exp\bigl(i\,\alpha_0\,\tau_1\bigr),
\qquad
\cos\alpha_0=\frac{m_\pi^2}{\mu_I^2},
\qquad
\sin\alpha_0=\sqrt{1-\frac{m_\pi^4}{\mu_I^4}}.
\]
In the basis
\[
\phi_1=(\sigma+i\pi^0),\qquad \phi_2=\pi_2+i\pi_1,
\]
the condensates satisfy
\[
\langle|\phi_2|\rangle=\sin\alpha_0,\qquad
\langle|\phi_1|\rangle=\cos\alpha_0>0
\qquad (\mu_I>m_\pi).
\]
Thus both neutral and charged pions have nonzero expectation values in the bulk once $\mu_I$ exceeds $m_\pi$ [2509.20844].

This coexistence is structurally important because it allows two different topological sectors to coexist in the same medium. The charged-pion condensate breaks $U(1)_{\rm EM}$ spontaneously and supports local flux-carrying vortices. Simultaneously, the nonzero neutral-pion condensate leaves a phase modulus that supports global vortices, with explicit breaking by the pion mass converting the neutral sector into a vortex–domain-wall system. A plausible implication is that the baryonic lattice is not an auxiliary add-on to the pion-condensed phase but an instability of that phase once $\mu_B$ becomes sufficiently large.

## 3. Local, global, and linked vortex configurations

In the phase $\mu_I>m_\pi$, the charged condensate supports an Abrikosov–Nielsen–Olesen-like local vortex. For a straight unit-winding vortex along the $z$-axis, the standard cylindrical ansatz is
\[
\phi_2(\rho,\theta)=f(\rho)e^{i\theta},
\qquad
A_\theta(\rho)=\frac{a(\rho)}{\rho},
\]
with boundary conditions
\[
f(0)=0,\quad a(0)=0,\qquad
f(\infty)=\sin\alpha_0,\quad a(\infty)=-1.
\]
In the simplest pure ANO solution one holds the neutral-pion phase $\beta$ constant, so that the azimuthal baryon current vanishes, $j_B^\theta=0$ [2509.20844].

Because $\phi_1=\cos\alpha_0\,e^{i\beta}$ remains nonzero in the bulk, the neutral sector supports a global vortex of winding $m\in\mathbb Z$,
\[
\phi_1(\rho,\theta)=|\phi_1(\rho)|\,e^{i\,m\theta},
\qquad
|\phi_1(\infty)|=\cos\alpha_0,
\qquad
|\phi_1(0)|=0.
\]
At finite pion mass, the global vortex line is attached to a two-dimensional domain wall, the chiral soliton, across which $\beta$ jumps by $2\pi$. Locally the wall is a sine-Gordon soliton solving
\[
\partial_z^2\beta=m_\pi^2\sin\beta,
\]
with a kink interpolating $\beta:0\mapsto2\pi$ [2509.20844].

The baryonic configuration of the linked-vortex proposal is obtained by superposing the local charged vortex with a closed neutral-pion vortex ring that encircles it and is capped by a domain wall. Earlier work described a related baryonic vortex phase from a different angle: a vortex with the same quantized magnetic flux as the conventional Abrikosov vortex carries baryon number once one takes into account a modulation of the neutral pion inside the vortex core, and Qiu and Nitta therefore dubbed the resulting object the baryonic vortex [2403.07433].

## 4. Topology, baryon number, and Skyrmion interpretation

The topological classification follows from $\pi_3(SU(2))\simeq \mathbb Z$ in the chiral formulation and, equivalently, from $\pi_3(S^3)$ in the $(\sigma,\pi)$ description. For any smooth configuration approaching a constant at spatial infinity, the baryon number is
\[
B=\int d^3x\,j_B^0
=\frac{1}{24\pi^2}\int d^3x\,\epsilon^{0ijk}
\mathrm{Tr}\!\left[
(U^{-1}\partial_iU)(U^{-1}\partial_jU)(U^{-1}\partial_kU)
\right]
\in\mathbb Z.
\]
In the linked-vortex construction, this integer is identified directly with the linking number,
\[
B={\rm Link}(\pi^\pm\hbox{-vortex},\,\pi^0\hbox{-loop})\in\mathbb Z.
\]
Because each linked pair carries $B=1$, it is interpreted as a single Skyrmion realized without any Skyrme term; the stability is attributed to the Wess–Zumino–Witten-induced coupling to $\mu_B$ [2509.20844].

In the vortex-Skyrmion formulation, a periodic modulation of the neutral pion along the vortex axis plays an analogous role. For the ansatz
\[
\Sigma(\rho,\phi,z)
=\begin{pmatrix}
e^{i\varphi(z)}\cos\alpha(\rho)&e^{-i\phi}\sin\alpha(\rho)\\[6pt]
e^{i\phi}\sin\alpha(\rho)&e^{-i\varphi(z)}\cos\alpha(\rho)
\end{pmatrix},
\]
imposing $\alpha(0)=0$, $\alpha(\infty)=\pi/2$, and $\varphi(0)=0$, $\varphi(L)=2\pi$ yields one unit of baryon number per period $0\le z\le L$ [2403.07433].

Two points are often conflated. First, the baryonic vortex lattice is not merely an Abrikosov lattice of ordinary $\pi^\pm$ flux tubes; the baryon number arises from nontrivial neutral-pion structure. Second, the 2025 linked-vortex construction explicitly emphasizes that the Skyrmion-type object is realized without the Skyrme term, whereas the 2024 effective-theory exposition notes that beyond $\mathcal O(p^2)$ one may add the usual Skyrme term for stability. This suggests that the central topological mechanism is not identical to a conventional Skyrme-model stabilization, even though the resulting objects carry integer baryon charge [2509.20844] [2403.07433].

## 5. Energetics, lattice formation, and phase structure

The energetics are organized around the competition between the cost of forming a vortex and the gain from the baryon chemical potential. In the linked-vortex formulation, the energy of one baryonic vortex is
\[
E_{\rm vortex}
=E_{\rm chiral}+E_{\rm EM}+E_{\rm WZW}
=E_{\rm core}+E_{\rm wall}+(-\mu_B)\times 1,
\]
and becomes negative when $\mu_B$ exceeds a critical value $\mu_B^c(\mu_I)$. Physically, this means that it becomes favorable to nucleate the linked pairs once $\mu_B$ is large enough. To construct a dense phase, one arranges many identical baryonic vortices on a two-dimensional lattice in the transverse plane and extends them uniformly in $z$. The unit-cell energy density is then computed as a function of $\mu_I$, $\mu_B$, $m_\pi$, and the lattice spacing $a$ or vortex density. Above $\mu_B>\mu_B^c(\mu_I)$, the unit-cell energy density is lower than that of a homogeneous charged-pion condensate or an Abrikosov lattice of pure $\pi^\pm$ vortices. Varying the lattice geometry, the lowest unit-cell energy is typically achieved by a triangular (hexagonal) lattice, and the optimal spacing $a_{\rm opt}$ decreases as $\mu_B$ grows [2509.20844].

The earlier baryonic-vortex-phase analysis expresses the same instability in terms of the string tension of a single vortex on top of the homogeneous $\pi^\pm$ condensate,
\[
T(k;\mu_I,\mu_B)=\frac{1}{L}\Bigl(H[\alpha,a,k]-H_0\Bigr),
\]
with the term $-\mu_B k$ arising from $-\mu_B\int j_B^0\,d^2x$. For $\mu_B$ below a certain $\mu_B^{\rm crit}$ the minimum tension remains positive, while above
\[
\mu_B^{\rm crit}(\mu_I)=\min_k \frac{T_{2+4+\mathrm{EM}}(k;\mu_I)}{k},
\]
one finds negative tension, so the homogeneous condensate is unstable. In the London approximation the free energy per unit volume of a periodic Abrikosov lattice with areal density $n_v$ is
\[
F/V=-|T|\,n_v+\frac12\,n_v\sum_{\mathbf R\neq 0}V_{\rm int}(|\mathbf R|),
\]
and minimization yields an optimal vortex density. For a triangular lattice the inter-vortex spacing is
\[
a(\mu_I,\mu_B)=\Bigl(\frac{2\Phi_0}{\sqrt3\,|T|}\Bigr)^{1/2},
\qquad
\Phi_0=2\pi/e,
\]
while a square lattice is slightly higher in energy by the standard factor $\simeq 1.06$. The triangular-array parameter is $\beta\approx1.1596$ [2403.07433].

These energetics lead to a simple low-temperature, zero-magnetic-field phase diagram. For $\mu_I<m_\pi$, the ground state is the trivial vacuum. For $\mu_I>m_\pi$ but $\mu_B<\mu_B^c(\mu_I)$, the ground state is the homogeneous charged-pion condensate, possibly interrupted by an Abrikosov lattice of pure $\pi^\pm$ vortices if an external magnetic field is present. Once $\mu_B$ crosses the critical line $\mu_B^c(\mu_I)$, it becomes energetically favorable to form a two-dimensional lattice of linked vortices, with each unit cell carrying one unit of baryon number [2509.20844].

## 6. Magnetic fields, rotation, and physical interpretation

One of the striking consequences of the baryonic vortex phase is that magnetic flux is intrinsic to the vortex lattice rather than externally imposed. In the 2024 analysis, when the tension becomes negative the system lowers its energy by creating vortices spontaneously without an external magnetic field. The magnetic field in the lattice phase is a superposition of single-vortex profiles,
\[
B(\mathbf x)=\sum_{m,n}B_{\rm v}\bigl(\mathbf x-m\mathbf a_1-n\mathbf a_2\bigr),
\]
with maximal magnitude near the core reaching
\[
B_{\max}\sim 0.1\,f_\pi^2 s^2\sim 10^{17}\,{\rm G}
\]
in the quoted numerics, while the spatial average
\[
\langle B\rangle=n_v\Phi_0
\]
may be of order $10^{13}$–$10^{15}\,{\rm G}$ for astrophysical densities. For $\mu_I\gtrsim m_\pi$ and $\mu_B$ in the range $0.74$–$0.94\,{\rm GeV}$, the ground state is described there as this vortex lattice rather than uniform nuclear matter or pion condensate [2403.07433].

A distinct extension concerns rotating nuclear matter. In a co-rotating frame with angular velocity $\Omega$ about the $z$-axis, the chiral action is formulated with the rotating metric
\[
g_{00}=1-\rho^2\Omega^2,\qquad
g_{0\phi}=g_{\phi0}=-\rho^2\Omega,\qquad
g_{\phi\phi}=-\rho^2,\qquad
g_{\rho\rho}=g_{zz}=-1.
\]
Within this framework two configurations are identified: a local vortex with charged-pion winding and a global vortex with neutral-pion winding, both carrying baryon number. The analysis emphasizes a point often treated as an objection to global vortices: in an infinite system their energy diverges logarithmically, but the finite-size constraint dictated by causality in a rotating frame regularizes the divergence physically, rendering the global vortex a viable excitation. The critical angular velocity is defined by
\[
\min_d T(d;\Omega_c)=0,
\]
and numerically one finds
\[
\Omega_c\sim O(0.01\,f_\pi \ldots 0.1\,f_\pi)
\]
for $f_\pi R\approx 5\ldots 7$ and $\mu\approx(3\ldots20)\,f_\pi$, corresponding to $\Omega_c\sim(10\ldots100)\,{\rm MeV}$ for $R\sim5\,{\rm fm}$. Once $\Omega>\Omega_c$, the vortex areal density obeys the Feynman relation
\[
n_v\approx \frac{2\Omega}{\kappa},
\qquad
\kappa=\frac{2\pi}{\mu_B},
\]
and the equilibrium lattice is triangular, with
\[
A_{\rm cell}=\frac{\sqrt3}{2}a^2=\frac{1}{n_v},
\qquad
a=\sqrt{\frac{2}{\sqrt3\,n_v}}.
\]
For heavy-ion-collision parameters $\Omega\sim50\,{\rm MeV}$ and $\mu_B\sim300\,{\rm MeV}$, this gives $a\sim5\ldots10\,{\rm fm}$ [2603.29325].

The present literature therefore places the baryonic vortex lattice at the intersection of several low-energy QCD mechanisms: pion condensation at finite isospin density, anomaly-induced baryon number from the Wess–Zumino–Witten or Goldstone–Wilczek current, Skyrmion topology, Abrikosov lattice energetics, and, in rotating systems, causality-limited finite-size regularization of global vortices. The astrophysical interpretation advanced in these works is that such a phase may be relevant to neutron-star interiors, neutron-star crusts, or phases of heavy-ion collisions where both isospin and baryon asymmetries are large; the rotational analysis further suggests that the previously overlooked global vortex can play a significant role in the topological structure of rotating dense QCD matter [2509.20844] [2603.29325].

Source: https://www.emergentmind.com/topics/baryonic-vortex-lattice