$νp$-process in Core-Collapse Supernovae: Imprints of General Relativistic Effects

This study demonstrates that incorporating General Relativistic effects into core-collapse supernova models significantly enhances the efficiency of the $\nu$p-process, enabling a sufficiently massive progenitor to reproduce the full range of solar system $p$-nuclide abundances up to $^{102}$Pd, whereas Newtonian calculations fail to do so.

The Gist

The Cosmic Puzzle: Where Do the "Rare Earth" Elements Come From?

Imagine the universe as a giant kitchen. Most of the ingredients (elements) we see are made by two main chefs:

1. The Slow Cook (s-process): Takes its time, building heavy elements slowly.
2. The Speed Chef (r-process): Works in a frenzy, slamming neutrons together to make heavy stuff fast.

But there's a problem. Some rare, proton-rich ingredients (called p-nuclides) are like "forbidden fruits." They are shielded by stable elements, so neither the Slow Cook nor the Speed Chef can reach them. For decades, scientists have been puzzled: How does the universe make these rare ingredients?

The Suspect: The "Neutrino-P" Chef

The leading theory is a process called the $\nu$p-process (nu-p process). It happens inside a Core-Collapse Supernova—the massive explosion of a dying star.

Think of the supernova as a pressure cooker. Inside, a tiny, super-dense ball called a Proto-Neutron Star (PNS) is cooling down. It shoots out a massive stream of ghostly particles called neutrinos.

• These neutrinos hit the gas around the star.
• They turn some protons into neutrons.
• These neutrons act as "keys" that unlock the door to the rare p-nuclides, allowing the star to cook them up.

The Missing Ingredient: General Relativity

For a long time, scientists modeled this cooking process using Newtonian physics (the rules of gravity that work for apples falling from trees). But the core of a supernova is so dense and hot that Einstein's General Relativity (GR) should actually apply.

The Big Question: Does Einstein's gravity change the recipe?

In this paper, the authors (a team from SLAC and other universities) decided to cook the meal twice: once with Newton's rules and once with Einstein's rules, to see if the taste (the amount of rare elements produced) changes.

The Analogy: The Gravity Well and the Blue-Shifted Flashlight

To understand what the authors found, imagine the Proto-Neutron Star is a deep, dark gravity well (like a very deep canyon).

1. The Flashlight (Neutrinos): The star shines a flashlight (neutrinos) up out of the canyon.
2. The Blue Shift: As the light climbs out of the deep gravity well, it gets "blueshifted." This means the light gains energy, becoming more intense and "hotter" as it leaves the bottom.
3. The Newtonian Mistake: In the old Newtonian models, scientists assumed the flashlight had the same brightness at the bottom of the canyon as it did at the top. They missed the fact that the light gets super-charged as it escapes the deep gravity.
4. The Einstein Correction: The authors used Einstein's math. They realized that because the gravity is so strong, the neutrinos get a massive energy boost as they escape.

The Results: A Faster, More Efficient Kitchen

When the authors turned on the "Einstein switch," the results were dramatic:

• The Wind Blows Faster: The extra energy from the blueshifted neutrinos pushes the gas out of the star much faster.
• Less "Seed" Production: Because the gas is moving so fast, it doesn't spend enough time in the "seed-making" zone (where heavy nuclei are built). It rushes past it.
• The Good News: This is actually good for making the rare p-nuclides! The $\nu$p-process needs a specific ratio of protons to seeds. By rushing past the seed-making phase, the gas ends up with too many protons and too few seeds. This is the perfect recipe for the rare elements to form.

The Analogy: Imagine trying to bake a cake.

• Newtonian Model: You let the batter sit too long. It turns into a dense, heavy loaf (too many seeds, not enough rare elements).
• Einstein Model: You blast the oven with extra heat (GR effects). The batter rises so fast it skips the "dense loaf" stage and instantly transforms into the perfect, fluffy rare cake.

The Specific Wins

The paper highlights three major victories for the Einstein model:

1. The "Mo" and "Ru" Jackpot: The model successfully reproduces the exact amounts of Molybdenum (Mo) and Ruthenium (Ru) found in our solar system. The Newtonian model failed to make enough of these.
2. The "Nb" Miracle: There is a rare isotope called Niobium-92. In the Newtonian model, it was almost non-existent. In the Einstein model, the production of this element jumped by a factor of 25. This is huge because Niobium-92 is a "cosmic clock" used to date the age of the solar system.
3. The Mass Matters: The model works best for stars that are "just right" (around 18 times the mass of our Sun). If the star is too light, the explosion is too violent, and the gas flies away too fast (supersonic), ruining the recipe. If it's too heavy, the gravity is too strong. The 18-solar-mass star is the "Goldilocks" zone where Einstein's gravity makes the recipe perfect.

