---
title: O–C Shell Mergers in Massive Stars
url: https://www.emergentmind.com/topics/o-c-shell-mergers
type: topic
---

# O–C Shell Mergers in Massive Stars

O–C shell mergers, also termed C–O shell mergers and, in some contexts, O–Ne–C shell interactions, are late-stage convective-reactive events in massive stars in which the convective oxygen-burning shell ingests material from adjacent carbon- and neon-rich layers and the two shells effectively merge into a single extended mixed region. They occur after central He burning, during advanced C/Ne/O shell burning, typically hours to days before core collapse, and they can strongly modify pre-supernova structure, odd-\(Z\) nucleosynthesis, p-nucleus production, and the radioactive inventory ejected by the subsequent core-collapse supernova [2504.18867].

## 1. Definition and stellar setting

In the relevant massive-star regime, successive H, He, C, Ne, O, and Si burning stages build an onion-like interior consisting of an Fe core surrounded by distinct burning shells. Near collapse, the star may contain an inner O-burning shell at temperatures \(T \sim 2.6{-}2.8\) GK and an overlying C-burning shell at \(T \sim 1\) GK. An O–C shell merger occurs when the convective O-burning shell penetrates into the overlying C- and Ne-rich layers, ingests that material, and forms a large mixed convective region in which C-rich fuel is exposed to O-burning temperatures [2504.18867].

The physical trigger is weakening of the entropy and composition barrier between the shells. Once that barrier is sufficiently eroded, convection in the O shell entrains C/Ne-rich material upward and downward across the interface, and burning and mixing become strongly coupled. In the broader literature, these events are variously described as “C–O shell mergers,” “O–C shell mergers,” and “O–Ne–C shell interactions,” reflecting the fact that the O shell often interacts not only with C-rich material but also with Ne-rich layers in the same region [2509.13749].

The phenomenon is intrinsically convective-reactive. A standard diagnostic is the Damköhler number,
\[
Da=\frac{\tau_{\mathrm{mix}}}{\tau_{\mathrm{reac}}},
\]
with \(Da \sim 1\) indicating that convective transport and nuclear burning operate on comparable timescales. In this regime, the assumption of passive diffusive mixing is no longer adequate, because the flow changes the burning and the burning feeds back on the flow [1704.05985].

## 2. Occurrence and early diagnostics

A systematic survey of 209 published 1D stellar models from FRANEC, MESA, KEPLER, and GENEC identified 41 models with a C–O shell merger. Those models span \(9.2 \leq M_{\rm ini}/M_\odot \leq 25\), \(0 \leq Z_{\rm ini} \leq 0.05\), and \(0 \leq v_{\rm ini} \leq 800\,\mathrm{km\,s^{-1}}\). The key result is that merger occurrence can be approximated from the outcome of central He burning using two quantities: the CO core mass \(M_{\rm CO}\), defined as the helium-exhausted core mass, and the central \(^{12}\mathrm C\) mass fraction \(X_{\rm C12}\) at the end of He-core burning [2504.18867].

The empirical diagnostic is
\[
X_{\rm C12}<0.277,\qquad M_{\rm CO}<4.90\,M_\odot.
\]
Ninety percent of merger models satisfy these bounds, with average values
\[
\langle M_{\rm CO}\rangle \approx 4.02\,M_\odot,\qquad \langle X_{\rm C12}\rangle \approx 0.176.
\]
These conditions are necessary but not sufficient. Within that region of the \((M_{\rm CO},X_{\rm C12})\) plane, 79 models are present, of which 41 merge and 38 do not, implying a merger likelihood of about 50% once a model enters the favored region [2504.18867].

The physical interpretation is twofold. A smaller CO core places later C, Ne, O, and Si burning shells closer together in mass coordinate, which makes interaction more likely. A lower central carbon abundance reduces subsequent Ne production and weakens the entropy jumps at the C–Ne–O interfaces, thereby making upward penetration by the O shell easier. Rotation and metallicity modify \(M_{\rm CO}\) and \(X_{\rm C12}\) indirectly, but rotation is not a direct trigger, and mergers occur at both solar and very low metallicity once the two structural diagnostics fall in the critical range [2504.18867].

## 3. Hydrodynamics, convection, and modeling

