pith. sign in

arxiv: 2605.30430 · v1 · pith:AHQ4TJE5new · submitted 2026-05-28 · 🌀 gr-qc · astro-ph.HE

Exact Mass Conservation in Binary Neutron Star Merger Simulations

Boris Daszuta , Sebastiano Bernuzzi , Joan Fontbut\'e , Ruocheng Zhai , Alan Tsz-Lok Lam , Jacob Fields , David Radice This is my paper

Pith reviewed 2026-06-29 05:43 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HE
keywords neutron star mergersmass conservationartificial atmospherenumerical relativitygravitational wavesejectahydrodynamics
0
0 comments X

The pith

A rescaling algorithm applied to the artificial atmosphere in neutron star merger simulations guarantees exact conservation of baryon mass and electron number to round-off precision.

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

The paper addresses the long-standing issue of baryon-mass violation introduced by artificial low-density atmospheres in Eulerian hydrodynamics simulations of neutron star spacetimes. It proposes a simple local rescaling algorithm that restores exact conservation of mass and electron number to round-off precision. The scheme is tested in binary neutron star merger simulations spanning multiple orbits and the postmerger phase with a microphysical equation of state, and can be combined with flux correction and a pseudo-vacuum treatment. If the claim holds, it removes a source of non-conservation in predictions of gravitational waves and ejected material, making simulations more reliable for interpreting observations.

Core claim

The authors show that their rescaling algorithm for the artificial atmosphere guarantees mass and electron number conservation to round-off precision. The pseudo-vacuum treatment shows slightly larger but approximately constant violations and improves computation of fast tail ejecta while providing convergent gravitational waves of quality comparable to the standard atmosphere. Results suggest that current computations of gravitational waves and dynamical ejecta are robust provided conservative adaptive mesh refinement with flux correction is employed.

What carries the argument

The local rescaling algorithm that adjusts the density in the artificial atmosphere to enforce global conservation constraints.

If this is right

  • Baryon mass and electron number conservation holds to round-off precision throughout the simulation.
  • The pseudo-vacuum option improves accuracy for fast tail ejecta.
  • Gravitational wave signals converge at quality levels comparable to standard treatments.
  • Merger simulations remain robust across different atmosphere treatments when using conservative AMR and flux correction.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This rescaling could be adapted for other Eulerian simulations that use artificial atmospheres, such as those involving black holes or accretion disks.
  • The approximately constant violations in the pseudo-vacuum case suggest that error accumulation is predictable and may not grow with simulation time.
  • Better ejecta modeling might lead to more accurate predictions for electromagnetic counterparts like kilonovae.

Load-bearing premise

That locally rescaling the artificial atmosphere density does not introduce unphysical effects that change the evolution of the neutron star matter or observable quantities.

What would settle it

A simulation test where the total baryon mass after rescaling differs from the initial mass by an amount larger than round-off error would falsify the exact conservation guarantee.

Figures

Figures reproduced from arXiv: 2605.30430 by Alan Tsz-Lok Lam, Boris Daszuta, David Radice, Jacob Fields, Joan Fontbut\'e, Ruocheng Zhai, Sebastiano Bernuzzi.

