Proto-Neutron Star Wind Models

• Proto-neutron star wind models describe quasi-steady, neutrino-heated outflows that drive nucleosynthesis in supernovae.
• They incorporate advanced general relativistic equations and detailed microphysics to predict mass-loss rates, entropy, and electron fraction conditions.
• Wave-induced shock heating modulates nucleosynthetic yields, shifting outcomes between νp-process enhancement and fast-outflow r-process signatures.

A proto-neutron star (PNS) wind model describes the quasi-steady, mass-loaded outflow driven by intense neutrino heating in the seconds following core collapse. This wind phase is central to the theory of nucleosynthesis in supernovae, as it sets the physical conditions for the synthesis of trans-iron nuclei by the νp-process, the α-process, and potentially the r-process. Modern wind models incorporate general relativistic effects, sophisticated treatments of microphysics (neutrino interactions, equation of state, charged-current rates), convection-driven instabilities, rotation and magnetization, and secondary energy deposition via gravito-acoustic waves. These ingredients combine to regulate key diagnostic quantities—mass-loss rate, entropy per baryon, electron fraction, expansion timescale—that control the assembled abundances of heavy nuclei.

1. Steady-State General Relativistic Wind Equations

PNS wind models are governed by a set of coupled differential equations describing the outflow in spherical symmetry and full general relativity. The line element is

$$ds^2 = -e^{2\Lambda(r)}c^2dt^2 + e^{-2\Lambda(r)}dr^2 + r^2d\Omega^2$$

with

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$

for a neutron star mass $M_{NS}$. The key equations (Nevins et al., 2024) are:

• Mass Conservation:

$\dot{M} = 4\pi r^2 e^\Lambda W \rho v$

where $v(r)$ is the radial velocity, $\rho(r)$ the mass density, and $W \equiv 1/\sqrt{1-v^2/c^2}$ the Lorentz factor.

• Momentum Conservation (including wave stresses):

$\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$

$f_1 = [1 - v^2/c_s^2] + \delta f_1$

$f_2 = \text{(standard wind terms)} + \delta f_2$

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$0 is the local sound speed, and $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$1 encode extra momentum deposition by gravito-acoustic waves.

• Entropy Evolution:

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$2

with $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$3—the sum of neutrino and wave heating per unit mass.

• Electron Fraction Evolution:

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$4

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$5

where $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$6, $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$7 are absorption rates.

• Wave Action Evolution (gravito-acoustic waves):

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$8

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$9 encapsulates the local wave energy. The dissipation length $M_{NS}$0 specifies over which scale waves shock.

Typical parameters include $M_{NS}$1, $M_{NS}$2 erg/s, $M_{NS}$3 with $M_{NS}$4, $M_{NS}$5 s$M_{NS}$6, $M_{NS}$7 km.

2. Boundary Conditions and Solution Techniques

Solutions are integrated from the neutrinosphere ($M_{NS}$8, $M_{NS}$9 MeV, $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$0 g/cm$\dot{M} = 4\pi r^2 e^\Lambda W \rho v$1), imposing boundary values for $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$2, mean neutrino energy, and equilibrium $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$3. Wave luminosity is fixed as a fraction of $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$4.

A shooting method is employed: an initial guess for mass-flux $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$5 is iteratively refined, integrating the ODE system through the sonic point (critical point where $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$6), enforcing transonic regularity to machine precision. High-resolution grids (200–500 log-spaced radial zones) are necessary; the system is converged when the critical condition is met to $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$7.

3. Wind Thermodynamics, Regime Classification, and Wave Effects

The inclusion of convection-driven wave luminosity $\dot{M} = 4\pi r^2 e^\Lambda W \rho v$8 reorganizes wind dynamics into three distinct regimes (Nevins et al., 2024):

• Regime I ($\dot{M} = 4\pi r^2 e^\Lambda W \rho v$9):

Mild wind acceleration, modest entropy enhancement ($v(r)$0–20$v(r)$1), $v(r)$2 at equilibrium, expansion timescale shortened, resulting in enhanced $v(r)$3-process nucleosynthesis up to $v(r)$4-140.

• Regime II ($v(r)$5):

Early acceleration reduces exposure to neutrino heating, lowers $v(r)$6 and entropy; seed production increases—impeding $v(r)$7-process and stifling nucleosynthesis near the iron peak.

• Regime III ($v(r)$8):

Shocks form at small radii ($v(r)$9 km), injecting heat ($\rho(r)$01 MeV/baryon), entropy rises above 100$\rho(r)$1, very rapid outflow ($\rho(r)$2 ms), $\rho(r)$3 recombination disrupted. An $\rho(r)$4-driven “fast-outflow r-process” commences, proceeding up to $\rho(r)$5 despite equilibrium $\rho(r)$6.

The wind response is strongly nonmonotonic: a dip in $\rho(r)$7 and maximum seed formation at intermediate $\rho(r)$8 suppresses heavy-element yields; higher $\rho(r)$9 correlates with heavier $W \equiv 1/\sqrt{1-v^2/c^2}$0-processing up to $W \equiv 1/\sqrt{1-v^2/c^2}$1.

4. Nuclear Reaction Network and Nucleosynthetic Outcomes

Post-processing employs a large network (SkyNet, $W \equiv 1/\sqrt{1-v^2/c^2}$28000 isotopes) spanning strong/electromagnetic (n,γ), (p,γ), (α,γ), (α,n), (α,p), weak (β$W \equiv 1/\sqrt{1-v^2/c^2}$3, e$W \equiv 1/\sqrt{1-v^2/c^2}$4 capture), and neutrino-induced channels (notably $W \equiv 1/\sqrt{1-v^2/c^2}$5 and $W \equiv 1/\sqrt{1-v^2/c^2}$6). Fission for $W \equiv 1/\sqrt{1-v^2/c^2}$7 is included.

