TY - JOUR

T1 - HARM3D+NUC

T2 - A New Method for Simulating the Post-merger Phase of Binary Neutron Star Mergers with GRMHD, Tabulated EOS, and Neutrino Leakage

AU - Murguia-Berthier, Ariadna

AU - Noble, Scott C.

AU - Roberts, Luke F.

AU - Ramirez-Ruiz, Enrico

AU - Werneck, Leonardo R.

AU - Kolacki, Michael

AU - Etienne, Zachariah B.

AU - Avara, Mark

AU - Campanelli, Manuela

AU - Ciolfi, Riccardo

AU - Cipolletta, Federico

AU - Drachler, Brendan

AU - Ennoggi, Lorenzo

AU - Faber, Joshua

AU - Fiacco, Grace

AU - Giacomazzo, Bruno

AU - Gupte, Tanmayee

AU - Ha, Trung

AU - Kelly, Bernard J.

AU - Krolik, Julian H.

AU - Armengol, Federico G. Lopez

AU - Ben Margalit, null

AU - Moon, Tim

AU - O'Shaughnessy, Richard

AU - Rueda-Becerril, Jesus M.

AU - Schnittman, Jeremy

AU - Zenati, Yossef

AU - Zlochower, Yosef

PY - 2021/9/29

Y1 - 2021/9/29

N2 - 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 gamma-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.

AB - 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 gamma-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.

KW - COOLED ACCRETION DISKS

KW - R-PROCESS NUCLEOSYNTHESIS

KW - GAMMA-RAY BURSTS

KW - MASS EJECTION

KW - DRIVEN WINDS

KW - BLACK-HOLES

KW - EVOLUTION

KW - TRANSPORT

KW - MODELS

KW - KILONOVA

U2 - 10.3847/1538-4357/ac1119

DO - 10.3847/1538-4357/ac1119

M3 - Journal article

VL - 919

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 2

M1 - 95

ER -