We continue our study of the binary neutron star parameter space by investigating the effect of the spin orientation on the dynamics, gravitational wave emission, and mass ejection during the binary neutron star coalescence. We simulate seven different configurations using multiple resolutions to allow a reasonable error assessment. Due to the particular choice of the setups, five configurations show precession effects, from which two show a precession (`wobbling') of the orbital plane, while three show a `bobbing' motion, i.e., the orbital angular momentum does not precess, while the orbital plane moves along the orbital angular momentum axis. Considering the ejection of mass, we find that precessing systems can have an anisotropic mass ejection, which could lead to a final remnant kick of about $\sim 40 \rm km/s$ for the studied systems. Furthermore, for the chosen configurations, anti-aligned spins lead to larger mass ejecta than aligned spins, so that brighter electromagnetic counterparts could be expected for these configurations. Finally, we compare our simulations with the precessing, tidal waveform approximant IMRPhenomPv2_NRTidalv2 and find good agreement between the approximant and our numerical relativity waveforms with phase differences below 1.2 rad accumulated over the last $\sim$ 16 gravitational wave cycles.

Accurate and fast gravitational waveform (GW) models are essential to extract information about the properties of compact binary systems that generate GWs. Building on previous work, we present an extension of the NRTidal model for binary neutron star (BNS) waveforms. The upgrades are: (i) a new closed-form expression for the tidal contribution to the GW phase which includes further analytical knowledge and is calibrated to more accurate numerical relativity data than previously available; (ii) a tidal correction to the GW amplitude; (iii) an extension of the spin-sector incorporating equation-of-state-dependent finite size effects at quadrupolar and octupolar order; these appear in the spin-spin tail terms and cubic-in-spin terms, both at 3.5PN. We add the new description to the precessing binary black hole waveform model IMRPhenomPv2 to obtain a frequency-domain precessing binary neutron star model. In addition, we extend the SEOBNRv4_ROM and IMRPhenomD aligned-spin binary black hole waveform models with the improved tidal phase corrections. Focusing on the new IMRPhenomPv2_NRTidalv2 approximant, we test the model by comparing with numerical relativity waveforms as well as hybrid waveforms combining tidal effective-one-body and numerical relativity data. We also check consistency against a tidal effective-one-body model across large regions of the BNS parameter space.

As current gravitational wave (GW) detectors increase in sensitivity, and particularly as new instruments are being planned, there is the possibility that ground-based GW detectors will observe GWs from highly eccentric neutron star binaries. We present the first detailed study of highly eccentric BNS systems with full (3+1)D numerical relativity simulations using consistent initial conditions, i.e., setups which are in agreement with the Einstein equations and with the equations of general relativistic hydrodynamics in equilibrium. Overall, our simulations cover two different equations of state (EOSs), two different spin configurations, and three to four different initial eccentricities for each pairing of EOS and spin. We extract from the simulated waveforms the frequency of the f-mode oscillations induced during close encounters before the merger of the two stars. The extracted frequency is in good agreement with f-mode oscillations of individual stars for the irrotational cases, which allows an independent measure of the supranuclear equation of state not accessible for binaries on quasi-circular orbits. The energy stored in these f-mode oscillations can be as large as $10^{-3}M_\odot \sim 10^{51}$ erg, even with a soft EOS. In order to estimate the stored energy, we also examine the effects of mode mixing due to the stars' offset from the origin on the f-mode contribution to the GW signal. While in general (eccentric) neutron star mergers produce bright electromagnetic counterparts, we find that the luminosity decreases when the eccentricity becomes too large, due to a decrease of the ejecta mass. Finally, the use of consistent initial configurations also allows us to produce high-quality waveforms for different eccentricities which can be used as a testbed for waveform model development of highly eccentric binary neutron star systems.

We construct closed-form gravitational waveforms (GWs) with tidal effects for the coalescence and merger of binary neutron stars. The method relies on a new set of eccentricity-reduced and high-resolution numerical relativity (NR) simulations and is composed of three steps. First, tidal contributions to the GW phase are extracted from the time-domain NR data. Second, those contributions are employed to fix high-order coefficients in an effective and resummed post-Newtonian expression. Third, frequency-domain tidal approximants are built using the stationary phase approximation. Our tidal approximants are valid from the low frequencies to the strong-field regime and up to merger. They can be analytically added to any binary black hole GW model to obtain a binary neutron star waveform, either in the time or in the frequency domain. This work provides simple, flexible, and accurate models ready to be used in both searches and parameter estimation of binary neutron star events.

We present the Computational Relativity CoRe collaboration's public database of gravitational waveforms from binary neutron star mergers. The database currently contains 367 waveforms from numerical simulations that are consistent with general relativity and that employ constraint satisfying initial data in hydrodynamical equilibrium. It spans 164 physically distinct configuration with different binary parameters (total binary mass, mass-ratio, initial separation, eccentricity, and stars' spins) and simulated physics. Waveforms computed at multiple grid resolutions and extraction radii are provided for controlling numerical uncertainties. We also release an exemplary set of 18 hybrid waveforms constructed with a state-of-art effective-one-body model spanning the frequency band of advanced gravitational-wave detectors. We outline present and future applications of the database to gravitational-wave astronomy.

