Shell Merger in Galaxies and Stars

• Shell merger is the process by which distinct, stratified shells—whether from galactic interactions or stellar burning zones—merge through dynamical, hydrodynamical, or nuclear mechanisms.
• Analytical models and simulations show that in galaxies, nearly radial mergers lead to phase-wrapping, suppressing angular momentum and producing sharp shell caustics that constrain gravitational potentials.
• In stellar evolution, shell mergers mix burning layers to alter core structure and nucleosynthetic yields, impacting supernova explosions and chemical evolution with enhanced odd-Z and radioactive isotopes.

A shell merger is a process in both galaxy dynamics and stellar evolution wherein distinct, energetically and/or compositionally stratified structures—the so-called "shells"—interact, overlap, or combine via dynamical, hydrodynamical, or nuclear processes. While the terminology is rooted in different physical domains (galactic versus stellar), the unifying feature is the merger or interpenetration of sharp boundaries separated by significant phase-space, angular momentum, or compositional differences. Shell mergers serve as powerful diagnostic probes of assembly history in galaxies and the late evolutionary pathways in massive stars, profoundly affecting observable signatures, internal structure, explosive outcomes, and nucleosynthetic yields.

1. Shell Merger in Galaxy Dynamics: Morphology and Kinematic Signature

In extragalactic astronomy, "shell merger" describes the formation of large-scale stellar shell systems in the outskirts of galaxies, predominantly early-type galaxies (ETGs), as a direct consequence of radial or near-radial mergers. Shells appear as concentric, sharp-edged arcs or ripples, usually symmetric about the galaxy center, and are morphologically distinct from tidal streams or tails. Yoon et al. (2024) established that shell detection is far more prevalent in slow-rotating ETGs, with nearly half of shell-bearing ETGs classified as slow rotators ($f_\mathrm{slow}\approx0.46\pm0.06$), in contrast to systems with tails ($f_\mathrm{slow}\approx0.22$), streams ($f_\mathrm{slow}\approx0.20$), or no features ($f_\mathrm{slow}\approx0.14$) (Yoon et al., 2024).

The physical mechanism is phase-wrapping: during a nearly head-on ("radial") merger, a disrupting satellite's stars are deposited onto low-angular-momentum orbits within the host. As these stars reach orbital apocenter, their velocity slows and they accumulate into sharp density caustics—observable as shells. The detection of shells in more massive ETGs is attributed to dynamical friction, which both drags massive satellites towards radial orbits and erodes their angular momentum, establishing a strong connection between shell incidence and merger orbit geometry.

Specific stellar angular momentum, quantified by $\lambda_{R_e}$ (the flux-weighted ratio of projected rotation to total support within the effective radius), is sharply suppressed in shell hosts: median $\lambda_{R_e}\approx0.13$ for shells vs. $\approx0.22$ for tails/streams and $\approx0.31$ for relaxed ETGs. This deficit is robust in all mass and velocity dispersion bins, indicating that radial orbit geometry, rather than mass ratio alone, dictates post-merger kinematics (Yoon et al., 2024).

2. Phase-Wrapping: Theoretical Framework and Timescales

The phase-wrapping model underpins shell formation dynamics. Stripped stars released near the host's center with a distribution of binding energies oscillate radially, with each energy corresponding to a unique radial period $T(E)$. For integer $n$, the locus where stars have completed $f_\mathrm{slow}\approx0.22$0 half-oscillations forms a caustic—i.e., a shell—at radius $f_\mathrm{slow}\approx0.22$1 determined by $f_\mathrm{slow}\approx0.22$2 in a spherical potential (Dong-Páez et al., 2021, Paudel et al., 2016, Pop et al., 2017, Ebrova, 2013).

Shell edges propagate outward at a phase velocity $f_\mathrm{slow}\approx0.22$3, scaling as $f_\mathrm{slow}\approx0.22$4, where $f_\mathrm{slow}\approx0.22$5 is the time since stripping (in the small radius limit). Multiple stripping episodes (typically at each pericentric passage) produce distinct shell systems interleaved in radius, as shown via action–angle techniques in six-dimensional phase space (Dong-Páez et al., 2021). Typical visibility timescales for sharp shells are $f_\mathrm{slow}\approx0.22$6 Gyr after coalescence. Phase mixing erases the shell signature beyond $f_\mathrm{slow}\approx0.22$7–$f_\mathrm{slow}\approx0.22$8 radial periods post-stripping, restricting their utility as markers of very recent mergers (Paudel et al., 2016, Cooper et al., 2011).

