The Lagrangian Numerical Relativity code SPHINCS_BSSN_v1.0

Stephan Rosswog,Francesco Torsello,Peter Diener
DOI: https://doi.org/10.3389/fams.2023.1236586
2023-10-07
Abstract:We present version 1.0 of our Lagrangian Numerical Relativity code SPHINCS_BSSN. This code evolves the full set of Einstein equations, but contrary to other Numerical Relativity codes, it evolves the matter fluid via Lagrangian particles in the framework of a high-accuracy version of Smooth Particle Hydrodynamics (SPH). The major new elements introduced here are: i) a new method to map the stress--energy tensor (known at the particles) to the spacetime mesh, based on a local regression estimate; ii) additional measures that ensure the robust evolution of a neutron star through its collapse to a black hole; and iii) further refinements in how we place the SPH particles for our initial data. The latter are implemented in our code SPHINCS_ID which now, in addition to LORENE, can also couple to initial data produced by the initial data library FUKA. We discuss several simulations of neutron star mergers performed with SPHINCS_BSSN_v1.0, including irrotational cases with and without prompt collapse and a system where only one of the stars has a large spin ($\chi = 0.5)$.
General Relativity and Quantum Cosmology
What problem does this paper attempt to address?
The paper attempts to address the development of a Lagrangian numerical relativity code named SPHINCS BSSN v1.0, which is used to simulate the general relativistic effects and hydrodynamics in neutron star mergers. Specifically, the paper introduces the following points: 1. **New Mapping Method**: A new method is proposed to map the stress-energy tensor on particles to the spacetime grid, based on Local Regression Estimate (LRE). 2. **Robustness of Neutron Star Collapse**: Additional measures are introduced to ensure stable evolution of neutron stars when collapsing into black holes. 3. **Improvement of Initial Data**: Further optimization of the initial data setup for SPH particles is conducted, which can be coupled with the LORENE and FUKA initial data libraries. With these improvements, SPHINCS BSSN v1.0 is capable of simulating various scenarios of neutron star mergers with high precision, including prompt collapse in non-rotating cases, partially rotating neutron star systems, and more. These simulations are significant for understanding gravitational wave signals, the equation of state of neutron star matter, and electromagnetic radiation in multi-messenger astronomy.