Flavor-Wave Synchronization

• Flavor-Wave Synchronization is the collective phase locking of neutrino flavor waves, establishing a steady, collisionless equilibrium in dense neutrino gases.
• It employs Fourier analysis and nonlinear coupling to demonstrate how phase alignment prevents secular evolution of the system’s mean field.
• This mechanism has astrophysical implications, influencing neutrino distributions and energy transfer in core-collapse supernovae and neutron-star mergers.

Searching arXiv for the specified paper and closely related work on fast neutrino flavor conversion. Flavor-wave synchronization is the phenomenon by which the spatially inhomogeneous “flavor waves” that develop in a dense neutrino gas lock their relative phases so as to support a time-independent, or collisionless, flavor equilibrium of the whole system. In the analysis of fast flavor conversion (FFC), it is presented alongside nonlinear Landau damping and collisionless equilibria as an analogue of time-honored nonlinear phenomena in plasmas and self-gravitating systems, and as a mechanism likely important for the unsolved puzzle of neutrino oscillations in core-collapse supernovae and neutron star mergers (Liu et al., 30 Sep 2025).

1. Definition in the two-beam fast flavor-conversion model

The setting is a one-dimensional gas of neutrinos streaming to the right (“R”) and to the left (“L”) in a periodic box $r\in[0,L]$. At each point $r$ and time $t$, the local flavor state of each beam is described by a polarization (Bloch) vector $\bm P_I(r,t)$, with $I=R,L$. In FFC, these vectors evolve under a purely mean-field self-interaction,

$$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$

where $J\neq I$, and the overall interaction strength is set to $\mu=1$ (Liu et al., 30 Sep 2025).

A Fourier decomposition,

$$\bm P_I(r,t) =\sum_{n\in\mathbb Z}\bm P_{I,n}(t)\,e^{i k_n r}, \quad k_n\equiv\frac{2\pi n}{L},$$

identifies the non-zero modes $n\neq0$ as flavor waves of wavenumber $r$0. The spatially averaged quantities are

$r$1

Flavor-wave synchronization is the special situation in which, despite the existence of these inhomogeneous modes, $r$2 and $r$3 remain strictly constant in time.

Two conditions must be met. First, the zeroth modes are in mixing equilibrium,

$r$4

so that the $r$5 dynamics by itself would be static. Second, the flavor-wave modes are phase aligned in just the right way that they do not feed back on $r$6. Specifically,

$r$7

and this vanishes if all modes satisfy

$r$8

In that case, each Fourier pair is circularly polarized and the two beams share the same phase, or are exactly anti-phase. The inhomogeneous flavor waves then “synchronize,” and the box average is truly steady.

2. Fourier-space equations and the dispersion problem

Starting from the real-space equations of motion,

$r$9

one Fourier-transforms to obtain, for each integer $t$0,

$t$1

$t$2

Because $t$3 is real, $t$4. From the $t$5 component, one finds

$t$6

Since $t$7 in mixing equilibrium, the first term vanishes and only the sum over inhomogeneous modes remains (Liu et al., 30 Sep 2025).

For small-amplitude flavor waves, the analysis introduces $t$8 and linearizes about the background $t$9. This yields a Schrödinger-like linear system

$\bm P_I(r,t)$0

where the non-Hermitian “Hamiltonian” $\bm P_I(r,t)$1 depends on $\bm P_I(r,t)$2. The standard dispersion relation,

$\bm P_I(r,t)$3

yields four eigenvalues $\bm P_I(r,t)$4. Two of them come in complex-conjugate pairs when $\bm P_I(r,t)$5 lies within an unstable band, signaling exponential growth at rate $\bm P_I(r,t)$6. In the two-beam model, instability occurs for

$\bm P_I(r,t)$7

and the maximal growth rate is $\bm P_I(r,t)$8 in these units.

3. Single-wave synchronization and the flavor-wave pendulum

In the single-wave regime, the box size $\bm P_I(r,t)$9 is chosen so that only one mode $I=R,L$0 lies in the unstable band. The nonlinear coupling is then effectively between $I=R,L$1 and $I=R,L$2 alone. In this case the dynamics reduces to a flavor-wave pendulum with two degrees of freedom,

$I=R,L$3

with

$I=R,L$4

and

$I=R,L$5

Equation (4) admits librating (trapped) or circulating orbits in the $I=R,L$6 phase plane, all periodic and regular (Liu et al., 30 Sep 2025).

The fixed points occur exactly at $I=R,L$7, that is, at phase alignment. These fixed points define neutrino BGK modes, described as collisionless equilibria supported by phase-locked flavor waves. The stable fixed points of Eq. (4), satisfying $I=R,L$8 at $I=R,L$9, are also exact stationary solutions of the full nonlinear equations of motion. Each such solution has a self-consistent wave amplitude and phase such that $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$0 co-precess about $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$1 with no secular drift.

The numerical example for a box of size $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$2, where only $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$3 is unstable, compares the full QKE evolution with the reduced pendulum equation. Both show perfectly periodic oscillations of $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$4 and $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$5. The corresponding phase-space orbits are regular, with librating and circulating trajectories around the centers $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$6, identified as synchronized BGK modes.

4. Many-wave synchronization, dephasing, and collisionless equilibration

