O–C Shell Merger Nucleosynthesis

• O–C shell merger nucleosynthesis is a process where the convective oxygen-burning shell merges with the carbon-burning shell, leading to convective-reactive burning and altered nucleosynthetic pathways.
• 3D hydrodynamic simulations reveal rapid mixing and high convective velocities that significantly boost the production of odd-Z elements, p-nuclei, and radioactive isotopes in the final pre-supernova stages.
• Uncertainties in turbulent entrainment, nuclear reaction rates, and mixing physics demand advanced simulations for accurate predictions of stellar yields and their impact on Galactic chemical evolution.

Oxygen–carbon (O–C) shell merger nucleosynthesis refers to the suite of nuclear transmutations, mixing physics, and resulting elemental yields triggered when the convective O-burning shell in an evolved massive star merges with the neighboring C-burning shell late in stellar evolution. This phenomenon radically alters the internal structure and composition in the final hours to years before core collapse, with profound implications for the origin of odd-Z elements, p-nuclei, radioactive isotopes, and pre-supernova asphericity. O–C shell mergers are now recognized as critical drivers of chemical enrichment in massive stars, with unique signatures observable in supernova remnants and the Galactic chemical inventory.

1. Physical Conditions and Onset of O–C Shell Mergers

O–C shell mergers occur when the convective boundaries of the inner O-burning region (T ≈ 2–3 GK, ρ ≈ 10^5–10^6 g cm⁻³) and the overlying C-burning shell (T ≈ 1–1.5 GK) are eroded by turbulent entrainment and reduced entropy gradients, often accelerated by late stage core contraction (Roberti et al., 26 Apr 2025, Rizzuti et al., 2024, Andrassy et al., 2018). Three-dimensional hydrodynamic simulations confirm that non-spherical entrainment across convective boundaries (entrainment rate $$\dot{M}_\text{ent} \approx A_{CB} \rho_{CB} v_{\text{conv,CB}}$$) merges the two shells, forming a single, large, vigorously-convective zone with turnover times $\tau_{c} \sim 10^{2}$–$10^{3}$ s—substantially shorter, and with velocities 5–30× higher than predicted by 1D MLT (Rizzuti et al., 2024, Andrassy et al., 2018, Issa et al., 23 Sep 2025).

The physical conditions in the merged shell are highly dynamic:

• Temperatures rise to $T_\mathrm{base} \sim 2.6$–$2.85$ GK.
• Entrained C and Ne from the disrupted interface are advected deep into the O-burning zone.
• Convective velocities can reach $v_\text{conv} \sim 5 \times 10^7$ cm s⁻¹ immediately following merger (Rizzuti et al., 2024).
• Multiple burning phases (O, Ne, C, and sometimes He) operate concurrently due to multi-scale mixing (Rizzuti et al., 2024, Sato et al., 10 Jul 2025).

2. Nucleosynthetic Pathways in Merged Shells

A defining feature of O–C shell mergers is “convective-reactive” burning: the simultaneous ingestion of unburned fuel into a high-entropy convective shell, where nuclear timescales become comparable to mixing timescales ($\tau_\text{nuc} \sim \tau_\text{mix}$) (Ritter et al., 2017). Dominant nuclear flows include:

• $\alpha$- and proton-capture chains: $^{16}$O$(\alpha,\gamma)^{20}$Ne, $^{20}$Ne$(\alpha,\gamma)^{24}$Mg, $^{24}$Mg$(\alpha,\gamma)^{28}$Si.
• Heavy-ion fusions: $^{12}$C+$^{12}$C$\to^{20}$Ne$+\alpha$, $^{12}$C+$^{16}$O$\to^{28}$Si$+\alpha$.
• The “SPAr process” (proton captures on S, P, Ar): e.g., $^{32}$S$(p,\gamma)^{33}$Cl, $^{34}$S$(p,\gamma)^{35}$Cl, $^{38}$Ar$(p,\gamma)^{39}$K, $^{31}$P$(p,\alpha)^{28}$Si, with energy output dominating the nuclear budget by factors $\gtrsim400$ over C+O fusion for $X_p\gtrsim10^{-9}$ (Roberti et al., 17 Sep 2025).

