{"id":"01cece02-34a0-46b6-815e-51a7a81bb4d2","arxiv_id":"2602.02507","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"A validated GPU-accelerated astrophysical simulation code now includes Newtonian self-gravity and a tabulated nuclear EOS, and delivers an eleven-progenitor adiabatic core-collapse benchmark whose shock speed anti-correlates with compactness.","lead":"GenASiS, a code for simulating supernova explosions, is extended with validated Newtonian self-gravity and GPU-accelerated hydrodynamics using a nuclear equation of state, checked against known analytic collapse solutions. The authors run an adiabatic collapse-bounce-explosion benchmark across eleven pre-supernova models and propose it to the community as a standard code-comparison test.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline anti-correlation rests on two hand-chosen measurement definitions — compactness at m=2.5 M_sun and a two-point shock-speed estimate truncated at the outer grid boundary — with no sensitivity tests reported.","rationale":"The reader identified the same weakest assumption I do: the headline physical result depends on two hand-chosen measurements without sensitivity tests. This is genuinely load-bearing because the abstract presents the inverse correlation as a finding, and the proposed benchmark's value depends on that correlation being a robust property of the adiabatic collapse dynamics, not an artifact of how shock speed and compactness are quantified. I independently weighed other potential concerns — e.g., the Poisson solver's multipole truncation at L=12, the missing gravitational mass beyond the grid, the phase-transition entropy artifact — but none undercuts the central claim as directly. The L=12 truncation is plausibly controlled by the convergence tests and the quasi-spherical nature of collapse; missing mass outside r_out contributes only a constant to the interior potential under spherical symmetry and thus does not affect the shock dynamics inside the grid; the phase-transition artifact is acknowledged and the kinetic energy appears well converged. The reader's conditions (1)–(3) are appropriate and consistent with my read. I therefore agree with the CONDITIONAL verdict and see no reason to move it.","tokens_in":21823,"tokens_out":9898,"duration_ms":101482,"concrete_test":"Using the public dataset (DOI 10.13139/OLCF/3011754), recompute for each of the 11 models: (a) the average shock speed with a fixed, reproducible shock detection criterion (e.g., maximum of the entropy gradient) over a fixed post-bounce window (e.g., 50–300 ms), and (b) compactness ξ_m for m = 1.5, 2.0, 3.0, 3.5 M_sun from the progenitor models. Then compute Spearman rank correlations between shock speed / kinetic energy and collapse mass / compactness. If the anti-correlation remains significant (|ρ| ≳ 0.8, p < 0.05) across all choices, the claim is robust; otherwise the abstract's statement is definition-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The only new physical claim in the abstract — that average shock expansion speed and total kinetic energy are inversely correlated with progenitor mass at collapse and with compactness — is not load-bearing for the solver validation, but it is the headline result. Its robustness rests on two definitional choices made in §3.5. First, the average shock speed is computed from exactly two points: the position/time when the shock is 'first detected' and the position/time 'just before the shock passes the grid boundary.' The detection criterion is unspecified, and the outer boundary at r_out ≈ 1.09×10^4 km (in the oxygen shell) truncates the measured trajectory, so the average may omit late-time deceleration that could differ across progenitors. Second, compactness (Eq. 67) is evaluated at a single mass scale, m = 2.5 M_sun, following O'Connor & Ott (2011), but no sensitivity to m is shown. If the 'especially tight' anti-correlation with compactness (Fig. 20) shifts under a different shock-detection definition or a different m, the proposed benchmark's baseline expectation would be overstated. The solver validation in §3.1–3.4 (spheroid potential, dust collapse, polytropic collapse) is credible and not the source of concern; the concern is specifically the generality of the headline trend.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This methods paper (GenASiS II) adds Newtonian self-gravity and updated finite-volume hydrodynamics to the GenASiS framework. The Poisson solver uses a multipole expansion with GPU offloading, and the fluid solver uses parabolic reconstruction, tabulated baryonic equations of state, and angular coarsening to handle spherical-coordinate singularities. The new solvers are validated against four external test problems: the analytic potential of homogeneous spheroids, the collapse of a homogeneous dust sphere, the collapse of a homogeneous dust spheroid, and the self-similar polytropic collapse solution of Yahil and collaborators. The paper then presents adiabatic collapse, bounce, and explosion simulations of 11 pre-supernova progenitors, proposes this problem as a community benchmark, and reports that the average shock speed and total kinetic energy are inversely correlated with the progenitor mass at collapse and with the compactness parameter.","tokens_in":22140,"tokens_out":13342,"duration_ms":123067,"significance":"The solver-validation component is strong and is the main scientific contribution: the four external benchmarks are genuine, with independently derived reference solutions, and the reported convergence behavior (Figures 5, 7, 9, 12) is consistent with the expected orders. The public dataset and released code are valuable assets, and the proposed adiabatic benchmark is well timed given known code-to-code differences in core-collapse supernova modeling. The headline physical trend in §3.5 is plausible and physically motivated, but its robustness is not yet established because it rests on several definitional choices that are not tested. If those choices are demonstrated to be inconsequential, or if the claims are suitably qualified, the paper would be a solid contribution to the CCSN methods literature.","major_comments":[{"comment":"The headline physical claim—average shock speed and total kinetic energy inversely correlated with collapse mass and compactness—rests on two hand-chosen definitions. The average shock speed is a two-point estimate: 'first detected' and 'just before the shock passes the grid boundary' at r_out≈1.09×10^4 km, but the detection criterion is not specified and the trajectory is truncated. Compactness is evaluated only at m=2.5 M_sun (Eq. 67), and the time at which the total kinetic energy is evaluated is not stated. No sensitivity tests are reported for these choices. Since this trend is proposed as a benchmark baseline expectation, please specify the shock-detection algorithm, test alternative detection criteria and outer-boundary placements, show compactness at several m values (e.g., 1.5, 2.0, 3.0 M_sun), and report correlation coefficients. Without these, the 'especially tight' anti-corre","section":"§3.5, Eq. (67), Figures 19–20"},{"comment":"Near the coordinate singularities the code coarsens angular blocks and, in the collapse runs, zeroes lateral momenta in the first polar cell and within 0.1 r_core. These operations can suppress non-radial degrees of freedom by construction, yet the paper claims the 3D runs remain spherically symmetric and are 'visually indistinguishable' from 1D/2D. A quantitative comparison (e.g., L1 or L∞ differences between 1D, 2D, and 3D radial profiles) and a test with the lateral-momentum zeroing disabled or with varied coarsening thresholds should be reported to show that the benchmark result is not an artifact of this regularization.","section":"§2.3 and §3.5"}],"minor_comments":[{"comment":"The summation over r∈{R,I} is easy to misread as requiring both regular and irregular terms at every radius; in fact the two terms are mutually exclusive Green's-function pieces. Please rewrite the expression with interior and exterior contributions clearly labeled.","section":"§2.2, Eq. (21)"},{"comment":"The symbol A^a_lm is used both for angular kernel functions and for angular moments of the source (Eq. 14). Different notation for the moments (e.g., calligraphic A) would avoid confusion.","section":"§2.2, Eqs. (8)–(16)"},{"comment":"Please state how D(0)=1.75 is selected from Yahil (1983) and how the error in the numerically integrated ODE reference solution is controlled; as written, the reference solution is not fully reproducible from the text.","section":"§3.4"},{"comment":"The abstract's 'multidimensional computations' covers 2D for all eleven progenitors but only S12, S25, and S40 in 3D. Please state this explicitly in the abstract or conclusions to avoid overgeneralization.","section":"§3.5 and Abstract"},{"comment":"The evaluation time and definition of 'total kinetic energy' used in Figure 19 should be stated in the caption or text; currently the reader must infer it from context.","section":"§3.5, Figure 19"},{"comment":"Typo: 'neighbhors' should be 'neighbors'.","section":"§2.3"}],"recommendation":"major_revision","confidential_remarks":"The solver validation is credible and the dataset release is a genuine strength. My recommendation is driven by the robustness of the §3.5 benchmark trend, not by the validation sections. If the authors add the requested sensitivity tests and qualify the abstract's physical claim accordingly, I would support publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read this as a straightforward code paper, and the core validation holds up. The Poisson solver and finite-volume updates are tested against real external benchmarks — Binney & Tremaine spheroids, Lin et al. dust collapse, Yahil & Lattimer polytropic collapse — with convergence behavior that looks right (first-order for the discontinuous dust problems, near-second-order for the smooth polytrope). The GPU offload details are concrete, including the OpenMP target reductions. The modification to Skinner-style parabolic reconstruction is clearly documented and motivated by an actual instability, which is the kind of thing that needs to be in the literature.\n\nThe most useful contribution is the proposed adiabatic benchmark: eleven Woosley & Heger progenitors, a tabulated EOS, and a public dataset DOI. That gives the CCSN community a cheap, shared test problem for the code-to-code differences Janka keeps pointing at. The paper is also honest about its own warts — the phase-transition entropy bump near 10 km, the resolution dependence of central temperature and entropy, the absence of neutrinos. That honesty is real credit.\n\nThe stress-test concern lands, though not as a fatal strike. The abstract's anti-correlation claim rests on two specific measurement choices: compactness at m = 2.5 M_sun and a two-point shock speed estimate truncated at the outer boundary. No sensitivity study is reported, and the shock detection criterion is not defined. So the 'especially tight' compactness correlation in Figure 20 is, right now, an observation plus a plausible physical story, not a robust benchmark result. I would want the authors to add correlation statistics and test m = 1.5 and 3.0 M_sun, plus a different definition of average shock speed (e.g. shock radius at a fixed post-bounce time). If those change the ranking, the benchmark still works, but the stated baseline expectation would need rewording.\n\nMinor things: the dataset DOI should be verified as actually public and figure-reproducing before I'd rely on it, and the ~15x GPU speedup is mentioned without much detail — fine for a methods paper, but I'd like a footnote on how the proportional resource test was done. The phase-transition artifact's effect on explosion energies is acknowledged but not quantified; a brief estimate would help.\n\nOverall: the solver validation is the load-bearing part and it passes. The benchmark proposal is worth taking seriously. This should go to peer review — ApJS is the right venue — with the requested sensitivity analysis and dataset verification as conditions, not as grounds for rejection.","headline":"Solid methods paper with credible solver validation; the headline anti-correlation is under-quantified but not fatal.","tokens_in":22720,"tokens_out":2080,"would_cite":true,"duration_ms":23015,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["85A15","85-08","65M08"],"pacs":["95.30.Lz","97.60.Bw"],"model":"deepseek-v4-flash","headline":"A new GPU-accelerated Poisson and finite-volume solver passes known Newtonian collapse tests, and adiabatic core-collapse supernova models show shock speed and kinetic energy anti-correlated with progenitor compactness.","keywords":["core-collapse supernovae","Newtonian self-gravity","Poisson solver","finite-volume methods","GPU acceleration","tabular equations of state","compactness parameter","adiabatic collapse benchmark"],"falsifier":"Run the adiabatic explosion for the same eleven progenitors with the outer boundary extended well beyond 10^4 km (or record the full shock-radius time series) and recompute average shock speed; similarly vary the compactness mass cut in the compactness definition over a range such as 1.5–3.5 solar masses. If the anti-correlation with compactness becomes noisy or reverses, the claimed trend depends on the measurement window and the single mass scale.","tokens_in":21612,"feed_emoji":"💥","tokens_out":10114,"duration_ms":89057,"temperature":0.7,"pith_summary":"This paper introduces the Newtonian self-gravity and updated finite-volume capabilities of GenASiS, a simulation system aimed at core-collapse supernovae. The authors test a multipole Poisson solver and a shock-capturing fluid solver on homogeneous spheroids, on dust collapse, and on self-similar polytropic collapse, all of which match known solutions. They then simulate the adiabatic collapse, bounce, and prompt explosion of eleven pre-supernova progenitors with a tabulated equation of state, finding that the average shock speed and total kinetic energy of the explosion are inversely correlated with the mass at the onset of collapse and with the compactness parameter—most tightly with compactness. Because the same solvers run efficiently on GPUs, the paper proposes this adiabatic benchmark as a reproducible standard for comparing core-collapse supernova codes before neutrino transport is added.","feed_headline":"Shock speed in supernova models tracks compactness, not birth mass","feed_subtitle":"A validated GPU Poisson solver plus adiabatic supernova runs show explosion speed and energy track core compactness.","key_machinery":"The central machinery is a multipole Poisson solver: it expands the Green's function 1/|r−r′| in real spherical harmonics, computes angular moments of the density in each radial shell, and performs outward and inward radial integrations to obtain regular and irregular moments—so the cost scales with the number of radial shells times a modest number of multipole terms. The fluid dynamics rests on a finite-volume method with parabolic reconstruction constrained to be monotone, explicit time stepping, a tabulated baryonic equation of state, and a spherical-coordinate mesh that coarsens angular cell blocks near the origin and polar axis so that the explicit time-step limit is set by the radial c","core_discovery":"On its own terms, the paper claims that GenASiS now handles Newtonian self-gravity and baryonic fluid dynamics reliably: a multipole Poisson solver converges to the analytic potential of homogeneous spheroids, and the fluid solver reproduces the collapse of a dust sphere, a dust spheroid, and a self-similar polytropic sphere. In the adiabatic collapse, bounce, and prompt explosion of eleven pre-supernova progenitors, the paper finds that the explosions are prompt, remain spherically symmetric in 1D, 2D, and 3D, and that the average shock expansion speed and total kinetic energy are inversely correlated with the mass at the onset of collapse and with the compactness parameter evaluated at 2.5","pith_inferences":["If the compactness–vigour anti-correlation survives the addition of neutrinos and multidimensional stochasticity, a single pre-collapse measurement (the 2.5-solar-mass compactness) could serve as a quick predictor of explosion kinetic energy in supernova models.","Because the average shock speed is a two-point estimate truncated by the outer boundary at about 10^4 km, the headline correlation may depend on boundary placement and shock-detection criteria; moving the boundary outward or using a shock-radius time series would test this.","The coarsening-block technique for spherical-coordinate singularities could transfer to other central or polar problems—self-gravitating cloud collapse, accretion disks, or any spherical-grid simulation—regardless of the application domain.","The released dataset would allow re-analysis with different compactness mass cuts (say 1.5–3.5 solar masses) or integral density measures, which might reveal whether the single mass scale of 2.5 solar masses is the optimal predictor."],"forward_implications":["If the validation holds, GenASiS becomes a production-capable Newtonian self-gravity code for problems with tabulated equations of state, including relativistic approximations that fit into a Poisson-like framework.","The proposed adiabatic collapse-bounce-explosion benchmark could become a community standard: persistent code-to-code differences on this simplified problem could be traced to mesh, reconstruction, gravity solver, or equation-of-state handling before neutrino transport is added.","The tight inverse correlation between shock speed and compactness gives a quantitative baseline: in adiabatic models, more compact cores explode more slowly and with less kinetic energy, informing expectations for 'explodability' studies with more complete physics.","The reported speedup of roughly 15x for GPU over CPU-only runs indicates the solvers are practical on current accelerator-based supercomputers, making the benchmark reproducible on such machines."],"fun_headline_variants":["Supernova shock speed tied to compactness, not mass","Compactness, not mass, sets shock speed in supernova runs","GenASiS shows prompt explosions track core compactness","Simulation: shock speed inversely linked to compactness"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The headline result rests on two hand-picked measurement definitions—the average shock speed from first detection to just before the outer boundary at about 10^4 km, and the compactness at a single mass cut of 2.5 solar masses—rather than on the solver validation itself; if these definitions are not representative, the inverse correlation could be an artifact of how the explosion and the progenitor structure were measured.","fun_headline_variants_meta":{"raw":{"variants":["Supernova shock speed tied to compactness, not mass","Compactness, not mass, sets shock speed in supernova runs","GenASiS shows prompt explosions track core compactness","Simulation: shock speed inversely linked to compactness"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000146,"raw_usage":{"total_tokens":1007,"prompt_tokens":723,"completion_tokens":284,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":467,"completion_tokens_details":{"reasoning_tokens":216}},"tokens_in":467,"tokens_out":284,"duration_ms":5589,"temperature":1.0,"reasoning_tokens":216,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T08:40:21.942272+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the adiabatic explosion for the same eleven progenitors with the outer boundary extended well beyond 10^4 km (or record the full shock-radius time series) and recompute average shock speed; similarly vary the compactness mass cut in the compactness definition over a range such as 1.5–3.5 solar masses. If the anti-correlation with compactness becomes noisy or reverses, the claimed trend depends on the measurement window and the single mass scale.","supporting_citations":[],"review_version":1}