The Conclusion

The paper concludes that General Relativity is not just a tiny detail; it is a crucial ingredient.

Without Einstein's corrections, our understanding of how the universe makes rare elements is broken. With them, the $\nu$p-process in a supernova becomes a "unifying explanation" that can account for almost all the rare, proton-rich elements we see in the solar system today.

In short: The universe is a relativistic kitchen. To understand how it cooks the rarest ingredients, we must use Einstein's recipe, not Newton's.

Technical Summary

1. Problem Statement

The origin of proton-rich isotopes (p-nuclides) in the solar system, particularly those in the mass range $74 \le A \le 102$ (such as $^{92,94}\text{Mo}$, $^{96,98}\text{Ru}$, and the shielded isotope $^{92}\text{Nb}$), remains a significant astrophysical puzzle. The $\nu$p-process (neutrino-driven proton capture process) is a leading candidate mechanism, occurring in the neutrino-driven outflows (NDOs) of core-collapse supernovae (CCSNe).

While previous studies have analyzed various physical inputs (neutrino luminosities, entropy, progenitor mass), the impact of General Relativity (GR) on the $\nu$p-process has remained unexplored. Previous GR hydrodynamic studies often yielded conflicting results regarding entropy enhancements (some reporting increases of 30–40 units, others negligible) and relied on assumptions (e.g., fixed neutrino spectra at the neutrinosphere) that may not reflect modern simulation outputs. This paper aims to systematically quantify how GR effects influence the hydrodynamics of NDOs and, consequently, the integrated yields of p-nuclides.

2. Methodology

The authors developed a comprehensive framework combining fully relativistic hydrodynamics with time-dependent nucleosynthesis calculations.

• General Relativistic Hydrodynamics:

  • Derived and solved steady-state, spherically symmetric hydrodynamic equations for NDOs in Schwarzschild spacetime.
  • Formulated the equations using covariant variables to transparently identify relativistic corrections: gravitational redshift/blueshift, relativistic velocity corrections ($\mu$ factor), and the contribution of radiation enthalpy to inertia (radiation mass).
  • Implemented a Boundary Value Problem (BVP) approach. Instead of fixing initial conditions at the proto-neutron star (PNS) surface, they matched the outflow to the homologously expanding material behind the supernova shock front. This allows the flow regime (subsonic vs. transonic) to be determined self-consistently based on progenitor properties.
  • Included variable relativistic degrees of freedom (RDF) in the equation of state (EoS) to account for the transition of electron-positron pairs as temperature drops.
  • Applied fully relativistic Rankine-Hugoniot conditions for termination shocks when the flow becomes supersonic.

• Neutrino Physics Inputs:

  • Used neutrino luminosities and spectra predicted by modern CCSN simulations (reported at a far radius, $R_{\text{far}} = 500$ km) rather than at the neutrinosphere.
  • Applied gravitational blueshift to neutrino energies as they travel from the far radius to the PNS surface, and redshift for the reverse.
  • Included neutrino-electron scattering and weak magnetism/recoil corrections to determine the electron fraction ($Y_e$).

• Nucleosynthesis Calculation:

  • Constructed tracer particle trajectories by "gluing" the steady-state outflow solution to the homologous expansion phase.
  • Post-processed these trajectories using the open-source reaction network SkyNet.
  • Calculated time-integrated yields by summing contributions from mass elements launched at different times ($t_{\text{launch}} = 1$ to $8$ s) during the cooling phase, weighted by the mass outflow rate ($\dot{M}$).

• Models:

  • Benchmark: $18 M_\odot$ progenitor with a $1.8 M_\odot$ PNS.
  • Variations: $12.75 M_\odot$ and $9 M_\odot$ progenitors with varying front shock velocities ($v_{\text{FS}}$) to test sensitivity to flow regimes (subsonic vs. transonic).

3. Key Contributions

1. Resolution of Entropy Discrepancies: The paper clarifies why previous studies reported large GR-induced entropy enhancements. The authors demonstrate that these large values ($\Delta S \sim 30$) arise from prescribing identical neutrino spectra at the neutrinosphere for both GR and Newtonian cases. By prescribing spectra at a far radius (consistent with simulations) and applying blueshift only to the GR case, the entropy enhancement is modest ($\Delta S \lesssim 5$).
2. Identification of Dominant GR Mechanisms: The study isolates two primary GR effects:
   • Neutrino Blueshift: Increases neutrino energy near the PNS, boosting heating rates and accelerating the outflow.
   • Effective Gravitational Potential: Deepens the potential well, affecting the entropy profile.
   • The authors show that the acceleration caused by blueshift is the dominant factor, reducing the time matter spends in the seed-formation window.
