Search

For proposed third-generation gravitational-wave detectors such as the Einstein Telescope and Cosmic Explorer, whose sensitive bands are proposed to extend down to 5 Hz or below, the signals of low-mass compact binaries such as binary neutron stars remain in the detector's sensitive band long enough (up to a few days for the smallest proposed low-frequency cutoff of 1 Hz) that one cannot neglect the effects of the Earth's rotation on the detector's response and the changing Doppler shift of the signal. In the latter case, one also needs to consider the effects of the Earth's orbital motion, which is currently only included in analyses of compact binary signals using continuous wave techniques. These effects are also relevant for current detectors and signals from putative subsolar-mass binaries. Here we present simple Fourier series methods for computing these effects in the frequency domain, giving explicit expressions for the Earth's orbital motion in terms of low-order Fourier series, which will be sufficiently accurate for all compact binary signals except for those from very low-mass subsolar-mass binaries. The expression for the effects of the Earth's rotation on the antenna pattern functions does not use the stationary phase approximation (SPA), so we are able to show that the SPA is indeed quite accurate in these situations and present a Fourier series expression equivalent to it which is an order of magnitude faster. We also provide illustrations of these effects on detector sensitivity and the accumulation of information about various binary parameters with frequency., Comment: 31 pages, 11 figures

The observation of gravitational waves from compact binary coalescences is a promising tool to test the validity of general relativity (GR) in a highly dynamical strong-field regime. There are now a variety of tests of GR performed on the observed compact binary signals. In this paper, we propose a new test of GR that compares the results of these individual tests. This meta inspiral-merger-ringdown consistency test (IMRCT) involves inferring the final mass and spin of the remnant black hole obtained from the analyses of two different tests of GR and checking for consistency. If there is a deviation from GR, we expect that different tests of GR will recover different values for the final mass and spin, in general. We check the performance of the meta IMRCT using a standard set of null tests used in various gravitational-wave analyses: the original IMRCT, parameterized phasing tests (TIGER and FTI) and the modified dispersion test. However, the meta IMRCT is applicable to any tests of GR that infer the initial masses and spins or the final mass and spin, including ones that are applied to binary neutron star or neutron star--black hole signals. We apply the meta IMRCT to simulated quasi-circular GR and non-GR binary black hole (BBH) signals as well as to eccentric BBH signals in GR (analyzed with quasicircular waveforms). We find that the meta IMRCT gives consistency with GR for the quasi-circular GR signals and picks up a deviation from GR in the other cases, as do other tests. In some cases, the meta IMRCT finds a significant GR deviation for a given pair of tests (and specific testing parameters) while the individual tests do not, showing that it is more sensitive than the individual tests to certain types of deviations. In addition, we also apply this test to a few selected real compact binary signals and find them consistent with GR., Comment: 27 pages, 30 figures

General relativity (GR) has proven to be a highly successful theory of gravity since its inception. The theory has thrivingly passed numerous experimental tests, predominantly in weak gravity, low relative speeds, and linear regimes, but also in the strong-field and very low-speed regimes with binary pulsars. Observable gravitational waves (GWs) originate from regions of spacetime where gravity is extremely strong, making them a unique tool for testing GR, in previously inaccessible regions of large curvature, relativistic speeds, and strong gravity. Since their first detection, GWs have been extensively used to test GR, but no deviations have been found so far. Given GR's tremendous success in explaining current astronomical observations and laboratory experiments, accepting any deviation from it requires a very high level of statistical confidence and consistency of the deviation across GW sources. In this paper, we compile a comprehensive list of potential causes that can lead to a false identification of a GR violation in standard tests of GR on data from current and future ground-based GW detectors. These causes include detector noise, signal overlaps, gaps in the data, detector calibration, source model inaccuracy, missing physics in the source and in the underlying environment model, source misidentification, and mismodeling of the astrophysical population. We also provide a rough estimate of when each of these causes will become important for tests of GR for different detector sensitivities. We argue that each of these causes should be thoroughly investigated, quantified, and ruled out before claiming a GR violation in GW observations., Comment: Review article; 1 figure; 1 table; comments welcome

