Tabulated equation of state support within IllinoisGRMHD

Presently, the best tools for understanding BNS mergers in their most extreme stages are detailed general relativistic magneto-hydrodynamic (GRMHD) simulations. Although recent years have seen the development and improvement of accurate GRMHD codes, only a few of these are open-source. Among these is the widely used IllinoisGRMHD code, which in its present open-source state lacks the ability to simulate fluids with the use of multi-dimensional, finite temperature equation of state tables. We present new developments to IllinoisGRMHD which allow for its use with realistic, finite temperature equations of state. Our adaptation of the code includes advection of the electron fraction, which is crucial in understanding the nucleosynthetic processes associated with BNS mergers, and new conservative-to-primitive recovery routines which are based on modern solution algorithms. We discuss stringent tests of our version of IllinoisGRMHD and present detailed comparisons to other open-source GRMHD codes. We consider the evolution and amplification of magnetic fields in BNS systems described using finite temperature equations of state to showcase the full capability of our code. We find that our version of IllinoisGRMHD behaves as expected, and similar to other open-source GRMHD codes, when considering many systems of astrophysical interest.