Pith. sign in

REVIEW 3 major objections 2 minor 40 references

Towards high-precision inspiral gravitational waveforms from binary neutron star mergers in numerical relativity

T0 review · 3 major / 2 minor · reviewed 2026-08-05 · deepseek-v4-flash

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

Pith's one-line read Upgrading the hydrodynamics Riemann solver to fourth order reduces the continuum-limit phase error of binary neutron star inspiral waveforms to about 0.27 radians at merger.

desk verdict The supplied full text is a different paper (ASAS-SN supernova rates), so the numerical-relativity claims in the abstract have no supporting evidence in the manuscript. read the letter →

arxiv 2508.10981 v2 pith:EGLOVPQB submitted 2025-08-14 astro-ph.HE gr-qc

Kenta Kiuchi This is my paper · ORCID
classification astro-ph.HEgr-qc
keywords gravitationalwavesbinaryneutronstarsnumericalrelativityHLLCRiemannsolverfourth-orderfinitevolumeinspiralwaveformaccuracyconvergenceorderadaptivemeshrefinement
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Numerical relativity simulations of binary neutron star mergers must supply very accurate waveforms for gravitational-wave detectors, and the hydrodynamics solver is often the dominant source of phase error. This paper reports that a fourth-order accurate finite-volume HLLC Riemann solver, implemented in the SACRA-MPI code, halves that error: the continuum-limit residual phase at merger drops from $0.58\pm0.22$ rad to $0.27\pm0.07$ rad out of a total phase of $\approx 176$ rad. The measured convergence order of the inspiral phase is about $2.1$ to $2.4$ for the fourth-order solver, compared with about $2.0$ for the second-order solver. If this is right, the hydrodynamics truncation error is no longer the limiting ingredient, and further waveform accuracy gains must come from other parts of the numerical setup.

What carries the argument

The central object is the fourth-order accurate finite-volume HLLC Riemann solver, an approximate Riemann solver that computes numerical fluxes through cell interfaces using a three-wave structure (two acoustic waves and a contact discontinuity), upgraded to fourth-order spatial accuracy. It carries the argument by removing the hydrodynamics truncation error that previously dominated the inspiral phase; the convergence-order measurements and the continuum extrapolation are both made possible by this solver upgrade.

What would settle it

Add a fifth and sixth finer resolution (for example about 60 m and 50 m) to the same binary configuration and recompute the residual phase error at merger; if the new residuals do not follow the same power-law trend and the extrapolated residual moves by more than the quoted $0.07$ rad, the claimed convergence order and continuum-limit error would be falsified. A complementary check is to change the mesh-refinement layout or the outer-boundary extraction radius and see whether the phase error shifts by a comparable amount.

Watch

Extended reading notes

Core claim

The authors claim that replacing the second-order finite-volume Riemann solver with a fourth-order solver in full dynamical-spacetime binary neutron star simulations improves the gravitational-wave phase accuracy in a quantified way. In a resolution study with grid spacings of about 78, 94, 118, and 135 m, the inspiral phase error converges with order $2.1\pm0.05$ to $2.4\pm0.27$ for the fourth-order solver, while the second-order solver stays near order $2.0$. A Richardson-style extrapolation to the continuum limit gives a residual phase error at merger of $0.27\pm0.07$ rad for the fourth-order solver and $0.58\pm0.22$ rad for the second-order solver, out of a total accumulated phase of $\a

Load-bearing premise

The four chosen grid spacings (about 78, 94, 118, and 135 m) all lie in the asymptotic convergent regime, so the phase error follows a single power law in grid spacing and the continuum extrapolation to $0.27$ rad is trustworthy.

Editorial extensions

