Tidal effects in binary neutron star coalescence
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.
💡 Research Summary
The paper presents a rigorous comparison between long‑duration numerical‑relativity (NR) simulations of binary neutron‑star (BNS) coalescence and the tidal effective‑one‑body (EOB) model that incorporates analytical tidal corrections up to second post‑Newtonian (2PN) order. The NR simulations span roughly ten orbital cycles before the stars touch, providing a high‑fidelity dataset that captures the cumulative influence of tidal interactions on the orbital dynamics and the emitted gravitational‑wave (GW) signal.
On the analytical side, the authors employ the state‑of‑the‑art tidal EOB framework. This model builds on the standard EOB description of binary black‑hole (BBH) systems, adding tidal terms derived from post‑Newtonian theory. The tidal sector includes the leading‑order (LO) contribution as well as the next‑to‑leading 2PN correction, which has been recently derived from high‑order analytical calculations. Importantly, the EOB model retains the same set of free parameters that were calibrated to BBH waveforms; no additional “tidal‑specific” calibration constants are introduced.
The comparison shows that the tidal‑completed EOB model reproduces the NR dynamics and waveforms up to the moment of contact with remarkable accuracy. Phase differences remain within the estimated NR uncertainty (≈0.5 rad), and amplitude discrepancies are below the 1 % level across the entire inspiral. The inclusion of the 2PN tidal term systematically reduces the residuals, confirming that higher‑order analytical information improves the model. However, the statistical analysis indicates that the improvement over the LO tidal term is modest; the current NR error budget does not allow a clear discrimination between the LO and 2PN contributions.
A notable outcome is that the EOB model requires no extra calibration of tidal parameters beyond those already fixed by BBH data. This demonstrates that the present analytical knowledge encoded in the tidal EOB Hamiltonian is already sufficient to capture the essential physics of BNS inspiral up to contact. The authors also investigate the possibility of a strong amplification of tidal effects, which some earlier studies suggested. Their analysis attributes any apparent amplification to numerical artefacts—such as insufficient grid resolution, imperfect initial data, or boundary‑condition effects—rather than to genuine physical mechanisms.
The paper emphasizes the importance of extensive error‑budget studies for NR simulations. It argues that future progress will rely on (i) higher‑resolution and longer‑duration NR runs, (ii) systematic exploration of different equations of state, and (iii) possibly new numerical strategies (e.g., adaptive mesh refinement tailored to tidal regions) to reduce systematic uncertainties.
In conclusion, the authors assert that the tidal‑completed EOB model, as it stands, provides the most advanced and accurate tool for modelling GW signals from BNS inspirals up to the point of contact. This model can be directly employed in GW data‑analysis pipelines to extract astrophysical information—such as neutron‑star radii and the equation of state—from observed BNS events. The work also sets a clear agenda for future improvements: refining NR simulations to tighten error bars and extending analytical tidal knowledge beyond 2PN to further enhance waveform fidelity.