The construction of accurate and consistent initial data for various binary parameters is a critical ingredient for numerical relativity simulations of the compact binary coalescence. In this article, we present an upgrade of the pseudospectral SGRID code, which enables us to access even larger regions of the binary neutron star parameter space. As a proof of principle, we present a selected set of first simulations based on initial configurations computed with the new code version. In particular, we simulate two millisecond pulsars close to their breakup spin, highly compact neutron stars with masses at about 98% of the maximum supported mass of the employed equation of state, and an unequal mass systems with mass ratios even outside the range predicted by population synthesis models (q = 2.03). The discussed code extension will help us to simulate previously unexplored binary configurations. This is a necessary step to construct and test new gravitational wave approximants and to interpret upcoming binary neutron star merger observations. When we construct initial data, one has to specify various parameters, such as a rotation parameter for each star. Some of these parameters do not have direct physical meaning, which makes comparisons with other methods or models difficult. To facilitate this, we introduce simple estimates for the initial spin, momentum, mass, and center of mass of each individual star.

We present the second data release of gravitational waveforms from binary neutron star merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different binary neutron star configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to l=m=4. They span a significant portion of the mass, mass-ratio,spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

We present the second data release of gravitational waveforms from binary neutron star merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different binary neutron star configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to = m = 4. They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data

The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numericalrelativity simulations for four different physical setups with mass ratios q = 1.25, 1.50, 1.75, 2.00, and total gravitational mass M = 2.7M . Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately 2nd order converging results for the dominant (2,2), but also subdominant (2,1), (3,3), (4,4) modes, while, generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current NRTidal model to describe tidal effects is a valid description for high-mass ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs.

No abstract