{"id":"d393a8a9-f133-455b-af0b-40700b48c586","arxiv_id":"2412.18226","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Companion-interaction wakes in Type Ia supernova remnants drive off-center reverse-shock convergence, pull interstellar medium to the remnant center, and keep the remnant X-ray asymmetric for thousands of years.","lead":"This paper simulates what happens when the blast wave of a Type Ia supernova slams into its companion star, carving a low-density wake in the ejected gas. The wake survives into the remnant phase, skewing the reverse shock, pulling interstellar gas to the center, and keeping the remnant lopsided in X-rays for thousands of years.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Axial density floor and wake extrapolation to r=0 are the load-bearing input for the ISM-to-center claim; a sensitivity test is needed.","rationale":"The reader correctly identified the ad hoc density floor as a weakness, and I agree that it is the least secure input. My concern sharpens it: the fitting procedure not only chooses an arbitrary floor value, it also extends the low-density wake to r = 0, where the physical wake should not exist because the companion is located at a = 0.269R_sun and the wake is formed downstream of that radius. The central region should instead contain unperturbed or bow-shocked ejecta at much higher density. Because the paper's headline result—ISM reaching the exact center—depends on the reverse shock propagating fastest along this axis, both the floor value and the radial extent of the wake are load-bearing. The paper gives the Sprout and Athena++ code versions, a 1D/semi-analytic validation in Section 5.4, and an η = 3 spherical-limit test in Section 6, all of which support the broader framework. However, none of these validate the central axial density structure, and there is no resolution study for the 268M-cell run. The proposed test directly isolates the sensitivity of the center-reaching claim. If it passes, the conclusion is substantially strengthened; if it fails, the long-lived X-ray asymmetry from ISM reaching the center would need to be revised to a near-center or off-center effect. Therefore the reader's CONDITIONAL verdict is appropriate, with the added condition of a wake-density sensitivity study.","tokens_in":17114,"tokens_out":11637,"duration_ms":110780,"concrete_test":"Re-run the Sprout DD remnant evolution with the same 268M-cell setup but replace the Appendix A axial density prescription with two variants: (i) floor density 0.1ρG(r) instead of 0.01ρG(r); (ii) a wake truncated at the homologously-expanded donor position r_wake(t0) = 0.269R_sun, i.e., unperturbed Gaussian density for r < r_wake at the initial time and the low-density cone only for r > r_wake. In each run, record the time and radius of reverse-shock convergence and the earliest time at which the ISM passive tracer reaches a sphere of radius 0.1 pc. If both variants still show ISM reaching r < 0.1 pc by 3000 yr and the convergence time changes by less than roughly 30%, the center-reaching claim is robust; if either fails, the claim is an artifact of the ad hoc axial density structure.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim of Section 5.3—that the reverse shock traverses the wake, converges off-center, and draws ISM all the way to the remnant center—rests on the density structure imposed along the symmetry axis in the Sprout initial condition. The Athena++ data are not mapped directly into Sprout; instead, the Appendix A fits are used and 'we also enforce ρ(r, θ) ≥ 0.01ρG(r) so that the density does not drop too low at θ = 0.' This floor is not derived from the collision simulation, and no sensitivity study is presented. The issue is more than the floor value: the same fitting procedure extrapolates the low-density cone to r = 0, although physically the wake is carved downstream of the donor at a = 0.269R_sun and should not exist interior to the homologously-expanded donor radius (≈0.002 pc at t = 10 yr, ≈0.2 pc at t ≈ 1000 yr). The unperturbed ejecta in that central region is roughly 100 times denser than the imposed floor at late times. Since the reverse-shock propagation speed and the advection of ISM into the center are controlled precisely by this axial density contrast, either an increased floor or a physically motivated truncation of the wake at the donor radius could delay or prevent the claimed center-reaching of ISM material. This is a concrete, testable vulnerability in the central claim, not a numerical detail.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the late-time evolution of Type Ia supernova remnants that contain a low-density wake carved by the interaction of the supernova ejecta with a surviving binary companion. The authors first run Athena++ simulations of the ejecta-companion collision for a double-degenerate (double-detonation) model, construct analytic fits to the resulting density structure, and compare with the Kasen (2010) single-degenerate fit. These structures are mapped into the expanding-grid code Sprout and evolved from 10 yr to 3000 yr. The reported results are that the forward shock becomes spherical by roughly 1000 yr owing to transverse flow, the reverse shock crosses the low-density wake, converges off-center, and draws interstellar medium into the remnant center, and that Rayleigh-Taylor plumes at the wake edge create a toroidal ejecta structure with observable thermal X-ray asymmetries persisting for thousands of years.","tokens_in":17406,"tokens_out":6826,"duration_ms":64284,"significance":"If the central claims hold, the paper provides a concrete observational discriminant between remnants of Type Ia supernovae with a surviving companion and remnants from mergers, at ages of order 10^3 yr, tied to off-center reverse-shock convergence and deep ISM ingress. The work is a forward simulation rather than a fit to observed remnants, uses standard publicly available codes, and is transparent about its free parameters and numerical choices; the code forks are made available for reproduction. The main value is the falsifiable prediction of a long-lived X-ray asymmetry and its viewing-angle dependence. The strength of the claim is currently limited by the sensitivity of the initial wake structure to unverified assumptions, as detailed below.","major_comments":[{"comment":"The ad hoc density floor and the extrapolation of the fitted low-density cone to r = 0 are load-bearing for the central claim of §5.3 that the reverse shock traverses the wake and draws ISM to the remnant center. Physically, the wake is carved downstream of the donor at a = 0.269 R_sun and should not extend interior to the homologously expanded donor radius (about 0.002 pc at 10 yr and 0.2 pc at 1000 yr); in that interior region, unperturbed ejecta is roughly 100 times denser than the imposed floor. Because the reverse-shock speed and the advection of ISM into the center are controlled by this axial density contrast, the center-reaching behavior is not robust until at least one sensitivity test is provided: for example, a higher floor, a wake truncated at the donor radius, or a direct mapping of the Athena++ profile instead of the fitted form.","section":"Appendix A, §5.3"},{"comment":"No resolution study is presented for either code. The Athena++ collision run uses only 10 azimuthal cells, and the Sprout remnant runs use a single 512x512x1024 grid. The claims of large Rayleigh-Taylor plumes around the wake edge and a toroidal ejecta structure (§5.3 and Figs. 6, 7, 11) are morphological results that could depend on resolution and on the suppression of non-axisymmetric modes. Please add a resolution test (for example, at least one coarser and one finer Sprout run, and ideally a check of the Athena++ azimuthal resolution) or otherwise bound the numerical uncertainty of these features.","section":"§2, §4"}],"minor_comments":[{"comment":"The title contains a stray space: 'Ejecta W akes' should be 'Ejecta Wakes'.","section":"Title"},{"comment":"The sentence defining 'unperturbed' appears to have reversed wording: it says the term refers to the portion of the ejecta altered by the companion interaction, but the subsequent usage refers to ejecta unaffected by the companion.","section":"§4"},{"comment":"The text contains a typo, 'morpohology', which should be 'morphology'.","section":"§5.3"},{"comment":"Equation (6) is followed by an empty equation number (7); the numbering should be checked.","section":"§2, Eq. (6)"},{"comment":"The statement that the semi-analytic model gives the correct reverse-shock convergence time 'only due to our choice of aRS' is an important limitation and should be flagged more prominently in the discussion, since Figure 10 could otherwise be read as independent validation.","section":"§5.4"}],"recommendation":"major_revision","confidential_remarks":"No concerns about citation practice or novelty. The main risk is the density-floor sensitivity of the initial conditions; the requested sensitivity runs appear feasible and would substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you the short version: this is a good simulation paper, worth reading carefully, but the flashiest claim—ISM drawn all the way to the center through the companion wake—rides on an assumed axial density floor and an extrapolation of the wake to r=0 that the paper doesn't test.\n\nWhat's actually new: Prust et al. take a 3D Athena++ collision simulation for a double-degenerate Ia, fit the wake structure, then evolve both DD and SD ejecta to 3000 yr on the expanding grid Sprout. The late-time behaviors—off-center reverse-shock convergence, toroidal vortex of ejecta, forward shock sphericalizing around 1000 yr—are genuine extensions of earlier SD simulations. They also compare with 1D hydro and a Truelove-McKee semi-analytic model, and they're transparent that the latter's agreement on convergence time hinges on the chosen reverse-shock acceleration. The X-ray emission measure maps give a concrete observable for XRISM-era work.\n\nCredit where due: the paper is clearly written, the codes and forks are public, and the authors distinguish what they derive from what they choose. That honesty matters.\n\nThe soft spot is concentrated in Section 5.3. The claim that ISM reaches the remnant center depends on the axial density in the wake being very low all the way in. But the fitted wake profile enforces ρ ≥ 0.01 ρ_G at θ=0 (Appendix A), and the same fit extends the low-density cone to r=0. Physically, the wake is carved by the donor at orbital separation a=0.269 R_sun; interior to the homologously expanded donor radius there shouldn't be an empty cone. The unperturbed ejecta there is roughly 100× denser at late times. So the reverse shock propagation and ISM advection are set by an imposed, not derived, density contrast. The stress-test note is right that this is a concrete testable vulnerability, not a numerical quibble.\n\nMinor issues: no resolution study anywhere (Athena++ runs with 10 azimuthal cells; Sprout at 512×512×1024, no refinement test), and the SD comparison uses the Kasen filling factor fh=0.1 without sensitivity. These are smaller because the global morphology likely doesn't hinge on them, but they should be addressed.\n\nOverall: the paper deserves a serious referee. The central dynamical picture—wake stays asymmetric for thousands of years, reverse shock converges off-center—is plausible and well supported. But the ISM-to-center result should be conditioned on a sensitivity study of the density floor and a wake truncated at the donor radius. If that holds, this becomes an important paper for the subfield.\n\nRecommendation: send it to peer review, but ask for the sensitivity test in revision.","headline":"Solid 3D study of companion wakes in Ia remnants, but the flashiest claim—ISM drawn to the center—rests on an untested density floor and a questionable wake extrapolation to r=0.","tokens_in":17951,"tokens_out":2004,"would_cite":true,"duration_ms":19230,"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 reverse shock races down the low-density wake carved by a surviving companion, converges off-center, and pulls interstellar gas into the remnant's core, leaving an X-ray asymmetry that lasts thousands of years.","keywords":["Type Ia supernovae","supernova remnants","companion interaction","reverse shock","interstellar medium entrainment","X-ray emission","hydrodynamic simulation","white dwarf binaries"],"falsifier":"High-resolution X-ray spectroscopy of the central few parsecs of a young type Ia remnant, such as SNR 0509-67.5, should detect ISM line emission if the claim holds; if no interstellar material is found at the center and the reverse-shock ring is centered on the explosion point instead of offset, the claim would be ruled out.","tokens_in":16888,"feed_emoji":"💫","tokens_out":9044,"duration_ms":78973,"temperature":0.7,"pith_summary":"This paper tries to show that the conical wake a surviving companion carves out of type Ia supernova ejecta controls the remnant's late-time appearance. In hydrodynamical simulations, the wake is underdense, so the reverse shock crosses it early, converges off-center, and draws interstellar medium all the way to the remnant's center. That inward channel plus a ring vortex around the wake keeps the remnant's X-ray emission measurably asymmetric for thousands of years, even after the forward shock itself has become spherical. If correct, this gives an observational route to identifying which type Ia remnants had a surviving companion rather than a merger.","feed_headline":"Companion wake sends interstellar gas to a supernova remnant's center","feed_subtitle":"Simulations show the reverse shock dives down the wake, twists the remnant's core, and keeps the X-ray asymmetry visible for millennia.","key_machinery":"The central object is the companion-interaction wake: the cone of low-density, higher-velocity ejecta left behind after the supernova blast collides with the surviving donor star. Its underdensity is the mechanism that lets the reverse shock travel quickly into the remnant, converge off-center, and entrain interstellar gas toward the center. A second mechanism, transverse flow of shocked ejecta driven by pressure gradients, refills the wake on roughly a sound-crossing time and restores sphericity to the forward shock without erasing the asymmetry in the core. The numerical machinery is a two-stage calculation: Athena++ for the ejecta-companion collision, followed by the uniformly expanding-grid code Sprout for the 3-D remnant evolution to 3000 years.","core_discovery":"The paper's central discovery is that the companion-interaction wake, not the forward shock, determines the long-lived asymmetry of a type Ia supernova remnant. In both the double-degenerate model simulated with Athena++ and the single-degenerate model built from the Kasen (2010) wake fit, the wake is a cone of low-density, shock-heated, slightly accelerated ejecta. Because it is underdense, the reverse shock traverses the wake sooner than elsewhere, converges off-center (colliding with the reverse shock from the other side roughly 2 pc from the explosion), rebounds, and drives high-velocity gas back up the wake. A passive tracer shows interstellar gas is pulled through the wake and reaches the center of the remnant in both models. The forward shock becomes spherical by about a thousand years because transverse pressure gradients push shocked ejecta into the wake, but the reverse shock, the ejecta distribution, and the X-ray emission remain asymmetric.","pith_inferences":["If interstellar gas really reaches the center through the wake, then spatially resolved X-ray line-ratio maps of the cores of nearby type Ia remnants could directly test the surviving-companion scenario, because a merger remnant should lack this inward ISM channel.","The finding that the reverse shock converges off-center predicts a central hot region offset from the explosion point; high-resolution Doppler mapping of X-ray lines could locate the companion's original direction.","A natural numerical extension is to repeat the remnant evolution with the density floor in the wake fit raised or lowered, since the ISM-to-center result depends on how empty the wake is allowed to become.","Analogous wake dynamics may appear in other transients with a surviving companion, although the different ejecta structure would change the timing of reverse-shock convergence."],"forward_implications":["Thermal X-ray maps of type Ia remnants can remain measurably asymmetric for more than a thousand years, so the absence of an early spherical remnant does not by itself require a merger.","Interstellar gas can be transported to the center of a young type Ia remnant through the wake, so central X-ray spectra may show ISM-like abundances mixed with ejecta.","A forward shock that is spherical at late times does not imply the underlying ejecta are spherically distributed; the contact discontinuity and reverse shock stay off-center.","The X-ray asymmetry is strongest when the remnant is viewed perpendicular to the wake, but it is still visible from other angles, so viewing geometry must be folded into any comparison with observations.","Wider binaries with smaller donor solid angles produce increasingly spherical X-ray morphologies, meaning a symmetric remnant does not rule out a surviving companion."],"supporting_citations":[{"why":"Supplies the Gaussian ejecta density profile, donor radius and orbital separation, and companion-response modeling used to set up the Athena++ collision simulation.","marker":"Wong et al. (2024)"},{"why":"Provides the single-degenerate wake density fit that is used as the SD initial condition for the remnant evolution.","marker":"Kasen (2010)"},{"why":"Introduces Sprout, the uniformly expanding-grid code whose moving-mesh flux treatment is used to evolve the remnant from 10 to 3000 years.","marker":"Mandal & Duffell (2023)"},{"why":"Supplies the inferred ISM density for SNR 0509-67.5 adopted as the ambient medium in the remnant simulations.","marker":"Arunachalam et al. (2022)"},{"why":"Earlier axisymmetric SPH simulations whose wake-edge Rayleigh-Taylor plumes and early remnant morphology are reproduced and extended by the present 3-D runs.","marker":"García-Senz et al. (2012)"},{"why":"Provided the flat-edge X-ray asymmetry seen at 100 and 300 years that this paper's later-time evolution builds on.","marker":"Gray et al. (2016)"},{"why":"Earlier double-degenerate remnant simulation with a different code, used to compare protrusion and reverse-shock behavior in the wake.","marker":"Ferrand et al. (2022)"},{"why":"Semi-analytic forward/reverse shock trajectory framework that the paper extends to a Gaussian ejecta profile to validate the hydrodynamic results.","marker":"Truelove & McKee (1999)"}],"fun_headline_variants":["Companion's wake funnels interstellar gas to remnant's center","Reverse shock dives wake, twists supernova remnant core","Wake-shaped remnant stays asymmetric for thousands of years","How a companion's wake leaves a lasting X-ray mark"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the fitted wake density never falls below 1% of the unperturbed ejecta density in the double-degenerate model (and 10% in the single-degenerate model); if a real wake is denser than that floor, the reverse shock would slow down and interstellar gas might not reach the remnant's center.","fun_headline_variants_meta":{"raw":{"variants":["Companion's wake funnels interstellar gas to remnant's center","Reverse shock dives wake, twists supernova remnant core","Wake-shaped remnant stays asymmetric for thousands of years","How a companion's wake leaves a lasting X-ray mark"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1553,"prompt_tokens":969,"completion_tokens":584,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":585,"completion_tokens_details":{"reasoning_tokens":517}},"tokens_in":585,"tokens_out":584,"duration_ms":5520,"temperature":1.0,"reasoning_tokens":517,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:54:34.967295+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"High-resolution X-ray spectroscopy of the central few parsecs of a young type Ia remnant, such as SNR 0509-67.5, should detect ISM line emission if the claim holds; if no interstellar material is found at the center and the reverse-shock ring is centered on the explosion point instead of offset, the claim would be ruled out.","supporting_citations":[{"cited_title":"Shocking and Mass Loss of Compact Donor Stars in Type Ia Supernovae","cited_arxiv_id":"2408.00125","evidence_quote":"Supplies the Gaussian ejecta density profile, donor radius and orbital separation, and companion-response modeling used to set up the Athena++ collision simulation."},{"cited_title":"P., Hovey , L., & Eriksen , K","cited_arxiv_id":null,"evidence_quote":"Supplies the inferred ISM density for SNR 0509-67.5 adopted as the ambient medium in the remnant simulations."},{"cited_title":"2012, , 745, 75, 10.1088/0004-637X/745/1/75","cited_arxiv_id":null,"evidence_quote":"Earlier axisymmetric SPH simulations whose wake-edge Rayleigh-Taylor plumes and early remnant morphology are reproduced and extended by the present 3-D runs."},{"cited_title":"J., Raskin , C., & Owen , J","cited_arxiv_id":null,"evidence_quote":"Provided the flat-edge X-ray asymmetry seen at 100 and 300 years that this paper's later-time evolution builds on."}],"review_version":1}