If this is right

  • If the reported accuracy holds, numerical-relativity waveforms for binary neutron star inspirals can be used as templates with a phase error below a third of a radian at merger, easing the accuracy burden on the inspiral part of waveform models.
  • The measured convergence order of about 2.1 to 2.4 implies that the hydrodynamics solver is no longer the dominant source of phase error; further gains must come from metric evolution, mesh-refinement interface treatment, or gravitational-wave extraction.
  • The same fourth-order solver, validated on shock-tube and smooth-flow tests and applied through a short post-merger phase, is positioned for longer post-merger and remnant-disk simulations where small numerical errors accumulate over many dynamical times.
  • Preserving the $\pi$-symmetry without imposing it means simulations of non-spinning equal-mass binaries can be checked for symmetry-breaking artifacts, strengthening confidence in the waveform.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural next step, not reported in the paper, would be to add a fifth resolution level near 60 m and verify that the phase-error curve follows the same power law; this would directly test whether the quoted 0.27 rad residual is a true continuum limit or a fit artifact.
  • Because the phase error is no longer dominated by hydrodynamics, pairing this solver with a weak-form or spectral metric evolution may push inspiral accuracy below 0.1 rad, which would matter for next-generation detectors.
  • The fourth-order solver may also reduce spurious numerical angular momentum transport in the post-merger disk, changing predictions for ejecta masses and kilonova light curves; this is a testable extension outside the paper's stated scope.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 2 minor

Summary. The submitted manuscript, arXiv:2508.10981, has a title and abstract describing a numerical-relativity study: a newly implemented fourth-order HLLC Riemann solver in the SACRA-MPI code, 1D special-relativistic hydrodynamics tests, binary neutron star inspiral simulations, and a convergence study reporting inspiral phase-error orders of ~2.1–2.4 and continuum-limit residuals of 0.27±0.07 rad and 0.58±0.22 rad. However, the full text provided is an entirely different paper: Pessi et al., 'Supernova rates and luminosity functions from ASAS-SN II: 2014–2017 core-collapse supernovae and their subtypes' (arXiv:2508.10985v2). The body contains no equations or discussion of SACRA-MPI, HLLC solvers, finite-volume methods, gravitational waveforms, binary neutron star simulations, or any of the numerical results cited in the abstract. The claims made in the abstract are therefore not supported by any evidence in the submitted manuscript.

Significance. If the abstract's claims were substantiated, the work would be significant for numerical relativity and gravitational-wave template construction: a robust fourth-order finite-volume hydrodynamics solver in a production AMR code, with quantified continuum-limit phase errors, would be a useful development for high-precision inspiral waveforms from binary neutron star mergers. The reported factor-of-two reduction in phase error and the statement that hydrodynamics errors are no longer dominant would be valuable and would merit careful evaluation. However, the submitted manuscript contains none of the apparatus needed to assess these claims: no method description, no validation tests, no convergence tables, no figures, and no reproducibility artifacts. Significance cannot be granted to an abstract alone.

major comments (3)
  1. [Full text (entire manuscript)] The manuscript body is not the numerical-relativity paper described by the title and abstract. It is a supernova-rates paper by Pessi et al. (arXiv:2508.10985v2), with no mention of SACRA-MPI, HLLC, Riemann solvers, special-relativistic hydrodynamics tests, binary neutron star mergers, gravitational waveforms, or a resolution study. None of the central claims of the abstract are supported by any equation, figure, table, or description in the submitted text.
  2. [Abstract (headline numerical claims)] The quantitative results quoted in the abstract — e.g., convergence order ≈2.1±0.05–2.4±0.27, residual continuum phase errors of 0.27±0.07 rad and 0.58±0.22 rad over ≈176 rad, and the claimed difference between fourth-order and second-order solvers — are presented without definitions of the phase-error estimator, the Richardson-extrapolation procedure, the resolution grid, or the error bars. In the present manuscript these are not results, because no methodology is given. This is a missing-evidence failure, not a technical flaw in a described analysis.
  3. [References and bibliography] The reference list is entirely consistent with the ASAS-SN supernova-rate paper and contains no citations to numerical-relativity methods, Riemann solvers, SACRA-MPI, or binary neutron star waveform accuracy studies. This independently confirms that the submitted text is a different paper, and it leaves no route for a referee to trace or verify the claimed numerical methods.
minor comments (2)
  1. [Header/front matter] The PDF header displays 'arXiv:2508.10985v2' and the A&A manuscript number aa56799-25, along with a title, author list, and affiliations that all correspond to the ASAS-SN supernova paper, not to the submitted title arXiv:2508.10981.
  2. [General] If the intended submission is the numerical-relativity paper, the current file is the wrong manuscript. At minimum, the correct file must be submitted before any review can proceed.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found; manuscript body is a different paper than the abstract, but no derivation reduces to its inputs.

