Numerical Simulations of Realistic Binary Neutron Stars

The new field of Multi-Messenger Astrophysics adds gravitational waves to the detection of electromagnetic signals, enabling us to gain previously unavailable insights about our Universe. It started with the detection in multiple wavelengths of the neutron star merger (GW/GRB 170817A), and proved its extraordinary potential in answering fundamental questions about the characteristics of neutron stars. The promise of MMA can be realized only with accurate theoretical models to decode the many unknown mechanisms leading to powerful electromagnetic signals such as fast radio bursts, cosmic jets, and gamma ray bursts that might be present during the merger. This work seeks to discern the role played by the magnetic field geometry and strength, equation of state and mass ratio during the chaotic dynamics of the merger, in generating or shutting off detectable multi-messenger signals from binary neutron star collisions. Specifically, we carry out numerical simulations of realistic binary neutron stars, to investigate the key factors responsible for steering the production of detectable electromagnetic counterparts to the gravitational wave signals. This work addresses also the reproducibility of results by using publicly available, open-source software and documenting the workflow, which will enable researchers and students to regenerate the same findings and produce new results.