As the gravitational wave detector network is upgraded and the sensitivity of the detectors improves, novel scientific avenues open for exploration. For example, tests of general relativity will become more accurate as smaller deviations can be probed. Additionally, the detection of lensed gravitational waves becomes more likely. However, these new avenues could also interact with each other, and a gravitational wave event presenting deviations from general relativity could be mistaken for a lensed one. Here, we explore how phenomenological deviations from general relativity or binaries of exotic compact objects could impact those lensing searches focusing on a single event. We consider strong lensing, millilensing, and microlensing and find that certain phenomenological deviations from general relativity may be mistaken for all of these types of lensing. Therefore, our study shows that future candidate lensing events would need to be carefully examined to avoid a false claim of lensing where instead a deviation from general relativity has been seen.

The spin orientations of spinning binary black hole (BBH) mergers detected by ground-based gravitational wave detectors such as LIGO and Virgo can provide important clues about the formation of such binaries. However, these spin tilts, i.e., the angles between the spin vector of each black hole and the binary's orbital angular momentum vector, can change due to precessional effects as the black holes evolve from a large separation to their merger. The tilts inferred at a frequency in the sensitive band of the detectors by comparing the signal with theoretical waveforms can thus be significantly different from the tilts when the binary originally formed. These tilts at the binary's formation are well approximated in many scenarios by evolving the BBH backwards in time to a formally infinite separation. Using the tilts at infinite separation also places all binaries on an equal footing in analyzing their population properties. In this paper, we perform parameter estimation for simulated BBHs and investigate the differences between the tilts one infers directly close to merger and those obtained by evolving back to infinite separation. We select simulated observations such that their configurations show particularly large differences in their orientations close to merger and at infinity. While these differences may be buried in the statistical noise for current detections, we show that in future plus-era (A$+$ and Virgo$+$) detectors, they can be easily distinguished in some cases. We also consider the tilts at infinity for BBHs in various spin morphologies and at the endpoint of the up-down instability. In particular, we find that we are able to easily identify the up-down instability cases as such from the tilts at infinity., Comment: 18 pages, 10 figures

Detections of gravitational waves emitted from binary black hole coalescences allow us to probe the strong-field dynamics of general relativity (GR). One can compare the observed gravitational-wave signals with theoretical waveform models to constrain possible deviations from GR. Any physics that is not included in these waveform models might show up as apparent GR deviations. The waveform models used in current tests of GR describe binaries on quasicircular orbits, since most of the binaries detected by ground-based gravitational-wave detectors are expected to have negligible eccentricities. Thus, a signal from an eccentric binary in GR is likely to show up as a deviation from GR in the current implementation of these tests. We study the response of four standard tests of GR to eccentric binary black hole signals with the forecast O4 sensitivity of the LIGO-Virgo network. Specifically, we consider two parameterized tests (TIGER and FTI), the modified dispersion relation test, and the inspiral-merger-ringdown consistency test. To model eccentric signals, we use non-spinning numerical relativity simulations from the SXS catalog with three mass ratios $(1,2,3)$, which we scale to a redshifted total mass of $80M_\odot$ and luminosity distance of $400$ Mpc. For each of these mass ratios, we consider signals with eccentricities of $\sim0.05$ and $\sim 0.1$ at $17$ Hz. We find that signals with larger eccentricity lead to very significant false GR deviations in most tests while signals having smaller eccentricity lead to significant deviations in some tests. For the larger eccentricity cases, one would even get a deviation from GR with TIGER at $\sim 90\%$ credibility at a distance of $\gtrsim 1.5$ Gpc. Thus, it will be necessary to exclude the possibility of an eccentric binary in order to make any claim about detecting a deviation from GR., Comment: 16 pages, 6 figures, version accepted by PRD

The precessional motion of binary black holes can be classified into one of three morphologies, based on the evolution of the angle between the components of the spins in the orbital plane: Circulating, librating around 0, and librating around $\pi$. These different morphologies can be related to the binary's formation channel and are imprinted in the binary's gravitational wave signal. In this paper, we develop a Bayesian model selection method to determine the preferred spin morphology of a detected binary black hole. The method involves a fast calculation of the morphology which allows us to restrict to a specific morphology in the Bayesian stochastic sampling. We investigate the prospects for distinguishing between the different morphologies using gravitational waves in the Advanced LIGO/Advanced Virgo network with their plus-era sensitivities. For this, we consider fiducial high- and low-mass binaries having different spin magnitudes and signal-to-noise ratios (SNRs). We find that in the cases with high spin and high SNR, the true morphology is strongly favored with $\log_{10}$ Bayes factors $\gtrsim 4$ compared to both alternative morphologies when the binary's parameters are not close to the boundary between morphologies. However, when the binary parameters are close to the boundary between morphologies, only one alternative morphology is strongly disfavored. In the low-spin, high-SNR cases, the true morphology is still favored with a $\log_{10}$ Bayes factor $\sim 2$ compared to one alternative morphology. We also consider the gravitational wave signal from GW200129_065458 that has some evidence for precession (modulo data quality issues) and find that there is no preference for a specific morphology. Our method for restricting the prior to a given morphology is publicly available through an easy-to-use Python package called bbh_spin_morphology_prior. (Abridged), Comment: 14 pages, 5 figures, version accepted by PRD, including the correction from the erratum