full rationale

The abstract accompanying arXiv:2508.10981 describes a convergence study of a fourth-order HLLC solver in SACRA-MPI, with convergence orders and residual phase errors measured via standard Richardson-style extrapolation across resolutions. No parameter is fitted to a target result, no equation is defined in terms of the quantity it purports to predict, and no self-citation chain is invoked; the reported convergence orders and phase errors are empirical measurements, not outputs of an optimization. The provided full text, however, is an entirely different manuscript (ASAS-SN supernova rates, arXiv:2508.10985), so the abstract's numerical-relativity derivation chain is absent from the submitted body. This is a severe missing-evidence or manuscript-mismatch problem, but it is not circularity: there is no quoted equation or fitting step that reduces to its own input. Similarly, the ASAS-SN paper, if taken as the manuscript, uses standard completeness corrections from injection-recovery simulations and empirical rate calculations, with no load-bearing self-citation or fitted-input-called-prediction pattern. Therefore, the appropriate finding is no significant circularity (score 0), while noting the manuscript/abstract inconsistency as a separate integrity concern.

Assumptions & free parameters 1 free parameters · 3 assumptions · 0 invented entities

Ledger for an abstract-only review. No physical free parameters are fitted: the only chosen numbers are the four grid spacings defining the convergence study, whose selection directly shapes the quoted orders and residual errors. The axioms are the modeling and analysis assumptions the abstract silently relies on: faithful GRHD implementation, transferability of 1D analytic tests to 3D performance, and single-power-law behavior across the resolution band. The residual 0.27/0.58 rad numbers are extrapolated estimates, not measured values at any single resolution, which is why the third axiom is flagged ad hoc. No invented physical entities are proposed.

free parameters (1)
  • Grid spacings of the resolution study (78, 94, 118, 135 m)
    Chosen by hand to define the convergence study. The quoted convergence orders and the extrapolated 0.27/0.58 rad residual phase errors are direct outputs of this choice; with only a 1.73:1 resolution span, the results are sensitive to it.
assumptions (3)
  • domain assumption The general-relativistic hydrodynamics and Einstein equations implemented in SACRA-MPI faithfully model binary neutron star mergers.
    Implicit in every claim in the abstract; the validation tests probe the solver, but the transfer of 1D test performance to 3D merger accuracy is assumed.
  • domain assumption The one-dimensional analytic solutions (simple wave, shock tube) exercise the same solver behavior relevant to the 3D inspiral problem.
    The abstract states validation in one-dimensional special relativistic hydrodynamics tests. The relevance of these tests to 3D BNS waveform accuracy is asserted, not demonstrated.
  • ad hoc to paper All four resolutions lie in the asymptotic convergent regime so the phase error obeys a single power law; the continuum-limit residual (0.27/0.58 rad) is then well-defined.
    The headline phase-error numbers depend on this unstated analysis assumption. From the abstract alone there is no evidence the 78-135 m range is in the clean power-law regime, e.g., no fifth resolution check is reported.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Towards high-precision inspiral gravitational waveforms from binary neutron star mergers in numerical relativity." pith.science (2026). https://pith.science/paper/EGLOVPQB

@misc{pith2026250810981,
  author       = {Pith},
  title        = {Pith review of: Towards high-precision inspiral gravitational waveforms from binary neutron star mergers in numerical relativity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EGLOVPQB}},
  note         = {Machine review of arXiv:2508.10981}
}
abstract

