The thermodynamic relation between pressure and density (i.e. the equation of state) of cold supranuclear matter is critical in describing neutron stars, yet it remains one of the largest uncertainties in nuclear physics. The extraction of tidal deformabilities from the gravitational waves emitted in the coalescence of neutron star binaries, such as GW170817, is a promising tool to probe this thermodynamic relation. Equation-of-state insensitive relations between symmetric and antisymmetric combinations of individual tidal deformabilities, the so-called "binary Love relations", have proven important to infer the radius of neutron stars, and thus constrain the equation of state, from such gravitational waves. A similar set of relations between the moment of inertia, the tidal deformability, the quadrupole moment, and the compactness of neutron stars, the so-called `I-Love-Q' and `C-Love' relations, allow for future tests of General Relativity in the extreme gravity regime. But even the most insensitive of such relations still presents some degree of equation-of-state variability that could introduce systematic uncertainties in parameter extraction and in model selection. We here reduce this variability by more than 50% by imposing a prior on the allowed set of equations of state, derived from the posteriors generated from the analysis of GW170817. The resulting increase in insensitivity reduces systematic uncertainties in the extraction of the tidal deformability from future gravitational wave observations, although statistical uncertainties currently dominate the error budget, and will continue to do so until the era of Voyager-class detectors.

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing non-trivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

The parameterized post-Einsteinian formalism was developed to search for generic deviations from general relativity with gravitational waves. We here present constraints on this framework using Bayesian analysis of a set of binary pulsar observations. In particular, we use measurements of the Keplerian and post-Keplerian parameters of six different binary pulsar systems, and Markov-Chain Monte-Carlo exploration to calculate posteriors on the parameterized post-Einsteinian parameters and derive robust constraints. We find improvements of 1-2 orders of magnitude in the strength of constraints when combining all six observations, relative to what one can achieve when using only the double binary pulsar. We also find that the constraints are robust to covariances with the binary's component masses. The bounds on the parameterized post-Einsteinian framework derived here could be used as a prior in future Bayesian tests of general relativity with gravitational wave observations.

The Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface eccentricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by seven orders of magnitude.

Gravitational-wave astronomy offers not only new vistas into the realm of astrophysics, but it also opens an avenue for probing, for the first time, general relativity in its strong-field, nonlinear, and dynamical regime, where the theory's predictions manifest themselves in their full glory. We present a study of whether the gravitational wave events detected so far by the LIGO-Virgo scientific collaborations can be used to probe higher-curvature corrections to general relativity. In particular, we focus on two examples: Einstein-dilaton-Gauss-Bonnet and dynamical Chern-Simons gravity. We find that the two events with a low-mass $m \approx 7 M_{\odot}$ BH (GW151226 and GW170608) place stringent constraints on Einstein-dilaton-Gauss-Bonnet gravity, ${\alpha}^{1/2}_{\rm EdGB} \lesssim 5.6$ km, whereas dynamical Chern-Simons gravity remains unconstrained by the gravitational-wave observations analyzed.

Dynamical Chern-Simons gravity cannot be strongly constrained with current experiments because it reduces to General Relativity in the weak-field limit. This theory, however, introduces modifications in the non-linear, dynamical regime, and thus, it could be greatly constrained with gravitational waves from the late inspiral of black hole binaries. We complete the first self-consistent calculation of such gravitational waves in this theory. For favorable spin-orientations, advanced ground-based detectors may improve existing solar-system constraints by 6 orders of magnitude.

We connect two formalisms that describe deformations away from general relativity, one valid in the strong-field regime of neutrons stars and another valid in the radiative regime of gravitational waves: the post-Tolman-Oppenheimer-Volkoff and the parametrized-post-Einsteinian formalisms respectively. We find that post-Tolman-Oppenheimer-Volkoff deformations of the exterior metric of an isolated neutron star induce deformations in the orbital binding energy of a neutron star binary. Such a modification to the binding energy then percolates into the gravitational waves emitted by such a binary, with the leading-order post-Tolman-Oppenheimer-Volkoff modifications introducing a second post-Newtonian order correction to the gravitational wave phase. The lack of support in gravitational wave data for general relativity deformations at this post-Newtonian order can then be used to place constraints post-Tolman-Oppenheimer-Volkoff parameters. As an application, we use the binary neutron star merger event GW170817 to place the constraint $-2.4 \leq \chi \leq 44$ (at 90% credibility) on a combination of post-Tolman-Oppenheimer-Volkoff parameters. We also explore the implications of this result to the possible deformations of the mass-radius relation of neutron stars allowed within this formalism. This work opens the path towards theory-independent tests of gravity, combining astronomical observations of neutron stars and gravitational wave observations.

The observation of the inspiral and merger of compact binaries by the LIGO-Virgo collaboration has allowed for new tests of Einstein's theory in the extreme gravity regime, where gravitational interactions are simultaneously strong, non-linear, and dynamical. Theories beyond Einstein's can also be constrained by detecting the polarization modes of gravitational waves. In this paper, we show that dynamical Chern-Simons and Einstein-dilaton-Gauss-Bonnet gravity cannot be differentiated from general relativity based on the detection of polarization modes alone. To prove this result, we use the Newman-Penrose method and an irreducible decomposition method to find that only the tensorial modes can be detected in both these theories.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

Electromagnetic observations have been used over the past decades to understand the nature of black holes and the material around them. Our ability to learn about the fundamental physics relies on our understanding of two key ingredients in the modeling of these electromagnetic observations: the gravity theory that describes the black hole, and the astrophysics that produces the observed radiation. In this work we study our current ability to constrain and detect deviations from General Relativity using the accretion disk spectrum of stellar-mass black holes in binary systems. Our analysis combines relativistic ray-tracing and Markov-Chain Monte-Carlo sampling techniques to determine how well such tests of General Relativity can be carried out in practice. We show that even when a very simple astrophysical model for the accretion disk is assumed a priori, the uncertainties and covariances between the parameters of the model and the parameters that control the deformation from General Relativity make any test of General Relativity very challenging with accretion disk spectrum observations. We also discuss the implications of assuming that General Relativity is correct a priori on the estimation of parameters of the astrophysical model when the data is not described by Einstein's theory, which can lead to a fundamental systematic bias.

The recent detection of gravitational waves by the LIGO/VIRGO collaboration has allowed for the first tests of Einstein's theory in the extreme gravity regime, where the gravitational interaction is simultaneously strong, non-linear and dynamical. One such test concerns the rate at which the binaries inspiral, or equivalently the rate at which the gravitational wave frequency increases, which can constrain the existence of hairy black holes. This is because black holes with scalar hair typically excite dipole radiation, which in turn leads to a faster decay rate and frequency chirping. In this paper, we present a mathematical proof that scalar hair is not sourced in vacuum, spherically symmetric spacetimes when considering extensions of Einstein's theory that break parity in gravity, focusing on dynamical Chern-Simons theory as a particular toy model. This result implies that the observational confirmation of the baldness of black holes cannot be used to constrain parity violation in gravity, unless the black holes observed are also spinning.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

The observation of the x-ray pulse profile emitted by hotspots on the surface of neutron stars offers a unique tool to measure the bulk properties of these objects, including their masses and radii. The x-ray emission takes place at the star's surface, where the gravitational field is strong, making these observations an incise probe to examine the curvature of spacetime generated by these stars. Motivated by this and the upcoming data releases by x-ray missions, such as NICER (Neutron star Interior Composition Explorer), we present a complete toolkit to model pulse profiles of rotating neutron stars in scalar-tensor gravity. We find that in this class of theories the presence of the scalar field affects the pulse profile's overall shape, producing strong deviations from the General Relativity expectation. This finding opens the possibility of potentially using x-ray pulse profile data to obtain new constraints on scalar-tensor gravity, if the pulse profile is found to be in agreement with General Relativity.

Gravitational wave astrophysics has only just begun, and as current detectors are upgraded and new detectors are built, many new, albeit faint, features in the signals will become accessible. One such feature is the presence of time-dependent Doppler shifts, generated by the acceleration of the center of mass of the gravitational-wave emitting system. We here develop a generic method that takes a frequency-domain, gravitational-wave model devoid of Doppler shifts and introduces modifications that incorporate them. Building upon a perturbative expansion that assumes the Doppler-shift velocity is small relative to the speed of light, the method consists of the inclusion of a single term in the Fourier phase and two terms in the Fourier amplitude. We validate the method through matches between waveforms with a Doppler shift in the time domain and waveforms constructed with our method for two toy problems: constant accelerations induced by a distant third body and Gaussian accelerations that resemble a kick profile. We find mismatches below $\sim\!10^{-6}$ for all of the astrophysically relevant cases considered, and improve further at smaller velocities. The work presented here will allow for the use of future detectors to extract new, faint features in the signal from the noise.

Future detectors could be used as a gravitational microscope to probe the horizon structure of merging black holes with gravitational waves. But can this microscope probe the quantum regime? We study this interesting question and find that (i) the error in the distance resolution is exponentially sensitive to errors in the Love number, and (ii) the uncertainty principle of quantum gravity forces a fundamental resolution limit. Thus, although the gravitational microscope can distinguish between black holes and other exotic objects, it is resolution limited well above the Planckian scale.

One of largest uncertainties in nuclear physics is the relation between the pressure and density of supranuclear matter: the equation of state. Some of this uncertainty may be removed through future gravitational wave observations of neutron star binaries by extracting the tidal deformabilities (or Love numbers) of neutron stars. Previous studies showed that only a certain combination of the individual deformabilities of each body (chirp tidal deformability) can be measured with second-generation gravitational wave interferometers, such as Adv. LIGO, due to correlations between the individual deformabilities. To overcome this, we search for approximately universal (or equation-of-state independent) relations between two combinations of the individual tidal deformabilities, such that once one of them has been measured, the other can be automatically obtained and the individual ones decoupled through these relations. We find an approximately universal relation between the symmetric and the anti-symmetric combination of the individual tidal deformabilities that is equation-of-state-insensitive to $20\%$ for binaries with masses less than $1.7M_\odot$. We show that these relations can be used to eliminate a combination of the tidal parameters from the list of model parameters, thus breaking degeneracies and improving the accuracy in parameter estimation. A simple study shows that the universal relations can improve the accuracy in the extraction of the symmetric combination of tidal parameters by as much as an order of magnitude, making the overall accuracy in the extraction of this parameter slightly better than that of the chirp tidal deformability. These new universal relations and the improved measurement accuracy on tidal parameters not only are important to astrophysics and nuclear physics, but also impact our ability to probe extreme gravity with gravitational waves and cosmology.

Neutron stars exhibit a set of universal relations independent of their equation of state that bears semblance to the black hole no hair relations. Motivated by this, we analytically and numerically explore other relations that connect neutron star and black hole universality. By analyzing two different measures, we find that certain rescaled entropies possess a nearly universal behavior. We also discover that when the compactness of neutron stars approaches the black hole limit, the rescaled entropy approaches that of a black hole, and the thermodynamic entropy scales with the stellar surface area in an ever more universal way.

