Toward More Realistic Machine-Learning Inference of the Dense-Matter Equation of State from Supernova Gravitational Waves

Aknur Sakan; Alejandro Torres-Forn\'e; Alisher Zhunuskanov; Almat Akhmetali; Ernazar Abdikamalov; Jos\'e Antonio Font; Marat Zaidyn; Nurzhan Ussipov; Y. Sultan Abylkairov

REVIEW 2 major objections 4 minor 3 cited by

Real noise, varied progenitors, and bounce-time uncertainty do not break machine-learning classification of the dense-matter equation of state from supernova gravitational waves.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-13 16:48 UTC pith:T6R67TSL

load-bearing objection Solid incremental robustness check: real O4a noise, four progenitors, and 20 ms bounce jitter do not kill discrete EOS classification; dataset size helps more than progenitor diversity hurts. the 2 major comments →

arxiv 2603.27680 v2 pith:T6R67TSL submitted 2026-03-29 astro-ph.HE

Toward More Realistic Machine-Learning Inference of the Dense-Matter Equation of State from Supernova Gravitational Waves

Almat Akhmetali , Y. Sultan Abylkairov , Marat Zaidyn , Aknur Sakan , Alisher Zhunuskanov , Nurzhan Ussipov , Jos\'e Antonio Font , Alejandro Torres-Forn\'e

show 1 more author

Ernazar Abdikamalov

This is my paper

classification astro-ph.HE

keywords core-collapse supernovaegravitational wavesnuclear equation of statemachine learningsupport vector machineLIGO noiseprogenitor diversitybounce-time uncertainty

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Gravitational waves emitted at the bounce of a rotating stellar core carry information about the nuclear equation of state (EOS). Earlier machine-learning studies showed that EOS models can be classified from those signals, but only under idealized conditions: simulated detector noise, a single progenitor star, and a perfectly known bounce time. This paper removes those three idealizations. Signals are injected into real LIGO Hanford O4a noise, the training set is expanded to four progenitors (12–40 solar masses) with many rotation rates, and the bounce epoch is allowed to wander by as much as 20 ms. None of the three changes substantially lowers classification accuracy. The larger multi-progenitor data set actually raises it, reaching roughly 92 percent at high signal-to-noise ratio when the frequency-domain representation is used. The work therefore argues that EOS inference from core-collapse gravitational waves remains feasible under more realistic observing conditions.

Core claim

When gravitational-wave signals from rotating core bounce are classified by a linear support-vector machine, neither real O4a detector noise, nor progenitor diversity across 12–40 solar masses, nor bounce-time uncertainty of up to 20 ms significantly degrades accuracy relative to the idealized single-progenitor, simulated-noise case. The expanded data set instead improves training, so that time-domain accuracy at SNR = 200 rises from about 84 percent (single progenitor) to 91.6 ± 3.1 percent (all four progenitors).

What carries the argument

A linear-kernel support-vector machine trained on 30 ms windows of whitened, band-pass-filtered strain, using either the raw time series or its Fourier amplitude spectrum after the simulated waveform is injected into real or colored Gaussian noise at a controlled signal-to-noise ratio.

Load-bearing premise

The short axisymmetric waveforms around bounce (roughly -2 to +6 ms) already contain the EOS-dependent features the classifier needs, so that later multi-dimensional effects such as prompt convection can be neglected or only roughly included.

What would settle it

Repeat the identical classification pipeline on a large set of high-resolution three-dimensional simulations that include realistic prompt convection and anisotropic neutrino emission; if accuracy collapses once those later components enter the analysis window, the claimed robustness fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

Solid incremental robustness check: real O4a noise, four progenitors, and 20 ms bounce jitter do not kill discrete EOS classification; dataset size helps more than progenitor diversity hurts.

read the letter

This is a clean, controlled extension of the authors’ own earlier ML-EOS papers. The new result is empirical and useful: injecting into real O4a Hanford noise, expanding from one to four progenitors (12–40 M⊙) with many T/|W| values, and allowing up to 20 ms bounce-time uncertainty leaves classification accuracy essentially intact. Frequency-domain SVM stays high (~85 % at SNR=200 even with 20 ms jitter); time-domain collapses under the same jitter, as expected. Larger multi-progenitor training sets actually improve accuracy (time-domain from ~84 % single-progenitor to 91.6 %±3.1 % with all four), and the balanced-dataset control shows that the gain is mostly from sample size rather than from any magical invariance to progenitor structure.

They do the experiments properly: 50 random train–test splits with error bars, real vs simulated noise side-by-side, time vs frequency, balanced vs unbalanced, and an appendix that re-runs on Richers et al. waveforms that include prompt convection. The convection test drops accuracy only ~11 % and the restricted-window check confirms the classifier is still mostly using the bounce/ring-down segment. Math and citation pattern look solid; the self-citations are the natural baseline.

