Inspiralling Compact Binaries Research Articles

The quadrupole moment of compact binaries to the fourth post-Newtonian order: II. Dimensional regularization and renormalization

The regularization and renormalization of the radiative mass-type quadrupole moment of inspiralling compact binaries (without spins) is investigated at the fourth post-Newtonian (4PN) approximation of general relativity. As clear from the conservative 4PN equations of motion, a dimensional regularization has to be implemented in order to properly treat the non-linear interactions experienced by gravitational waves during their propagation towards future null infinity. By implementing such procedure, we show that the poles coming from the source moment (computed in a companion paper) are exactly cancelled in the radiative moment, as expected for a physical quantity. We thus define and obtain a ‘renormalized’ source quadrupole, three-dimensional by nature, which is an important step towards the computation of the gravitational-wave flux with 4PN accuracy. Furthermore, we explicitly prove the equivalence between the dimensional regularization and the previously used Hadamard partie finie scheme up to the 3PN order.

Linear momentum flux from inspiralling compact binaries in quasielliptical orbits at 2.5 post-Newtonian order

Emission of anisotropic gravitational radiation from compact binary system leads to a flux of linear momentum. This results in the recoil of the system. We investigate the rate of loss of Linear momentum flux in the far zone of the source using various mass type and current type multipole moments for inspiralling compact binary mergers in quasi-elliptical orbits at 2.5 Post Newtonian order. We compute the linear momentum flux accurate up to $\mathcal{O}(e_t)$ in harmonic coordinate. A 2.5 Post Newtonian Quasi-Keplarian representation of the parametric solution to the Post Newtonian equation of motion for the compact binary system has been adopted here. We also provide a closed-form expression for the accumulated linear momentum from the remote past through the binary evolution.

The mass quadrupole moment of compact binary systems at the fourth post-Newtonian order

The mass-type quadrupole moment of inspiralling compact binaries (without spins) is computed at the fourth post-Newtonian (4PN) approximation of general relativity. The multipole moments are defined by matching between the field in the exterior zone of the matter system and the PN field in the near zone, following the multipolar-post-Minkowskian (MPM)-PN formalism. The matching implies a specific regularization for handling infra-red (IR) divergences of the multipole moments at infinity, based on the Hadamard finite part procedure. On the other hand, the calculation entails ultra-violet (UV) divergences due to the modelling of compact objects by delta-functions, that are treated with dimensional regularization (DR). In future work we intend to systematically study the IR divergences by means of dimensional regularization as well. Our result constitutes an important step in the goal of obtaining the gravitational wave templates of inspiralling compact binary systems with 4PN/4.5PN accuracy.

Post-Newtonian Lagrangian of an N -body system with arbitrary mass and spin multipoles

The present paper derives the post-Newtonian Lagrangian of translational motion of N arbitrary-structured bodies with all mass and spin multipoles in a scalar-tensor theory of gravity. The multipoles depend on time and evolve in accordance with their own dynamic equations of motion. The Lagrangian is retrieved from the post-Newtonian equations of motion by solving the inverse problem of the Lagrangian mechanics and generalizes a well-known Lagrangian of pole-dipole-quadrupole massive particles to the particles of higher multipolarity. Analytic treatment of the higher-order multipole contributions is important for more rigorous computation of gravitational waveform of inspiralling compact binaries at the latest stage of their orbital evolution before merger when tidal and rotational deformations of stars are no longer small and rapidly change in time. The Lagrangian of an N-body system with arbitrary mass and spin multipoles is instrumental for formulation of the post-Newtonian conservation laws of energy, momenta and the integrals of the center of mass.

Second post-Newtonian order radiative dynamics of inspiralling compact binaries in the effective field theory approach

We use the Effective Field Theory (EFT) framework to compute the mass quadrupole moment, the equation of motion, and the power loss of inspiralling compact binaries at the second order in the Post-Newtonian (PN) approximation. We present expressions for the stress-energy pseudo-tensor components of the binary system in higher PN orders. The 2PN correction to the mass quadrupole moment as well as to the acceleration computed in the linearized harmonic gauge presented here are the ingredients needed for the calculation of the next-to-next-to leading order radiation reaction force, which will be presented elsewhere. While this paper reproduces known results, it supplies the building blocks necessary for future higher order calculations in the EFT methodology.