3. Self-Consistent Flow Regime Determination: Unlike prior works that assumed subsonic or transonic flows, this study solves for the flow regime dynamically. It reveals that GR effects can push near-critical flows into the supersonic regime, forming termination shocks, which drastically alters nucleosynthesis outcomes.

4. Key Results

Hydrodynamic Effects

• Subsonic Flows ($18 M_\odot$): GR corrections lead to a modest entropy increase ($\sim 2$ units) but a significant increase in outflow velocity (factor of $\sim 3$) due to blueshift. This reduces the time spent in the seed-formation temperature window ($6 \text{ GK} > T > 3 \text{ GK}$), suppressing seed nuclei production ($^{56}\text{Ni}$).
• Transonic Flows ($12.75 M_\odot$): For lighter progenitors or higher shock velocities, GR effects can trigger a transition to supersonic flow earlier than in Newtonian cases. This leads to the formation of termination shocks, which suppresses the $\nu$p-process efficiency.

Nucleosynthesis Yields

• Enhanced p-Nuclide Production: In the subsonic $18 M_\odot$ model, GR effects significantly boost the production of p-nuclides.
  • $^{92,94}\text{Mo}$ and $^{96,98}\text{Ru}$: Production is enhanced by factors of 4 to 9 compared to Newtonian calculations.
  • $^{102}\text{Pd}$: Enhanced by a factor of $\sim 30$.
  • Mechanism: The suppression of seed nuclei (due to faster expansion) increases the proton-to-seed ratio ($Y_p/Y_{\text{seed}}$). Combined with a favorable ratio of neutron production time to seed formation time ($\tau_2/\tau_1$), this allows the $\nu$p-process to bypass waiting points more efficiently.
• The Case of $^{92}\text{Nb}$:
  • $^{92}\text{Nb}$ is a shielded isotope produced in Stage III (late-time neutron captures, $T < 1.5 \text{ GK}$).
  • GR effects boost $^{92}\text{Nb}$ production by a factor of 25.
  • This is because GR-driven acceleration leads to higher proton-to-seed ratios at late times, and the optimal production window for $^{92}\text{Nb}$ occurs later ($t \sim 5$ s) when GR effects on the contracting PNS are strongest.
• Solar Abundance Reproduction: The GR model for the $18 M_\odot$ progenitor successfully reproduces both the relative and absolute solar system abundances for the entire set of p-nuclides in the range $74 \le A \le 102$. The Newtonian model fails to produce sufficient absolute yields for these isotopes.

Dependence on Progenitor Mass

• $9 M_\odot$ Progenitor: Outflows are supersonic almost immediately. GR corrections enhance yields but not enough to match solar observations; the $\nu$p-process is inefficient in very light progenitors.
• $12.75 M_\odot$ Progenitor: The outcome is highly sensitive to the front shock velocity ($v_{\text{FS}}$).
  • At $v_{\text{FS}} = 6000$ km/s, flows remain subsonic, and yields are marginally compatible with observations.
  • At $v_{\text{FS}} = 8000$ km/s, GR induces a supersonic transition at $t \sim 2.3$ s, causing a sharp drop in yields, making them incompatible with solar data.

5. Significance

• Unifying Explanation: The paper provides a robust, unified explanation for the origin of all solar system p-nuclides up to $^{102}\text{Pd}$, provided the progenitor is sufficiently massive ($\sim 18 M_\odot$) and GR effects are included.
• Critical Role of GR: It establishes that General Relativity is not a minor correction but a decisive factor in determining the efficiency of the $\nu$p-process. The interplay between neutrino blueshift and gravitational potential fundamentally alters the thermodynamic history of the ejecta.
• Methodological Advancement: By resolving the discrepancy in entropy calculations and implementing a self-consistent BVP with variable RDF and realistic neutrino inputs, the authors set a new standard for modeling neutrino-driven outflows.
• Cosmochronometry: The successful reproduction of $^{92}\text{Nb}$ yields (a key cosmochronometer) validates the $\nu$p-process as the origin of this extinct radionuclide, resolving a long-standing debate about its production mechanism.

In conclusion, the study demonstrates that General Relativistic effects are essential for the $\nu$p-process to function as a viable source of heavy p-nuclides in the universe, particularly by suppressing seed formation and enhancing the neutron-to-seed ratio in massive progenitor explosions.