Chiral Plasma Instability Overview

• Chiral plasma instability is a helical gauge-field instability in parity-violating plasmas driven by a chiral imbalance and axial anomaly.
• It is analyzed using Berry-curvature kinetic theory and magnetohydrodynamic models to characterize unstable modes and determine growth rates.
• Its nonlinear evolution leads to turbulence, inverse cascade, and modified transport properties that impact astrophysical and cosmological magnetic fields.

Chiral plasma instability (CPI) is a long-wavelength helical gauge-field instability of a parity-violating plasma with nonzero chirality imbalance. In the relevant literature, the imbalance is parameterized either by $\mu_5 \equiv \mu_R-\mu_L$ or by $\mu_5 \equiv (\mu_R-\mu_L)/2$, depending on convention. The defining mechanism is anomaly-driven: a chiral asymmetry induces a parity-odd current or polarization term, one circular polarization of the electromagnetic, hypercharge, or non-Abelian gauge field becomes unstable, and fermionic chirality is converted into magnetic helicity or Chern–Simons number. CPI is therefore simultaneously a collective-mode instability, a helicity-transfer process, and a transport phenomenon with nonlinear consequences for turbulence, viscosity, heating, and magnetic-field generation in astrophysical and cosmological plasmas (Akamatsu et al., 2013, Akamatsu et al., 2014, Mace et al., 2020).

1. Anomalous origin and helicity transfer

The microscopic origin of CPI is the axial anomaly. In chiral gauge plasmas, axial charge is not conserved independently of gauge topology; instead, fermionic chirality can be exchanged continuously with gauge-field helicity. In QED language this is encoded in the global relation

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$

while in the chiral limit of the lattice QED simulations the conserved quantity is

$N_5(t)+2N_fN_h(t),$

with $N_h$ equal to magnetic helicity in QED (Ohnishi et al., 2014, Mace et al., 2020).

This conservation structure explains the energetic tendency behind the instability. A plasma with $\mu_5\neq 0$ can lower its free energy by transferring chirality stored in fermions into helical gauge-field structure. In kinetic descriptions this appears as a parity-odd component of the induced current or polarization tensor; in magnetohydrodynamic descriptions it appears as a current parallel to $\mathbf B$, i.e. the chiral magnetic effect. CPI is the dynamical channel through which that anomalous coupling reorganizes the plasma into helical magnetic configurations (Akamatsu et al., 2013).

A recurring point in the literature is that chirality imbalance is the driver, not magnetic helicity itself. The instability vanishes when the parity-odd coefficient vanishes. In Berry-curvature kinetic theory, for example, the relevant coefficient is

$C_E=\frac{e^2\mu_5}{4\pi^2},$

so $\mu_5=0$ removes the parity-odd response and the unstable branch disappears (Kumar et al., 2016).

2. Linear theory and unstable mode structure

The canonical linear analysis uses Berry-curvature-modified kinetic theory. In the isotropic collisionless case, the spatial self-energy is decomposed as

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$

and the transverse modes obey

$\mu_5 \equiv (\mu_R-\mu_L)/2$0

In the quasistatic regime $\mu_5 \equiv (\mu_R-\mu_L)/2$1, the parity-even sector gives the familiar Landau-damped response while the parity-odd term splits the two circular polarizations. One branch admits a purely imaginary frequency, $\mu_5 \equiv (\mu_R-\mu_L)/2$2, and grows exponentially when $\mu_5 \equiv (\mu_R-\mu_L)/2$3 (Akamatsu et al., 2013, Kumar et al., 2016).

For the Abelian isotropic instability, the growth rate may be written as

$\mu_5 \equiv (\mu_R-\mu_L)/2$4

or, in the notation of the turbulent-transport analysis,

$\mu_5 \equiv (\mu_R-\mu_L)/2$5

The unstable band is therefore long-wavelength, and the fastest growth occurs at

$\mu_5 \equiv (\mu_R-\mu_L)/2$6