Gravitational waves from inspiralling compact binaries in conformal gravity

We investigate the production of gravitational waves during the inspiral of compact binaries close to their merger in the context of a conformal gravity model. The model incorporates five massive polarization degrees of freedom, besides the two massless gravitational wave polarizations of general relativity. For small graviton mass, we find that the amplitude of the gravitational waves is strongly suppressed as compared to general relativity and decreases as coalescence is approached, which contradicts the observational fact. We conclude that this model with small graviton mass, including a regime that can explain galaxy rotation curves without dark matter, cannot describe the observed gravitational wave events. For a large graviton mass, the modifications to the waveform, compared to the one from general relativity, are negligible on the relevant distance scales and hence a conformal gravity model with a large graviton mass is in agreement with LIGO/VIRGO observations and leads to chirp mass and distance estimates that agree with those from general relativity.

The eccentric behavior of inspiralling compact binaries

Even if there is no gauge invariant definition of eccentricity, it has an important impact on the observed gravitational wave signal of such systems, generating power in all possible harmonics of the orbital period. We here clarify the possible discrepancies between different eccentricity parameters used to describe the orbital dynamics of binary systems across different approximations, specifically the post-Newtonian approximation, the self-force approximation, and numerical relativity. To this end, we highlight disparities between the typically used orbit averaged method of evolving binary systems under radiation reaction, and more direct techniques of solving the two-body problem in post-Newtonian theory. We show, both numerically and analytically, that the orbit averaged method breaks down in the late inspiral, failing to capture a strong secular growth in the Keplerian eccentricity parameter and producing a orbital de-phasing relative to direct integration of the two-body equations of motion. We show that the secular growth and de-phasing affect the observed gravitational wave signal, which could bias how accurately we may recover parameters for systems with signal-to-noise ratios . We further develop a frequency domain post-adiabatic waveform model to capture these effects, and study the precision to which we may estimate parameters with this model through a Fisher information matrix analysis.

Precise peculiar velocities from gravitational waves accompanied by electromagnetic signals and cosmological applications

Peculiar velocities are a precious tool to study the large-scale distribution of matter in the local universe and test cosmological models. However, present measurements of peculiar velocities are based on empirical distance indicators, which introduce large error bars. Here we present a new method to measure the peculiar velocities, by directly estimating luminosity distances through waveform signals from inspiralling compact binaries and measuring redshifts from electromagnetic (EM) counterparts. In the future, with the distance uncertainty of GW events reducing to $0.1$ per cent by future GW detectors, the uncertainty of the peculiar velocity can be reduced to $10$ km/s at 100 mega parsecs. We find that dozens of GW events with EM counterparts can provide a Hubble constant $H_0$ uncertainty of $0.5\%$ and the growth rate of structure with a $0.6\%$ precision in the third-generation ground-base GW detectors, which can reconcile the $H_0$ tension and determine the origins for cosmic accelerated expansion.

Testing Brans-Dicke gravity using the Einstein telescope

Gravitational radiation is an excellent field for testing theories of gravity in strong gravitational fields. The current observations on the gravitational-wave (GW) bursts by LIGO have already placed various constraints on the alternative theories of gravity. In this paper, we investigate the possible bounds which could be placed on the Brans-Dicke gravity using GW detection from inspiralling compact binaries with the proposed Einstein Telescope, a third-generation GW detector. We first calculate in details the waveforms of gravitational radiation in the lowest post-Newtonian approximation, including the tensor and scalar fields, which can be divided into the three polarization modes, i.e. "plus mode", "cross mode" and "breathing mode". Applying the stationary phase approximation, we obtain their Fourier transforms, and derive the correction terms in amplitude, phase and polarization of GWs, relative to the corresponding results in General Relativity. Imposing the noise level of Einstein Telescope, we find that the GW detection from inspiralling compact binaries, composed of a neutron star and a black hole, can place stringent constraints on the Brans-Dicke gravity. The bound on the coupling constant $\omega_{\rm BD}$ depends on the mass, sky-position, orbital angle, polarization angle, luminosity distance, redshift distribution and total observed number $N_{\rm GW}$ of the binary systems. Taking into account all the burst events up to redshift $z=5$, we find that the bound could be $\omega_{\rm BD}\gtrsim 10^{6}\times(N_{\rm GW}/10^4)^{1/2}$. Even for the conservative estimation with $10^{4}$ observed events, the bound is still more than one order tighter than the current limit from Solar System experiments. So, we conclude that Einstein Telescope will provide a powerful platform to test alternative theories of gravity.

