Superfluid Pion Multi-Vortices

• Superfluid pion multi-vortices are quantized defects in a rotating pion condensate where vortex winding reveals key aspects of chiral transport in QCD.
• An effective SU(2) chiral Lagrangian with a radial profile ansatz rigorously defines vortex core structure and circulation quantization.
• Fermionic zero modes confined to vortex lines yield enhanced chiral currents that diverge from standard hydrodynamic anomaly predictions.

Superfluid pion multi-vortices are quantized topological defects arising in a rotating pion condensate in the presence of finite isospin chemical potential. These structures manifest as lattices or arrangements of vortex lines, each associated with nontrivial winding of the superfluid order parameter and supporting localized fermionic zero modes propagating at the speed of light. The study of multi-vortex configurations is central to understanding chiral transport phenomena in quantum chromodynamics (QCD) at low temperature and high isospin density, particularly the chiral vortical effect in superfluid media (Kirilin et al., 2012).

1. Effective Theory and Superfluid Pion Condensation

The low-energy sector of two-flavor QCD at chemical potential $\mu_I$ and vanishing temperature is effectively captured by the SU(2) chiral Lagrangian. The pertinent Lagrangian is

$$\mathcal{L} = \frac{f_\pi^2}{4} \operatorname{Tr}[D_\mu U (D^\mu U)^\dagger]$$

with the covariant derivatives $D_0 U = \partial_0 U - \frac{\mu_I}{2} [\tau_3,U]$ and $D_i U = \partial_i U$. Superfluidity stems from spontaneous breaking of $U(1)_{L+R}$ (isospin-$I_3$ rotations) via condensation of $\pi^\pm$ mesons. The order parameter is encapsulated in the nonlinear field $U(x) = e^{i \tau_3 \phi(x)}$, where $\phi$ is the Goldstone mode corresponding to the superfluid phase. The ground state at rest obeys $\phi = \mu_I t$, reflecting the Josephson relation $\partial_0\phi = \mu_I$.

2. Vortex Solutions: Ansatz, Quantization, and Core Structure

A superfluid vortex of winding number $n$ is characterized by a multi-valued phase $\phi$. In the London approximation (constant amplitude), the phase takes the form $\phi(t, r, \varphi) = \mu_I t + n\varphi$, with $(r,\varphi)$ denoting plane-polar coordinates orthogonal to the vortex. A more refined ansatz incorporates a radial profile $\alpha(r)$: $$U(r, \varphi, t) = \cos\alpha(r) + i \sin\alpha(r)\left[\tau_1\cos(n\varphi) + \tau_2\sin(n\varphi)\right] \; e^{i\tau_3 \mu_I t}$$ with boundary conditions $\alpha(0) = 0$ (uncondensed vortex core) and $\alpha(r\to\infty)=\pi/2$ (fully condensed far from the core). The profile satisfies

$$\alpha''(r) + \frac{1}{r}\alpha'(r) - \frac{n^2}{r^2}\sin\alpha\cos\alpha - \frac{\mu_I^2}{4}\sin(2\alpha) = 0.$$

The quantization of circulation follows from the integral of $\partial_i\phi$ along a closed path: $\oint \partial_i \phi\,dx^i = 2\pi n.$ The vortex core radius $a$ is set by the inverse radial (Higgs) mode mass, $a \sim m^{-1} \sim \mu_I^{-1}$.

3. Fermionic Zero Modes: Dirac Equation and Index Theorem

A massless Dirac fermion $\Psi$ couples to the superfluid phase as an axial gauge field, with Lagrangian

$$\mathcal{L}_\psi = \bar{\Psi} i\gamma^\mu (\partial_\mu + i\partial_\mu\phi) \Psi.$$

In the presence of a vortex, $\partial_i\phi$ acts as an effective Aharonov–Bohm flux. Decomposing into Weyl components and Fourier transforming along the vortex axis ($z$-direction), the right-handed sector's effective Hamiltonian is $H_R = p_3 \sigma_3 + H_\perp$, where $$H_\perp \equiv (-i\partial_a - \partial_a\phi)\sigma_a$$ ($a=1,2$). The anticommutator $\{\sigma_3, H_\perp\}=0$ ensures paired nonzero eigenvalues $\pm\lambda$, but zero modes persist and are determined via the index theorem: $$\mathrm{Index}\,H_\perp = N_+ - N_- = \frac{1}{2\pi} \oint \partial_i\phi\,dx^i = n,$$ where $N_+$ ($N_-$) counts zero modes with chirality $+1$ ($-1$) under $\sigma_3$. For $n=1$, one finds a single normalizable zero mode localized on the vortex.

