
Collision-induced flavor instability in dense neutrino gases with energy-dependent
scattering
Yu-Chia Lin1, 2, 3, ∗and Huaiyu Duan1, †
1Department of Physics and Astronomy, University of New Mexico, Albuquerque, New Mexico 87131, USA
2Department of Astronomy/Steward Observatory, University of Arizona, Tucson, Arizona 85721, USA
3Department of Physics, University of Arizona, Tucson, Arizona 85721, USA
(Dated: April 27, 2023)
We investigate the collision-induced flavor instability in homogeneous, isotropic, dense neutrino
gases in the two-flavor mixing scenario with energy-dependent scattering. We uncover a simple
expression of the growth rate of this instability in terms of the flavor-decohering collision rates and
the electron lepton number distribution of the neutrino. This growth rate is common to the neutrinos
and antineutrinos of different energies, and is independent of the mass-splitting and vacuum mixing
angle of the neutrino, the matter density, and the neutrino density, although the initial amplitude
of the unstable oscillation mode can be suppressed by a large matter density. Our results suggest
that neutrinos are likely to experience collision-induced flavor conversions deep inside a core-collapse
supernova even when both the fast and slow collective flavor oscillations are suppressed.
I. INTRODUCTION
Neutrinos help shape the physical and chemical evolu-
tion of the early universe, core-collapse supernovae, and
neutron star mergers where they are copiously produced.
Flavor oscillations, which alter the flavor composition of
the neutrinos (see, e.g., Ref. [1] for a review), can have a
significant impact on the physical conditions in these in-
teresting astrophysical environments. In addition to the
well-known vacuum oscillations [2, 3] and the Mikheyev-
Smirnov-Wolfenstein (MSW) effect [4, 5], the ambient
neutrinos further change the refraction of the neutrinos in
these extreme environments [6–8]. As a result, the dense
neutrino gas can experience a collective flavor transfor-
mation because of the tight coupling among the neutrinos
themselves [9–11]. (See, e.g., Ref. [12, 13] for reviews on
this topic and the references therein.) It has been shown
that the so-called fast flavor conversions, a special type of
collective flavor transformation that arises on the scales
of centimeters to meters, can take place near or even be-
low the neutrino-decoupling layer of a neutrino-emission
compact object [14, 15]. (See also Ref. [16] for a review
and the references therein.) An interesting recent devel-
opment in the research of collective neutrino oscillations
is that neutrino collisions, which are usually thought to
damp neutrino oscillations, are shown to be able to in-
duce flavor conversions [17–26].
In this work, we focus on the collision-induced flavor
conversions in homogeneous and isotropic neutrino gases.
We first follow the pioneering work of Johns [20] and
solve the flavor evolution in a mono-energetic neutrino
gas. By comparing the numerical and analytical solu-
tions, we demonstrate the dependence of the collision-
induced instability on the effective mixing angle and the
density of the neutrino (Sec. II). We then consider the ef-
∗yuchialin@arizona.edu
†duan@unm.edu
fect of energy-dependent neutrino scattering and derive
a simple expression for the exponential growth rate of
the collision-induced flavor instability for which we also
provide a few illustrative numerical examples (Sec. III).
We conclude by summarizing our results and discussing
their implications (Sec. IV).
II. MONO-ENERGETIC NEUTRINO GAS
A. Physics model
The flavor content of a dense neutrino gas can be rep-
resented by the flavor density matrix ρwhose diagonal
elements are the neutrino occupancies in different weak-
interaction states and the off-diagonal elements are the
coherences among those states [27]. For a homogeneous
and isotropic environment, the flavor evolution of the
neutrino gas can be solved from the following equation:
˙ρ=−i[H, ρ] + C,(1)
where Hand Care the flavor-evolution Hamiltonian and
the collision term for the neutrino, respectively. In the
two-flavor mixing scenario, say between the eand xfla-
vors, one can expand ρin terms of the 2 ×2 identity
matrix σ0and the Pauli matrices σi(i= 1,2,3) so that
ρ∝σ0P0+P·σ,(2)
where Pis the flavor polarization (Bloch) vector, and P0
is proportional to the total density of the eand xflavor
neutrinos. The corresponding quantities ¯ρ,¯
P0, and ¯
P
can be defined for the antineutrinos in a similar way.
We first consider a mono-energetic neutrino gas of a
single vacuum oscillation frequency ω= ∆m2/2E, where
∆m2and Eare the mass-squared difference and energy
of the neutrino, respectively. We assume a minimal col-
lisional model that damps the flavor coherence without
changing the particle numbers [28]. In this model, Eq. (1)
arXiv:2210.09218v3 [hep-ph] 26 Apr 2023