In 1D stellar evolution, O–C shell mergers are identified by following convective boundaries, burning rates, and composition profiles through the final stages of evolution. A merger is registered when the O-burning convective shell overtakes the C/Ne shell and a single mixed convective zone forms. Detailed nucleosynthesis is then usually computed with post-processing networks that include heavy-ion fusion, proton captures, \(\alpha\) captures, and photodisintegrations, while the final explosion is often represented by parameterized or semi-analytic prescriptions rather than full neutrino-radiation hydrodynamics [2504.18867].

Three-dimensional simulations show that the multi-D dynamics differ substantially from mixing-length theory. Shell merging is driven by turbulent entrainment and erosion of the stable layer between shells, convective velocities are much faster than MLT predicts, and merged shells develop strongly asymmetric, often dipolar, composition distributions. In a detailed 3D study of late shell merging in a \(20\,M_\odot\) model, entrainment and erosion of stable regions were identified as the main drivers of the merger, and multiple burning phases were observed within the same merged shell [2407.15544].

For O-shell convection with C ingestion, 3D simulations find robust scaling laws:
\[
\dot{M}_{\rm entr}\propto L_{\rm tot},\qquad v_{\rm conv}\propto L_{\rm tot}^{1/3},
\]
and show that carbon burning can provide \(14\%\)–\(33\%\) of the total luminosity, depending on the adopted reaction set. Most runs settle into a quasi-stationary state dominated by a few large-scale convective cells, but an experimental high-feedback case develops global oscillations with Mach numbers in excess of \(0.2\), representing a violent regime akin to a full O–C shell merger [1808.04014].

This multi-D behavior motivated later “3D-inspired” post-processing studies that replace MLT mixing with diffusion profiles featuring a downturn near the boundary, boosted convective velocities, quenched regions, and explicit C-shell ingestion rates. Those studies demonstrate that merger nucleosynthesis depends sensitively on the assumed macro-physics of entrainment and boundary mixing, not only on the underlying stellar structure [2512.17705].

## 4. Energetics and nucleosynthetic products

When an O–C shell merger occurs, the O shell expands by ingesting C/Ne-rich material, the former O and C/Ne shells are replaced by a single extended convective region, and composition profiles are flattened by convection. In the affected layers, C-burning products such as \(^{23}\mathrm{Na}\), \(^{24}\mathrm{Mg}\), and \(^{27}\mathrm{Al}\) are depleted, while O-burning products such as \(^{28}\mathrm{Si}\), \(^{31}\mathrm{P}\), and \(^{32}\mathrm{S}\) are enhanced. Odd-\(Z\) isotopes including Cl, K, and Sc are produced efficiently in the same zone [2504.18867].

A major revision of the standard picture is that, under typical merger conditions in 1D models, the dominant energy source is not classical O fusion but proton captures on S, P, and Ar. This set of reactions is termed the SPAr process and consists chiefly of
\[
^{32}\mathrm S(p,\gamma)^{33}\mathrm{Cl},\quad
^{34}\mathrm S(p,\gamma)^{35}\mathrm{Cl},\quad
^{31}\mathrm P(p,\gamma)^{32}\mathrm S,\quad
^{38}\mathrm{Ar}(p,\gamma)^{39}\mathrm K,
\]
plus the corresponding reverse photodisintegration channels. In a representative \(15\,M_\odot\) model, SPAr proton captures produce approximately 400 times more energy than classical O fusion under the same conditions, and about 330 times more than the sum of C, O, and C+O fusion channels [2509.13749].

The merger environment strongly enhances light-particle abundances. At the bottom of the O shell, the proton mass fraction rises from \(1.54\times10^{-10}\) in classical O-shell burning to \(9.39\times10^{-9}\) at the beginning of the merger and \(7.21\times10^{-9}\) at the end; \(\alpha\)-particle and neutron abundances increase by similar factors, and the neutron density reaches \(n_n\sim 2\times10^{16}\,\mathrm{cm^{-3}}\). These light particles drive the odd-\(Z\) reaction flows. Examples include \(^{34}\mathrm S(p,\gamma)^{35}\mathrm{Cl}\), \(^{38}\mathrm{Ar}(p,\gamma)^{39}\mathrm K\), and the Sc-producing sequence \(^{38}\mathrm{Ar}(\alpha,\gamma)^{42}\mathrm{Ca}(\alpha,p)^{45}\mathrm{Sc}\) [2509.13749].