3. Quantitative Shell Kinematics and Gravitational Potential Constraints

Spatially resolved kinematics of shells provide independent constraints on the gravitational potential, especially in the outer regions of massive galaxies. Analytical work and test-particle simulations demonstrate that the line-of-sight velocity distribution (LOSVD) of shell stars near a propagating shell adopts a quadruple-peaked structure—the splitting of the classic double-peaked profile due to shell expansion velocity $f_\mathrm{slow}\approx0.22$9 (Ebrova et al., 2012, Ebrova, 2013). For projected radius $f_\mathrm{slow}\approx0.20$0 inside a shell at $f_\mathrm{slow}\approx0.20$1, the four LOSVD peaks occur at

$f_\mathrm{slow}\approx0.20$2

where $f_\mathrm{slow}\approx0.20$3 is the circular velocity at $f_\mathrm{slow}\approx0.20$4. Measuring the $f_\mathrm{slow}\approx0.20$5-dependence of the LOSVD peak separations allows simultaneous inference of $f_\mathrm{slow}\approx0.20$6 and $f_\mathrm{slow}\approx0.20$7, and hence a direct probe of the dark matter distribution out to $f_\mathrm{slow}\approx0.20$8 kpc in early-type galaxies (Ebrova et al., 2012, Ebrova, 2013).

4. Shell Merger Statistics, Cosmological Context, and Environmental Dependence

Cosmological simulations (e.g., Illustris, Aquarius) robustly predict that shell systems arise primarily from relatively major ($f_\mathrm{slow}\approx0.20$9) mergers on low-angular-momentum orbits. At $f_\mathrm{slow}\approx0.14$0, 18% $f_\mathrm{slow}\approx0.14$1 3% of the 220 most massive Illustris galaxies ($f_\mathrm{slow}\approx0.14$2) exhibit shells, rising to $f_\mathrm{slow}\approx0.14$3 for $f_\mathrm{slow}\approx0.14$4; this incidence declines to $f_\mathrm{slow}\approx0.14$5 by $f_\mathrm{slow}\approx0.14$6 due to shorter dynamical times and rapid phase-mixing at earlier epochs (Pop et al., 2017). Shell-producing satellites predominantly have mass ratios $f_\mathrm{slow}\approx0.14$7 (i.e., major/minor mergers) and infall radialities $f_\mathrm{slow}\approx0.14$8. The temporal window for visible shells coincides with mergers accreted $f_\mathrm{slow}\approx0.14$94–8 Gyr prior to observation, whose stars are typically stripped $\lambda_{R_e}$02 Gyr ago.

Environmental complexity, such as multiple shell-forming progenitors, satellite-of-satellite infall, or rapid successive mergers, can produce overlapping, misaligned, or prematurely erased shell systems. Deep imaging surveys confirm the presence of shell features with $\lambda_{R_e}$1–29 mag arcsec$\lambda_{R_e}$2 in both giant and dwarf early-type galaxies, with layering, radial extent, and alignment matching phase-wrapping plus dynamical friction predictions (Paudel et al., 2016, Cooper et al., 2011).

5. Shell Merger Phenomena in Stellar Evolution: Internal Mixing and Nucleosynthesis

In stellar contexts, a "shell merger" refers to the hydrodynamic or convective merging of distinct nuclear-burning shells during the late phases of massive star evolution (M ≳ 12 M⊙). Major types include C–O, O–Ne, Si–C, and Ne–O shell mergers. Such events are typically mediated by turbulent entrainment, convective boundary erosion, and associated mixing driven by rapid core contraction and burning luminosity escalation.

3D hydrodynamic simulations reveal that the merger proceeds via the downward migration of the convective boundary ("entrainment"), mixing fresh fuel (e.g., C, Ne) into hotter burning zones, inducing convective-reactive burning and complex, multi-zone nuclear energy generation within an enlarged single convective region (Rizzuti et al., 2024, Yadav et al., 2019). 3D simulations find much faster convective velocities (by factors of 5–10 over MLT), strong large-scale asymmetries (e.g., dipolar abundance distributions), and only partial homogenization prior to collapse. The net effect is profound alteration of the presupernova core structure (entropy, density, composition) and mixing of nuclear species, impacting subsequent shock propagation and explosion dynamics (Roberti et al., 26 Apr 2025, Côté et al., 2019).

