An efficient approach to resistive GRMHD simulations of binary neutron star mergers
Abstract
Astrophysical plasmas are often highly conducting, but finite resistivity can nevertheless influence the evolution of strongly magnetised systems. Capturing these effects in general relativistic magnetohydrodynamics is computationally challenging because high conductivity leads to stiff source terms in the evolution equations. Standard numerical approaches rely on implicit-explicit (IMEX) time integration to maintain stability, but this substantially increases the complexity and computational cost of the simulations. Building upon the previously developed Resistive Extension upGrade for Ideal MagnEtohydrodynamics (REGIME), we present a new fluid-frame formulation that simplifies the incorporation of leading-order resistive effects by avoiding matrix inversions required in previous approaches. The method is implemented in a general relativistic magnetohydrodynamics code and validated through one-dimensional current-sheet tests and simulations of isolated magnetised neutron stars and binary neutron star mergers. We find that the new formulation reproduces the expected resistive behaviour with accuracy comparable to IMEX-based methods, while substantially reducing the computational cost and complexity of the simulations. In binary neutron star mergers, finite resistivity leads to measurable differences in the post-merger gravitational-wave signal and remnant structure, showing that resistive effects can produce measurable changes in the post-merger dynamics even in the high-conductivity regime.
Related
- Charged-current neutrino opacity within the relativistic Hartree-Fock framework for astrophysical simulations of core-collapse supernovae and binary neutron star mergers
- Nuclear Physics of Binary Neutron Star Mergers
- The imprint of individual neutron star spins on ejecta, r -process nucleosynthesis, and kilonovae in binary neutron star mergers
- Magnetic Eruption and Nucleosynthesis in GRνMHD Simulations of Spinning Neutron Star Mergers
- Detectability of bulk viscosity effects on post-merger gravitational wave signals from binary neutron star mergers
BibTeX
@misc{miravettenes2026efficient,
title = {{An efficient approach to resistive GRMHD simulations of binary neutron star mergers}},
author = {Miquel Miravet-Tenés and Ian Hawke and Rahime Matur and Alexander J. Wright},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.12998},
eprint = {2609.12998},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.12998},
}
RIS
TY - UNPB
TI - An efficient approach to resistive GRMHD simulations of binary neutron star mergers
AU - Miquel Miravet-Tenés
AU - Ian Hawke
AU - Rahime Matur
AU - Alexander J. Wright
PY - 2026
DA - 2026/09/11
PB - arXiv
DO - 10.48550/arxiv.2609.12998
UR - https://arxiv.org/abs/2609.12998
AB - Astrophysical plasmas are often highly conducting, but finite resistivity can nevertheless influence the evolution of strongly magnetised systems. Capturing these effects in general relativistic magnetohydrodynamics is computationally challenging because high conductivity leads to stiff source terms in the evolution equations. Standard numerical approaches rely on implicit-explicit (IMEX) time integration to maintain stability, but this substantially increases the complexity and computational cost of the simulations. Building upon the previously developed Resistive Extension upGrade for Ideal MagnEtohydrodynamics (REGIME), we present a new fluid-frame formulation that simplifies the incorporation of leading-order resistive effects by avoiding matrix inversions required in previous approaches. The method is implemented in a general relativistic magnetohydrodynamics code and validated through one-dimensional current-sheet tests and simulations of isolated magnetised neutron stars and binary neutron star mergers. We find that the new formulation reproduces the expected resistive behaviour with accuracy comparable to IMEX-based methods, while substantially reducing the computational cost and complexity of the simulations. In binary neutron star mergers, finite resistivity leads to measurable differences in the post-merger gravitational-wave signal and remnant structure, showing that resistive effects can produce measurable changes in the post-merger dynamics even in the high-conductivity regime.
ER -