Gravitational-wave phasing for low-eccentricity inspiralling compact binaries to 3PN order

[abridged] Although gravitational radiation causes inspiralling compact binaries to circularize, a variety of astrophysical scenarios suggest that binaries might have small but nonnegligible orbital eccentricities when they enter the low-frequency bands of ground and space-based gravitational-wave detectors. If not accounted for, even a small orbital eccentricity can cause a potentially significant systematic error in the mass parameters of an inspiralling binary. Gravitational-wave search templates typically rely on the quasi-circular approximation, which provides relatively simple expressions for the gravitational-wave phase to 3.5 post-Newtonian (PN) order. The quasi-Keplerian formalism provides an elegant but complex description of the post-Newtonian corrections to the orbits and waveforms of inspiralling binaries with any eccentricity. Here we specialize the quasi-Keplerian formalism to binaries with low eccentricity. In this limit the non-periodic contribution to the gravitational-wave phasing can be expressed explicitly as simple functions of frequency or time, with little additional complexity beyond the well-known formulas for circular binaries. These eccentric phase corrections are computed to 3PN order and to leading order in the eccentricity for the standard PN approximants. For a variety of systems these eccentricity corrections cause significant corrections to the number of gravitational wave cycles that sweep through a detector's frequency band. This is evaluated using several measures, including a modification of the useful cycles. We also evaluate the role of periodic terms that enter the phasing and discuss how they can be incorporated into some of the PN approximants. While the eccentric extension of the PN approximants is our main objective, this work collects a variety of results that may be of interest to others modeling eccentric relativistic binaries.

Estimating effective higher order terms in the post-Newtonian binding energy and gravitational-wave flux: Nonspinning compact binary inspiral

In the adiabatic post-Newtonian (PN) approximation, the phase evolution of gravitational waves (GWs) from inspiralling compact binaries in quasicircular orbits is computed by equating the change in binding energy with the GW flux. This energy balance equation can be solved in different ways, which result in multiple approximants of the PN waveforms. Due to the poor convergence of the PN expansion, these approximants tend to differ from each other during the late inspiral. Which of these approximants should be chosen as templates for detection and parameter estimation of GWs from inspiraling compact binaries is not obvious. In this paper, we present estimates of the effective higher order (beyond the currently available 4PN and 3.5PN) non-spinning terms in the PN expansion of the binding energy and the GW flux that minimize the difference of multiple PN approximants (TaylorT1, TaylorT2, TaylorT4, TaylorF2) with effective one body waveforms calibrated to numerical relativity (EOBNR). We show that PN approximants constructed using the effective higher order terms show significantly better agreement (as compared to 3.5PN) with the inspiral part of the EOBNR. For non-spinning binaries with component masses $m_{1,2} \in [1.4 M_\odot, 15 M_\odot]$, most of the approximants have a match (faithfulness) of better than 99% with both EOBNR and each other.

Effectual template bank for the detection of gravitational waves from inspiralling compact binaries with generic spins

We report the construction of a three-dimensional template bank for the search for gravitational waves from inspiralling binaries consisting of spinning compact objects. The parameter space consists of two dimensions describing the mass parameters and one "reduced-spin" parameter, which describes the secular (non-precessing) spin effects in the waveform. The template placement is based on an efficient stochastic algorithm and makes use of the semi-analytical computation of a metric in the parameter space. We demonstrate that for "low-mass" ($m_1 + m_2 \lesssim 12\,M_\odot$) binaries, this template bank achieves effective fitting factors $\sim0.92$--$0.99$ towards signals from generic spinning binaries in the advanced detector era over the entire parameter space of interest (including binary neutron stars, binary black holes, and black hole-neutron star binaries). This provides a powerful and viable method for searching for gravitational waves from generic spinning low-mass compact binaries. Under the assumption that spin magnitudes of black-holes [neutron-stars] are uniformly distributed between 0--0.98 [0 -- 0.4] and spin angles are isotropically distributed, the expected improvement in the average detection volume (at a fixed signal-to-noise-ratio threshold) of a search using this reduced-spin bank is $\sim20-52\%$, as compared to a search using a non-spinning bank.

