TY - JOUR
T1 - Quantifying uncertainties in general relativistic magnetohydrodynamic codes
AU - Espino, Pedro L.
AU - Bozzola, Gabriele
AU - Paschalidis, Vasileios
N1 - Publisher Copyright:
© 2023 American Physical Society.
PY - 2023/5/15
Y1 - 2023/5/15
N2 - In this paper, we show that similar open-source codes for general relativistic (magneto)hydrodynamic [GR(M)HD] produce different results for key features of binary neutron star mergers. First, we present a new open-source version of the publicly available illinoisgrmhd code that provides support for realistic, finite temperature equations of state. After stringent tests of our upgraded code, we perform a code comparison between grhydro, illinoisgrmhd, spritz, and whiskythc, which implement the same physics, but slightly different computational methods. The benefit of the comparison is that all codes are embedded in the einstein toolkit suite, hence their only difference is algorithmic. We find similar convergence properties, fluid dynamics, and gravitational waves, but different merger times, remnant lifetimes, and gravitational wave phases. Such differences must be resolved before the postmerger dynamics modeled with such simulations can be reliably used to infer the properties of nuclear matter especially in the era of precision gravitational wave astronomy.
AB - In this paper, we show that similar open-source codes for general relativistic (magneto)hydrodynamic [GR(M)HD] produce different results for key features of binary neutron star mergers. First, we present a new open-source version of the publicly available illinoisgrmhd code that provides support for realistic, finite temperature equations of state. After stringent tests of our upgraded code, we perform a code comparison between grhydro, illinoisgrmhd, spritz, and whiskythc, which implement the same physics, but slightly different computational methods. The benefit of the comparison is that all codes are embedded in the einstein toolkit suite, hence their only difference is algorithmic. We find similar convergence properties, fluid dynamics, and gravitational waves, but different merger times, remnant lifetimes, and gravitational wave phases. Such differences must be resolved before the postmerger dynamics modeled with such simulations can be reliably used to infer the properties of nuclear matter especially in the era of precision gravitational wave astronomy.
UR - http://www.scopus.com/inward/record.url?scp=85161141707&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85161141707&partnerID=8YFLogxK
U2 - 10.1103/PhysRevD.107.104059
DO - 10.1103/PhysRevD.107.104059
M3 - Article
AN - SCOPUS:85161141707
SN - 2470-0010
VL - 107
JO - Physical Review D
JF - Physical Review D
IS - 10
M1 - 104059
ER -