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$

$c_s$ is the local sound speed, and $\delta f_{1,2}$ encode extra momentum deposition by gravito-acoustic waves.

• Entropy Evolution:

$$\frac{ds}{dr} = \frac{r}{v}\frac{\dot{q}_{\text{tot}}}{e^\Lambda W T}$$

with $\dot{q}_{\text{tot}} = \dot{q}_\nu + \dot{q}_w$—the sum of neutrino and wave heating per unit mass.

• Electron Fraction Evolution:

$$\frac{dY_e}{dr} = \frac{r}{v}\frac{\dot{Y}_e}{e^\Lambda W}$$

$$\dot{Y}_e = \lambda_{\nu_e}(1-Y_e) - \lambda_{\bar{\nu}_e}Y_e$$

where $\lambda_{\nu_e}$, $\lambda_{\bar{\nu}_e}$ are absorption rates.

• Wave Action Evolution (gravito-acoustic waves):

$$\frac{dS}{dr} = -S\left(\frac{2}{r} + \frac{1}{l_d} + \frac{1}{v_g}\frac{dv_g}{dr}\right)$$

$S(r) = \frac{L_w}{4\pi r^2 c_s \omega}$ encapsulates the local wave energy. The dissipation length $l_d$ specifies over which scale waves shock.

Typical parameters include $M_{NS}=1.4-2.1\,M_\odot$, $L_\nu=3\times10^{51}-1.2\times10^{53}$ erg/s, $L_w=\epsilon_w L_\nu$ with $\epsilon_w=10^{-5}-10^{-2}$, $\omega=10^2-10^4$ s$^{-1}$, $R_{NS}\approx 12$ km.

2. Boundary Conditions and Solution Techniques

Solutions are integrated from the neutrinosphere ($r\approx R_{NS}$, $T=3-4$ MeV, $\rho \sim 10^{11}-10^{12}$ g/cm$^3$), imposing boundary values for $L_\nu$, mean neutrino energy, and equilibrium $Y_{e}$. Wave luminosity is fixed as a fraction of $L_\nu$.

A shooting method is employed: an initial guess for mass-flux $\dot M$ is iteratively refined, integrating the ODE system through the sonic point (critical point where $v=c_s$), 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 $<10^{-8}$.

3. Wind Thermodynamics, Regime Classification, and Wave Effects

The inclusion of convection-driven wave luminosity $\epsilon_w$ reorganizes wind dynamics into three distinct regimes (Nevins et al., 2024):

• Regime I ($\epsilon_w\lesssim 2\times10^{-4}$):

Mild wind acceleration, modest entropy enhancement ($\Delta s \sim 10$–20$k_B$), $Y_e$ at equilibrium, expansion timescale shortened, resulting in enhanced $\nu p$-process nucleosynthesis up to $A\sim100$-140.

• Regime II ($2\times10^{-4}\lesssim \epsilon_w \lesssim 10^{-3}$):

Early acceleration reduces exposure to neutrino heating, lowers $Y_e$ and entropy; seed production increases—impeding $\nu p$-process and stifling nucleosynthesis near the iron peak.

• Regime III ($\epsilon_w \gtrsim 10^{-3}$):

Shocks form at small radii ($r_s\lesssim30$ km), injecting heat ($\sim$1 MeV/baryon), entropy rises above 100$k_B$, very rapid outflow ($\tau<1$ ms), $\alpha$ recombination disrupted. An $(n,\gamma)$-driven “fast-outflow r-process” commences, proceeding up to $A\sim200$ despite equilibrium $Y_e=0.6$.

The wind response is strongly nonmonotonic: a dip in $s^2/\tau$ and maximum seed formation at intermediate $\epsilon_w$ suppresses heavy-element yields; higher $\epsilon_w$ correlates with heavier $\nu p$-processing up to $A\sim200$.

4. Nuclear Reaction Network and Nucleosynthetic Outcomes

Post-processing employs a large network (SkyNet, $\sim$8000 isotopes) spanning strong/electromagnetic (n,γ), (p,γ), (α,γ), (α,n), (α,p), weak (β$^\pm$, e$^\pm$ capture), and neutrino-induced channels (notably $p(\bar{\nu}_e,e^+)n$ and $n(\nu_e,e^-)p$). Fission for $A>240$ is included.

Key reaction rates:

• Triple-$\alpha$/α-capture rates set seed formation.
• $p(\bar{\nu}_e,e^+)n$ determines free-neutron availability for the νp-process (at $T\sim1.5$–3 GK).

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

$$(n/\text{seed}) \propto Y_p n_\nu \exp(-\tau_{\exp}/\cdots) / Y_{\text{seed}}$$

where $n_\nu = \int \lambda_{\bar{\nu}_e} dt$ from $T=3$ GK down.

Results:

• For $\epsilon_w \lesssim 2\times10^{-4}$, classic νp-process signatures peak at $A\sim90$–120, endpoint $A_{max}$ correlated with $s^2/\tau$.
• For $\epsilon_w \gtrsim 10^{-3}$, shock heating yields a suppressed, r-process-like pattern with peaks near $A\sim130$ and 200, albeit lower abundances than full solar r-process.

Modulating the wind termination radius $r_{ts}$ affects nucleosynthetic yields: in the $\nu p$-process regime, a smaller $r_{ts}$ prolongs high $T$ 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 $S \propto P_\star^{-5/6}$, with favorable conditions for third-peak r-process ($S\gtrsim150-200$, $\tau_{\exp}\lesssim5\times10^{-4}$ s, $Y_e\lesssim0.48$). For $B_0=3\times10^{15}$ G, $M_r\sim1$–$5\times10^{-5} M_\odot$ synthesized in $\sim1$–$2$ s.

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

Rapid rotation focuses outflows equatorially, increases mass-loss rates by $>10\times$ but lowers entropy and $Y_e$, 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 $L_w/L_\nu \sim 10^{-5}$. Three distinct yield regimes are established as $\epsilon_w$ rises: extended νp-processing, seed-dominated suppression, and shock-driven, fast-outflow r-process. The transition points ($\epsilon_w \approx 2\times10^{-4}$, $10^{-3}$) are robust under varying $M_{NS}$ and $L_\nu$.

Proto-neutron star convection should excite gravity waves with $L_w/L_\nu \sim 10^{-5}$–$10^{-2}$ and $\omega \sim 10^2$–$10^4$ s$^{-1}$ (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

$\epsilon_w (=L_w/L_\nu)$	Dominant Process	Entropy $s$ ($k_B$)	Expansion Time $\tau$ (ms)	Nucleosynthetic Endpoint $A_{max}$	Notes
$\lesssim 2\times10^{-4}$	Enhanced νp-process	$80$–$100$	$5$–$10$	$100$–$140$	Higher $M_{NS}$ shifts $A_{max}$ up
$2\times10^{-4}–10^{-3}$	Seed overproduction	$80$–$100$	$5$–$10$	$60$–$90$	νp-process stifled, iron-peak
$\gtrsim 10^{-3}$	Shock-driven r-process	$>100$	$<1$	$130$–$200$	Early shock, "fast r-process"

The termination radius $r_{ts}$ and wind microphysics further modulate yields, especially at high $\epsilon_w$.

• "Proto-Neutron Star Convection and the Neutrino-Driven Wind: Implications for the $ν$p-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)