The observation of the shadows cast by the event horizon of black holes on the light emitted in its neighborhood is the target of current very-long-baseline-interferometric observations. When considering supermassive black holes, the light source is the black hole's accretion disk, and therefore, the observation of the shadow may reveal information about the black hole and the accretion flow. We here study the shadows cast by stellar-mass black holes that are illuminated not by an accretion disk but by a stellar companion in a wide binary orbit. We calculate the shadows produced in such a configuration for the first time and show snapshots of the time-dependent shadow "movie" that is generated. We also study the minimal criteria for detecting and resolving such shadows with very-long-baseline-interferometric observations. We find that one would need telescopes capable of resolving apparent magnitudes greater than $33$ with baselines larger than $10^{6}$ km. Therefore, although very-long-baseline-interferometric efforts such as the Event Horizon Telescope would not be able to detect these shadows, their detection is possible with future telescopes in the next few decades.

Observations of the black hole shadow of supermassive black holes, such as Sagittarius A* at the center of our Milky Way galaxy, allow us to study the properties of black holes and the nature of strong-field gravity. According to the Kerr hypothesis, isolated, stationary, and axisymmetric astrophysical black holes are described by the Kerr metric. The Kerr hypothesis holds in General Relativity and in some modified gravity theories, but there are others in which it is violated. In principle, black hole shadow observations can be used to determine if the Kerr metric is the correct description for black holes, and in turn, they could be used to place constraints on modified gravity theories that do not admit the Kerr solution. We here investigate whether black hole shadow observations can constrain deviations from general relativity, focusing on two well-motivated modified quadratic gravity theories: Einstein-dilaton-Gauss-Bonnet gravity and dynamical Chern-Simons gravity. We show that current constraints on Einstein-dilaton-Gauss-Bonnet gravity are stronger than any that could be placed with black hole shadow observations of supermassive black holes. We also show that the same holds for dynamical Chern-Simons gravity through a systematic bias and a likelihood analysis when considering slowly-rotating supermassive black holes. However, observations of more rapidly-rotating black holes, with dimensionless spins $|\vec{J}|/M^{2}\simeq0.5$, could be used to better constrain dynamical Chern-Simons gravity.

The recent gravitational wave observations provide insight into the extreme gravity regime of coalescing binaries, where gravity is strong, dynamical and non-linear. The interpretation of these observations relies on the comparison of the data to a gravitational wave model, which in turn depends on the orbital evolution of the binary, and in particular on its orbital energy and angular momentum decay. In this paper, we calculate the latter in the inspiral of a non-spinning compact binary system within Einstein-\AE{}ther theory. From the theory's gravitational wave stress energy tensor and a balance law, we compute the angular momentum decay both as a function of the fields in the theory and as a function of the multipole moments of the binary. We then specialize to a Keplerian parameterization of the orbit to express the angular momentum decay as a function of the binary's orbital elements. We conclude by combining this with the orbital energy decay to find expressions for the decay of the semi-major axis and the orbital eccentricity of the binary. We find that these rates of decay are typically faster in Einstein-\AE{}ther theory than in General Relativity due to the presence of dipole radiation. Such modifications will imprint onto the chirp rate of gravitational waves, leaving a signature of Einstein-\AE{}ther theory that if absent in the data could be used to stringently constrain it.

The LIGO/Virgo detections of binary black hole mergers marked a watershed moment in astronomy, ushering in the era of precision tests of Kerr dynamics. We review theoretical and experimental challenges that must be overcome to carry out black hole spectroscopy with present and future gravitational wave detectors. Among other topics, we discuss quasinormal mode excitation in binary mergers, astrophysical event rates, tests of black hole dynamics in modified theories of gravity, parameterized "post-Kerr" ringdown tests, exotic compact objects, and proposed data analysis methods to improve spectroscopic tests of Kerr dynamics by stacking multiple events.

The observation of the inspiral and merger of compact binaries by the LIGO/Virgo collaboration ushered in a new era in the study of strong-field gravity. We review current and future tests of strong gravity and of the Kerr paradigm with gravitational-wave interferometers, both within a theory-agnostic framework (the parametrized post-Einsteinian formalism) and in the context of specific modified theories of gravity (scalar-tensor, Einstein-dilaton-Gauss-Bonnet, dynamical Chern-Simons, Lorentz-violating, and extra dimensional theories). In this contribution we focus on (i) the information carried by the inspiral radiation, and (ii) recent progress in numerical simulations of compact binary mergers in modified gravity.

Is gravity parity violating? Given the recent observations of gravitational waves from coalescing compact binaries, we develop a strategy to find an answer with current and future detectors. We identify the key signatures of parity violation in gravitational waves: amplitude birefringence in their propagation and a modified chirping rate in their generation. We then determine the optimal binaries to test the existence of parity violation in gravity, and prioritize the research in modeling that will be required to carry out such tests before detectors reach their design sensitivity.

The exterior gravitational field of a slowly-rotating neutron star can be characterized by its multipole moments, the first few being the neutron star mass, moment of inertia, and quadrupole moment to quadratic order in spin. In principle, all of these quantities depend on the neutron star's internal structure, and thus, on unknown nuclear physics at supra-nuclear energy densities. We here find relations between the moment of inertia, the Love numbers and the quadrupole moment (I-Love-Q relations) that do not depend sensitively on the neutron star's internal structure. Three important consequences derive from these I-Love-Q relations. On an observational astrophysics front, the measurement of a single member of the I-Love-Q trio would automatically provide information about the other two, even when the latter may not be observationally accessible. On a gravitational wave front, the I-Love-Q relations break the degeneracy between the quadrupole moment and the neutron-star spins in binary inspiral waveforms, allowing second-generation ground-based detectors to determine the (dimensionless) averaged spin to $\mathcal{O}(10)%$, given a sufficiently large signal-to-noise ratio detection. On a fundamental physics front, the I-Love-Q relations allow for tests of General Relativity in the neutron-star strong-field that are both theory- and internal structure-independent. As an example, by combining gravitational-wave and electromagnetic observations, one may constrain dynamical Chern-Simons gravity in the future by more than 6 orders of magnitude more stringently than Solar System and table-top constraints.

The coalescence of compact objects is one of the most promising sources of gravitational waves for ground-based interferometric detectors, such as advanced LIGO and Virgo. Generically, com- pact objects in binaries are expected to be spinning with spin angular momenta misaligned with the orbital angular momentum, causing the orbital plane to precess. This precession adds rich structure to the gravitational waves, introducing such complexity that an analytic closed-form description has been unavailable until now. We here construct the first closed-form frequency- domain gravitational waveforms that are valid for generic spin-precessing quasicircular compact binary inspirals. We first construct time-domain gravitational waves by solving the post-Newtonian precession equations of motion with radiation reaction through multiple scale analysis. We then Fourier transform these time-domain waveforms with the method of shifted uniform asymptotics to obtain closed-form expressions for frequency-domain waveforms. We study the accuracy of these analytic, frequency-domain waveforms relative to waveforms obtained by numerically evolving the post-Newtonian equations of motion and find that they are suitable for unbiased parameter estimation for 99.2%(94.6%) of the binary configurations we studied at a signal-to-noise ratio of 10(25). These new frequency-domain waveforms could be used for detection and parameter estimation studies due to their accuracy and low computational cost.

Neutron stars are ideal to probe, not only nuclear physics, but also strong-field gravity. Approximate universal relations insensitive to the star's internal structure exist among certain observables and are useful in testing General Relativity, as they project out the uncertainties in the equation of state. One such set of universal relations between the moment of inertia $(I)$, the tidal Love number and the quadrupole moment $(Q)$ has been studied both in General Relativity and in modified theories. In this paper, we study the relations in dynamical Chern-Simons gravity, a well-motivated, parity-violating effective field theory, extending previous work in various ways. First, we study how projected constraints on the theory using the I-Love relation depend on the measurement accuracy of $I$ with radio observations and that of the Love number with gravitational-wave observations. Provided these quantities can be measured with future observations, we find that the latter could place bounds on dynamical Chern-Simons gravity that are six orders of magnitude stronger than current bounds. Second, we study the I-Q and Q-Love relations in this theory by constructing slowly-rotating neutron star solutions to quadratic order in spin. We find that the approximate universality continues to hold in dynamical Chern-Simons gravity, and in fact, it becomes stronger than in General Relativity, although its existence depends on the normalization of the dimensional coupling constant of the theory. Finally, we study the variation of the eccentricity of isodensity contours inside a star and its relation to the degree of universality. We find that, in most cases, the eccentricity variation is smaller in dynamical Chern-Simons gravity than in General Relativity, providing further support to the idea that the approximate self-similarity of isodensity contours is responsible for universality.

Certain bulk properties of neutron stars, in particular their moment of inertia, rotational quadrupole moment and tidal Love number, when properly normalized, are related to one another in a nearly equation of state independent way. The goal of this paper is to test these relations with extreme equations of state at supranuclear densities constrained to satisfy only a handful of generic, physically sensible conditions. By requiring that the equation of state be (i) barotropic and (ii) its associated speed of sound be real, we construct a piecewise function that matches a tabulated equation of state at low densities, while matching a stiff equation of state parametrized by its sound speed in the high-density region. We show that the I-Love-Q relations hold to 1 percent with this class of equations of state, even in the extreme case where the speed of sound becomes superluminal and independently of the transition density. We also find further support for the interpretation of the I-Love-Q relations as an emergent symmetry due to the nearly constant eccentricity of isodensity contours inside the star. These results reinforce the robustness of the I-Love-Q relations against our current incomplete picture of physics at supranuclear densities, while strengthening our confidence in the applicability of these relations in neutron star astrophysics.

The population of stellar-mass, compact object binaries that merge with non-negligible eccentricity may be large enough to motivate searches with ground-based gravitational wave detectors. Such events could be exceptional laboratories to test General Relativity in the dynamical, strong-field regime, as a larger fraction of the energy is emitted at high-velocities, compared to quasi-circular inspirals. A serious obstacle here, however, is the challenge of computing theoretical waveforms for eccentric systems with the requisite accuracy for use in a matched-filter search. The corresponding waveforms are more a sequence of concentrated bursts of energy emitted near periapse than a continuous waveform. Based on this, an alternative approach, stacking excess power over the set of time-frequency tiles coincident with the bursts, was recently suggested as a more practical (though sub-optimal) detection strategy. The leading-order "observable" that would be inferred from such a detection would be a sequence of discrete numbers characterizing the position and size of each time-frequency tile. In General Relativity, this (possibly large) sequence of numbers is uniquely determined by the small set of parameters describing the binary at formation. In this work, following the spirit of the parameterized post-Einsteinian framework developed for quasi-circular inspiral, we propose a simple, parameterized deformation of the baseline general relativistic burst algorithm for eccentric inspiral events that would allow for model-independent tests of Einstein's theory in this high-velocity, strong-field regime.

Binary pulsars are ideal to test the foundations of General Relativity, such as Lorentz symmetry, which requires that experiments produce the same results in all free-falling (i.e.inertial) frames. We here break this symmetry in the gravitational sector by specifying a preferred time direction, and thus a preferred frame, at each spacetime point. We then examine the consequences of this gravitational Lorentz symmetry breaking in the orbital evolution of binary pulsars, focusing on the dissipative effects. We find that Lorentz symmetry breaking modifies these effects, and thus the orbital dynamics, in two different ways. First, it generically causes the emission of dipolar radiation, which makes the orbital separation decrease faster than in General Relativity. Second, the quadrupole component of the emission is also modified. The orbital evolution depends critically on the sensitivities of the stars, which measure how their binding energies depend on the motion relative to the preferred frame. We calculate the sensitivities numerically and compute the predicted orbital decay rate of binary pulsars in Lorentz-violating gravity. By testing these predictions against observations, we place very stringent constraints on gravitational Lorentz violation.

