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
SN - 0004-637X
VL - 919
JO - Astrophysical Journal
JF - Astrophysical Journal
IS - 2
M1 - 95
ER -