Charged-current neutrino opacity within the relativistic Hartree-Fock framework for astrophysical simulations of core-collapse supernovae and binary neutron star mergers

This paper introduces a relativistic Hartree-Fock framework with momentum-dependent nuclear interactions to improve the calculation of charged-current neutrino opacities for astrophysical simulations, revealing significant discrepancies and substantial shifts in medium-dependent modifications compared to commonly used relativistic mean-field models.

The Gist

Imagine the heart of a dying star or the violent crash of two neutron stars as a cosmic pressure cooker. Inside this furnace, particles called neutrinos are born in massive numbers. These are ghostly particles that rarely interact with anything, but in these extreme environments, they act like the star's lifeblood: they carry away heat, transport energy, and help decide what new elements are forged in the fire.

To understand how these stars explode or merge, scientists run computer simulations. A critical part of these simulations is calculating how easily neutrinos can move through the dense soup of protons and neutrons inside the star. This "ease of movement" is called opacity. If the opacity is high, neutrinos get stuck (like trying to walk through a crowded concert); if it's low, they zip right through.

The Old Map vs. The New Map

For a long time, scientists used a standard map to calculate this opacity, called the Relativistic Mean-Field (RMF) model. Think of this model as a simplified map where every particle in the star is treated as if it's moving in a smooth, average ocean. It assumes that the "water" (the nuclear medium) affects all particles the same way, regardless of how fast they are swimming.

In this new paper, the authors say, "That map is too simple." They introduce a more detailed map called the Relativistic Hartree-Fock (RHF) approach.

The Analogy of the Traffic Jam:

• The RMF Model (Old Way): Imagine a highway where every car feels the same average traffic pressure. It doesn't matter if you are driving a sports car or a truck; the road treats you the same.
• The RHF Model (New Way): This model realizes that traffic is messy. A fast car feels the air differently than a slow truck. It accounts for the fact that particles have specific speeds and that their interactions depend on exactly how fast they are moving and in what direction. It's like realizing that in a real traffic jam, your experience depends heavily on your specific speed and the cars immediately around you.

What They Found

When the authors applied this new, more detailed "traffic-aware" model to calculate neutrino opacity, they found some surprising differences compared to the old model:

1. The "Ghost" vs. The "Wall": For certain types of neutrinos (electron neutrinos), the new model suggests they get stuck much more easily in the star's core than the old model predicted. It's as if the old map said the road was clear, but the new map reveals a hidden wall.
2. The Reverse for Anti-Neutrinos: For the opposite type of particle (anti-neutrinos), the new model suggests they can actually move more freely than the old model thought. The "wall" is less of a barrier for them.
3. The Speed Matters: The biggest difference comes from the fact that in the new model, the "density" of the star changes depending on how fast the particles are moving. In the old model, the density was static. This speed-dependence shifts the energy levels where neutrinos can be absorbed, effectively changing the "rules of the game" for how the star evolves.

Why This Matters for the Simulation

The authors didn't just change the math for the sake of it; they showed that these changes are huge.

• In the old simulations, the difference between how neutrinos and anti-neutrinos behave was exaggerated.
• In the new simulations, the behavior of these two particle types is actually more similar to each other than previously thought, but the magnitude of their interaction with the star's matter is different.

Think of it like tuning a musical instrument. The old model was slightly out of tune, making the "notes" (the energy and flow of neutrinos) sound too different from each other. The new model tightens the strings, bringing the pitch closer to what the physics of the nuclear medium actually dictates.

The Bottom Line

This paper doesn't claim to have solved how stars explode or how neutron stars merge. Instead, it provides a more accurate instrument for the scientists who do those simulations. By including the fact that particles interact differently based on their speed (momentum), the authors have created a more realistic description of the "nuclear soup" inside these cosmic events.

They found that the old, simpler models were missing a crucial detail: the "personality" of the particles changes based on how fast they are moving. Ignoring this leads to significant errors in predicting how much heat is trapped or released, which is vital for understanding the life and death of stars.

In short: The authors built a better microscope to look at the tiny interactions inside a dying star, and they found that the view is much more complex—and different from what we thought—than the old, blurry picture allowed.

Technical Summary

Technical Summary: Charged-current neutrino opacity within the relativistic Hartree-Fock framework

Problem Statement
Neutrinos and their weak interactions are fundamental to the physics of core-collapse supernovae and binary neutron star mergers, governing shock revival, nucleosynthesis, and the transport of lepton number and entropy. Accurate astrophysical simulations require precise descriptions of weak rates, particularly charged-current (CC) processes involving the nuclear medium. Previous simulations have largely relied on the Relativistic Mean-Field (RMF) approximation (Hartree level), which treats nucleons with momentum-independent self-energies. However, this approach neglects explicit momentum-dependent nuclear interactions arising from exchange (Fock) terms. The authors identify a need to improve the description of the nuclear medium to calculate CC weak rates more accurately, specifically by incorporating momentum-dependent nucleon self-energies.

Methodology
The authors develop a formalism for charged-current neutrino opacity within the Relativistic Hartree-Fock (RHF) framework.

