arXiv preprint ·

RAKI: High-Order Alternative Finite Differences for General Relativistic Magnetohydrodynamics

Rahime Matur, Beyhan Karakas, Ian Hawke and Nils Andersson

Published
Source
arXiv preprint, posted before peer review
Original
Read it at arXiv
Topics
High Energy Astrophysical PhenomenaGeneral Relativity & Quantum CosmologyHigh Energy Physics - PhenomenologyComputational Physics

Abstract

Binary neutron star (BNS) and black hole--neutron star mergers are among the most important multimessenger sources. Observations of these events could provide valuable insights into several unresolved phenomena, ranging from the efficiency of magnetic field amplification during and after mergers to the formation of magnetars and jets. Understanding these phenomena requires accurate modelling of magnetic field evolution, which remains a challenge for state-of-the-art general relativistic magnetohydrodynamics (GRMHD) simulations. In this work, we present \texttt{RAKI} (Relativistic Astrophysics Code Implementation of the Alternative Finite Difference Approach for GRMHD), the first GRMHD code to simulate BNS mergers using high-order alternative finite-difference methods with adaptive-order WENO (AFD-WENO-AO) reconstruction and a fully non-staggered, node-centred vector-potential formulation. The code achieves up to ninth-order convergence in smooth test problems while accurately capturing discontinuities. In magnetised Tolman--Oppenheimer--Volkoff (TOV) star tests, \texttt{RAKI} retains $99.2$ per cent of the initial maximum magnetic field strength over ${\sim}2$ ms, exhibiting an order of magnitude less magnetic field loss than a reference GRMHD code at the same resolution. In BNS merger simulations, fifth-order AFD-WENO-AO reconstruction yields approximately third-order convergence of the gravitational-wave phase in both magnetised and unmagnetised models. The discrete no-monopole constraint is preserved to machine precision without staggered grids, even during BNS mergers with moving refinement regions. These results establish \texttt{RAKI} as a high-order, non-staggered GRMHD framework for accurate modelling of magnetic field evolution in neutron star mergers and its implications for multimessenger astrophysics.

The authors’ abstract, as arXiv lists it. arXiv’s listings, abstracts included, are free for anyone to reuse.

Read it on ReadSearch

Each opens after you create an account, and brings you back to this paper.

Already a reader? Sign in

All new papers in High Energy Astrophysical Phenomena

Get the few papers worth your time in High Energy Astrophysical Phenomena every morning.

ReadSearch reads the new papers in your field each day and brings you the few worth reading, each with its key finding in a sentence.

Papers belong to their authors and publishers; every page here links to the original. All research topics