Abstract:The observation of gravitational waves from compact objects has now become an active part of observational astronomy. For a sound interpretation, one needs to compare such observations against detailed Numerical Relativity simulations, which are essential tools to explore the dynamics and physics of compact binary mergers. To date, essentially all simulation codes that solve the full set of Einstein's equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations, we focus on perfectly symmetric binary systems that are irrotational and have 1.3 M⊙ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle–mesh mapping), and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the Ψ4 formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times 10−3as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity (∼0.5...0.7c) ejecta component of ∼10−4 that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission, and it could also potentially cause a kilonova afterglow.

Numerical Relativity Multimodal Waveforms using Absorbing Boundary Conditions

Fixing the BMS Frame of Numerical Relativity Waveforms

Simulating binary black hole mergers using discontinuous Galerkin methods

Formulating the complete initial boundary value problem in numerical relativity to model black hole echoes

Numerical Tests of Evolution Systems, Gauge Conditions, and Boundary Conditions for 1D Colliding Gravitational Plane Waves

Gravitational Waves from Binary Neutron Star Mergers with a Spectral Equation of State

Improved Moving-Puncture Techniques for Compact Binary Simulations

Accurate waveforms for eccentric, aligned-spin binary black holes: The multipolar effective-one-body model SEOBNRv5EHM

Thinking Outside the Box: Numerical Relativity with Particles

Eccentric binary black hole surrogate models for the gravitational waveform and remnant properties: comparable mass, nonspinning case

Numerical relativity surrogate model with memory effects and post-Newtonian hybridization

Numerical-relativity validation of effective-one-body waveforms in the intermediate-mass-ratio regime

Reflecting boundary conditions in numerical relativity as a model for black hole echoes

A phenomenological gravitational waveform model of binary black holes incorporating horizon fluxes

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

Accuracy limitations of existing numerical relativity waveforms on the data analysis of current and future ground-based detectors

Gravitational waveforms from SpEC simulations : neutron star-neutron star and low-mass black hole-neutron star binaries

Relieving scale disparity in binary black hole simulations

Effective-one-body waveform model for noncircularized, planar, coalescing black hole binaries II: high accuracy by improving logarithmic terms in resummations

Binary Black Holes in Modified Gravity