Soft spots are real but already flagged by the authors and do not undermine the reported numbers. Everything is still axisymmetric CoCoNuT, truncated to −2 to +6 ms, four discrete EOSs, and no public code or waveforms. Continuous EOS regression and full 3-D turbulence remain future work. Those are limitations of scope, not hidden flaws in the robustness claim itself.

This is for people already working on CCSN GW parameter estimation who need to know whether the earlier idealizations were load-bearing. It is not a new method or a nuclear-physics breakthrough. I would send it to peer review; a serious referee can push on the remaining idealizations without desk-rejecting a carefully executed incremental study. Worth reading if you care about the practical pipeline; skip if you only want first-principles advances.

Referee Report

2 major / 4 minor

Summary. The manuscript tests whether three realism upgrades degrade machine-learning classification of four nuclear EOSs (SFHo, LS220, HSDD2, GShenFSU2.1) from rotating core-collapse supernova gravitational-wave signals. Using axisymmetric CoCoNuT waveforms for four progenitors (12–40 M⊙) spanning T/|W| ≈ 0.02–0.18 (886 total waveforms), the authors inject signals into real O4a LIGO Hanford noise and into simulated Gaussian noise colored by the detector PSD, allow bounce-time uncertainty Δtb up to 20 ms inside 30 ms analysis windows, and train a linear SVM (C=10) on time-domain and frequency-domain representations. Accuracy is reported as mean ± 1σ over 50 random train–test splits. Real and simulated noise give statistically consistent performance (e.g., ~91–92 % time-domain accuracy at SNR=200). Expanding from one to four progenitors raises accuracy when the dataset grows (Table I) but modestly lowers it when the training set is size-balanced (Table II). Frequency-domain features remain robust under bounce-time shifts while time-domain accuracy collapses; enlarging the frequency-domain training set by denser time shifts recovers ~91.6 % accuracy even at Δtb=20 ms. An appendix using longer Richers et al. waveforms that include prompt convection shows only a modest accuracy drop and confirms that the classifier still relies primarily on the bounce/ring-down segment.

Significance. If the reported robustness holds under the stated assumptions, the work supplies a concrete, controlled demonstration that real detector noise, multi-progenitor diversity, and realistic bounce-time uncertainty need not destroy discrete EOS classification for rapidly rotating CCSNe. The systematic real-versus-simulated noise comparison, the balanced-dataset control that isolates size from diversity, and the frequency-domain resilience under time shifts are useful methodological results for the growing literature on ML inference of nuclear physics from supernova GWs. The study is explicitly incremental—a step toward more realistic conditions—and correctly flags remaining limitations (axisymmetry, discrete EOS set, optimal orientation). These controlled experiments therefore constitute a solid intermediate contribution rather than a definitive inference pipeline.

major comments (2)

Section II.A and Appendix A: the central claim of robustness rests on the assertion that the −2 to +6 ms bounce/ring-down window (or its approximate extension with unphysical 2-D convection) captures the EOS-dependent features the classifier learns. While the appendix test with Richers et al. waveforms is a useful control and shows only an ~11 % accuracy drop, the simulations remain axisymmetric and the convection is acknowledged to be poorly resolved. A clearer quantitative statement of how much the reported accuracies could change under high-resolution 3-D bounce-plus-convection signals is needed before the robustness claim can be taken as fully load-bearing for realistic Galactic events.
Section III and Tables I–II: classification is performed on a fixed discrete set of four EOSs. The conclusion (Section IV) correctly notes that a continuous, parameterized EOS family and regression would be preferable, yet the abstract and main results still present high classification accuracies as evidence that EOS inference is robust. The manuscript should more explicitly bound the claim: the numbers demonstrate robustness of four-way classification under the tested systematics, not yet continuous EOS parameter recovery.

minor comments (4)

Figure 3 caption: the quoted accuracies for real versus simulated noise appear swapped relative to the body text of Section III.A; please align caption and text.
Section II.B: the precise definition of the 30 ms window placement relative to the random bounce-time offset Δtb could be stated more formally (e.g., injection start uniform on [0, Δtb]) so that the experiment is fully reproducible.
Throughout: a few typographical slips remain (“asses”, “constrainghts”, “intoroduction” in the Tukey reference). A light copy-edit pass would clean them.
Section II.C: the choice of linear SVM with C=10 is justified by prior work, but a one-sentence reminder of why a linear kernel remains adequate once real noise and multi-progenitor diversity are added would help readers who have not read the earlier papers.

Circularity Check

1 steps flagged

Empirical ML robustness study; prior self-citations supply baselines and method choice but do not force the reported accuracies by construction.

specific steps

self citation load bearing [Introduction / Section II.C (method choice)]
"In previous studies [75–77], we demonstrated the feasibility of using machine learning (ML) to infer the nuclear EOS from gravitational waves produced during rotating core bounce. ... To classify the EOS models using GW signals, we employ a support vector machine (SVM). This supervised ML algorithm demonstrated the best performance for similar tasks among the tested classical and deep learning models conducted by Abylkairov et al. [77]."

