{"id":"1ecde371-ac16-4bfb-979f-225f1ad3ddcb","arxiv_id":"2505.07479","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":7,"one_line_summary":"XRISM's high-resolution spectra of SN 1987A show 1500-1700 km/s Doppler-broadened Si, S, and Fe lines, indicating that the observed X-rays come from kinematic, non-metal-rich shocked ejecta.","lead":"New X-ray spectra from the XRISM satellite show that the hot plasma in SN 1987A now glows mostly from outer ejecta gas, not just the ring of material thrown off before the explosion. The broadened atomic lines reveal fast bulk motion in the ejecta, giving the clearest spectral sign yet that the reverse shock has begun plowing into the exploded star's outer layers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section 6.3 asserts electron-ion temperature equilibrium without providing a Coulomb equilibration calculation; if ions are hotter than electrons, thermal broadening could be a substantial part of the 1500-1700 km/s line widths.","rationale":"The reader's weakest_assumption correctly identifies the Coulomb-equilibration premise in Section 6.3 as the weakest link. My stress-test agrees that this is the single most load-bearing concern: the entire kinematic interpretation and the reverse-shock location argument rest on negligible thermal broadening, which in turn rests on T_i ≈ T_e. The paper provides no quantitative support for this equilibration claim, and the standard Spitzer timescale can be comparable to or longer than the remnant age depending on the actual density of the shocked ejecta. Other issues raised by the reader (e.g., 'first clear evidence' wording, Fe abundance excluding 1.0 LMC) are secondary wording or consistency matters that do not threaten the physical interpretation. The concern does not by itself overturn the paper, because the MHD simulation's a priori agreement with the line profiles is a genuine independent support; however, it does warrant keeping the verdict at CONDITIONAL because the central conclusion would need a dedicated check of the equilibrium assumption. Therefore I do not change the reader's verdict, but I emphasize that the equilibration statement is not just a weak justification but a quantitatively untested premise. My proposed concrete test — either a direct Spitzer-timescale calculation or a spectral fit allowing T_i ≠ T_e — would settle the issue.","tokens_in":30505,"tokens_out":18634,"duration_ms":181966,"concrete_test":"Derive the electron density from the best-fit emission measure normalization of the pshock_wide2 model in Table 2, then compute the Spitzer electron-ion equilibration time for Fe, Si, and S at kT = 2.84 keV and compare it with the 37-year remnant age. If t_eq is longer than the age, fit the six measured lines in Table 4 with a model that adds ion-temperature-dependent thermal broadening (σ ∝ sqrt(kT_i/m_i)) as a free parameter alongside the kinematic Gaussian width; if the fit requires T_i >> T_e, the kinematic interpretation of the line widths is weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the observed 1500-1700 km/s line widths (Tables 2 and 4) are purely kinematic depends on the statement in Section 6.3: 'the electron and ion temperatures are expected to have reached equilibrium at τ ∼ 10^10-11 cm^-3 s, even for Fe.' No derivation or reference supports this timescale. The bvpshock model used throughout has a single temperature parameter (k T ≈ 2.84 keV) and a separate Gaussian Doppler broadening; it cannot distinguish thermal ion broadening from bulk motion. If ions are heated preferentially at the reverse shock and have not equilibrated with electrons, the thermal width for Fe at a few keV already approaches ~100 km/s, and for T_i substantially above T_e — as is common in collisionless shocks — the Fe, Si, and S lines could acquire thermal widths of hundreds of km/s, comparable to the measured widths. The LMC-like abundances constrain composition, not the electron-ion equilibration state, which depends on the product of density and time via the Spitzer relaxation rate. The paper's supporting MHD simulation (Sapienza et al. 2024a; Orlando et al. 2020) likely assumes a single plasma temperature, so its successful reproduction of the line profiles does not independently rule out a hot ion population. The claim is plausible, but the equilibrium assumption is load-bearing and is asserted rather than established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first high-resolution X-ray spectroscopic analysis of SN 1987A with XRISM Resolve, using a 290.5 ks observation from June 2024. The 1.7–10 keV spectrum is well reproduced by a single plane-parallel shock (bvpshock) model with kT = 2.84 keV, ionization parameter τ = 2.64×10^11 s cm^-3, and abundances consistent with LMC values for Si, S, Ar, Ca, and Fe. Doppler-broadened line widths of roughly 1,500–1,700 km/s are measured, and the authors argue that thermal broadening is negligible in the LMC-abundance plasma, so the widths represent bulk kinematic motion of the ejecta. The 3D MHD simulation of Orlando et al. (2020) reproduces the observed spectrum and line profiles, with ejecta dominating the line wings and circumstellar material the line cores, supporting the shocked-ejecta interpretation. The paper also reports a 90% upper limit on pulsar wind nebula emission and a 1σ upper limit on the 44Sc K line, corresponding to an initial 44Ti mass below about 2.6×10^-4 M_sun.","tokens_in":30599,"tokens_out":18312,"duration_ms":166842,"significance":"If correct, this is a major observational result: it provides the first direct, high-resolution evidence that the reverse shock in SN 1987A is now heating the outer, LMC-composition ejecta envelope, and that the inner metal-rich layers have not yet been reached. The analysis is careful and reproducible in structure: the non-X-ray background is modeled and its normalization checked, LMC and solar-flare contamination are quantified and shown to be negligible, three energy bands are fitted independently and found consistent, and the results are robust to different abundance priors. The comparison with the previously published 3D MHD simulation is a genuine, parameter-free prediction that strengthens the ejecta interpretation. The upper limits on the pulsar wind nebula and 44Ti yield are consistent with NuSTAR and INTEGRAL results and are useful additions. The conclusion is falsifiable with continued XRISM monitoring, which should reveal the transition to metal-rich ejecta.","major_comments":[{"comment":"The claim that \"the electron and ion temperatures are expected to have reached equilibrium at τ ∼ 10^10–10^11 cm^-3 s, even for Fe\" is asserted without any derivation or reference. This assumption is load-bearing because the conclusion that the measured 1,500–1,700 km/s line widths are purely kinematic rests on thermal broadening being negligible. A quantitative estimate of the Coulomb equilibration timescale for Fe and Si at kT ≈ 2.8 keV and the plausible post-shock density should be provided, along with the resulting thermal line widths if T_i were a factor of 3–10 above T_e. Even if the final conclusion remains unchanged, the paper should demonstrate this explicitly rather than stating it. The supporting MHD simulation, which presumably assumes a single plasma temperature, does not by itself rule out a hotter ion population.","section":"Section 6.3"}],"minor_comments":[{"comment":"The stated line-width range of 1,500–1,700 km/s (abstract) and 1,500–1,800 km/s (Section 6.3) appears to refer to the full-band model values in Tables 2 and 3, but the line-by-line measurements in Table 4 span roughly 800–2,100 km/s. Please clarify which quantity is being quoted to avoid apparent inconsistency.","section":"Abstract and Section 6.3"},{"comment":"The axis labels in the top row of Figure 10 are garbled in the arXiv version (e.g., \"4J9***\"); please ensure the final published figure renders correctly.","section":"Figure 10"},{"comment":"The phrase \"The observation id of this observation was\" should be \"The observation ID was\" for style consistency.","section":"Section 2"},{"comment":"When introducing the Fe XXVI lines at 6.952 and 6.973 keV, please state explicitly that these energies are in the rest frame, given that the paper elsewhere discusses Doppler shifts.","section":"Section 4.1"},{"comment":"The -1,000 km/s shift applied to the 44Sc line should be explicitly identified as a blueshift relative to the systemic velocity, since the text describes it as the sum of intrinsic and recession velocities.","section":"Section 6.5"}],"recommendation":"major_revision","confidential_remarks":"This is a strong observational paper with a clear result. The main technical caveat is the electron-ion equilibrium assumption in Section 6.3, which is stated without support and is central to the kinematic interpretation. The MHD comparison is legitimate because the simulation was published before the observation and the authors clearly state that the agreement is not from an ad-hoc fit. With a proper equilibration calculation or a suitably caveated thermal-broadening discussion, the paper would be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is the first resolved X-ray spectroscopy of SN1987A, and the main conclusion—that the 2024 X-ray emission is dominated by non-metal-rich shocked ejecta with bulk velocities of 1500–1700 km/s—survives scrutiny. The analysis is careful, and the paper has a real predictive success behind it.\n\nWhat's genuinely new: resolved line profiles for Si, S, Ar, and Fe, measured individually rather than inferred from broad-band fits. The multi-band fitting is consistent, the abundance-prior checks don't move the results, and the NXB and solar-flare contamination are quantified. The strongest point is the comparison with the Orlando et al. (2020) MHD simulation: Sapienza et al. (2024a) predicted before the observation that Resolve would see broadened lines from shocked ejecta, and the synthetic spectra reproduce the observed profiles without refitting. That is a genuine prediction, not a postdictive fit. The shared co-authors with the simulation team don't bother me because the prediction was on record first.\n\nSoft spots, in approximate order. (1) Section 6.3 asserts electron–ion equilibrium at τ ~ 10^10–10^11 cm^-3 s without showing the Coulomb calculation or giving a reference. I'd want that added, but I don't think it's load-bearing: even at T_i = 30 keV, thermal broadening for Si and Fe is a few hundred km/s, far below the measured 1500–1700 km/s. The claim is plausible; it just needs support. (2) The abstract says abundances \"align with the LMC value,\" but Table 2 gives Fe = 0.72(+0.15/–0.11) LMC, whose 90% interval excludes 1.0. That's a wording mismatch, not a result-breaker. (3) \"First clear evidence of X-ray emission from shocked ejecta\" overstates what Ravi et al. (2024) and Sun et al. (2025) already argued—both are cited in the same section. \"First direct kinematic evidence\" would be accurate and still important.\n\nWho this is for: anyone working on SN 1987A, reverse shocks, or XRISM science. It deserves a serious referee and should be published after minor revision. I'd bring it to reading group too, because the prediction-then-confirmation arc is a good example of how to use MHD simulations honestly.","headline":"First resolved X-ray line profiles in SN 1987A make a strong, mostly well-supported case that the 2024 X-ray emission is dominated by kinematically broadened, non-metal-rich shocked ejecta; the remaining issues are wording and one underargued assumption.","tokens_in":31425,"tokens_out":3897,"would_cite":true,"duration_ms":40474,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The 1.7–10 keV spectrum of SN 1987A is dominated by reverse-shocked outer ejecta with LMC-like abundances, and its 1,500–1,700 km/s line widths are bulk motion, not thermal broadening.","keywords":["SN 1987A","supernova remnants","X-ray spectroscopy","XRISM","shock-heated ejecta","reverse shock","Doppler line broadening","plasma diagnostics"],"falsifier":"Measure Doppler widths of emission lines from ions of very different masses, such as O or Ne near 16–20 amu against Fe or Ni near 56–58 amu, in the same shocked plasma. Kinematic broadening predicts nearly the same velocity for all ions, while thermal broadening predicts widths proportional to $m^{-1/2}$, so Fe lines would be only about 0.6 times as wide as O lines; a high-signal spectrum deciding between those two scaling laws would settle whether the 1,500–1,700 km/s widths are bulk motion. The current XRISM Fe and Si measurements already point to Fe being as wide as or wider than Si, but with larger error bars than this test requires.","tokens_in":30037,"feed_emoji":"🔭","tokens_out":8172,"duration_ms":80718,"temperature":0.7,"pith_summary":"Using the Resolve microcalorimeter on XRISM, this paper measures the 1.7–10 keV spectrum of SN 1987A 37.3 years after explosion. It finds that a single plane-parallel shock component at $kT \\approx 2.84$ keV with Large Magellanic Cloud abundances fits the whole band, and that a 3D magnetohydrodynamic simulation of the remnant reproduces the spectrum with the dominant contribution coming from shocked ejecta rather than circumstellar ring material. The observed Si, S, Ar, and Fe lines are Doppler-broadened to 1,500–1,700 km/s, far larger than thermal broadening would be for a $\\sim$2.8 keV plasma with LMC composition, so the paper concludes the line widths trace the bulk expansion of the ejecta. The implication is that the reverse shock is currently heating the outer, unenriched envelope of the ejecta and has not yet reached the metal-rich inner layers. The same observation also places upper limits on a pulsar wind nebula and on the initial $^{44}$Ti yield that are consistent with earlier missions.","feed_headline":"SN 1987A's broad X-ray lines are ejecta motion, not heat","feed_subtitle":"XRISM's Resolve spectra trace 1,500–1,700 km/s Doppler broadening to LMC-composition, reverse-shocked outer ejecta.","key_machinery":"The load-bearing tool is the plane-parallel shock plasma model (bvpshock in XSPEC), which parameterizes the post-shock plasma by temperature $kT$ and an upper ionization timescale $\\tau_u$ and assigns each atomic line a Doppler broadening velocity. The paper fits this model over three sub-bands and the full 1.7–10 keV band, checks internal consistency, and additionally compares the observed spectrum with a synthetic spectrum built from a 3D MHD simulation of SN 1987A, which separates the contributions of circumstellar matter and ejecta to the line cores and wings. The diagnostic that turns measured widths into bulk velocities is the comparison with thermal broadening, using the known scaling of thermal line width with ion mass and the measured energy dependence $\\alpha \\simeq 1.25$ of the broadening. The observation itself is made with the Resolve microcalorimeter, whose roughly 4.5 eV resolution is what makes the line widths measurable.","core_discovery":"The paper's central claim is that in June 2024 the X-rays from SN 1987A were emitted by non-metal-rich, shock-heated ejecta, and that the enhanced width of the X-ray lines is kinematic. The full 1.7–10 keV Resolve spectrum is described by one plane-parallel shock plasma with $kT = 2.84^{+0.09}_{-0.08}$ keV and an ionization parameter $\\tau_u = 2.64^{+0.58}_{-0.45} \\times 10^{11}$ s cm$^{-3}$, and the metal abundances agree with LMC values rather than with metal-enriched ejecta. Doppler widths of Si XIII, Si XIV, S XV, S XVI, Ar XVII, and Fe XXV lines correspond to velocities of 1,500–1,700 km/s. Because at this temperature and composition the electrons and ions should be Coulomb-equilibrated, thermal broadening is negligible, so the line widths must reflect bulk motion of the ejecta. The paper states that these are the first clear evidence of X-ray emission from shocked ejecta in SN 1987A, and that the reverse shock has not reached the inner metal-rich region of the ejecta.","pith_inferences":["A direct test of the kinematic interpretation is to compare line widths of low- and high-mass ions: bulk motion gives nearly equal velocities, while thermal broadening gives widths $\\propto m^{-1/2}$; a longer XRISM exposure or an imaging calorimeter could make that test decisive.","If the reverse shock is still in the LMC-composition envelope, SN 1987A offers a rare clean laboratory for measuring the composition and velocity structure of the outer layers of a core-collapse supernova before metal-rich ejecta arrive, and the transition time will calibrate simulations of mixing and reverse-shock propagation.","The apparent excess in the Fe line wings may be an early sign of the onset of metal-rich ejecta emission, in which case the non-metal-rich phase could last only a few more years and repeat observations should see the line centers and widths change.","The same approach of fitting a single shock model and comparing with 3D MHD predictions could be applied to other young supernova remnants with high-resolution spectra, where line-width ratios between different ions would discriminate thermal from kinematic broadening even when equilibration assumptions are less secure."],"forward_implications":["If the central claim is right, continued monitoring should show the X-ray spectrum gradually becoming metal-rich as the reverse shock reaches the inner ejecta, and the timing of that transition would constrain the density and asymmetry of the progenitor's envelope.","The measured 1,500–1,700 km/s velocities are slower than the early 3,000–4,000 km/s expansion of the blast wave, implying that the X-ray-emitting ejecta have been decelerated while the forward shock has left the equatorial ring.","Because the Fe XXV line appears broader than the model predicts and its wings are under-reproduced, the data hint that some Fe-rich plasma may already be contributing or has its own velocity component, a point future MHD comparisons should settle.","The 44Sc K line upper limit of $1.0 \\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ translates to an initial $^{44}$Ti mass of roughly $2 \\times 10^{-4} M_\\odot$, consistent with NuSTAR and INTEGRAL, so the finding introduces no new nucleosynthesis discrepancy.","The pulsar wind nebula upper limit of $4.3 \\times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV band aligns with NuSTAR, leaving the presence of a compact object in SN 1987A unconfirmed."],"supporting_citations":[{"why":"Supplies the 3D MHD simulation whose synthetic spectra reproduce the Resolve data and separate the ejecta and circumstellar-medium contributions.","marker":"Orlando et al. (2020)"},{"why":"Predicts the XRISM-Resolve spectrum from that simulation and anticipates broad line emission from newly shocked ejecta.","marker":"Sapienza et al. (2024a)"},{"why":"Provides the plane-parallel shock plasma model (bvpshock) used for all spectral fits.","marker":"Borkowski et al. (2001)"},{"why":"Supplies the adopted lighter-element abundances and the earlier measurement of a line-broadening energy power of approximately 2 to compare against.","marker":"Sun et al. (2021)"},{"why":"Provides an alternative abundance set and the long-term Chandra temperature trend that the measured $kT$ continues.","marker":"Ravi et al. (2024)"},{"why":"Defines the LMC abundance reference against which the plasma abundances are judged non-metal-rich.","marker":"Russell & Dopita (1992)"},{"why":"Provides the NuSTAR $^{44}$Ti flux used to predict the expected $^{44}$Sc line intensity at the XRISM epoch.","marker":"Boggs et al. (2015)"},{"why":"Supplies the pulsar wind nebula model and flux limits that the XRISM upper limits are compared with.","marker":"Greco et al. (2022)"},{"why":"Measures thermal broadening below 450 km/s in the earlier CSM-dominated phase, used as the benchmark for why thermal broadening cannot explain the new widths.","marker":"Miceli et al. (2019)"}],"fun_headline_variants":["SN 1987A: Ejecta motion, not heat, widens X-ray lines","XRISM: SN 1987A's X-rays from unmixed outer ejecta","SN 1987A's X-ray spectrum: Kinematics dominate thermal","LMC-like metals: SN 1987A's reverse shock misses core","SN 1987A: 1,700 km/s ejecta seen in X-ray line widths"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument's load-bearing premise is that electrons and ions in the non-metal-rich shocked plasma are Coulomb-equilibrated, so that at $kT \\approx 2.8$ keV the true thermal line broadening is negligible; if the collisionless shock instead heats ions preferentially and they remain much hotter than the electrons, part or all of the 1,500–1,700 km/s widths could be thermal and the bulk-motion conclusion would weaken.","fun_headline_variants_meta":{"raw":{"variants":["SN 1987A: Ejecta motion, not heat, widens X-ray lines","XRISM: SN 1987A's X-rays from unmixed outer ejecta","SN 1987A's X-ray spectrum: Kinematics dominate thermal","LMC-like metals: SN 1987A's reverse shock misses core","SN 1987A: 1,700 km/s ejecta seen in X-ray line widths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00037,"raw_usage":{"total_tokens":2129,"prompt_tokens":1239,"completion_tokens":890,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":855,"completion_tokens_details":{"reasoning_tokens":779}},"tokens_in":855,"tokens_out":890,"duration_ms":8054,"temperature":1.0,"reasoning_tokens":779,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:16:56.108921+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure Doppler widths of emission lines from ions of very different masses, such as O or Ne near 16–20 amu against Fe or Ni near 56–58 amu, in the same shocked plasma. Kinematic broadening predicts nearly the same velocity for all ions, while thermal broadening predicts widths proportional to $m^{-1/2}$, so Fe lines would be only about 0.6 times as wide as O lines; a high-signal spectrum deciding between those two scaling laws would settle whether the 1,500–1,700 km/s widths are bulk motion. The current XRISM Fe and Si measurements already point to Fe being as wide as or wider than Si, but with larger error bars than this test requires.","supporting_citations":[{"cited_title":"2020, A&A, 636, A22 ,","cited_arxiv_id":null,"evidence_quote":"Supplies the 3D MHD simulation whose synthetic spectra reproduce the Resolve data and separate the ejecta and circumstellar-medium contributions."},{"cited_title":"J., Lyerly, W","cited_arxiv_id":null,"evidence_quote":"Provides the plane-parallel shock plasma model (bvpshock) used for all spectral fits."},{"cited_title":"2021, ApJ, 916, 41,","cited_arxiv_id":null,"evidence_quote":"Supplies the adopted lighter-element abundances and the earlier measurement of a line-broadening energy power of approximately 2 to compare against."},{"cited_title":"C., & Dopita, M","cited_arxiv_id":null,"evidence_quote":"Defines the LMC abundance reference against which the plasma abundances are judged non-metal-rich."},{"cited_title":"E., Harrison, F","cited_arxiv_id":null,"evidence_quote":"Provides the NuSTAR $^{44}$Ti flux used to predict the expected $^{44}$Sc line intensity at the XRISM epoch."},{"cited_title":"N., et al","cited_arxiv_id":null,"evidence_quote":"Measures thermal broadening below 450 km/s in the earlier CSM-dominated phase, used as the benchmark for why thermal broadening cannot explain the new widths."}],"review_version":1}