skip to main content

Title: Beyond the linear tide: impact of the non-linear tidal response of neutron stars on gravitational waveforms from binary inspirals

Tidal interactions in coalescing binary neutron stars modify the dynamics of the inspiral and hence imprint a signature on their gravitational wave (GW) signals in the form of an extra phase shift. We need accurate models for the tidal phase shift in order to constrain the supranuclear equation of state from observations. In previous studies, GW waveform models were typically constructed by treating the tide as a linear response to a perturbing tidal field. In this work, we incorporate non-linear corrections due to hydrodynamic three- and four-mode interactions and show how they can improve the accuracy and explanatory power of waveform models. We set up and numerically solve the coupled differential equations for the orbit and the modes and analytically derive solutions of the system’s equilibrium configuration. Our analytical solutions agree well with the numerical ones up to the merger and involve only algebraic relations, allowing for fast phase shift and waveform evaluations for different equations of state over a large parameter space. We find that, at Newtonian order, non-linear fluid effects can enhance the tidal phase shift by $\gtrsim 1\, {\rm radian}$ at a GW frequency of 1000 Hz, corresponding to a $10{{\%}}-20{{\%}}$ correction to the linear theory. The scale of the more » additional phase shift near the merger is consistent with the difference between numerical relativity and theoretical predictions that account only for the linear tide. Non-linear fluid effects are thus important when interpreting the results of numerical relativity and in the construction of waveform models for current and future GW detectors.

« less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 4325-4343
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different BNS 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.CoResimulations 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.

  2. We study the possible occurrence of the hadron-quark phase transition (PT) during the merging of neutron star binaries by hydrodynamical simulations employing a set of temperature dependent hybrid equations of state (EoSs). Following previous work we describe an unambiguous and measurable signature of deconfined quark matter in the gravitational-wave (GW) signal of neutron star binary mergers including equal-mass and unequal-mass systems of different total binary mass. The softening of the EoS by the PT at higher densities, i.e. after merging, leads to a characteristic increase of the dominant postmerger GW frequency f_peak relative to the tidal deformability Lambda inferred during the premerger inspiral phase. Hence, measuring such an increase of the postmerger frequency provides evidence for the presence of a strong PT. If the postmerger frequency and the tidal deformability are compatible with results from purely baryonic EoS models yielding very tight relations between f_peak and Lambda, a strong PT can be excluded up to a certain density. We find tight correlations of f_peak and Lambda with the maximum density during the early postmerger remnant evolution. These GW observables thus inform about the density regime which is probed by the remnant and its GW emission. Exploiting such relations we devisemore »a directly applicable, concrete procedure to constrain the onset density of the QCD PT from future GW measurements. We point out two interesting scenarios: if no indications for a PT are inferred from a GW detection, our procedure yields a lower limit on the onset density of the hadron quark PT. On the contrary, if a merger event reveals evidence for the occurrence of deconfined quark matter, the inferred GW parameters set an upper limit on the PT onset density. (abridged)« less
  3. Abstract

    The effects of strong magnetic fields on the deconfinement phase transition expected to take place in the interior of massive neutron stars are studied in detail for the first time. For hadronic matter, the very general density-dependent relativistic mean field model is employed, while the simple, but effective vector-enhanced bag model is used to study quark matter. Magnetic-field effects are incorporated into the matter equation of state and in the general-relativity solutions, which also satisfy Maxwell’s equations. We find that for large values of magnetic dipole moment, the maximum mass, canonical mass radius, and dimensionless tidal deformability obtained for stars using spherically symmetric Tolman–Oppenheimer–Volkoff (TOV) equations and axisymmetric solutions attained through the LORENE library differ considerably. The deviations depend on the stiffness of the equation of state and on the star mass being analyzed. This points to the fact that, unlike what was assumed previously in the literature, magnetic field thresholds for the approximation of isotropic stars and the acceptable use of TOV equations depend on the matter composition and interactions.

  4. ABSTRACT Compact white dwarf (WD) binaries are important sources for space-based gravitational-wave (GW) observatories, and an increasing number of them are being identified by surveys like Extremely Low Mass (ELM) and Zwicky Transient Facility (ZTF). We study the effects of non-linear dynamical tides in such binaries. We focus on the global three-mode parametric instability and show that it has a much lower threshold energy than the local wave-breaking condition studied previously. By integrating networks of coupled modes, we calculate the tidal dissipation rate as a function of orbital period. We construct phenomenological models that match these numerical results and use them to evaluate the spin and luminosity evolution of a WD binary. While in linear theory the WD’s spin frequency can lock to the orbital frequency, we find that such a lock cannot be maintained when non-linear effects are taken into account. Instead, as the orbit decays, the spin and orbit go in and out of synchronization. Each time they go out of synchronization, there is a brief but significant dip in the tidal heating rate. While most WDs in compact binaries should have luminosities that are similar to previous traveling-wave estimates, a few per cent should be about 10 times dimmermore »because they reside in heating rate dips. This offers a potential explanation for the low luminosity of the CO WD in J0651. Lastly, we consider the impact of tides on the GW signal and show that the Laser Interferometer Space Antenna (LISA) and TianGO can constrain the WD’s moment of inertia to better than $1{{\ \rm per\ cent}}$ for centi-Hz systems.« less

    We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with ourmore »older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

    « less