Abstract:One of the main challenges in the numerical modeling of binary neutron-star mergers are long-term simulations of the postmerger remnant over timescales of the order of seconds. When this modeling includes all the aspects of the complex physics accompanying the remnant, the computational costs can easily become enormous. To address this challenge in part, we have developed a novel hybrid approach in which the solution from a general-relativistic magnetohydrodynamics (GRMHD) code solving the full set of the Einstein equations in Cartesian coordinates is coupled with another GRMHD code in which the Einstein equations are solved under the conformally flat condition (CFC). The latter approximation has a long history and has been shown to provide an accurate description of compact objects in nonvacuum spacetimes. An important aspect of the CFC approximation is that the elliptic equations need to be solved only for a fraction of the steps needed for the underlying hydrodynamical/magnetohydrodynamical evolution, thus allowing for a gain in computational efficiency that can be up to a factor of ∼6(230) in three-dimensional (two-dimensional) simulations. We present here the basic features of the new code, the strategies necessary to interface it when importing both two- and three-dimensional data, and a novel and robust approach to the recovery of the primitive variables. To validate our new framework, we have carried out a number of tests with various coordinates systems and different numbers of spatial dimensions, involving a variety of astrophysical scenarios, including the evolution of the postmerger remnant of a binary neutron-star merger over a timescale of one second. Overall, our results show that the new code, bhac +, is able to accurately reproduce the evolution of compact objects in nonvacuum spacetimes and that, when compared with the evolution in full general relativity, the CFC approximation reproduces accurately both the gravitational fields and the matter variables at a fraction of the computational costs. This opens the way for the systematic study of the secular matter and electromagnetic emission from binary-merger remnants. https://doi.org/10.1103/PhysRevD.109.064061 © 2024 American Physical Society

Improved moving-puncture techniques for compact binary simulations

Improved Moving-Puncture Techniques for Compact Binary Simulations

Insights into Binary Neutron Star Merger Simulations: A Multi-Code Comparison

Improved methods for simulating nearly extremal binary black holes

Numerical Relativity Multimodal Waveforms using Absorbing Boundary Conditions

Simulating binary black hole mergers using discontinuous Galerkin methods

New public code for initial data of unequal-mass, spinning compact-object binaries

Fast and faithful interpolation of numerical relativity surrogate waveforms using a meshfree approximation

Accurate and Efficient Waveform Model for Precessing Binary Black Holes

Thinking Outside the Box: Numerical Relativity with Particles

Initial Data and Eccentricity Reduction Toolkit for Binary Black Hole Numerical Relativity Waveforms

Adaptive mesh refinement in binary black holes simulations

Relieving scale disparity in binary black hole simulations

A hybrid approach to long-term binary neutron-star simulations

Advances in Simulations of Generic Black-Hole Binaries

Formulation Improvements for Critical Collapse Simulations

SEOBNRv5PHM: Next generation of accurate and efficient multipolar precessing-spin effective-one-body waveforms for binary black holes

Hybrid approach to long-term binary neutron-star simulations

Closed-form solutions of spinning, eccentric binary black holes at 1.5 post-Newtonian order

Binary Black Holes in Modified Gravity