We present a method to construct conformally curved initial data for charged black hole binaries with spin on arbitrary orbits. We generalize the superposed Kerr-Schild, extended conformal thin sandwich construction from [Lovelace et al., Phys. Rev. D {78}, 084017 (2008)] to use Kerr-Newman metrics for the superposed black holes and to solve the electromagnetic constraint equations. We implement the construction in the pseudospectral code SGRID. The code thus provides a complementary and completely independent excision-based construction, compared to the existing charged black hole initial data constructed using the puncture method [Bozzola and Paschalidis, Phys. Rev. D {99}, 104044 (2019)]. It also provides an independent implementation (with some small changes) of the Lovelace et al. vacuum construction. We construct initial data for different configurations of orbiting binaries, e.g., with black holes that are highly charged or rapidly spinning (90 and 80 percent of the extremal values, respectively, for this initial test, though the code should be able to produce data with even higher values of these parameters using higher resolutions), as well as for generic spinning, charged black holes. We carry out exploratory evolutions with the finite difference, moving punctures codes BAM (in the vacuum case) and HAD (for head-on collisions including charge), filling inside the excision surfaces. In the charged case, evolutions of these initial data provide a proxy for binary black hole waveforms in modified theories of gravity. Moreover, the generalization of the construction to Einstein-Maxwell-dilaton theory should be straightforward., Comment: 17 pages, 7 figures

We construct a template to model the post-merger phase of a binary black hole coalescence in the presence of a remnant $U(1)$ charge. We include the quasi-normal modes typically dominant during a binary black hole coalescence, $(\ell,m,n) = \{(2,2,0), (2,2,1)\}$ and also present analytical fits for the quasinormal mode frequencies of a Kerr-Newman black hole in terms of its spin and charge, here also including the $(3,3,0)$ mode. Aside from astrophysical electric charge, our template can accommodate extensions of the Standard Model, such as a dark photon. Applying the model to LIGO-Virgo detections, we find that we are unable to distinguish between the charged and uncharged hypotheses from a purely post-merger analysis of the current events. However, restricting the mass and spin to values compatible with the analysis of the full signal, we obtain a 90th percentile bound $\bar{q} < 0.33$ on the black hole charge-to-mass ratio, for the most favorable case of GW150914. Under similar assumptions, by simulating a typical loud signal observed by the LIGO-Virgo network at its design sensitivity, we assess that this model can provide a robust measurement of the charge-to-mass ratio only for values $\bar{q} \gtrsim 0.5$; here we also assume that the mode amplitudes are similar to the uncharged case in creating our simulated signal. Lower values, down to $\bar{q} \sim 0.3$, could instead be detected when evaluating the consistency of the pre-merger and post-merger emission., Comment: 21 pages, 11 figures, 4 tables. Matches published version

Gravitational wave observations of compact binary coalescences provide precision probes of strong-field gravity. There is thus now a standard set of null tests of general relativity (GR) applied to LIGO-Virgo detections and many more such tests proposed. However, the relation between all these tests is not yet well understood. We start to investigate this by applying a set of standard tests to simulated observations of binary black holes in GR and with phenomenological deviations from GR. We consider four types of tests: residuals, inspiral-merger-ringdown consistency, parameterized phasing (two varieties), and modified dispersion relation. We also check the consistency of the unmodeled reconstruction of the waveforms with the waveform recovered using GR templates. These tests are applied to simulated observations similar to GW150914 with both large and small deviations from GR and similar to GW170608 just with small deviations from GR. We find that while very large deviations from GR are picked up with high significance by almost all tests, more moderate deviations are picked up by only a few tests, and some deviations are not recognized as GR violations by any test at the moderate signal-to-noise ratios we consider. Moreover, the tests that identify various deviations with high significance are not necessarily the expected ones. In particular, the parameterized tests recover PN test parameters much closer to zero than their true values in some cases. Additionally, we find that of the GR deviations we consider, the residuals test is only able to detect extreme deviations from GR that no longer look like binary black hole coalescences in GR. The reconstruction comparison shows more promise for detecting relatively small GR deviations in an unmodeled framework, at least for high-mass systems., Comment: 29 pages, 16 figures. Accepted manuscript