These flows synthesize large abundances of intermediate-mass $\alpha$-nuclei (Mg, Si, S), “odd-Z” species (P, Cl, K, Sc), and light radioactive isotopes (e.g., $^{40}$K, $^{44}$Ti, $^{41}$Ca) (Ritter et al., 2017, Issa et al., 23 Sep 2025, Issa et al., 19 Dec 2025).

3. Quantitative Yield Predictions and Uncertainties

The yield of key isotopes is highly sensitive to mixing physics and ingestion rates, as confirmed by both 1D post-processing and 3D hydro-informed models (Issa et al., 22 Jul 2025, Issa et al., 23 Sep 2025, Issa et al., 19 Dec 2025). Examples include:

Isotope	Pre-SN Yield Range (M$_\odot$; 15 M$_\odot$ model, various mixing)	Overproduction Factor (1D/3D)
$^{44}$Ti	$1.94\times10^{-7}$—$6.69\times10^{-3}$ (4.8 dex spread) (Issa et al., 19 Dec 2025)	up to $>1.5$ dex
$^{40}$K	$>10^3$ variation across 3D mixing regimes (Issa et al., 23 Sep 2025)	up to 3 dex
$^{39,41}$K, Sc	up to $10$–$100\times$ enhancement relative to standard O-shell yields (Ritter et al., 2017, Roberti et al., 26 Apr 2025)	$>1$ dex
p-nuclei ($A>100$)	$F_i \sim 10$–$40\times$ solar (e.g., $^{130,132}$Ba, $^{144}$Sm) (Roberti et al., 2024)	up to 1.5 dex

Yields for p-nuclei and odd-Z elements are non-monotonic in mixing speed and ingestion rate, with non-linearities arising from feedback between burning and flow (Issa et al., 22 Jul 2025). Increased ingestion of C/Ne enhances α and p densities, strongly boosting heavy-element production via photodisintegration and $(p,\gamma)$ or $(\alpha,\gamma)$ chains (Roberti et al., 26 Apr 2025, Roberti et al., 17 Sep 2025).

Macro-physics uncertainties (boundary downturn, quenching, velocity boosts) produce yield variations for odd-Z and p-nuclei on par with those from order-of-magnitude uncertainties in experimental photo-disintegration rates. The sensitivity of $^{44}$Ti yield to diffusion coefficient $D$ is super-linear, with $\partial\log Y/\partial \log D \sim2$–$4$ (Issa et al., 19 Dec 2025).

4. Observational Evidence and Remnant Signatures

O–C shell merger nucleosynthetic outputs are directly observed in several contexts:

• Cassiopeia A SNR: Inhomogeneous mixing signatures (Ne-rich downflows, Si-rich upflows) and spatially segregated O–Ne–rich versus O–Si–rich ejecta, confirming rapid, multi-scale mixing hours before core collapse (Sato et al., 10 Jul 2025).
• SNR G359.0–0.9: X-ray spectra exhibiting anomalously high Mg/Ne and elevated Si/Mg, consistent with Ne-shell intrusion and O–C shell merger models for progenitors of $>15$ M$_\odot$ (Matsunaga et al., 2024).
• EMP star abundances: Galactic [K/Mg] and [Sc/Mg] ratios reproduced only if O–C shell merger yields are included, with merger models populating the observed [K/Mg] > 0, [Sc/Mg] > 0 domain unreachable by standard nucleosynthesis and neutrino-processed models (Roberti et al., 26 Apr 2025, Ritter et al., 2017).

Gamma-ray line observations of $^{44}$Ti and $^{22}$Na in young SNRs yield total masses and morphologies best explained by pre-explosive merger yields in 3D-informed models, rather than by purely explosive nucleosynthesis (Issa et al., 19 Dec 2025).

5. Impact on Galactic Chemical Evolution and Astrophysical Consequences

O–C shell mergers solve the longstanding underproduction problem for P, Cl, K, Sc, and certain p-nuclei in Galactic chemical evolution models (Ritter et al., 2017, Roberti et al., 26 Apr 2025). Their yields, when included at even ~10–50% frequency among massive stars, elevate [K/Fe], [Sc/Fe], and [P/Fe] to match Milky Way trends; a higher merger fraction is needed for solar-like Cl/Fe. The merger regime also provides a robust p-nuclei production floor independent of explosion energy, with pre-explosive γ-process dominating for $A\gtrsim 100$ nuclei in O–C merger progenitors (Roberti et al., 2024).