Optimizing gravitational-wave searches for a population of coalescing binaries: Intrinsic parameters

We revisit the problem of searching for gravitational waves from inspiralling compact binaries in Gaussian coloured noise. For binaries with quasicircular orbits and non-precessing component spins, considering dominant mode emission only, if the intrinsic parameters of the binary are known then the optimal statistic for a single detector is the well-known two-phase matched filter. However, the matched filter signal-to-noise ratio is /not/ in general an optimal statistic for an astrophysical population of signals, since their distribution over the intrinsic parameters will almost certainly not mirror that of noise events, which is determined by the (Fisher) information metric. Instead, the optimal statistic for a given astrophysical distribution will be the Bayes factor, which we approximate using the output of a standard template matched filter search. We then quantify the possible improvement in number of signals detected for various populations of non-spinning binaries: for a distribution of signals uniformly distributed in volume and with component masses distributed uniformly over the range $1\leq m_{1,2}/M_\odot\leq 24$, $(m_1+m_2) /M_\odot\leq 25$ at fixed expected SNR, we find $\gtrsim 20\%$ more signals at a false alarm threshold of $10^{-6}\,$Hz in a single detector. The method may easily be generalized to binaries with non-precessing spins.

Gravitational Radiation from Post-Newtonian Sources and Inspiralling Compact Binaries.

To be observed and analyzed by the network of gravitational wave detectors on ground (LIGO, VIRGO, etc.) and by the future detectors in space (eLISA, etc.), inspiralling compact binaries — binary star systems composed of neutron stars and/or black holes in their late stage of evolution — require high-accuracy templates predicted by general relativity theory. The gravitational waves emitted by these very relativistic systems can be accurately modelled using a high-order post-Newtonian gravitational wave generation formalism. In this article, we present the current state of the art on post-Newtonian methods as applied to the dynamics and gravitational radiation of general matter sources (including the radiation reaction back onto the source) and inspiralling compact binaries. We describe the post-Newtonian equations of motion of compact binaries and the associated Lagrangian and Hamiltonian formalisms, paying attention to the self-field regularizations at work in the calculations. Several notions of innermost circular orbits are discussed. We estimate the accuracy of the post-Newtonian approximation and make a comparison with numerical computations of the gravitational self-force for compact binaries in the small mass ratio limit. The gravitational waveform and energy flux are obtained to high post-Newtonian order and the binary’s orbital phase evolution is deduced from an energy balance argument. Some landmark results are given in the case of eccentric compact binaries — moving on quasi-elliptical orbits with non-negligible eccentricity. The spins of the two black holes play an important role in the definition of the gravitational wave templates. We investigate their imprint on the equations of motion and gravitational wave phasing up to high post-Newtonian order (restricting to spin-orbit effects which are linear in spins), and analyze the post-Newtonian spin precession equations as well as the induced precession of the orbital plane.