We derive a stationary and axisymmetric black hole solution to quadratic order in the spin angular momentum. The previously found, linear-in-spin terms modify the odd-parity sector of the metric, while the new corrections appear in the even-parity sector. These corrections modify the quadrupole moment, as well as the (coordinate-dependent) location of the event horizon and the ergoregion. Although the linear-in-spin metric is of Petrov type D, the quadratic order terms render it of type I. The metric does not possess a second-order Killing tensor or a Carter-like constant. The new metric does not possess closed timelike curves or spacetime regions that violate causality outside of the event horizon. The new, even-parity modifications to the Kerr metric decay less rapidly at spatial infinity than the leading-order in spin, odd-parity ones, and thus, the former are more important when considering black holes that are rotating moderately fast. We calculate the modifications to the Hamiltonian, binding energy and Kepler's third law. These modifications are crucial for the construction of gravitational wave templates for black hole binaries, which will enter at second post-Newtonian order, just like dissipative modifications found previously.

Observations of the continuum spectrum emitted by accretion disks around black holes allows us to infer their properties, including possibly whether black holes are described by the Kerr metric. Some modified gravity theories do not admit the Kerr metric as a solution, and thus, continuum spectrum observations could be used to constrain these theories. We here investigate whether current and next generation X-Ray observations of the black hole continuum spectrum can constrain such deviations from Einstein's theory, focusing on two well-motivated modified quadratic gravity theories: dynamical Chern-Simons gravity and Einstein-dilaton-Gauss-Bonnet gravity. We do so by determining whether the non-Kerr deviations in the continuum spectrum introduced by these theories are larger than the observational error intrinsic to the observations. We find that dynamical Chern-Simons gravity cannot be constrained better than current bounds with current or next generation continuum spectrum observations. Einstein-dilaton-Gauss-Bonnet gravity, however, may be constrained better than current bounds with next generation telescopes, as long as the systematic error inherent in the accretion disk modeling is decreased below the predicted observational error.

Binary systems of two compact objects circularize and spiral toward each other via the emission of gravitational waves. The coupling of the spins of each object with the orbital angular momentum causes the orbital plane to precess, which leads to modulation of the gravitational wave signal. Until now, generating frequency-domain waveforms for fully precessing systems for use in gravitational wave data analysis meant numerically integrating the equations of motion, then Fourier transforming the result, which is very computationally intensive for systems that complete hundreds or thousands of cycles in the sensitive band of a detector. Previously, analytic solutions were only available for certain special cases or for simplified models. Here we describe the construction of closed-form, frequency-domain waveforms for fully-precessing, quasi-circular binary inspirals.

Rapidly rotating black hole solutions in theories beyond general relativity play a key role in experimental gravity, as they allow us to compute observables in extreme spacetimes that deviate from the predictions of general relativity. Such solutions are often difficult to find in beyond-general-relativity theories due to the inclusion of additional fields that couple to the metric non-linearly and non-minimally. In this paper, we consider rotating black hole solutions in one such theory, dynamical Chern-Simons gravity, where the Einstein-Hilbert action is modified by the introduction of a dynamical scalar field that couples to the metric through the Pontryagin density. We treat dynamical Chern-Simons gravity as an effective field theory and work in the decoupling limit, where corrections are treated as small perturbations from general relativity. We perturb about the maximally-rotating Kerr solution, the so-called extremal limit, and develop mathematical insight into the analysis techniques needed to construct solutions for generic spin. First we find closed-form, analytic expressions for the extremal scalar field, and then determine the trace of the metric perturbation, giving both in terms of Legendre decompositions. Retaining only the first three and four modes in the Legendre representation of the scalar field and the trace, respectively, suffices to ensure a fidelity of over 99% relative to full numerical solutions. The leading-order mode in the Legendre expansion of the trace of the metric perturbation contains a logarithmic divergence at the extremal Kerr horizon, which is likely to be unimportant as it occurs inside the perturbed dynamical Chern-Simons horizon. The techniques employed here should enable the construction of analytic, closed-form expressions for the scalar field and metric perturbations on a background with arbitrary rotation.

The gravitational wave observations GW150914 and GW151226 by Advanced LIGO provide the first opportunity to learn about physics in the extreme gravity environment of coalescing binary black holes. The LIGO Scientific Collaboration and the Virgo Collaboration have verified that this observation is consistent with General Relativity. This paper expands their analysis to a larger class of anomalies, highlighting the inferences that can be drawn on non-standard theoretical physics mechanisms. We find that these events constrain a plethora of mechanisms associated with the generation and propagation of gravitational waves, including the activation of scalar fields, gravitational leakage into large extra dimensions, the variability of Newton's constant, a modified dispersion relation, gravitational Lorentz violation and the strong equivalence principle. Though other observations limit many of these mechanisms already, GW150914 and GW151226 are unique in that they are direct probes of dynamical strong-field gravity and of gravitational wave propagation. We also show that GW150914 constrains inferred properties of exotic compact object alternatives to Kerr black holes. We argue, however, that the true potential for GW150914 to both rule out exotic objects and constrain physics beyond General Relativity is severely limited by the lack of understanding of the merger regime in almost all relevant modified gravity theories. This event thus significantly raises the bar that these theories have to pass, both in terms of having a sound theoretical underpinning, and being able to solve the equations of motion for binary merger events. We conclude with a discussion of the additional inferences that can be drawn if the lower-confidence observation of an electromagnetic counterpart to GW150914 holds true; this would provide dramatic constraints on the speed of gravity and gravitational Lorentz violation.

Neutron stars and quark stars are ideal laboratories to study fundamental physics at supra nuclear densities and strong gravitational fields. Astrophysical observables, however, depend strongly on the star's internal structure, which is currently unknown due to uncertainties in the equation of state. Universal relations, however, exist among certain stellar observables that do not depend sensitively on the star's internal structure. One such set of relations is between the star's moment of inertia ($I$), its tidal Love number (Love) and its quadrupole moment ($Q$), the so-called I-Love-Q relations. Similar relations hold among the star's multipole moments, which resemble the well-known black hole no-hair theorems. Universal relations break degeneracies among astrophysical observables, leading to a variety of applications: (i) X-ray measurements of the nuclear matter equation of state, (ii) gravitational wave measurements of the intrinsic spin of inspiraling compact objects, and (iii) gravitational and astrophysical tests of General Relativity that are independent of the equation of state. We here review how the universal relations come about and all the applications that have been devised to date.

When in a tight binary, the mutual tidal deformations of neutron stars imprint onto observables, encoding information about their internal structure at supranuclear densities and gravity in the extreme-gravity regime. Gravitational wave observations of their late binary inspiral may serve as a tool to extract the individual tidal deformabilities, but this is made difficult by degeneracies between them in the gravitational wave model. We here resolve this problem by discovering approximately universal relations between dimensionless combinations of the individual tidal deformabilities. We show that these relations break degeneracies in the gravitational wave model, allowing for the accurate extraction of both deformabilities. Such measurements can be used to better differentiate between equation-of-state models, and improve tests of General Relativity and cosmology.

The relations between most observables associated with a compact star, such as the mass and radius of a neutron star or a quark star, typically depend strongly on their unknown internal structure. The I-Love-Q relations (between the moment of inertia, the tidal deformability and the quadrupole moment) are however approximately insensitive to this structure. These relations become exact for stationary black holes in General Relativity as shown by the no-hair theorems. In this paper, we take the first steps toward studying how the approximate I-Love-Q relations become exact in the limit as compact stars become black holes. To do so, we consider a toy model, i.e. incompressible stars with anisotropic pressure, which allows us to model an equilibrium sequence of stars with their compactness approaching the black hole limit arbitrarily closely. We extract the I-Love-Q trio by numerically constructing such a sequence in the slow-rotation and small-tide approximations. We find that the I-Love-Q relations approach the black hole limit in a nontrivial way, with the quadrupole moment and the tidal deformability changing sign as the compactness and the amount of anisotropy are increased. Generalizing Maclaurin spheroids to anisotropic stars, we show that the multipole moments also change sign in the Newtonian limit as the amount of anisotropy is increased. We also prove analytically that the stellar moment of inertia reaches the black hole limit as the compactness reaches the black hole value in the strongly anisotropic limit. Modeling the black hole limit through a sequence of anisotropic stars, however, fails when considering other theories of gravity. We calculate the scalar dipole charge and the moment of inertia in a parity-violating modified theory and find that these quantities do not tend to their black hole counterparts as the anisotropic stellar sequence approaches the black hole limit.

The aLIGO detection of the black-hole binary GW150914 opened a new era for probing extreme gravity. Many gravity theories predict the emission of dipole gravitational radiation by binaries. This is excluded to high accuracy in binary pulsars, but entire classes of theories predict this effect predominantly (or only) in binaries involving black holes. Joint observations of GW150914-like systems by aLIGO and eLISA will improve bounds on dipole emission from black-hole binaries by six orders of magnitude relative to current constraints, provided that eLISA is not dramatically descoped.

Future electromagnetic observations of black holes may allow us to test General Relativity in the strong-field regime. Such tests, however, require knowledge of rotating black hole solutions in modified gravity theories, a class of which does not admit the Kerr metric as a solution. Several rotating black hole solutions in modified theories have only been found in the slow-rotation approximation (i.e. assuming the spin angular momentum is much smaller than the mass squared). We here investigate whether the systematic error due to the approximate nature of these black hole metrics is small enough relative to the observational error to allow their use in electromagnetic observations to constrain deviations from General Relativity. We address this by considering whether electromagnetic observables constructed from a slow-rotation approximation to the Kerr metric can fit observables constructed from the full Kerr metric with systematic errors smaller than current observational errors. We focus on black hole shadow and continuum spectrum observations, as these are the least influenced by accretion disk physics, with current observational errors of about 10%. We find that the fractional systematic error introduced by using a second-order, slowly-rotating Kerr metric is at most 2% for shadows created by black holes with dimensionless spins $\chi\leq0.6$. We also find that the systematic error introduced by using the slowly-rotating Kerr metric as an exact metric when constructing continuum spectrum observables is negligible for black holes with dimensionless spins of $\chi \lesssim 0.9$. Our results suggest that the modified gravity solutions found in the slow-rotation approximation may be used to constrain realistic deviations from General Relativity with continuum spectrum and black hole shadow observations.

We study isolated and binary neutron stars in dynamical Chern-Simons gravity. This theory modifies the Einstein-Hilbert action through the introduction of a dynamical scalar field coupled to the Pontryagin density. We here treat this theory as an effective model, working to leading order in the Chern-Simons coupling. We first construct isolated neutron star solutions in the slow-rotation expansion to quadratic order in spin. We find that isolated neutron stars acquire a scalar dipole charge that corrects its spin angular momentum to linear order in spin and corrects its mass and quadrupole moment to quadratic order in spin, as measured by an observer at spatial infinity. We then consider neutron stars binaries that are widely separated and solve for their orbital evolution in this modified theory. We find that the evolution of post-Keplerian parameters is modified, with the rate of periastron advance being the dominant correction at first post-Newtonian order. We conclude by applying these results to observed pulsars with the aim to place constraints on dynamical Chern-Simons gravity. We find that the modifications to the observed mass are degenerate with the neutron star equation of state, which prevents us from testing the theory with the inferred mass of the millisecond pulsar J1614-2230. We also find that the corrections to the post-Keplerian parameters are too small to be observable today even with data from the double binary pulsar J0737-3039. Our results suggest that pulsar observations are not currently capable of constraining dynamical Chern-Simons gravity, and thus, gravitational-wave observations may be the only path to a stringent constraint of this theory in the imminent future.