Two important parameters inferred from the gravitational wave signals of binaries of precessing black holes are the spin tilt angles, i.e., the angles at which the black holes' spin axes are inclined with respect to the binary's orbital angular momentum. The LIGO-Virgo parameter estimation analyses currently provide spin tilts at a fiducial reference frequency, often the lowest frequency used in the data analysis. However, the most astrophysically interesting quantities are the spin tilts when the binary was formed, which can be significantly different from those at the reference frequency for strongly precessing binaries. The spin tilts at formally infinite separation are a good approximation to the tilts at formation in many formation channels and can be computed efficiently for binary black holes using precession-averaged evolution. Here, we present a new code for computing the tilts at infinity that combines the precession-averaged evolution with orbit-averaged evolution at high frequencies and illustrate its application to GW190521 and other binary black hole detections from O3a. We have empirically determined the transition frequency between the orbit-averaged and precession-averaged evolution to produce tilts at infinity with a given accuracy. We also have regularized the precession-averaged equations in order to obtain good accuracy for the very close-to-equal-mass binary parameters encountered in practice. This additionally allows us to investigate the singular equal-mass limit of the precession-averaged expressions, where we find an approximate scaling of $1/(1 - q)$ with the mass ratio $q$., Comment: 25 pages, 16 figures

(Abridged) While the gravitational-wave (GW) signal GW170817 was accompanied by a variety of electromagnetic (EM) counterparts, sufficiently high-mass binary neutron star (BNS) mergers are expected to be unable to power bright EM counterparts. The putative high-mass binary BNS merger GW190425, for which no confirmed EM counterpart has been identified, may be an example of such a system. It is thus important to understand how well we will be able to distinguish high-mass BNSs and low-mass binary black holes (BBHs) solely from their GW signals. To do this, we consider the imprint of the tidal deformability of the neutron stars on the GW signal for systems undergoing prompt black hole formation after merger. We model the BNS signals using hybrid numerical relativity -- tidal effective-one-body waveforms. Specifically, we consider a set of five nonspinning equal-mass BNS signals with masses of 2.7, 3.0, 3.2 Msun and with three different equations of state, as well as the analogous BBH signals. We perform parameter estimation on these signals in three networks: Advanced LIGO-Advanced Virgo and Advanced LIGO-Advanced Virgo-KAGRA with sensitivities similar to O3 and O4, respectively, and a 3G network of two Cosmic Explorers (CEs) and one Einstein Telescope, with a CE sensitivity similar to Stage 2. Our analysis suggests that we cannot distinguish the signals from high-mass BNSs and BBHs at a 90% credible level with the O3-like network even at 40 Mpc. However, we can distinguish all but the most compact BNSs that we consider in our study from BBHs at 40 Mpc at a >= 95% credible level using the O4-like network, and can even distinguish them at a > 99.2% (>= 97%) credible level at 369 (835) Mpc using the 3G network. Additionally, we present a simple method to compute the leading effect of the Earth's rotation on the response of a gravitational wave detector in the frequency domain., Comment: 18 pages, 6 figures. Accepted manuscript

The "Shapiro delay" experienced by an astronomical messenger traveling through a gravitational field has been used to place constraints on possible deviations from the equivalence principle. The standard Shapiro delay used to obtain these constraints is not itself an observable in general relativity, but is rather obtained by comparing with a fiducial Euclidean distance. There is not a mapping between the constraints obtained in this manner and alternative theories that exhibit equivalence principle violations. However, even assuming that the comparison with the fiducial Euclidean distance is carried out in a way that is useful for some class of alternative theories, we show that the standard calculation of these constraints cannot be applied on cosmological scales, as is often done. Specifically, we find that the Shapiro delay computed in the standard way (taking the Newtonian potential to vanish at infinity) diverges as one includes many remote sources. We use an infinite homogeneous lattice model to illustrate this divergence, and also show how the divergence can be cured by using Fermi coordinates associated with an observer. With this correction, one finds that the Shapiro delay is no longer monotonic with the number of sources. Thus, one cannot compute a conservative lower bound on the Shapiro delay using a subset of the sources of the gravitational field without further assumptions and/or observational input. As an illustration, we compute the Shapiro delay by applying the Fermi coordinate expression to two catalogs of galaxy clusters, illustrating the dependence of the result on the completeness of the catalogue and the mass estimates., Comment: 15 pages, 7 figures. Accepted manuscript in Phys. Rev. D