Erratum: 2.5PN linear momentum flux from inspiralling compact binaries in quasicircular orbits and associated recoil: Nonspinning case [Phys. Rev. D85, 044021 (2012)

Erratum: 2.5PN linear momentum flux from inspiralling compact binaries in quasicircular orbits and associated recoil: Nonspinning case [Phys. Rev. D<b>85</b>, 044021 (2012)

The third and a half-post-Newtonian gravitational wave quadrupole mode for quasi-circular inspiralling compact binaries

We compute the quadrupole mode of the gravitational waveform of inspiralling compact binaries at the third and half-post-Newtonian (3.5PN) approximation of general relativity. The computation is performed using the multipolar post-Newtonian formalism, and restricted to binaries without spins moving on quasi-circular orbits. The new inputs mainly include the 3.5PN terms in the mass quadrupole moment of the source, and the control of required subdominant corrections to the contributions of hereditary integrals (tails and nonlinear memory effect). The result is given in the form of the quadrupolar mode (2, 2) in a spin-weighted spherical harmonic decomposition of the waveform, and may be used for comparison with the counterpart quantity computed in numerical relativity. It is a step towards the computation of the full 3.5PN waveform, whose knowledge is expected to reduce the errors on the location parameters of the source.

Generic bounds on dipolar gravitational radiation from inspiralling compact binaries

Various alternative theories of gravity predict dipolar gravitational radiation in addition to quadrupolar radiation. We show that gravitational wave (GW) observations of inspiralling compact binaries can put interesting constraints on the strengths of the dipole modes of GW polarizations. We put forward a physically motivated gravitational waveform for dipole modes, in the Fourier domain, in terms of two parameters: one which captures the relative amplitude of the dipole mode with respect to the quadrupole mode (α) and the other a dipole term in the phase (β). We then use this two-parameter representation to discuss typical bounds on their values using GW measurements. We obtain the expected bounds on the amplitude parameter α and the phase parameter β for advanced LIGO (AdvLIGO) and Einstein telescope (ET) noise power spectral densities using Fisher information matrix. AdvLIGO and ET may at best bound α to an accuracy of ∼10−2 and ∼10−3 and β to an accuracy of ∼10−5 and ∼10−6.

2.5PN linear momentum flux from inspiralling compact binaries in quasicircular orbits and associated recoil: Nonspinning case

Anisotropic emission of gravitational waves (GWs) from inspiralling compact binaries leads to the loss of linear momentum and hence gravitational recoil of the system. The loss rate of linear momentum in the far-zone of the source (a nonspinning binary system of black holes in quasicircular orbit) is investigated at the 2.5 post-Newtonian (PN) order and used to provide an analytical expression in harmonic coordinates for the 2.5PN accurate recoil velocity of the binary accumulated in the inspiral phase. We find that the recoil velocity at the end of the inspiral phase (i.e at the innermost stable circular orbit (ISCO)) is maximum for a binary with symmetric mass ratio of \nu~0.2 and is roughly about ~4.58 km/s. Going beyond inspiral, we also provide an estimate of the more important contribution to the recoil velocity from the plunge phase. Again the recoil velocity at the end of the plunge, involving contributions both from inspiral and plunge phase, is maximum for a binary with \nu~0.2 and is of the order of ~180 km/s.

Addressing the spin question in gravitational-wave searches: Waveform templates for inspiralling compact binaries with nonprecessing spins

This paper presents a post-Newtonian (PN) template family of gravitational waveforms from inspiralling compact binaries with nonprecessing spins, where the spin effects are described by a single ``reduced-spin'' parameter. This template family, which reparametrizes all the spin-dependent PN terms in terms of the leading-order (1.5PN) spin-orbit coupling term in an approximate way, has very high overlaps (fitting factor $&gt;0.99$) with nonprecessing binaries with arbitrary mass ratios and spins. We also show that this template family is ``effectual'' for the detection of a significant fraction of generic spinning binaries in the comparable-mass regime (${m}_{2}/{m}_{1}\ensuremath{\lesssim}10$), providing an attractive and feasible way of searching for gravitational waves from spinning low-mass binaries. We also show that the secular (nonoscillatory) spin-dependent effects in the phase evolution (which are taken into account by the nonprecessing templates) are more important than the oscillatory effects of precession in the comparable-mass (${m}_{1}\ensuremath{\simeq}{m}_{2}$) regime. Hence the effectualness of nonspinning templates is particularly poor in this case, as compared to non-precessing-spin templates. For the case of binary neutron stars observable by Advanced LIGO, even moderate spins (${\stackrel{^}{\mathbf{L}}}_{N}\ifmmode\cdot\else\textperiodcentered\fi{}\mathbf{S}/{m}^{2}\ensuremath{\simeq}0.015--0.1$) will cause considerable mismatches ($\ensuremath{\sim}3%--25%$) with nonspinning templates. This is contrary to the expectation that neutron-star spins may not be relevant for gravitational wave detection.

Spin induced multipole moments for the gravitational wave flux from binary inspirals to third Post-Newtonian order

Using effective field theory techniques we calculate the source multipole moments needed to obtain the spin contributions to the power radiated in gravitational waves from inspiralling compact binaries to third Post-Newtonian order (3PN). The multipoles depend linearly and quadratically on the spins and include both spin(1)spin(2) and spin(1)spin(1) components. The results in this paper provide the last missing ingredient required to determine the phase evolution to 3PN including all spin effects which we will report in a separate paper.