Key reaction rates:

• Triple-$W \equiv 1/\sqrt{1-v^2/c^2}$8/α-capture rates set seed formation.
• $W \equiv 1/\sqrt{1-v^2/c^2}$9 determines free-neutron availability for the νp-process (at $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$0–3 GK).

The main diagnostic is the neutron-to-seed ratio:

$\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$1

where $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$2 from $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$3 GK down.

Results:

• For $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$4, classic νp-process signatures peak at $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$5–120, endpoint $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$6 correlated with $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$7.
• For $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$8, shock heating yields a suppressed, r-process-like pattern with peaks near $\frac{dv}{dr} = \frac{v}{r}\frac{f_2}{f_1}$9 and 200, albeit lower abundances than full solar r-process.

Modulating the wind termination radius $f_1 = [1 - v^2/c_s^2] + \delta f_1$0 affects nucleosynthetic yields: in the $f_1 = [1 - v^2/c_s^2] + \delta f_1$1-process regime, a smaller $f_1 = [1 - v^2/c_s^2] + \delta f_1$2 prolongs high $f_1 = [1 - v^2/c_s^2] + \delta f_1$3 exposure and increases heavy-element output.

5. Comparative Model Context: Magnetized and Rotating Winds

Proto-neutron star winds under rapid rotation or strong magnetization further modify nucleosynthetic regimes.

• 2D MHD Models (Prasanna et al., 2024):

Magnetized, rapidly rotating winds eject high-entropy plasmoids quasi-periodically. The maximum entropy $f_1 = [1 - v^2/c_s^2] + \delta f_1$4, with favorable conditions for third-peak r-process ($f_1 = [1 - v^2/c_s^2] + \delta f_1$5, $f_1 = [1 - v^2/c_s^2] + \delta f_1$6 s, $f_1 = [1 - v^2/c_s^2] + \delta f_1$7). For $f_1 = [1 - v^2/c_s^2] + \delta f_1$8 G, $f_1 = [1 - v^2/c_s^2] + \delta f_1$9–$f_2 = \text{(standard wind terms)} + \delta f_2$0 synthesized in $f_2 = \text{(standard wind terms)} + \delta f_2$1–$f_2 = \text{(standard wind terms)} + \delta f_2$2 s.

• Non-rotating vs. Rotating Models (Desai et al., 2022):

Rapid rotation focuses outflows equatorially, increases mass-loss rates by $f_2 = \text{(standard wind terms)} + \delta f_2$3 but lowers entropy and $f_2 = \text{(standard wind terms)} + \delta f_2$4, suppressing heavy r-process, but possibly powering light neutron-rich element production (LEPP).

• Nucleosynthetic Impact:

The occurrence rate, field strength, and birth spin of magnetars fundamentally constrain their contribution to Galactic r-process inventories (Vincenzo et al., 2021).

6. Astrophysical Implications and Future Directions

Wave effects—convection-driven gravito-acoustic fluxes—alter NDW nucleosynthetic endpoints even at $f_2 = \text{(standard wind terms)} + \delta f_2$5. Three distinct yield regimes are established as $f_2 = \text{(standard wind terms)} + \delta f_2$6 rises: extended νp-processing, seed-dominated suppression, and shock-driven, fast-outflow r-process. The transition points ($f_2 = \text{(standard wind terms)} + \delta f_2$7, $f_2 = \text{(standard wind terms)} + \delta f_2$8) are robust under varying $f_2 = \text{(standard wind terms)} + \delta f_2$9 and $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$00.

Proto-neutron star convection should excite gravity waves with $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$01–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$02 and $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$03–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$04 s$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$05 (Nevins et al., 2024). Consequently, realistic NDW models must self-consistently integrate these effects to accurately predict p-nuclei and r-process contributions. Observational comparisons—meteoritic isotopic ratios, Galactic chemical evolution—require multi-dimensional simulations coupling time-dependent convection and wave transport.

7. Summary Table: Wind Regimes and Nucleosynthetic Outcomes

$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$06	Dominant Process	Entropy $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$07 ($e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$08)	Expansion Time $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$09 (ms)	Nucleosynthetic Endpoint $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$10	Notes
$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$11	Enhanced νp-process	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$12–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$13	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$14–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$15	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$16–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$17	Higher $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$18 shifts $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$19 up
$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$20	Seed overproduction	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$21–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$22	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$23–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$24	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$25–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$26	νp-process stifled, iron-peak
$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$27	Shock-driven r-process	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$28	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$29	$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$30–$e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$31	Early shock, "fast r-process"

The termination radius $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$32 and wind microphysics further modulate yields, especially at high $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$33.

• "Proto-Neutron Star Convection and the Neutrino-Driven Wind: Implications for the $e^\Lambda \equiv \sqrt{1 - 2GM_{NS}/(rc^2)}$34p-Process" (Nevins et al., 2024)
• "Favorable conditions for heavy element nucleosynthesis in rotating proto-magnetar winds" (Prasanna et al., 2024)
• "Three-Dimensional General-Relativistic Simulations of Neutrino-Driven Winds from Rotating Proto-Neutron Stars" (Desai et al., 2022)
• "Nucleosynthesis signatures of neutrino-driven winds from proto-neutron stars: a perspective from chemical evolution models" (Vincenzo et al., 2021)