This nucleosynthesis is largely primary, because the seeds are generated within the O shell itself rather than inherited from the initial metallicity. In yield space, merger models occupy the high-\([\mathrm K/\mathrm{Mg}]\), high-\([\mathrm{Sc}/\mathrm{Mg}]\) region: \([\mathrm K/\mathrm{Mg}] \gtrsim 0\) and \([\mathrm{Sc}/\mathrm{Mg}] \gtrsim 0\), whereas non-merger models, especially at low metallicity, remain sub-solar in both ratios. This is the basis for the claim that C–O shell mergers resolve the longstanding underproduction of K and Sc in standard massive-star yields [2504.18867].

O–C shell mergers are also pre-explosive \(\gamma\)-process sites. In a \(15\,M_\odot\), \(Z=0.02\) post-processing study, the average spread in p-nucleus overproduction across MLT and 3D-inspired mixing prescriptions is \(0.96\) dex, ingestion-rate uncertainties produce spreads of \(1.22\)–\(1.84\) dex depending on the mixing profile, and nuclear-rate variations for unstable neutron-deficient isotopes from Se to Po contribute \(0.56\)–\(0.79\) dex. This places macro-physics and nuclear-physics uncertainties on a comparable footing for quantitative p-nucleus predictions [2507.16965].

## 5. Structural consequences and the approach to explosion

Beyond nucleosynthesis, O–C shell mergers alter the pre-supernova structure itself. They weaken entropy and density discontinuities between shells, can create a density drop at the Si/O interface, and introduce non-radial asymmetries that are expected to be larger in reality than in 1D models. Because the merged region often sits at relatively large radius, the supernova shock modifies but does not erase the hydrostatic signature; in simplified shock arguments, the peak shock temperature scales as \(T_{\rm peak}\propto R_*^{-3/4}\), so only moderate explosive reprocessing occurs in much of the merged shell [2504.18867].

These structural changes are directly relevant to explodability. Closely related studies of strong C-shell burning adjacent to the O shell show that extra heating in the C shell can prevent contraction of outer layers, decrease the compactness parameter at \(M_r=2.5\,M_\odot\), and steepen the density and pressure gradients at the O-burning shell. In that framework, the pressure gradient is described by \(V/U\) and the density drop by \(1/U\), where
\[
U=\frac{4\pi r^3\rho}{M_r},\qquad
V=\frac{G M_r \rho}{rP}.
\]
The proposed critical values are \(\log(V/U_{\rm max})=1.205^{+0.070}_{-0.055}\) and \(\log U_{\rm min}=-0.745^{+0.053}_{-0.032}\), above and below which, respectively, models are more likely to explode and form neutron stars rather than black holes [2502.11012].

Three-dimensional shell-merger simulations support the same qualitative picture. In a \(4\pi\) simulation of a violent oxygen–neon shell merger just before collapse, the merged shell attained convective Mach numbers of \(\approx 0.1\), erased the interface between the Ne layer and the Si-enriched O layer, and developed large-scale \(\ell=1\) and \(\ell=2\) asymmetries. Although that case is an O–Ne merger with partial C-shell involvement rather than a pure O–C merger, it demonstrates that shell mergers can seed the large-scale perturbations known to favor shock revival in core-collapse simulations [1905.04378].

## 6. Observational and chemical-evolution evidence

Stellar archaeology supplies one line of evidence. A comparison to JINAbase abundances selected 121 metal-poor stars with K and Sc measurements in the range \(-4.12 \leq [\mathrm{Fe}/\mathrm{H}] \leq -1.51\), including 37 extremely metal-poor stars with \([\mathrm{Fe}/\mathrm{H}] \leq -3\). The observations occupy a junction region between three model populations: non-merger CCSNe with low K and Sc, C–O merger CCSNe with high K and Sc, and neutrino-enhanced CCSNe with elevated Sc but not K. No single progenitor type reproduces the full observed distribution, but a significant contribution from C–O shell mergers is required to explain near-solar \([\mathrm K/\mathrm{Mg}]\) and moderately sub-solar \([\mathrm{Sc}/\mathrm{Mg}]\) at very low metallicity [2504.18867].

Galactic chemical-evolution modeling reaches a similar conclusion. Earlier work showed that, if merger-like O–C shell interaction yields are taken as representative and such events occur in more than 50% of all stars, one-zone chemical-evolution models reproduce the observed Galactic trends of the odd-\(Z\) elements [1704.05985]. A later analysis linked the improvement at low metallicity specifically to rotating massive-star yields that frequently undergo C–O shell mergers, reinforcing the idea that the mergers are chemically important rather than exceptional curiosities [2504.18867].