Accurate and fast gravitational waveform (GW) models are essential to extract information about the properties of compact binary systems that generate GWs. Building on previous work, we present an extension of the NRTidal model for binary neutron star (BNS) waveforms. The upgrades are: (i) a new closed-form expression for the tidal contribution to the GW phase which includes further analytical knowledge and is calibrated to more accurate numerical relativity data than previously available; (ii) a tidal correction to the GW amplitude; (iii) an extension of the spin-sector incorporating equation-of-state-dependent finite size effects at quadrupolar and octupolar order; these appear in the spin-spin tail terms and cubic-in-spin terms, both at 3.5PN. We add the new description to the precessing binary black hole waveform model IMRPhenomPv2 to obtain a frequency-domain precessing binary neutron star model. In addition, we extend the SEOBNRv4_ROM and IMRPhenomD aligned-spin binary black hole waveform models with the improved tidal phase corrections. Focusing on the new IMRPhenomPv2_NRTidalv2 approximant, we test the model by comparing with numerical relativity waveforms as well as hybrid waveforms combining tidal effective-one-body and numerical relativity data. We also check consistency against a tidal effective-one-body model across large regions of the BNS parameter space., Comment: Accepted manuscript

We summarize the four tests of general relativity carried out using the ten binary black hole signals detected by Advanced LIGO and Advanced Virgo and included in their first catalog, GWTC-1. These events provide unprecedented opportunities for testing the predictions of general relativity for the gravitational waveforms from these highly dynamical, strong-field events. The first two tests check the consistency of the residuals with noise and the consistency of the low- and high-frequency parts of the signal. The other two tests check that parameterized deviations in the waveform model---including in the post-Newtonian coefficients---are consistent with zero, and that the propagation of the waves is nondispersive. These tests reveal no evidence for deviations from general relativity, and the combined constraints improve previous results by factors of up to 2.5. We also check that the binary black hole signals observed by all three detectors do not give stronger constraints on alternative polarizations than those obtained from GW170817., Comment: 8 pages, 3 figures, contribution to the 2019 Gravitation session of the 54th Rencontres de Moriond; summary of arXiv:1903.04467; v2: small updates to some results (includes some updates that did not make it into the published version)

As gravitational-wave (GW) observations of binary black holes are becoming a precision tool for physics and astronomy, several subdominant effects in the GW signals need to be accurately modeled. Previous studies have shown that neglecting subdominant modes in the GW templates causes an unacceptable loss in detection efficiency and large systematic errors in the estimated parameters for binaries with large mass ratios. Our recent work [Mehta et al., Phys. Rev. D 96, 124010 (2017)] constructed a phenomenological gravitational waveform family for nonspinning black-hole binaries that includes subdominant spherical harmonic modes $(\ell = 2, m = \pm 1)$, $(\ell = 3, m = \pm 3)$, and $(\ell = 4, m = \pm 4)$ in addition to the dominant quadrupole mode, $(\ell = 2, m=\pm 2)$. In this article, we construct analytical models for the ($\ell = 3, m = \pm 2$) and ($\ell = 4, m = \pm 3$) modes and include them in the existing waveform family. Accurate modeling of these modes is complicated by the mixing of multiple spheroidal harmonic modes. We develop a method for accurately modeling the effect of mode mixing, thus producing an analytical waveform family that has faithfulness >99.6%., Comment: 12 pages, 11 figures, Accepted manuscript

As current gravitational wave (GW) detectors increase in sensitivity, and particularly as new instruments are being planned, there is the possibility that ground-based GW detectors will observe GWs from highly eccentric neutron star binaries. We present the first detailed study of highly eccentric BNS systems with full (3+1)D numerical relativity simulations using consistent initial conditions, i.e., setups which are in agreement with the Einstein equations and with the equations of general relativistic hydrodynamics in equilibrium. Overall, our simulations cover two different equations of state (EOSs), two different spin configurations, and three to four different initial eccentricities for each pairing of EOS and spin. We extract from the simulated waveforms the frequency of the f-mode oscillations induced during close encounters before the merger of the two stars. The extracted frequency is in good agreement with f-mode oscillations of individual stars for the irrotational cases, which allows an independent measure of the supranuclear equation of state not accessible for binaries on quasi-circular orbits. The energy stored in these f-mode oscillations can be as large as $10^{-3}M_\odot \sim 10^{51}$ erg, even with a soft EOS. In order to estimate the stored energy, we also examine the effects of mode mixing due to the stars' offset from the origin on the f-mode contribution to the GW signal. While in general (eccentric) neutron star mergers produce bright electromagnetic counterparts, we find that the luminosity decreases when the eccentricity becomes too large, due to a decrease of the ejecta mass. Finally, the use of consistent initial configurations also allows us to produce high-quality waveforms for different eccentricities which can be used as a testbed for waveform model development of highly eccentric binary neutron star systems., Comment: (26 pages, 17 figures), "CoRe webpage, see http://www.computational-relativity.org/"

