Quantifying uncertainties in general relativistic magnetohydrodynamic codes

Pedro L. Espino, Gabriele Bozzola, Vasileios Paschalidis

Research output: Contribution to journalArticlepeer-review

2 Scopus citations


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.

Original languageEnglish (US)
Article number104059
JournalPhysical Review D
Issue number10
StatePublished - May 15 2023

ASJC Scopus subject areas

  • Nuclear and High Energy Physics


Dive into the research topics of 'Quantifying uncertainties in general relativistic magnetohydrodynamic codes'. Together they form a unique fingerprint.

Cite this