- The paper demonstrates that steep radial lepton gradients can suppress fast neutrino-flavor conversion by driving non-adiabatic drift of the unstable dispersion-relation branch.
- The simulations show that stronger Hamiltonian attenuation increases apparent suppression, while high matter potentials can also hide conversion through unresolved high-frequency modes.
- The authors propose an adiabaticity diagnostic based on growth rates, velocity, and radial matter gradients to distinguish physical suppression from numerical artifacts in global transport models.
Overview
Fast neutrino-flavor conversion (FFC), triggered by fast flavor instabilities (FFI) arising from ELN-XLN angular crossings, has been studied extensively in local boxes, but its coupling to global background geometry remains a central open problem in quantum kinetic neutrino transport. In "Fast Neutrino-Flavor Conversion with Attenuation and Global Lepton Gradient" (2604.13617), Zaizen and Nagakura perform global, spherically symmetric quantum kinetic simulations with radially varying matter backgrounds and an attenuated oscillation Hamiltonian, and show that steep radial lepton gradients can suppress FFC entirely. Crucially, they demonstrate that this suppression is highly sensitive to the attenuation parameter ξ, implying that attenuation—introduced purely as a numerical device—can artificially overestimate the impact of background variation and produce spurious stabilization.
Global simulation setup
The authors solve the QKE for the neutrino density matrix in spherical symmetry using their new GPU-accelerated code {\tt GANTS-QK}, which employs fifth-order WENO reconstruction on non-uniform Gauss-Legendre angular grids, SSP-RK(5,4) time integration with CFL number 0.4, and multi-GPU MPI parallelization. The oscillation Hamiltonian is rescaled by an attenuation parameter ξ∈{10−4,4×10−4,10−2,1}, a standard trick to relax the scale disparity between oscillation and transport scales (2604.13617). The matter potential follows a power law λ(r)=λ0(r/Rin)m with m∈{0,−1,−3} and λ0/μ0 up to 30, where μ0 is the self-interaction strength at the inner boundary. Dirichlet boundary conditions inject monochromatic (Eν=12 MeV) νe and νˉe distributions whose crossing at v=0.5 is progressively forward-peaked by momentum advection as neutrinos propagate outward, narrowing the ELN-XLN crossing and weakening the local instability at larger radii.
Suppression of flavor conversion by radial gradients
The central numerical result is stark: with strong attenuation (ξ∈{10−4,4×10−4,10−2,1}0) and a flat matter profile (ξ∈{10−4,4×10−4,10−2,1}1), robust flavor conversion develops across the domain, whereas with a steep gradient (ξ∈{10−4,4×10−4,10−2,1}2) flavor coherence fails to build any coherent wave pattern beyond ξ∈{10−4,4×10−4,10−2,1}3 km and conversion is completely suppressed. The degree of suppression depends strongly on ξ∈{10−4,4×10−4,10−2,1}4: at ξ∈{10−4,4×10−4,10−2,1}5, conversion occurs in all background models, with onset radii below 50.1 km, indicating that the physical effect of the gradient is modest when growth is fast enough. This immediately establishes the paper's cautionary conclusion: attenuation artificially inflates the propagation distance required for linear saturation, thereby increasing exposure to background variation before saturation, so strongly attenuated global simulations can misreport gradient-induced suppression as physical.
The authors also identify a purely numerical suppression channel: for large matter potentials (ξ∈{10−4,4×10−4,10−2,1}6), the unstable branch migrates to high frequencies that approach the maximum resolvable frequency set by the CFL time step, ξ∈{10−4,4×10−4,10−2,1}7. Increasing radial resolution restores conversion, confirming that under-resolution—not physics—caused the missing growth in those models.
Adiabaticity of flavor waves
To interpret these results, the authors develop a local dispersion-relation analysis of the unstable branch ξ∈{10−4,4×10−4,10−2,1}8 evaluated at each radius from the classical steady state. Because the background lepton density shifts the branch primarily along the ξ∈{10−4,4×10−4,10−2,1}9 axis while the group velocity λ(r)=λ0(r/Rin)m0 and growth rate λ(r)=λ0(r/Rin)m1 are nearly frame-independent and insensitive to the phase shift, a propagating flavor wave remains unstable only if it tracks the drifting branch adiabatically. They define the adiabaticity parameter
λ(r)=λ0(r/Rin)m2
with conversion possible only when λ(r)=λ0(r/Rin)m3. Under attenuation, λ(r)=λ0(r/Rin)m4 scales as λ(r)=λ0(r/Rin)m5 because λ(r)=λ0(r/Rin)m6 carries one power of λ(r)=λ0(r/Rin)m7 while λ(r)=λ0(r/Rin)m8 carries two—a clean analytic explanation of why stronger attenuation makes gradient-induced suppression more likely. The computed adiabaticity profiles reproduce the simulated onset radii: models violating λ(r)=λ0(r/Rin)m9 before the observed onset radius show suppressed or delayed conversion, and vice versa.
Notably, even with a flat matter profile (m∈{0,−1,−3}0), spherical advection itself drives m∈{0,−1,−3}1 at large radii through the combined effects of decreasing neutrino density and sharpening angular crossings. Global geometry alone can therefore break adiabaticity, independent of any imposed matter gradient.
The authors contrast this measure with the comoving-frame treatment of Bhattacharyya et al., which uses the width of the unstable branch m∈{0,−1,−3}2 in place of m∈{0,−1,−3}3 in the denominator. They argue the two measures are not interchangeable: the comoving-frame picture implicitly assumes the dispersion relation connects across radii, which fails when the branch drifts in the m∈{0,−1,−3}4–m∈{0,−1,−3}5 plane, whereas their formulation explicitly tracks m∈{0,−1,−3}6.
Approximate diagnostic
Since full linear stability analysis requires complete momentum-space information unavailable in classical CCSN/BNSM transport data, the authors derive a practical estimator requiring only three ingredients: the empirical two-beam growth rate built from the positive and negative ELN-XLN angular integrals, a representative velocity taken as either the angular-crossing direction m∈{0,−1,−3}7 or the flavor-averaged flux factor m∈{0,−1,−3}8, and the approximation m∈{0,−1,−3}9. The resulting dimensionless estimate,
λ0/μ00
can be evaluated directly from hydrodynamical simulation output. The authors state it is accurate to within a factor of a few, not order-of-magnitude, though they concede it is a rough diagnostic. They also note a caveat: the angular-crossing direction happens to track the group velocity in their models but does not universally—for example, in CCSN preshock regions the group velocity substantially exceeds λ0/μ01—so the flux-factor proxy is generally preferable.
Limitations and open questions
Several limitations are stated plainly by the authors. First, the initial neutrino distributions are fixed throughout the parametric study; realistic CCSN environments feature narrower angular crossings than those adopted here, weakening the instability and potentially altering the adiabaticity competition, so generalization to realistic backgrounds remains untested. Second, the fluid bulk velocity is neglected, restricting the analysis to static backgrounds. Third, the approximate formula inherits the empirical two-beam growth-rate estimate and the assumption that matter dominates λ0/μ02; both approximations would require validation against dynamical models. Finally, the finding that attenuation artificially enhances suppression raises a question the paper leaves open: whether collisional and slow flavor instabilities—which have longer growth timescales and may exhibit intrinsically global growth—can be reliably studied with attenuated Hamiltonians at all, or whether their apparent absence in global simulations reflects the same artificial mechanism identified here for FFC.
Conclusion
This work establishes that radial lepton gradients can physically suppress fast neutrino-flavor conversion via non-adiabatic drift of the unstable dispersion-relation branch, provides a quantitative adiabaticity criterion computable from classical transport quantities, and delivers a clear methodological warning: attenuation and finite resolution can each manufacture artificial suppression of collective flavor conversion in global simulations. Any global quantum kinetic study reporting absent flavor conversion under spatially varying backgrounds should verify, via the proposed adiabaticity diagnostic and resolution tests, that the suppression is physical rather than numerical.