LIGO and Virgo have recently observed a number of gravitational wave (GW) signals that are fully consistent with being emitted by binary black holes described by general relativity. However, there are theoretical proposals of exotic objects that can be massive and compact enough to be easily confused with black holes. Nevertheless, these objects differ from black holes in having nonzero tidal deformabilities, which can allow one to distinguish binaries containing such objects from binary black holes using GW observations. Using full Bayesian parameter estimation, we investigate the possibility of constraining the parameter space of such "black hole mimickers" with upcoming GW observations. Employing perfect fluid stars with a polytropic equation of state as a simple model that can encompass a variety of possible black hole mimickers, we show how the observed masses and tidal deformabilities of a binary constrain the equation of state. We also show how such constraints can be used to rule out some simple models of boson stars., Comment: 13 pages, 5 figures; v3: Accepted manuscript

Advanced LIGO's recent observations of gravitational waves (GWs) from merging binary black holes have opened up a unique laboratory to test general relativity (GR) in the highly relativistic regime. One of the tests used to establish the consistency of the first LIGO event with a binary black hole merger predicted by GR was the inspiral-merger-ringdown consistency test. This involves inferring the mass and spin of the remnant black hole from the inspiral (low-frequency) part of the observed signal and checking for the consistency of the inferred parameters with the same estimated from the post-inspiral (high-frequency) part of the signal. Based on the observed rate of binary black hole mergers, we expect the advanced GW observatories to observe hundreds of binary black hole mergers every year when operating at their design sensitivities, most of them with modest signal to noise ratios (SNRs). Anticipating such observations, this paper shows how constraints from a large number of events with modest SNRs can be combined to produce strong constraints on deviations from GR. Using kludge modified GR waveforms, we demonstrate how this test could identify certain types of deviations from GR if such deviations are present in the signal waveforms. We also study the robustness of this test against reasonable variations of a variety of different analysis parameters., Comment: 23 pages, 10 figures

The coalescences of stellar-mass black-hole binaries through their inspiral, merger, and ringdown are among the most promising sources for ground-based gravitational-wave (GW) detectors. If a GW signal is observed with sufficient signal-to-noise ratio, the masses and spins of the black holes can be estimated from just the inspiral part of the signal. Using these estimates of the initial parameters of the binary, the mass and spin of the final black hole can be uniquely predicted making use of general-relativistic numerical simulations. In addition, the mass and spin of the final black hole can be independently estimated from the merger--ringdown part of the signal. If the binary black hole dynamics is correctly described by general relativity (GR), these independent estimates have to be consistent with each other. We present a Bayesian implementation of such a test of general relativity, which allows us to combine the constraints from multiple observations. Using kludge modified GR waveforms, we demonstrate that this test can detect sufficiently large deviations from GR, and outline the expected constraints from upcoming GW observations using the second-generation of ground-based GW detectors., Comment: 5 pages, 2 figs

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., Comment: 11 pages, 5 figures

Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spins, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $\sim10^{-2}M_\odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged], Comment: (37pages, 26 figures)

It is now possible to compute linear in mass-ratio terms in the post-Newtonian (PN) expansion for compact binaries to very high orders using black hole perturbation theory applied to various invariants. For instance, a computation of the redshift invariant of a point particle in a circular orbit about a black hole in linear perturbation theory gives the linear-in-mass-ratio portion of the binding energy of a circular binary with arbitrary mass ratio. This binding energy, in turn, encodes the system's conservative dynamics. We give a method for extracting the analytic forms of these PN coefficients from high-accuracy numerical data using experimental mathematics techniques, notably an integer relation algorithm. Such methods should be particularly important when the calculations progress to the considerably more difficult case of perturbations of the Kerr metric. As an example, we apply this method to the redshift invariant in Schwarzschild. Here we obtain analytic coefficients to 12.5PN, and higher-order terms in mixed analytic-numerical form to 21.5PN, including analytic forms for the complete 13.5PN coefficient, and all the logarithmic terms at 13PN. At these high orders, an individual coefficient can have over 30 terms, including a wide variety of transcendental numbers, when written out in full. We are still able to obtain analytic forms for such coefficients from the numerical data through a careful study of the structure of the expansion. The structure we find also allows us to predict certain "leading logarithm"-type contributions to all orders. The additional terms in the expansion we obtain improve the accuracy of the PN series for the redshift observable, even in the very strong-field regime inside the innermost stable circular orbit, particularly when combined with exponential resummation., Comment: 25 pages, 3 figures; version accepted by PRD, including various small additions and corrections

Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipse's inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero., Comment: 23 pages, 10 figures