The recent detection of gravitational waves and electromagnetic counterparts emitted during and after the collision of two neutron stars marks a breakthrough in the field of multi-messenger astronomy. Numerical relativity simulations are the only tool to describe the binary's merger dynamics in the regime when speeds are largest and gravity is strongest. In this work we report state-of-the-art binary neutron star simulations for irrotational (non-spinning) and spinning configurations. The main use of these simulations is to model the gravitational-wave signal. Key numerical requirements are the understanding of the convergence properties of the numerical data and a detailed error budget. The simulations have been performed on different HPC clusters, they use multiple grid resolutions, and are based on eccentricity reduced quasi-circular initial data. We obtain convergent waveforms with phase errors of 0.5-1.5 rad accumulated over approximately 12 orbits to merger. The waveforms have been used for the construction of a phenomenological waveform model which has been applied for the analysis of the recent binary neutron star detection. Additionally, we show that the data can also be used to test other state-of-the-art semi-analytical waveform models.

We present the first set of numerical relativity simulations of binary neutron mergers that include spin precession effects and are evolved with multiple resolutions. Our simulations employ consistent initial data in general relativity with different spin configurations and dimensionless spin magnitudes $\sim 0.1$. They start at a gravitational-wave frequency of $\sim392$~Hz and cover more than $1$ precession period and about 15 orbits up to merger. We discuss the spin precession dynamics by analyzing coordinate trajectories, quasi-local spin measurements, and energetics, by comparing spin aligned, antialigned, and irrotational configurations. Gravitational waveforms from different spin configuration are compared by calculating the mismatch between pairs of waveforms in the late inspiral. We find that precession effects are not distinguishable from nonprecessing configurations with aligned spins for approximately face-on binaries, while the latter are distinguishable from a nonspinning configurations. Spin precession effects are instead clearly visible for approximately edge-on binaries. For the parameters considered here, precession does not significantly affect the characteristic postmerger gravitational-wave frequencies nor the mass ejection. Our results pave the way for the modeling of spin precession effects in the gravitational waveform from binary neutron star events.

We present new (3+1) dimensional numerical relativity simulations of the binary neutron star (BNS) mergers that take into account the NS spins. We consider different spin configurations, aligned or antialigned to the orbital angular momentum, for equal and unequal mass BNS and for two equations of state. All the simulations employ quasiequilibrium circular initial data in the constant rotational velocity approach, i.e. they are consistent with Einstein equations and in hydrodynamical equilibrium. We study the NS rotation effect on the energetics, the gravitational waves (GWs) and on the possible electromagnetic (EM) emission associated to dynamical mass ejecta. For dimensionless spin magnitudes of $\chi\sim0.1$ we find that spin-orbit interactions and also spin-induced-quadrupole deformations affect the late-inspiral-merger dynamics. The latter is, however, dominated by finite-size effects. Spin (tidal) effects contribute to GW phase differences up to 5 (20) radians accumulated during the last eight orbits to merger. Similarly, after merger the collapse time of the remnant and the GW spectrogram are affected by the NSs rotation. Spin effects in dynamical ejecta are clearly observed in unequal mass systems in which mass ejection originates from the tidal tail of the companion. Consequently kilonovae and other EM counterparts are affected by spins. We find that spin aligned to the orbital angular momentum leads to brighter EM counterparts than antialigned spin with luminosities up to a factor of two higher.