Young supernova remnants now provide a second, more direct observational route. A 2026 comparison of eight CCSN model sets to Cassiopeia A identified large \(\mathrm{Ar}/\mathrm{Ne}\), \(\mathrm{Si}/\mathrm{Ne}\), and \(\mathrm{S}/\mathrm{Ne}\) as robust diagnostics of a C–O shell merger. Merger models produce \(\mathrm{Ar}/\mathrm{Ne}\gtrsim 0.1\), \(\mathrm{Si}/\mathrm{Ne}\gtrsim 1\), \(\mathrm{S}/\mathrm{Ne}\gtrsim 1\), and typically \(\mathrm{Ca}/\mathrm{Ne}\gtrsim 0.2\), whereas non-merger models remain at lower ratios. For Cas A, only models with a C–O shell merger match the combined X-ray and infrared abundance constraints [2603.24758].

The same remnant constrains merger-produced \(^{44}\mathrm{Ti}\). Merger-compatible models predict \(5\times10^{-6}\) to \(5\times10^{-5}\,M_\odot\) of \(^{44}\mathrm{Ti}\) synthesized in the C–O shell and eventually ejected, corresponding to roughly \(5\)–30% of the total Cas A \(^{44}\mathrm{Ti}\). Because this component originates above the Si/Si–O interface, it should lie outside the reverse shock and remain below current NuSTAR and COSI detection limits, though it may be accessible to future instruments such as ASCENT. By contrast, for SN 1987A the observed redshift of the \(^{44}\mathrm{Ti}\) line argues against a dominant C–O-merger origin for the measured \(^{44}\mathrm{Ti}\) [2603.24758].

## 7. Uncertainties and future directions

The frequency of O–C shell mergers remains uncertain. The sample of 41 mergers among 209 models is not weighted by realistic distributions in mass, metallicity, and rotation, so it is not itself a cosmic merger fraction. At the same time, the presence of mergers across FRANEC, MESA, KEPLER, and GENEC suggests that they are a generic possible feature of \(12\)–\(25\,M_\odot\) massive-star evolution rather than a code-specific artifact [2504.18867].

The dominant theoretical uncertainties are convective-boundary physics, dimensionality, and nuclear rates. One-dimensional models treat convection with MLT plus diffusive mixing, whereas 3D calculations show more efficient entrainment, higher convective velocities, and possible quenching or global oscillatory states. This suggests that mergers may be more common in nature than in current 1D grids, but also that the detailed extent of mixing and burning in any given event may differ substantially from the 1D prediction [2504.18867].

Nuclear uncertainties remain central. The \(^{12}\mathrm C(\alpha,\gamma)^{16}\mathrm O\) rate is relatively well constrained, with an uncertainty of about 20%, but shifts in that rate change the post-He-burning \(X_{\rm C12}\) and therefore the location of models in merger-favored parameter space. The heavy-ion fusion rates \(^{12}\mathrm C+^{12}\mathrm C\) and \(^{12}\mathrm C+^{16}\mathrm O\) are more uncertain, and they directly affect shell energetics and the onset of convective-reactive behavior. Future measurements at LUNA, JUNA, and LNS-INFN are expected to tighten those constraints [2504.18867].

Quantitatively, uncertainties in 3D macro-physics are already as large as the nuclear uncertainties in some channels. For p nuclei, the spread from different mixing prescriptions rivals the spread from varying photodisintegration rates; for light odd-\(Z\) products in 24 mixing cases, the pre-explosive yields of \([\mathrm P/\mathrm{Fe}]\), \([\mathrm{Cl}/\mathrm{Fe}]\), \([\mathrm K/\mathrm{Fe}]\), and \([\mathrm{Sc}/\mathrm{Fe}]\) vary over the intervals \([-0.33,0.23]\), \([-0.84,0.64]\), \([-0.78,1.48]\), and \([-0.36,1.29]\) dex, respectively, and the inferred \(^{40}\mathrm K\) yield varies by more than three orders of magnitude [2507.16965, 2509.19240].

The immediate research program implied by these results is well defined: larger systematic grids that vary single stellar-physics inputs, more multi-D hydrodynamic simulations of late burning shells with realistic nuclear energetics, improved low-metallicity K and Sc observations, and chemical-evolution models that disentangle the roles of rotation, neutrino physics, and O–C shell mergers. A plausible implication is that the field is moving from asking whether O–C shell mergers can occur to asking how often they occur, under what structural conditions, and how strongly their multi-D dynamics reshape the final hours of massive-star evolution.

Source: https://www.emergentmind.com/topics/o-c-shell-mergers