Deciphering the Remnants of Core-Collapse Supernovae: Reconstructing Progenitor Star Properties and Explosion Mechanisms

Technical Summary: Deciphering the Remnants of Core-Collapse Supernovae

Problem Statement
Recent high-resolution observations of the Cassiopeia A (Cas A) supernova remnant (SNR) by the James Webb Space Telescope (JWST) have revealed unprecedented structural details that challenge current theoretical understanding. Two primary features require explanation: (1) an intricate, web-like network of unshocked, oxygen-rich (O-rich) filaments in the remnant's interior, and (2) the "Green Monster" (GM), an enigmatic region on the blueshifted (near) side characterized by a pockmarked appearance with nearly circular holes and rings in the shocked circumstellar medium (CSM). Traditional models often treat the supernova (SN) explosion, the evolution of the ejecta, and the interaction with the circumstellar environment as separate phases. This paper addresses the need for a unified physical narrative that links the earliest instabilities of the core-collapse explosion to the specific morphological features observed centuries later, bridging the gap between the explosion engine and the progenitor's mass-loss history.

Methodology
The authors employ a unified numerical framework utilizing high-resolution, fully three-dimensional (3D) hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations. The approach follows the system continuously from the moment of core collapse to an age of approximately 1,000 years.

• Initial Conditions: The simulations begin with the W15-2-cw-IIb model, a 3D neutrino-driven explosion of a 15 $M_{\odot}$ zero-age main sequence progenitor. The progenitor was artificially stripped to a Type IIb configuration (final radius $\sim$21.5 $R_{\odot}$) to match Cas A, retaining only $\sim$0.3 $M_{\odot}$ of hydrogen. The explosion was calibrated to a kinetic energy of $\sim$1.5 $\times$ 10$^{51}$ erg, ejecting 3.3 $M_{\odot}$ of material and synthesizing $\sim$0.1 $M_{\odot}$ of radioactive $^{56}$Ni.
• Physical Processes: The models incorporate neutrino-driven explosion dynamics, hydrodynamic instabilities (specifically Rayleigh-Taylor and Kelvin-Helmholtz), the "Ni-bubble" effect (expansion of radioactive ejecta compressing surrounding layers), and optically thin radiative cooling. Crucially, the simulations account for non-equilibrium ionization and electron-proton temperature equilibration.
• Resolution and Scope: Using a grid of 2048$^3$ cells and a remapping technique, the simulations resolve structures down to $\sim$0.01 pc, comparable to JWST's spatial resolution. Chemical elements are tracked via passive tracers to reconstruct spatial distributions of O, Si, Fe, and Ti.
• CSM Interaction: To model the Green Monster, the ejecta interact with a thin, dense, asymmetric circumstellar shell located at $\sim$1.5 pc, representing a major mass-loss episode occurring $10^4$–$10^5$ years prior to the explosion. The shell is modeled with enhanced density on the blueshifted side.
• Observational Alignment: Simulations are rotated to match the observed distribution of Fe- and Ti-rich ejecta in Cas A, allowing for direct comparison with multiwavelength data (JWST, Chandra, NuSTAR).

Key Contributions and Results

1. Origin of the Filamentary O-Rich Network:
   The simulations demonstrate that the complex web of unshocked O-rich filaments is not a product of late-time hydrodynamic evolution but a direct "fossil" imprint of the explosion mechanism established within the first second after core collapse.

   • Mechanism: Neutrino-driven plumes break spherical symmetry immediately. As the shock propagates, these plumes drive vigorous mixing. The subsequent formation and expansion of hot, buoyant, Ni-rich bubbles (the Ni-bubble effect) compress and stretch the surrounding O-rich layers.
   • Evolution: Rayleigh-Taylor (RT) instabilities rapidly fragment these compressed layers into a filamentary pattern. By the time the remnant is a few years old, a highly organized, interconnected web of dense O-rich filaments occupies the unshocked interior.
   • Persistence: At the current age ($\sim$350 years), the simulated network matches JWST images in geometry, hierarchy, and thickness ($\sim$0.01 pc). The model predicts this structure is transient; as Cas A approaches $\sim$700 years, the reverse shock will erode the filaments via RT and Kelvin-Helmholtz instabilities, erasing the JWST-scale morphology. This implies that such networks in older Galactic remnants have likely already been destroyed.

2. Origin of the Green Monster (GM):
   The peculiar morphology of the GM is reproduced by the interaction of fast, clumpy SN ejecta with the pre-existing asymmetric circumstellar shell.

   • Mechanism: Dense, radiatively cooled RT fingers of ejecta, accelerated to several thousand km/s, strike the circumstellar shell. These fingers act as ballistic streams, punching through the shell and creating narrow perforations.
   • Morphology: The displaced shell material is swept laterally, accumulating around the penetration sites to form near-circular rims. When multiple fingers strike in close proximity, they create compound or merged rings, reproducing the "Swiss-cheese" appearance.
   • Role of Physics: Radiative cooling is decisive, allowing the fingers to maintain high density contrasts and coherence long enough to create clean holes rather than dispersing. Magnetic fields further confine these knots, enhancing their ability to burrow through the shell.
   • Validation: The simulated holes (typically $\sim$5 arcseconds) match JWST measurements. The chemical composition within the holes is predicted to be enriched in O, C, and Ne, reflecting the intermediate-mass layers of the progenitor.

Significance and Claims
The paper claims that these studies provide a unified physical interpretation linking the explosion engine, ejecta evolution, and progenitor mass-loss history.

• Archaeological Record: The authors posit that young SNRs like Cas A function as "archaeological records." The internal filaments preserve the hydrodynamic fingerprints of the explosion's first seconds (neutrino heating, mixing instabilities), while the GM records the final centuries of the progenitor's life (mass-loss geometry).
• Unified Framework: The work highlights the necessity of modeling stellar evolution, SN dynamics, and remnant development in a single, self-consistent framework. Treating these phases separately fails to explain the specific morphological features now visible.
• Broader Impact: By connecting the physics of the progenitor, the explosion, and the remnant, these models allow for the extraction of physical information previously inaccessible, such as the nature of neutrino-driven turbulence and the details of massive star mass loss immediately before collapse. The authors suggest this approach will become increasingly vital as new observatories resolve young remnants like SN 1987A, transforming them into precise diagnostic tools for understanding how massive stars end their lives.