Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger remnants. However, most general relativistic studies have been performed for nearly circular binaries, with the exception of head-on collisions. We present the first numerical relativity investigation of mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and circular. In addition to gravitational waves generated by the orbital motion, we find that the signal also contains a strong component due to stellar oscillations (f-modes) induced by tidal forces, extending a classical result for Newtonian binaries. The merger can lead to rather massive disks on the order of 10% of the total initial mass.

We present numerical relativity simulations of nine-orbit equal-mass binary neutron star covering the quasicircular late inspiral and merger. The extracted gravitational waveforms are analyzed for convergence and accuracy. Second order convergence is observed up to contact, i.e. about 3-4 cycles to merger; error estimates can be made up to this point. The uncertainties on the phase and the amplitude are dominated by truncation errors and can be minimized to 0.13 rad and less then 1%, respectively, by using several simulations and extrapolating in resolution. In the latter case finite-radius extraction uncertainties become a source of error of the same order and have to be taken into account. The waveforms are tested against accuracy standards for data analysis. The uncertainties on the waveforms are such that accuracy standards are generically not met for signal-to-noise ratios relevant for detection, except for some best cases using extrapolation from several runs. A detailed analysis of the errors is thus imperative for the use of numerical relativity waveforms from binary neutron stars in quantitative studies. The waveforms are compared with the post-Newtonian Taylor T4 approximants both for point-particle and including the analytically known tidal corrections. The T4 approximants accumulate significant phase differences of 2 rad at contact and 4 rad at merger, underestimating the influence of finite size effects. Tidal signatures in the waveforms are thus important at least during the last six orbits of the merger process.

We compare dynamics and waveforms from binary neutron star coalescence as computed by new long-term ($\sim 10 $ orbits) numerical relativity simulations and by the tidal effective-one-body (EOB) model including analytical tidal corrections up to second post-Newtonian order (2PN). The current analytical knowledge encoded in the tidal EOB model is found to be sufficient to reproduce the numerical data up to contact and within their uncertainties. Remarkably, no calibration of any tidal EOB free parameters is required, beside those already fitted to binary black holes data. The inclusion of 2PN tidal corrections minimizes the differences with the numerical data, but it is not possible to significantly distinguish them from the leading-order tidal contribution. The presence of a relevant amplification of tidal effects is likely to be excluded, although it can appear as a consequence of numerical inaccuracies. We conclude that the tidally-completed effective-one-body model provides nowadays the most advanced and accurate tool for modelling gravitational waveforms from binary neutron star inspiral up to contact. This work also points out the importance of extensive tests to assess the uncertainties of the numerical data, and the potential need of new numerical strategies to perform accurate simulations.

We present a new numerical relativity code designed for simulations of compact binaries involving matter. The code is an upgrade of the BAM code to include general relativistic hydrodynamics and implements state-of-the-art high-resolution-shock-capturing schemes on a hierarchy of mesh refined Cartesian grids with moving boxes. We test and validate the code in a series of standard experiments involving single neutron star spacetimes. We present test evolutions of quasi-equilibrium equal-mass irrotational binary neutron star configurations in quasi-circular orbits which describe the late inspiral to merger phases. Neutron star matter is modeled as a zero-temperature fluid; thermal effects can be included by means of a simple ideal-gas prescription. We analyze the impact that the use of different values of damping parameter in the Gamma-driver shift condition has on the dynamics of the system. The use of different reconstruction schemes and their impact in the post-merger dynamics is investigated. We compute and characterize the gravitational radiation emitted by the system. Self-convergence of the waves is tested, and we consistently estimate error-bars on the numerically generated waveforms in the inspiral phase.