This is not a Weibel instability of an anisotropic distribution; it exists already in an isotropic plasma provided the chiral imbalance is nonzero (Akamatsu et al., 2013, Kumar et al., 2016).

The same physics can be written directly at the MHD level. For a relativistic plasma with magnetic diffusivity $\mu_5 \equiv (\mu_R-\mu_L)/2$7,

$\mu_5 \equiv (\mu_R-\mu_L)/2$8

and for circularly polarized Fourier modes

$\mu_5 \equiv (\mu_R-\mu_L)/2$9

Hence one helicity is unstable for

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$0

with maximal growth near $$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$1 and characteristic time

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$2

For $$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$3, the unstable modes are right-circularly polarized; for $$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$4, they are left-circularly polarized (Armstrong et al., 7 May 2026).

3. Nonlinear evolution, turbulence, and transport

Linear growth is only the first stage of CPI. In classical-statistical real-time lattice simulations of strongly coupled QED with $$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$5, a chirality imbalance in the fermion sector produces exponential growth of magnetic helicity in a narrow band around

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$6

The linear phase terminates at approximately

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$7

when the right-handed magnetic occupation reaches the nonperturbative regime,

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$8

The subsequent evolution is not simple damping: magnetic helicity is redistributed toward the infrared through a self-similar inverse cascade,

$$\frac{d}{dt}\left(N_5+\frac{\alpha}{\pi}{\cal H}\right)=0, \qquad {\cal H}=\int d{\bm x}\,{\bm A}\cdot{\bm B},$$9

with extracted exponents

$N_5(t)+2N_fN_h(t),$0

This is the direct numerical realization of chirality transfer into long-range helical magnetic fields (Mace et al., 2020).

The same nonlinear sector alters transport coefficients. In the turbulent-transport analysis, unstable chiral modes generate fluctuating electromagnetic fields that enhance scattering through a decorrelation rate. Resonance broadening replaces the sharp factor

$N_5(t)+2N_fN_h(t),$1

by

$N_5(t)+2N_fN_h(t),$2

with decorrelation time $N_5(t)+2N_fN_h(t),$3 determined self-consistently. Evaluated at the most unstable mode, the anomalous shear viscosity is estimated as

$N_5(t)+2N_fN_h(t),$4

In the regime $N_5(t)+2N_fN_h(t),$5, where $N_5(t)+2N_fN_h(t),$6, this implies $N_5(t)+2N_fN_h(t),$7, while in the low-temperature limit $N_5(t)+2N_fN_h(t),$8 (Kumar et al., 2016).

Energy redistribution in the nonlinear stage is also thermodynamic rather than purely magnetic. A recent analysis shows that corrected chiral MHD equations conserve energy locally and that the chiral sector does not transfer all of its free energy into the growing magnetic field. In the regime $N_5(t)+2N_fN_h(t),$9, the plasma is heated by

$N_h$0

and, in the monochromatic unstable range $N_h$1, the thermal-energy gain exceeds the magnetic-energy gain,

$N_h$2

Ohmic dissipation is therefore part of the endpoint of CPI, not merely a perturbation of it (Armstrong et al., 7 May 2026).

A further nonlinear qualification comes from fully backreacting classical-statistical field theory. There, chirality pumping in parallel $N_h$3 and $N_h$4 is self-limited because the electric field is screened by produced on-shell fermions, and for initially chirally imbalanced states the relation between inverse cascade and axial-charge decay is not always one-to-one. In particular, very clear inverse cascade can occur without appreciable axial-charge decay, which suggests that the simplest anomalous-Maxwell description is incomplete at the level of full fermion backreaction (Buividovich et al., 2015).

4. Anisotropy, non-Abelian plasmas, and related instability families

CPI does not remain isolated in more general plasma states. In anisotropic chiral plasmas, Berry-curvature kinetic theory yields a mixed chiral-Weibel spectrum. For weak anisotropy, the chiral and Weibel growth rates become comparable near

$N_h$5