The choice of SVM and the claim of prior promising accuracy rest on the authors' own earlier numerical experiments. This is ordinary self-citation of independent simulation results, not a definitional loop that forces the new multi-progenitor / real-noise accuracies. It is therefore only a minor, non-load-bearing circularity.

full rationale

The paper's central claims are measured classification accuracies under controlled relaxations of three assumptions (real O4a noise, multi-progenitor diversity, bounce-time uncertainty). These are obtained by training/evaluating an SVM on an expanded library of CoCoNuT waveforms injected into detector noise (Sections II–III, Tables I–II, Figures 3–4). No quantity is defined in terms of the target accuracy, no free parameter is fitted to the same data that is then called a prediction, and no uniqueness theorem is invoked. Self-citations to the authors' earlier works [75–77] establish the prior single-progenitor baseline and the choice of linear SVM; those earlier results are independent numerical experiments, not algebraic identities that make the present numbers inevitable. The balanced-dataset control (Table II) further isolates dataset-size effects from progenitor diversity, confirming that the improvement is empirical rather than definitional. Appendix A tests the effect of approximate prompt convection and finds only modest degradation, again an empirical measurement. Consequently the derivation chain is self-contained against its own simulation library; circularity is limited to ordinary, non-load-bearing self-citation of prior numerical work.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The claim rests on standard numerical-relativity and ML practices plus a handful of domain modeling choices that are not independently verified inside the paper. No new physical entities are postulated; free parameters are ordinary hyper-parameters and analysis-window choices.

free parameters (3)

SVM regularization C = 10
Fixed at C=10 after grid-search cross-validation; the reported accuracies depend on this choice.
analysis window length and placement = 30 ms / Δtb ≤ 20 ms
30 ms windows with signal injected randomly up to Δtb; the 8 ms signal duration (−2 to +6 ms) is a modeling cut.
Tukey taper α and filter bands = α=0.1
α=0.1, 20–2000 Hz bandpass, notches at 60/120/240 Hz; standard but still free analysis choices.

axioms (4)

domain assumption Axisymmetric CoCoNuT simulations with Ye(ρ) deleptonization and leakage/heating adequately capture the EOS-dependent bounce and early ring-down GW signal.
Stated in Section II.A; justified by prior literature but not re-validated in 3-D here.
domain assumption The four chosen EOS models (SFHo, LS220, HSDD2, GShenFSU2.1) are sufficiently representative that classification accuracy among them generalizes to the broader EOS landscape.
Section II.A; the paper notes that a continuous parameterized family would be preferable.
domain assumption Optimal source and detector orientations can be assumed without loss of the robustness conclusion.
Explicitly listed as a remaining limitation in Section IV.
domain assumption Linear-kernel SVM is an adequate classifier for the task; more complex architectures would not reverse the robustness findings.
Chosen because it performed best in the authors’ prior comparison (Abylkairov et al. 2024).

pith-pipeline@v1.1.0-grok45 · 24329 in / 2809 out tokens · 33815 ms · 2026-07-13T16:48:20.484916+00:00 · methodology

0 comments

read the original abstract

Gravitational waves from core-collapse supernovae offer a unique probe of the equation of state (EOS) of dense nuclear matter. For rapidly rotating stars, previous machine-learning studies demonstrated promising EOS classification accuracy. However, these analyses relied on several simplifying assumptions. In this work, we relax three key assumptions. First, we include real detector noise. Second, we expand the analysis from a single progenitor model to four models spanning 12 to 40 solar masses, and for each mass we consider multiple rotational configurations, from slow to rapid. Third, we introduce uncertainty in the core bounce time of up to 20 ms, rather than assuming it is known precisely. We find that none of these effects significantly degrades EOS classification performance. Instead, the larger dataset associated with multiple progenitor models and noise realizations improves training and classification accuracy. This study is a step in a broader effort to progressively incorporate more realistic conditions into gravitational-wave inference for core-collapse supernovae.

discussion (0)

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Contrastive self-supervised convolutional autoencoder for core-collapse supernova gravitational-wave detection
gr-qc 2026-05 unverdicted novelty 7.0

A contrastive self-supervised convolutional autoencoder detects core-collapse supernova gravitational waves with performance comparable to supervised CNNs, better generalization to unseen waveforms, and ~120 kpc sensi...
The Generalization Gap in Machine Learning EoS Inference from Core-Collapse Supernova Gravitational Waves
astro-ph.HE 2026-07 conditional novelty 6.5

LightGBM and other regressors achieve R^{2}≈0.6–0.7 under random CV on CCSN GW catalogues but collapse to worse-than-mean performance under Leave-One-EoS-Out validation, exposing a generalisation gap for unseen EoS families.
Parameter Estimation Horizon of Core-Collapse Supernovae with Current and Next-Generation Gravitational-Wave Detectors
astro-ph.HE 2026-05 unverdicted novelty 5.0

Machine learning extracts core rotation and signal properties from CCSN gravitational waves, with next-generation detectors constraining rotation beyond 100 kpc for favorable orientations despite some uncertainties.