We present new (3+1)D numerical relativity simulations of the binary neutron star (BNS) merger and postmerger phase. We focus on a previously inaccessible region of the binary parameter space spanning the binary's mass-ratio $q\sim1.00-1.75$ for different total masses and equations of state, and up to $q\sim2$ for a stiff BNS system. We study the mass-ratio effect on the gravitational waves (GWs) and on the possible electromagnetic emission associated to dynamical mass ejecta. We compute waveforms, spectra, and spectrograms of the GW strain including all the multipoles up to $l=4$. The mass-ratio has a specific imprint on the GW multipoles in the late-inspiral-merger signal, and it affects qualitatively the spectra of the merger remnant. The multipole effect is also studied by considering the dependency of the GW spectrograms on the source's sky location. Unequal mass BNSs produce more ejecta than equal mass systems with ejecta masses and kinetic energies depending almost linearly on $q$. We estimate luminosity peaks and light curves of macronovae events associated to the mergers using a simple approach. For $q\sim2$ the luminosity peak is delayed for several days and can be up to four times larger than for the $q=1$ cases. The macronova emission associated with the $q\sim2$ BNS is more persistent in time and could be observed for weeks instead of few days ($q=1$) in the near infrared. Finally, we estimate the flux of possible radio flares produced by the interaction of relativistic outflows with the surrounding medium. Also in this case a large $q$ can significantly enhance the emission and delay the peak luminosity. Overall, our results indicate that BNS merger with large mass ratio have EM signatures distinct from the equal mass case and more similar to black hole - neutron star binaries.

Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spins, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $\sim10^{-2}M_\odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged]

In general neutron stars in binaries are spinning. Recently, a new quasi-equilibrium approximation that includes a rotational velocity piece for each star has been proposed to describe binary neutron stars with arbitrary rotation states in quasi-circular orbits. We have implemented this approximation numerically for the first time, to generate initial data for neutron star binaries with spin. If we choose the rotational velocity piece such that it equals the Newtonian rigid rotation law, we obtain stars with fluid 4-velocities that have expansion and shear of approximately zero, as one would expect for quasi-equilibrium configurations. We also use the new approach to construct and study initial data sequences for irrotational, corotating and fixed rotation binaries.

Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipse's inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero.

Simulations of binary neutron stars have seen great advances in terms of physical detail and numerical quality. However, the spin of the neutron stars, one of the simplest global parameters of binaries, remains mostly unstudied. We present the first, fully nonlinear general relativistic dynamical evolutions of the last three orbits for constraint satisfying initial data of spinning neutron star binaries, with astrophysically realistic spins aligned and anti-aligned to the orbital angular momentum. The initial data is computed with the constant rotational velocity approach. The dynamics of the systems is analyzed in terms of gauge-invariant binding energy vs. orbital angular momentum curves. By comparing to a binary black hole configuration we can estimate the different tidal and spin contributions to the binding energy for the first time. First results on the gravitational wave forms are presented. The phase evolution during the orbital motion is significantly affected by spin-orbit interactions, leading to delayed or early mergers. Furthermore, a frequency shift in the main emission mode of the hyper massive neutron star is observed. Our results suggest that a detailed modeling of merger waveforms requires the inclusion of spin, even for the moderate magnitudes observed in binary neutron star systems.

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]

In general neutron stars in binaries are spinning. Due to the existence of millisecond pulsars we know that these spins can be substantial. We argue that spins with periods on the order a few dozen milliseconds could influence the late inspiral and merger dynamics. Thus numerical simulations of the last few orbits and the merger should start from initial conditions that allow for arbitrary spins. We discuss quasi-equilibrium approximations one can make in the construction of binary neutron star initial data with spins. Using these approximations we are able to derive two new matter equations. As in the case of irrotational neutron star binaries one of these equations is algebraic and the other elliptic. If these new matter equations are solved together with the equations for the metric variables following the Wilson-Mathews or conformal thin sandwich approach one can construct neutron star initial data. The spin of each star is described by a rotational velocity that can be chosen freely so that one can create stars in arbitrary rotation states. Our new matter equations reduce to the well known limits of both corotating and irrotational neutron star binaries.

We consider black holes resulting from binary black hole mergers. By fitting to numerical results we construct analytic formulas that predict the mass and spin of the final black hole. Our formulas are valid for arbitrary initial spins and mass ratios and agree well with available numerical simulations. We use our spin formula in the context of two common merger scenarios for supermassive galactic black holes. We consider the case of isotropically distributed initial spin orientations (when no surrounding matter is present) and also the case when matter closely aligns the spins with the orbital angular momentum. The spin magnitude of black holes resulting from successive generations of mergers (with symmetric mass ratio $\eta$) has a mean of $1.73\eta + 0.28$ in the isotropic case and 0.94 for the closely aligned case.

We present results from several simulations of equal mass black holes with spin. The spin magnitudes are $S/m^2=0.8$ in all cases, but we vary the spin orientations arbitrarily, in and outside the orbital plane. We find that in all but one case the final merged black hole acquires a kick of more than 1000 km/s, indicating that kicks of this magnitude are likely to be generic and should be expected for mergers with general spin orientations. The maximum kick velocity we find is 2500 km/s and occurs for initial spins which are anti-aligned in the initial orbital plane.