Gradient-Produced Neutrinos

Abstract: Sufficiently strong electric fields can produce charged-particle pairs via the Schwinger effect. We argue that steep matter-density gradients, as can arise in neutron star interiors, would analogously produce neutrino-antineutrino pairs. We then discuss observational signatures of these gradient-produced (anti)neutrinos and how they could provide new probes of neutron-star structure and baryon-dense QCD.

• The paper introduces a generalized Schwinger effect where steep matter-potential gradients in neutron stars produce neutrino-antineutrino pairs.
• It employs a Sauter profile and Bogoliubov transformations to derive production rates that scale with the matter-potential jump, distinguishing sharp and gradual interfaces.
• The study reveals that gradient-produced neutrinos modify neutron star cooling curves by inducing heating plateaus or accelerated cooling, offering observable astrophysical signatures.

Gradient-Induced Neutrino Pair Production in Neutron Stars

Generalized Schwinger Effect for Neutrinos

The paper "Gradient-Produced Neutrinos" (2604.21968) explores a generalization of the Schwinger effect, extending it from gauge fields to inhomogeneous matter density backgrounds, particularly relevant in neutron star (NS) interiors. In the standard Schwinger mechanism, sufficiently strong electric fields destabilize the vacuum, resulting in spontaneous electron-positron pair production. The authors consider the analogous case where steep gradients in the effective matter potential $V(x)$ produce neutrino-antineutrino pairs via mixing of positive- and negative-frequency modes in the Dirac equation.

The interaction Lagrangian formulated is $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$, where $j_{\rm ext}^\mu$ is a static external current, e.g., the matter potential from neutrino forward scattering in baryon-rich matter. For Dirac neutrinos, $c_V = c_A = 1/2$. In steep density/composition gradients, $V(x)$ exhibits critical jumps, analogously to phase boundaries in 1st-order phase transitions within NSs.

Calculation of Gradient-Induced Pair Production

The spatial profile of the matter potential is modeled as a Sauter (tanh) profile: $V(x) = \frac{1}{2}\Delta V\tanh(z/l)$, with the jump $\Delta V$ and width $l$. A necessary condition for unsuppressed pair production is that $\Delta V$ overcomes the mass threshold ($\Delta V > 2 m_\nu$). The pair production rate per unit area is derived using Bogoliubov transformations and mapped directly onto the known formalism for Schwinger pair production by spatially varying vector potentials [Fedotov:2022ely, Chervyakov:2009bq]. In the sharp (supercritical) limit, $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$0, the rate scales cubically with the jump:

$$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$1

For gradual transitions, $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$2, the rate is linearly suppressed by $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$3. Numerical results, as well as analytic approximations, confirm the crossover from cubic to quadratic scaling with $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$4 as the interface becomes less sharp.

Figure 1: Pair-production rate per unit area $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$5 versus transition width $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$6, exhibiting higher rates for sharper matter-potential jumps.

Consequences in Neutron Star Structure

The application to neutron stars centers on the scenario where a sharp density/composition jump exists (e.g., hadronic-to-quark matter phase boundary), inducing a supercritical matter-potential gradient. The benchmark parameters adopt $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$7, $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$8 km, $$-\mathcal{L}_{\rm int}=j_{\rm ext}^\mu\bar{\psi}\gamma_\mu(c_V-c_A\gamma_5)\psi$$9, $j_{\rm ext}^\mu$0, and $j_{\rm ext}^\mu$1 km. The resulting neutrino production rate is highly efficient:

$j_{\rm ext}^\mu$2

Neutrinos are confined within the jump radius—accumulating and becoming degenerate as Pauli blocking saturates the phase space, with a total degenerate number scaling as $j_{\rm ext}^\mu$3.

Figure 2: Schematic of pair production at a sharp interface in NS, showing bound neutrino formation and thermal upscattering escape.

Heating and Cooling Channels Due to Gradient-Produced Neutrinos

Bound neutrinos can be depleted by two primary mechanisms: absorption in matter and upscattering, which enables escape. The absorption mean free path ($j_{\rm ext}^\mu$4) is controlled by available phase space and exothermic transitions, sensitive to chemical potential imbalances. If upscattered above the potential barrier (when $j_{\rm ext}^\mu$5), neutrinos escape, carrying thermal energy and cooling the core.

The paper constructs a thermal evolution model with terms for photon luminosity ($j_{\rm ext}^\mu$6), standard thermal-neutrino emission ($j_{\rm ext}^\mu$7), gradient-produced scattering cooling ($j_{\rm ext}^\mu$8), and absorption heating ($j_{\rm ext}^\mu$9):

$c_V = c_A = 1/2$0

Luminosities $c_V = c_A = 1/2$1 and $c_V = c_A = 1/2$2 can exceed $c_V = c_A = 1/2$3 and $c_V = c_A = 1/2$4 at relevant NS surface temperatures ($c_V = c_A = 1/2$5).

Figure 3: Heating and cooling luminosities as a function of $c_V = c_A = 1/2$6, highlighting dominance over conventional photon and neutrino channels in specific temperature ranges.

NS Cooling Curve Modification and Observational Signatures

The gradient-neutrino mechanisms can produce either heating plateaus or accelerated cooling in NS temperature-age curves, depending on which depletion mechanism dominates. The heating from absorption manifests as a late-time temperature plateau ($c_V = c_A = 1/2$7 yr), whereas enhanced upscattering induces a characteristic knee feature in the cooling history ($c_V = c_A = 1/2$8 yr).

Figure 4: NS cooling curves plotted as $c_V = c_A = 1/2$9 vs age $V(x)$0, illustrating the effect of gradient-produced neutrino channels.

Surface-temperature measurements of old NSs—particularly within $V(x)$1 K—by IR telescopes (JWST, ELT, TMT) can test these predictions. The anomalously warm NSs may be explained by enhanced gradient-induced heating, contingent on substantial chemical imbalances.

Spectral Properties of Gradient-Produced Neutrinos

In the sharp limit, the energy spectra of produced particles peak near $V(x)$2, below the potential barrier, implying confinement except in rare cases of upscattering. This spectral function is symmetric in $V(x)$3, arising from the pair-production formalism.

Figure 5: Kinetic-energy spectra for gradient-produced particles in the sharp interface limit, demonstrating confinement by NS matter potential.

Robustness and Implications

Absorption-heating signals are robust under changes to interface width and can persist even for relatively small $V(x)$4, provided at least one neutrino mass eigenstate is sufficiently light. Observations of cold NSs can, under a reliable model for chemical potential imbalance, set lower bounds on the lightest neutrino mass. The effect persists even for smooth density profiles, with irreducible gradient production, although rates are greatly enhanced by abrupt phase transitions.

The theoretical implications extend to broader classes of backgrounds (not limited to gauge or matter potentials): any sufficiently steep inhomogeneous vector field can drive pair production. This has potential applications in probing BSM fields, new fermion species, and relic neutrino backgrounds.

Conclusion

This study establishes the existence of Schwinger-like neutrino-antineutrino pair production from steep matter-potential gradients as generic in neutron star interiors. For supercritical jumps, the effect yields heating and cooling channels capable of dominating NS late-stage thermal evolution and may provide a probe of baryon-dense QCD. The predicted observational signatures—plateaus and knees in NS cooling curves—are accessible by IR telescopes and radio surveys. The mechanism is theoretically robust and generalizable, prompting further work into NS phenomenology, flavor/chiral effects, interface modeling, and the search for new physics associated with low-energy neutrino matter interactions.