and for QED with $N_h$6, $N_h$7. At $N_h$8, Weibel modes dominate for $N_h$9 when $\mu_5\neq 0$0, whereas sufficiently large $\mu_5\neq 0$1 or special propagation angles can restore chiral-mode dominance. Parallel propagation enhances the instability; perpendicular propagation suppresses it (Kumar et al., 2016, Kumar et al., 2014).

In high-temperature non-Abelian plasmas, the instability shifts from magnetohydrodynamics to the magnetic scale

$\mu_5\neq 0$2

Because this lies below the mean-free-path scale, conventional anomalous hydrodynamics is not the appropriate effective theory. The relevant description is a stochastic chiral Langevin theory,

$\mu_5\neq 0$3

with the corresponding instability time

$\mu_5\neq 0$4

Here the anomaly transfers chiral charge into non-Abelian Chern–Simons number rather than ordinary magnetic helicity (Akamatsu et al., 2014).

The same anomaly-driven logic also generates related instability classes beyond standard CME-driven CPI. In the chiral magnetovortical instability, the chiral vortical effect induces a coupled magnetic-vortical mode in the presence of a background magnetic field; for the pure CVE-driven case, the linear instability appears when

$\mu_5\neq 0$5

and the chirality transfer proceeds primarily into cross helicity rather than magnetic helicity (Wang et al., 2023). In rotating chiral MHD, rigid rotation splits the Alfvén mode into fast and slow magneto-Coriolis branches, and the slow branch can be destabilized by arbitrarily weak chiral vortical effect in an appropriate kinematic window, making rotation a catalyst of the chiral magnetovortical instability (Wang et al., 11 Feb 2026).

5. Astrophysical and cosmological realizations

Compact-star applications motivated much of the early CPI literature. In proto-neutron-star cores, the parity-violating capture process

$\mu_5\neq 0$6

preferentially removes left-handed electrons and can generate $\mu_5\neq 0$7. Using $\mu_5\neq 0$8 with $\mu_5\neq 0$9, the characteristic instability time is estimated as

$\mathbf B$0

and full conversion of the chiral free energy yields an order-of-magnitude core field

$\mathbf B$1

Because the same process generates magnetic helicity,

$\mathbf B$2

the amplified field is not merely strong but topologically stabilized in the sense emphasized in magnetar models (Ohnishi et al., 2014).

A distinct supernova mechanism replaces electron chiral imbalance by nonequilibrium left-handed neutrino backreaction. In that case the matter sector acquires a current

$\mathbf B$3

even with no electron chiral imbalance, with an estimated upper bound $\mathbf B$4 under typical supernova conditions. The corresponding CPI dispersion relation,

$\mathbf B$5

again yields $\mathbf B$6, and conversion of the effective chiral free energy gives

$\mathbf B$7

Three-dimensional chiral MHD simulations then show magnetic-helicity growth and inverse cascade that is mainly CPI-driven when

$\mathbf B$8

(Matsumoto et al., 2022).

Magnetospheric applications are more selective. In gap regions of pulsar and black-hole magnetospheres, $\mathbf B$9 generates chiral charge via

$C_E=\frac{e^2\mu_5}{4\pi^2},$0

For supermassive black holes in AGN, the resulting $C_E=\frac{e^2\mu_5}{4\pi^2},$1 is negligibly small. For magnetars, however, the gap-induced chiral asymmetry can be substantial; the fastest-growing mode has

$C_E=\frac{e^2\mu_5}{4\pi^2},$2

with an estimated instability time

$C_E=\frac{e^2\mu_5}{4\pi^2},$3

for $C_E=\frac{e^2\mu_5}{4\pi^2},$4 MeV. The unstable helical modes imply transient circularly polarized emission over

$C_E=\frac{e^2\mu_5}{4\pi^2},$5

spanning radio to near-infrared frequencies (Gorbar et al., 2021).

Cosmological CPI is formulated in the unbroken electroweak phase in terms of hypermagnetic fields. There the control parameter is the hypercharge-weighted asymmetry $C_E=\frac{e^2\mu_5}{4\pi^2},$6, with unstable modes satisfying