Corrections to general relativity that introduce long-ranged scalar fields which are non-minimally coupled to curvature typically predict that neutron stars possess a non-trivial scalar field profile. An observer far from a star is most sensitive to the spherically-symmetric piece of this profile that decays linearly with the inverse of the distance, the so-called scalar charge, which is related to the emission of dipolar radiation from compact binaries. The presence of dipolar radiation has the potential to very strongly constrain extended theories of gravity. These facts may lead people to believe that gravitational theories with long-ranged scalar fields have already been constrained strongly from binary pulsar observations. Here we challenge this "lore" by investigating the decoupling limit of Gauss-Bonnet gravity as an example, in which the scalar field couples linearly to the Gauss-Bonnet density in the action. We prove a theorem that neutron stars in this theory cannot possess a scalar charge. Thus Gauss-Bonnet gravity evades the strong binary pulsar constraints on dipole radiation. We discuss the astrophysical systems yielding the best constraints on Gauss-Bonnet gravity and related quadratic gravity theories. To achieve this we explicitly compute the scalar charge for slowly-rotating neutron stars in quadratic gravity theories. In generic case, either neutron star-binary or neutron star-black hole systems can be used to constrain the theory, but Gauss-Bonnet gravity evades the neutron star-binary constraints. However, black holes in Gauss-Bonnet gravity do anchor scalar charge. The best constraints on Gauss-Bonnet gravity will thus come from black hole observations, for example via pulsar-black hole binaries. We estimate these constraints to be ten times better than the current estimated bound, and also include estimated constraints on generic quadratic gravity theories.

We consider a general class of quantum gravity-inspired, modified gravity theories, where the Einstein-Hilbert action is extended through the addition of all terms quadratic in the curvature tensor coupled to scalar fields with standard kinetic energy. This class of theories includes Einstein-Dilaton-Gauss-Bonnet and Chern-Simons modified gravity as special cases. We analytically derive and solve the coupled field equations in the post-Newtonian approximation, assuming a comparable-mass, spinning black hole binary source in a quasi-circular, weak-field/slow-motion orbit. We find that a naive subtraction of divergent piece associated with the point-particle approximation is ill-suited to represent compact objects in these theories. Instead, we model them by appropriate effective sources built so that known strong-field solutions are reproduced in the far-field limit. In doing so, we prove that black holes in Einstein-Dilaton-Gauss-Bonnet and Chern-Simons theory can have hair, while neutron stars have no scalar monopole charge, in diametrical opposition to results in scalar-tensor theories. We then employ techniques similar to the direct integration of the relaxed Einstein equations to obtain analytic expressions for the scalar field, metric perturbation, and the associated gravitational wave luminosity measured at infinity. We find that scalar field emission mainly dominates the energy flux budget, sourcing electric-type (even-parity) dipole scalar radiation and magnetic-type (odd-parity) quadrupole scalar radiation, correcting the General Relativistic prediction at relative -1PN and 2PN orders. Such modifications lead to corrections in the emitted gravitational waves that can be mapped to the parameterized post-Einsteinian framework. Such modifications could be strongly constrained with gravitational wave observations.

Gravitational waves from neutron star binary inspirals contain information about the equation of state of supranuclear matter. In the absence of definitive experimental evidence that determines the correct equation of state, a number of diverse models that give the pressure in a neutron star as function of its density have been proposed. These models differ not only in the approximations and techniques they use to solve the many-body Schr\"odinger equation, but also in the neutron star composition they assume. We study whether gravitational wave observations of neutron star binaries in quasicircular inspirals will allow us to distinguish between equations of state of differing internal composition, thereby providing important information about the properties of extremely high density matter. We carry out a Bayesian model selection analysis, and find that second generation gravitational wave detectors can heavily constrain equations of state that contain only quark matter, but hybrid stars containing both normal and quark matter are harder to distinguish from normal matter stars. A gravitational wave detection with a signal-to-noise ratio of 30 and masses around $1.4M_{\odot}$ could either detect or rule out strange quark stars with a 20 to 1 confidence. The presence of kaon condensates or hyperons in neutron star inner cores cannot be easily confirmed. For example, for the equations of state studied in this paper, even a gravitational wave signal with a signal-to-noise ratio as high as 60 would not allow us to claim a detection of kaon condensates or hyperons with confidence greater than 5 to 1. On the other hand, if kaon condensates and hyperons do not form in neutron stars, a gravitational wave signal with similar signal-to-noise ratio would be able to constrain their existence with an 11 to 1 confidence for high-mass systems.

The opportunity to study physics at supra-nuclear densities through X-ray observations of neutron stars has led to in-depth investigations of certain approximately universal relations that can remove degeneracies in pulse profile models. One such set of relations determines all of the multipole moments of a neutron star just from the first three (the mass monopole, the current dipole and the mass quadrupole moment) approximately independently of the equation of state. These three-hair relations were found to hold in neutron stars that rotate rigidly, as is the case in old pulsars, but neutron stars can also rotate differentially, as is the case for proto-neutron stars and hypermassive transient remnants of binary mergers. We here extend the three-hair relations to differentially rotating stars for the first time with a generic rotation law using two approximations: a weak-field scheme (an expansion in powers of the neutron star compactness) and a perturbative differential rotation scheme (an expansion about rigid rotation). These approximations allow us to analytically derive approximately universal relations that allow us to determine all of the multipole moments of a (perturbative) differentially rotating star in terms of only the first four moments. These new four-hair relations for differentially rotating neutron stars are found to be approximately independent of the equation of state to a higher degree than the three-hair relations for uniformly rotating stars. Our results can be instrumental in the development of four-hair relations for rapidly differentially rotating stars in full General Relativity using numerical simulations.

Certain physical quantities that characterize neutron stars and quark stars (e.g. their mass, spin angular momentum and quadrupole moment) are interrelated in a way that is approximately insensitive to their internal structure. Such approximately universal relations are useful to break degeneracies in data analysis for future radio, X-ray and gravitational wave observations. Although the pressure inside compact stars is most likely nearly isotropic, certain scenarios have been put forth that suggest otherwise, for example due to phase transitions. We here investigate whether pressure anisotropy affects the approximate universal relations and whether it prevents their use in future observations. We achieve this by numerically constructing slowly-rotating and tidally-deformed, anisotropic, compact stars in General Relativity to third order in spin. We find that anisotropy affects the universal relations only weakly; the relations become less universal by a factor of 1.5-3 relative to the isotropic case, but remain approximately universal to 10%. We succeed in explaining this increase in variability as an increase in the eccentricity variation of isodensity contours, which provides further support for the emergent approximate symmetry explanation of universality. Anisotropy does not affect the universal relations to a sufficient level to prevent their use in gravitational wave astrophysics or in experimental relativity. We provide an explicit example of the latter in dynamical Chern-Simons gravity. The increase in variability of the universal relations due to pressure anisotropy could affect their use in future X-ray observations. Given expected observational uncertainties, however, the relations remain sufficiently universal for use in such observations if the anisotropic modifications to the moment of inertia and the quadrupole moment are less than 10% of their isotropic values.

Compact stars satisfy certain no-hair relations through which their multipole moments are given by their mass, spin and quadrupole moment. These relations are approximately independent of their equation of state, relating pressure to density. Such relations are similar to the black hole no-hair theorems, but these possess event horizons inside which information that led to their formation can hide. Compact stars do not possess horizons, so whether their no-hair relations are related to the black hole ones is unclear. We investigate how the two relations are related by studying relations among multipole moments for compact stars with anisotropic pressure as a toy model, which allows such stars to be more compact than those with isotropic pressure. We here show numerically that the compact star no-hair relations approach the black hole ones as the compactness approaches that of a black hole. We also prove analytically that the current dipole moment exactly reaches the black hole limit quadratically in compactness as strongly-anisotropic stars approach the black hole limit. We moreover show that compact stars become progressively oblate in this limit, even if prolate at low compactness due to strong anisotropies.

No-hair like relations between the multipole moments of the exterior gravitational field of neutron stars have recently been found to be approximately independent of the star's internal structure. This approximate, equation-of-state universality arises after one adimensionalizes the multipole moments appropriately, which then begs the question of whether there are better ways to adimensionalize the moments to obtain stronger universality. We here investigate this question in detail by considering slowly-rotating neutron stars both in the non-relativistic limit and in full General Relativity. We find that there exist normalizations that lead to stronger equation-of-state universality in the relations among the moment of inertia and the quadrupole, octopole and hexadecapole moments of neutron stars. We determine the optimal normalization that minimizes the equation-of-state dependence in these relations. The results found here may have applications in the modeling of X-ray pulses and atomic line profiles from millisecond pulsars with NICER and LOFT.

Mergers of compact stellar remnant are prime targets for the LIGO/Virgo gravitational wave detectors. One hopes that the gravitational wave signals from these merger events can be used to study the mass and spin distribution of stellar remnants, and provide information about black hole horizons and the material properties of neutron stars. However, it has been suggested that degeneracies in the way that the star's mass and spin are imprinted in the waveforms may make it impossible to distinguish between black holes and neutron stars. Here we show that the precession of the orbital plane due to spin-orbit coupling breaks the mass-spin degeneracy, and allows us to distinguish between standard neutron stars and alternative possibilities, such as black holes or exotic neutron stars with large masses and spins.

The detection and characterization of gravitational wave signals from compact binary coalescence events relies on accurate waveform templates in the frequency domain. The stationary phase approximation (SPA) can be used to compute closed-form frequency-domain waveforms for non-precessing, quasi-circular binary inspirals. However, until now, no fast frequency-domain waveforms have existed for generic, spin-precessing quasi-circular compact binary inspirals. Templates for these systems have had to be computed via a discrete Fourier transform of finely sampled time-domain waveforms, which is far more computationally expensive than those constructed directly in the frequency-domain, especially for those systems that are dominated by the inspiral part. There are two obstacles to deriving frequency-domain waveforms for precessing systems: (i) the spin-precession equations do not admit closed-form solutions for generic systems; (ii) the SPA fails catastrophically. Presently there is no general solution to the first problem, so we must resort to numerical integration of the spin precession equations. This is not a significant obstacle, as numerical integration on the slow precession timescale adds very little to the computational cost of generating the waveforms. Our main result is to solve the second problem, by providing an alternative to the SPA that we call the method of Shifted Uniform Asymptotics, or SUA, that cures the divergences in the SPA. The construction of frequency-domain templates using the SUA can be orders of magnitude more efficient than the time-domain ones obtained through a discrete Fourier transform. Moreover, this method is very faithful to the discrete Fourier transform, with mismatches on the order of $10^{-5}$.

Modified gravity theories generically predict a violation of Lorentz invariance, which may lead to a modified dispersion relation for propagating modes of gravitational waves. We construct a parametrized dispersion relation that can reproduce a range of known Lorentz-violating predictions and investigate their impact on the propagation of gravitational waves. A modified dispersion relation forces different wavelengths of the gravitational wave train to travel at slightly different velocities, leading to a modified phase evolution observed at a gravitational-wave detector. We show how such corrections map to the waveform observable and to the parametrized post-Einsteinian framework, proposed to model a range of deviations from General Relativity. Given a gravitational-wave detection, the lack of evidence for such corrections could then be used to place a constraint on Lorentz violation. The constraints we obtain are tightest for dispersion relations that scale with small power of the graviton's momentum and deteriorate for a steeper scaling.

