{"id":"e0451025-c590-4175-ba04-e44b236a8d45","arxiv_id":"2607.06338","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The superluminous Type Ia supernova SN 2022erq was powered by ejecta crashing into a ~3-solar-mass circumstellar shell, implying the progenitor's mass-loss rate escalated from ~0.04 to ~0.6 solar masses per year in its final decades.","lead":"SN 2022erq is a superluminous thermonuclear supernova whose brightness is powered mostly by the explosion slamming into a dense shell of gas lost by its companion star. The data suggest the mass-loss rate increased tenfold in the final decades before explosion, from about 0.04 to 0.6 solar masses per year.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hα and bolometric mass-loss diagnostics appear internally inconsistent when the Hα emitting radius is set to the shock radius; the claimed 0.04→0.6 M_sun/yr escalation may be an artifact of the bolometric inversion.","rationale":"The reader's CONDITIONAL verdict is appropriate: the central power-source claim (ejecta–CSM interaction) is well supported by the early Fe-dominated spectra, persistent narrow Hα, slow decline, and the Arnett-law argument. The quantitative mass-loss history, however, is less secure. The reader flagged the constant-wind and ε assumptions; I identify a sharper, testable internal inconsistency: the Hα-derived mass-loss rates at the two observed epochs are consistent with a shallow r^-2 profile when evaluated at the actual shock radii, while the bolometric inversion yields a much steeper r^-3.66 profile. If this inconsistency is confirmed by the proposed re-analysis, the abstract's central quantitative claim—the order-of-magnitude escalation—would not be independently supported by the Hα data and would depend entirely on the bolometric model. This does not overturn the paper's classification or interaction-powered interpretation, but it does require a careful reanalysis before the mass-loss history can be accepted. Hence I keep the reader's CONDITIONAL verdict rather than moving to ACCEPT or REJECT.","tokens_in":40311,"tokens_out":12647,"duration_ms":123549,"concrete_test":"Re-derive the Hα-based Mdot at τ=19 and 77 d using the physical emitting radius r_sh(t)=∫v_sh dt (≈6×10^14 cm and ≈3×10^15 cm, respectively) instead of a fixed outer radius, and overlay these as density points on Figure 14. If the points fall on the ρ∝r^-3.66 power law (Mdot≈0.2–0.6), the consistency claim survives; if they remain near Mdot≈0.04, the two probes are mutually inconsistent and the bolometric inversion (Eq. 2) must be revised—e.g., lower ε, different v_sh(t), or non-spherical/clumpy CSM—before the escalation can be accepted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline quantitative claim is the escalating mass-loss history. Section 5.2 derives Mdot≈0.04 M_sun/yr from narrow Hα at τ≈19 and 77 d using the Ofek et al. (2013) steady-wind formula with v_wind=180 km/s. Section 5.3 inverts Eq. (2) to obtain ρ∝r^-3.66, which corresponds to Mdot=4πr²ρv_wind ∝ r^-1.66: at r≈3×10^15 cm (τ≈77 d) this gives Mdot≈0.2–0.25 M_sun/yr, not 0.04. But the narrow Hα at τ≈77 d arises from unshocked CSM just ahead of the forward shock, at r_sh≈∫v_sh dt≈3×10^15 cm. If the Hα formula is evaluated at this radius, the derived Mdot is ≈0.04, consistent with ρ∝r^-2, not r^-3.66. The paper's claim that the two diagnostics are 'broadly consistent' (Section 5.3, Figure 14) appears to hold only if the Hα emission is referred to the outer CSM radius (~3.5×10^16 cm) rather than the instantaneous shock radius. Thus the Hα data, as interpreted, do not independently support the steep density profile; the escalation rests entirely on the bolometric inversion with its assumed ε≈50% and two-epoch v_sh(t) parameterization.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"SN 2022erq is presented as the earliest spectroscopically confirmed Type Ia-CSM supernova, discovered about 1.8 days after explosion and followed for about 1350 days. From multi-band photometry and spectroscopy the authors derive a peak bolometric luminosity of about 8e43 erg/s, an unusually slow post-peak decline, and an IGE-rich/IME-weak spectrum with persistent narrow H-alpha. They argue that Arnett's law would require roughly 6 Msun of 56Ni, ruling out radioactive decay as the dominant power source, and that ejecta-CSM interaction supplies most of the luminosity. Using the Ofek et al. (2013) narrow-H-alpha formula they infer a wind velocity of about 180 km/s and a mass-loss rate near 0.04 Msun/yr; inverting the bolometric light curve they derive rho proportional to r^-3.66, which they translate into an escalating mass-loss rate reaching about 0.6 Msun/yr in the final years and a total CSM mass of about 3 Msun. They conclude that the progenitor is a white dwarf with an intermediate-mass companion that underwent escalating late-stage mass loss.","tokens_in":1938,"tokens_out":1679,"duration_ms":138732,"significance":"If the quantitative reconstruction holds, the paper is a valuable addition to the rare Ia-CSM class: it has the earliest spectra, dense photometry, a clean host environment, and one of the few attempts to reconstruct a pre-explosion mass-loss history for this class. The central claim that CSM interaction, not 56Ni decay, dominates the luminosity is well supported by three independent observational arguments: the Arnett-limit inconsistency, the persistent narrow Balmer lines, and the slow post-peak decline. The extensive dataset and the machine-readable tables and spectra are concrete assets. The quantitative mass-loss history is the least secure part because it depends on a single-zone inversion with several hand-tuned parameters, and the H-alpha data do not independently corroborate the steep density profile as presented. This weakness does not undermine the classification or the interaction-dominance claim, but it does affect the headline 'escalating mass loss' result.","major_comments":[{"comment":"The claimed 'broad consistency' between the H-alpha and bolometric diagnostics is not supported when the H-alpha emission is referred to the instantaneous shock radius. The Ofek et al. formula used in Section 5.2 is the recombination luminosity of a wind outside an inner radius r; for a local mass-loss rate Mdot_eff = 4 pi r^2 rho v_wind it gives L_Halpha proportional to Mdot_eff^2 / (r v_wind^2). At tau about 77 d, the adopted v_sh(t) = 2800 t^0.15 gives r_sh about 3e15 cm, so the H-alpha point at this epoch measures Mdot_eff about 0.04 Msun/yr at that radius, i.e. rho about 1.2e-14 g/cm3. The bolometric profile rho proportional to r^-3.66 normalized to the quoted inner value (about 0.6 Msun/yr near 1.5e15 cm) implies Mdot_eff proportional to r^-1.66 and predicts Mdot_eff(3e15 cm) about 0.1-0.2 Msun/yr, a factor 3-5 above the H-alpha value. Moreover, the two H-alpha epochs (4.5e-2 and 3","section":"Sections 5.2 and 5.3"},{"comment":"The derivation of s = 3.59 from the post-peak slope alpha = -1.46 is called an 'independent check,' but it is not independent: both quantities are derived from the same bolometric light curve, and the Moriya et al. (2013) relation is an analytic approximation of the same inversion with an assumed ejecta index n about 10. It verifies only internal consistency of the power-law model. In addition, the inversion of Eq. (2) assumes a single, time-independent epsilon = 0.5. That value is calibrated from the total radiated energy divided by a canonical 1.4 foe, but the instantaneous conversion efficiency is not measured; the two-epoch v_sh(t) parameterization and the 40-330 d phase boundaries are also unquantified inputs. The resulting 0.04 to 0.6 Msun/yr escalation therefore has no formal error budget. The authors should either propagate these systematic uncertainties or present the mass-loss","section":"Section 5.3"},{"comment":"The total CSM mass and outer radius are derived with two different velocity assumptions. The H-alpha mass of about 2.4 Msun assumes a constant 180 km/s wind over about 60 yr, while the bolometric mass of about 3 Msun comes from the r^-3.66 profile; the near-agreement of the two numbers is therefore not an independent confirmation. Also, the outer CSM radius of about 3.5e16 cm is obtained from an 'ejecta velocity' of 10,000 km/s, whereas the CDS velocities used in Section 5.3 are 5,400-5,900 km/s; the latter gives r about 2e16 cm over 400 d, closer to the 1.6e16 cm bolometric extent. A single, clearly defined interaction radius should be used consistently in both diagnostics.","section":"Sections 5.2 and 5.3"}],"minor_comments":[{"comment":"The parameter beta in the Ofek et al. formula is not defined; please state its assumed value or give the reference for the adopted normalization.","section":"Section 5.2"},{"comment":"Please specify the phase range over which the L proportional to t^-1.46 fit is performed; 'post-peak' is ambiguous for a light curve with such a long, nearly flat evolution.","section":"Figure 13"},{"comment":"The SED-derived stellar metallicity log(Z/Zsun) = -2.13 (about 0.7% solar) is remarkably low, far below the gas-phase R23 value; a sentence discussing whether this result is robust to the BayeSED prior choices would be helpful.","section":"Section 5.1"},{"comment":"Minor typographical issues: 'NaiD' should read 'Na I D'; in Section 4.1, 'overlay Hbeta and Hgamma absorption' would read more clearly as 'overlaid by Hbeta and Hgamma absorption.'","section":"Section 2.3"}],"recommendation":"major_revision","confidential_remarks":"The observational case for an interacting Ia-CSM is strong and the dataset is a major asset. My concern is entirely with the quantitative mass-loss reconstruction. If the authors can either resolve the H-alpha/bolometric radius inconsistency documented in Section 5.3/Figure 14 or explicitly reframe the escalation as a model-dependent illustration with a full systematic budget, I would support publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives us the earliest spectrum of a confirmed Ia-CSM event (τ≈1.8 d) plus an exceptionally dense optical/NIR campaign that runs to ~1350 d. That dataset alone makes it a reference object for the class. The central claim—that SN 2022erq's luminosity and slow decline are powered by ejecta–CSM interaction, not 56Ni decay—holds up. The Arnett-law argument (would require ~6 M_sun of Ni), the persistent narrow Hα, the slow multi-band decline, and the IGE-rich spectra with weak IME features all point the same way. The qualitative picture of a massive, hydrogen-rich CSM around an efficient-burning thermonuclear explosion is well supported.\n\nWhat is genuinely new is the reconstructed mass-loss escalation: ~0.04 to ~0.6 M_sun/yr over the final decades, with a ~3 M_sun CSM shell. This is the kind of claim that matters for progenitor channels and for Ia cosmology. But the quantitative reconstruction is the soft spot. The bolometric inversion (Eq. 2) depends on three weakly anchored assumptions: ε≈50% from a single energy-budget argument, v_sh(t) fitted to only two epochs, and a constant v_wind=180 km/s. These are not unreasonable, but no systematic sensitivity analysis is presented. The paper's own caveat that a time-varying wind would shift the timing is honest, but it doesn't quantify how much the mass-loss history or total mass would change.\n\nThe stress-test note raises a sharper issue, and I think it lands. The Hα-derived Mdot≈0.04 M_sun/yr at τ≈77 d is evaluated at the shock radius (r≈3×10^15 cm). At that radius, the bolometrically inverted density profile (ρ∝r^-3.66) predicts Mdot≈0.2–0.25 M_sun/yr, not 0.04. The agreement shown in Figure 14 appears to rely on placing the Hα point at the outer CSM radius (~3.5×10^16 cm) rather than at the instantaneous shock radius. If that is correct, the Hα data do not independently confirm the steep density profile; the escalation rests almost entirely on the bolometric inversion. The authors call the s=3.59 from the post-peak slope an independent check, but it uses the same bolometric light curve that produced the numerical density profile, so it is a consistency check at best.\n\nThe environment analysis (young ~100 Myr host, subsolar metallicity) is plausible but comes from a single SED fit with weak constraints; it is a minor point and not load-bearing. The spectral comparison to 1991T/2007if is heuristic, and the authors are appropriately careful not to overclaim a specific subtype link.\n\nFor whom is this paper? Anyone working on Ia-CSM, superluminous SNe, or SN Ia progenitor channels. It deserves a serious referee: the observations are excellent, the central interaction claim is solid, and the mass-loss reconstruction is interesting even if it needs more work. My recommendation is to send it to peer review, with the expectation that the authors will need to (a) re-examine the Hα radius assumption, (b) provide systematic uncertainty estimates on ε, v_sh, and the density profile, and (c) soften or better defend the quantitative escalation claim.","headline":"SN 2022erq is the earliest confirmed Ia-CSM with dense multi-band data to 1350 d; the central interaction-powered claim is solid, but the headline 0.04→0.6 M_sun/yr escalation rests on bolometric inversion assumptions that need systematic sensitivity tests.","tokens_in":41496,"tokens_out":2072,"would_cite":false,"duration_ms":20553,"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 2022erq's extraordinary luminosity and slow decline are powered by ejecta–CSM interaction, not radioactive decay, and its progenitor's mass-loss rate escalated from about 0.04 to 0.6 solar masses per year in the final decades before expl","keywords":["supernovae","Type Ia supernovae","circumstellar medium","mass loss","superluminous supernovae","H-alpha emission","white dwarf progenitors","light curves"],"falsifier":"A direct measurement of the circumstellar density profile, for example from late-time radio or X-ray emission of the shock, would test the claim: if that profile does not match the steep ρ ∝ r^{-3.7} power law inferred from the bolometric light curve, the escalating-path-to-explosion interpretation would be wrong. Alternatively, detecting resolved absorption features from the wind at multiple epochs could reveal whether the wind speed was constant.","tokens_in":40206,"feed_emoji":"💥","tokens_out":4993,"duration_ms":47689,"temperature":0.7,"pith_summary":"This paper reports on SN 2022erq, a thermonuclear supernova whose peak luminosity is roughly ten times that of a normal Type Ia. It argues that this extreme brightness does not come from radioactive nickel but from the shock wave of the explosion converting kinetic energy into light as it plows through a dense hydrogen-rich shell of circumstellar material. By combining the strength of narrow hydrogen emission with the evolution of the bolometric light curve, the authors reconstruct the pre-explosion mass-loss history: the mass-loss rate rose from about 0.04 to about 0.6 solar masses per year over the final decades, building a shell of about 3 solar masses extending to about 3.5×10^16 cm. This matters because it directly implicates a white dwarf with a nondegenerate companion as the progenitor and shows that some of the brightest thermonuclear supernovae owe their luminosity to environment rather than to an unusually powerful explosion.","feed_headline":"Exploding star's wind surged 15-fold before blast","feed_subtitle":"A superluminous supernova's glow comes from shock-heated gas, not radioactive decay, changing how such events are typed.","key_machinery":"The central device is the inversion of the bolometric light curve, L = ε dEkin/dt = 2πε ρ(r) r² v³, which converts the observed luminosity into a circumstellar density profile. The shock velocity is taken from the width of the broad Hα component, parameterized as v_sh ∝ t^{0.15}, and the kinetic-to-radiative efficiency ε is set to about 50 percent, consistent with the ratio of radiated energy to the canonical ejecta kinetic energy of about 1.4 foe. Independent density estimates come from the narrow Hα flux, which probes the outer, unshocked wind. The steep power-law density profile (index ~3.6) is then translated into a mass-loss history by assuming a constant wind speed of 180 km/s, so that","core_discovery":"The central claim is that SN 2022erq is an Ia-CSM event, a thermonuclear supernova interacting with hydrogen-rich circumstellar material, and that its light curve is dominated by long-lived ejecta–circumstellar interaction rather than radioactive decay. Applying the standard peak-luminosity–nickel-mass relation would require an implausible nickel mass of about 6 solar masses, ruling out radioactive decay as the power source. Instead, the authors model the post-peak decline as the shock's kinetic energy being converted into radiation with roughly 50 percent efficiency. Inverting the bolometric light curve yields a circumstellar density profile that steeply declines with radius (ρ ∝ r^{-3.7}),","pith_inferences":["If the wind speed varied with time, the absolute timing of the mass-loss escalation would shift, so the 0.04→0.6 M_sun/yr history is model-dependent; a direct measurement of velocity stratification could resolve this.","The steep density profile could also be produced by an asymmetric or clumpy shell, which would change the mass estimate; polarimetric observations during the interaction phase could test sphericity.","The same inversion method could be applied to other Ia-CSM events to see whether escalating mass loss is common or unique to this event.","If the underlying explosion is super-Chandrasekhar, the assumed ejecta kinetic energy and efficiency would change, which would alter the absolute density scale."],"forward_implications":["If true, the brightest Type Ia-like supernovae can be powered mainly by circumstellar interaction, so luminosity-based cosmological calibrations must account for such contamination.","The derived mass-loss escalation of about a factor of 15 over decades places strong constraints on binary evolution models, favoring a brief, violent ejection episode shortly before explosion.","The young host environment (about 100 Myr) and the massive shell point to a white dwarf with an intermediate-mass companion, a channel that current population synthesis may underestimate.","The event demonstrates that early spectroscopic classification (within days) can identify Ia-CSM events, so future wide-field surveys should find more such objects.","The success of this analytic inversion suggests that bolometric light curves of interacting supernovae can be used to reconstruct progenitor mass-loss histories."],"fun_headline_variants":["Supernova light from wind shock, not radioactive decay","Stellar wind surged 15-fold before supernova blast","Rare supernova reveals escalating mass loss before explosion","Thermonuclear supernova's glow from gas collision, not decay","SN 2022erq: interaction, not radioactivity, powers light curve"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The reconstruction of the mass-loss history rests on assuming that the circumstellar wind had a constant speed of 180 km/s and was spherically symmetric; if the wind velocity evolved with time or the shell is clumpy or asymmetric, the derived escalation and total mass could change substantially.","fun_headline_variants_meta":{"raw":{"variants":["Supernova light from wind shock, not radioactive decay","Stellar wind surged 15-fold before supernova blast","Rare supernova reveals escalating mass loss before explosion","Thermonuclear supernova's glow from gas collision, not decay","SN 2022erq: interaction, not radioactivity, powers light curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000289,"raw_usage":{"total_tokens":1563,"prompt_tokens":813,"completion_tokens":750,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":665}},"tokens_in":557,"tokens_out":750,"duration_ms":9048,"temperature":1.0,"reasoning_tokens":665,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T08:16:36.759233+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the circumstellar density profile, for example from late-time radio or X-ray emission of the shock, would test the claim: if that profile does not match the steep ρ ∝ r^{-3.7} power law inferred from the bolometric light curve, the escalating-path-to-explosion interpretation would be wrong. Alternatively, detecting resolved absorption features from the wind at multiple epochs could reveal whether the wind speed was constant.","supporting_citations":[],"review_version":3}