General Relativistic Magneto-Hydrodynamic Simulations with BAM: Implementation and Code Comparison

Anna Neuweiler,Tim Dietrich,Bernd Brügmann,Edoardo Giangrandi,Kenta Kiuchi,Federico Schianchi,Philipp Mösta,Swapnil Shankar,Bruno Giacomazzo,Masaru Shibata

2024-07-31

Abstract:Binary neutron star mergers are among the most energetic events in our Universe, with magnetic fields significantly impacting their dynamics, particularly after the merger. While numerical-relativity simulations that correctly describe the physics are essential to model their rich phenomenology, the inclusion of magnetic fields is crucial for realistic simulations. For this reason, we have extended the BAM code to enable general relativistic magneto-hydrodynamic (GRMHD) simulations employing a hyperbolic `divergence cleaning' scheme. We present a large set of standard GRMHD tests and compare the BAM code to other GRMHD codes, SPRITZ, GRaM-X, and SACRA$_{\rm KK22}$, which employ different schemes for the evolution of the magnetic fields. Overall, we find that the BAM code shows a good performance in simple special-relativistic tests. In addition, we find good agreement and consistent results when comparing GRMHD simulation results between BAM and SACRA$_{\rm KK22}$.

General Relativity and Quantum Cosmology,High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is the numerical simulation of magnetohydrodynamic effects during the merger of double neutron stars. Specifically, the focus of the paper is as follows:

1. **Extension of the bam code**: In order to conduct general - relativistic magnetohydrodynamic (GRMHD) simulations, the researchers extended the existing bam code and introduced the hyperbolic "divergence cleaning" scheme to handle the evolution of the magnetic field. This extension is crucial for achieving more realistic simulations because the magnetic field plays an important role in the dynamic process after the merger of double neutron stars.

2. **Evaluation of code performance**: Through a series of standard GRMHD tests, the researchers evaluated the performance of the bam code in simple special - relativistic tests and compared it with other GRMHD codes (such as spritz, GRaM - X, and sacra KK22) that use different magnetic field evolution schemes. These tests not only verified the accuracy of the new implementation but also evaluated the performance differences with other mature codes.

3. **Handling of the divergence - free condition of the magnetic field**: The paper discussed in detail how to maintain the divergence - free condition of the magnetic field in numerical simulations, which is a key challenge in GRMHD simulations. The researchers adopted the hyperbolic divergence cleaning method. Although this method cannot completely avoid divergence errors, it can effectively reduce the influence of these errors while maintaining the simplicity and high - order accuracy of the algorithm.

4. **Simulation of double - neutron - star mergers**: Finally, the paper presented the results of double - neutron - star merger simulations using the bam code and compared them with the results of the sacra KK22 code. These simulations not only verified the effectiveness of the code but also provided important data support for understanding the physical phenomena during the double - neutron - star merger process, such as matter ejection and the generation of gravitational wave signals.

In summary, this paper aims to improve and verify the GRMHD functions of the bam code to enhance the understanding and prediction ability of double - neutron - star mergers and their subsequent physical processes. This is not only of great significance for astrophysics research but also provides strong support for the development of multi - messenger astronomy.