4. Chiral Current and its Lattice Extension

Zero modes correspond to 1+1D chiral fermions propagating along the vortex at light speed. Their dispersion is $\epsilon(p_3) = +p_3$ for $N_+$ zero modes and $-p_3$ for $N_-$. States with $|p_3|\le \mu_I$ are occupied. The axial current carried by a single vortex is

$$J^3_5 = \int d^2x\, \langle \bar{\Psi} \gamma^3\gamma^5 \Psi \rangle = (N_+ - N_-)\int_{|p_3| \le \mu_I} \frac{dp_3}{2\pi} = n\, \frac{\mu_I}{\pi}.$$

In contrast, standard hydrodynamic anomaly-based arguments predict a chiral current density in uniform rotation

$$j^3_5 = \frac{\mu_I^2}{2\pi^2} \Omega \qquad \text{(per unit volume)},$$

where $\Omega$ is the angular velocity.

When multiple ($N$) vortices are present, typically arranged in a triangular lattice via global rotation, the area density is $\rho_v = N/V \simeq \mu_I \Omega/\pi$. The total axial current and current density generalize to

$$J^3_5(\mathrm{total}) = N\, \frac{\mu_I}{\pi} \quad \Longrightarrow \quad j^3_5 = \frac{\mu_I^2 \Omega}{\pi^2}.$$

5. Hydrodynamic Anomaly Versus Zero-Mode Mechanism: Factor-of-Two Discrepancy

A comparison of the vortex-based and continuum approaches reveals a numerical discrepancy: the zero-mode calculation yields an axial current twice that of the standard hydrodynamic result. The origin of the factor-two difference is traced to the underlying anomaly mechanism. In the superfluid vortex case, both vertices of the relevant diagram correspond to identical insertions of $\mu_I u^\mu \bar{\Psi}\gamma_\mu\Psi$. In contrast, the magnetic case involves distinct vertices for the chemical potential and the external gauge field. Physically, the zero modes propagate at the speed of light along the vortex—unlocked from the local superfluid four-velocity $u^\mu$. Hydrodynamic treatments that assume all chiral charge is carried by a single velocity field thus undercount this contribution by a factor of $1/2$ (Kirilin et al., 2012).

6. Microscopic Realization and Physical Implications

In a rotating superfluid pion condensate, the chiral vortical effect is realized via 1+1D fermionic zero modes confined to vortex lines. The chiral current scales linearly with the total vortex winding and, in a lattice configuration, reflects an enhancement over conventional anomaly-induced currents. Each unit-winding vortex supports one chiral zero mode; the total current is obtained by summing over all vortices. The result is

$j^3_5 = \frac{\mu_I^2\Omega}{\pi^2}$

for a uniformly rotating sample, explicitly demonstrating the departure from hydrodynamic expectations.

Significance arises both in confirming nontrivial topological transport in QCD superfluids and in the elucidation of anomaly-related enhancements due to light-like zero-mode propagation. This robust mechanism illustrates the limitations of naive hydrodynamic approaches and underscores the necessity of microscopic treatments when addressing quantized vortices in quantum fluids (Kirilin et al., 2012).

7. Summary Table: Key Quantities in Superfluid Pion Multi-vortices

Quantity	Formula	Context
Vortex circulation	$2\pi n$	$n$ is winding number
Zero modes per vortex	$n$	Index theorem (chirality $\pm 1$)
Axial current (per vortex)	$J^3_5 = n\, \frac{\mu_I}{\pi}$	Occupied up to $|p_3| = \mu_I$
Vortex lattice density	$\rho_v = \frac{\mu_I \Omega}{\pi}$	Triangular lattice, uniform rotation
Total chiral current	$J^3_5(\mathrm{total}) = N\, \frac{\mu_I}{\pi}$	$N$ unit-winding vortices
Current density	$j^3_5 = \frac{\mu_I^2\Omega}{\pi^2}$	Per unit volume, lattice configuration

The collective axial current in superfluid pion multi-vortex systems manifests the nonlocal, light-like transport properties intrinsic to fermionic zero modes—a distinctive phenomenon beyond the reach of local hydrodynamic anomaly calculations (Kirilin et al., 2012).