{"id":"2d59d332-a3b6-4fa9-9f74-8851bd8d3bbf","arxiv_id":"2608.06461","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A 3D simulation of nine repeated nova eruptions shows that nested shell collisions produce a slowly fading soft X-ray component and episodic hard X-ray flares, roughly matching the extended X-rays seen around RS Ophiuchi.","lead":"Repeated nova eruptions build nested shells that collide and emit X-rays for over a century. A three-dimensional simulation of a nine-eruption sequence reproduces the order of magnitude of RS Ophiuchi's diffuse X-ray glow.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NEI ionization-clock assignment via the pressure threshold in Appendix A is the least secure link; the paper never quantifies how Pth changes the band luminosities, so the RS Oph comparison may be threshold-dependent.","rationale":"The reader's weakest-assumption analysis correctly identifies the NEI clock. The central quantitative claim—the order-of-magnitude match of the fiducial soft-band luminosity to RS Oph's extended emission—is the bridge from simulation to observation, and it sits directly on the pressure-mask timescale assignment. The hydrodynamics itself has independent support: the nested-shell morphology is visible in the maps, the resolution comparison bounds the numerical error at about 0.14 dex (soft) and 0.33 dex (hard), and the internal light-curve behavior (soft secular decline, hard episodic flares) is physically plausible. What is not supported is the absolute normalization of the X-ray emissivities: vvpshock depends on net, and the only calibration of that parameter in this multi-eruption flow is an empirical pressure threshold validated against morphology, not against a tracer or observations. The paper's own caveat about mixed/repeatedly shocked cells and its failure to quantify luminosity sensitivity to Pth make the concern concrete rather than speculative. A purely post-processing sweep of Pth—or a passive-scalar tracer—would settle it without changing the hydrodynamics. If the luminosities are stable to about 0.5 dex, the RS Oph comparison stands; if they move by more than that, the conclusion should be rephrased as qualitative. The absence of code/data and the missing baseline run are secondary; they do not change the verdict but reinforce CONDITIONAL. No adjustment to the reader's CONDITIONAL verdict is needed.","tokens_in":11684,"tokens_out":9906,"duration_ms":88063,"concrete_test":"Recompute Fig. 5 from the stored snapshots with Pth swept over 1e-7, 3.16e-7, and 1e-6 dyn/cm2, and with a passive-scalar tracer of the most recent ejecta replacing the pressure mask (or, failing that, a rerun of one cycle with a passive scalar). Record the soft and hard luminosities at 3.43 and 5.28 yr after the t=108 yr eruption, the epochs of the Montez et al. (2022) points. If the soft luminosity varies by more than about 0.5 dex across the sweep, the RS Oph comparison is threshold-dominated; if stable, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix A assigns NEI ages by a single pressure cut: cells with P > 3.16e-7 dyn/cm2 are 'fresh' and get t_char = t - t_last; all others get one mean elapsed time over all prior eruptions. This is the pivot for the band luminosities in Fig. 5. vvpshock spectra are strongly nonlinear in net, so a wrong clock in a small but dense volume can shift soft/hard fluxes by factors. The calibration is only morphological: the paper says modest Pth variations do not change the identified post-shock layer, but it never reports the corresponding change in the 0.5-2.0 or 2.0-10.0 keV luminosities, let alone at the RS Oph epochs. Because the same threshold separates the 'young' hard-emitting gas from the 'old' soft-emitting gas, the claimed order-of-magnitude agreement with Montez et al. (2022) and the episodic hard flares could both be artifacts of the mask rather than of the hydrodynamics. The appendix's 'averaged time since all prior eruptions' for old cells is also unmotivated: most old gas was shocked at specific earlier times and then expanded, so a single mean age does not represent its ionization state. A tracer-based or swept-threshold test is needed before the quantitative comparison is accepted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents three-dimensional AMR hydrodynamical simulations of a nine-eruption recurrent-nova sequence spanning 130 years, using RS Oph-like parameters and a cadence derived from the documented RS Oph outburst record. The gas is evolved with CIE cooling plus on-the-fly dust cooling, and the density and temperature snapshots are post-processed with the XSPEC vvpshock NEI model to compute intrinsic luminosities in the 0.5-2.0 keV and 2.0-10.0 keV bands. The authors find that repeated eruptions excavate a bipolar cavity bounded by nested shell structures; soft X-rays trace dense shell rims and decline secularly as the remnant expands, while hard X-rays arise from the hottest shocked gas and appear as episodic flares following individual eruptions. The simulated soft-band luminosity is within an order of magnitude of the extended RS Oph luminosities reported by Montez et al. (2022). The central quantitative results rely on an empirically calibrated pressure threshold (Appendix A) that assigns NEI ionization timescales to freshly shocked versus older mixed gas.","tokens_in":12052,"tokens_out":3767,"duration_ms":33009,"significance":"If the results are robust, the paper establishes that diffuse shell-shock X-ray emission from nested nova remnants is a natural, observable long-lived component of recurrent novae, with implications for interpreting extended X-ray emission around symbiotic recurrent novae and for the CSM structure encountered by any eventual supernova explosion. The work moves beyond single-eruption modeling and produces concrete, falsifiable predictions of band-dependent variability on decade timescales. Strengths include the 3D AMR setup with physically motivated wind and ejecta inputs, a numerical convergence check shown as shaded bands in Fig. 5, and a transparent statement of the post-processing assumptions and their limitations.","major_comments":[{"comment":"The NEI timescale assignment through the pressure threshold Pth = 3.16e-7 dyn/cm2 is load-bearing for the computed band luminosities, but the paper only verifies that modest variations do not change the identified post-shock morphology; it does not report how the 0.5-2.0 keV and 2.0-10.0 keV light curves, or the RS Oph comparison in Fig. 5, respond to Pth variations. Since vvpshock fluxes are strongly nonlinear in the ionization timescale, and since the same threshold separates young hard-emitting gas from old soft-emitting gas, the claimed order-of-magnitude agreement and the episodic hard flares could be artifacts of the mask rather than of the hydrodynamics. I request a sensitivity study over at least a factor of a few in Pth, with the resulting band luminosities shown, or a tracer-based alternative that does not rely on a single pressure cut.","section":"Appendix A, Fig. 5"},{"comment":"The hydrodynamics uses the CIE cooling curve of Schure et al. (2009), while the X-ray post-processing uses the NEI vvpshock model. The paper acknowledges this inconsistency and argues that dust-grain cooling dominates at T ≳ 2×10^6 K, reducing its impact at high temperatures, but the diffuse gas that dominates the soft band is often at lower temperatures where gas-phase CIE cooling is the only cooling channel. The paper does not quantify the resulting error in the thermal structure or in the soft-band luminosity. A quantitative estimate, or a test run with NEI-aware cooling, is needed to show that the central light curves are not strongly affected by this inconsistency.","section":"§2.2"},{"comment":"The RS Oph comparison rests on a single fiducial parameter set (Mej, Ek, wind mass-loss rate, recurrence cadence), with no exploration of the plausible ranges of these quantities. Because the RS Oph eruption cadence is used as an input and the comparison is made to the same system, the 'broadly consistent' claim would be substantially strengthened by a small parameter-space survey (for example, varying Mej or the wind mass-loss rate by factors of a few) showing that the soft-band secular decline and the order-of-magnitude luminosity level are robust.","section":"§3.2, Table 1"}],"minor_comments":[{"comment":"The velocity profile v(r) = v∞ (r/Rw) is written without an explicit piecewise specification; it would be clearer to state that this form holds for r < Rw and that v = v∞ for r ≥ Rw.","section":"Eq. (1)"},{"comment":"The color-bar labels use 'log fx' with units erg s−1 cm−2, but these are projections through the computational box; please clarify whether these are projected surface-brightness maps or per-cell flux values, and define the line-of-sight integration if used.","section":"Captions of Figs. 3 and 4"},{"comment":"The inset labels 1–9 identify the eruption cycles, but the caption does not fully explain which curve corresponds to which eruption or how the labels map to the post-eruption time axis; please expand the caption for readability.","section":"Inset of Fig. 5"},{"comment":"The 'no subsequent wind mass or energy injection' simplification is important for the late-time density in the excavated cavity and for the soft-band emission measure; it would be helpful to discuss in the conclusions how ongoing wind replenishment could alter the secular decline.","section":"§2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A and makes a useful contribution to recurrent-nova remnant modeling. The central concern is the sensitivity of the quantitative X-ray predictions to the NEI pressure threshold; this is fixable with additional post-processing tests. The use of the RS Oph cadence as an input and then comparing with RS Oph makes the comparison less independent, but this is not circular in the technical sense and should be stated more explicitly in the text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is the full multi-cycle simulation: 130 years, nine eruptions, 3D hydrodynamics with adaptive mesh, post-processed into soft and hard X-ray bands. That is real progress. Prior work mostly modeled isolated outbursts or applied the nested-shell idea to supernova remnants. The paper shows clearly that repeated eruptions build a bipolar cavity bounded by nested shells, that soft X-rays trace dense shell interfaces, and that hard X-rays flare when fresh ejecta catch up with older shells. The secular decline in the soft band and the order-of-magnitude agreement with RS Oph's extended luminosity are plausible, and the convergence check in Figure 5 is a credit to the authors. They are also upfront about missing magnetic fields, conduction, and the single fiducial setup.\n\nWhere the paper is soft is exactly where the reader's report puts the finger: the ionization clock. The quantitative X-ray output depends on a pressure threshold in Appendix A that separates 'fresh' from 'old' ejecta. Cells above the threshold get the time since the last eruption; all other cells get one mean elapsed time over all prior eruptions. That is a crude prescription. The authors verify that modest threshold variations do not change the morphology of the identified post-shock layer, but they never report how the band luminosities change when the threshold moves. Since vvpshock is strongly nonlinear in the ionization timescale, a wrong clock in a small but dense volume can shift soft and hard fluxes by factors. The 'mean elapsed time' for old gas is also unmotivated, because most of it was shocked at specific earlier times and then expanded. These are fixable with sensitivity runs and tracer-based assignments, but until then the RS Oph comparison is not as tight as the abstract implies. It is a single fiducial run with several empirically tuned inputs, including the eruption cadence itself. That does not invalidate the result, but it does make the quantitative agreement less convincing.\n\nSmaller points: the hydrodynamics cools via CIE while the spectra use NEI, a mismatch the authors justify by dust-dominated cooling; that is reasonable but not fully settled. No code or data are released, which hurts reproducibility given the number of free parameters.\n\nOverall, this is a serious paper with a new simulation and a clear qualitative result. It deserves peer review. A good referee should push for a threshold sensitivity study and a better-motivated clock for the older gas before the quantitative claims are accepted. I would cite it if I worked on shocks or cataclysmic variables, and I would bring it to our reading group.","headline":"A genuinely new 130-year, nine-eruption 3D simulation of a recurrent-nova remnant with X-ray post-processing; the RS Oph comparison is plausible but rests on an unquantified NEI pressure-threshold clock.","tokens_in":12537,"tokens_out":1885,"would_cite":true,"duration_ms":17034,"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":"Repeated eruptions build nested shells that keep a recurrent nova's extended X-ray glow alive for more than a century.","keywords":["recurrent novae","nova remnants","X-ray emission","shock waves","circumstellar medium","hydrodynamic simulations","RS Ophiuchi","nested shells"],"falsifier":"Observe RS Oph's extended emission with a high-resolution X-ray spectrum several years after the 2021 eruption and measure the ionization timescale of the fresh shell; if the ionization age does not match the time since 2021 while the diffuse 0.5-2 keV luminosity continues to fade along the simulated two-order-of-magnitude decline, or if a hard flare appears at a time not tied to the eruption cadence, the nested-shell explanation is falsified.","tokens_in":11478,"feed_emoji":"🔭","tokens_out":7719,"duration_ms":61643,"temperature":0.7,"pith_summary":"This paper tries to show that the extended X-ray glow seen around recurrent novae like RS Ophiuchi is not a one-time blast artifact but a long-lived structure built by many eruptions. Using three-dimensional hydrodynamics of nine eruptions over 130 years, the authors find that each eruption deposits a fast bipolar shell into a red-giant wind, and these shells pile up as nested, cavity-bounded walls. Soft X-rays come mostly from the dense compressed shell interfaces and fade slowly as the whole remnant expands, while hard X-rays come from the hottest gas and flare within a few years of each new eruption. The simulated diffuse soft X-ray luminosity lands within an order of magnitude of the value inferred for RS Oph, which is what connects the simulation to observation. If right, diffuse shell-shock X-rays should be counted as a real, persistent component of recurrent-nova remnants, not just a short post-outburst phase.","feed_headline":"Repeated nova blasts keep X-ray glow alive for over a century","feed_subtitle":"Simulations match RS Oph's extended soft X-ray brightness and tie its fading to nested shell expansion.","key_machinery":"The load-bearing mechanism is the nested-shell prescription: each nova is inserted as a thin, clumpy shell with a latitude-dependent expansion speed (faster toward the poles, set by shape parameters $\\alpha$ and $\\beta$) that expands into a wind-shaped density profile before the next eruption deposits another shell on top. The resulting remnant is post-processed with a non-equilibrium-ionization (NEI) shock model, and the ionization clock for each cell is assigned by a pressure threshold: cells above $3.16\\times10^{-7}$ dyn cm$^{-2}$ are treated as freshly shocked ejecta with time since the last eruption, while all other cells get one mean elapsed time. That assignment is what converts the hydrodynamics into soft and hard band luminosities.","core_discovery":"The central claim is that repeated eruptions in a symbiotic recurrent nova assemble a bipolar cavity bounded by nested shells, and that this structure itself radiates the extended X-ray emission seen between outbursts. In the fiducial nine-eruption sequence, the soft 0.5--2.0 keV luminosity peaks near $4\\times10^{32}$ erg/s about 11 years after the first eruption, then declines secularly by roughly two orders of magnitude by 128 years, while the hard 2--10 keV luminosity rises and falls episodically with peaks after each eruption, from about $1.25\\times10^{30}$ erg/s after the first event down to $10^{27}$--$10^{28}$ erg/s in later cycles. The brightest simulated cycles reach the same order of magnitude and decline rate as the extended 0.5--1.8 keV luminosities inferred from RS Oph at 3.4 and 5.3 years after the 2006 eruption. The paper concludes that nested-shell evolution can sustain long-lived diffuse soft X-ray emission between recurrent-nova eruptions.","pith_inferences":["If this picture holds, late-time X-ray spectra of RS Oph should show non-equilibrium ionization signatures whose ionization age tracks the time since the most recent eruption; measuring those ages directly would test the pressure-threshold assignment used here.","The nine-cycle sequence is likely a lower bound for systems like RS Oph, and extending the same physics to many more eruptions would predict an even smoother, fainter, more volume-filling remnant that could be compared with the reported large-scale super-remnant around RS Oph.","The pressure threshold used to tag fresh ejecta is the least physically anchored step; an alternative tracer such as shock velocity or a self-consistent ionization-evolution calculation would reveal whether the band luminosities are stable or threshold-sensitive.","The same nested-shell mechanism may apply to other symbiotic recurrent novae with inferred bipolar structure, giving a uniform prediction that their extended X-ray luminosities should fade secularly between eruptions."],"forward_implications":["Diffuse soft X-ray emission can persist between eruptions for more than a century, so recurrent-nova remnants should be visible as extended soft X-ray sources even when no outburst is in progress.","Hard-band X-ray flares follow each eruption within a few years, with amplitudes that decline as the shell system grows, giving a predicted decade-scale variability pattern tied to eruption cadence.","The nested-shell remnant becomes smoother and more volume-filling over time, so late-time X-ray morphology is expected to look diffuse rather than shell-like.","The accumulated shell system sets the density structure that any future explosion, including a Type Ia supernova, would expand into.","The soft-band secular decline is driven mainly by the falling emission measure of the expanding remnant, not by the details of individual eruptions."],"supporting_citations":[{"why":"Supplies the documented RS Oph eruption epochs that define the fiducial recurrence cadence.","marker":"Schaefer 2010"},{"why":"Adds the 2021 RS Oph outburst to the cadence and links the sequence to modern observations.","marker":"Page et al. 2022"},{"why":"Provides the extended X-ray luminosity measurements of RS Oph that the simulated soft-band luminosities are compared against.","marker":"Montez et al. 2022"},{"why":"Provides the nested-shell ejecta prescription and the wind-injection setup adapted here to recurrent novae.","marker":"Serrano-Hernández et al. 2025"},{"why":"Defines the NEI plane-parallel shock model used for X-ray post-processing.","marker":"Borkowski et al. 2001"},{"why":"Gives the RS Oph elemental enhancement and depletion factors used in the emissivity calculation.","marker":"Orio et al. 2023"},{"why":"Supplies the collisional-ionization-equilibrium cooling curve and the electron-to-hydrogen density ratio used in the hydrodynamics and post-processing.","marker":"Schure et al. 2009"},{"why":"Supplies the clumpy-ejecta perturbation method and the dust-induced cooling treatment used in the simulations.","marker":"Martínez-González et al. 2018"}],"fun_headline_variants":["Nested nova shells sustain X-ray glow for decades","Recurrent nova blasts create X-ray-lit nested shells","Simulations link nova shell stacking to RS Oph's X-rays","Repeated eruptions build X-ray-emitting nova remnant shells","Nova shell structure drives long-term X-ray variability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation leans on a pressure threshold that labels which gas is freshly shocked and which is old, but that threshold is calibrated by matching shock morphology rather than by a physical measurement; if it misplaces fresh ejecta, the soft and hard X-ray fluxes are systematically wrong even though the gas dynamics may be right.","fun_headline_variants_meta":{"raw":{"variants":["Nested nova shells sustain X-ray glow for decades","Recurrent nova blasts create X-ray-lit nested shells","Simulations link nova shell stacking to RS Oph's X-rays","Repeated eruptions build X-ray-emitting nova remnant shells","Nova shell structure drives long-term X-ray variability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000212,"raw_usage":{"total_tokens":1494,"prompt_tokens":1099,"completion_tokens":395,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":715,"completion_tokens_details":{"reasoning_tokens":315}},"tokens_in":715,"tokens_out":395,"duration_ms":4090,"temperature":1.0,"reasoning_tokens":315,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:33:33.780262+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe RS Oph's extended emission with a high-resolution X-ray spectrum several years after the 2021 eruption and measure the ionization timescale of the fresh shell; if the ionization age does not match the time since 2021 while the diffuse 0.5-2 keV luminosity continues to fade along the simulated two-order-of-magnitude decline, or if a hard flare appears at a time not tied to the eruption cadence, the nested-shell explanation is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the documented RS Oph eruption epochs that define the fiducial recurrence cadence."},{"cited_title":"L., Beardmore, A","cited_arxiv_id":null,"evidence_quote":"Adds the 2021 RS Oph outburst to the cadence and links the sequence to modern observations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the extended X-ray luminosity measurements of RS Oph that the simulated soft-band luminosities are compared against."},{"cited_title":"2023, ApJ, 955, 37","cited_arxiv_id":null,"evidence_quote":"Gives the RS Oph elemental enhancement and depletion factors used in the emissivity calculation."},{"cited_title":"M., Kosenko, D., Kaastra, J","cited_arxiv_id":null,"evidence_quote":"Supplies the collisional-ionization-equilibrium cooling curve and the electron-to-hydrogen density ratio used in the hydrodynamics and post-processing."}],"review_version":2}