Compact binaries that emit gravitational waves in the sensitivity band of ground-based detectors can have non-negligible eccentricities just prior to merger, depending on the formation scenario. We develop a purely analytic, frequency-domain model for gravitational waves emitted by compact binaries on orbits with small eccentricity, which reduces to the quasi-circular post-Newtonian approximant TaylorF2 at zero eccentricity and to the post-circular approximation of Yunes et al. (2009) at small eccentricity. Our model uses a spectral approximation to the (post-Newtonian) Kepler problem to model the orbital phase as a function of frequency, accounting for eccentricity effects up to ${\cal{O}}(e^8)$ at each post-Newtonian order. Our approach accurately reproduces an alternative time-domain eccentric waveform model for eccentricities $e\in [0, 0.4]$ and binaries with total mass less than 12 solar masses. As an application, we evaluate the signal amplitude that eccentric binaries produce in different networks of existing and forthcoming gravitational waves detectors. Assuming a population of eccentric systems containing black holes and neutron stars that is uniformly distributed in co-moving volume, we estimate that second generation detectors like Advanced LIGO could detect approximately 0.1-10 events per year out to redshift $z\sim 0.2$, while an array of Einstein Telescope detectors could detect hundreds of events per year to redshift $z \sim 2.3$.

Astrophysical charge-free black holes are known to satisfy no-hair relations through which all multipole moments can be specified in terms of just their mass and spin angular momentum. We here investigate the possible existence of no-hair-like relations among multipole moments for neutron stars and quark stars that are independent of their equation of state. We calculate the multipole moments of these stars up to hexadecapole order by constructing uniformly-rotating and unmagnetized stellar solutions to the Einstein equations. For slowly-rotating stars, we construct stellar solutions to quartic order in spin in a slow-rotation expansion, while for rapidly-rotating stars, we solve the Einstein equations numerically with the LORENE and RNS codes. We find that the multipole moments extracted from these numerical solutions are consistent with each other. We confirm that the current-dipole is related to the mass-quadrupole in an approximately equation of state independent fashion, which does not break for rapidly rotating neutron stars or quark stars. We further find that the current-octupole and the mass-hexadecapole moments are related to the mass-quadrupole in an approximately equation of state independent way to $\sim 10%$, worsening in the hexadecapole case. All of our findings are in good agreement with previous work that considered stellar solutions to leading-order in a weak-field expansion. The quartic in spin, slowly-rotating solutions found here allow us to estimate the systematic errors in the measurement of the neutron star's mass and radius with future X-ray observations, such as NICER and LOFT. We find that the effect of these quartic-in-spin terms on the quadrupole and hexadecapole moments and stellar eccentricity may dominate the error budget for very rapidly-rotating neutron stars. The new universal relations found here should help to reduce such systematic errors.

Black holes are said to have no hair because all of their multipole moments can be expressed in terms of just their mass, charge and spin angular momentum. The recent discovery of approximately equation-of-state-independent relations among certain multipole moments in neutron stars suggests that they are also approximately bald. We here explore the yet unknown origin for this universality. First, we investigate which region of the neutron star's interior and of the equation of state is most responsible for the universality. We find that the universal relation between the moment of inertia and the quadrupole moment is dominated by the star's outer-core, a shell of width (50-95)% of the total radius, which corresponds to the density range $(10^{14}$-$10^{15})$g/cm$^3$. Second, we study the impact on the universality of approximating stellar isodensity contours as self-similar ellipsoids. An analytical calculation in the Newtonian limit reveals that the shape of the ellipsoids does not affect the universal relations, but relaxing the self-similarity assumption can completely destroy it. Third, we investigate the eccentricity profiles of rotating relativistic stars and find that the ellipticity is roughly constant, with variations of roughly (20-30)% in the region that matters to the universal relations. Fourth, we repeat the above analysis for non-compact, regular stars and find that the ellipticity is not constant, with variations that easily exceed 100% and universality is lost. These findings suggest that universality arises as an emergent approximate symmetry: as one flows in the stellar-structure phase space from non-compact star region to the relativistic star region, the eccentricity variation inside stars decreases, leading to approximate self-similarity in their isodensity contours, which then leads to the universal behavior observed in their exterior multipole moments.

The gravitational field outside of astrophysical black holes is completely described by their mass and spin frequency, as expressed by the no-hair theorems. These theorems assume vacuum spacetimes, and thus they apply only to black holes and not to stars. Despite this, we analytically find that the gravitational potential of arbitrarily rapid rigidly rotating stars can still be described completely by only their mass, spin angular momentum, and quadrupole moment. Although these results are obtained in the nonrelativistic limit (to leading order in a weak-field expansion of general relativity, GR), they are also consistent with fully relativistic numerical calculations of rotating neutron stars. This description of the gravitational potential outside the source in terms of just three quantities is approximately universal (independent of equation of state). Such universality may be used to break degeneracies in pulsar and future gravitational wave observations to extract more physics and test GR in the strong-field regime.

We lay the foundations for the construction of analytic expressions for Fourier-domain gravitational waveforms produced by eccentric, inspiraling compact binaries in a post-circular or small-eccentricity approximation. The time-dependent, "plus" and "cross" polarizations are expanded in Bessel functions, which are then self-consistently re-expanded in a power series about zero initial eccentricity to eighth order. The stationary phase approximation is then employed to obtain explicit analytic expressions for the Fourier transform of the post-circular expanded, time-domain signal. We exemplify this framework by considering Newtonian-accurate waveforms, which in the post-circular scheme give rise to higher harmonics of the orbital phase and amplitude corrections both to the amplitude and the phase of the Fourier domain waveform. Such higher harmonics lead to an effective increase in the inspiral mass reach of a detector as a function of the binary's eccentricity e_0 at the time when the binary enters the detector sensitivity band. Using the largest initial eccentricity allowed by our approximations (e_0 < 0.4), the mass reach is found to be enhanced up to factors of approximately 5 relative to that of circular binaries for Advanced LIGO, LISA, and the proposed Einstein Telescope at a signal-to-noise ratio of ten. A post-Newtonian generalization of the post circular scheme is also discussed, which holds the promise to provide "ready-to-use" Fourier-domain waveforms for data analysis of eccentric inspirals.

Binary pulsars are excellent laboratories to test the building blocks of Einstein's theory of General Relativity. One of these is Lorentz symmetry which states that physical phenomena appear the same for all inertially moving observers. We study the effect of violations of Lorentz symmetry in the orbital evolution of binary pulsars and find that it induces a much more rapid decay of the binary's orbital period due to the emission of dipolar radiation. The absence of such behavior in recent observations allows us to place the most stringent constraints on Lorentz violation in gravity, thus verifying one of the cornerstones of Einstein's theory much more accurately than any previous gravitational observation.

The main goal of a scientific workshop is to bring together experts in a specific field or related fields to collaborate, to discuss, and to creatively make progress in a particular area. The organizational aspects of such a meeting play a critical role in achieving these goals. We here present suggestions from scientists to scientists that hopefully help in organizing a successful scientific workshop that maximizes collaboration and creativity.

Waveform templates are a powerful tool for extracting and characterizing gravitational wave signals, acting as highly restrictive priors on the signal morphologies that allow us to extract weak events buried deep in the instrumental noise. The templates map the waveform shapes to physical parameters, thus allowing us to produce posterior probability distributions for these parameters. However, there are attendant dangers in using highly restrictive signal priors. If strong field gravity is not accurately described by General Relativity (GR), then using GR templates may result in fundamental bias in the recovered parameters, or even worse, a complete failure to detect signals. Here we study such dangers, concentrating on three distinct possibilities. First, we show that there exist modified theories compatible with all existing tests that would fail to be detected by the LIGO/Virgo network using searches based on GR templates, but which would be detected using a one parameter post-Einsteinian extension. Second, we study modified theories that produce departures from GR that turn on suddenly at a critical frequency, producing waveforms that do not naively fit into the simplest parameterized post-Einsteinian (ppE) scheme. We show that even the simplest ppE templates are still capable of picking up these strange signals and diagnosing a departure from GR. Third, we study whether using inspiral-only ppE waveforms for signals that include merger and ringdown can lead to problems in misidentifying a GR departure. We present an easy technique that allows us to self-consistently identify the inspiral portion of the signal, and thus remove these potential biases, allowing GR tests to be performed on higher mass signals that merge within the detector band. We close by studying a parameterized waveform model that may allow us to test GR using the full inspiral-merger-ringdown signal.

The waveform phase for a neutron star binary can be split into point-particle terms and finite-size terms (characterized by the Love number) that account for equation of state effects. The latter first enter at 5 post-Newtonian (PN) order (i.e. proportional to the tenth power of the orbital velocity), but the former are only known completely to 3.5 PN order, with higher order terms only known to leading-order in the mass-ratio. We here find that not including point-particle terms at 4PN order to leading- and first-order in the mass ratio in the template model can severely deteriorate our ability to measure the equation of state. This problem can be solved if one uses numerical waveforms once their own systematic errors are under control.

The Newman-Janis algorithm has been widely used to construct rotating black hole solutions from non-rotating counterparts. While this algorithm was developed within General Relativity, it has more recently been applied to non-rotating solutions in modified gravity theories. We find that the application of the Newman-Janis algorithm to an arbitrary non-GR spherically-symmetric solution introduces pathologies in the resulting axially-symmetric metric. This then establishes that, in general, the Newman-Janis algorithm should not used to construct rotating black hole solutions outside of General Relativity.

We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.

Neutron stars and quark stars are not only characterized by their mass and radius, but also by how fast they spin, through their moment of inertia, and how much they can be deformed, through their Love number and quadrupole moment. These depend sensitively on the star's internal structure, and thus on unknown nuclear physics. We find universal relations between the moment of inertia, the Love number and the quadrupole moment that are independent of the neutron star's and quark star's internal structure. These can be used to learn about the deformability of these compact objects through observations of the moment of inertia, break degeneracies in gravitational wave detection to measure spin in binary inspirals and test General Relativity in a nuclear-structure independent fashion.

This review is focused on tests of Einstein's theory of General Relativity with gravitational waves that are detectable by ground-based interferometers and pulsar timing experiments. Einstein's theory has been greatly constrained in the quasi-linear, quasi-stationary regime, where gravity is weak and velocities are small. Gravitational waves will allow us to probe a complimentary, yet previously unexplored regime: the non-linear and dynamical strong-field regime. Such a regime is, for example, applicable to compact binaries coalescing, where characteristic velocities can reach fifty percent the speed of light and compactnesses can reach a half. This review begins with the theoretical basis and the predicted gravitational wave observables of modified gravity theories. The review continues with a brief description of the detectors, including both gravitational wave interferometers and pulsar timing arrays, leading to a discussion of the data analysis formalism that is applicable for such tests. The review ends with a discussion of gravitational wave tests for compact binary systems.

We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a given direction splits into two modes, one subluminal and one superluminal. The width of the splitting depends on the azimuthal components of the propagation vector, is linearly proportional to the mass of the black hole, and decreases with the third inverse power of the distance from the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable amplitude is excited, and vice-versa.

