The Lagrangian numerical relativity code SPHINCS_BSSN_v1.0

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 (χ = 0.5).