We report the performance of a newly implemented fourth-order accurate finite-volume HLLC Riemann solver in the adaptive-mesh-refinement numerical relativity code {\tt SACRA-MPI}. First, we validate our implementation in one-dimensional special relativistic hydrodynamics tests, i.e., a simple wave and shock tube test, which have analytic solutions. We demonstrate that the fourth-order convergence is achieved for the smooth flow, which cannot be achieved in our original second-order accurate finite-volume Riemann solver. We also show that our new solver is robust for the strong shock wave emergence problem. Second, we validate the implementation in a dynamical spacetime by demonstrating that {\tt SACRA-MPI} perfectly preserves the $\pi$-symmetry without imposing the $\pi$-symmetry in a short-term ($\sim 20~{\rm ms}$ in the inspiral and subsequent post-merger phase) non-spinning equal-mass binary neutron star merger simulations. Finally, we quantify the accuracy of $\approx 28$ cycles inspiral gravitational waveforms from binary neutron star mergers by conducting a resolution study with $\approx 78, 94$, $118$, and $135$ m. We find that the fourth-order accurate Riemann solver achieves the convergence order $\approx 2.1\pm{0.05}$--$2.4\pm{0.27}$, i.e., slightly evolving with time, in the inspiral gravitational wave phase, while the second-order accurate Riemann solver achieves the convergence order $\approx 2.0\pm{0.5}$. The residual phase error towards the continuum limit at the merger is $0.27\pm 0.07$ rad and $0.58\pm 0.22$ rad out of a total phase of $\approx 176$ rad, respectively, for the fourth- and second-order accurate Riemann solver.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 24 canonical work pages

  1. [1]

    D., Bahcall, J

    Arnett, W. D., Bahcall, J. N., Kirshner, R. P., & Woosley, S. E. 1989, ARA&A, 27, 629 Astropy Collaboration, Price-Whelan, A. M., Sip˝ocz, B. M., et al. 2018, AJ, 156, 123 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

    work page 1989

  2. [2]

    K., Driver, S

    Baldry, I. K., Driver, S. P., Loveday, J., et al. 2012, MNRAS, 421, 621

    work page 2012

  3. [3]

    2009, A&A, 499, 653

    Bazin, G., Palanque-Delabrouille, N., Rich, J., et al. 2009, A&A, 499, 653

    2009

  4. [4]

    C., Kulkarni, S

    Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002

    2019

  5. [5]

    Bessell, M. S. 1990, PASP, 102, 1181

    work page 1990

  6. [6]

    A., Brown, G

    Bethe, H. A., Brown, G. E., Applegate, J., & Lattimer, J. M. 1979, Nucl. Phys. A, 324, 487

    work page 1979

  7. [7]

    T., Riello, M., Cappellaro, E., et al

    Botticella, M. T., Riello, M., Cappellaro, E., et al. 2008, A&A, 479, 49

    work page 2008

  8. [8]

    S., Stanek, K

    Brown, J. S., Stanek, K. Z., Holoien, T. W. S., et al. 2019, MNRAS, 484, 3785

    work page 2019

