REVIEW 2 major objections 9 minor 53 references
Reviewed by Pith at T0; open to challenge.
T0 review · glm-5.2
Supernova 1987A's X-ray ring grows energy-dependent after 2012
2026-07-08 10:15 UTC pith:AEUXEO5V
load-bearing objection Energy-dependent torus broadening in SN 1987A is a real new result; the alignment concern is worth raising but probably does not sink the signal the 2 major comments →
Chandra X-Ray Imaging and Spatially Resolved Spectroscopy of SN 1987A: Energy-Dependent Morphology of the Equatorial Ring
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central finding is that the X-ray ring of SN 1987A has developed an energy-dependent radial structure: after the early 2010s, soft X-rays (0.5–1.5 keV) trace a broader distribution that extends inward to smaller radii than hard X-rays (1.5–7.0 keV). This inward extension of soft emission, with the inner boundary reaching roughly 0.5 arcseconds by the 2020s, signals a growing contribution from interior or high-latitude plasma — likely reverse-shock-heated ejecta — marking the remnant's transition away from a purely equatorial-ring-dominated phase.
What carries the argument
The analysis uses Chandra ACIS and HETG observations binned into 4-year intervals, subpixel event repositioning (EDSER) at quarter-pixel scale, Richardson–Lucy deconvolution with MARX-simulated point-spread functions, and an elliptical Gaussian torus model fitted to the deconvolved images to extract ring radius and width as a function of energy band and epoch. Spatially resolved spectroscopy divides the remnant into four quadrants (East, West, North, South) and fits an absorbed power-law plus Gaussian model in the 5–8 keV band to characterize Fe K line emission.
Load-bearing premise
The internal alignment strategy uses the remnant itself as a reference for registering multi-epoch images rather than external point sources, which assumes that the centroid of the X-ray emission does not shift between epochs in a way that would bias the merged image. If the emission centroid migrates — for example, due to the east-west brightness asymmetry reversal reported around 9500 days — the alignment could introduce systematic radial structure differences between bands
What would settle it
If the energy-dependent broadening of the soft-band torus width disappears or reverses when a different center-determination or alignment method is used, the result would be called into question.
If this is right
- If the inward soft-band extension traces reverse-shock-heated ejecta, continued monitoring should show the soft X-ray morphology becoming increasingly distinct from the hard-band ring as more ejecta is shock-heated.
- The eastern enhancement of Fe K emission, if confirmed to predate 2018, may point to an intrinsic asymmetry in the explosion or in the circumstellar environment that predates the current epoch.
- The transition from equatorial-ring-dominated to interior-emission-dominated X-ray morphology provides an observational benchmark for testing 3D hydrodynamic and magnetohydrodynamic simulations of supernova remnant evolution at 30–40 years post-explosion.
- Future X-ray observatories with higher angular resolution could resolve whether the inward soft emission is structured (e.g., bubble-like, as suggested by JWST) or smoothly distributed, distinguishing between reverse-shock and high-latitude origin scenarios.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a systematic Chandra imaging and spatially resolved spectroscopic study of SN 1987A spanning 1999–2025, combining ACIS and HETG zeroth-order data in 4-year bins. The authors employ EDSER subpixel reconstruction, MARX-simulated PSFs, and Richardson–Lucy deconvolution to characterize the projected torus morphology in soft (0.5–1.5 keV) and hard (1.5–7.0 keV) bands. The central new result is that the soft-band torus width becomes systematically broader than the hard-band width after the early 2010s, and when radius and width are considered together, the soft-band emission extends inward to r_ell ~ 0.5 arcsec in the 2020s. The authors interpret this as evidence for an increasing contribution from reverse-shock-heated and/or high-latitude material. A complementary spectral analysis tracks the Fe K line flux evolution and its east–west asymmetry. The analysis is methodologically careful, with Monte Carlo error estimation (300 realizations), cross-checks against individual observations (Appendix C), and a model-independent radial-profile analysis (Appendix D).
Significance. The paper provides a valuable, uniformly analyzed long-term Chandra dataset that quantifies the emergence of energy-dependent radial structure in SN 1987A as it transitions beyond the ER-dominated phase. The finding that the soft-band emission extends inward in the 2020s is physically well-motivated and consistent with HD/MHD predictions (Orlando et al. 2015–2025) and JWST reverse-shock observations (Larsson et al. 2023). The Fe K spatial analysis, while not a new discovery, provides an independent Chandra-based confirmation of the east–west asymmetry and extends its possible detection to earlier epochs. The multiwavelength comparison in Section 3.3 provides useful context. The analysis pipeline (RL deconvolution with MARX PSFs, Monte Carlo uncertainties, multiple cross-checks) is thorough and reproducible in principle.
major comments (2)
- Section 2.3–2.4: The center position for each observation is determined exclusively from the 1.5–7.0 keV band and then applied unchanged to the 0.5–1.5 keV images. If the soft-band emission centroid migrates relative to the hard-band centroid over time — which is physically plausible given that the two bands trace increasingly different components and the east–west asymmetry reverses around 9500 days (Section 4) — the soft-band torus fit would be systematically broadened because the model center is offset from the true soft-band centroid. The width difference in Figure 3b is ~0.05–0.10 arcsec in the latest epochs; a soft–hard centroid offset of comparable size (within the scatter seen in Table 3 center coordinates, which ranges over ~0.1–0.3 arcsec across observations) could account for a substantial fraction of this signal. Critically, all three analysis pathways — merged torus fitting,
- individual-observation fitting (Appendix C), and radial-profile FWHM (Appendix D) — use the same hard-band-derived center positions, so a shared systematic bias would not be caught by any of these cross-checks. The paper should either (a) demonstrate quantitatively that a soft–hard centroid offset of the plausible magnitude cannot reproduce the observed width difference (e.g., by injecting simulated offsets and re-fitting), or (b) repeat the soft-band torus fitting with centers determined independently from the soft-band images and show that the width evolution persists. Without one of these tests, the central claim of energy-dependent width evolution remains vulnerable to this systematic.
minor comments (9)
- Section 2.1: The pileup fraction is stated as '≲10%' for most observations, but several bright epochs between 2002 and 2004 reached higher values. Please quantify the peak pileup fraction in those epochs and briefly discuss whether it could bias the 2000–2003 or 2004–2007 torus measurements.
- Table 1: The 2008–2011 NONE observations (~20 ks) are excluded because contemporaneous HETG data provide deeper coverage. Please state explicitly whether any HETG data from this interval exist or whether there is a temporal gap, as this affects the continuity of the time series.
- Figure 3: Error bars are not visible or not described in the caption. Please confirm that Monte Carlo uncertainties (Appendix B) are plotted and state this in the caption.
- Section 3.2, Table 2: Several early-epoch fits are marked 'N/C' (not constrained). Please define this abbreviation in the table notes (it currently appears only in the footnote marker).
- Section 3.2: The Gaussian width is fixed at sigma = 0.04 keV based on XRISM results. Since XRISM has much higher spectral resolution than ACIS, please briefly justify why this narrow-line assumption is appropriate for the ACIS data and whether the Fe K flux is sensitive to this choice.
- Figure 4 caption: 'Sourth' should be 'South' in panel (e) label.
- Appendix D: The FWHM is described as measured from RL-deconvolved profiles, but the method for determining the half-maximum level (especially in the presence of the central depression or asymmetric wings) is not specified. Please add a brief description.
- Section 3.3: The multiwavelength torus radii are compared, but the X-ray values come from a single ObsID (25514) while the radio/optical/NIR values come from different instruments and modeling approaches. A sentence acknowledging that the absolute values are not directly comparable across wavelengths would strengthen this section.
- Table 3: The center coordinate scatter for HETG observations in 2018 (e.g., ObsIDs 21037, 21038, 21042) shows RA offsets of ~0.3–0.5 arcsec from the mean. Please comment on whether these outliers affect the merged image quality for the 2015–2018 bin.
Simulated Author's Rebuttal
We thank the referee for a careful and constructive report. The referee raises one major concern: the center position used for torus fitting is determined exclusively from the hard band (1.5–7.0 keV) and applied unchanged to the soft band (0.5–1.5 keV), which could systematically broaden the soft-band torus if a soft–hard centroid offset exists. We agree this is a valid and important concern. We will address it by repeating the soft-band torus fitting with centers determined independently from the soft-band images and by performing a simulation-based test injecting known centroid offsets. We expect the width evolution to persist, but if it does not, we will revise the central claim accordingly.
read point-by-point responses
-
Referee: Section 2.3–2.4: The center position for each observation is determined exclusively from the 1.5–7.0 keV band and then applied unchanged to the 0.5–1.5 keV images. If the soft-band emission centroid migrates relative to the hard-band centroid over time, the soft-band torus fit would be systematically broadened. All three analysis pathways use the same hard-band-derived centers, so a shared systematic bias would not be caught by cross-checks. The paper should either (a) demonstrate quantitatively that a soft–hard centroid offset cannot reproduce the observed width difference, or (b) repeat the soft-band torus fitting with independently determined soft-band centers and show the width evolution persists.
Authors: We agree that this is a legitimate and important concern. The referee is correct that all three analysis pathways (merged torus fitting, individual-observation fitting, and radial-profile FWHM) share the same hard-band-derived center positions, so they cannot independently catch a systematic bias from a soft–hard centroid offset. We will address this in the revised manuscript by implementing both tests the referee suggests. First, we will repeat the soft-band torus fitting using centers determined independently from the 0.5–1.5 keV images for each epoch. We note that at early epochs (before ~2012), the soft and hard morphologies are similar, so the centroid offset should be minimal; the critical test is whether the width difference persists at late epochs when the soft-band center is allowed to differ. Second, we will perform a simulation-based test: we will inject known centroid offsets (ranging from 0.05 to 0.3 arcsec, covering the scatter seen in Table 3) into synthetic soft-band images, re-fit the torus model with the hard-band center, and quantify the resulting artificial broadening. This will allow us to determine whether a plausible offset can account for the observed ~0.05–0.10 arcsec width difference. We expect the width evolution to persist with independently determined soft-band centers, because the radial-profile analysis in Appendix D already shows that the soft-band profiles have more extended inner wings (not just a symmetric broadening that a pure centroid offset would produce), but we will present the quantitative results transparently and revise the central claim if the tests indicate the signal is not robust. We will add a new subsection or appendix describing these tests and their outcomes. revision: yes
Circularity Check
No circularity: observational paper with self-contained measurements and no derivation chain that reduces to its inputs.
full rationale
This is an observational astronomy paper measuring Chandra X-ray morphology and spectroscopy of SN 1987A. There is no first-principles derivation chain that could be circular. The torus model parameters (radius r_ell, width sigma_r) are independently fitted to observed images in each energy band; they are not defined in terms of each other or in terms of the quantities they claim to measure. The Fe K line flux and centroid energy are extracted from spectral fits to observed photon counts, not predicted from fitted parameters. The center-determination procedure (Section 2.3) uses the 1.5-7.0 keV band to fix (x0, y0), then applies this to the 0.5-1.5 keV band; this is a methodological choice that could introduce systematic bias (as the skeptic notes), but it is not circularity — the soft-band width is still independently fitted to the soft-band image data, not defined by the hard-band fit. Self-citations (Sakai et al. 2024, Sato et al. 2018) appear only in Section 2.2 as examples of WCS-based alignment methods that the authors explicitly state they cannot use, so these citations are not load-bearing for any result. The inclination and position angle are fixed from external HST measurements (Tegkelidis et al. 2024), not from the authors' own prior work. The Fe K Gaussian width is fixed at sigma=0.04 keV from XRISM measurements (an external result), not from a self-citation. No step in the analysis reduces to its inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (9)
- Torus radius r =
varies by epoch and band, ~0.5–0.9 arcsec
- Torus width sigma_r =
varies by epoch and band, ~0.15–0.35 arcsec
- Torus normalization N_r =
varies
- Background N_0 =
varies
- RL deconvolution iterations =
30 (merged), 10 (individual/center)
- Fe K Gaussian width sigma =
0.04 keV
- N_H =
2.17e21 cm^-2
- Inclination i =
42.85 deg
- Position angle PA =
-6.24 deg
axioms (4)
- domain assumption The projected X-ray emission can be adequately modeled as an elliptical Gaussian torus with uniform emissivity
- domain assumption RL deconvolution with 30 iterations does not introduce energy-dependent artifacts in the torus width measurement
- domain assumption The internal alignment using the remnant centroid is sufficient for sub-arcsecond image registration
- domain assumption Pileup fractions of ~10% do not significantly bias the morphological measurements
read the original abstract
We present a systematic imaging and spatially resolved spectral study of SN 1987A using Chandra observations obtained between 1999 and 2025. By combining multiepoch ACIS and HETG data, we investigate the long-term evolution of the remnant in both the soft and hard X-ray bands. To characterize the radial structure, we model the projected emission with a torus profile and derive its radius and width on the image plane. We find an energy dependence in the ring morphology: while the soft and hard bands exhibit similar structures at early epochs, the soft-band emission becomes systematically broader than the hard-band emission after the early 2010s. Furthermore, when considering the radius and width together, the soft-band emission shows an inward extension, suggesting an increasing contribution from interior and/or high-latitude emission components. The flux evolution of the Fe K line is consistent with previous XMM-Newton results, and we detect its presence already in earlier epochs (~2007-2009) using combined Chandra spectra. Spatially resolved analysis further indicates that the Fe K emission is enhanced in the eastern region. These results provide a unified view of the long-term morphological and spectral evolution of SN 1987A and highlight the emergence of energy-dependent radial structure as a key feature in its late-time evolution.
Figures
Reference graph
Works this paper leans on
-
[1]
, keywords =. doi:10.1086/340591 , eprint =
-
[2]
doi:10.1088/0004-637X/692/2/1190 , eprint =
, keywords =. doi:10.1088/0004-637X/692/2/1190 , eprint =
-
[3]
doi:10.3847/1538-4357/acd9b3 , eid =. arXiv , author =:2306.13355 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/acd9b3
-
[4]
Chandra observations of SN 1987A: the soft X-ray light curve revisited
doi:10.1088/0004-637X/764/1/11 , eid =. arXiv , author =:1212.2664 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/764/1/11
-
[5]
doi:10.1086/157084 , journal =
-
[6]
The matter beyond the ring: the recent evolution of SN 1987A observed by the Hubble Space Telescope
doi:10.3847/1538-4357/ab4ff2 , eid =. arXiv , author =:1910.09582 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ab4ff2 1910
-
[7]
Chandra Observations of Supernova 1987A
Supernova 1987A: 20 Years After: Supernovae and Gamma-Ray Bursters , date-added =. doi:10.1063/1.3682881 , editor =. arXiv , author =:0704.0209 , keywords =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1063/1.3682881
-
[8]
, keywords =. doi:10.1086/505023 , eprint =
-
[9]
, keywords =. doi:10.1086/498848 , eprint =
-
[10]
, keywords =. doi:10.1086/421701 , eprint =
-
[11]
Tracing the ejecta structure of SN 1987A: Insights and diagnostics from 3D MHD simulations
doi:10.1051/0004-6361/202554862 , eid =. arXiv , author =:2504.19896 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/202554862
-
[12]
doi:10.1051/0004-6361/201834487 , eid =. arXiv , author =:1812.00021 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201834487
-
[13]
Supernova 1987A: a Template to Link Supernovae to their Remnants
doi:10.1088/0004-637X/810/2/168 , eid =. arXiv , author =:1508.02275 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/810/2/168
-
[14]
doi:10.1051/0004-6361/201220595 , eid =. arXiv , author =:1211.1220 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201220595
-
[15]
Hydrodynamic simulations unravel the progenitor-supernova-remnant connection in SN 1987A
doi:10.1051/0004-6361/201936718 , eid =. arXiv , author =:1912.03070 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201936718 1912
-
[16]
JWST NIRSpec observations of Supernova 1987A -- from the inner ejecta to the reverse shock
doi:10.3847/2041-8213/acd555 , eid =. arXiv , author =:2302.03576 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/2041-8213/acd555 2041
-
[17]
Evolution of the Radio Remnant of Supernova 1987A: Morphological Changes from Day 7000
doi:10.1088/0004-637X/777/2/131 , eid =. arXiv , author =:1310.3913 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/777/2/131
-
[18]
doi:10.3847/1538-4357/ac052d , eid =. arXiv , author =:2103.02612 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ac052d
-
[19]
High Resolution X-Ray Spectroscopy of SN 1987A: Monitoring with XMM-Newton
doi:10.1051/0004-6361/200913317 , eid =. arXiv , author =:1002.1865 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/200913317
-
[20]
doi:10.3847/1538-4357/adaea7 , eid =. arXiv , author =:2501.18090 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/adaea7
-
[21]
doi:10.1088/0004-637X/706/1/L100 , eprint =
, keywords =. doi:10.1088/0004-637X/706/1/L100 , eprint =
-
[22]
doi:10.1088/0004-637X/703/2/1752 , eprint =
, keywords =. doi:10.1088/0004-637X/703/2/1752 , eprint =
-
[23]
, keywords =. doi:10.1086/321338 , eprint =
-
[24]
New Century of X-ray Astronomy , date-added =
-
[25]
, keywords =. doi:10.1086/374967 , eprint =
-
[26]
doi:10.1093/pasj/psaf056 , eprint =
, keywords =. doi:10.1093/pasj/psaf056 , eprint =
-
[27]
, keywords =. doi:10.1086/430396 , eprint =
-
[28]
, keywords =. doi:10.1086/338492 , eprint =
-
[29]
doi:10.1086/309714 , journal =
-
[30]
doi:10.1086/176339 , journal =
-
[31]
, keywords =. doi:10.1086/317271 , eprint =
-
[32]
CASA, the Common Astronomy Software Applications for Radio Astronomy
doi:10.1088/1538-3873/ac9642 , eid =. arXiv , author =:2210.02276 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/1538-3873/ac9642
-
[33]
doi:10.1093/mnras/stae1032 , eprint =
, keywords =. doi:10.1093/mnras/stae1032 , eprint =
-
[34]
arXiv , author =:2511.07394 , journal =
doi:10.1051/0004-6361/202555479 , eid =. arXiv , author =:2511.07394 , journal =
-
[35]
doi:10.1146/annurev.aa.31.090193.001135 , journal =
-
[36]
doi:10.1146/annurev.aa.27.090189.003213 , journal =
-
[37]
doi:10.1146/annurev-astro-082615-105405 , journal =
-
[38]
Spectral Evolution of the X-Ray Remnant of SN 1987A: A High-Resolution $Chandra$ HETG Study
doi:10.3847/1538-4357/ac249a , eid =. arXiv , author =:2109.02881 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ac249a
-
[39]
doi:10.1086/156922 , journal =
-
[40]
doi:10.3847/1538-4357/ad812e , eid =. arXiv , author =:2410.00625 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ad812e
-
[41]
The Re-Acceleration of the Shock Wave in the Radio Remnant of SN 1987A
doi:10.3847/1538-4357/aae261 , eid =. arXiv , author =:1809.02364 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/aae261
-
[42]
The Post-impact Evolution of the X-ray Emitting Gas in SNR 1987A Viewed by XMM-Newton
doi:10.3847/1538-4357/ac033d , eid =. arXiv , author =:2103.03844 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ac033d
-
[43]
doi:10.1117/12.926937 , editor =
Space Telescopes and Instrumentation 2012: Ultraviolet to Gamma Ray , date-added =. doi:10.1117/12.926937 , editor =
-
[44]
Chandra Observes the End of an Era SN 1987A
doi:10.3847/0004-637X/829/1/40 , eid =. arXiv , author =:1608.02160 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/0004-637x/829/1/40
-
[45]
Latest Evolution of the X-Ray Remnant of SN 1987A: Beyond the Inner Ring
doi:10.3847/1538-4357/ad3800 , eid =. arXiv , author =:2403.07990 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ad3800
-
[46]
Jingqiang Li and Joel H. Kastner and Gregory Y. Prigozhin and Norbert S. Schulz and Eric D. Feigelson and Konstantin V. Getman , date-added =. Chandra ACIS Subpixel Event Repositioning: Further Refinements and Comparison between Backside- and Frontside-illuminated X-Ray CCDs , url =. 2004 , bdsk-url-1 =. doi:10.1086/421866 , journal =
-
[47]
Astronomical Data Analysis Software and Systems V , date-added =
-
[48]
Fruscione, Antonella and McDowell, Jonathan C. and Allen, Glenn E. and Brickhouse, Nancy S. and Burke, Douglas J. and Davis, John E. and Durham, Nick and Elvis, Martin and Galle, Elizabeth C. and Harris, Daniel E. and Huenemoerder, David P. and Houck, John C. and Ishibashi, Bish and Karovska, Margarita and Nicastro, Fabrizio and Noble, Michael S. and Nowa...
-
[49]
Journal of the Optical Society of America (1917-1983) , month = jan, number =
work page 1917
-
[50]
doi:10.1086/111605 , journal =
-
[51]
Global deceleration and inward movements of X-ray knots and rims of RCW 103
doi:10.3847/1538-4357/ad00ae , eid =. arXiv , author =:2310.04177 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ad00ae
-
[52]
doi:10.3847/1538-4357/ad739f , eid =. arXiv , author =:2410.15072 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ad739f
-
[53]
X-ray Measurements of the Particle Acceleration Properties at Inward Shocks in Cassiopeia A
doi:10.3847/1538-4357/aaa021 , eid =. arXiv , author =:1710.06992 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/aaa021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.