Flavor Pendulum in Neutrino Physics

• Flavor pendulum is a mechanical analogy that maps collective neutrino oscillations onto the motion of a pendulum or gyroscopic top, clarifying complex flavor dynamics.
• It captures key mechanisms such as energy exchange, instability growth, and nonlinear saturation in dense neutrino gases found in supernovae and mergers.
• The model integrates both Hamiltonian and collisional frameworks to provide analytic predictions for flavor conversion and coherence decay.

In modern neutrino physics, the term "flavor pendulum" refers to the mapping of nonlinear collective oscillations in dense neutrino gases onto the equations of motion of a mechanical pendulum or gyroscopic top. This analogy encompasses several regimes: the "slow" pendulum driven by vacuum oscillations, the "fast" pendulum driven solely by neutrino-neutrino interactions in the massless (fast) limit, and the more recent "collisional flavor pendulum" where decoherence is also active. The flavor pendulum provides analytic and conceptual closure for understanding nonlinear flavor evolution, saturation regimes, and the impact of environmental perturbations in supernovae, neutron-star mergers, and cosmology.

1. Mechanical Mapping and Equations of Motion

The core of the flavor pendulum analogy is the representation of collective flavor evolution as spin-precession in SU(2) flavor space. In the fast-flavor, axisymmetric case, each momentum or angular mode $v = \cos\theta$ is associated with a polarization vector $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$, where $P^z_v = D_v$ (diagonal lepton-number difference) and $P^x_v + iP^y_v = \Psi_v$ (coherence).

The single-wave equations (for a spatially homogeneous or monochromatic spatial mode) are (Fiorillo et al., 21 Jan 2026, Fiorillo et al., 26 Jan 2026): $$\dot D_v = \frac{i}{2}[(\Psi_0 - v\Psi_1)\Psi_v^* - (\Psi_0^* - v\Psi_1^*)\Psi_v],$$

$$\dot \Psi_v = i [ (D_0 - v D_1 - vK)\Psi_v - (\Psi_0 - v\Psi_1) D_v ].$$

Integrated over $v$, these yield conservation of total lepton number and oscillatory dynamics in $\Psi_0$ (collective coherence).

Each $\mathbf{P}_v$ precesses about a Hamiltonian field $$\mathbf{H}_v = (\text{Re} (\Psi_0 - v\Psi_1), -\text{Im} (\Psi_0 - v\Psi_1), D_0 - vD_1 - vK)$$, with the full set of EoMs taking the flavor–isospin cross-product form,

$\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$0

These equations exhibit a global spin–pendulum analogy, where the system's angular degrees of freedom and collective invariants mirror those of a gyroscopic top.

2. Dispersion Relations, Instabilities, and the Pendulum Reduction

The stability and nonlinear behavior of neutrino ensembles are set by the eigenstructure of the linearized dispersion relation. In the fast-flavor regime, the commonly used lepton-number distribution $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$1 (with a "crossing" at some $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$2) is essential.

The generic unstable mode satisfies (Padilla-Gay et al., 2021, Fiorillo et al., 2024): $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$3 where $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$4 is the ELN spectrum, $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$5 the wave number, and $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$6 the eigenfrequency (complex; its imaginary part is the instability growth rate $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$7).

For shallow crossings and weak instability ($\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$8), the nonlinear development is dominated by phase-locked resonant neutrinos satisfying $\mathbf{P}_v(t) = (P^x_v, P^y_v, P^z_v)$9. Their collective evolution isolates a two-variable subsystem,

$P^z_v = D_v$0

where $P^z_v = D_v$1, $P^z_v = D_v$2. The elimination of $P^z_v = D_v$3 yields a canonical pendulum equation,

$P^z_v = D_v$4

This describes periodic full flavor reversal and back-transfer of power to the unstable wave—an explicit nonlinear analog of Landau resonance in plasma physics (Fiorillo et al., 26 Jan 2026).

3. Energy, Integrability, and Taxonomy of Pendulum Solutions

In mechanical terms, the flavor pendulum conserves "energy" $P^z_v = D_v$5 and oscillates in the potential $P^z_v = D_v$6 For an inverted initial state, the analytic "soliton" solution is

$P^z_v = D_v$7

with small-amplitude oscillation frequency $P^z_v = D_v$8.

Integrability is critically dependent on spatial and angular homogeneity. In single-angle cases (slow or fast, homogeneous), Gaudin invariants guarantee complete integrability and coherent multi-mode pendula (Fiorillo et al., 21 Jan 2026). In the single-wave case, no Gaudin invariants exist, restricting exact pendulum solutions to two-beam systems; there is no extension to continuous angular spectra except for weak instabilities. This distinguishes “fast” homogeneous, integrable systems from “single-wave” inhomogeneous, non-integrable regimes.

Taxonomy of solutions:

Regime	Pendulum Mapping	Integrable Multi-mode	Nonlinear Extension
Slow/homogeneous	Yes	Yes (Gaudin)	Solitons, waves
Fast/homogeneous	Yes	Yes (Gaudin)	Solitons, coprecession
Single-wave/inhomogeneous	Yes (2-beam)	Only 2-beam	Approx. pendulum for weak instability

4. Collisional, Damped, and Asymmetric Pendulum Dynamics

In supernovae and mergers, collisional damping impacts flavor coherence. The mean-field kinetic EOMs (Padilla-Gay et al., 2022, Johns et al., 2023) with a common damping rate $P^z_v = D_v$9 (and, possibly, asymmetric $P^x_v + iP^y_v = \Psi_v$0) yield

$P^x_v + iP^y_v = \Psi_v$1

driving decay of transverse coherence and, hence, pendulum "length" over time.