(Abridged) High-order terms in the post-Newtonian (PN) expansions of various quantities for compact binaries exhibit a combinatorial increase in complexity, including ever-increasing numbers of transcendentals. Here we consider the gravitational wave energy flux at infinity from a point particle in a circular orbit around a Schwarzschild black hole, which is known to 22PN beyond the lowest-order Newtonian prediction, at which point each order has over 1000 terms. We introduce a factorization that considerably simplifies the spherical harmonic modes of the energy flux (and thus also the amplitudes of the spherical harmonic modes of the gravitational waves); it is likely that much of the complexity this factorization removes is due to curved-space wave propagation (e.g., tail effects). For the modes with azimuthal number l of 7 or greater, this factorization reduces the expressions for the modes that enter the 22PN total energy flux to pure integer PN series with rational coefficients, which amounts to a reduction of up to a factor of ~150 in the total number of terms in a given mode. The reduction in complexity becomes less dramatic for smaller l, due to the structure of the expansion, though the factorization is still able to remove all the half-integer PN terms. For the 22PN l = 2 modes, this factorization still reduces the total number of terms (and size) by a factor of ~10 and gives purely rational coefficients through 8PN. This factorization also improves the convergence of the series, though we find the exponential resummation introduced for the full energy flux by Isoyama et al. to be even more effective at improving the convergence of the individual modes, producing improvements of over four orders of magnitude over the original series for some modes. However, the exponential resummation is not as effective at simplifying the series, particularly for the higher-order modes., Comment: 18 pages, 4 figures, version accepted by PRD, including various minor corrections

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters., Comment: 51 pages, 10 figures; published version

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients., Comment: 33 pages, 2 figures

We show that there is a direct relation between upper limits on (or potential future measurements of) the m = 2 quadrupole moments of slowly rotating neutron stars and the l = m = 2 deformation of the star's surface, in full general relativity, to first order in the perturbation. This relation only depends on the star's structure through its mass and radius. All one has to assume about the star's constituents is that the stress-energy tensor at its surface is that of a perfect fluid, which will be true with good accuracy in almost all the situations of interest, and that the magnetic field configuration there is force-free, which is likely to be a good approximation. We then apply this relation to the stars which have direct LIGO/Virgo bounds on their m = 2 quadrupole moment, below the spin-down limit, and compare with the expected surface deformations due to rotation. In particular, we find that LIGO observations have constrained the Crab pulsar's l = m = 2 surface deformation to be smaller than its l = 2, m = 0 deformation due to rotation, for all the causal equations of state we consider, a statement that would not have been able to be made just using the upper bounds on the l = m = 2 deformation from electromagnetic observations., Comment: 13 pages, 3 figures, version accepted by PRD, including expanded discussion of the calculation and a few corrections

We present a method for calculating the maximum elastic quadrupolar deformations of relativistic stars, generalizing the previous Newtonian, Cowling approximation integral given by [G. Ushomirsky et al., Mon. Not. R. Astron. Soc. 319, 902 (2000)]. (We also present a method for Newtonian gravity with no Cowling approximation.) We apply these methods to the m = 2 quadrupoles most relevant for gravitational radiation in three cases: crustal deformations, deformations of crystalline cores of hadron-quark hybrid stars, and deformations of entirely crystalline color superconducting quark stars. In all cases, we find suppressions of the quadrupole due to relativity compared to the Newtonian Cowling approximation, particularly for compact stars. For the crust these suppressions are up to a factor ~6, for hybrid stars they are up to ~4, and for solid quark stars they are at most ~2, with slight enhancements instead for low mass stars. We also explore ranges of masses and equations of state more than in previous work, and find that for some parameters the maximum quadrupoles can still be very large. Even with the relativistic suppressions, we find that 1.4 solar mass stars can sustain crustal quadrupoles of a few times 10^39 g cm^2 for the SLy equation of state or close to 10^40 g cm^2 for equations of state that produce less compact stars. Solid quark stars of 1.4 solar masses can sustain quadrupoles of around 10^44 g cm^2. Hybrid stars typically do not have solid cores at 1.4 solar masses, but the most massive ones (~2 solar masses) can sustain quadrupoles of a few times 10^41 g cm^2 for typical microphysical parameters and a few times 10^42 g cm^2 for extreme ones. All of these quadrupoles assume a breaking strain of 0.1 and can be divided by 10^45 g cm^2 to yield the fiducial "ellipticities" quoted elsewhere., Comment: 21 pages, 11 figures, version accepted by PRD, including the corrected maximum hybrid star quadrupoles (from the erratum to the shear modulus calculation) and the corrected binding energy computation