Several model-independent parameterizations of deviations from General Relativity have been developed to test Einstein's theory. Although these different parameterizations were developed for different gravitational observables, they ultimately all test the same underlying physics. In this paper, we develop connections between the parameterized post-Newtonian, parameterized post-Keplerian, and the parameterized post-Einsteinian frameworks, developed to carry out tests of General Relativity with Solar System, binary pulsar, and gravitational wave observations respectively. These connections allow us to use knowledge gained from one framework to inform and guide tests using the others. Relating these parameterizations and combining the results from each approach strengthens our tests of General Relativity.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We calculate analytical gravitational waveforms in the time- and frequency-domain for precessing quasi-circular binaries with spins of arbitrary magnitude, but nearly aligned with the orbital angular momentum. We first derive an analytical solution to the precession equations by expanding in the misalignment angle and using multiple scale analysis to separate timescales. We then use uniform asymptotic expansions to analytically Fourier transform the time-domain waveform, thus extending the stationary-phase approximation, which fails when precession is present. The resulting frequency-domain waveform family has a high overlap with numerical waveforms obtained by direct integration of the post-Newtonian equations of motion and discrete Fourier transformations. Such a waveform family lays the foundations for the accurate inclusion of spin precession effects in analytical gravitational waveforms, and thus, it can aid in the detection and parameter estimation of gravitational wave signals from the inspiral phase of precessing binary systems.

We study generic tests of strong-field General Relativity using gravitational waves emitted during the inspiral of compact binaries. Previous studies have considered simple extensions to the standard post-Newtonian waveforms that differ by a single term in the phase. Here we improve on these studies by (i) increasing the realism of injections and (ii) determining the optimal waveform families for detecting and characterizing such signals. We construct waveforms that deviate from those in General Relativity through a series of post-Newtonian terms, and find that these higher-order terms can affect our ability to test General Relativity, in some cases by making it easier to detect a deviation, and in some cases by making it more difficult. We find that simple single-phase post-Einsteinian waveforms are sufficient for detecting deviations from General Relativity, and there is little to be gained from using more complicated models with multiple phase terms. The results found here will help guide future attempts to test General Relativity with advanced ground-based detectors.

A stellar-mass compact object spiraling into a supermassive black hole, an extreme-mass-ratio inspiral (EMRI), is one of the targets of future gravitational-wave detectors and it offers a unique opportunity to test General Relativity (GR) in the strong-field. We study whether generic scalar-tensor (ST) theories can be further constrained with EMRIs. We show that in the EMRI limit, all such theories universally reduce to massive or massless Brans-Dicke theory and that black holes do not emit dipolar radiation to all orders in post-Newtonian (PN) theory. For massless theories, we calculate the scalar energy flux in the Teukolsky formalism to all orders in PN theory and fit it to a high-order PN expansion. We derive the PN ST corrections to the Fourier transform of the gravitational wave response and map it to the parameterized post-Einsteinian framework. We use the effective-one-body framework adapted to EMRIs to calculate the ST modifications to the gravitational waveform. We find that such corrections are smaller than those induced in the early inspiral of comparable-mass binaries, leading to projected bounds on the coupling that are worse than current Solar System ones. Brans-Dicke theory modifies the weak-field, with deviations in the energy flux that are largest at small velocities. For massive theories, superradiance can lead to resonances in the scalar energy flux that can lead to floating orbits outside the innermost stable circular orbit and that last until the supermassive black hole loses enough mass and spin-angular momentum. If such floating orbits occur in the frequency band of LISA, they would lead to a large dephasing (~1e6 rads), preventing detection with GR templates. A detection that is consistent with GR would then rule out floating resonances at frequencies lower than the lowest observed frequency, allowing for the strongest constraints yet on massive ST theories.

As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.

We describe a new kludge scheme to model the dynamics of generic extreme-mass-ratio inspirals (EMRIs; stellar compact objects spiraling into a spinning supermassive black hole) and their gravitational-wave emission. The Chimera scheme is a hybrid method that combines tools from different approximation techniques in General Relativity: (i) A multipolar, post-Minkowskian expansion for the far-zone metric perturbation (the gravitational waveforms) and for the local prescription of the self-force; (ii) a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and (iii) a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The EMRI trajectory is made out of Kerr geodesic fragments joined via the method of osculating elements as dictated by the multipolar post-Minkowskian radiation-reaction prescription. We implemented the proper coordinate mapping between Boyer-Lindquist coordinates, associated with the Kerr geodesics, and harmonic coordinates, associated with the multipolar post-Minkowskian decomposition. The Chimera scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme to intermediate mass ratios, and hence, it can provide valuable information for future space-based gravitational-wave observatories, like LISA, and even for advanced ground detectors. The local character in time of our multipolar post-Minkowskian self-force makes this scheme amenable to study the possible appearance of transient resonances in generic inspirals.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

We introduce a new kludge scheme to model the dynamics of generic extreme mass-ratio inspirals (stellar compact objects spiraling into a spinning supermassive black hole) and to produce the gravitational waveforms that describe the gravitational-wave emission of these systems. This scheme combines tools from different techniques in General Relativity: It uses a multipolar, post-Minkowskian (MPM) expansion for the far-zone metric perturbation (which provides the gravitational waveforms, here taken up to mass hexadecapole and current octopole order) and for the local prescription of the self-force (since we are lacking a general prescription for it); a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The orbital evolution is thus equivalent to solving the geodesic equations with time-dependent orbital elements, as dictated by the MPM radiation-reaction prescription. To complete the scheme, both the orbital evolution and wave generation require to map the Boyer-Lindquist coordinates of the orbits to the harmonic coordinates in which the different MPM quantities have been derived, a mapping that we provide explicitly in this paper. This new kludge scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme mass ratios to systems with more moderate mass ratios, and hence can provide valuable information for future space-based gravitational-wave observatories like LISA and even for advanced ground detectors. Finally, due to the local character in time of our MPM self-force, this scheme can be used to perform studies of the possible appearance of transient resonances in generic inspirals.

Gravitational wave observations will probe non-linear gravitational interactions and thus enable strong tests of Einstein's theory of general relativity. We present a numerical relativity study of the late inspiral and merger of binary black holes in scalar-tensor theories of gravity. We consider black hole binaries in an inhomogeneous scalar field, specifically binaries inside a scalar field bubble, in some cases with a potential. We calculate the emission of dipole radiation. We also show how these configurations trigger detectable differences between gravitational waves in scalar-tensor gravity and the corresponding waves in general relativity. We conclude that, barring an external mechanism to induce dynamics in the scalar field, scalar-tensor gravity binary black holes alone are not capable of awaking a dormant scalar field, and are thus observationally indistinguishable from their general relativistic counterparts.

The future detection of gravitational wave forces us to consider the many ways in which astrophysics, gravitational wave theory and fundamental theory will interact. In this paper, I summarize some recent work done to develop such an interface. In particular, I concentrate on how non-vacuum astrophysical environments can modify the gravitational wave signal emitted by compact binary inspirals, and whether signatures from the former are detectable by current and future gravitational wave detectors. I also describe the interface between gravitational wave modeling and fundamental theory, focusing on the status of the parameterized post-Einsteinian framework (a general framework to detect deviations away from General Relativity in future gravitational wave data) and its current data analysis implementation.

Gravitational wave astronomy has tremendous potential for studying extreme astrophysical phenomena and exploring fundamental physics. The waves produced by binary black hole mergers will provide a pristine environment in which to study strong field, dynamical gravity. Extracting detailed information about these systems requires accurate theoretical models of the gravitational wave signals. If gravity is not described by General Relativity, analyses that are based on waveforms derived from Einstein's field equations could result in parameter biases and a loss of detection efficiency. A new class of "parameterized post-Einsteinian" (ppE) waveforms has been proposed to cover this eventuality. Here we apply the ppE approach to simulated data from a network of advanced ground based interferometers (aLIGO/aVirgo) and from a future spaced based interferometer (LISA). Bayesian inference and model selection are used to investigate parameter biases, and to determine the level at which departures from general relativity can be detected. We find that in some cases the parameter biases from assuming the wrong theory can be severe. We also find that gravitational wave observations will beat the existing bounds on deviations from general relativity derived from the orbital decay of binary pulsars by a large margin across a wide swath of parameter space.

We study the effect of black hole spin on the accuracy of the post-Newtonian approximation. We focus on the gravitational energy flux for the quasicircular, equatorial, extreme mass-ratio inspiral of a compact object into a Kerr black hole of mass M and spin J. For a given dimensionless spin a=J/M^2 (in geometrical units), the energy flux depends only on the orbital velocity v or (equivalently) on the Boyer-Lindquist orbital radius r. We investigate the formal region of validity of the Taylor post-Newtonian expansion of the energy flux (which is known up to order v^8 beyond the quadrupole formula), generalizing previous work by two of us. The "error function" used to determine the region of validity of the post-Newtonian expansion can have two qualitatively different kinds of behavior, and we deal with these two cases separately. We find that, at any fixed post-Newtonian order, the edge of the region of validity (as measured by v/v_{ISCO}, where v_{ISCO} is the orbital velocity at the innermost stable circular orbit) is only weakly dependent on a. Unlike in the nonspinning case, the lack of sufficiently high order terms does not allow us to determine if there is a convergent to divergent transition at order v^6. Independently of a, the inclusion of angular multipoles up to and including l=5 in the numerical flux is necessary to achieve the level of accuracy of the best-known (N=8) PN expansion of the energy flux.

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW-driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk leaves an imprint on the GW through its angular momentum exchange with the binary, the mass increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that observations of EMRIs may place new sensitive constraints on the physics of accretion disks.

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

We study two large classes of alternative theories, modifying the action through algebraic, quadratic curvature invariants coupled to scalar fields. We find one class that admits solutions that solve the vacuum Einstein equations and another that does not. In the latter, we find a deformation to the Schwarzschild metric that solves the modified field equations in the small coupling approximation. We calculate the event horizon shift, the innermost stable circular orbit shift, and corrections to gravitational waves, mapping them to the parametrized post-Einsteinian framework.

Extreme mass ratio inspirals, in which a stellar-mass object merges with a supermassive black hole, are prime sources for space-based gravitational wave detectors because they will facilitate tests of strong gravity and probe the spacetime around rotating compact objects. In the last few years of such inspirals, the total phase is in the millions of radians and details of the waveforms are sensitive to small perturbations. We show that one potentially detectable perturbation is the presence of a second supermassive black hole within a few tenths of a parsec. The acceleration produced by the perturber on the extreme mass-ratio system produces a steady drift that causes the waveform to deviate systematically from that of an isolated system. If the perturber is a few tenths of a parsec from the extreme-mass ratio system (plausible in as many as a few percent of cases) higher derivatives of motion might also be detectable. In that case, the mass and distance of the perturber can be derived independently, which would allow a new probe of merger dynamics.