$C_E=\frac{e^2\mu_5}{4\pi^2},$7

Three-dimensional simulations find maximally helical hypermagnetic growth, inverse cascade with

$C_E=\frac{e^2\mu_5}{4\pi^2},$8

and a relic gravitational-wave signal with

$C_E=\frac{e^2\mu_5}{4\pi^2},$9

An earlier gravitational-wave estimate for $\mu_5=0$0 GeV and $\mu_5=0$1 gave a peak around $\mu_5=0$2 Hz, illustrating the model dependence of the redshifted signal (Brandenburg et al., 2023, Anand et al., 2018).

6. Energy accounting, suppression mechanisms, and viability

A persistent controversy concerns whether CPI can generate large astrophysical fields under realistic microphysics. An energy-conservation analysis of a homogeneous neutral plasma shows that the effective potential for a maximally helical mode has negative curvature at the origin only when there is a chiral asymmetry in electron Fermi energy, not merely in number density. The potential rises at large field amplitude, so there is no unlimited runaway, and the true ground state has zero magnetic helicity. In that framework, finite electron mass allows an excited helical state to relax quickly, leading to the conclusion that weak-interaction-triggered CPI is not a viable mechanism for stellar magnetic fields except possibly when the dynamics is driven far from equilibrium (Kaplan et al., 2016).

Far-from-equilibrium driving is precisely where some later work reopens the question. Abrupt density and temperature perturbations generated by shocks in core-collapse supernovae and neutron-star mergers can sustain a chiral imbalance despite electron-mass damping. The practical criterion is

$\mu_5=0$3

A representative strong-shock case gives

$\mu_5=0$4

corresponding to amplification by

$\mu_5=0$5

In an already magnetized medium, the same shock-induced CME can produce Joule heating

$\mu_5=0$6

and in the strong example the heating exceeds the shock thermal deposition,

$\mu_5=0$7

This suggests that CPI and CME heating may coexist as microphysical shock effects rather than as equilibrium dynamos (Harris et al., 24 Feb 2026).

Cosmological viability is likewise sensitive to source and washout terms. A sourced early-universe model introduces

$\mu_5=0$8

which allows CPI to operate below the nominal $\mu_5=0$9 TeV chirality-washout threshold. In that setting the generalized total-helicity production law is

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$0

and the saturated magnetic helicity is estimated as

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$1

This suggests that continuous chirality injection changes CPI from one-time relaxation into a sustained magnetogenesis mechanism (Gurgenidze et al., 9 Dec 2025).

A more pessimistic analysis of finite-time pumping argues that realistic source durations strongly suppress the dynamo. If $$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$2 is the ideal instantaneous CPI rate and

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$3

then the actual growth rate satisfies

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$4

Chiral flipping hinders the dynamo when

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$5

and completely suppresses it when

$$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$6

Because realistic $$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$7 often gives $$\Pi^{ij}(K)=\Pi_T P_T^{ij}+\Pi_L P_L^{ij}+\Pi_A P_A^{ij}, \qquad P_A^{ij}=i\epsilon^{ijk}\hat{k}^k,$$8, the chiral dynamo becomes highly vulnerable to suppression. The result is strongly unfavorable in protoneutron stars and only marginally avoidable near the electroweak transition (Skoutnev et al., 8 Mar 2026).

Taken together, these results define the modern status of CPI. Its linear existence is well established; its helicity-transfer mechanism is fixed by the anomaly; and its nonlinear evolution can generate inverse cascade, anomalous transport, and plasma heating. Its quantitative efficacy, however, depends sensitively on backreaction, conductivity, chirality-flip rates, source duration, anisotropy, and departures from equilibrium. CPI is therefore best understood not as a universal magnetic dynamo, but as a family of anomaly-driven instabilities whose realization is highly regime dependent (Buividovich et al., 2015, Armstrong et al., 7 May 2026).