REVIEW 7 cited by
GR-Athena++: General-relativistic magnetohydrodynamics simulations of neutron star spacetimes
Not yet reviewed by Pith; the record is open.
This paper has not been read by Pith yet. Machine review is queued; the pith claim, tier, and objections will appear here once it completes.
SPECIMEN: schema-true, not a live event
T0 review · schema-true
One-sentence machine reading of the paper's core claim.
pith:XXXXXXXX · record.json · timestamp
abstract
We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in non-rotating stars with total relative violation of the divergence-free constraint of $10^{-16}$. It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin-Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above $80\%$ in excess of $10^5$ CPU cores and excellent weak scaling is shown up to $\sim 5 \times 10^5$ CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exascale computers.
Forward citations
Cited by 7 Pith papers
-
3D Binary Neutron Star Merger Ejecta Evolution up to Seconds Timescale: Dynamics, Element Distribution, and Light Curves
Long 3D simulations of neutron star merger ejecta show radioactive heating keeps reshaping heavy-element outflows, and 3D light curves are dimmer but no closer to AT2017gfo than 2D ones.
-
Transport Properties of the MRI in Differentially Rotating Neutron Stars
Saturated MRI turbulence in differentially rotating neutron stars yields ℓ_mix ≈ (0.01–0.1) λ_MRI, largely independent of density, so standard GRLES mixing-length prescriptions overestimate transport by about an order...
-
GRACE: An Open-Source Framework for GPU-Accelerated Numerical Relativity
GRACE is a validated, open-source, Kokkos+p4est GPU-portable framework that evolves ideal GRMHD with constrained transport self-consistently coupled to Z4c Einstein equations on fixed or adaptive meshes.
-
Magnetic Field Configurations in Binary Neutron Star Mergers II: Inspiral, Merger and Ejecta
Initial magnetic field topology, especially anti-aligned poloidal fields, strongly controls post-merger field amplification and ejecta magnetisation in neutron star merger simulations.
-
Extending the infrastructure of the BAM code towards resistive general-relativistic magnetohydrodynamics: tests and first applications
BAM gains a validated resistive GRMHD module; BNS runs show the parallel electric field reaching up to 10% of total strength in low-density disks, unlike ideal MHD.
-
Magnetic field dynamics in isolated neutron stars with an external dipole field
Long-term numerical relativity simulations find that neutron star magnetic fields relax to stable mixed configurations with toroidal energy fraction ≲10% within one Alfvén time after Tayler instability saturation.
-
Impact of rotation on magnetic field stability and orientation in isolated neutron stars
In 3D GRMHD simulations, neutron star rotation delays Tayler, kink, and Parker instabilities, so faster-spinning models retain more magnetic energy over 10 Alfven times.
Discussion (0). Continue with ORCID to comment.