The inspiral of binary systems in vacuum is controlled by the stress-energy of gravitational radiation and any other propagating degrees of freedom. For gravitational waves, the dominant contribution is characterized by an effective stress-energy tensor at future null infinity. We employ perturbation theory and the short-wavelength approximation to compute this stress-energy tensor in a wide class of alternative theories. We find that this tensor is generally a modification of that first computed by Isaacson, where the corrections can dominate over the general relativistic term. In a wide class of theories, however, these corrections identically vanish at asymptotically flat, future, null infinity, reducing the stress-energy tensor to Isaacson's. We exemplify this phenomenon by first considering dynamical Chern-Simons modified gravity, which corrects the action via a scalar field and the contraction of the Riemann tensor and its dual. We then consider a wide class of theories with dynamical scalar fields coupled to higher-order curvature invariants, and show that the gravitational wave stress-energy tensor still reduces to Isaacson's. The calculations presented in this paper are crucial to perform systematic tests of such modified gravity theories through the orbital decay of binary pulsars or through gravitational wave observations.

The emergent area of gravitational wave astronomy promises to provide revolutionary discoveries in the areas of astrophysics, cosmology, and fundamental physics. One of the most exciting possibilities is to use gravitational-wave observations to test alternative theories of gravity. In this contribution we describe how to use observations of extreme-mass-ratio inspirals by the future Laser Interferometer Space Antenna to test a particular class of theories: Chern-Simons modified gravity.

One of the most interesting high-energy, astrophysical phenomena are relativistic jets emitted from highly localized sky location. Such jets are common in Nature, observed to high redshift and in a range of wavelengths. Their precise generation mechanism remains a bit of a mystery, but they are generically believed to be powered by black holes. We here summarize the recent simulations of Palenzuela, Lehner and Liebling that shed light on the jet generation mechanism. These authors studied the merger of two non-spinning black holes in the presence of a magnetic field, perpendicular to the orbital plane and anchored by a circumbinary accretion disk, in the "force-free" approximation. They found that each black hole essentially acts as a "straw" that stirs the magnetic field lines around the center of mass as the black holes inspiral. The twisting of the magnetic field lines then generates jets around each black hole, even though these are not spinning. Their simulations show the formation of such a dual jet geometry and how it transitions to a single jet one, as the black holes merge due to gravitational wave emission.

We constrain the parameterized post-Einsteinian framework with binary pulsar observations of orbital period decay due to gravitational wave emission. This framework proposes to enhance the amplitude and phase of gravitational waveform templates through post-Einsteinian parameters to search for generic deviations from General Relativity in gravitational wave data. Such enhancements interpolate between General Relativity and alternative theory predictions, but their magnitude must be such as to satisfy all current experiments and observations. The data that currently constrains the parameterized post-Einsteinian framework the most is the orbital period decay of binary pulsars. We use such observations to place upper limits on the magnitude of post-Einsteinian parameters, which will be critical when gravitational waves are detected and this framework is implemented.

Gravitational parity violation is a possibility motivated by particle physics, string theory and loop quantum gravity. One effect of it is amplitude birefringence of gravitational waves, whereby left and right circularly-polarized waves propagate at the same speed but with different amplitude evolution. Here we propose a test of this effect through coincident observations of gravitational waves and short gamma-ray bursts from binary mergers involving neutron stars. Such gravitational waves are highly left or right circularly-polarized due to the geometry of the merger. Using localization information from the gamma-ray burst, ground-based gravitational wave detectors can measure the distance to the source with reasonable accuracy. An electromagnetic determination of the redshift from an afterglow or host galaxy yields an independent measure of this distance. Gravitational parity violation would manifest itself as a discrepancy between these two distance measurements. We exemplify such a test by considering one specific effective theory that leads to such gravitational parity-violation, Chern-Simons gravity. We show that the advanced LIGO-Virgo network and all-sky gamma-ray telescopes can be sensitive to the propagating sector of Chern-Simons gravitational parity violation to a level roughly two orders of magnitude better than current stationary constraints from the LAGEOS satellites.

Gravitational waves are messengers carrying valuable information about their sources. For sources at cosmological distances, the waves will contain also the imprint left by the intervening matter. The situation is in close analogy with cosmic microwave photons, for which the large-scale structures the photons traverse contribute to the observed temperature anisotropies, in a process known as the integrated Sachs-Wolfe effect. We derive the gravitational wave counterpart of this effect for waves propagating on a Friedman-Robertson-Walker background with scalar perturbations. We find that the phase, frequency and amplitude of the gravitational waves experience Sachs-Wolfe type integrated effects, this in addition to the magnification effects on the amplitude from gravitational lensing. We show that for supermassive black hole binaries, the integrated effects could account for measurable changes on the frequency, chirp mass and luminosity distance of the binary, thus unveiling the presence of inhomogeneities, and potentially dark energy, in the Universe.

Neutron stars are sensitive laboratories for testing general relativity, especially when considering deviations where velocities are relativistic and gravitational fields are strong. One such deviation is described by dynamical, Chern-Simons modified gravity, where the Einstein-Hilbert action is modified through the addition of the gravitational parity-violating Pontryagin density coupled to a field. This four-dimensional effective theory arises naturally both in perturbative and non-perturbative string theory, loop quantum gravity, and generic effective field theory expansions. We calculate here Chern-Simons modifications to the properties and gravitational fields of slowly spinning neutron stars. We find that the Chern-Simons correction affects only the gravitomagnetic sector of the metric to leading order, thus introducing modifications to the moment of inertia but not to the mass-radius relation. We show that an observational determination of the moment of inertia to an accuracy of 10%, as is expected from near-future observations of the double pulsar, will place a constraint on the Chern-Simons coupling constant of \xi^{1/4} < 5 km, which is at least three-orders of magnitude stronger than the previous strongest bound.

Space-borne gravitational wave detectors, such as the proposed Laser Interferometer Space Antenna, are expected to observe black hole coalescences to high redshift and with large signal-to-noise ratios, rendering their gravitational waves ideal probes of fundamental physics. The promotion of Newton's constant to a time-function introduces modifications to the binary's binding energy and the gravitational wave luminosity, leading to corrections in the chirping frequency. Such corrections propagate into the response function and, given a gravitational wave observation, they allow for constraints on the first time-derivative of Newton's constant at the time of merger. We find that space-borne detectors could indeed place interesting constraints on this quantity as a function of sky position and redshift, providing a {\emph{constraint map}} over the entire range of redshifts where binary black hole mergers are expected to occur. A LISA observation of an equal-mass inspiral event with total redshifted mass of 10^5 solar masses for three years should be able to measure $\dot{G}/G$ at the time of merger to better than 10^(-11)/yr.

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Nicol\'as Yunes (Princeton University) to extend in high detail his recent work on EMRI waveforms for our highlight article.

We present the first models of extreme-mass-ratio inspirals within the effective-one-body (EOB) formalism, focusing on quasi-circular orbits into non-rotating black holes. We show that the phase difference and (Newtonian normalized) amplitude difference between analytical EOB and numerical Teukolsky-based gravitational waveforms can be reduced to less than 10^(-1) rad and less than 2 x 10^(-3), respectively, after a 2-year evolution. The inclusion of post-Newtonian self-force terms in the EOB approach leads to a phase disagreement of roughly 6-27 rad after a 2-year evolution. Such inclusion could also allow for the EOB modeling of waveforms from intermediate-mass ratio, quasi-circular inspirals.

We consider the concept of fundamental bias in gravitational wave astrophysics as the assumption that general relativity is the correct theory of gravity during the entire wave-generation and propagation regime. Such an assumption is valid in the weak field, as verified by precision experiments and observations, but it need not hold in the dynamical strong-field regime where tests are lacking. Fundamental bias can cause systematic errors in the detection and parameter estimation of signals, which can lead to a mischaracterization of the universe through incorrect inferences about source event rates and populations. We propose a remedy through the introduction of the parameterized post-Einsteinian framework, which consists of the enhancement of waveform templates via the inclusion of post-Einsteinian parameters. These parameters would ostensibly be designed to interpolate between templates constructed in general relativity and well-motivated alternative theories of gravity, and also include extrapolations that follow sound theoretical principles, such as consistency with conservation laws and symmetries. As an example, we construct parameterized post-Einsteinian templates for the binary coalescence of equal-mass, non-spinning compact objects in a quasi-circular inspiral. The parametrized post-Einsteinian framework should allow matched filtered data to select a specific set of post-Einsteinian parameters without a priori assuming the validity of the former, thus either verifying general relativity or pointing to possible dynamical strong-field deviations.

Testing deviation of GR is one of the main goals of the proposed {\emph{Laser Interferometer Space Antenna}}, a space-based gravitational-wave observatory. For the first time, we consistently compute the generation of gravitational waves from extreme-mass ratio inspirals (stellar compact objects into supermassive black holes) in a well-motivated alternative theory of gravity, that to date remains weakly constrained by double binary pulsar observations. The theory we concentrate on is Chern-Simons (CS) modified gravity, a 4-D, effective theory that is motivated both from string theory and loop-quantum gravity, and which enhances the Einstein-Hilbert action through the addition of a dynamical scalar field and the parity-violating Pontryagin density. We show that although point particles continue to follow geodesics in the modified theory, the background about which they inspiral is a modification to the Kerr metric, which imprints a CS correction on the gravitational waves emitted. CS modified gravitational waves are sufficiently different from the General Relativistic expectation that they lead to significant dephasing after 3 weeks of evolution, but such dephasing will probably not prevent detection of these signals, but instead lead to a systematic error in the determination of parameters. We end with a study of radiation-reaction in the modified theory and show that, to leading-order, energy-momentum emission is not CS modified, except possibly for the subdominant effect of scalar-field emission. The inclusion of radiation-reaction will allow for tests of CS modified gravity with space-borne detectors that might be two orders of magnitude larger than current binary pulsar bounds.

The low-energy limit of string theory contains an anomaly-canceling correction to the Einstein-Hilbert action, which defines an effective theory: Chern-Simons (CS) modified gravity. The CS correction consists of the product of a scalar field with the Pontryagin density, where the former can be treated as a background field (non-dynamical formulation) or as an evolving field (dynamical formulation). Many solutions of general relativity persist in the modified theory; a notable exception is the Kerr metric, which has sparked a search for rotating black hole solutions. Here, for the first time, we find a solution describing a rotating black hole within the dynamical framework, and in the small-coupling/slow-rotation limit. The solution is axisymmetric and stationary, constituting a deformation of the Kerr metric with dipole scalar "hair," whose effect on geodesic motion is to weaken the frame-dragging effect and shift the location of the inner-most stable circular orbit outwards (inwards) relative to Kerr for co-rotating (counter-rotating) geodesics. We further show that the correction to the metric scales inversely with the fourth power of the radial distance to the black hole, suggesting it will escape any meaningful bounds from weak-field experiments. For example, using binary pulsar data we can only place an initial bound on the magnitude of the dynamical coupling constant of $\xi^{1/4} \lesssim 10^{4} {\textrm{km}}$. More stringent bounds will require observations of inherently strong-field phenomena.

Chern-Simons modified gravity is an effective extension of general relativity that captures leading-order, gravitational parity violation. Such an effective theory is motivated by anomaly cancelation in particle physics and string theory. In this review, we begin by providing a pedagogical derivation of the three distinct ways such an extension arises: (1) in particle physics, (2) from string theory and (3) geometrically. We then review many exact and approximate, vacuum solutions of the modified theory, and discuss possible matter couplings. Following this, we review the myriad astrophysical, solar system, gravitational wave and cosmological probes that bound Chern-Simons modified gravity, including discussions of cosmic baryon asymmetry and inflation. The review closes with a discussion of possible future directions in which to test and study gravitational parity violation.

