The error budget of binary neutron star merger simulations for configurations with high spin

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Kuan, Hao-Jui, Markin, Ivan, Ujevic, Maximiliano, Dietrich, Tim, Kiuchi, Kenta, Shibata, Masaru, Tichy, Wolfgang
Format: Preprint
Veröffentlicht: 2025
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866913871982231552
author Kuan, Hao-Jui
Markin, Ivan
Ujevic, Maximiliano
Dietrich, Tim
Kiuchi, Kenta
Shibata, Masaru
Tichy, Wolfgang
author_facet Kuan, Hao-Jui
Markin, Ivan
Ujevic, Maximiliano
Dietrich, Tim
Kiuchi, Kenta
Shibata, Masaru
Tichy, Wolfgang
contents Numerical-relativity simulations offer a unique approach to investigating the dynamics of binary neutron star mergers and provide the most accurate predictions of waveforms in the late inspiral phase. However, the numerical predictions are prone to systematic biases originating from the construction of initial quasi-circular binary configurations, the numerical methods used to evolve them, and to extract gravitational signals. To assess uncertainties arising from these aspects, we analyze mergers of highly spinning neutron stars with dimensionless spin parameter $χ=0.5$. The initial data are prepared by two solvers, \textsc{FUKA} and \textsc{SGRID}, which are then evolved by two independent codes, \textsc{SACRA} and \textsc{BAM}. We assess the impact of numerical discretizations, finite extraction radii, and differences in numerical frameworks on the resulting gravitational waveforms. Our analysis reveals that the primary source of uncertainty in numerical waveforms is the evolution code, while the initial data solver has a smaller impact. We also compare our numerical-relativity waveforms with state-of-the-art analytical models, finding that the discrepancies between them exceed the estimated numerical uncertainties. Few suggestions are offered: (i) the analytic waveform becomes an inadequate approximation after the two neutron stars come into contact and the binary enters the essentially-one-body phase, (ii) the analytical models may not capture finite-size effects beyond quadrupole moment, and (iii) the inconsistent use of the binary black hole baseline in the analytical models may also be contributing to these discrepancies. The presented results benchmark the error budget for numerical waveforms of binary neutron star mergers, and provide information for the analytic models to explore further the high spin parameter space of binary neutron star mergers.
format Preprint
id arxiv_https___arxiv_org_abs_2506_02115
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle The error budget of binary neutron star merger simulations for configurations with high spin
Kuan, Hao-Jui
Markin, Ivan
Ujevic, Maximiliano
Dietrich, Tim
Kiuchi, Kenta
Shibata, Masaru
Tichy, Wolfgang
General Relativity and Quantum Cosmology
High Energy Astrophysical Phenomena
Numerical-relativity simulations offer a unique approach to investigating the dynamics of binary neutron star mergers and provide the most accurate predictions of waveforms in the late inspiral phase. However, the numerical predictions are prone to systematic biases originating from the construction of initial quasi-circular binary configurations, the numerical methods used to evolve them, and to extract gravitational signals. To assess uncertainties arising from these aspects, we analyze mergers of highly spinning neutron stars with dimensionless spin parameter $χ=0.5$. The initial data are prepared by two solvers, \textsc{FUKA} and \textsc{SGRID}, which are then evolved by two independent codes, \textsc{SACRA} and \textsc{BAM}. We assess the impact of numerical discretizations, finite extraction radii, and differences in numerical frameworks on the resulting gravitational waveforms. Our analysis reveals that the primary source of uncertainty in numerical waveforms is the evolution code, while the initial data solver has a smaller impact. We also compare our numerical-relativity waveforms with state-of-the-art analytical models, finding that the discrepancies between them exceed the estimated numerical uncertainties. Few suggestions are offered: (i) the analytic waveform becomes an inadequate approximation after the two neutron stars come into contact and the binary enters the essentially-one-body phase, (ii) the analytical models may not capture finite-size effects beyond quadrupole moment, and (iii) the inconsistent use of the binary black hole baseline in the analytical models may also be contributing to these discrepancies. The presented results benchmark the error budget for numerical waveforms of binary neutron star mergers, and provide information for the analytic models to explore further the high spin parameter space of binary neutron star mergers.
title The error budget of binary neutron star merger simulations for configurations with high spin
topic General Relativity and Quantum Cosmology
High Energy Astrophysical Phenomena
url https://arxiv.org/abs/2506.02115