Numerical and analytical study reveals that the steady-state flavor conversion fraction obeys an empirical universal law: $P^x_v + iP^y_v = \Psi_v$2 with $P^x_v + iP^y_v = \Psi_v$3 from linear analysis. This final value is independent of angular details, seed amplitude, and even $P^x_v + iP^y_v = \Psi_v$4 (for $P^x_v + iP^y_v = \Psi_v$5).

When $P^x_v + iP^y_v = \Psi_v$6, new "collisional instabilities" can arise even without ELN crossings; the instability criterion reads $P^x_v + iP^y_v = \Psi_v$7, permitting flavor conversion in ELN-stable systems (Padilla-Gay et al., 2022, Johns et al., 2023).

5. Practical Implications, Diagnostic Tools, and Physical Interpretation

The flavor pendulum paradigm provides both theoretical and practical benefits:

• Predictive diagnostics: The depth of flavor conversion ($P^x_v + iP^y_v = \Psi_v$8) is extractable directly from the real and imaginary parts ($P^x_v + iP^y_v = \Psi_v$9, $$\dot D_v = \frac{i}{2}[(\Psi_0 - v\Psi_1)\Psi_v^* - (\Psi_0^* - v\Psi_1^*)\Psi_v],$$0) of the linear eigenfrequency without solving the full nonlinear equations (Padilla-Gay et al., 2021). This enables rapid assessment in multidimensional supernova/merger simulations.
• Energy-exchange mechanism: Only resonant neutrinos (with velocities matching the wave phase velocity) participate in energy exchange, driving and damping the flavor wave cyclically. This is the flavor analog of O’Neil’s nonlinear wave-particle saturation mechanism in plasmas.
• Relation to solitons and waves: The flavor pendulum is the homogeneous, infinite-speed limit of the full fast-flavor nonlinear wave equation. "Solitons" emerge as limiting cases ("one-swing pendulum")—either temporal or spatial—when the seed amplitude vanishes (Fiorillo et al., 2023).
• Fragility to perturbations: Imposed homogeneity or axisymmetry is a delicate symmetry; arbitrarily small spatial inhomogeneities or angular asymmetries damp the pendulum via mode-cascade, multipole diffusion, or wave-number mixing (Mangano et al., 2014, Bhattacharyya et al., 2022).

In environments with two-beam or multi-azimuth-angle (MAA) structure, spontaneous symmetry breaking is generic: even in the normal-hierarchy case, arbitrary seeds away from perfect symmetry project onto unstable pendulum modes (Raffelt et al., 2013).

6. Limits, Extensions, and Generalizations

• Weak vs. strong instability: In weakly unstable regimes ($$\dot D_v = \frac{i}{2}[(\Psi_0 - v\Psi_1)\Psi_v^* - (\Psi_0^* - v\Psi_1^*)\Psi_v],$$1), the pendulum model remains quantitatively accurate for many cycles ($$\dot D_v = \frac{i}{2}[(\Psi_0 - v\Psi_1)\Psi_v^* - (\Psi_0^* - v\Psi_1^*)\Psi_v],$$2) before dephasing. Strongly unstable regimes activate large portions of momentum space (beyond the resonant band), so pendular dynamics give way to non-pendulum multi-mode saturation (Fiorillo et al., 26 Jan 2026, Fiorillo et al., 21 Jan 2026).
• Spectral diversity in collisional conversion: For multi-energy neutrino gases, collision-induced flavor conversion exhibits two branches: one with flavor equipartition for high-energy neutrinos, another where low-energy modes "swap" flavor before being damped. The branch realized depends on the sign of the average collision rate difference (α) (Zaizen, 13 Feb 2025).
• Three-flavor generalization: The two-flavor pendulum maps onto the "Bloch triangle" for three flavors; depolarization and flavor swapping recipes can be adapted with weighted projection factors for realistic supernova flux modeling (Bhattacharyya et al., 2022).
• Relation to coffee-sloshing and parametric excitation: Nonlinear planar pendula subject to parametric excitation (as in sloshing or mixing) display resonance-enhanced amplitudes near twice the natural frequency (Guarín-Zapata, 2021). The "flavor pendulum" metaphor in such systems provides analytic tools for optimizing mixing and excitation amplitude.

7. Summary Table: Flavor Pendulum Regimes and Key Properties

Flavor Pendulum Type	Driving Mechanism	Integrability	Saturation/Resting Behavior	Astrophysical Relevance
Slow/homogeneous	vacuum mixing (ω)	Yes (Gaudin)	periodic exchanges, solitons	early universe, SN
Fast/homogeneous	ν–ν refraction (μ)	Yes (Gaudin)	periodic/soliton oscillations	SN, mergers
Single-wave (2-beam)	spatial Fourier K	Partial	exact pendulum for 2-beam	idealized models
Single-wave (multi-v)	spatial Fourier K	No	approximate (weak γ) pendulum	SN, merger edge
Collisional	decoherence (Γ)	Yes/No	damped pendulum, steady swap	SN, mergers

The flavor pendulum, in its various mechanistic realizations, remains a central analytic structure for understanding collective neutrino flavor evolution, nonlinear saturation, and the impact of environmental perturbations. Its scope now includes both Hamiltonian and non-Hamiltonian channels ("collisional pendulum"), exact and approximate integrable sectors, and direct tie-ins to observable supernova/merger neutrino signal phenomenology (Fiorillo et al., 21 Jan 2026, Fiorillo et al., 26 Jan 2026, Padilla-Gay et al., 2022, Johns et al., 2023, Padilla-Gay et al., 2021, Fiorillo et al., 2024, Zaizen, 13 Feb 2025, Bhattacharyya et al., 2022, Mangano et al., 2014, Raffelt et al., 2013).