In the many-wave regime, $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$7 is chosen so that tens or hundreds of modes lie in the unstable band. Each mode $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$8 then undergoes its own linear growth and later nonlinear saturation. Different modes acquire different frequencies and growth rates $$(\partial_t + v_I\,\partial_r)\,\bm P_I(r,t) \;=\; 2\,\bm P_J(r,t)\,\times\,\bm P_I(r,t)\,, \quad (v_R=+1,\;v_L=-1),$$9, so their mutual couplings lead to dephasing and an effectively irreversible relaxation of $J\neq I$0 (Liu et al., 30 Sep 2025).

The numerical examples with box sizes $J\neq I$1 and $J\neq I$2 correspond to 11 and 450 unstable modes respectively, while $J\neq I$3 corresponds to 1 unstable mode. As $J\neq I$4 grows, exact recurrences give way to effectively irreversible relaxation, but $J\neq I$5 nonetheless at late times. This distinguishes single-wave synchronization, which is exact and periodic, from many-wave synchronization, which is statistical and quasi-steady.

In the many-wave case, the system settles into a quasi-steady state in which the time-averaged distributions of each phase difference

$J\neq I$6

become sharply peaked at $J\neq I$7 and $J\neq I$8. During the nonequilibrium window, those distributions are broad and dephased. During the equilibrium window, they peak sharply at $J\neq I$9 and $\mu=1$0, revealing the many-wave flavor-wave synchronization that underpins the collisionless equilibrium. The resulting phase-aligned many-wave state sustains $\mu=1$1 even though each $\mu=1$2 is small but coherent.

5. Relation to inverse and nonlinear Landau damping

The early linear phase is described as linear Landau (inverse) damping. Each unstable eigenmode grows as

$\mu=1$3

driven by the crossing of the electron and positron, here $\mu=1$4 and $\mu=1$5, dispersion curves (Liu et al., 30 Sep 2025).

The nonlinear stage is described as nonlinear Landau damping. In the single-wave pendulum, once $\mu=1$6 grows enough to backreact on $\mu=1$7, the growth slows and reverses, with energy sloshing back and forth between the mean field $\mu=1$8 and the wave $\mu=1$9. This bidirectional energy exchange is identified as the analogue of O’Neil’s trapped-particle pendulum. It explains why flavor waves do not simply damp to zero but instead can execute sustained periodic oscillations.

Within this interpretation, flavor-wave synchronization is not merely a kinematic coincidence of phases. It is the phase structure associated with dynamical equilibria of the self-interacting neutrino field. If the modes co-precess and preserve the phase relation required by $$\bm P_I(r,t) =\sum_{n\in\mathbb Z}\bm P_{I,n}(t)\,e^{i k_n r}, \quad k_n\equiv\frac{2\pi n}{L},$$0 or $$\bm P_I(r,t) =\sum_{n\in\mathbb Z}\bm P_{I,n}(t)\,e^{i k_n r}, \quad k_n\equiv\frac{2\pi n}{L},$$1, then the nonlinear coupling does not drive secular evolution of the box average. If they dephase, the feedback on $$\bm P_I(r,t) =\sum_{n\in\mathbb Z}\bm P_{I,n}(t)\,e^{i k_n r}, \quad k_n\equiv\frac{2\pi n}{L},$$2 drives damping or growth. A plausible implication is that synchronization provides the phase-space organization corresponding to collisionless neutrino equilibria in the same sense that BGK-type structures organize collisionless equilibria in other mean-field kinetic systems.

6. Physical interpretation and astrophysical implications

Flavor-wave synchronization is the collective locking of the relative oscillation phases of the self-interacting neutrino field across different spatial harmonics. Physically, neutrino self-interactions produce a nonlinear coupling among all Fourier modes. If those modes co-precess, nothing in the coupling drives the box-averaged occupation, and the system is in a dynamical equilibrium. If they dephase, the feedback on the mean field $$\bm P_I(r,t) =\sum_{n\in\mathbb Z}\bm P_{I,n}(t)\,e^{i k_n r}, \quad k_n\equiv\frac{2\pi n}{L},$$3 drives secular evolution (Liu et al., 30 Sep 2025).

For core-collapse supernovae and neutron-star mergers, the mechanism implies that once fast flavor instabilities have grown and saturated across a spectrum of modes, the system will tend toward a phase-synchronized state rather than incoherent decoherence. The paper states that such synchronization can alter the angular and spectral distributions of neutrinos, with consequences for neutrino heating behind the stalled shock in a supernova, the electron-fraction ($$\bm P_I(r,t) =\sum_{n\in\mathbb Z}\bm P_{I,n}(t)\,e^{i k_n r}, \quad k_n\equiv\frac{2\pi n}{L},$$4) of outflows and hence nucleosynthesis yields, the predicted neutrino signal in terrestrial detectors, and the large-scale evolution of the neutrino radiation field in merger disks.

The same perspective also identifies a target for subgrid and coarse-grained modeling of neutrino quantum kinetics: one need not track all phase-space filaments, only the collective synchronized modes and their approach to BGK-type equilibrium. This suggests a modeling strategy in which the physically relevant late-time state is characterized not by complete incoherence, but by collisionless equilibria supported by phase-locked or statistically phase-aligned flavor waves.