Robust arguments predict that a hadron-quark mixed phase may exist in the cores of some "neutron" stars. Such a phase forms a crystalline lattice with a shear modulus higher than that of the crust due to the high density and charge separation, even allowing for the effects of charge screening. This may lead to strong continuous gravitational-wave emission from rapidly rotating neutron stars and gravitational-wave bursts associated with magnetar flares and pulsar glitches. We present the first detailed calculation of the shear modulus of the mixed phase. We describe the quark phase using the bag model plus first-order quantum chromodynamics corrections and the hadronic phase using relativistic mean-field models with parameters allowed by the most massive pulsar. Most of the calculation involves treating the "pasta phases" of the lattice via dimensional continuation, and we give a general method for computing dimensionally continued lattice sums including the Debye model of charge screening. We compute all the shear components of the elastic modulus tensor and angle average them to obtain the effective (scalar) shear modulus for the case where the mixed phase is a polycrystal. We include the contributions from changing the cell size, which are necessary for the stability of the lower-dimensional portions of the lattice. Stability also requires a minimum surface tension, generally tens of MeV/fm^2 depending on the equation of state. We find that the shear modulus can be a few times 10^33 erg/cm^3, two orders of magnitude higher than the first estimate, over a significant fraction of the maximum mass stable star for certain parameter choices., Comment: 22 pages, 12 figures, version accepted by Phys. Rev. D, with the corrections to the shear modulus computation and Table I given in the erratum

We generalize the standard Poisson summation formula for lattices so that it operates on the level of theta series, allowing us to introduce noninteger dimension parameters (using the dimensionally continued Fourier transform). When combined with one of the proofs of the Jacobi imaginary transformation of theta functions that does not use the Poisson summation formula, our proof of this generalized Poisson summation formula also provides a new proof of the standard Poisson summation formula for dimensions greater than 2 (with appropriate hypotheses on the function being summed). In general, our methods work to establish the (Voronoi) summation formulae associated with functions satisfying (modular) transformations of the Jacobi imaginary type by means of a density argument (as opposed to the usual Mellin transform approach). In particular, we construct a family of generalized theta series from Jacobi theta functions from which these summation formulae can be obtained. This family contains several families of modular forms, but is significantly more general than any of them. Our result also relaxes several of the hypotheses in the standard statements of these summation formulae. The density result we prove for Gaussians in the Schwartz space may be of independent interest., Comment: 12 pages, version accepted by JFAA, with various additions and improvements

(Abridged) By asymptotically matching a post-Newtonian (PN) metric to two tidally perturbed Schwarzschild metrics, we generate approximate initial data (in the form of a 4-metric) for a nonspinning black hole binary in a circular orbit. We carry out this matching through O(v^4) in the binary's orbital velocity v, so the resulting data are conformally curved. Far from the holes, we use the appropriate PN metric that accounts for retardation, which we construct using the highest-order PN expressions available to compute the binary's past history. The data set's uncontrolled remainders are thus O(v^5) throughout the timeslice; we also generate an extension to the data set that has uncontrolled remainders of O(v^6) in the purely PN portion of the timeslice (i.e., not too close to the holes). The resulting data are smooth, since we join all the metrics together by smoothly interpolating between them. We perform this interpolation using transition functions constructed to avoid introducing excessive additional constraint violations. Due to their inclusion of tidal deformations and outgoing radiation, these data should substantially reduce the initial spurious ("junk") radiation observed in current simulations that use conformally flat initial data. Such reductions in the nonphysical components of the initial data will be necessary for simulations to achieve the accuracy required to supply Advanced LIGO and LISA with the templates necessary for parameter estimation., Comment: 48 pages, 23 figures