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 (BH) has become nearly stationary, the evolution of the surrounding accretion disk on cartesian grids over long timescales (∼1s) is suboptimal, as the plasma flows obliquely across coordinate lines and numerical errors quickly dissipate angular momentum 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 BH 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.