{"id":"39f7d77d-6280-4ce0-9d62-120b8225b6e5","arxiv_id":"2411.16857","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"The critical neutrino luminosity for a stalled supernova shock depends on the pre-shock Mach number, and the antesonic ratio is the most stable explosion criterion across the tested 1D parameter space.","lead":"Researchers ran 1D simulations of stalled supernova shocks and mapped the neutrino luminosity needed to trigger explosion across many conditions. They found that the thermal content of the gas falling onto the shock matters as much as its density, helping explain why some stars explode.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The low-M normalization shift in Eq. 22 may be partly a numerical artifact because the free-nucleon suppression factor chi_N (Eq. 17) fails to suppress pre-shock heating for M=1.3, as the paper itself states in Section 3.1 and the Figure 17 caption.","rationale":"The reader's weakest-assumption analysis and my independent read converge on the same load-bearing point: the paper's new physics claim about M-dependence of L_crit hinges on the pre-shock flow being isentropic between the outer boundary and the shock, with entropy set by the boundary Mach number. The paper itself provides textual evidence that this condition fails precisely in the regime (M=1.3) where the largest effect is claimed. The statement in Section 3.1 that 'this treatment is less effective for M=1.3, where Qdot slowly approaches zero in the pre-shock region' and the Figure 17 caption's explicit admission of an inadequate suppression are internal flags that cannot be dismissed. A numerical artifact here would not only change the normalization by a few percent; it could systematically bias the 44% shift that is the central new result and the physical connection to compositional interfaces. My concrete test is designed to isolate this effect by removing the contamination entirely. I do not believe this concern is strong enough to reject the paper or to demand a more severe verdict than CONDITIONAL, because the contamination is acknowledged, the affected regime is clearly identified, and the main scaling relation (Eq. 22) is anchored at M=2.0 where the suppression works. However, the concern is real and must be quantified before the M-dependence is used for progenitor-level conclusions. The reader's conditional verdict with moderate confidence remains appropriate; no verdict change is needed.","tokens_in":75,"tokens_out":3991,"duration_ms":51488,"concrete_test":"Re-run the M=1.3 series at two accretion rates (e.g., Mdot=0.5 and 1.0 Msun/s) with the free-nucleon suppression factor chi_N in Eq. 17 replaced by an exact step function that sets Qdot=0 for all r>R_shock, while keeping all other inputs and the tau=2/3 feedback unchanged. If L_crit for these models rises by more than ~5% toward the M=2.0 values, the claimed M-dependence of the critical luminosity is partly contaminated by the pre-shock heating artifact. As a supporting diagnostic, compare the entropy at R_shock with the entropy at the outer boundary: if the difference for M=1.3 exceeds the ~18% variation quoted in Figure 17, the outer-boundary Mach number does not uniquely determine the thermal content that reaches the shock.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most novel quantitative claim is that L_nu^crit drops by roughly 44% as the outer boundary Mach number decreases from 2.0 to 1.3 (Table 1), and this is interpreted as higher-entropy accreted material, i.e., compositional interfaces, lowering the critical luminosity (Section 3.6). This interpretation relies entirely on the pre-shock thermal content being controlled by the single outer-boundary parameter M via Eq. 21. However, the same section acknowledges that for M=1.3 the suppression factor chi_N (Eq. 17) does not drive the pre-shock neutrino heating to zero; instead Qdot 'slowly approaches zero in the pre-shock region' (Figure 1 caption), and Figure 17 shows a nonzero pre-shock entropy gradient for M=1.3 because 'Qdot is inadequately suppressed by equation 17 in this parameter space.' In these low-M models, spurious heating between the outer boundary and the shock adds entropy to the infalling flow, so the material that actually crosses the shock does not have the thermal content implied by the boundary condition. If this contamination lowers the effective ram pressure or pre-heats the accreted gas, then part of the measured 44% decrease in L_crit is a numerical artifact rather than a physical effect of pressurized inflow. This would weaken the direct connection drawn to compositional-interface accretion and would also affect the stated comparison to multi-dimensional instability reductions (~30%) made in Section 3.2.1.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents spherically symmetric, time-dependent Athena++ accretion models of the stalled core-collapse supernova shock with a general equation of state and neutrino heating/cooling. Varying the accretion rate, neutrino luminosity, average neutrino energy, PNS radius and mass, and the pre-shock Mach number at fixed neutrino optical depth tau = 2/3, the authors compute the critical neutrino luminosity L_nu^crit separating accretion from explosion. They fit the numerical critical curves to a power-law scaling relation (Eq. 22) with Mdot-dependent exponents (Eqs. 23-25), and they test three proposed explosion diagnostics: the antesonic condition, the advection/heating timescale criterion, and the force explosion condition. The main qualitative claims are that low pre-shock Mach number lowers the normalization of L_nu^crit and that, across the explored model space, the antesonic ratio shows the smallest relative variation, with a near-critical value of about 0.215 +/- 0.01. The low-M behavior is interpreted as evidence that accretion of higher-entropy material across compositional interfaces, such as the Si/O interface, can promote explosion in addition to the well-known decrease in Mdot.","tokens_in":23856,"tokens_out":7134,"duration_ms":73781,"significance":"If the results hold, this is a useful quantitative map of the 1D critical condition and a natural baseline for the planned 2D and 3D extension. The paper has real strengths: the simulations enforce tau = 2/3 through a feedback loop, use stabilized checkpoint models to step in L_nu, check resolution convergence of L_nu^crit to about 1%, and test several published explosion criteria on the same time-dependent data set. The comparison of the antesonic, timescale, and force conditions is informative even where it shows that none of the criteria is exactly constant. However, the central scaling relation is a least-squares fit with no reported uncertainties, and the low-M branch of the study is acknowledged by the authors themselves to be affected by imperfect suppression of pre-shock neutrino heating. The significance is therefore conditional: the framework and diagnostics are valuable, but the quantitative low-M shift and the fit parameters need additional support before the results can be used as calibrated predictions.","major_comments":[{"comment":"The central product of the paper is the power-law scaling relation for L_nu^crit, but the text reports no fit uncertainties, no goodness-of-fit statistic, and no number of fitted model points. Equations (23)-(25) claim weak Mdot dependence of the exponents based on curves like Figure 8, yet no error bars are shown on those fits, and Table 1 calibrates the normalization for each Mach number using a single point. In addition, Eq. (22) contains no M term at all, even though M is a headline input parameter; the M dependence is handled by separately rescaling the prefactor at three discrete values. I request that the authors report the fit covariances and residuals, state the fitting range and number of points, and either extend Eq. (22) to include M explicitly or clearly restrict its validity to M = 2.0.","section":"Sec. 3.2.1, Eqs. (22)-(25)"},{"comment":"The manuscript explicitly acknowledges that for M = 1.3 the free-nucleon suppression factor chi_N in Eq. (17) does not drive the pre-shock Qdot to zero; Figure 1 shows Qdot slowly approaching zero in the pre-shock region, and the Figure 17 caption states that the M = 1.3 profile has a non-zero entropy gradient because Qdot is inadequately suppressed. This contamination is load-bearing for the paper's main quantitative claim, because the 44% decrease in L_nu^crit between M = 2.0 and M = 1.3 (Table 1, Figure 7) is the basis for the compositional-interface interpretation in Section 3.6. Spurious pre-shock heating can raise the entropy and lower the ram pressure of material entering the shock, thereby lowering L_nu^crit artificially. I request a control experiment in which the heating and cooling rates are artificially set to zero for r > R_shock for all M, together with a quantitative estimate of the integrated spurious pre-shock heating relative to the accretion enthalpy flux, so that the physical part of the low-M shift can be separated from the numerical artifact.","section":"Sec. 3.1 and Fig. 17 caption"},{"comment":"The abstract's statement that the antesonic ratio shows the least variation across the model space is based on Figure 9, but the comparison mixes time-averaged maxima for oscillatory models with final-time maxima for non-oscillatory models, and the bars show that instantaneous values can exceed the critical value without leading to explosion. The quoted range 0.215 +/- 0.01 also includes the M = 1.3 models affected by the pre-shock heating problem discussed above. I ask for a version of Figure 9 and the associated statistics restricted to M >= 2.0, and for a statement of how the time-averaging window affects the 3-5% variation quoted in the conclusions. Without this, the 'least variation' claim is not cleanly separated from the grid of models and diagnostics used to define it.","section":"Sec. 3.3, Figure 9"}],"minor_comments":[{"comment":"There are several typos, including 'relativistc' in Section 2.2 and 'timscles' in Section 3.4; a careful proofreading pass is needed.","section":"General"},{"comment":"The notation \\dot{M}^{1.0} is ambiguous; it should be written as \\dot{M}/(1 M_sun/s) or defined explicitly before use, as is done in the text but not in the equation.","section":"Eq. (34)"},{"comment":"The caption states that R_shock stabilizes at ~150-200 km and ~600 km for models that are not the fiducial Rstar = 30 km models; please clarify in the caption which parameter sets produce these radii, since the range quoted in the text is much smaller.","section":"Figure 5 caption"},{"comment":"The distinction between the first exploding model L_nu^crit and the near-critical stable model L_nu^{crit,n} is clear, but the reader would benefit from an explicit statement of why different sections use one or the other; currently the transition between Sections 3.2 and 3.3 is abrupt.","section":"Section 2.4 and Section 3.2"},{"comment":"The data availability statement says the implementation and data are available upon request; depositing the problem generator and model outputs in a persistent repository would strengthen reproducibility, especially since readers cannot otherwise check the fit values in Eq. (22).","section":"Data availability"},{"comment":"The caption attributes the black profile to Eqs. (28)-(29) and the blue profile to Eq. (30), but the latter is defined together with Eq. (31); please correct the cross-reference.","section":"Figure 12 caption"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the low-M normalization is not a manufactured issue: the authors themselves report in the Figure 1 and Figure 17 captions that chi_N fails to suppress pre-shock heating for M = 1.3, and that is exactly the branch carrying the largest claimed effect. I would therefore ask the editor to require a control calculation before the 44% statement is accepted. I do not see grounds for rejection, because the numerical framework and the qualitative comparison of explosion criteria are solid and the contamination is, in principle, quantifiable. The main revisions are: report fit uncertainties, include M in the scaling relation or restrict its validity, and isolate the numerical contamination in the low-M models."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nThis is a careful 1D parameter study of the critical neutrino luminosity for core-collapse supernova shock revival. The genuinely new pieces are the systematic variation of pre-shock Mach number M and the time-dependent comparison of the antesonic, advection-heating timescale, and force explosion conditions across that parameter space. The simulations are well constructed: resolution convergence to about 1 percent, tau held fixed at 2/3 via ghost-zone density adjustment, and checkpoint stepping through L_nu. The central claim, that pressurized pre-shock flow lowers L_nu^crit (with roughly a 44 percent normalization shift between M=1.3 and 2.0) and that the antesonic ratio is the least variable explosion criterion (~0.215±0.01), is supported by the figures, and the fitted scaling relation (Eq. 22) is a useful summary of the numerics.\n\nThe main soft spot is the one the paper itself flags: for M=1.3 the free-nucleon suppression factor chi_N does not drive pre-shock heating to zero. Figure 17 shows a nonzero pre-shock entropy gradient, and the text says Qdot is inadequately suppressed. That means part of the low-M lowering of L_crit could be a boundary-condition artifact rather than a pure physical effect of higher-entropy inflow. The authors are upfront about it, but they do not quantify how much of the 44 percent shift survives if the spurious heating is removed. That matters because the compositional-interface argument in Section 3.6 leans on the size of the shift. I would not call the central argument wrong—the trend is monotonic and consistent with the earlier Pejcha & Thompson result—but the magnitude is less secure than the abstract suggests.\n\nOther soft spots are minor. Equation 22 is a least-squares fit with no reported uncertainties, and the exponents themselves depend on Mdot, so the relation is descriptive rather than predictive. The data and code are \"available upon request\" rather than public, which limits independent checking. The comparison to ~30 percent multidimensional instability reductions is suggestive but not a quantitative match.\n\nWho gets value: anyone working on critical luminosity phenomenology, explosion conditions, or explodability of progenitors. It is a solid reference result for 1D parameterized models. It deserves a serious referee; I would send it to review and expect the referee to ask for quantified pre-shock heating tests and error bars on the fit.","headline":"A careful 1D parameter study that usefully extends the critical luminosity framework to finite pre-shock Mach number, with a real but openly acknowledged numerical caveat at low M that should temper the headline 44% shift.","tokens_in":24289,"tokens_out":2061,"would_cite":true,"duration_ms":19447,"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":"A fitted scaling relation now predicts the critical neutrino luminosity for supernova explosions from four inputs.","keywords":["core-collapse supernovae","critical neutrino luminosity","antesonic condition","accretion shock","neutrino heating","proto-neutron star","pre-shock Mach number","Si/O interface"],"falsifier":"Run the $\\mathcal{M}=1.3$ series again with pre-shock neutrino heating forcibly set to zero; if $L_\\nu^{\\mathrm{crit}}$ rises back toward the $\\mathcal{M}=2.0$ values, the low-Mach lowering is a numerical artifact rather than the pressurized-inflow effect.","tokens_in":23236,"feed_emoji":"💥","tokens_out":10285,"duration_ms":81564,"temperature":0.7,"pith_summary":"This paper asks what sets the exact boundary between a stalled accretion shock and a successful supernova explosion in the neutrino-heating mechanism. Using time-dependent one-dimensional accretion models with a general equation of state, neutrino heating and cooling, and a controlled Mach number for the pre-shock flow, it derives a scaling relation for the critical neutrino luminosity $L_\\nu^{\\mathrm{crit}}$ as a function of accretion rate, proto-neutron star radius and mass, and average neutrino energy. It finds that pressurized pre-shock inflow lowers the critical luminosity by up to roughly 44 percent relative to free-fall models, connecting this to entropy jumps across compositional interfaces like the Si/O layer. Among three proposed explosion diagnostics, the antesonic ratio varies least along the critical curve, staying near $0.215 \\pm 0.01$. If correct, the relation gives a fast quantitative way to rank progenitors by explodability and prepares the same measurements for multi-dimensional models.","feed_headline":"One scaling relation predicts supernova explosion thresholds","feed_subtitle":"Pre-shock pressure shifts critical luminosity by ~44%, and the antesonic ratio stays nearly constant at explosion.","key_machinery":"The central object is the critical curve in the $(L_\\nu, \\dot{M})$ plane, whose normalization is mapped by raising $L_\\nu$ in small steps until the stalled shock advects outward. The controlling device is the outer-boundary Mach number $\\mathcal{M}$, which sets the pre-shock pressure through $P_{\\mathrm{go}}\\simeq \\rho v^2/(\\Gamma_{\\mathrm{go}}\\mathcal{M}^2)$ and thereby the thermal content of the accreted matter, while the neutrino optical depth is held fixed at $\\tau=2/3$ by adjusting ghost-zone densities. The main diagnostic is the antesonic ratio $\\max(c_s^2/v_{\\mathrm{esc}}^2)$, the squared ratio of post-shock sound speed to escape velocity, whose near-constant value along the critical curve anchors the paper's comparison of explosion conditions.","core_discovery":"Across a grid of spherically symmetric, time-dependent accretion models with a general equation of state and optically thin neutrino heating and cooling, the paper finds that the explosion threshold is a critical neutrino luminosity $L_\\nu^{\\mathrm{crit}}$ and fits it as $$L_\\$nu^{{\\mathrm{crit}}$} = 32.8 \\left(\\frac{\\dot{M}}{0.5\\,M_\\odot\\,\\mathrm{s}^{-1}}\\right)^{1.29} \\left(\\frac{R_\\star}{30\\,\\mathrm{km}}\\right)^{n_{R_\\star}} \\left(\\frac{M_\\star}{1.4\\,M_\\odot}\\right)^{n_{M_\\star}} \\left(\\frac{\\langle\\epsilon_{\\nu_e}\\rangle}{12.6\\,\\mathrm{MeV}}\\right)^{n_\\epsilon} $10^{{51}}$\\,\\mathrm{erg\\,s}^{-1},$$ with $n_{R_\\star}=-2.58(\\dot{M}/0.5)^{-0.16}$, $n_{M_\\star}=1.97(\\dot{M}/0.5)^{-0.15}$, and $n_\\epsilon=-1.69(\\dot{M}/0.5)^{0.06}$ at $\\mathcal{M}=2.0$. The key physical claim is that the thermal content of the pre-shock flow, parameterized by the outer-boundary Mach number $\\mathcal{M}$, changes the normalization: lowering $\\mathcal{M}$ from 2.0 to 1.3 reduces $L_\\nu^{\\mathrm{crit}}$ by about 44 percent, while raising it from 2.0 to 3.0 changes it by about 12 percent, with convergence to pressureless free-fall at high $\\mathcal{M}$. The paper further argues that accretion of compositional interfaces such as the Si/O layer therefore promotes explosion two ways at once: a drop in $\\dot{M}$ and a rise in entropy that lowers the critical curve. When the antesonic, advection-to-heating timescale, and force-explosion conditions are compared, the antesonic ratio shows the least variation along the critical curve, staying near $\\max(c_s^2/v_{\\mathrm{esc}}^2)\\simeq 0.215\\pm0.01$.","pith_inferences":["Extension: if the roughly 44 percent normalization shift survives in 2D and 3D, then entropy jumps at compositional interfaces matter as much as $\\dot{M}$ drops for predicting which progenitors explode, and progenitor-to-explosion mapping should carry both.","Extension: a single pre-shock Mach number may be too coarse; a two-parameter family that varies entropy and infall speed separately would test whether the effect is truly pressurized inflow or a more general thermal-content dependence.","Extension: the near-constant antesonic ratio suggests a practical explosion probe for simulations and eventually observations: record $\\max(c_s^2/v_{\\mathrm{esc}}^2)$ behind the shock and flag values above roughly 0.215, with oscillations time-averaged."],"forward_implications":["At fixed accretion rate, pressurized pre-shock inflow lowers the neutrino luminosity required to explode, so models that assume pressureless free-fall overestimate how hard it is to explode realistic progenitors.","Accretion of an Si/O compositional interface pushes the shock toward explosion through two simultaneous effects: $\\dot{M}$ decreases and higher-entropy material lowers $L_\\nu^{\\mathrm{crit}}$ at fixed $\\dot{M}$.","Along the critical curve, the antesonic ratio stays near $0.215\\pm0.01$, varying only about 3-5 percent depending on $\\mathcal{M}$, making it the most stable of the tested explosion diagnostics.","Shock oscillations can temporarily push any diagnostic past its nominal critical value without producing an explosion, so time-averaged rather than instantaneous thresholds are needed near criticality.","The paper anticipates that the fitted power-law exponents may also depend on $R_\\star$, $M_\\star$, and $\\langle\\epsilon_\\nu\\rangle$ themselves, so the relation is a reference point rather than a universal scaling."],"supporting_citations":[{"why":"Establishes the critical-luminosity concept and the nonexistence of steady accretion solutions above a threshold.","marker":"Burrows & Goshy 1993"},{"why":"Supplies the antesonic condition and the pressureless-free-fall scaling relation the paper extends and tests.","marker":"Pejcha & Thompson 2012"},{"why":"Provides the advection and heating timescale criterion and the multi-dimensional normalization shift used for context.","marker":"Murphy & Burrows 2008"},{"why":"Defines the force explosion condition $\\tilde{\\Psi}$ that the paper evaluates along the critical curve.","marker":"Murphy & Dolence 2017"},{"why":"Provides the neutrino reaction rates and heating and cooling expressions used in the microphysics.","marker":"Scheck et al. 2006"},{"why":"Supplies the charged-current rates and the free-nucleon fraction $\\chi_N$ used to suppress pre-shock heating.","marker":"Qian & Woosley 1996"},{"why":"Provides the realistic progenitor entropy profiles used to connect pre-shock pressure to Si/O interface accretion.","marker":"Woosley & Heger 2007"},{"why":"Documents oscillatory but stable accretion models that the paper's stability classification builds on.","marker":"Gabay et al. 2015"}],"fun_headline_variants":["One scaling rule ties supernova explosion to neutrino luminosity","Pre-shock pressure shifts supernova explosion threshold by 44%","Antesonic ratio predicts supernova explosion onset","One formula captures supernova explosion threshold","Critical neutrino luminosity sets supernova explosions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that one outer-boundary Mach number captures the thermal state of the infalling material while the free-nucleon suppression factor completely stops pre-shock neutrino heating; for the lowest-Mach models that suppression fails, so spurious heating may contaminate the pre-shock entropy and partly produce the lowered critical luminosity the paper attributes to pressurized inflow.","fun_headline_variants_meta":{"raw":{"variants":["One scaling rule ties supernova explosion to neutrino luminosity","Pre-shock pressure shifts supernova explosion threshold by 44%","Antesonic ratio predicts supernova explosion onset","One formula captures supernova explosion threshold","Critical neutrino luminosity sets supernova explosions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000745,"raw_usage":{"total_tokens":3511,"prompt_tokens":1326,"completion_tokens":2185,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":942,"completion_tokens_details":{"reasoning_tokens":2113}},"tokens_in":942,"tokens_out":2185,"duration_ms":14302,"temperature":1.0,"reasoning_tokens":2113,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:48:01.392434+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the $\\mathcal{M}=1.3$ series again with pre-shock neutrino heating forcibly set to zero; if $L_\\nu^{\\mathrm{crit}}$ rises back toward the $\\mathcal{M}=2.0$ values, the low-Mach lowering is a numerical artifact rather than the pressurized-inflow effect.","supporting_citations":[{"cited_title":"W., Dolence J","cited_arxiv_id":null,"evidence_quote":"Defines the force explosion condition $\\tilde{\\Psi}$ that the paper evaluates along the critical curve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents oscillatory but stable accretion models that the paper's stability classification builds on."}],"review_version":1}