{"id":"65aea117-b8d0-4b89-bf50-85bd75a3137a","arxiv_id":"2607.26189","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In SN 1987A's inner ejecta, HCO+ is co-spatial with CO, and its two-line RADEX mass of 3–9×10^-6 solar masses implies hydrogen was mixed into the carbon-oxygen zones.","lead":"New ALMA maps of supernova 1987A's debris reveal HCO+ (a hydrogen ion attached to CO) sitting where CO sits, carrying a few millionths of a solar mass. Finding HCO+ means hydrogen reached the carbon- and oxygen-rich core layers, evidence of mixing inside the exploded star.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"J=4-3 line, which anchors the RADEX mass, is only ~2σ after propagating the paper's own continuum error; quoted ±1.0e-20 is inconsistent with Table 2, so the derived 3–9e-6 M_sun mass range is not robust.","rationale":"The reader's weakest_assumption focuses on the H2 collisional density, which is a valid concern about the mass normalization. However, I find the more load-bearing weakness to be the significance of the J=4-3 line after proper propagation of the continuum uncertainty. This line is essential to break the T-N degeneracy in the RADEX fit. The paper reports ±1.0e-20 for the J=4-3 intensity, but Table 2 lists a continuum error of ±2.54e-20, and the MC SED fit gives ±9.2e-20. Quadrature yields σ≈2.6e-20, making the line ~2.2σ; with the MC error it is insignificant. Under these conditions the RADEX solution is not constrained, and the mass range is not robust. The H2 density assumption changes the mass by a factor of ~3 between the two tested values, but the J=4-3 problem can change it by an order of magnitude or more and is an internal inconsistency. Both concerns point to CONDITIONAL, which the reader has already assigned, so I recommend UNCHANGED with the explicit condition that the error budget be corrected and the RADEX fit re-run with proper J=4-3 significance.","tokens_in":24821,"tokens_out":13513,"duration_ms":132964,"concrete_test":"Recompute the J=4-3 line intensity and its uncertainty from the ALMA data using the same aperture as the paper, but with the continuum level and uncertainty taken from the Monte Carlo SED fit (including the ±0.92×10^-19 uncertainty) instead of the single best-fit value. Evaluate the detection significance. Then re-run the RADEX χ^2 analysis twice: (a) with the J=4-3 line treated as a 3σ upper limit, and (b) with the full propagated error bars on both lines. If the allowed range of N_HCO+ and T_kin at 1σ expands so that the derived mass spans more than an order of magnitude (e.g., <1×10^-6 or >3×10^-5 M_sun), the quoted 3–9×10^-6 M_sun central value is not robust and any 'moderate mixing' quantification is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim (HCO+ mass 3–9×10^-6 M_sun and the inferred 'moderate amount' of hydrogen mixing) rests on a two-transition RADEX analysis. The J=3-2 line alone cannot break the T_kin–N_HCO+ degeneracy; the J=4-3 line provides the necessary anchor. The paper reports J=4-3 continuum-subtracted intensity I = (5.69 ± 1.0)×10^-20 W/m^2, but Table 2 lists uncertainties of ±0.40 (calibration), ±0.24 (RMS), and ±2.54×10^-20 (continuum). Propagating in quadrature gives σ ≈ 2.6×10^-20, so the line is only ~2.2σ. Furthermore, the paper's own Monte-Carlo SED fit gives a continuum uncertainty of ±0.92×10^-19, larger than the continuum level itself, under which the line is entirely sub-significant. The adopted 30% continuum error in Table 2 contradicts the MC result. If the J=4-3 line is not a secure detection, the RADEX fit has no meaningful constraint on T_kin; the allowed N_HCO+ and hence mass become very broad or unbounded. The qualitative detection of HCO+ and its co-spatiality with CO stand, but the quantitative '3–9×10^-6 M_sun' and the implied 'moderate amount' would not be supported by the data. This is distinct from the H2 collisional-density assumption (which changes the mass by only a factor ~3 between the two tested values); the continuum problem can change the mass by an order of magnitude or more and is an internal inconsistency in the error budget.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports ALMA observations of HCO+ J=3-2 and J=4-3 emission from the ejecta of SN 1987A. It finds the J=3-2 emission to be clumpy and co-spatial with CO, with a Spearman correlation of 0.72, and less similar to SiO and Hα. Using RADEX non-LTE radiative transfer on the two line intensities, the authors derive an HCO+ column density and mass of 3–9×10^-6 M_sun (depending on the adopted H2 collisional density, 1e5 or 1e6 cm^-3). They compare this to the CO mass and conclude that the fractional abundance M_HCO+/M_CO = 3×10^-6–3×10^-4 implies that a moderate amount of hydrogen was mixed into the carbon- and oxygen-rich zones. A chemical feasibility estimate for formation via CO + H3+ yields ~1e-7 M_sun, an order of magnitude below the claimed mass unless the H2 ionization rate is raised by an order of magnitude.","tokens_in":25116,"tokens_out":9428,"duration_ms":88873,"significance":"If the quantitative result were robust, this would be an important observational constraint: HCO+ is a direct tracer of hydrogen transport into the metal-rich core of an SN remnant, and no other remnant currently offers a comparable measurement. The paper's qualitative contributions are solid and should be credited: the spatially resolved ALMA detection of HCO+, the morphological comparison with CO/SiO/Hα/H2, and the explicit attempt to connect the observed molecular gas to hydrodynamic mixing scenarios. However, the quantitative anchor for the central claim—the HCO+ mass and the inferred 'moderate' hydrogen mixing—is currently too fragile. The J=4-3 detection is marginal under the paper's own error budget, and the formation-rate estimate does not independently validate the mass. The paper is worth pursuing, but the load-bearing quantitative steps need substantial revision.","major_comments":[{"comment":"The J=4-3 line, which anchors the RADEX mass, is not a secure detection under the paper's own error budget. Table 2 lists σ_cal=0.40, σ_RMS=0.24, and σ_cont=2.54 (10^-20 W/m^2). Propagating in quadrature gives σ≈2.58, so the continuum-subtracted line (5.69×10^-20) is only ≈2.2σ, not the ±1.0 quoted in the text. Moreover, §3.1 gives the MC SED continuum uncertainty at 356 GHz as ±0.92×10^-19, roughly 90% of the continuum; if that uncertainty is used, the line is sub-significant. Since the J=3-2 line alone cannot break the T_kin–N_HCO+ degeneracy, the derived N_HCO+ and mass (3–9×10^-6 M_sun) are not robust. The error budget and the quoted line intensities must be reconciled before the central quantitative claim can be accepted.","section":"§3.4 / Table 2"},{"comment":"The formation-rate 'feasibility test' yields ~1×10^-7 M_sun, an order of magnitude below the adopted 3–9×10^-6 M_sun, and the gap is closed by assuming ζ=3×10^-16 s^-1. Because ζ is not independently constrained for SN ejecta, this does not demonstrate that the observed HCO+ can be formed; it simply adjusts a free parameter. The calculation also assumes spatially uniform CO and H3+, which §5.2 later concedes is unrealistic. Consequently, the abstract's statement that the fractional abundance 'suggests a moderate amount of hydrogen was mixed' is not quantitatively supported by this section. Either an independent constraint on ζ must be provided, or the chemistry should be presented as a qualitative plausibility argument only.","section":"§5.1"},{"comment":"The H2 collisional density input to RADEX is derived by spreading the full 6 M_sun hydrogen envelope uniformly through the ejecta volume. This presumes the very H–C/O co-location the paper aims to infer. The two tested densities (1e5 and 1e6 cm^-3) change the derived mass by only a factor ~3, but they do not sample the possibility that H2 is absent or clumpy in the HCO+ region; if the local n_H2 differs, the excitation solution and N_HCO+ shift. The paper should state explicitly that the mass is conditional on H2 being co-located with HCO+ and should quantify how an H2 filling factor changes the allowed mass. As written, the RADEX mass cannot independently support the 'moderate mixing' conclusion.","section":"§4.1"}],"minor_comments":[{"comment":"Errorc for J=4-3 is listed as 2.54×10^-20, but 30% of Cdust=10.17×10^-20 is 3.05×10^-20. Reconcile the table, the text, and the quoted total uncertainty ±1.0×10^-20.","section":"Table 2"},{"comment":"The continuum contour levels are inconsistent: the figure caption lists 4.5e-5, 6e-5, 8e-5 Jy/beam, while the text quotes 7e-5 and 9e-4 Jy/beam. Please check and unify.","section":"Fig. 1 / §3.1"},{"comment":"The J=4-3 Gaussian fit fixes the FWHM to the J=3-2 value (1906 km/s). Given the limited spectral coverage of the J=4-3 data, state explicitly how this assumption affects the integrated intensity and the subsequent RADEX fit.","section":"§3.4"},{"comment":"The paper first assumes a uniform distribution of CO and H3+ for the formation estimate and then argues that uniform mixing is unrealistic. This tension should be acknowledged earlier and the quantitative conclusions adjusted accordingly.","section":"§5.1 / §5.2"},{"comment":"The mass range 3–9×10^-6 M_sun and the fractional abundance range spanning two orders of magnitude are conditional on two fixed H2 densities and f=1. The abstract should present these as conditional estimates, not as a measured range.","section":"Abstract / §4.2"}],"recommendation":"major_revision","confidential_remarks":"The qualitative detection and morphological analysis are sound and likely publishable, but the quantitative mass and the 'moderate hydrogen mixing' conclusion rest on the marginal J=4-3 detection and an inconsistent error budget. These are fixable by reanalysis and reframing. If the J=4-3 line cannot be established at higher significance, the paper should present an upper limit on the HCO+ mass and explicitly temper the mixing conclusion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — worth a look, but with eyes open. The paper gives the first resolved HCO+ J=3-2 map of SN 1987A's ejecta, plus a J=4-3 measurement, and argues that HCO+ requires hydrogen mixed into the C/O zones. The qualitative case is solid: the J=3-2 emission is clearly detected, clumpy, co-spatial with CO (Spearman 0.72), and sits where H-alpha is fainter. That is a genuine new observational result, and it does imply some hydrogen transport into the metal-rich layers, given HCO+ needs H, C and O in the same gas.\n\nThe quantitative anchor is the problem. The RADEX mass of 3–9e-6 M_sun rests on the J=4-3 line, but that line's significance collapses under the paper's own uncertainties. Table 2 gives continuum error 2.54e-20 for J=4-3, calibration 0.40, RMS 0.24; combining them gives ~2.6e-20, so the 5.69e-20 line is ~2 sigma. The text quotes ±1.0e-20, inconsistent with the table. The MC SED fit gives an even larger continuum uncertainty (0.92e-19). So the two-line RADEX constraint on T_kin and column is not secure; the mass range is likely much broader and not robust.\n\nAlso note the H2 collisional density in RADEX is derived by spreading the whole 6 M_sun hydrogen envelope uniformly through the ejecta. That assumes hydrogen is everywhere, which is the conclusion. Switching 1e6 vs 1e5 cm^-3 changes the mass by only a factor ~3, but the deeper problem is unquantified: if hydrogen is only in clumps, the local density and excitation change. The chemistry feasibility test is order-of-magnitude and needs a factor-ten boost in ionization rate to match observations, so the 'moderate mixing' language overreaches.\n\nNone of this kills the main point. The data show HCO+ exists where CO is, and that demands some hydrogen mixing. But the specific mass and mixing fraction should be treated as illustrative, not measured. The paper needs a corrected error budget and a clearer statement of assumptions before the numbers go into the literature.\n\nWho's it for? SN remnant and astrochemistry people. The resolved map and the CO-HCO+ correlation will be useful benchmarks. I'd send it to a referee, but the referee should insist on fixing the error propagation and softening the quantitative claims.","headline":"First resolved HCO+ map in SN 1987A is a real advance, but the mass and mixing conclusions lean on a J=4-3 detection that the paper's own error budget makes marginal.","tokens_in":25969,"tokens_out":2881,"would_cite":true,"duration_ms":27549,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"SN 1987A's HCO+ emission shows hydrogen was mixed into the remnant's carbon–oxygen core.","keywords":["HCO+","SN 1987A","supernova remnants","molecular ejecta","chemical mixing","astrochemistry","ALMA","radiative transfer"],"falsifier":"High-spatial-resolution maps of H2 emission in the same velocity channels as the HCO+ clumps: if the H2 there is far below 10^5 cm^-3, or absent, the adopted collision densities collapse and the inferred HCO+ mass and mixing conclusion no longer hold. A cleaner test would be detection of the HCO+ J=1–0 line: its intensity, combined with J=3–2 and J=4–3, would break the column-density/temperature degeneracy and settle the mass without fixing the collision density.","tokens_in":24493,"feed_emoji":"🔭","tokens_out":4797,"duration_ms":41273,"temperature":0.7,"pith_summary":"SN 1987A is a rare supernova remnant where HCO+ is detected in the ejecta itself, and this paper uses that fact to read the explosion's mixing history. The authors establish that the HCO+ J=3–2 emission is co-spatial with CO, quantify the HCO+ mass at 3–9 × 10^-6 solar masses from two rotational transitions, and argue that the abundance ratio to CO can only be reached if hydrogen from the star's envelope was transported inward to the carbon- and oxygen-rich nuclear zones. This makes HCO+ a tracer of mixing that other remnant observations do not provide. A simple formation-rate estimate falls short of the observed mass unless the ionisation rate is higher than the canonical interstellar value or additional hydrogen sits with the CO — either way, hydrogen must be co-located with CO.","feed_headline":"Mixing pulled hydrogen into SN 1987A's metal core","feed_subtitle":"ALMA maps of the ion put 3–9 millionths of a solar mass of HCO+ in the ejecta, tracing how deep hydrogen sank.","key_machinery":"The load-bearing object is HCO+ itself, a molecular ion that requires carbon, oxygen, and hydrogen to meet: the proposed formation route CO + H3+ → HCO+ + H2. Two ALMA lines (J=3–2 and J=4–3) feed a non-local-thermodynamic-equilibrium radiative-transfer calculation that turns observed intensities into column density and kinetic temperature, with the H2 collision-partner density fixed at 10^6 or 10^5 cm^-3. Spatial comparison with CO, SiO, Hα, and H2 isolates where the ion forms, and a timescale-integrated formation-rate estimate tests whether the reaction can account for the observed mass. The electron destruction channel (dissociative recombination) sets the condition that HCO+ survives onl","core_discovery":"The central claim is that HCO+ in SN 1987A exists because the ejecta are not radially stratified: hydrogen reached the CO-rich nuclear zones before and during the explosion. High-resolution ALMA maps show the J=3–2 HCO+ emitting region overlaps the CO J=2–1 emission (Spearman correlation 0.72), while the brighter HCO+ peaks sit apart from Hα and H2, placing HCO+ in compact, mildly ionised gas inside the ejecta. Using the J=3–2 and J=4–3 line intensities with a non-LTE radiative-transfer calculation, the paper derives column density and kinetic temperature for two assumed H2 collision densities, giving an HCO+ mass of 3–9 × 10^-6 solar masses and a fractional abundance relative to CO of 3×10^","pith_inferences":["If the true H2 distribution is clumpy rather than uniform, the density used in the radiative-transfer calculation could be locally higher or lower; high-resolution H2 maps would test whether the assumed 10^5–10^6 cm^-3 densities coincide with the HCO+ clumps.","The same technique applied to other young core-collapse remnants with detectable molecular ejecta could show whether deep hydrogen mixing is a generic feature of the explosion mechanism or peculiar to SN 1987A's blue supergiant progenitor.","A full chemical network including CH+, OH+, H2O, and CO+ would likely change the inferred required hydrogen fraction; the paper's single-reaction estimate is deliberately minimal, so the mixing conclusion is more robust than the specific mass of hydrogen inferred."],"forward_implications":["HCO+ becomes a quantitative tracer of hydrogen transport into metal-rich ejecta, giving models of Rayleigh–Taylor and smaller-scale mixing a new observable to match.","The HCO+/CO ratio of 3×10^-6 to 3×10^-4 sets a floor on how much hydrogen must be mixed into the C/O zones; reproducing it requires either an elevated H2 ionisation rate (about 3×10^-16 s^-1) or additional co-located H2.","Because HCO+ tracks CO spatially and forms from CO, CO maps can be used to predict where HCO+ emission should appear in other remnants.","The non-coincidence of HCO+ with Hα and H2 implies HCO+ marks shielded, low-ionisation gas and supports UV irradiation, not electron collisions, as the H2 excitation mechanism in SN 1987A's ejecta."],"fun_headline_variants":["SN 1987A's HCO+ signals hydrogen mixing into metal core","ALMA spots HCO+ in SN 1987A, mapping ejecta mixing","Hydrogen reached CO zones in SN 1987A, HCO+ shows","HCO+ mass in SN 1987A traces deep hydrogen mixing","ALMA links HCO+ to hydrogen mixing in SN 1987A's core"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The HCO+ mass rests on an assumed H2 collisional-partner density that is obtained by spreading the star's entire 6-solar-mass hydrogen envelope uniformly through the ejecta volume — effectively assuming hydrogen is already everywhere, which is exactly the mixing the paper sets out to infer.","fun_headline_variants_meta":{"raw":{"variants":["SN 1987A's HCO+ signals hydrogen mixing into metal core","ALMA spots HCO+ in SN 1987A, mapping ejecta mixing","Hydrogen reached CO zones in SN 1987A, HCO+ shows","HCO+ mass in SN 1987A traces deep hydrogen mixing","ALMA links HCO+ to hydrogen mixing in SN 1987A's core"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000354,"raw_usage":{"total_tokens":1807,"prompt_tokens":836,"completion_tokens":971,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":881}},"tokens_in":580,"tokens_out":971,"duration_ms":8076,"temperature":1.0,"reasoning_tokens":881,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T00:33:49.305619+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"High-spatial-resolution maps of H2 emission in the same velocity channels as the HCO+ clumps: if the H2 there is far below 10^5 cm^-3, or absent, the adopted collision densities collapse and the inferred HCO+ mass and mixing conclusion no longer hold. A cleaner test would be detection of the HCO+ J=1–0 line: its intensity, combined with J=3–2 and J=4–3, would break the column-density/temperature degeneracy and settle the mass without fixing the collision density.","supporting_citations":[],"review_version":1}