Figure 1
Figure 1. Figure 1: FIG. 1. Rest-mass density [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Evolution of the relative error in baryon (rest) mass and electron number conservation for the four vacuum algorithms [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Evolution of the [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Evolution of ejecta mass (top) and mass-averaged [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Mass-weighted histograms of ejecta angular distri [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Gravitational waves strain and instantaneous fre [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: All the atmosphere treatments give consistent [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Self convergence of the phase of the gravitational [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Gravitational waves phase (top) and amplitude (bot [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
read the original abstract

A long-standing problem in the simulation of neutron star spacetimes is the treatment of vacuum regions outside the stars. The use of an artificial low-density atmosphere is a common robust approach within Eulerian hydrodynamics that, however, introduces baryon-mass violation even with conservative numerical schemes. We propose a simple numerical algorithm that ensures exact mass conservation by means of an appropriate local rescaling of the atmosphere. The scheme is combined with a low-order flux correction and it can be further augmented by a pseudo-vacuum treatment that enforces strict vacuum in the outer regions far from the central objects. We demonstrate the effectiveness of these vacuum treatments with binary neutron star mergers simulations spanning multiple orbits and the postmerger phase, and including a microphysical equation of state. The rescaling algorithm guarantees mass and electron number conservation to round-off precision. The pseudo-vacuum treatment shows slightly larger but approximately constant violations and can improve the computation of fast tail ejecta as well as provide convergent gravitational waves of quality comparable to the standard atmosphere. Overall, results from different atmosphere treatments and a two-code comparison suggest that current computations of gravitational waves and (dynamical) ejecta in the presence of an artifical atmosphere are robust, provided that conservative adaptive mesh refinement with flux correction is employed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript presents a numerical algorithm for treating artificial low-density atmospheres in Eulerian hydrodynamics simulations of binary neutron star mergers. It proposes local rescaling of the atmosphere density after the conservative update, combined with low-order flux correction, to enforce exact baryon mass and electron number conservation to round-off precision. An optional pseudo-vacuum treatment is also introduced for outer regions. The scheme is tested on multi-orbit BNS mergers with a microphysical equation of state, claiming that the rescaling achieves the stated conservation while the pseudo-vacuum yields slightly larger but constant violations, and that results for gravitational waves and dynamical ejecta are comparable across treatments, indicating robustness when conservative AMR with flux correction is used.

Significance. If the central claims hold, the work addresses a persistent practical issue in numerical relativity by providing a simple, exact-conservation fix for atmosphere treatments without apparent loss of accuracy in key observables. The use of microphysical EOS and multi-orbit runs, plus the two-code comparison, supports broader applicability. The explicit demonstration of round-off conservation is a clear strength for reproducibility in the field.

major comments (1)
  1. [Abstract] Abstract: the claim that results from different atmosphere treatments are comparable (and thus that computations are robust) is invoked to support the overall conclusion, but the manuscript provides no quantitative bound on the magnitude of interface perturbations from local density rescaling before they would affect GW strains or ejecta properties at the level of the reported conservation gains. This assumption is load-bearing for the robustness statement.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of the work and the constructive comment. We address the major comment below and will incorporate revisions as indicated.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that results from different atmosphere treatments are comparable (and thus that computations are robust) is invoked to support the overall conclusion, but the manuscript provides no quantitative bound on the magnitude of interface perturbations from local density rescaling before they would affect GW strains or ejecta properties at the level of the reported conservation gains. This assumption is load-bearing for the robustness statement.

    Authors: We agree that the current manuscript does not supply an explicit quantitative bound on the size of interface perturbations from the rescaling that would begin to affect GW strains or ejecta at the scale of the reported conservation improvements. The robustness statement rests on the empirical finding that GW and ejecta quantities remain comparable across atmosphere treatments and between independent codes when conservative AMR with flux correction is used. While this provides practical evidence of robustness, a formal a priori bound would require additional analysis of the rescaling operator's effect on the solution that is not present. We will revise the abstract to moderate the language on robustness and add a short discussion quantifying the typical magnitude of the density adjustments introduced by the rescaling (which are localized and at the level of the atmosphere floor) together with their observed impact on the reported observables. revision: yes

Circularity Check

0 steps flagged

Numerical conservation scheme is self-contained with external verification

full rationale

The paper presents an explicit numerical algorithm that enforces exact mass and electron-number conservation via post-update local rescaling of the artificial atmosphere density, combined with low-order flux correction. This guarantee to round-off precision is a direct consequence of the rescaling construction itself rather than a derived prediction or first-principles result that reduces to fitted inputs. No load-bearing self-citations, uniqueness theorems, or ansatzes appear in the provided text; results are instead validated against independent external benchmarks including conservation laws, multi-orbit BNS simulations, microphysical EOS, and cross-code comparisons. The contribution therefore remains self-contained against external requirements without circular reduction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard framework of conservative Eulerian relativistic hydrodynamics plus the new rescaling step; the atmosphere density itself remains a tunable but conventional parameter.

free parameters (1)
  • atmosphere density floor
    Standard tunable low-density value used to fill vacuum regions; its specific choice is not altered by the rescaling but still affects where the algorithm is applied.
axioms (1)
  • domain assumption Eulerian hydrodynamics schemes remain conservative away from the artificial atmosphere
    The paper builds on the assumption that the underlying numerical scheme is conservative except where the atmosphere treatment intervenes.

pith-pipeline@v0.9.1-grok · 5778 in / 1384 out tokens · 36795 ms · 2026-06-29T05:43:19.125002+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 35 canonical work pages · 18 internal anchors

  1. [1]

    J. A. Fontet al., Phys. Rev.D65, 084024 (2002), arXiv:gr-qc/0110047

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  2. [2]

    Relativistic simulations of rotational core collapse. I. Methods, initial models, and code tests

    H. Dimmelmeier, J. A. Font, and E. M¨ uller, Astron. 10 Astrophys.388, 917 (2002), arXiv:astro-ph/0204288

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  3. [3]

    Baiotti, I

    L. Baiotti, I. Hawke, P. J. Montero, F. Loffler, L. Rez- zolla,et al., Phys.Rev.D71, 024035 (2005), arXiv:gr- qc/0403029 [gr-qc]

    work page arXiv 2005

  4. [4]

    Numerical relativity simulations of binary neutron stars

    M. Thierfelder, S. Bernuzzi, and B. Br¨ ugmann, Phys.Rev.D84, 044012 (2011), arXiv:1104.4751 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  5. [5]

    W. Cook, B. Daszuta, J. Fields, P. Hammond, S. Al- banesi, F. Zappa, S. Bernuzzi, and D. Radice, Astrophys. J. Suppl.277, 3 (2025), arXiv:2311.04989 [gr-qc]

    work page arXiv 2025

  6. [6]

    Kastaun, J

    W. Kastaun, J. V. Kalinani, and R. Ciolfi, Phys. Rev. D103, 023018 (2021), arXiv:2005.01821 [gr-qc]

    work page arXiv 2021

  7. [7]

    Hammond, I

    P. Hammond, I. Hawke, and N. Andersson, Phys. Rev. D104, 103006 (2021), arXiv:2108.08649 [astro-ph.HE]

    work page arXiv 2021

  8. [8]

    E. F. Toro,Riemann Solvers and Numerical Methods for Fluid Dynamics, 2nd ed. (Springer-Verlag, 1999)

    1999

  9. [9]

    R. J. LeVeque,Finite Volume Methods for Hyperbolic Problems(Cambridge University Press, 2002)

    2002

  10. [10]

    Chandrashekar and C

    P. Chandrashekar and C. Klingenberg, SIAM Jour- nal on Scientific Computing37, B382 (2015), https://doi.org/10.1137/140984373

    work page doi:10.1137/140984373 2015

  11. [11]

    J. P. Berberich, R. K¨ appeli, P. Chandrashekar, and C. Klingenberg, arXiv e-prints , arXiv:2005.01811 (2020), arXiv:2005.01811 [math.NA]

    work page arXiv 2005

  12. [12]

    P. V. F. Edelmann, L. Horst, J. P. Berberich, R. An- drassy, J. Higl, G. Leidi, C. Klingenberg, and F. K. R¨ opke, Astron. and Astrophys.652, A53 (2021), arXiv:2102.13111 [astro-ph.SR]

    work page arXiv 2021

  13. [13]

    W. E. East, F. Pretorius, and B. C. Stephens, Phys.Rev. D85, 124010 (2012), arXiv:1112.3094 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  14. [14]

    Three-Dimensional General-Relativistic Hydrodynamic Simulations of Binary Neutron Star Coalescence and Stellar Collapse with Multipatch Grids

    C. Reisswig, R. Haas, C. D. Ott, E. Abdikamalov, P. M¨ osta, D. Pollney, and E. Schnetter, Phys. Rev.D87, 064023 (2013), arXiv:1212.1191 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  15. [15]

    Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement

    T. Dietrich, S. Bernuzzi, M. Ujevic, and B. Br¨ ugmann, Phys. Rev.D91, 124041 (2015), arXiv:1504.01266 [gr- qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  16. [16]

    Poudel, W

    A. Poudel, W. Tichy, B. Br¨ ugmann, and T. Dietrich, Phys. Rev. D102, 104014 (2020), arXiv:2009.06617 [gr- qc]

    work page arXiv 2020

  17. [17]

    Tichy, L

    W. Tichy, L. Ji, A. Adhikari, A. Rashti, and M. Pirog, Class. Quant. Grav.40, 025004 (2023), arXiv:2212.06340 [gr-qc]

    work page arXiv 2023

  18. [18]

    Calibration of Moving Puncture Simulations

    B. Br¨ ugmann, J. A. Gonzalez, M. Hannam, S. Husa, U. Sperhake,et al., Phys.Rev.D77, 024027 (2008), arXiv:gr-qc/0610128 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  19. [19]

    M. N. Lemaster and J. M. Stone, Astrophys. J.691, 1092 (2009), arXiv:0809.4005 [astro-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  20. [20]

    J. M. Stone, P. D. Mullen, D. Fielding, P. Grete, M. Guo, P. Kempski, E. R. Most, C. J. White, and G. N. Wong, arXiv e-prints , arXiv:2409.16053 (2024), arXiv:2409.16053 [astro-ph.IM]

    work page arXiv 2024

  21. [21]

    Fields, H

    J. Fields, H. Zhu, D. Radice, J. M. Stone, W. Cook, S. Bernuzzi, and B. Daszuta, Astrophys. J. Suppl.276, 35 (2025), arXiv:2409.10384 [astro-ph.HE]

    work page arXiv 2025

  22. [22]

    GR-Athena++: Magnetohy- drodynamical Evolution with Dynamical Space-Time,

    B. Daszuta and W. Cook, “GR-Athena++: Magnetohy- drodynamical Evolution with Dynamical Space-Time,” (2025) arXiv:2406.05126 [gr-qc]

    work page arXiv 2025

  23. [23]

    Daszuta, W

    B. Daszuta, W. Cook, P. Hammond, J. Fields, E. M. Guti´ errez, S. Bernuzzi, and D. Radice, Phys. Rev. D 112, 103006 (2025), arXiv:2406.09139 [gr-qc]

    work page arXiv 2025

  24. [24]

    Banyuls, J

    F. Banyuls, J. A. Font, J. M. A. Ibanez, J. M. A. Marti, and J. A. Miralles, Astrophys. J.476, 221 (1997)

    1997

  25. [25]

    A. Neumaier, ZAMM - Journal of Applied Math- ematics and Mechanics / Zeitschrift f¨ ur Ange- wandte Mathematik und Mechanik54, 39 (1974), https://onlinelibrary.wiley.com/doi/pdf/10.1002/zamm.19740540106

    work page doi:10.1002/zamm.19740540106 1974

  26. [26]

    Daszuta, S

    B. Daszuta, S. Bernuzzi, M. Jacobi, E. M. Guti´ errez, P. Hammond, W. Cook, and D. Radice, (2026), arXiv:2602.18290 [gr-qc]

    work page arXiv 2026

  27. [27]

    Borges, M

    R. Borges, M. Carmona, B. Costa, and W. S. Don, Jour- nal of Computational Physics227, 3191 (2008)

    2008

  28. [28]

    L. D. Zanna and N. Bucciantini, Astron. Astrophys.390, 1177 (2002), arXiv:astro-ph/0205290

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  29. [29]

    van Leer, Journal of Computational Physics14, 361 (1974)

    B. van Leer, Journal of Computational Physics14, 361 (1974)

    1974

  30. [30]

    Quasiequilibrium sequences of synchronized and irrotational binary neutron stars in general relativity. I. Method and tests

    E. Gourgoulhon, P. Grandclement, K. Taniguchi, J.- A. Marck, and S. Bonazzola, Phys.Rev.D63, 064029 (2001), arXiv:gr-qc/0007028 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2001

  31. [31]

    New Equations of State in Simulations of Core-Collapse Supernovae

    M. Hempel, T. Fischer, J. Schaffner-Bielich, and M. Liebendorfer, Astrophys. J.748, 70 (2012), arXiv:1108.0848 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  32. [32]

    Gottlieb and C.-W

    S. Gottlieb and C.-W. Ketcheson, David I.and Shu, Jour- nal of Scientific Computing38, 251 (2009)

    2009

  33. [33]

    Perego, S

    A. Perego, S. Bernuzzi, and D. Radice, Eur. Phys. J. A55, 124 (2019), arXiv:1903.07898 [gr-qc]

    work page arXiv 2019

  34. [34]

    Nedora, S

    V. Nedora, S. Bernuzzi, D. Radice, A. Perego, A. En- drizzi, and N. Ortiz, Astrophys. J.886, L30 (2019), arXiv:1907.04872 [astro-ph.HE]

    work page arXiv 2019

  35. [35]

    Constraint violation in free evolution schemes: comparing BSSNOK with a conformal decomposition of Z4

    S. Bernuzzi and D. Hilditch, Phys. Rev.D81, 084003 (2010), arXiv:0912.2920 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  36. [36]

    Compact binary evolutions with the Z4c formulation

    D. Hilditch, S. Bernuzzi, M. Thierfelder, Z. Cao, W. Tichy, and B. Bruegmann, Phys. Rev.D88, 084057 (2013), arXiv:1212.2901 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  37. [37]

    Constraint damping for the Z4c formulation of general relativity

    A. Weyhausen, S. Bernuzzi, and D. Hilditch, Phys. Rev. D85, 024038 (2012), arXiv:1107.5539 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  38. [38]

    Binary Neutron Star Mergers: Mass Ejection, Electromagnetic Counterparts and Nucleosynthesis

    D. Radice, A. Perego, K. Hotokezaka, S. A. Fromm, S. Bernuzzi, and L. F. Roberts, Astrophys. J.869, 130 (2018), arXiv:1809.11161 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  39. [39]

    Properties of hypermassive neutron stars formed in mergers of spinning binaries

    W. Kastaun and F. Galeazzi, Phys.Rev.D91, 064027 (2015), arXiv:1411.7975 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  40. [40]

    Notes on the integration of numerical relativity waveforms

    C. Reisswig and D. Pollney, Class.Quant.Grav.28, 195015 (2011), arXiv:1006.1632 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  41. [41]

    High-Order Fully General-Relativistic Hydrodynamics: new Approaches and Tests

    D. Radice, L. Rezzolla, and F. Galeazzi, Class.Quant.Grav.31, 075012 (2014), arXiv:1312.5004 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2014