{"id":"70b84408-8a08-430f-bda3-37941c9c52ec","arxiv_id":"2506.21667","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Supernova blast waves can resupply planetary debris to old white dwarfs by kicking exo-Oort cloud material inward, with metre-sized boulders delivered at least once over a cooling age.","lead":"Supernova blast waves can fling distant debris around old white dwarfs into the inner planetary region, giving a possible new source for the rocky pollution seen in very old white dwarfs. The paper derives the geometry, kick sizes, and debris-size limits, and predicts that metre-sized boulders should be delivered at least once over a white dwarf's cooling age.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (5)'s energy-coupling kick likely overestimates boulder Δv by orders of magnitude; gamma=1 is not merely uncertain but momentum-inconsistent for large bodies.","rationale":"The central claim is the resupply of metre-sized boulders. The reader correctly identified Section 3's energy-coupling assumption as the weakest point, but the issue is more acute than an uncertain gamma. Equation (4) treats the body as a sail that converts intercepted blast energy into its own kinetic energy with efficiency gamma. Momentum conservation forbids this for a body that intercepts a tiny gas mass: the blast's energy is carried by gas whose density is ~1 cm^-3, so a 1 m boulder at 20 pc sweeps only ~10^-4 kg of gas; even 100% elastic momentum transfer gives Δv ~ 0.1 m/s, whereas the paper's gamma=1 formula gives ~200 m/s and the orbital threshold is ~75 m/s. In the standard ram-pressure treatment, the kick scales as (Σ_gas v_shock)/(ρ_body R), not sqrt(ESN/(ρ_body R D^2)); using Sedov parameters yields Rmax ~ 1 mm for the Table 1 'Yes' configurations, not 7.5 m. This does not invalidate the paper's geometric probability derivations or its conclusions for micron dust and mm pebbles, which are accelerated efficiently by either model. But the headline boulder claim is unsupported unless the coupling model is replaced with one that respects momentum conservation. The proposed concrete test—comparing momentum-based Rmax against Eq. (47) or running a blast-wave hydro simulation—would settle this. Until then, CONDITIONAL remains appropriate; if the test confirms the drag-based estimate, the boulder claim should be removed or drastically revised, reducing the paper's main conclusion to dust/pebble resupply.","tokens_in":23369,"tokens_out":21883,"duration_ms":229767,"concrete_test":"Compute the momentum-based kick for radii R=1 μm to 10 m using the standard ram-pressure impulse Δv ≈ 3 Σ_gas v_shock / (4 ρ_body R), with Σ_gas = (1/3) ρ_ISM D_SN and v_shock from the Sedov-Taylor solution for E_SN=10^44 J. Overlay these curves on Fig. 6 and compare with the minimum kick Eq. (32) for each Table 1 row. If the momentum-based Rmax is below 1 m for D_SN ≥ 20 pc while Eq. (5) gives ≥1 m, the boulder conclusion fails. An independent check is to run a 1D or 3D hydrodynamic simulation of a 1 m sphere at 20 pc from a 10^44 J blast and measure its velocity change.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The load-bearing weak point is the Section 3 kick prescription. Equation (4) assumes a fraction γ of the blast energy intercepted by the body's geometric cross-section becomes translational kinetic energy, with γ=1 adopted by fiat. For a macroscopic solid body in a blast wave, this is not just an uncertain parameter but the wrong scaling: momentum is transferred by ram pressure from the tenuous shocked gas, giving Δv ~ C_d Σ_gas v_shock / ρ_body (per unit radius), where Σ_gas ~ (1/3)ρ_ISM D_SN is the shell column density. For the fiducial D_SN=20 pc, v_shock ≈ 350 km/s and Σ_gas ≈ 3×10^-4 kg/m^2, so a 1 m radius boulder (ρ≈1500 kg/m^3) receives Δv ≈ 0.1 m/s—about three orders of magnitude smaller than the ~75 m/s threshold from Eq. (32) for the Table 1 'Yes' rows. The energy formula Eq. (5) instead gives Δv ≈ 200 m/s because it assumes the body absorbs all incident kinetic energy, violating momentum conservation for a body whose swept-up gas mass is tiny. Equation (47)'s Rmax ∝ γ is therefore an artifact of the assumed coupling; a physically derived Rmax from direct gas drag is millimetre-scale, not metre-scale. The 'metre-sized boulders resupplied at least once' headline claim thus rests on an unphysical coupling model, even before the value of γ is debated. Dust and mm pebbles remain plausibly resupplied by either model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript develops an analytic framework for the impulsive action of a supernova blast wave on small bodies in white-dwarf exo-Oort clouds. It derives post-blast orbital elements, the fraction of blast geometries that lower a body's pericentre into the ~100 au perturbation zone, the minimum kicks needed for bounded, leaking, and escaping post-blast orbits, and a 23 per cent upper bound on the initial true anomalies that permit repeated pericentre passages. It then combines the kick threshold with an energy-coupling model for the blast to obtain a maximum deliverable body size Rmax and, using the local supernova rate, estimates the number of inward pericentre thrusts over a white dwarf cooling age, concluding that dust and millimetre pebbles are always affected and that metre-sized boulders are resupplied at least once to very old white dwarfs.","tokens_in":23689,"tokens_out":10837,"duration_ms":126961,"significance":"The orbital-mechanics core is a genuine strength: the analytic geometry predictions agree with three Monte Carlo suites to within a few per cent (Section 8), the 23 per cent bound is derived rather than calibrated, and the only fitted element is the correction factor k≈0.87. If the kick-coupling assumption were physically justified, the mechanism would offer a quantitative route for maintaining pollution in old and very old white dwarfs without long-lived inner reservoirs, and the paper's analytic expressions would be reusable. The main weakness is that the size and rate conclusions are controlled by an energy-coupling prescription in Section 3 that is explicitly unvalidated and, as argued below, has the wrong physical scaling for macroscopic bodies.","major_comments":[{"comment":"The energy-coupling prescription is load-bearing and, as written, physically problematic for the size regime the paper emphasizes. Equation (4) assumes that a fraction γ of the blast energy intercepted by the body's geometric cross-section is converted into translational kinetic energy of the body, and Eq. (5) then gives Δv ∝ [γ E_SN/(ρ R D_SN^2)]^{1/2}. For a macroscopic solid body in a blast wave, however, momentum is transferred by the ram pressure of the shocked gas, giving Δv ~ C_d Σ_gas v_shock/(ρ_body R), with Σ_gas the shell column density. For the fiducial parameters of Section 9 (D_SN=20 pc, R=1 m, ρ=1500 kg/m^3, v_shock≈350 km/s) this yields Δv ~ 0.05 m/s, about three orders of magnitude below the ~75 m/s threshold implied by Eq. (32), whereas Eq. (5) with γ=1 yields Δv≈200 m/s. Equation (47) and Table 1 inherit this problem: Rmax scales linearly with γ, and the 'metre-sized boulders' conclusion is an artifact of an energy-coupling model that violates momentum conservation for bodies much more massive than the intercepted gas mass. The manuscript itself notes the lack of physical justification for γ=1; I recommend replacing Eq. (5) with a momentum-drag-based kick and recomputing Rmax, Table 1, and the abstract's size claims.","section":"Section 3, Eqs. (4)-(5); Section 9.1, Eq. (47); Table 1"},{"comment":"Equation (48) defines N as the number of inward pericentre thrusts per small body over a cooling age, and the Table 1 entries for D_SN=20 pc give N ≈ 0.025-0.15 t_cool Γ. With Γ≈0.4 Gyr^-1 (one 20 pc supernova per 2.5 Gyr) and t_cool=10 Gyr, N is below unity for every row, so the statement in the abstract that 'metre-sized boulders [are] resupplied at least once' is a statement about a population of many boulders rather than about an individual body. Because the manuscript does not specify the exo-Oort cloud boulder population or its size distribution, the population-level probability of at least one resupply event is not computed; the abstract overstates what Eq. (48) and Table 1 can support. I recommend either adding a population model or qualifying the claim to per-object probabilities.","section":"Section 9.2 and Table 1"}],"minor_comments":[{"comment":"The phrase 'A then relevant question' should read 'A relevant question'.","section":"Section 5.2"},{"comment":"There are typographical errors: 'repleneshed' in Section 1 should be 'replenished', and 'exo-Oort cloulds' in Section 10 should be 'exo-Oort clouds'.","section":"Section 1 and Section 10"},{"comment":"The sentence 'I assume that D_SN does not vary across the white dwarf cooling age' is in tension with the preceding paragraph, which emphasizes that Γ depends on D_SN(t); the assumption should be explicitly labelled as a simplification for the estimates in Table 1.","section":"Section 9.2"},{"comment":"The right-hand y-axis label 'Number of supernovae needed' is the reciprocal of the plotted probability, but this is not stated; please clarify that the number of supernovae is 1/P and assumes independent blasts.","section":"Figure 1"},{"comment":"The word 'prove' is strong for a result obtained with series expansions and the small-qf/ai approximation; consider using 'show' or 'demonstrate'.","section":"Abstract and Section 7.6"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of MNRAS and the orbital-mechanics part is solid. The main obstacle is the unvalidated and physically implausible energy-coupling model for macroscopic bodies; if the author replaces it with a momentum-transfer prescription and recomputes the size and rate claims, a revised version would likely be publishable. I do not see citation or novelty concerns beyond what is stated in the report."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first analytic treatment of supernova blast waves as a way to resupply debris to old white dwarfs, and the orbital mechanics part is genuinely solid. The size-dependent headline—metre-sized boulders resupplied at least once—rests on a gamma=1 energy coupling that is not physically justified, and the stress-test note is right that this is not just an uncertain parameter but the wrong scaling for a macroscopic boulder. For dust and mm pebbles the mechanism survives; for boulders, probably not.\n\nWhat's new: the probability integrals for blast geometries (eqs 27/30), the kick thresholds for inward thrust, leaking and breaking (eqs 32-34), the 23% bound on true anomalies that allow repeated pericentre passages, and the Rmax formula (eq 47). The three Monte Carlo tests check out to within a few per cent, which gives confidence in the analytic machinery. The paper also does a good job of linking to the local supernova rate and tabulating expected numbers of events. This is a real contribution to the white dwarf pollution literature.\n\nWhere it's soft: the energy coupling in Sec 3. Adopting gamma=1 by fiat, with no physical justification, converts a small-body blast interaction into pure energy capture. The stress-test calculation shows that momentum transfer from the shocked gas alone gives a couple of orders of magnitude less delta-v for a 1 m boulder than eq (5). That makes the 'metre-sized boulders resupplied at least once' claim look like an artifact of the coupling model. The 23% bound is an approximate numerical maximum, not a theorem, though the author does phrase it 'within this formalism', so that's a minor complaint. Also, the paper does not consider radiation-driven effects, but those would mostly reduce the coupling further.\n\nBottom line: the geometry and kick-threshold derivations are worth refereeing and likely worth publishing after the coupling issue is addressed. The author should do a sensitivity analysis over gamma or a physically motivated coupling model, and the boulder claim should be downgraded until then. Send to a serious referee, not desk reject.","headline":"A genuinely new analytic treatment of supernova debris resupply for old white dwarfs, with a solid orbital-mechanics core, but the headline boulder claim rests on an unjustified gamma=1 energy-coupling assumption.","tokens_in":24244,"tokens_out":3678,"would_cite":false,"duration_ms":37336,"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":"Supernova blast waves can resupply debris to old white dwarfs from their exo-Oort clouds, delivering metre-sized boulders at least once over a cooling age.","keywords":["planets and satellites: dynamical evolution and stability","planet-star interactions","stars: white dwarfs","stars: supernovae: general","ISM: supernova remnants","celestial mechanics","exo-Oort clouds","white dwarf pollution"],"falsifier":"Run a radiation-hydrodynamics simulation of a supernova blast wave sweeping over a metre-sized boulder at $10^4$ au and measure the fraction of intercepted energy that becomes bulk translation; a value of $\\gamma$ much below 1 would push $R_{\\mathrm{max}}$ in equation (47) below the metre scale and reduce the expected boulder delivery rate below once per cooling age.","tokens_in":23103,"feed_emoji":"💥","tokens_out":9524,"duration_ms":96124,"temperature":0.7,"pith_summary":"Old and very old white dwarfs show signs of polluting debris, but their inner reservoirs (analogues of the asteroid and Kuiper belts) are depleted long before a 10 Gyr cooling age. This paper argues that supernova blast waves can resupply those systems from the outside: an impulsive kick from a nearby supernova can thrust small bodies orbiting at $10^4$ au in an exo-Oort cloud inward to within $100$ au of the white dwarf, where the usual planet-driven delivery can take over. The author derives the blast geometries, kick magnitudes, and resulting orbit types analytically, and connects them to debris size and the local supernova rate. The bottom line is a size ladder: micron dust and millimetre pebbles are ejected or redirected by essentially every relevant blast, objects above $10$ km are almost never delivered, and metre-sized boulders should be resupplied at least once to very old white dwarfs over their cooling ages. If correct, this provides a quantitative route for sustaining observable white-dwarf pollution without a long-lived inner reservoir.","feed_headline":"Supernova blasts fling metre boulders to old white dwarfs","feed_subtitle":"Blast waves sweep Oort-cloud debris inward, explaining how 10-billion-year-old white dwarfs stay polluted.","key_machinery":"The central machinery is the impulse approximation for a supernova blast, adopted from Jackson et al. (2014): the blast is treated as an instantaneous velocity kick $\\Delta v = \\sqrt{3\\gamma E_{\\mathrm{SN}}/(4\\pi \\rho R D_{\\mathrm{SN}}^2)}$ applied to a small body on an eccentric orbit around a $0.6\\,M_\\odot$ white dwarf, with $\\gamma$ the fraction of intercepted blast energy converted to translational kinetic energy (taken as 1). The kick direction is parameterised by polar angle $\\theta$ and azimuthal angle $\\phi$, and the post-blast semi-major axis, eccentricity, pericentre and apocentre are given by closed-form expressions. The argument then hinges on comparing three minimum kick thresholds — the kick needed to push the pericentre into the perturbation zone, the kick that makes the orbit leak past the Hill ellipsoid, and the kick that breaks the ellipse into a hyperbola — which together bound the geometries and true anomalies that permit repeated inner passages and set the maximum deliverable size $R_{\\mathrm{max}}$.","core_discovery":"The central claim is that a supernova blast wave acts as an impulsive velocity kick, $\\Delta v$, on small bodies in an exo-Oort cloud, and that this kick can shrink their pericentre into the inner $\\sim 100$ au perturbation zone where debris can eventually be accreted. The paper proves, within its impulse formalism, that for a post-blast elliptical orbit to stay entirely inside the white dwarf's Hill ellipsoid and yield repeated pericentre passages, the pre-blast true anomaly must lie in a restricted interval around apocentre; the maximum fraction of true anomalies that allow this is about 23 per cent. It further derives the maximum debris radius that can be delivered, $R_{\\mathrm{max}} \\propto \\gamma E_{\\mathrm{SN}} a_i (1-e_i^2)/[\\rho M_\\star D_{\\mathrm{SN}}^2 (1+e_i \\cos f_i)^2]$, and, combining this with a local supernova rate, concludes that micron dust and millimetre sand and pebbles are redirected or ejected in nearly every blast, metre boulders are resupplied at least once over a 10 Gyr cooling age, and asteroids larger than about 10 km are essentially never delivered unless the supernova goes off inside the cloud itself.","pith_inferences":["The paper's size ladder implies a time-ordering of pollution that is not spelled out: dust and pebbles are processed quickly by blasts, while metre boulders arrive on longer timescales, so younger-old versus very-old white dwarfs should show different characteristic grain sizes in their debris — a testable prediction against observed debris-disc spectral energy distributions.","Because the maximum deliverable radius scales linearly with the energy-coupling fraction $\\gamma$, the metre-boulder conclusion is the least robust part of the paper: if real blast coupling is an order of magnitude below unity, the delivered objects are decimetre-scale, not metre-scale.","The 23 per cent bound on true anomalies suggests that many resupplied objects make only a single pass through the inner system; observable pollution from this channel may therefore be episodic or transient rather than steady, unless multiple blasts or replenished clouds keep resupplying the perturbation zone.","A white dwarf's Galactic trajectory controls the local supernova rate, so the mechanism predicts that pollution incidence among old white dwarfs should correlate with kinematics (thin-disc versus thick-disc or halo populations); Gaia-style kinematic samples could test this correlation statistically."],"forward_implications":["Micron-sized dust and millimetre-sized sand and pebbles in exo-Oort clouds around old white dwarfs are ejected or have their orbits significantly altered by essentially every local supernova blast; their size distribution in those clouds will become top-heavy unless continuously replenished.","Metre-sized boulders are, more likely than not, thrust into the inner perturbation zone at least once over the cooling age of a very old white dwarf.","Objects larger than about 10 km can only be delivered by a supernova occurring inside the exo-Oort cloud itself; for typical blast distances they are effectively unaffected.","Repeated pericentre passages — the condition for sustained delivery — require the small body to be near apocentre before the blast, and even then only at most about 23 per cent of true anomalies allow this; otherwise the post-blast orbit leaks or becomes hyperbolic and the debris passes through the inner region at most once.","Old white dwarfs can therefore maintain observable pollution without relying on long-lived inner debris belts, as long as an exo-Oort cloud exists and the system experiences sufficiently frequent nearby supernovae along its Galactic path."],"supporting_citations":[{"why":"Supplies the impulse formalism and the post-blast orbital element equations (7), (10)-(13) that the paper's kick thresholds are built on.","marker":"Jackson et al. (2014)"},{"why":"Numerically demonstrates that supernovae can perturb Oort cloud material into the inner solar system and sets the blast duration (~0.05 Myr) used to fix the impulse regime, with gamma = 1/2 as a comparison case.","marker":"Smith et al. (2024)"},{"why":"Provides the census of core-collapse supernova rates in the Solar neighbourhood used to convert kick thresholds into the expected number of inward pericentre thrusts over a cooling age.","marker":"Quintana et al. (2025)"},{"why":"Gives evidence that nearly two dozen supernovae occurred within 300 pc of the Sun in the last 300 Myr, supporting the assumed frequency of relevant nearby blasts.","marker":"Firestone (2014)"},{"why":"Sets the Galactic Hill ellipsoid scale (~149,000 au for a Solar analogue) from which the escape boundary resc = 125,000 au is adopted.","marker":"Veras & Evans (2013)"},{"why":"Documents white dwarfs that have been cooling for ~10 Gyr and still show pollution, defining the very-old population the resupply mechanism targets.","marker":"Elms et al. (2022)"}],"fun_headline_variants":["Supernova blast waves kick boulders to ancient white dwarfs","Supernovae resupply boulders to old white dwarfs","Blast waves fling boulders to pollute ancient white dwarfs","Supernova kick delivers boulders to ancient white dwarfs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire size and rate argument assumes that all ($\\gamma = 1$) of the supernova blast energy intercepted by a small body's geometric cross-section is converted into translational kinetic energy; if the real coupling is even a tenth of that, the maximum deliverable radius shrinks by an order of magnitude and the metre-boulder conclusion can fail.","fun_headline_variants_meta":{"raw":{"variants":["Supernova blast waves kick boulders to ancient white dwarfs","Supernovae resupply boulders to old white dwarfs","Blast waves fling boulders to pollute ancient white dwarfs","Supernova kick delivers boulders to ancient white dwarfs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001649,"raw_usage":{"total_tokens":6650,"prompt_tokens":1147,"completion_tokens":5503,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":763,"completion_tokens_details":{"reasoning_tokens":5430}},"tokens_in":763,"tokens_out":5503,"duration_ms":39435,"temperature":1.0,"reasoning_tokens":5430,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:22:01.157809+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a radiation-hydrodynamics simulation of a supernova blast wave sweeping over a metre-sized boulder at $10^4$ au and measure the fraction of intercepted energy that becomes bulk translation; a value of $\\gamma$ much below 1 would push $R_{\\mathrm{max}}$ in equation (47) below the metre scale and reduce the expected boulder delivery rate below once per cooling age.","supporting_citations":[{"cited_title":"P., Wyatt, M","cited_arxiv_id":null,"evidence_quote":"Supplies the impulse formalism and the post-blast orbital element equations (7), (10)-(13) that the paper's kick thresholds are built on."},{"cited_title":"A., & Fields, B","cited_arxiv_id":null,"evidence_quote":"Numerically demonstrates that supernovae can perturb Oort cloud material into the inner solar system and sets the blast duration (~0.05 Myr) used to fix the impulse regime, with gamma = 1/2 as a comparison case."},{"cited_title":"& Evans, N","cited_arxiv_id":null,"evidence_quote":"Sets the Galactic Hill ellipsoid scale (~149,000 au for a Solar analogue) from which the escape boundary resc = 125,000 au is adopted."},{"cited_title":"K., Tremblay, P.-E., G¨ ansicke, B","cited_arxiv_id":null,"evidence_quote":"Documents white dwarfs that have been cooling for ~10 Gyr and still show pollution, defining the very-old population the resupply mechanism targets."}],"review_version":1}