No Arabic abstract
The combined observation of gravitational and electromagnetic waves from the coalescence of two neutron stars marks the beginning of multi-messenger astronomy with gravitational waves (GWs). The development of accurate gravitational waveform models is a crucial prerequisite to extract information about the properties of the binary system that generated a detected GW signal. In binary neutron star systems (BNS), tidal effects also need to be incorporated in the modeling for an accurate waveform representation. Building on previous work [Phys.Rev.D96 121501], we explore the performance of inspiral-merger waveform models that are obtained by adding a numerical relativity (NR) based approximant for the tidal part of the phasing (NRTidal) to existing models for nonprecessing and precessing binary black hole systems (SEOBNRv4, PhenomD and PhenomPv2), as implemented in the LSC Algorithm Library Suite. The resulting BNS waveforms are compared and contrasted to target waveforms hybridizing NR waveforms, covering the last approx. 10 orbits up to merger and extending through the postmerger phase, with inspiral waveforms calculated from 30Hz obtained with TEOBResumS. The latter is a state-of-the-art effective-one-body waveform model that blends together tidal and spin effects. We probe that the combination of the PN-based self-spin terms and of the NRTidal description is necessary to obtain minimal mismatches (< 0.01) and phase differences (< 1 rad) with respect to the target waveforms. However, we also discuss possible improvements and drawbacks of the NRTidal approximant in its current form, since we find that it tends to overestimate the tidal interaction with respect to the TEOBResumS model during the inspiral.
Spinning neutron stars acquire a quadrupole moment due to their own rotation. This quadratic-in-spin, self-spin effect depends on the equation of state (EOS) and affects the orbital motion and rate of inspiral of neutron star binaries. We incorporate the EOS-dependent self-spin (or monopole-quadrupole) terms in the spin-aligned effective-one-body (EOB) waveform model TEOBResumS at next-to-next-to-leading (NNLO) order, together with other (bilinear, cubic and quartic) nonlinear-in-spin effects (at leading order, LO). The structure of the Hamiltonian of TEOBResumS is such that it already incorporates, in the binary black hole case, the recently computed quartic-in-spin LO term. Using the gauge-invariant characterization of the phasing provided by the function $Q_omega=omega^2/dot{omega}$ of $omega=2pi f$ , where $f$ is the gravitational wave frequency, we study the EOS dependence of the self-spin effects and show that: (i) the next-to-leading order (NLO) and NNLO monopole-quadrupole corrections yield increasingly phase-accelerating effects compared to the corresponding LO contribution; (ii) the standard TaylorF2 post-Newtonian (PN) treatment of NLO (3PN) EOS-dependent self-spin effects makes their action stronger than the corresponding EOB description; (iii) the addition to the standard 3PN TaylorF2 post-Newtonian phasing description of self-spin tail effects at LO allows one to reconcile the self-spin part of the TaylorF2 PN phasing with the corresponding TEOBResumS one up to dimensionless frequencies $Momegasimeq 0.04-0.06$. By generating the inspiral dynamics using the post-adiabatic approximation, incorporated in a new implementation of TEOBResumS, one finds that the computational time needed to obtain a typical waveform (including all multipoles up to $ell=8$) from 10 Hz is of the order of 0.4 sec.
Gravitational wave (GW) astronomy has consolidated its role as a new observational window to reveal the properties of compact binaries in the Universe. In particular, the discovery of the first binary neutron star coalescence, GW170817, led to a number of scientific breakthroughs as the possibility to place constraints on the equation of state of cold matter at supranuclear densities. These constraints and all scientific results based on them require accurate models describing the GW signal to extract the source properties from the measured signal. In this article, we study potential systematic biases during the extraction of source parameters using different descriptions for both, the point-particle dynamics and tidal effects. We find that for the considered cases the mass and spin recovery show almost no systematic bias with respect to the chosen waveform model. However, the extracted tidal effects can be strongly biased, where we find generally that Post-Newtonian approximants predict neutron stars with larger deformability and radii than numerical relativity tuned models. Noteworthy, an increase in the Post-Newtonian order in the tidal phasing does not lead to a monotonic change in the estimated properties. We find that for a signal with strength similar to GW170817, but observed with design sensitivity, the estimated tidal parameters can differ by more than a factor of two depending on the employed tidal description of the waveform approximant. This shows the current need for the development of better waveform models to extract reliably the source properties from upcoming GW detections.
Gravitational waves radiated by the coalescence of compact-object binaries containing a neutron star and a black hole are one of the most interesting sources for the ground-based gravitational-wave observatories Advanced LIGO and Advanced Virgo. Advanced LIGO will be sensitive to the inspiral of a $1.4, M_odot$ neutron star into a $10,M_odot$ black hole to a maximum distance of $sim 900$ Mpc. Achieving this sensitivity and extracting the physics imprinted in observed signals requires accurate modeling of the binary to construct template waveforms. In a NSBH binary, the black hole may have significant angular momentum (spin), which affects the phase evolution of the emitted gravitational waves. We investigate the ability of post-Newtonian (PN) templates to model the gravitational waves emitted during the inspiral phase of NSBH binaries. We restrict the black holes spin to be aligned with the orbital angular momentum and compare several approximants. We examine restricted amplitude waveforms that are accurate to 3.5PN order in the orbital dynamics and complete to 2.5PN order in the spin dynamics. We also consider PN waveforms with the recently derived 3.5PN spin-orbit and 3PN spin-orbit tail corrections. We compare these approximants to the effective-one-body model. For all these models, large disagreements start at low to moderate black hole spins, particularly for binaries where the spin is anti-aligned with the orbital angular momentum. We show that this divergence begins in the early inspiral at $v sim 0.2$ for $chi_{BH} sim 0.4$. PN spin corrections beyond those currently known will be required for optimal detection searches and to measure the parameters of neutron star--black hole binaries. While this complicates searches, the strong dependence of the gravitational-wave signal on the spin dynamics will make it possible to extract significant astrophysical information.
We study the effect of superfluidity on the tidal response of a neutron star in a general relativistic framework. In this work, we take a dual-layer approach where the superfluid matter is confined in the core of the star. Then, the superfluid core is encapsulated with an envelope of ordinary matter fluid which acts effectively as the low-density crustal region of the star. In the core, the matter content is described by a two-fluid model where only the neutrons are taken as superfluid and the other fluid consists of protons and electrons making it charge neutral. We calculate the values of various tidal love numbers of a neutron star and discuss how they are affected due to the presence of entrainment between the two fluids in the core. We also emphasize that more than one tidal parameter is necessary to probe superfluidity with the gravitational wave from the binary inspiral.
We reanalyze gravitational waves from binary-neutron-star mergers GW170817 and GW190425 using a numerical-relativity (NR) calibrated waveform model, the TF2+_Kyoto model, which includes nonlinear tidal terms. For GW170817, by imposing a uniform prior on the binary tidal deformability $tilde{Lambda}$, the symmetric $90%$ credible interval of $tilde{Lambda}$ is estimated to be $481^{+436}_{-359}$ and $402^{+465}_{-279}$ for the case of $f_mathrm{max}=1000$ and $2048~mathrm{Hz}$, respectively, where $f_mathrm{max}$ is the maximum frequency in the analysis. We also reanalyze the event with other waveform models: two post-Newtonian waveform models (TF2_PNTidal and TF2+_PNTidal), the TF2+_NRTidal model that is another NR calibrated waveform model, and its upgrade, the TF2+_NRTidalv2 model. While estimates of parameters other than $tilde{Lambda}$ are broadly consistent among various waveform models, our results indicate that estimates of $tilde{Lambda}$ depend on waveform models. However, the difference is smaller than the statistical error. For GW190425, we can only obtain little information on the binary tidal deformability. The systematic difference among the NR calibrated waveform models will become significant to measure $tilde{Lambda}$ as the number of detectors and events increase and sensitivities of detectors are improved.