Show all 40 references
  1. [9]

    Byrne, R. A. & Fraser, M. 2022, MNRAS, 514, 1188

    2022

  2. [10]

    T., Pignata, G., et al

    Cappellaro, E., Botticella, M. T., Pignata, G., et al. 2015, A&A, 584, A62

    2015

  3. [11]

    1999, A&A, 351, 459

    Cappellaro, E., Evans, R., & Turatto, M. 1999, A&A, 351, 459

    1999

  4. [12]

    2005, A&A, 430, 83

    Cappellaro, E., Riello, M., Altavilla, G., et al. 2005, A&A, 430, 83

    2005

  5. [13]

    Chevalier, R. A. & Irwin, C. M. 2011, ApJ, 729, L6

    2011

  6. [14]

    & Nelemans, G

    Chruslinska, M. & Nelemans, G. 2019, MNRAS, 488, 5300

    2019

  7. [15]

    & Hjorth, J

    Cold, C. & Hjorth, J. 2023, A&A, 670, A48

    2023

  8. [16]

    2012, Nature, 491, 228

    Cooke, J., Sullivan, M., Gal-Yam, A., et al. 2012, Nature, 491, 228

    2012

  9. [17]

    2014, A&A, 572, A23

    Corbelli, E., Thilker, D., Zibetti, S., Giovanardi, C., & Salucci, P. 2014, A&A, 572, A23

    2014

  10. [18]

    Coughlin, E. R. & Armitage, P. J. 2018, MNRAS, 474, 3857

    2018

  11. [19]

    G., et al

    Dahlen, T., Strolger, L.-G., Riess, A. G., et al. 2012, ApJ, 757, 70

    2012

  12. [20]

    K., Kasliwal, M

    Das, K. K., Kasliwal, M. M., Fremling, C., et al. 2025, PASP, 137, 044203

    2025

  13. [21]

    W., Taggart, K., Tinyanont, S., et al

    Davis, K. W., Taggart, K., Tinyanont, S., et al. 2023, MNRAS, 523, 2530 de Jaeger, T., Anderson, J. P., Galbany, L., et al. 2018, MNRAS, 476, 4592

    2023

  14. [22]

    D., Kochanek, C

    Desai, D. D., Kochanek, C. S., Shappee, B. J., et al. 2024, MNRAS, 530, 5016

    2024

  15. [23]

    J., Waldman, R., Livne, E., & Blondin, S

    Dessart, L., Hillier, D. J., Waldman, R., Livne, E., & Blondin, S. 2012, MNRAS, 426, L76

    2012

  16. [24]

    J., Prieto, J

    Dong, S., Shappee, B. J., Prieto, J. L., et al. 2016, Science, 351, 257

    2016

  17. [25]

    H., Matthews, K., Neugebauer, G., & Persson, S

    Elias, J. H., Matthews, K., Neugebauer, G., & Persson, S. E. 1985, ApJ, 296, 379

    1985

  18. [26]

    Filippenko, A. V . 1997, ARA&A, 35, 309

    1997

  19. [27]

    V ., Li, W

    Filippenko, A. V ., Li, W. D., Treffers, R. R., & Modjaz, M. 2001, in Astronom- ical Society of the Pacific Conference Series, V ol. 246, IAU Colloq. 183: Small Telescope Astronomy on Global Scales, ed. B. Paczynski, W.-P. Chen, & C. Lemme, 121

    2001

  20. [28]

    D., Khandrika, H., Rubin, D., et al

    Fox, O. D., Khandrika, H., Rubin, D., et al. 2021, MNRAS, 506, 4199

    2021

  21. [29]

    D., Brennan, S

    Fraser, M., Stritzinger, M. D., Brennan, S. J., et al. 2021, arXiv e-prints, arXiv:2108.07278

    2021 arXiv

  22. [30]

    A., Sharma, Y ., et al

    Fremling, C., Miller, A. A., Sharma, Y ., et al. 2020, ApJ, 895, 32

    2020

  23. [31]

    R., Vincenzi, M., et al

    Frohmaier, C., Angus, C. R., Vincenzi, M., et al. 2021, MNRAS, 500, 5142

    2021

  24. [32]

    2012, Science, 337, 927

    Gal-Yam, A. 2012, Science, 337, 927

    2012

  25. [33]

    2017, in Handbook of Supernovae, ed

    Gal-Yam, A. 2017, in Handbook of Supernovae, ed. A. W. Alsabti & P. Murdin, 195

    2017

  26. [34]

    2022, Nature, 601, 201

    Gal-Yam, A., Bruch, R., Schulze, S., et al. 2022, Nature, 601, 201

    2022

  27. [35]

    J., Wetzel, A., Hopkins, P

    Gandhi, P. J., Wetzel, A., Hopkins, P. F., et al. 2022, MNRAS, 516, 1941

    2022

  28. [36]

    2012, A&A, 542, A29

    Georgy, C., Ekström, S., Meynet, G., et al. 2012, A&A, 542, A29

    2012

  29. [37]

    Z., Kochanek, C

    Godoy-Rivera, D., Stanek, K. Z., Kochanek, C. S., et al. 2017, MNRAS, 466, 1428 González-Gaitán, S., Gutiérrez, C. P., Martins, G., et al. 2025, A&A, 700, A119

    2017

  30. [38]

    B., Huang, S., et al

    Graur, O., Bianco, F. B., Huang, S., et al. 2017, ApJ, 837, 120

    2017

  31. [39]

    B., & Modjaz, M

    Graur, O., Bianco, F. B., & Modjaz, M. 2015, MNRAS, 450, 905

    2015

  32. [40]

    J., Hey, D., et al

    Hart, K., Shappee, B. J., Hey, D., et al. 2023, arXiv e-prints, arXiv:2304.03791

    2023 arXiv

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.