{"id":"77164ef1-ec9a-4a3f-b192-549dc187cfed","arxiv_id":"2507.06412","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The wake carved by a surviving companion in a Type Ia supernova produces detectable, viewing-angle dependent distortions in the [Co III] and [Ar III] nebular line profiles.","lead":"Companion stars in Type Ia supernova progenitors leave a cone-shaped wake in the explosion debris, and this paper predicts how that wake bends the shapes of late-time infrared emission lines. The predicted asymmetries are strong and angle-dependent enough to be searched for in JWST spectra, giving a new probe of how these supernovae form.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Full-sphere extrapolation at 80.2° sets the entire unperturbed background for every predicted line profile and viewing-angle trend; if the omitted polar cap is not spherical, the claimed JWST-identifiable horns and first-moment slopes are not yet established.","rationale":"I read the paper as establishing a new observable diagnostic: companion wakes produce viewing-angle-dependent line-profile asymmetries in nebular [Co III] and [Ar III] lines, strong enough for JWST. The argument is mostly sound: the hydrodynamics is standard Athena++, the passive-scalar analysis documents the wake and bow shock, the line-shape code is validated against analytic Gaussian and thin-shell profiles, and the wake-free comparison to SN 2021aefx via the Gaussian model is independent supporting evidence. The weakest load-bearing step is not the local wake physics but the construction of the full sphere. Section 4 explicitly extrapolates the density at 80.2° to all larger polar angles. Every synthetic line profile and every first-moment trend is computed on this assumed sphere. The strongest claim is framed as a quantitative detectability statement, which requires the baseline to be right; the validation case does not test it because it contains no wake and no extrapolated cap. I therefore agree with the reader's identification of this assumption. A full-polar-angle simulation is the natural, cheap, decisive check: same resolution, same parameters, extend θ to 180°. If the horns and moment slopes survive, the extrapolation was benign; if they shift, the current quantitative predictions are conditional on the cap. This does not change the verdict: conditional acceptance is appropriate, and the proposed simulation should be a stated condition.","tokens_in":13989,"tokens_out":9305,"duration_ms":134317,"concrete_test":"Rerun the same Athena++ simulation with the polar domain extended to at least 120° (ideally 180°) at matched resolution, keeping the identical donor, injection, and physical parameters. Recompute the [Co III] 11.89 µm and [Ar III] 8.99 µm profiles at viewing angles 0°, 40°, and 90° from the extended density field, and compare them with the 80°-extrapolated predictions. If the horns, double-horned shapes, and the sign and slope of the first-moment curves survive, and the extrapolated baseline matches the extended run within the 100 km/s line-shape resolution, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is in Section 4: 'We span the full sphere by extrapolating our empirical density-radius relation at our largest simulated angle of approximately 80 degrees around the rest of the sphere.' The simulation's polar domain ends at 80.2°, so all line profiles in Figure 5 and all first-moment curves in Figure 6 are computed on a sphere whose unobserved polar cap is assumed to be identical to the density at 80°. This matters because the strongest claim is quantitative ('distinctive enough to be identified in JWST observations'), and that assessment is made from normalized profiles and moment slopes built on that extrapolated baseline. The independent validation against SN 2021aefx uses a wake-free Gaussian profile and so gives no constraint on the polar cap. If the real D6 ejecta beyond 80° carries any global asymmetry—off-center ignition, a large-scale density gradient, or orbital structure—the baseline line shape, the normalization, and the viewing-angle trends would all change. The local wake features are probably robust, but the paper's central quantitative claim has not been tested against that possibility. This is an addressable modeling choice, not an observed constraint.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper simulates the interaction of Type Ia supernova ejecta with a surviving Roche-lobe-filling companion in the D6 double-degenerate double-detonation scenario using Athena++, showing that the companion leaves a low-density wake and a bow-shock overdensity. The authors introduce a lightweight line-shape code based on Jerkstrand (2017) that computes optically thin nebular line profiles from arbitrary 3D density fields, using an emissivity that interpolates between j proportional to rho and j proportional to rho^2 according to ne/ncrit. They validate the wake-free limit against analytic profiles and the 1D non-LTE calculations of Blondin et al. (2023), and reproduce the [Co III] 11.89 micron line of SN 2021aefx with an electron-density-sampled Gaussian model. With the simulated wake density structure, they predict viewing-angle-dependent horns and double-horned profiles for [Co III] and [Ar III] and compute first-velocity-moment versus cos(theta) trends, concluding that these features are distinctive enough to be identified in JWST observations.","tokens_in":14283,"tokens_out":8083,"duration_ms":100143,"significance":"This is a timely and useful contribution. If the predictions hold, the wake signatures provide a new, physically motivated observational diagnostic for surviving-companion SNe Ia in the JWST era, complementing the hypervelocity-WD evidence for the D6 channel. The line-shape code is simple, fast, and validated against analytic and 1D non-LTE benchmarks; the comparison to SN 2021aefx is a genuine non-calibrated test. The treatment of the critical-density transition with ne*rho sampling is a nice improvement over pure rho or rho^2 scalings. The main caveats are the unconstrained polar-cap extrapolation and the absence of a quantitative detectability threshold.","major_comments":[{"comment":"The full-sphere line profiles and moment trends rely on the assumption stated in Section 4: 'We span the full sphere by extrapolating our empirical density-radius relation at our largest simulated angle of approximately 80 degrees around the rest of the sphere.' Since the simulation domain extends only to a maximum polar angle of 80.2 degrees, every profile in Figure 5 and every first-moment curve in Figure 6 is computed on a sphere whose unobserved polar cap is assumed to be identical to the density at 80.2 degrees. The wake-free validation in Section 3.4 uses a spherically symmetric Gaussian model and therefore provides no constraint on this polar cap. I request an explicit statement of the empirical density-radius relation and a robustness test in which the polar cap is filled with at least one alternative prescription (for example, the undisturbed Gaussian density profile, or the 80-degree profile with a mild large-scale density gradient). The test should show whether the horn/double-horn features and the slope of the first-moment curve survive plausible variations of the cap; without such a test, the Section 5 claim that the features are 'distinctive enough to be identified in JWST observations' is not yet established.","section":"Section 4, first paragraph; Section 2.1"},{"comment":"The central claim of JWST detectability is stated without a quantitative detection criterion. The profiles in Figure 5 are normalized and shown without a noise model, spectral resolution, or exposure-time/sensitivity estimate, and the text does not specify how large the wake-induced deviations must be relative to the observational uncertainty to be identified. I request either a more cautious wording of the detectability claim or a simple estimate of the required S/N and spectral resolution, for example by convolving the synthetic profiles with a JWST-like line-spread function and adding representative noise. This would also help separate the intrinsic wake signature from the unconstrained polar-cap contribution discussed in the previous comment.","section":"Section 5 and Figure 5"},{"comment":"The emissivity model in Eq. (8) assumes a density-only dependence with a uniform isothermal, constant-ionization conversion ne proportional to rho, following the spherically averaged D6 model of Blondin et al. (2023). However, Section 2.2 and Figure 1 show that the wake is shock-heated and radiation-pressure dominated, with thermodynamic properties that differ from the ambient ejecta. If the temperature or ionization state in the wake differs from the surrounding ejecta at nebular epochs, the density contrasts may not translate directly into emissivity contrasts, and the predicted horns could be suppressed or enhanced. I request a sensitivity test that assigns to the wake region a plausible temperature or ionization offset (for example, a factor of two change in the effective emissivity normalization) and recomputes the [Co III] and [Ar III] profiles, or, alternatively, a discussion of why the uniform-ionization assumption remains valid inside the wake.","section":"Sections 2.2, 3.1, 3.2"}],"minor_comments":[{"comment":"The notation in Eq. (7) is ambiguous: the integral Vmax to V(nu) is not clearly connected to the discrete sum over j, and the quantities dA and V(nu) are not fully defined. Please clarify the mapping between velocity bins and planar slices and state whether the profiles are normalized before comparison.","section":"Section 3.3, Eq. (7)"},{"comment":"The observed peak flux ratio of approximately 7.3 is quoted without uncertainties or a description of the continuum subtraction and epoch choice, and it is compared with a predicted ratio of approximately 6.6; the discussion would benefit from error bars on the observed ratio and a note on how the result depends on the assumed density scaling.","section":"Section 3.4"},{"comment":"The statement that the first-moment behavior shows a 'strong dipolar asymmetry' should be qualified because the sign of the slope flips for rho^2 sampling and because the calculation inherits the polar-cap extrapolation uncertainty; including error bars or a band showing the sensitivity to the cap prescription would strengthen the figure.","section":"Section 4.3 and Figure 6"},{"comment":"The abstract says the paper presents a tool to quickly calculate line shapes, but no code repository or availability statement is provided. For reproducibility, please include a link to the code or a statement of availability.","section":"Abstract and Section 3.3"},{"comment":"The low-velocity extrapolation test in Figure 7 is shown only for the [Co III] line at 270 days with the ne*rho sampling scheme; a brief statement on whether the conclusion holds for the other sampling schemes and for [Ar III] would be useful.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid theoretical contribution, and the requested robustness tests are feasible within the scope of the manuscript. The main concern is the unconstrained polar-cap extrapolation, which affects the paper's central quantitative claim; this is addressable and does not require new observations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a genuinely useful paper. Its core contribution is explicit: it takes the low-density wake and bow-shock overdensity that companion interaction leaves in a D6 double-detonation ejecta and computes what those do to mid-IR nebular line profiles as a function of viewing angle. Prior work had noted wakes or studied ejecta structure, but not turned them into line-shape predictions for the specific JWST lines. The predicted horns, double-horned profiles, and first-moment shifts are new and testable.\n\nThe paper earns credit for being careful and honest. The emissivity derivation is standard physics, the code is validated against analytic Gaussian and thin-shell profiles and against the 1D non-LTE output of Blondin et al. (2023), and the comparison to SN 2021aefx gives a real sanity check. The authors also flag their own limitations in the appendix: the low-velocity interior is not modeled, and they show the effect of alternative inner-density choices. That is the right way to handle a known boundary of the model.\n\nThe soft spots are real but proportionate. The full-sphere extrapolation at 80.2 degrees is an assumption, not a result, and the stress-test note is correct that the first-moment slopes and the normalized baseline profiles inherit that uncertainty. I do not think it sinks the qualitative claim: the horns and the double-horned shape are local features of the wake region, which the simulation does resolve, so those should survive changes in the polar cap. What changes is the quantitative strength of the features and the slope of the moment-versus-angle curves. The paper is upfront about the extrapolation, but a referee should push for a short sensitivity test (e.g., symmetric vs. asymmetric polar-cap density) before publication. The hard-sphere donor is a separate approximation, but the authors justify it with prior work and discuss stripping in the appendix; I would not treat that as a fatal flaw. The code itself is not public, which is a minor disappointment given the paper advertises it as a tool, but the method is described well enough to reproduce.\n\nWho is this for? People working on Type Ia progenitor constraints with nebular spectroscopy, especially mid-IR JWST follow-up. They will get a concrete set of predictions and a fast way to test other density fields. This deserves a serious referee and, assuming the sensitivity test is added or the caveat is framed more carefully, should be published. I would take it to a reading group only if we were doing a session on SN Ia observables; otherwise it is a cite-and-move-on paper.","headline":"A solid, honest paper that converts companion-wake density perturbations into concrete JWST-era predictions for [Co III] and [Ar III] line shapes; the main caveat is the unseen polar cap, which mostly affects quantitative moment slopes rather than the qualitative horns.","tokens_in":14818,"tokens_out":2647,"would_cite":true,"duration_ms":31601,"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":"A surviving companion's wake shapes Type Ia supernova spectra into horns and notches that vary with viewing angle.","keywords":["Type Ia supernovae","white dwarf stars","nebular phase spectroscopy","companion interaction","ejecta wake","forbidden emission lines","double-detonation progenitors","JWST mid-infrared spectroscopy"],"falsifier":"Search a sample of a dozen nebular-phase JWST spectra of Type Ia supernovae at $t > 200$ days for the [Ar III] 8.99 micron line: if none of them shows the predicted single-horned or double-horned morphology, and the [Co III] 11.89 micron first velocity moments show no correlation with line morphology, the companion-wake signatures would fail to appear. A more direct test would measure the predicted linear slope of the first velocity moment versus $\\cos \\theta$ across randomly oriented objects; a flat distribution with no such trend would rule out the wake as the source of the predicted asymmetries.","tokens_in":13801,"feed_emoji":"🔭","tokens_out":11046,"duration_ms":102439,"temperature":0.7,"pith_summary":"This paper argues that the surviving companion star in a double-degenerate Type Ia supernova leaves a permanent scar on the exploding ejecta: a low-density wake flanked by a bow-shock overdensity. The authors simulate this ejecta-companion collision and then compute nebular-phase forbidden line profiles for the innermost iron-group material and for a shell of intermediate-mass elements. They find that the wake removes high-velocity material along the companion axis and piles up emitting material at the shock edge, producing strong peaks, single horns, and double-horned profiles whose shape changes with the viewing angle. If correct, these predicted morphologies give JWST a direct way to identify which Type Ia supernovae had a Roche-lobe-filling companion and to constrain the binary orientation at the time of explosion.","feed_headline":"Companion wakes leave horn-shaped marks in supernova spectra","feed_subtitle":"JWST late-time spectra can test double-degenerate Type Ia progenitors by reading line shapes.","key_machinery":"The load-bearing mechanism is the shock-heated low-density wake and its flanking bow-shock overdensity left behind when the ejecta runs into the Roche-lobe-filling donor. On top of that density field the paper builds a fast line-shape code that slices the homologous ejecta into planes perpendicular to the line of sight, assigns each plane a projected velocity, and integrates a density-based emissivity. The emissivity is sampled in three regimes: linearly in density when the electron density is far below the transition's critical density $n_{\\rm crit}$, quadratically when far above, and through the full expression $j \\propto \\rho\\, n_e/(1+n_e/n_{\\rm crit})$ when the ejecta density is near threshold, which is the case for [Co III] at 270 days. This critical-density sampling is what makes the predicted horns and notches angle-dependent rather than a simple rescaling of the spherically symmetric profile.","core_discovery":"The central discovery is that companion interaction permanently modifies the density, velocity, and composition structure of Type Ia ejecta in a way that survives to homologous expansion and shapes forbidden emission lines. In the simulation, the fastest-moving ejecta is trapped behind the donor and ends up as the slowest material inside the wake, inverting the usual velocity-composition ordering, while the bow shock creates an overdensity at roughly 40 degrees from the symmetry axis. When line emission is computed with the electron density tracked against the critical density of the transition, the [Co III] 11.89 micron line (tracing 56Ni in the innermost ejecta) loses high-velocity flux along the wake axis and gains a central bump viewed perpendicular to it, and the shell-like [Ar III] 8.99 micron line develops single-horned profiles along the axis and double-horned profiles perpendicular to it. The authors conclude that these angle-dependent asymmetries are distinctive enough to be identified in JWST nebular-phase spectra.","pith_inferences":["If the predicted horned profiles are found and are correlated with hypervelocity runaway white dwarfs, the line-shape test would connect individual supernovae to the double-detonation channel rather than only to the progenitor population as a whole.","The angle-dependent first-velocity-moment slope offers a statistical way to measure the orbital axis orientation of the progenitor binary from a sample of nebular spectra, without resolving the system.","The same plane-slicing tool could be applied to other 3D asymmetries, such as off-center ignition, large-scale 56Ni clumping, or circumstellar interaction, to separate their line-shape fingerprints from the companion wake.","Extending the simulation beyond 80 degrees with a self-consistent global ejecta model would test whether the predicted full-sphere line profiles survive; until that is done, the horned morphologies should be treated as a robust but not final prediction."],"forward_implications":["Late-time JWST spectra of Type Ia supernovae should show wake signatures: an axis-on view of [Co III] 11.89 microns missing high-velocity emission, and an edge-on view of [Ar III] 8.99 microns showing a double-horned shape.","The first velocity moment of the wake-perturbed lines follows a linear trend with $\\cos\\theta$, so a measured pattern of blueshifts and redshifts can be used to infer the binary orientation along the line of sight.","Monitoring a forbidden line across the epoch when the electron density passes $n_{\\rm crit}$ should reveal a measurable change in line shape, because the emissivity shifts from roughly $\\rho^2$ to linear in $\\rho$.","The same wake physics extends to other progenitor channels with surviving companions, including thick helium-shell and quadruple-detonation scenarios, widening the applicability beyond the double-degenerate case."],"supporting_citations":[{"why":"Supplies the hydrodynamical simulation setup, equation of state, donor treatment, and original wake calculation that this paper reuses.","marker":"Prust et al. (2025)"},{"why":"Provides the Gaussian ejecta density model and its mass and energy parameters used for the injected ejecta and the undisturbed baseline line shapes.","marker":"Wong et al. (2024)"},{"why":"Provides the 3D double-detonation model whose spherically averaged density profile underlies the Blondin et al. radiative-transfer comparison.","marker":"Gronow et al. (2021)"},{"why":"Supplies the non-LTE treatment, isothermal temperature and ionization profiles, and the 1D double-detonation spectrum against which the line-shape code is tested.","marker":"Blondin et al. (2023)"},{"why":"Gives the planar-slice formalism and the analytical Gaussian and thin-shell line shapes on which the line-shape code is built and validated.","marker":"Jerkstrand (2017)"},{"why":"Establishes that [Co III] 11.89 microns is optically thin in the nebular phase and supplies observed line behavior used to set the critical density argument.","marker":"Gerardy et al. (2007)"},{"why":"Supplies the Einstein A coefficients and collision strengths used to compute the critical density for [Co III].","marker":"Storey & Sochi (2016)"},{"why":"Provides the JWST spectrum of SN 2021aefx used to validate the wake-free line shape and the peak flux ratio between epochs.","marker":"DerKacy et al. (2023)"},{"why":"Provides additional JWST nebular-phase spectra of Type Ia supernovae against which the predicted wake asymmetries can be searched for.","marker":"Kwok et al. (2023)"}],"fun_headline_variants":["Wake inversion shapes supernova spectral horns","Companion wake flips velocity order in Type Ia ejecta","JWST can spot companion wake in nebular line shapes","Thermonuclear supernova wake leaves spectral signature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole-sphere line profiles assume that the ejecta beyond a polar angle of about 80 degrees looks exactly like the ejecta at 80 degrees, because the simulation only extends to 80.2 degrees and the rest of the sphere is filled by extrapolating that density-radius relation; if the true outer ejecta geometry differs, every predicted full-sphere line shape and viewing-angle trend changes.","fun_headline_variants_meta":{"raw":{"variants":["Wake inversion shapes supernova spectral horns","Companion wake flips velocity order in Type Ia ejecta","JWST can spot companion wake in nebular line shapes","Thermonuclear supernova wake leaves spectral signature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000397,"raw_usage":{"total_tokens":2122,"prompt_tokens":1033,"completion_tokens":1089,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":649,"completion_tokens_details":{"reasoning_tokens":1026}},"tokens_in":649,"tokens_out":1089,"duration_ms":8674,"temperature":1.0,"reasoning_tokens":1026,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:05:36.006778+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Search a sample of a dozen nebular-phase JWST spectra of Type Ia supernovae at $t > 200$ days for the [Ar III] 8.99 micron line: if none of them shows the predicted single-horned or double-horned morphology, and the [Co III] 11.89 micron first velocity moments show no correlation with line morphology, the companion-wake signatures would fail to appear. A more direct test would measure the predicted linear slope of the first velocity moment versus $\\cos \\theta$ across randomly oriented objects; a flat distribution with no such trend would rule out the wake as the source of the predicted asymmetries.","supporting_citations":[{"cited_title":"J., Kumar, G., & Bildsten, L","cited_arxiv_id":null,"evidence_quote":"Supplies the hydrodynamical simulation setup, equation of state, donor treatment, and original wake calculation that this paper reuses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian ejecta density model and its mass and energy parameters used for the injected ejecta and the undisturbed baseline line shapes."},{"cited_title":"E., Sim, S","cited_arxiv_id":null,"evidence_quote":"Provides the 3D double-detonation model whose spherically averaged density profile underlies the Blondin et al. radiative-transfer comparison."},{"cited_title":"J., Ramsbottom, C","cited_arxiv_id":null,"evidence_quote":"Supplies the non-LTE treatment, isothermal temperature and ionization profiles, and the 1D double-detonation spectrum against which the line-shape code is tested."},{"cited_title":"L., Meikle, W., Kotak, R., et al","cited_arxiv_id":null,"evidence_quote":"Establishes that [Co III] 11.89 microns is optically thin in the nebular phase and supplies observed line behavior used to set the critical density argument."},{"cited_title":"2016, Monthly Notices of the Royal Astronomical Society, 459, 2558","cited_arxiv_id":null,"evidence_quote":"Supplies the Einstein A coefficients and collision strengths used to compute the critical density for [Co III]."},{"cited_title":"2023, The Astrophysical Journal Letters, 945, L2","cited_arxiv_id":null,"evidence_quote":"Provides the JWST spectrum of SN 2021aefx used to validate the wake-free line shape and the peak flux ratio between epochs."},{"cited_title":"A., Jha, S","cited_arxiv_id":null,"evidence_quote":"Provides additional JWST nebular-phase spectra of Type Ia supernovae against which the predicted wake asymmetries can be searched for."}],"review_version":1}