-
Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD
Authors:
Leonardo R. Werneck,
Zachariah B. Etienne,
Ariadna Murguia-Berthier,
Roland Haas,
Federico Cipolletta,
Scott C. Noble,
Lorenzo Ennoggi,
Federico G. Lopez Armengol,
Bruno Giacomazzo,
Thiago Assumpção,
Joshua Faber,
Tanmayee Gupte,
Bernard J. Kelly,
Julian H. Krolik
Abstract:
We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which imp…
▽ More
We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit.
△ Less
Submitted 14 December, 2023; v1 submitted 30 August, 2022;
originally announced August 2022.
-
Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations
Authors:
Federico G. Lopez Armengol,
Zachariah B. Etienne,
Scott C. Noble,
Bernard J. Kelly,
Leonardo R. Werneck,
Brendan Drachler,
Manuela Campanelli,
Federico Cipolletta,
Yosef Zlochower,
Ariadna Murguia-Berthier,
Lorenzo Ennoggi,
Mark Avara,
Riccardo Ciolfi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Julian H. Krolik,
Vassilios Mewes,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril,
Jeremy Schnittman
Abstract:
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg…
▽ More
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach.
△ Less
Submitted 31 October, 2022; v1 submitted 17 December, 2021;
originally announced December 2021.
-
HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage
Authors:
Ariadna Murguia-Berthier,
Scott C. Noble,
Luke F. Roberts,
Enrico Ramirez-Ruiz,
Leonardo R. Werneck,
Michael Kolacki,
Zachariah B. Etienne,
Mark Avara,
Manuela Campanelli,
Riccardo Ciolfi,
Federico Cipolletta,
Brendan Drachler,
Lorenzo Ennoggi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Bernard J. Kelly,
Julian H. Krolik,
Federico G. Lopez Armengol,
Ben Margalit,
Tim Moon,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril
, et al. (3 additional authors not shown)
Abstract:
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge…
▽ More
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.
△ Less
Submitted 31 January, 2022; v1 submitted 9 June, 2021;
originally announced June 2021.
