We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from the Web site stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth, which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

• Received 18 March 2014

DOI:https://doi.org/10.1103/PhysRevD.89.104029