• Formalism: They utilize the real-time formalism and Kobes-Semenoff rules to derive the neutrino absorption rate from the neutrino self-energy. The interaction is modeled via $W^{\pm}$ boson exchange between leptons and nucleons.
• Nuclear Medium: Unlike the standard RMF approach, the RHF treatment explicitly includes momentum-dependent nucleon self-energies ($\Sigma^S, \Sigma^0, \Sigma^V$) derived from the exchange of $\sigma, \omega, \rho,$ and $\pi$ mesons. This leads to effective (dressed) nucleon masses, momenta, and dispersion relations that depend on the nucleon's momentum.
• Hadronic Current: The calculation incorporates the full hadronic vertex, including Leading Order (LO), Weak Magnetism (WM), Pseudoscalar (PS) terms, and momentum-dependent weak form factors (FF).
• Numerical Implementation: The authors evaluate the neutrino absorption rate for two representative conditions (A and B) relevant to supernovae and mergers, characterized by specific temperatures, baryon densities, and lepton fractions. They compare the full RHF results against:
  1. Common RMF parameterizations (TMA, NL3, DD2).
  2. A reduced RMF model ($RMF^*$) derived from the RHF parameters to isolate the effects of the momentum dependence and self-energy differences.
  3. Vacuum rates and semi-reduced RHF limits to validate the numerical implementation.
• Kinematics: The authors address the complexity introduced by momentum-dependent self-energies, which prevents the analytical removal of Dirac delta functions in the phase space integrals (a simplification possible in RMF). They employ numerical root-finding procedures (Newton-Raphson) to handle the energy conservation constraints within the integration.

Key Contributions

1. Derivation of RHF Opacity Expressions: The paper provides the first-principle expressions for charged-current weak rates within the RHF framework, explicitly accounting for momentum-dependent self-energies in the hadron polarization tensor and lepton-hadron tensor contraction.
2. Identification of RMF Limitations: The study demonstrates that the standard RMF approximation, which assumes identical scalar self-energies for neutrons and protons (in the absence of isovector-scalar $\delta$-mesons), fails to capture significant medium modifications present in the RHF approach.
3. Quantitative Discrepancies: The authors quantify the differences between RHF and RMF predictions, showing that the momentum dependence and the natural splitting of neutron and proton scalar self-energies in RHF lead to substantial shifts in the effective $Q$-values and the energy dependence of the opacity.

Results

• Neutrino vs. Antineutrino Asymmetry: The RHF approach consistently suppresses the opacity for electron and muon neutrinos ($\nu_e, \nu_\mu$) compared to RMF models, while enhancing the opacity for antineutrinos ($\bar{\nu}_e, \bar{\nu}_\mu$). This trend holds across various temperatures and densities.
• Effective $Q$-Value Shifts: The inclusion of momentum-dependent self-energies and the difference between neutron and proton scalar self-energies in RHF results in a shift of the effective medium-modified $Q$-value ($Q^*$). For example, under condition A, the $Q^*$ for $\nu_e$ absorption shifts from $\approx -16$ MeV (TMA) to $\approx -14.5$ MeV (RHF), while for $\bar{\nu}_e$ it shifts from $\approx +16$ MeV to $\approx +14.5$ MeV.
• Density Dependence: At intermediate densities, the RHF rates for $\nu_e$ are suppressed relative to RMF, while $\bar{\nu}_e$ rates are enhanced. At supranuclear densities, the RHF approach begins to dominate over all tested RMF models. The DD2 RMF model exhibits the strongest enhancement of $\nu_e$ opacity and suppression of $\bar{\nu}_e$ opacity among the RMF models, attributed to its large self-energies.
• Inverse Neutron Decay: The inverse neutron decay channel ($\bar{\nu}_e + p + e^- \to n$) is found to be largely suppressed within the RHF approximation compared to RMF.
• Role of Hadronic Current Terms: Higher-order contributions (Weak Magnetism, Pseudoscalar, Form Factors) significantly increase $\nu_e$ rates at high energies and densities, whereas their impact on $\bar{\nu}_e$ rates becomes negligible due to medium suppression.
• Reduced Model Validation: The $RMF^*$ model, which retains the RHF bulk properties but removes the momentum dependence and enforces equal scalar self-energies, qualitatively reproduces the RHF results. This confirms that the primary source of discrepancy between standard RMF and RHF is the difference in scalar self-energies between neutrons and protons and the momentum dependence of the interactions.

Significance
The paper claims that the inclusion of momentum-dependent nuclear interactions and the resulting differences in neutron and proton scalar self-energies are essential for accurate neutrino opacity calculations. The observed suppression of $\nu_e$ and enhancement of $\bar{\nu}_e$ opacities in the RHF framework, compared to standard RMF models, may have important consequences for the differences between neutrino and antineutrino fluxes and spectra in simulations of core-collapse supernovae. Specifically, these effects could influence the long-term deleptonization phase of the proto-neutron star. The authors conclude that these novel rates should be incorporated into future multi-dimensional simulations of supernovae and binary neutron star mergers to improve the reliability of predictions regarding nucleosynthesis and explosion dynamics.