Chern-Simons modified gravity is a string-theory and loop-quantum-gravity inspired effective theory that modifies General Relativity by adding a parity-violating Pontryagin density to the Einstein-Hilbert action multiplied by a coupling scalar. We strongly constrain non-dynamical Chern-Simons modified gravity with a timelike Chern-Simons scalar through observations of the double binary pulsar PSR J0737-3039A/B. We first calculate Chern-Simons corrections to the orbital evolution of binary systems. We find that the ratio of the correction to periastron precession to the general relativistic prediction scales quadratically with the semi-major axis and inversely with the square of the object's radius. Binary pulsar systems are thus ideal to test this theory, since periastron precession can be measured with sub-degree accuracies and the semi-major axis is millions of times larger than the stellar radius. Using data from PSR J0737-3039A/B we dramatically constrain the non-dynamical Chern-Simons coupling to $M_{\CS} := 1/|\dot{\vartheta}| > 33 {\textrm{meV}}$, approximately a hundred billion times better than current Solar System tests.

[abridged] Chern-Simons (CS) modified gravity is a 4D effective theory that descends both from string theory and loop quantum gravity, and that corrects the Einstein-Hilbert action by adding the product of a scalar field and the parity-violating, Pontryagin density. In this theory, the gravitational field of spinning black holes is described by a modified Kerr geometry whose multipole moments deviate from the Kerr ones only at the fourth multipole, l = 4. We investigate possible signatures of this theory in the gravitational wave emission produced in the inspiral of stellar compact objects into massive black holes, both for intermediate- and extreme-mass ratios. We use the semi-relativistic approximation, where the trajectories are geodesics of the massive black hole geometry and the gravitational waveforms are obtained from a multipolar decomposition of the radiative field. The main CS corrections to the waveforms arise from modifications to the geodesic trajectories, due to changes to the massive black hole geometry, and manifest themselves as an accumulating dephasing relative to the general relativistic case. We also explore the propagation and the stress-energy tensor of gravitational waves in this theory. We find that, although this tensor has the same form as in General Relativity, the energy and angular momentum balance laws are indeed modified through the stress-energy tensor of the CS scalar field. These balance laws could be used to describe the inspiral through adiabatic changes in the orbital parameters, which in turn would enhance the dephasing effect. Gravitational-wave observations of intermediate- or extreme-mass ratio inspirals with advanced ground detectors or with LISA could use such dephasing to test the dynamical theory to unprecedented levels.

We consider a loop-quantum gravity inspired modification of general relativity, where the Holst action is generalized by making the Barbero-Immirzi (BI) parameter a scalar field, whose value could be dynamically determined. The modified theory leads to a non-zero torsion tensor that corrects the field equations through quadratic first-derivatives of the BI field. Such a correction is equivalent to general relativity in the presence of a scalar field with non-trivial kinetic energy. This stress-energy of this field is automatically covariantly conserved by its own dynamical equations of motion, thus satisfying the strong equivalence principle. Every general relativistic solution remains a solution to the modified theory for any constant value of the BI field. For arbitrary time-varying BI fields, a study of cosmological solutions reduces the scalar field stress-energy to that of a pressureless perfect fluid in a comoving reference frame, forcing the scale factor dynamics to be equivalent to those of a stiff equation of state. Upon ultraviolet completion, this model could provide a natural mechanism for k-inflation, where the role of the inflaton is played by the BI field and inflation is driven by its non-trivial kinetic energy instead of a potential.

We study perturbations of a Schwarzschild black hole in Chern-Simons modified gravity. We begin by showing that Birkhoff's theorem holds for a wide family of Chern-Simons coupling functions, a scalar field present in the theory that controls the strength of the Chern-Simons correction to the Einstein-Hilbert action. After decomposing the perturbations in spherical harmonics, we study the linearized modified field equations and find that axial and polar modes are coupled, in contrast to general relativity. The divergence of the modified equations leads to the Pontryagin constraint, which forces the vanishing of the Cunningham-Price-Moncrief master function associated with axial modes. We analyze the structure of these equations and find that the appearance of the Pontryagin constraint yields an overconstrained system that does not allow for generic black hole oscillations. We illustrate this situation by studying the case characterized by a canonical choice of the coupling function and pure-parity perturbative modes. We end with a discussion of how to extend Chern-Simons modified gravity to bypass the Pontryagin constraint and the suppression of perturbations.

Extreme mass ratio bursts (EMRBs) have been proposed as a possible source for future space-borne gravitational wave detectors, such as the Laser Interferometer Space Antenna (LISA). These events are characterized by long-period, nearly-radial orbits of compact objects around a central massive black hole. The gravitational radiation emitted during such events consists of a short burst, corresponding to periapse passage, followed by a longer, silent interval. In this paper we investigate the impact of including relativistic corrections to the description of the compact object's trajectory via a geodesic treatment, as well as including higher-order multipole corrections in the waveform calculation. The degree to which the relativistic corrections are important depends on the EMRB's orbital parameters. We find that relativistic EMRBs (v_{max}}/c > 0.25) are not rare and actually account for approximately half of the events in our astrophysical model. The relativistic corrections tend to significantly change the waveform amplitude and phase relative to a Newtonian description, although some of this dephasing could be mimicked by parameter errors. The dephasing over several bursts could be of particular importance not only to gravitational wave detection, but also to parameter estimation, since it is highly correlated to the spin of the massive black hole. Consequently, we postulate that if a relativistic EMRB is detected, such dephasing might be used to probe the relativistic character of the massive black hole and obtain information about its spin.

We study Chern-Simons (CS) gravity in the parameterized post-Newtonian (PPN) framework through a weak-field solution of the modified field equations. We find that CS gravity possesses the same PPN parameters as general relativity, except for the inclusion of a new term, proportional to the CS coupling and the curl of the PPN vector potential. This new term leads to a modification of frame dragging and gyroscopic precession and we provide an estimate of its size. This correction might be used in experiments, such as Gravity Probe B, to bound CS gravity and test string theory.

During the inspiral and merger of a binary black hole, gravitational radiation is emitted anisotropically due to asymmetries in the merger configuration. This anisotropic radiation leads to a gravitational wave kick, or recoil velocity, as large as ~ 4000 km/sec. We investigate the effect gravitational recoil has on the retention of intermediate mass black holes (IMBH) within Galactic globular clusters. Assuming that our current understanding of IMBH-formation is correct and yields an IMBH-seed in every globular cluster, we find a significant problem retaining low mass IMBHs (1000 $\Msun$) in the typical merger-rich globular cluster environment. Given a uniform black hole spin distribution and orientation and a Kroupa IMF, we find that at most 3% of the globular clusters can retain an IMBH larger than 1000 $\Msun$ today. For a population of black holes that better approximates mass loss from winds and supernovae, we find that 16% of globulars can retain an IMBH larger than 1000 $\Msun$. Our calculations show that if there are black holes of mass $M > 60 \Msun$ in a cluster, repeated IMBH-BH encounters will eventually eject a 1000 $\Msun$ IMBH with greater than 30% probability. As a consequence, a large population of rogue black holes may exist in our Milky Way halo. We discuss the dynamical implications of this subpopulation, and its possible connection to ultraluminous X-ray sources (ULXs).

We investigate the weak-field, post-Newtonian expansion to the solution of the field equations in Chern-Simons gravity with a perfect fluid source. In particular, we study the mapping of this solution to the parameterized post-Newtonian formalism to 1 PN order in the metric. We find that the PPN parameters of Chern-Simons gravity are identical to those of general relativity, with the exception of the inclusion of a new term that is proportional to the Chern-Simons coupling parameter and the curl of the PPN vector potentials. We also find that the new term is naturally enhanced by the non-linearity of spacetime and we provide a physical interpretation for it. By mapping this correction to the gravito-electro-magnetic framework, we study the corrections that this new term introduces to the acceleration of point particles and the frame-dragging effect in gyroscopic precession. We find that the Chern-Simons correction to these classical predictions could be used by current and future experiments to place bounds on intrinsic parameters of Chern-Simons gravity and, thus, string theory.

[abridged] The coalescence of a binary black hole system is one of the main sources of gravitational waves that present and future detectors will study. Apart from the energy and angular momentum that these waves carry, for unequal-mass binaries there is also a net flux of linear momentum that implies a recoil velocity of the resulting final black hole in the opposite direction. We present a computation of the recoil velocity based on the close-limit approximation scheme, which gives excellent results for head-on and grazing collisions of black holes when compared to full numerical relativistic calculations. We obtain a maximum recoil velocity of ~ 64 km/s for a symmetric mass ratio eta = M_1 M_2/(M_1+M_2)^2 ~ 0.19 and an initial proper separation of 4 M, where M is the total ADM mass of the system. This separation is the maximum at which the close-limit approximation is expected to provide accurate results. If we supplement this estimate with PN calculations up to the innermost stable circular orbit, we obtain a lower bound for the recoil velocity, with a maximum around 84 km/s. This is a lower bound because it neglects the initial merger phase. We can however obtain a rough estimate by using PN methods or the close-limit approximation. Since both methods are known to overestimate the amount of radiation, we obtain in this way an upper bound for the recoil with maxima in the range of 220-265 km/s. We also provide non-linear fits to these estimated upper and lower bounds. These estimates are subject to uncertainties related to issues such as the choice of initial data and higher effects in perturbation theory. Nonetheless, our estimates are consistent with previous results in the literature and suggest a narrower range of possible recoil velocities.

The formation and growth of supermassive black holes is a key issue to unveil the secrets of galaxy formation. In particular, the gravitational recoil produced in the merger of unequal mass black hole binaries could have a number of astrophysical implications, such as the ejection of black holes from the host galaxy or globular cluster. We present estimates of the recoil velocity that include the effect of small eccentricities. The approach is specially suited for the last stage of the merger, where most of the emission of linear momentum in gravitational waves takes place. Supplementing our estimates with post-Newtonian approximations, we obtain lower and upper bounds that constrain previous recoil velocities estimates as well as a best estimate that agrees with numerical simulations in the quasi-circular case. For eccentricities e <= 0.1, the maximum recoil is found for mass ratios of M_1/M_2 ~ 0.38 with a best estimate of ~ 167 (1 + e) km/s and upper and lower bounds of 79 (1 + e) km/s and 216 (1 + e) km/s respectively.

Self-similarity induced by critical gravitational collapse is used as a paradigm to probe the mass distribution of subsolar objects. At large mass (solar mass and above) there is widespread agreement as to both the form and parameter values arising in the mass distribution of stellar objects. At subsolar mass there is still considerable disagreement as to the qualitative form of the mass distribution, let alone the specific parameter values characterizing that distribution. For the first time, the paradigm of critical gravitational collapse is applied to several concrete astrophysical scenarios to derive robust qualitative features of the subsolar mass distribution. We further contrast these theoretically derived ideas with the observational situation. In particular, we demonstrate that at very low mass the distribution is given by a power law, with an exponent opposite in sign to that observed in the high-mass regime. The value of this low-mass exponent is in principle calculable via dynamical systems theory applied to gravitational collapse. Qualitative agreement between theory, numerical experiments, and observational data is good, though quantitative issues remain troublesome.

Critical gravitational collapse and self similarity are used to probe the mass distribution of subsolar objects. We demonstrate that at very low mass the distribution is given by a power law, with an exponent opposite in sign to that observed at high-mass regime. We further show that the value of this low-mass exponent is in principle calculable via dynamical systems theory applied to gravitational collapse. Qualitative agreement between numerical experiments and observational data is good.