6. Nucleosynthetic and Chemical Evolution Consequences

Shell mergers in massive stars have decisive effects on the nucleosynthesis and expulsion of both stable and radioactive nuclei. In particular:

• C–O Shell Mergers: Ingestion of carbon into the O-burning zone at $\lambda_{R_e}$3 K boosts production of odd-Z isotopes (Sc, K) and p-nuclei (e.g., $\lambda_{R_e}$4Mo, $\lambda_{R_e}$5Ba), potentially solving longstanding Galactic chemical evolution deficiencies in these elements. Survey of 209 models finds that C–O merger incidence can be robustly predicted from post–He-burning CO core mass and $\lambda_{R_e}$6: 90% of mergers have $\lambda_{R_e}$7 and $\lambda_{R_e}$8 (Roberti et al., 26 Apr 2025).
• O–C Shell Mergers: Such events can lead to "pre-explosive" production of $\lambda_{R_e}$9Ti (up to 5 dex above standard predictions), as well as enhanced $\lambda_{R_e}\approx0.13$0Na, $\lambda_{R_e}\approx0.13$1Cr, $\lambda_{R_e}\approx0.13$2Fe, altering gamma-ray line diagnostics and the chemical signatures of pre-solar grains. The total yield can exceed that from the explosion itself if the entrainment rate is sufficiently high and convective mixing is strongly boosted (factors ≥ 10 above MLT) (Issa et al., 19 Dec 2025).
• Si–C Shell Mergers: These produce an outward mixing of incomplete Si-burning products (notably $\lambda_{R_e}\approx0.13$3V, $\lambda_{R_e}\approx0.13$4Cr), which can lead to anomalously high [Cr/Fe] abundance ratios in nucleosynthetic yield tables. However, observational constraints from Galactic disk stars suggest such events must be rare or less efficient than 1D models predict (Côté et al., 2019).
• Ne/O Shell Intrusions: Observations of supernova remnants with high Mg/Ne and Si/Mg ratios (e.g., G359.0−0.9: $\lambda_{R_e}\approx0.13$5, $\lambda_{R_e}\approx0.13$6) imply shell-merger induced mixing, allowing SNRs to be used as post-explosion diagnostics of shell merger occurrence and progenitor mass ($\lambda_{R_e}\approx0.13$7 for Ne-intrusion, $\lambda_{R_e}\approx0.13$8 for O-shell merger) (Matsunaga et al., 2024).

7. Broader Astrophysical Implications and Observational Signatures

Shell mergers, whether galactic or stellar, imprint testable observational signatures:

• Galactic shells: Serve as fossil records of recent, typically radial, major or minor mergers. Their morphology, kinematic profile, and persistence time constrain both the orbital parameters and mass assembly histories of ETGs. Shell kinematics enable mapping of host galaxy gravitational potentials well beyond optical radii and provide stringent tests for CDM-based merger rates and mass assembly (Yoon et al., 2024, Pop et al., 2017, Ebrova et al., 2012).
• Stellar shell mergers: Drastically reshape the structure of massive stars prior to core collapse, create nucleosynthetic anomalies (particularly in odd-Z and radioactive species) that propagate into supernova light curves, SNR abundance profiles, and gamma-ray emission. The extent, timing, and efficiency of such mergers, as well as their 3D versus 1D modeling, are essential uncertainties for predicting massive star fates and interpreting chemical evolution signatures (Roberti et al., 26 Apr 2025, Issa et al., 19 Dec 2025, Rizzuti et al., 2024).

A secondary form of "shell merger" arises in the context of mass transfer or accretion between evolved stars (e.g., the merger of a low-mass companion with a red supergiant envelope), which spins up the envelope and can create rapidly rotating supergiants like Betelgeuse without significantly impacting the inner core structure or magnetic field generation (Sullivan et al., 2020). In specific SNe Ia-CSM, shell ejection years to decades before explosion, for example from CO WD mergers, can lead to circumstellar interaction light-curve bumps, as seen in SN 2020aeuh (Chugai, 19 Oct 2025).

In sum, shell mergers are critical, multidomain phenomena that provide stringent constraints and diagnostic leverage across extragalactic, stellar evolutionary, and nucleosynthetic contexts. Their robust identification and precise modeling—especially in 3D, turbulence-resolved regimes—are imperative for unlocking their role in the assembly and chemical enrichment of galaxies and the evolution of massive stars.