Extending the infrastructure of the BAM code towards resistive general-relativistic magnetohydrodynamics: tests and first applications
Abstract
Many astrophysical phenomena, including pulsars and short gamma-ray bursts, are associated with the extremely strong magnetic fields present in neutron stars and neutron star mergers. While the ideal magnetohydrodynamic approximation, which assumes infinite conductivity, provides an excellent description of the neutron-star interior, it cannot capture non-ideal processes such as Ohmic dissipation and magnetic reconnection. To overcome this limitation, we present an extension of the numerical-relativity code BAM incorporating a resistive general-relativistic magnetohydrodynamic (GRMHD) description. We validate the new implementation through an extensive suite of special-relativistic magnetohydrodynamic benchmark tests and by performing stable simulations of isolated and binary neutron star systems. For the latter, we investigate the impact of finite conductivity on magnetic-field amplification, mass ejection, and non-ideal GRMHD effects. In particular, we find that the component of the electric field parallel to the magnetic field, which is zero in the ideal case, can reach up to 10% of the total electric field strength. This highlights the potential importance of non-ideal effects for accurately modeling the long-term evolution of post-merger remnants, particularly in low-density regions. Although the present study is restricted to simplified piecewise-polytropic equations of state, it demonstrates the capabilities of the new resistive GRMHD framework and paves the way for future investigations employing more realistic microphysics.
Versions
Related
- Chana: A GPU-Accelerated Adaptive-Mesh Code for General Relativistic Magnetohydrodynamics with Multigroup M1 Neutrino Transport
- Relativistic Magnetohydrodynamics from Myers-Pospelov Lorentz-Violating Electrodynamics
- FUGUE: A New GPU-accelerated General Relativistic Magnetohydrodynamics Code with Weak-Field Multipole Self-Gravity and Quasi-Stationary Kerr Black Hole Evolution
- Linear causality and stability constraints on relativistic second-order magnetohydrodynamics
- An HLLD Implementation for General Relativistic Magnetohydrodynamics in AthenaK
BibTeX
@misc{beaudoin2026extending,
title = {{Extending the infrastructure of the BAM code towards resistive general-relativistic magnetohydrodynamics: tests and first applications}},
author = {Matthew Beaudoin and Maximilano Ujevic and Ramon Jaeger and Anna Neuweiler and Henrique Gieg and Kenta Kiuchi and Tim Dietrich},
howpublished = {arXiv},
year = {2026},
month = jul,
doi = {10.48550/arxiv.2607.11670},
eprint = {2607.11670},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.11670},
note = {CC BY-SA 4.0},
}
RIS
TY - UNPB
TI - Extending the infrastructure of the BAM code towards resistive general-relativistic magnetohydrodynamics: tests and first applications
AU - Matthew Beaudoin
AU - Maximilano Ujevic
AU - Ramon Jaeger
AU - Anna Neuweiler
AU - Henrique Gieg
AU - Kenta Kiuchi
AU - Tim Dietrich
PY - 2026
DA - 2026/07/13
PB - arXiv
DO - 10.48550/arxiv.2607.11670
UR - https://arxiv.org/abs/2607.11670
AB - Many astrophysical phenomena, including pulsars and short gamma-ray bursts, are associated with the extremely strong magnetic fields present in neutron stars and neutron star mergers. While the ideal magnetohydrodynamic approximation, which assumes infinite conductivity, provides an excellent description of the neutron-star interior, it cannot capture non-ideal processes such as Ohmic dissipation and magnetic reconnection. To overcome this limitation, we present an extension of the numerical-relativity code BAM incorporating a resistive general-relativistic magnetohydrodynamic (GRMHD) description. We validate the new implementation through an extensive suite of special-relativistic magnetohydrodynamic benchmark tests and by performing stable simulations of isolated and binary neutron star systems. For the latter, we investigate the impact of finite conductivity on magnetic-field amplification, mass ejection, and non-ideal GRMHD effects. In particular, we find that the component of the electric field parallel to the magnetic field, which is zero in the ideal case, can reach up to 10% of the total electric field strength. This highlights the potential importance of non-ideal effects for accurately modeling the long-term evolution of post-merger remnants, particularly in low-density regions. Although the present study is restricted to simplified piecewise-polytropic equations of state, it demonstrates the capabilities of the new resistive GRMHD framework and paves the way for future investigations employing more realistic microphysics.
ER -