Mapping plasma properties of Cassiopeia A with XRISM/Resolve: a Bayesian analysis via UltraSPEX

Mapping the physical conditions of the shocked plasma of young supernova remnants (SNR) is crucial for understanding their explosion mechanisms, ejecta structure, and large-scale asymmetries. Using $>350$ ks of XRISM/Resolve high spectral resolution observations of Cassiopeia A (Cas A), the youngest known Galactic core-collapse SNR, we present the first microcalorimeter-based plasma parameter maps of any SNR. We tessellate Cas A into $1’\times1’$ regions and fit the broadband spectra as thermal emission from two pure-metal ejecta components – corresponding to intermediate-mass elements (IMEs) and iron-group elements (IGEs) – plus nonthermal synchrotron radiation. For robust inference, we introduce UltraSPEX, a Bayesian framework that couples the SPEX plasma code with the UltraNest nested-sampling algorithm, yielding full posterior distributions and exploration of parameter degeneracies. Key findings include enhanced Ar/Si and Ca/Si abundance ratios near the base of the Si-rich jets, and a high Ni/Fe mass ratio ($0.08\pm0.015$) in the base of NE jet. IGEs ejecta exhibit systematically higher Doppler velocities and broadenings than IMEs ejecta in most regions, with maximum differences of $\sim800$ km/s and $\sim1200$ km/s, respectively; Ca shows distinct (faster) kinematics from other IMEs in several SE regions. The ionization timescale and electron temperature show a robust anti-correlation, particularly for IGEs. This relation and measured parameter values could be explained by semi-analytical models with significant ejecta clumping (overdensities of $\sim10$ for IGEs and up to $\sim100$ for IMEs) and reduced historical reverse-shock velocities. Nonthermal emission accounts for a substantial fraction, with at least 47% of the 4–6 keV continuum and dominates in the western regions, where the spectrum hardens.

💡 Research Summary

This paper presents the first micro‑calorimeter‑based plasma property maps of any supernova remnant (SNR), using more than 350 ks of XRISM Resolve high‑resolution X‑ray observations of Cassiopeia A (Cas A), the youngest known Galactic core‑collapse SNR. The authors tessellate the remnant into 1′ × 1′ “super‑pixels” (each comprising 2 × 2 Resolve detector pixels) and fit the broadband (1.8–11.9 keV) spectra with a physically motivated model that includes two pure‑metal ejecta components—one representing intermediate‑mass elements (IMEs) and the other iron‑group elements (IGEs)—plus a non‑thermal synchrotron component.

To obtain robust, unbiased parameter estimates, they develop UltraSPEX, a Bayesian framework that couples the SPEX plasma code with the UltraNest nested‑sampling algorithm. UltraNest efficiently explores the high‑dimensional parameter space (electron temperature, ionization age, elemental mass fractions, Doppler shifts and broadenings, synchrotron power‑law index and normalization, etc.) and returns full posterior probability distributions, allowing a rigorous assessment of uncertainties and parameter degeneracies.

Key results are:

1. Abundance patterns – The Ar/Si and Ca/Si ratios are significantly enhanced at the bases of the Si‑rich jets, especially in the northeast (NE) jet. The Ni/Fe mass ratio in the same region is measured as 0.08 ± 0.015, higher than previous estimates, indicating either Ni over‑production in the explosion or selective heating of Ni‑rich material.

2. Kinematic differences – IGE ejecta show systematically larger bulk Doppler velocities (up to ~800 km s⁻¹) and broader line widths (up to ~1200 km s⁻¹) than IME ejecta across most tiles. Calcium, in several southeast (SE) tiles, displays a distinct, faster kinematic signature compared with the other IMEs, suggesting a separate dynamical component or mixing process.

3. Temperature–ionization anti‑correlation – Electron temperature (kTₑ) and ionization timescale (nₑt) are strongly anti‑correlated, particularly for the IGE component. The authors interpret this as evidence for strong clumping: IGE plasma appears to be overdense by a factor of ~10, while IME plasma may be overdense by up to ~100, combined with a reduced historical reverse‑shock velocity. Semi‑analytical models incorporating these clump factors reproduce the observed anti‑correlation.

4. Non‑thermal contribution – The synchrotron component accounts for at least 47 % of the 4–6 keV continuum, dominating the western side where the spectrum hardens (power‑law index ≈ 2.2). This fraction is higher than earlier CCD‑based estimates (25–50 %) and benefits from XRISM’s low instrumental background and sensitivity up to ~12 keV.

Methodologically, the paper demonstrates that a plane‑parallel shock (pshock) model, implemented as a linear combination of 200 NEI layers with a constant electron temperature and a logarithmically spaced ionization age distribution, better captures the complex ionization structure of Cas A than a single‑zone NEI model. The Bayesian approach avoids the local‑minimum traps common in χ² minimization and provides a transparent view of parameter covariances, which is crucial when fitting multi‑component spectra with overlapping line complexes.

The authors discuss the broader implications: the ability to map detailed plasma properties with high spectral resolution opens a new window on the three‑dimensional structure and explosion asymmetries of young SNRs. The observed Ni enrichment and velocity offsets support models where the explosion was highly anisotropic, possibly involving jet‑like outflows that carried distinct nucleosynthetic yields. The clumping factors inferred from the temperature–ionization relation align with high‑resolution optical and infrared studies that have identified dense knots and filaments within Cas A.

Finally, the paper argues that UltraSPEX, by delivering full posterior distributions for complex, multi‑component plasma models, will become a standard tool for analyzing forthcoming XRISM and Athena data sets, not only for SNRs but also for galaxy clusters, active galactic nuclei, and other extended X‑ray sources where spatially resolved spectroscopy is essential. Future work will aim to integrate these plasma maps with three‑dimensional hydrodynamic simulations to further constrain the progenitor’s explosion mechanism and the subsequent evolution of the remnant.