Axially Symmetric Fast Flavor Oscillations

• Axially symmetric fast flavor system is defined for neutrino ensembles using axially averaged density matrices and is governed by rapid flavor instabilities.
• The system’s instability criteria hinge on ELN crossings and Nyquist winding, allowing analytic determination of conditions for fast flavor conversion.
• Nonlinear dynamics map onto a gyroscopic pendulum and are explored via moment and PIC methods, with implications for supernovae and merger environments.

An axially symmetric fast flavor system refers to a many-body neutrino ensemble exhibiting collective flavor oscillations driven by neutrino-neutrino forward scattering, under the assumption of exact symmetry about a given axis in momentum space. In such systems, the angular distributions of neutrinos and antineutrinos, projected along the symmetry axis, are the central objects, and the system’s dynamics is governed by the so-called fast flavor instabilities. These instabilities, even in the absence of vacuum mixing, can induce rapid (nanosecond-scale) and potentially large-amplitude flavor conversions, especially when the electron lepton number (ELN) angular distribution features a crossing. The reduced system respects axial symmetry, thereby greatly facilitating analytic treatment while capturing the essential instability mechanisms, and serves as the principal context for investigating the homogeneous "fast flavor pendulum" and its nonlinear analogs.

1. Formulation and Equations of Motion

Axially symmetric fast flavor systems are defined by energy- and azimuth-averaged density matrices $\rho_v(z,t)$ for neutrinos and $\bar\rho_v(z,t)$ for antineutrinos, with $v\equiv \cos\theta$ the projection of the velocity on the axis of symmetry. In the fast regime (vacuum effects negligible), the evolution is described by the Liouville–von Neumann equations linearized around a flavor-diagonal state. The only nontrivial dynamics arises from the off-diagonal (flavor coherence) modes, characterized by the "flavor field" $\psi_v$. The linearized equation governing these modes under axial symmetry is

$$(\partial_t + v\partial_z)\psi_v(z,t) = i\mu\left[\psi_v(G_0 - v G_1) - G_v(\psi_0 - v\psi_1)\right],$$

where $\mu = \sqrt{2}(n_\nu + n_{\bar\nu})$ is the neutrino self-interaction potential, $G_v$ is the fixed lepton-number spectrum, and $G_n = \int_{-1}^{+1}dv\, G_v v^n$, $\psi_n = \int_{-1}^{+1}dv\, \psi_v v^n$ are angular moments (Fiorillo et al., 2024).

In the fully nonlinear regime, particularly in two-flavor models, the evolution can be cast into equations for polarization vectors or flavor-spin vectors (Bloch representation), which capture the geometric structure of the collective motion and allow exact conservation statements about the system's dynamics (Padilla-Gay et al., 2021, Fiorillo et al., 26 Jan 2026).

2. Instability Criteria and Dispersion Relations

The presence or absence of fast flavor instabilities in these systems is ultimately determined by the structure of the ELN distribution $G_v$. The linear normal-mode analysis seeks solutions of the form $\psi_v \propto e^{-i\Omega t + iKz}$, leading to the shifted frequency and wavenumber variables:

$\omega = \Omega + \mu G_0,\quad k = K + \mu G_1.$

The resulting dispersion relation in the axially symmetric case is

$$(I_0 - 1)(I_2 + 1) - I_1^2 = 0,\qquad I_n(\omega, k) = \mu \int_{-1}^{+1}\frac{G_v v^n}{\omega - k v + i\epsilon}dv.$$

Physical (causal) solutions are selected by the Landau prescription, which places the integration pole just below the real axis, automatically accounting for Landau damping and identifying modes contributing to late-time evolution (Fiorillo et al., 2024, Fiorillo et al., 26 Jan 2026).

The strictly homogeneous ($k=0$) mode, referred to as the fast flavor pendulum, admits a specialized form:

$$I_1(\omega, k=G_1) = \mu \int_{-1}^{+1} \frac{G_v v}{\omega - G_1 v + i\epsilon}dv = 0.$$

This equation determines the collective frequency for the pendulum mode.

3. Angular Crossing, Nyquist Criterion, and Growth Conditions

A necessary condition for instability is the existence of a crossing in the ELN distribution—a value $v_c$ such that $G_{v_c}=0$ within $(−1,1)$. Such a crossing guarantees two critical points in the subluminal domain $|u|=|ω/k|<1$ where modes can transition from damped to unstable. However, the existence of an unstable solution at a given $k$ requires a more stringent Nyquist criterion: the winding number $W$ of the complex argument of $$\Phi(u) = \int_{-1}^{+1}\frac{G_v v}{u-v + i\epsilon}dv$$ as $u$ traverses the real axis (skirting branch cuts at $u = \pm 1$). If $W - N_s/2 > 0$, where $N_s$ counts real superluminal roots $|u|>1$, then the system is unstable. In particular, instability for the homogeneous mode requires that no superluminal critical points "soak up" the Nyquist winding; otherwise, only inhomogeneous modes with $k \neq 0$ are unstable (Fiorillo et al., 2024).

The table below summarizes instability analysis steps:

Criterion	Physical Meaning	Diagnostic
ELN crossing	Necessary for instability	$G_v=0$ for some $v$ in $(-1,1)$
Nyquist winding $W$	Sufficient + counts unstable roots	$W = (1/2\pi)\Delta\arg\Phi(u)$
Superluminal zeros	Subtract possible N_s/2 for stability	Inspect $\Phi(u)$ for $|u|>1$

4. Nonlinear Dynamics and Pendulum Mapping

Nonlinear dynamics of the strictly homogeneous axially symmetric mode are integrable and can be mapped onto a gyroscopic pendulum in an abstract "flavor-spin" space. The full nonlinear equations admit several constants of motion, including the energy and "pendulum spin." The latter, defined as $S = 2λσ = 2\Omega/\mu$, is directly related to $\Omega$, the real part of the eigenfrequency of the linear normal mode, and encapsulates the depth of flavor conversion (Padilla-Gay et al., 2021).

In the regime $\sigma < 1$ (with $\sigma = \Omega/\sqrt{\Omega^2+\Gamma^2}$), the motion of the pendulum is periodic, with the maximal flavor conversion fixed by

$\cos\theta_{\rm max} = -1 + 2\sigma^2.$

This explicit connection allows prediction of the degree of flavor conversion solely from linear analysis without integrating the nonlinear system (Padilla-Gay et al., 2021).

5. Weak and Strong Instability Regimes: Resonant Behavior and Saturation

In the weak instability regime, such as when the ELN crossing is shallow, the growth rate $\gamma = \operatorname{Im} \Omega$ is small and the single-mode nonlinear solution becomes localized in velocity space: only the subset of neutrinos whose velocities are resonant with the wave ($v_{\mathrm{res}} = \omega_R/k$) are significantly converted. Their dynamics reduce to a single-pendulum Hamiltonian system:

$$\ddot{\Phi}_u = \gamma^2\sin\Phi_u, \qquad \mathcal{E} = \frac{1}{2}(\dot\Phi_u)^2 + \gamma^2\cos\Phi_u,$$

admitting full flavor swap and return cycles with characteristic period $2\pi/\gamma$ and maximal field amplitude $2\gamma/(1-Uu)$ (Fiorillo et al., 26 Jan 2026).

For broad crossings (strong instability, $\gamma \sim |\Omega_R|$), many angular modes interact, pendular regularity is lost unless the system reduces to a two-beam scenario, and the system exhibits complex multi-mode dynamics and rapid decoherence.

6. Numerical Approaches: Moment and PIC Methods

Direct simulation of the fast flavor system can employ angular-integrated moment schemes, typically closing at the level of the second (pressure) moment using maximum entropy (Minerbo) closures. These methods reproduce growth rates within $\sim$10–40% and capture the qualitative evolution, though with limitations for highly anisotropic or nonthermal ELN distributions (Grohs et al., 2023). A particle-in-cell (PIC) approach offers higher fidelity, tracking ensemble flavor evolution in phase space at substantial computational cost and resolving saturation and decoherence processes (Richers et al., 2021).

PIC simulations confirm:

• Quantitative agreement of linear growth rates for axially symmetric ($m=0$) and symmetry-breaking ($m=1$) modes with analytic predictions.
• Equipartition of flavors post-saturation for sufficiently strong crossings.
• The minimum between total $\nu$ and $\bar\nu$ sets the limiting change in flavor, with weak crossings resulting in suppressed conversion even after saturation.

7. Physical Role in Supernovae and Merger Environments

Fast flavor instabilities, mediated by the axially symmetric system, are ubiquitous in dense astrophysical environments such as core-collapse supernovae and neutron star mergers. The symmetry assumptions are locally justified by the geometry of emission regions. Regions exhibiting an ELN crossing are susceptible to rapid flavor equilibration, affecting nucleosynthesis, neutrino detection signatures, and overall dynamics. However, even when the homogeneous mode ($k=0$) is stable, inhomogeneous, higher-$k$ modes may drive flavor evolution. The rigorous distinction between the pendulum’s special symmetric, regular behavior and the generic, fragile character of most collective modes underlines the importance of symmetry assumptions and the presence of inhomogeneities, noise, or collisions that may break regularity and enhance flavor decoherence (Fiorillo et al., 2024, Grohs et al., 2023).

• Fiorillo, Goimil-García & Raffelt, "Fast Flavor Pendulum: Instability Condition" (Fiorillo et al., 2024)
• Capozzi, Dasgupta & Mirizzi, "Neutrino Flavor Pendulum Reloaded: The Case of Fast Pairwise Conversion" (Padilla-Gay et al., 2021)
• Grohs et al., "Two-Moment Neutrino Flavor Transformation with applications to the Fast Flavor Instability in Neutron Star Mergers" (Grohs et al., 2023)
• Richers et al., "Particle-in-cell Simulation of the Neutrino Fast Flavor Instability" (Richers et al., 2021)
• Duan, Raffelt & Tamborra, "Single-wave solutions of the neutrino fast flavor system. Part II. Weak instabilities and their resonant behavior" (Fiorillo et al., 26 Jan 2026)