Gravitational waves detected from well-localized inspiraling binaries would allow to determine, directly and independently, both binary luminosity and redshift. In this case, such systems could behave as standard candles providing an excellent probe of cosmic distances up to z < 0.1 and thus complementing other indicators of cosmological distance ladder.
It has been recently shown that quadruply lensed gravitational-wave (GW) events due to coalescing binaries can be localized to one or just a few galaxies, even in the absence of an electromagnetic counterpart. We discuss how this can be used to extract information on modified GW propagation, which is a crucial signature of modifications of gravity at cosmological scales. We show that, using quadruply lensed systems, it is possible to constrain the parameter $Xi_0$ that characterizes modified GW propagation, without the need of imposing a prior on $H_0$. A LIGO/Virgo/Kagra network at target sensitivity might already get a significant measurement of $Xi_0$, while a third generation GW detector such as the Einstein Telescope could reach a very interesting accuracy.
Waveforms of gravitational waves provide information about a variety of parameters for the binary system merging. However, standard calculations have been performed assuming a FLRW universe with no perturbations. In reality this assumption should be dropped: we show that the inclusion of cosmological perturbations translates into corrections to the estimate of astrophysical parameters derived for the merging binary systems. We compute corrections to the estimate of the luminosity distance due to velocity, volume, lensing and gravitational potential effects. Our results show that the amplitude of the corrections will be negligible for current instruments, mildly important for experiments like the planned DECIGO, and very important for future ones such as the Big Bang Observer.
Extending prior work by Pankow et al, we introduce RIFT, an algorithm to perform Rapid parameter Inference on gravitational wave sources via Iterative Fitting. We demonstrate this approach can correctly recover the parameters of coalescing compact binary systems, using detailed comparisons of RIFT to the well-tested LALInference software library. We provide several examples where the unique speed and flexibility of RIFT enables otherwise intractable or awkward parameter inference analyses, including (a) adopting either costly and novel models for outgoing gravitational waves; and (b) mixed approximations, each suitable to different parts of the compact binary parameter space. We demonstrate how RIFT{} can be applied to binary neutron stars, both for parameter inference and direct constraints on the nuclear equation of state.
We introduce a method based on the loudest event statistic to calculate an upper limit or interval on the astrophysical rate of binary coalescence. The calculation depends upon the sensitivity and noise background of the detectors, and a model for the astrophysical distribution of coalescing binaries. There are significant uncertainties in the calculation of the rate due to both astrophysical and instrumental uncertainties as well as errors introduced by using the post--Newtonian waveform to approximate the full signal. We catalog these uncertainties in detail and describe a method for marginalizing over them. Throughout, we provide an example based on the initial LIGO detectors.
We develop the foundations of an effective-one-body (EOB) model for eccentric binary coalescences that includes the conservative dynamics, radiation reaction, and gravitational waveform modes from the inspiral and the merger-ringdown signals. We use the same approach as is commonly employed in black-hole perturbation theory by introducing a relativistic parameterization of the dynamics that is defined by the orbital geometry and consists of a set of phase variables and quantities that evolve only due to gravitational radiation reaction. Specializing to nonspinning binaries, we derive the EOB evolution equations and compute the binarys radiative multipole moments that determine the gravitational waves through a decomposition into the fundamental frequencies of the motion. The major differences between our treatment and the quasi-Keplerian approach often used in post-Newtonian (PN) calculations are that the orbital parameters describe strong-field dynamics, and that expressing the multipole moments in terms of the frequencies simplifies the calculations and also results in an unambiguous orbit-averaging operation. While our description of the conservative dynamics is fully relativistic, we limit explicit derivations in the radiative sector to 1.5PN order for simplicity. This already enables us to establish methods for computing both instantaneous and hereditary contributions to the gravitational radiation in EOB coordinates that have straightforward extensions to higher PN order. The weak-field, small eccentricity limit of our results for the orbit-averaged fluxes of energy and angular momentum agrees with known PN results when expressed in terms of gauge-invariant quantities. We further address considerations for the numerical implementation of the model and the completion of the waveforms to include the merger and ringdown signals, and provide illustrative results.