The process injects large, spatially inhomogeneous quantities of $^{40}$K, with implications for early radiogenic heating and geodynamic variability in rocky exoplanets (Issa et al., 23 Sep 2025). $^{44}$Ti produced pre-explosively can dominate the total yield over any explosive channel for plausible 3D mixing and ingestion rates, affecting SNR light curves, radiogenic $^{44}$Ca in presolar grains, and observable $\gamma$-line fluxes (Issa et al., 19 Dec 2025).

The merger’s altered pre-supernova density profile (lower compactness) and large-scale compositional asymmetry may promote more energetic, asymmetric core-collapse supernovae, supporting enhanced neutron star kicks and facilitating shock revival in 3D explosion simulations (Rizzuti et al., 2024, Sato et al., 10 Jul 2025).

6. Theoretical Uncertainties, Open Problems, and Future Directions

The primary theoretical uncertainties in O–C shell merger nucleosynthesis stem from:

• The fidelity of 1D prescriptions for turbulent convective mixing, which underestimate both entrainment rates and flow velocities relative to 3D simulations (Andrassy et al., 2018, Rizzuti et al., 2024, Issa et al., 22 Jul 2025).
• Incomplete knowledge of the rate and global occurrence of O–C mergers in stellar populations. Observational constraints from SNRs and EMP stars strongly support a non-negligible frequency, but detailed statistics are model-dependent (Roberti et al., 26 Apr 2025).
• Nuclear physics uncertainties, particularly for proton and α-capture rates on S, P, Ar, and key photo-disintegration rates relevant to the p-process. These can induce $0.5$–$0.8$ dex variation in p-nuclei yields, comparable to mixing uncertainties (Issa et al., 22 Jul 2025, Roberti et al., 17 Sep 2025).
• Multi-dimensional effects on large-scale mixing, asphericity development, and feedback between burning and hydrodynamic flow remain incompletely captured, especially for dynamic shell mergers with global oscillations at high ingestion/burning rates (Andrassy et al., 2018).

Next-generation 3D hydrodynamic simulations with expanded nuclear networks and resolved boundary physics are required to calibrate mixing models, constrain nucleosynthetic yields, and enable accurate population-synthesis and GCE predictions (Rizzuti et al., 2024, Issa et al., 23 Sep 2025). Upcoming high-resolution X-ray and $\gamma$-ray spectroscopy (XRISM, Athena, COSI) will further test model predictions through abundance ratios, isotopic signatures, and the spatial structure of SNR ejecta (Matsunaga et al., 2024, Sato et al., 10 Jul 2025).

7. O-C Shell Mergers in He-CO White Dwarf Mergers

An analogous process occurs in post-merger remnants of He-CO white dwarf mergers, where the so-called shell-of-fire (SOF) forms under dynamic conditions ($T_\mathrm{SOF}\sim 1.2$–$2.4 \times 10^8$ K, $\rho_\mathrm{SOF}\sim 3$–$5\times10^4$ g cm$^{-3}$) (Menon et al., 2012). Key reactions include $^{14}\mathrm{N}(\alpha,\gamma)^{18}$F$(\beta^+\nu){}^{18}$O, $^{13}\mathrm{C}(\alpha,n)^{16}$O, and $^{18}$O$(\alpha,\gamma)^{22}$Ne. The evolution of $^{16}$O/$^{18}$O is highly sensitive to the duration and temperature of the SOF, with efficient envelope mixing required to match the extremely low $^{16}$O/$^{18}$O$\sim$1–10 observed in RCB stars. Post-merger rotational mixing and its subsequent decay $\sim10^6$ yr after merger determine final surface abundances (Menon et al., 2012).

---

In summary, O–C shell merger nucleosynthesis is a convective-reactive phenomenon essential for quantitative understanding of odd-Z species, p-nuclei, and key radioactive isotopes in massive stars. Accurately modeling these events requires both multi-dimensional hydrodynamic calibration and comprehensive nuclear networks. Their signatures are evident in both Galactic chemical trends and supernova remnant spectroscopy, and their macro-physical uncertainties now represent a limiting factor in predictive stellar nucleosynthesis across cosmic history (Issa et al., 19 Dec 2025, Roberti et al., 26 Apr 2025, Ritter et al., 2017, Issa et al., 22 Jul 2025).