{"id":"aad0f3c7-c327-46e2-b7e5-3f1fb38e293b","arxiv_id":"2608.11978","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"SNR 0509-67.5 shows a 1000 km/s bulk Doppler shift and a flattened ejecta edge, which the authors model as a companion shadow to infer a double-white-dwarf merger progenitor.","lead":"Deep MUSE observations of supernova remnant SNR 0509-67.5 reveal a [Fe XIV] emission line split into two Doppler components across the whole remnant, with a net redshift near 1000 km/s and a flattened northeastern edge. The authors interpret these as the orbital motion of a white dwarf and the shadow of a companion white dwarf, pointing to a double-degenerate merger as the supernova's progenitor.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Before accepting the double-degenerate merger claim, settle the internal 5303/5308 Å rest-wavelength conflict: a 5 Å error shifts the central -1000 km/s bulk velocity by ~280 km/s and breaks the claimed agreement with the inferred 990 km/s primary orbital velocity.","rationale":"I agree with the reader that the companion-shadow identification lacks a quantitative null test and depends on manually selected tracers. However, the even more load-bearing link in the chain is the bulk-velocity measurement itself: it is the independent observable used to anchor the orbital-velocity interpretation and to check the cone-derived solution. The manuscript contains an internal calibration conflict (5303 vs 5308 Å) that is not resolved. A 5 Å error is ~283 km/s, much larger than the quoted error bar. Since the paper's own figure caption undermines the zero-point, the central claim is not yet supported to the stated precision. This is a missing support that should have been caught before claiming agreement with the cone-derived orbital velocity. The proposed check is straightforward: recompute velocities with the correct rest wavelength from the public data. The reader's concern about the shadow remains valid and would keep the verdict conditional even if the calibration is fixed.","tokens_in":13270,"tokens_out":11222,"duration_ms":122578,"concrete_test":"Refit the [Fe XIV] line centroids from the archival MUSE cube using the laboratory rest wavelength of [Fe XIV] (5302.86 Å in air) and an explicit, stated LMC systemic velocity, reproducing the single/double-peak maps of Fig. 1. If the median bulk velocity changes by more than ~100 km/s from -1000 km/s, then the claimed agreement in Section 3.3 between the observed bulk velocity and the cone-derived 990 km/s primary orbital velocity is not established; the manuscript should be revised or the double-degenerate claim weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central observable is the bulk Doppler velocity of -1000±60 km/s (Section 3.1), interpreted as the line-of-sight orbital velocity of the primary white dwarf and later shown to agree with the 990 km/s primary velocity derived from the cone fit and Eggleton relation (Section 3.3). That observable is not calibrated securely. The abstract and methods state the [Fe XIV] line as λ5303, but Fig. 1's caption states that the rest wavelength used was 5308 Å. A 5 Å mismatch corresponds to ~283 km/s at 5303 Å, several times the quoted 60 km/s uncertainty and large enough to erase the agreement with the 990 km/s value. The quoted error budget also omits rest-wavelength and systemic-velocity (LMC ~293 km/s) uncertainties, so the 60 km/s precision is not a measure of accuracy. If 5308 Å was actually used, the bulk velocity is partly an artifact of the adopted zero-point; if 5303 Å was used, the caption is misleading. Either way, the manuscript does not currently demonstrate that the measured offset is real. Without a secure bulk velocity, the double-degenerate interpretation loses its quantitative anchor.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes deep MUSE integral-field observations of SNR 0509-67.5, using a supervised neural network to classify [Fe XIV] λ5303 line profiles as single- or double-peaked and to fit Gaussian components. From the resulting Doppler maps the authors measure a bulk Doppler velocity of -1000±60 km/s, which they interpret as the line-of-sight component of the primary white dwarf's orbital velocity in a double-degenerate merger. They also identify a flattened north-eastern edge in the red- and blue-shifted ejecta maps, interpret it as the shadow of a companion white dwarf, fit a cone to this feature, and use the Eggleton Roche-lobe relation to derive a secondary mass of ~0.6 Msun, an orbital separation of ~0.042 Rsun, and an orbital velocity of ~1730 km/s for the secondary, with a primary orbital velocity of ~990 km/s that is claimed to agree with the bulk Doppler shift. The paper concludes that SNR 0509-67.5 originated from a double-degenerate merger progenitor.","tokens_in":13431,"tokens_out":8235,"duration_ms":87292,"significance":"If correct, this would be one of the first direct kinematic and morphological constraints on the progenitor system of a Type Ia supernova remnant, providing a rare dynamical link between the explosion and a specific double-degenerate channel. The analysis has genuine strengths: the Doppler-split [Fe XIV] emission across the whole remnant is a novel dataset; the bulk velocity measurement and the cone-derived orbital velocity are independent observables, so their consistency is a real check rather than a constructed agreement; and the authors make their machine-learning code and reduced data publicly available. The paper also connects the observed morphology to forward models of double-detonation remnants, which is a useful physical template. However, the central quantitative claims rest on two points that are not yet secured: the rest-wavelength calibration of the bulk velocity, and the identification of the flattened edge as a companion shadow in the absence of any quantitative null test. These issues affect the derived binary parameters and the final double-degenerate conclusion.","major_comments":[{"comment":"The manuscript uses inconsistent rest wavelengths for [Fe XIV]: the abstract and §2.2 state λ5303, while the Fig. 1 caption states that the Doppler maps were calculated using a rest wavelength of 5308 Å. A 5 Å error corresponds to c×5/5303 ≈ 283 km/s, several times the quoted ±60 km/s uncertainty. If 5308 Å was actually used, the bulk velocity after LMC systemic correction would be shifted by roughly 283 km/s relative to the λ5303-based value, making the claimed agreement with the cone-derived primary orbital velocity (990 +310/−220 km/s) much weaker; if 5303 Å was used, the Fig. 1 caption is misleading. The authors must specify which rest wavelength was adopted, correct the caption if needed, and recompute or re-derive the central velocity with the correct value.","section":"§3.1, Fig. 1 caption, Abstract"},{"comment":"The identification of the flattening in the north-eastern rim as a companion shadow is based on visual comparison with a single D6 forward model. No quantitative null test is presented: the authors do not compare the observed morphology with a spherical shell plus noise, with an ellipsoidal or otherwise asymmetric explosion model, or with ambient-medium density variations. Because the tracer points for the cone fit are manually selected on this feature, all subsequent parameters (cone angles in §3.2, and the secondary mass, separation, and orbital velocity in §3.3 and Table 1) inherit this assumption. The authors should provide a quantitative measure of the flattening significance, for example by fitting a circle or ellipse to the outer isophotes and testing the residuals, and by demonstrating that the D6 shadow model is favored over a symmetric model using a model-comparison statistic.","section":"§3.2, Fig. 4"},{"comment":"The conversion from the fitted cone half-opening angle α = 36.7° to the secondary mass via Eggleton's formula assumes that the secondary filled its Roche lobe at the moment of explosion and that the cone opening angle directly equals the Roche-lobe angular radius as seen from the primary. This geometric/physical mapping is not derived or tested; a non-Roche-filling companion or a different shadow geometry would change q, the secondary mass, and the orbital velocity. The authors should justify this mapping or explore how the inferred parameters depend on it. In addition, the fixed positional uncertainty σ = 0.5 pc in the cone fit is ad hoc and likely does not capture the systematic uncertainty from tracer-point selection; the quoted posterior uncertainties are therefore probably optimistic.","section":"§3.3, Table 1"},{"comment":"The quoted uncertainty of ±60 km/s for the bulk Doppler velocity is computed from the standard error of the median and the scatter among three maps, but it omits the uncertainty in the MUSE wavelength calibration, the rest-wavelength error discussed above, and the uncertainty in the LMC systemic velocity (293 km/s, van der Marel et al. 2002). A fair error budget should include these terms, since the bulk velocity is the quantitative link between the observation and the orbital-motion interpretation.","section":"§3.1"}],"minor_comments":[{"comment":"The left panel is labeled 'Redshifted [Fe XIV] ejecta' but the color bar is -7000 to 0 km/s, while the right panel is labeled 'Blueshifted' with a color bar 0 to 4000 km/s; the sign convention is opposite to standard usage and should be clarified or corrected.","section":"Fig. 1"},{"comment":"The text refers to a 'bulk systematic red-shift velocity of ~1000 km/s' while the Doppler velocities and histograms are negative; the authors should define the sign convention for redshift and blueshift consistently throughout.","section":"§3.1 and Fig. S4"},{"comment":"The phrase 'after correcting for the systemic velocity of the LMC' is ambiguous about whether the LMC velocity was added or subtracted; please state the exact transformation used.","section":"§3.1"},{"comment":"There are typographical issues: 'matploblib' should be 'matplotlib', and several author names contain encoding artifacts such as 'R¨opke' and 'M¨unchen'.","section":"Acknowledgments"},{"comment":"The description of the neural network training omits the line-width, continuum, and noise model used to generate the synthetic spectra; please provide these details for reproducibility.","section":"§2.2"},{"comment":"Several references appear to be in-press items with unformatted article numbers (e.g., Das et al. 2026 'stag596', Pollin et al. 2026 'stag735'); please update these to the final journal identifiers.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The rest-wavelength inconsistency is a serious but potentially fixable issue; the authors should clarify whether 5303 Å or 5308 Å was used and propagate the corresponding systematic uncertainty. The lack of a null test for the flattened edge is the other load-bearing weakness; without it, the companion-shadow interpretation and the derived binary parameters are not yet established. I would also ask the editor to verify that the numerous 2026 references, several by the authors themselves, are all published or accepted, since some appear as in-press placeholders. The paper is potentially important, but it needs a careful revision before it can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper has a genuinely new observational result: the first global Doppler-split coronal-line map of an entire Type Ia supernova remnant in the optical, built from deep MUSE data. The map itself, plus the public code and forward-model comparison, is worth a close look. The clean consistency between the measured bulk redshift and the orbital velocity inferred from the cone geometry is real and encouraging — those two numbers are not derived from each other, so the agreement is the strongest point in the paper.\n\nThe soft spots are not minor. The rest-wavelength conflict is load-bearing. The abstract and text use 5303 Å, but the Fig. 1 caption says the map was made with 5308 Å. A 5 Å shift is ~280 km/s at this wavelength, several times the quoted 60 km/s uncertainty. If the wrong zero-point was used, the bulk velocity is partly an artifact; if the caption is wrong, the error budget still omits rest-wavelength and LMC systemic-velocity uncertainties. Either way, the -1000 km/s value is not securely calibrated, and the agreement with the cone-derived 990 km/s could be coincidental. The authors need to state the adopted rest wavelength in one place, justify it, and fold its uncertainty into the error budget.\n\nThe companion-shadow interpretation also needs more support. The flattened edge is identified by eye, with no null test against a spherical shell plus noise or against ambient density variations. The cone fit uses hand-selected tracer points on that feature, so all companion parameters inherit that selection. On top of that, the mapping from the cone half-opening angle (36.7 deg) to a Roche-lobe radius via Eggleton is not explained transparently; the implied angular radius of the secondary's Roche lobe for the derived mass ratio seems to be ~19 deg, about half the quoted opening angle. Unless there is a geometric factor I am missing, the companion mass and separation are not on solid footing.\n\nWhat holds up: the Doppler splitting across the remnant is likely real, the bulk motion is a plausible physical signal, and the overall double-degenerate interpretation is plausible and well contextualized. The authors are honest about the framework's limitations. But the current manuscript does not demonstrate the central claim at the level of \"reveal\".\n\nI would send this to a good referee with a request for major revision. The observational dataset and the core idea are important enough that a referee's time is well spent. The paper should provisionally be treated as a strong candidate, not a finished result.","headline":"A genuinely new Doppler-split [Fe XIV] map of SNR 0509, but the central progenitor claim rests on an uncalibrated rest wavelength and a visually identified shadow cone; deserves serious refereeing, not desk rejection.","tokens_in":14079,"tokens_out":2892,"would_cite":true,"duration_ms":31769,"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 paper claims that SNR 0509-67.5 exploded as a double-degenerate merger, with a ~0.6 solar-mass white dwarf companion casting a cone-shaped shadow on the ejecta and a bulk Doppler shift tracing the primary's orbital velocity.","keywords":["supernova remnants","Type Ia supernovae","white dwarf mergers","double-degenerate progenitors","integral field spectroscopy","coronal emission lines","Doppler kinematics","machine learning classification"],"falsifier":"A quantitative null test would settle this: generate synthetic [Fe XIV] shells with the same spaxel noise, binning, and background density variations but without a companion, and ask whether a north-eastern flat edge appears as often as observed; if it does, the shadow interpretation is not required by the data. A complementary test would map the rim's surface brightness or proper motion to check whether the supposed shadow region is genuinely depleted rather than merely viewed through dense gas.","tokens_in":12962,"feed_emoji":"🌌","tokens_out":5788,"duration_ms":57377,"temperature":0.7,"pith_summary":"Using deep integral-field spectroscopy of the young Type Ia supernova remnant SNR 0509-67.5, the paper maps Doppler-split coronal [Fe XIV] emission across the whole remnant and finds a bulk redshift of about 1000 km/s in the ejecta. It interprets this as the line-of-sight component of the exploding white dwarf's orbital velocity, a signature expected in a double-degenerate merger but not in single-degenerate or most other explosion channels. In the same velocity maps, the remnant's north-eastern rim is flattened; the paper reads this flat edge as the conical shadow cast by a companion white dwarf on the expanding ejecta and fits a cone to it with Bayesian inference. Using the cone geometry together with the Eggleton Roche-lobe relation, it infers a companion of roughly 0.6 solar masses with an orbital separation of about 0.042 solar radii and an orbital velocity consistent with the measured Doppler shift. If correct, this gives a self-consistent dynamical picture of the progenitor binary at the moment of explosion.","feed_headline":"Flat rim and 1000 km/s shift reveal a white-dwarf merger","feed_subtitle":"Deep spectroscopic maps of a Type Ia remnant yield a 0.6-solar-mass companion and a nearly edge-on orbit.","key_machinery":"The load-bearing objects are four linked pieces. First, the Doppler-split [Fe XIV] emission line profiles from MUSE integral-field data, classified and fitted per spaxel with a dense neural network followed by Bayesian Gaussian fitting. Second, the three-dimensional reconstruction of the ejecta shell under the assumption of homologous expansion, using cos(θ) = V_obs / V_max to convert line-of-sight velocity into depth. Third, the shadow cone: a conical underdensity carved in the ejecta by the companion star, parameterised by a half-opening angle and two orientation angles with the apex fixed at the dynamical centre, which converts a visual flat edge into a geometric fit. Fourth, the Eggleton Roche-lobe relation, which translates the cone's opening angle into the binary mass ratio and hence into companion mass, radius, and orbital speed.","core_discovery":"The central discovery is that SNR 0509-67.5's reverse-shocked ejecta, traced by the [Fe XIV] λ5303 coronal line, carries two independent kinematic imprints of a close double white dwarf binary at explosion. First, the bulk Doppler velocity of the whole shell, -1000 ± 60 km/s, matches the expected orbital velocity of a ~1.1 solar-mass primary if the system was a merger with a ~0.6 solar-mass companion. Second, the blue-shifted and red-shifted ejecta maps both show a flattened north-eastern edge; forward hydrodynamic models of dynamically driven double-detonation mergers produce exactly such a flat edge as the shadow of the companion on the ejecta. Fitting this edge as a cone anchored at the dynamical centre yields a half-opening angle of 36.7 degrees and an orientation tilted about 4 degrees from the plane of the sky, from which the companion mass, radius, separation, and orbital velocity follow. The inferred primary orbital velocity of about 990 km/s agrees with the measured bulk shift, supporting a nearly edge-on orbital plane.","pith_inferences":["A similar cone-shadow analysis could be applied to other young Type Ia remnants, such as Kepler, Tycho, or SNR 0519-69.0, with existing or future integral-field cubes; the method only needs a Doppler-split coronal line and a chosen dynamical centre.","The bulk-shift interpretation predicts that the geometric centre of the [Fe XIV] shell should be offset from the dynamical centre by a few tenths of a parsec along the line of sight, which is potentially measurable with future high-resolution X-ray or optical spectroscopy.","If the companion's shadow is real, the ejecta inside the cone should be less shocked and less luminous, so maps of [Fe XIV] surface brightness or X-ray emission should show a deficit coincident with the flat edge; testing this correlation would sharpen the claim beyond a visual comparison."],"forward_implications":["For this remnant, explosion scenarios that require no companion at the moment of detonation, such as core-degenerate mergers or head-on white-dwarf collisions, are disfavoured by the data.","The ~1000 km/s bulk Doppler shift means the entire [Fe XIV]-emitting shell is moving as a unit, so any kinematic model of SNR 0509-67.5 must include this systemic offset rather than assuming a spherical rest-frame expansion.","The inferred binary parameters imply an extremely close pair at explosion: a companion radius near 9800 km and an orbital separation near 0.042 solar radii, consistent with the companion filling its Roche lobe just before detonation.","Because the outer ejecta geometry would be nearly unchanged if the secondary also exploded, the remnant morphology alone cannot settle the secondary's fate; deeper searches for a surviving companion or for intermediate-mass-element enrichment in the inner ejecta are the next observational tests."],"supporting_citations":[{"why":"Supplies the deep MUSE observations, the discovery of the [Fe XIV] line splitting, and the data-reduction products that this paper re-analyses.","marker":"Das et al. 2025"},{"why":"Provides double-degenerate merger models in which the primary's orbital velocity is imprinted on the inner ejecta, the basis for reading the -1000 km/s shift as orbital motion.","marker":"Pakmor et al. 2022"},{"why":"Shows three-dimensional hydrodynamic remnant models of the D6 double-detonation scenario retaining a flattened 'shadow of the companion' edge at ~500 years, the template for the observed north-eastern rim.","marker":"Ferrand et al. 2022"},{"why":"Extends the companion-shadow prediction to quadruple-detonation merger scenarios, supporting the generality of a flat edge in the projected ejecta.","marker":"Ferrand et al. 2025"},{"why":"Simulates remnants with no companion, an exploding companion, and a surviving companion, demonstrating that the companion-induced asymmetry persists into the remnant phase.","marker":"Prust et al. 2026"},{"why":"Provides the dynamical centre coordinates used as the fixed apex of the fitted shadow cone.","marker":"Arunachalam et al. 2022"},{"why":"Supplies the Roche-lobe relation that converts the fitted cone half-opening angle into the binary mass ratio.","marker":"Eggleton 1983"},{"why":"Supplies the adopted primary white dwarf mass of ~1.1 solar masses for SNR 0509-67.5, used to scale the companion properties.","marker":"Mandal et al. 2026"}],"fun_headline_variants":["Companion shadow and Doppler shift expose white-dwarf merger","MUSE data reveal double-degenerate merger in Type Ia remnant","Flat ejecta edge and 1000 km/s shift pin down merger","Near edge-on orbit inferred from companion's shadow in SNR"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on the claim that the flattened north-eastern rim is the companion's shadow rather than an asymmetric explosion, an uneven ambient medium, or noise; if that identification is wrong, the cone fit and all inferred companion properties do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Companion shadow and Doppler shift expose white-dwarf merger","MUSE data reveal double-degenerate merger in Type Ia remnant","Flat ejecta edge and 1000 km/s shift pin down merger","Near edge-on orbit inferred from companion's shadow in SNR"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000377,"raw_usage":{"total_tokens":2061,"prompt_tokens":1050,"completion_tokens":1011,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":666,"completion_tokens_details":{"reasoning_tokens":939}},"tokens_in":666,"tokens_out":1011,"duration_ms":11280,"temperature":1.0,"reasoning_tokens":939,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:21:14.022133+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A quantitative null test would settle this: generate synthetic [Fe XIV] shells with the same spaxel noise, binning, and background density variations but without a companion, and ask whether a north-eastern flat edge appears as often as observed; if it does, the shadow interpretation is not required by the data. A complementary test would map the rim's surface brightness or proper motion to check whether the supposed shadow region is genuinely depleted rather than merely viewed through dense gas.","supporting_citations":[],"review_version":1}