The extreme-gravity collisions between black holes allow us to probe the underlying theory of gravity. We apply the theory-agnostic inspiral-merger-ringdown consistency test to an example theory beyond general relativity for the first time. Here we focus on the string-inspired Einstein-dilaton Gauss-Bonnet gravity and modify the inspiral, ringdown, and remnant black hole properties of the gravitational waveform. We found that future multiband observations allow us to constrain the theory stronger than current observations by an order of magnitude. The formalism developed here can easily be applied to other theories.

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.

Pulsar timing observations precisely test general relativity. Recently, the hierarchical triple system PSR J0337+1715 has placed new constraints on the existence of a fifth force from violation in the strong equivalence principle. Many alternative gravity theories exist with massive (pseudo-)scalar fields to explain a variety of phenomena from the accelerating expansion of the universe at large scales to the QCD strong CP problem at small scales. We here develop a generic formalism for the fifth force effect in theories involving massive scalar fields arising from e.g. string theory. With PSR J0337 measurements, we find the strongest bound on the simplest theory with a massive scalar field beyond general relativity and derive new constraints in other theories with axions, dark matter mediators, and higher-curvature corrections. These results show that the triple system J0337 provides a stringent test for massive scalar fields.

Recent observations of neutron stars (NS) with radio, X-rays and gravitational waves have begun to constrain the equation of state (EoS) for nuclear matter beyond the nuclear saturation density. To one's surprise, there exist approximate universal relations connecting certain bulk properties of NSs that are insensitive to the underlying EoS and having important applications on probing fundamental physics including nuclear and gravitational physics. To date, analytic works on universal relations for realistic NSs are lacking, which may lead to a better understanding of the origin of the universality. Here, we focus on the universal relations between the compactness(C), moment of inertia (I), and tidal deformability(related to the Love number), and derive analytic, approximate I-Love-C relations. To achieve this, we construct slowly-rotating/ tidally-deformed NS solutions analytically starting from an extended Tolman VII model that accurately describes non-rotating realistic NSs, allowing us to extract the moment of inertia and the tidal deformability on top of the compactness. We solve the field equations analytically by expanding them about the Newtonian limit and keeping up to 6th order in the stellar compactness. Based on these analytic solutions, we can mathematically demonstrate the O(10%) EoS variation in the I-C and Love-C relations and the O(1%) variation in the I-Love relation that have previously been found numerically. Our new analytic relations agree more accurately with numerical results for realistic NSs(especially the I-C and Love-C ones) than the analytic relations for constant density stars derived in previous work. Based on these analytic findings, we attribute a possible origin of the universality for the I-C and Love-C relations to the fact that the energy density of realistic NSs can be approximated as a quadratic function, as is the case for the Tolman VII solution.

Gravitational waves from extreme gravity events such as the coalescence of two black holes in a binary system fill our observable universe, bearing with them the underlying theory of gravity driving their process. One compelling alternative theory of gravity -- known as Einstein-dilaton Gauss-Bonnet gravity motivated by string theory -- describes the presence of an additional dilaton scalar field coupled directly to higher orders of the curvature, effectively describing a "fifth force" interaction and the emission of scalar dipole radiation between two scalarized black holes. Most previous studies focused on considering only the leading correction to the inspiral portion of the binary black hole waveforms. In our recent paper, we carried out inspiral-merger-ringdown consistency tests in this string-inspired gravity by including corrections to both the inspiral and ringdown portions, as well as those to the mass and spin of remnant black holes, valid to quadratic order in spin. We here extend the analysis by directly computing bounds on the theoretical coupling constant using the full inspiral-merger-ringdown waveform rather than treating the inspiral and merger-ringdown portions separately. We also consider the corrections valid to quartic order in spin to justify the validity of black hole's slow-rotation approximation. We find the quasinormal mode corrections to the waveform to be particularly important for high-mass events such as GW170729, in which the dilaton fields' small-coupling approximation fails without such effects included. We also show that future space-based and multiband gravitational-wave observations have the potential to go beyond existing bounds on the theory. The bounds presented here are comparable to those found in via the inspiral-merger-ringdown consistency tests.

Gravitational waves from the explosive merger of distant black holes are encoded with details regarding the complex extreme-gravity spacetime present at their source. Famously described by the Kerr spacetime metric for rotating black holes in general relativity, what if effects beyond this theory are present? One way to efficiently test this hypothesis is to first obtain a metric which parametrically deviates from the Kerr metric in a model-independent way. Given such a metric, one can then predict the ensuing corrections to both the inspiral and ringdown portions of the gravitational waveform for black holes present in the new spacetime. With these tools in hand, one can then test gravitational wave signals for such effects by two different methods, (i) inspiral-merger-ringdown consistency test, and (ii) parameterized test. In this paper, we demonstrate the exact recipe one needs to do just this. We first derive parameterized corrections to the waveform inspiral, ringdown, and remnant properties for a generic non-Kerr spacetime and apply this to two example beyond-Kerr spacetimes each parameterized by a single non-Kerr parameter. We then predict the beyond-Kerr parameter magnitudes required in an observed gravitational wave signal to be statistically inconsistent with the Kerr case in general relativity. We find that the two methods give very similar bounds. The constraints found with existing gravitational-wave events are comparable to those from x-ray observations, while future gravitational-wave observations using Cosmic Explorer (Laser Interferometer Space Antenna) can improve such bounds by two (three) orders of magnitude.

We study the effect of compact extra dimensions on the gravitational wave luminosity and waveform. We consider a toy model, with a compactified fifth dimension, and matter confined on a brane. We work in the context of five dimensional ($5d$) general relativity, though we do make connections with the corresponding Kaluza-Klein effective $4d$ theory. We show that the luminosity of gravitational waves emitted in $5d$ gravity by a binary with the same characteristics (same masses and separation distance) as a $4d$ binary is 20.8\% less relative to the $4d$ case, to leading post-Newtonian order. The phase of the gravitational waveform differs by 26\% relative to the $4d$ case, to leading post-Newtonian order. Such a correction arises mainly due to the coupling between matter and dilaton field in the effective $4d$ picture and agrees with previous calculations when we set black holes' scalar charges to be those computed from the Kaluza-Klein reduction. The above correction is inconsistent with the recent gravitational-wave observations and it thus effectively rules out the possibility of such a simple compactified higher dimensions scenario. We also comment on how our results change if there are several compactified extra dimensions, and show that the discrepancy with $4d$ general relativity only increases.

Binary pulsars allow us to carry out precision tests of gravity and have placed stringent bounds on a broad class of theories beyond general relativity. Current and future radio telescopes, such as FAST, SKA, and MeerKAT, may find a new astrophysical system, a pulsar orbiting around a black hole, which will provide us a new source for probing gravity. In this paper, we systematically study the prospects of testing general relativity with such black hole-pulsar binaries. We begin by finding a mapping between generic non-Einsteinian parameters in the orbital decay rate and theoretical constants in various modified theories of gravity and then summarize this mapping with a ready-to-use list. Theories we study here include scalar-tensor theories, varying $G$ theories, massive gravity theories, generic screening gravity and quadratic curvature-corrected theories. We next use simulated measurement accuracy of the orbital decay rate with FAST/SKA to derive projected upper bounds on the generic non-Einsteinian parameters. We find that such bounds from black hole-pulsars can be stronger than those from neutron star-pulsar and neutron star-white dwarf binaries by a few orders of magnitude when the correction enters at negative post-Newtonian orders. By mapping such bounds on generic parameters to those on various modified theories of gravity, we find that one can constrain the time variation in Newton's constant $G$ to be comparable and complementary to the current strongest bound from solar system experiments. We also study how well one can probe quadratic gravity from black hole quadrupole moment measurements of black hole-pulsars. We find that bounds on the parity-violating sector of quadratic gravity can be stronger than current bounds by six orders of magnitude. These results suggest that a new discovery of black hole-pulsars in the future will provide powerful ways to probe gravity further.

The NICER Collaboration recently reported the measurement of the mass and radius of a pulsar PSR J0030+0451. We here use this new measurement to constrain one of the higher-order nuclear matter parameters $K_{\mathrm{sym,0}}$. We further combine the tidal measurement of the binary neutron star merger GW170817 by LIGO/Virgo to derive a joint 1-$\sigma$ constraint as $K_{\mathrm{sym,0}} = -102^{+71}_{-72}$ MeV. We believe this is the most reliable bound on the parameter to date under the assumption that there is no new physics above the saturation density which impacts neutron star observations.

The tidal measurement of gravitational waves from the binary neutron star merger event GW170817 allows us to probe nuclear physics that suffers less from astrophysical systematics compared to neutron star radius measurements with electromagnetic wave observations. A recent work found strong correlation among neutron-star tidal deformabilities and certain combinations of nuclear parameters associated with the equation of state. These relations were then used to derive bounds on such parameters from GW170817 assuming that the relations and neutron star masses are known exactly. Here, we expand on this important work by taking into account a few new considerations: (1) a broader class of equations of state; (2) correlations with the mass-weighted tidal deformability that was directly measured with GW170817; (3) how the relations depend on the binary mass ratio; (4) the uncertainty from equation of state variation in the correlation relations; (5) adopting the updated tidal deformability measurement from GW170817. Upon these new considerations, we find GW170817 bounds on nuclear parameters (the incompressibility $K_0$, its slope $M_0$ and the curvature of symmetry energy $K_{\mathrm{sym},0}$ at nuclear saturation density) to be 81 MeV $\leq K_0 \leq$ 362 MeV, 1556 MeV $\leq M_0 \leq$ 4971 MeV, and -254 MeV $\leq K_\mathrm{sym,0} \leq$ 27 MeV, which are more conservative than previously found with systematic errors more properly taken into account.

The extreme-gravity collisions of binaries with one black hole and one neutron star provide for excellent tests of general relativity. We here study how well one can constrain theories beyond general relativity with additional scalar fields that allow for spontaneous scalarization of neutron stars, and those motivated from string theory. We find that existing bounds can be improved with current gravitational-wave detectors if the black hole mass is sufficiently small. Bounds will further improve by many orders of magnitude with future detections, especially by combining multiple events.

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.

The detection of gravitational waves has offered us the opportunity to explore the dynamical and strong-field regime of gravity. Because matched filtering is more sensitive to variations in the gravitational waveform phase than the amplitude, many tests of gravity with gravitational waves have been carried out using only the former. Such studies cannot probe the non-Einsteinian effects that may enter only in the amplitude. Besides, if not accommodated in the waveform template, a non-Einsteinian effect in the amplitude may induce systematic errors on other parameters such as the luminosity distance. In this paper, we derive constraints on a few modified theories of gravity (Einstein-dilaton-Gauss-Bonnet gravity, scalar-tensor theories, and varying-$G$ theories), incorporating both phase and amplitude corrections. We follow the model-independent approach of the parametrized post-Einsteinian formalism. We perform Fisher analyses with Monte-Carlo simulations using the LIGO/Virgo posterior samples. We find that the contributions from amplitude corrections can be comparable to the ones from the phase corrections in case of massive binaries like GW150914. Also, constraints derived by incorporating both phase and amplitude corrections differ from the ones with phase corrections only by 4% at most, which supports many of the previous studies that only considered corrections in the phase. We further derive reliable constraints on the time-evolution of a scalar field in a scalar-tensor theory for the first time with gravitational waves.

Despite the tremendous success of general relativity so far, modified theories of gravity have received increased attention lately, motivated from both theoretical and observational aspects. Gravitational wave observations opened new possibilities for testing the viability of such theories in the strong-field regime. One could test each theory against observed data one at a time, though a more efficient approach would be to first probe gravity in a theory-agnostic way and map such information to that on specific theories afterward. One example of such model-independent tests with gravitational waves is the parameterized post-Einsteinian formalism, where one introduces generic parameters in the amplitude and phase that capture non-Einsteinian effects. In this paper, we derive gravitational waveforms from inspiraling compact binaries in various modified theories of gravity that violate at least one fundamental pillar in general relativity. We achieve this by deriving relations between corrections to the waveform amplitude/phase and those to the frequency evolution and Kepler's third law, since the latter two have already been (or can easily be) derived in many theories. Such an analysis allows us to derive corrections to the waveform amplitude, which extends many of previous works that focused on phase corrections only. Moreover, we derive gravitational waveforms in theories with varying gravitational constant. We extend the previous work by introducing two gravitational constants (the conservative one entering in the binding energy and the dissipative one entering in the gravitational wave luminosity) and allowing masses of binary constituents to vary with time. We also correct some errors in previous literature. Our results can be used to improve current analyses of testing general relativity as well as to achieve new projected constraints on many modified theories of gravity.

The space-borne gravitational wave interferometer, LISA, is expected to detect signals from numerous binary white dwarfs. At small orbital separation, rapid rotation and large tidal bulges may allow for the stellar internal structure to be probed through such observations. Finite-size effects are encoded in quantities like the moment of inertia ($I$), tidal Love number (Love), and quadrupole moment ($Q$). The universal relations among them (I-Love-Q relations) can be used to reduce the number of parameters in the gravitational-wave templates. We here study I-Love-Q relations for more realistic white dwarf models than used in previous studies. In particular, we extend previous works by including (i) differential rotation and (ii) internal temperature profiles taken from detailed stellar evolution calculations. We use the publicly available stellar evolution code MESA to generate cooling models of both low- and high-mass white dwarfs. We show that differential rotation causes the I-Q relation (and similarly the Love-Q relation) to deviate from that of constant rotation. We also find that the introduction of finite temperatures causes the white dwarf to move along the zero-temperature mass sequence of I-Q values, moving towards values that suggest a lower mass. We further find that after only a few Myrs, high-mass white dwarfs are well-described by the zero-temperature model, suggesting that the relations with zero-temperature may be good enough in most practical cases. Low-mass, He-core white dwarfs with thick hydrogen envelopes may undergo long periods of H burning which sustain the stellar temperature and allow deviations from the I-Love-Q relations for longer times.

The violent collisions of black holes provide for excellent test-beds of Einstein's general relativity in the strong/dynamical gravity regime. We here demonstrate the resolving power one can gain upon the use of multi-band observations of gravitational waves from both ground- and space-based detectors. We find significant improvement in both generic parameterized tests of general relativity and consistency tests of inspiral-merger-ringdown parts of the waveform over single-band detections. Such multi-band observations are crucial for unprecedented probes of e.g. parity-violation in gravity.

Gravitational wave observations offer unique opportunities to probe gravity in the strong and dynamical regime, which was difficult to access previously. We here review two theory-agnostic ways to carry out tests of general relativity with gravitational waves, namely (i) parameterized waveform tests and (ii) consistency tests between the inspiral and merger-ringdown portions. For each method, we explain the formalism, followed by results from existing events, and finally we discuss future prospects with upgraded detectors, including the possibility of using multi-band gravitational-wave observations with ground-based and space-borne interferometers. We show that such future observations have the potential to improve upon current bounds on theories beyond general relativity by many orders of magnitude. We conclude by listing several open questions that remain to be addressed.

The recent detection of gravitational waves from merging neutron star events has opened a new window on the many unknown aspects of their internal dynamics. A key role in this context is played by the transition from baryon to quark matter described in the neutron star equation of state (EoS). In particular, the binary pulsar observation of heavy neutron stars requires appropriately stiff dense matter in order to counter gravitational collapse, at variance with the predictions of many phenomenological quark models. On the other side, the LIGO observations favor a softer EoS therefore providing a lower bound to the equation stiffness. We introduce a quantum chromodynamics (QCD) description of the neutron star's high baryon density regime where the pressure and energy density distributions are directly obtained from the matrix elements of the QCD energy momentum tensor. Recent ab initio calculations allow us to evaluate the energy-momentum tensor in a model independent way including both quark and gluon degrees of freedom. Our approach is a first effort to replace quark models and effective gluon interactions with a first principles, fully QCD-based description. Most importantly, the QCD energy momentum tensor matrix elements are connected to the Mellin moments of the generalized parton distributions which can be measured in deeply virtual exclusive scattering experiments. As a consequence, we establish a connection between observables from high energy experiments and from the analysis of gravitational wave events. Both can be used to mutually constrain the respective sets of data. In particular, the emerging QCD-based picture is consistent with the GW170817 neutron star merger event once we allow a first-order phase transition from a low-density nuclear matter EoS to the newly-constructed high-density quark-gluon one.

The gravitational wave emission from the merging binary neutron star system GW170817 arrived full of tidal information which can be used to probe the fundamental ultra-dense nuclear physics residing in these stars. In previous work, we used two-dimensional correlations between nuclear matter parameters and tidal deformabilities of neutron stars applying specifically to GW170817 to derive constraints on the former. Here, we extend this analysis by finding similar correlations for varying chirp masses, the dominant determining factor in the frequency evolution of the inspiral, such that one can apply the same method to future detections. We estimate how accurately one can measure nuclear parameters with future gravitational wave interferometers and show how such measurements can be improved by combining multiple events. We find that bounds on the nuclear parameters with future observations can improve from the current one with GW170817 only by $\sim 30\%$ due to the existence of systematic errors caused mainly by the remaining uncertainty in the equation of state near and just above the nuclear saturation density. We show that such systematic errors can be reduced by considering multidimensional correlations among nuclear parameters and tidal deformabilities with various neutron star masses.

Studies of neutron stars are extremely timely given the recent detection of gravitational waves from a binary neutron star merger GW170817, and an International Space Station payload NICER currently in operation that aims to determine radii of neutron stars to a precision better than 5%. In many cases, neutron star solutions are constructed numerically due to the complexity of the field equations with realistic equations of state. However, in order to relate observables like the neutron star mass and radius to interior quantities like central density and pressure, it would be useful to provide an accurate, analytic modeling of a neutron star interior. One such solution for static and isolated neutron stars is the Tolman VII solution characterized only by two parameters (e.g. mass and radius), though its agreement with numerical solutions is not perfect. We here introduce an improved analytic model based on the Tolman VII solution by introducing an additional parameter to make the analytic density profile agree better with the numerically obtained one. This additional parameter can be fitted in terms of the stellar mass, radius and central density in an equation-of-state-insensitive way. In most cases, we find that the new model more accurately describes realistic profiles than the original Tolman VII solution by a factor of 2-5. Our results are first-step calculations towards constructing analytic interior solutions for more realistic neutron stars under rotation or tidal deformation.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

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.

Gravitational-wave observations in the near future may allow us to measure tidal deformabilities of neutron stars, which leads us to the understanding of physics at nuclear density. In principle, the gravitational waveform depends on various tidal parameters, which correlate strongly. Therefore, it would be useful if one can express such tidal parameters with a single parameter. Here, we report on universal relations among various $\ell$-th (dimensionless) electric, magnetic and shape tidal deformabilities in neutron stars and quark stars that do not depend sensitively on the equation of state. Such relations allow us to break the degeneracy among the tidal parameters. In this paper, we focus on gravitational waves from non-spinning neutron-star binary inspirals. We first derive the leading contribution of the $\ell$-th electric and $\ell=2$ magnetic tidal deformabilities to the gravitational-wave phase, which enters at $2\ell +1$ and $6$ post-Newtonian orders relative to the leading Newtonian one respectively. We then calculate the useful number of gravitational-wave cycles and show that not only the $\ell=2$ but also $\ell=3$ electric tidal deformabilities are important for parameter estimation with third-generation gravitational-wave detectors such as LIGO III and Einstein Telescope. Although the correlation between the $\ell=2$ and $\ell=3$ electric tidal deformabilities deteriorate the measurement accuracy of the former deformability parameter, one can increase its measurement accuracy significantly by using the universal relation. We provide a fitting formula for the LIGO III noise curve in the Appendix.

Gravitational wave observations of GW170817 placed bounds on the tidal deformabilities of compact stars allowing one to probe equations of state for matter at supranuclear densities. Here we design new parametrizations for hybrid hadron-quark equations of state, that give rise to low-mass twin stars, and test them against GW170817. We find that GW170817 is consistent with the coalescence of a binary hybrid star--neutron star. We also test and find that the I-Love-Q relations for hybrid stars in the third family agree with those for purely hadronic and quark stars within $\sim 3\%$ for both slowly and rapidly rotating configurations, implying that these relations can be used to perform equation-of-state independent tests of general relativity and to break degeneracies in gravitational waveforms for hybrid stars in the third family as well.

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.

The detection of gravitational waves (GWs) from black hole (BH) mergers provides an inroad toward probing the interior of astrophysical BHs. The general-relativistic description of the BH interior is that of empty spacetime with a (possibly) singular core. Recently, however, the hypothesis that the BH interior does not exist has been gaining traction, as it provides a means for resolving the BH information-loss problem. Here, we propose a simple method for answering the following question: Does the BH interior exist and, if so, does it contain some distribution of matter or is it mostly empty? Our proposal is premised on the idea that, similar to the case of relativistic, ultra-compact stars, any BH-like object whose interior has some matter distribution should support fluid modes in addition to the conventional spacetime modes. In particular, the Coriolis-induced Rossby (r-) modes, whose spectrum is mostly insensitive to the composition of the interior matter, should be a universal feature of such BH-like objects. In fact, the frequency and damping time of these modes are determined by only the object's mass and speed of rotation. The r-modes oscillate at a lower frequency, decay at a slower rate and produce weaker GWs than do modes of the spacetime class. Hence, they imprint a model-insensitive signature of a non-empty interior in the GW spectrum resulting from a BH merger. We find that future GW detectors, such as Advanced LIGO with its design sensitivity, have the potential of detecting such r-modes if the amount of GWs leaking out quantum mechanically from the interior of a BH-like object is sufficiently large.

The recent discovery of gravitational wave events has offered us unique testbeds of gravity in the strong and dynamical field regime. One possible modification to General Relativity is the gravitational parity violation that gives rise to the amplitude birefringence in gravitational waves, where one of the circularly-polarized mode is amplified while the other one is suppressed during their propagation. In this paper, we study how well one can measure gravitational parity violation via the amplitude birefringence of gravitational waves from stellar-mass black hole binaries. We choose Chern-Simons gravity as an example and work within an effective field theory formalism. We consider gravitational waves from both individual sources and stochastic gravitational wave backgrounds. Regarding bounds from individual sources, we estimate them using a Fisher analysis and carry out Monte Carlo simulations by randomly distributing sources over their sky location and binary orientation. We find that the bounds on the scalar field evolution in Chern-Simons gravity from GW150914 are too weak to satisfy the weak Chern-Simons approximation, while aLIGO with its design sensitivity can place meaningful bounds. Regarding bounds from stochastic gravitational wave backgrounds, we set the threshold signal-to-noise ratio for detection of the parity-violation mode as 5 and estimate projected bounds with future detectors assuming that signals are consistent with no parity violation. We find that a network of two third-generation detectors is able to place bounds that are slightly stronger than current binary pulsar bounds. Since gravitational wave observations probe either the difference in parity violation between the source and the detector or the line-of-sight integration of the scalar field, such bounds are complementary to local measurements from solar system experiments and binary pulsar observations.

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.

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.

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.

The primary target for the planned space-borne gravitational wave interferometers DECIGO/BBO is a primordial gravitational wave background (PGWB). However there exist astrophysical foregrounds and among them, gravitational waves from neutron star (NS) binaries are the solid and strong component that must be identified and subtracted. In this paper, we discuss the geometry of detector configurations preferable for identifying the NS/NS binary signals. As a first step, we analytically estimate the minimum signal-to-noise ratios (SNRs) of the binaries for several static detector configurations that are characterized by adjustable geometrical parameters, and determine the optimal values for these parameters. Next we perform numerical simulations to take into account the effect of detector motions, and find reasonable agreements with the analytical results. We show that, with the standard network formed by 4 units of triangle detectors, the proposed BBO sensitivity would be sufficient in receiving gravitational waves from all the NS/NS binaries at $z\le 5$ with SNRs higher than 25. We also discuss the minimum sensitivity of DECIGO required for the foreground identification.

Two low mass neutron stars, J0737-3039B and the companion to J1756-2251, show strong evidence of being formed from the collapse of an ONeMg core in an electron capture supernova (ECSN) or in an ultra-stripped iron core collapse supernova (FeCCSN). Using three different systematically generated sets of equations of state we explore the relationship between the moment of inertia of J0737-3039A and the binding energy of the two low mass neutron stars. We find this relationship, a less strict variant of the recently discovered I-Love-Q relations, is nevertheless more robust than a previously explored correlation between the binding energy and the slope of the nuclear symmetry energy L. We find that, if either J0737-3039B or the J1756-2251 companion were formed in an ECSN, no more than 0.06 solar masses could have been lost from the progenitor core, more than four times the mass loss predicted by current supernova modeling. A measurement of the moment of inertia of J0737-3039A to within 10% accuracy from pulsar timing, possible within a decade, can discriminate between formation scenarios such as ECSN or ultra-stripped FeCCSN and, given current constraints on the predicted core mass loss, potentially rule them out. Using the I-Love-Q relations we find that an Advanced LIGO can potentially measure the neutron star tidal polarizability to equivalent accuracy in a neutron star-neutron star merger at a distance of 200 Mpc, thus obtaining similar constraints on the formation scenarios. Such information on the occurrence of ECSNe is important for population synthesis calculations, especially for estimating the rate of binary neutron star mergers and resulting electromagnetic and gravitational wave signals. Further progress needs to be made modeling the core collapse process that leads to low-mass neutron stars, particularly in making robust predictions for the mass loss from the progenitor core.

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.

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.

General relativity has passed all solar system experiments and neutron star based tests, such as binary pulsar observations, with flying colors. A more exotic arena for testing general relativity is in systems that contain one or more black holes. Black holes are the most compact objects in the universe, providing probes of the strongest-possible gravitational fields. We are motivated to study strong-field gravity since many theories give large deviations from general relativity only at large field strengths, while recovering the weak-field behavior. In this article, we review how one can probe general relativity and various alternative theories of gravity by using electromagnetic waves from a black hole with an accretion disk, and gravitational waves from black hole binaries. We first review model-independent ways of testing gravity with electromagnetic/gravitational waves from a black hole system. We then focus on selected examples of theories that extend general relativity in rather simple ways. Some important characteristics of general relativity include (but are not limited to) (i) only tensor gravitational degrees of freedom, (ii) the graviton is massless, (iii) no quadratic or higher curvatures in the action, and (iv) the theory is 4 dimensional. Altering a characteristic leads to a different extension of general relativity: (i) scalar-tensor theories, (ii) massive gravity theories, (iii) quadratic gravity, and (iv) theories with large extra dimensions. Within each theory, we describe black hole solutions, their properties, and current and projected constraints on each theory using black hole-based tests of gravity. We close this review by listing some of the open problems in model-independent tests and within each specific theory.

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.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

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.

Gravity theories beyond General Relativity typically predict dipolar gravitational emission by compact-star binaries. This emission is sourced by "sensitivity" parameters depending on the stellar compactness. We introduce a general formalism to calculate these parameters, and show that in shift-symmetric Horndeski theories stellar sensitivities and dipolar radiation vanish, provided that the binary's dynamics is perturbative (i.e. the post-Newtonian formalism is applicable) and cosmological-expansion effects can be neglected. This allows reproducing the binary-pulsar observed orbital decay.

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.

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.

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.

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.

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.

We parametrize a large class of corrections to general relativity which include a long-ranged gravitational scalar field as a dynamical degree of freedom in two ways: parametrizing the structure of the correction to the action, and parametrizing the scalar hair (multipole structure) that compact objects and black holes attain. The presence of this scalar hair violates the no-hair theorems present in general relativity, which leads to several important effects. The effects we consider are (i) the interaction between an isolated body and an external scalar field, (ii) the scalar multipole-multipole interaction between two bodies in a compact binary, (iii) the additional pericenter precession of a binary, (iv) the scalar radiation from a binary, and (v) the modification to the gravitational wave phase from a binary. We apply this framework to example theories including Einstein-dilaton-Gauss-Bonnet gravity and dynamical Chern-Simons gravity, and estimate the size of the effects. Finally, we estimate the bounds that can be placed on parameters of the theories from the precession of pulsar binaries and from gravitational waves.

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.

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.

DECIGO Pathfinder (DPF) has an ability to detect gravitational waves from galactic intermediate-mass black hole binaries. If the signal is detected, it would be possible to determine parameters of the binary components. Furthermore, by using future space-borne gravitational wave interferometers, it would be possible to test alternative theories of gravity in the strong field regime. In this review article, we first explain how the detectors like DPF and DECIGO/BBO work and discuss the expected event rates. Then, we review how the observed gravitational waveforms from precessing compact binaries with slightly eccentric orbits can be calculated both in general relativity and in alternative theories of gravity. For the latter, we focus on Brans-Dicke and massive gravity theories. After reviewing these theories, we show the results of the parameter estimation with DPF using the Fisher analysis. We also discuss a possible joint search of DPF and ground-based interferometers. Then, we show the results of testing alternative theories of gravity using future space-borne interferometers. DECIGO/BBO would be able to place 4--5 orders of magnitude stronger constraint on Brans-Dicke theory than the solar system experiment. This is still 1--2 orders of magnitude stronger than the future solar system mission such as ASTROD I. On the other hand, LISA should be able to put 4 orders of magnitude more stringent constraint on the mass of the graviton than the current solar system bound. DPF may be able to place comparable constraint on the massive gravity theories as the solar system bound. We also discuss the prospects of using eLISA and ASTROD-GW in testing alternative theories of gravity. The bounds using eLISA are similar to the LISA ones, but ASTROD-GW performs the best in constraining massive gravity theories among all the gravitational wave detectors considered in this article.

It was recently shown that a black hole (BH) is the only compact object that can acquire a scalar charge in scalar Gauss-Bonnet (sGB) theory under the small coupling approximation. This leads to the fact that scalar radiation is emitted from a binary containing at least one BH. In this letter, we find the constraints on this theory from BH low-mass X-ray binaries (BH-LMXBs). The main result of this letter is that from the orbital decay rate of A0620-00, we obtained a conservative bound that is six orders of magnitude stronger than the solar system bound. In addition to this, we look at XTE J1118+480, whose orbital decay rate has been recently measured with an excess compared to the theoretical prediction in GR due to the radiation reaction. The cause of this excess is currently unknown. Although it is likely that the cause is of astrophysical origin, here we investigate the possibility of explaining this excess with the additional scalar radiation in sGB theory. We find that there still remains a parameter range where the excess can be explained while also satisfying the constraint obtained from A0620-00. The interesting point is that for most of other alternative theories of gravity, it seems difficult to explain this excess with the additional radiation. This is because it would be difficult to evade the constraints from binary pulsars or they have already been constrained rather strongly from other observations such as solar system experiments. We propose several ways to determine whether the excess is caused by the scalar radiation in sGB gravity including future gravitational wave observations with space-borne interferometers, which can give a constraint three orders of magnitude stronger than that from A0620-00.

Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect $\sim10^6$ neutron-star (NS) binaries and determine the luminosity distances to the binaries with high precision. Combining the luminosity distances with cosmologically-induced phase corrections on the GWs, cosmological expansion out to high redshift can be measured without the redshift determinations of host galaxies by electromagnetic observation and can be a unique probe for dark energy. This article is based on the results obtained in [1] where we investigated constraining power of the GW standard siren without redshift information on the equation of state of dark energy with future space-based GW detectors. We also compare the results with those obtained with other instruments and methods.

If we assume that we live in the center of a spherical inhomogeneous universe, we can explain the apparent accelerating expansion of the universe without introducing the unknown dark energy or modifying gravitational theory. Direct measurement of the cosmic acceleration can be a powerful tool in distinguishing $\Lambda$CDM and the inhomogeneous models. If $\Lambda$CDM is the correct model, we have shown that DECIGO/BBO has sufficient ability to detect the positive redshift drift of the source by observing gravitational waves from neutron star binaries for 5-10 years. This enables us to rule out any Lema\^itre-Tolman-Bondi (LTB) void model with monotonically increasing density profile. Furthermore, by detecting the positive redshift drift at $z\sim 0$, we can even rule out generic LTB models unless we allow unrealistically steep density gradient at $z\sim 0$. We also show that the measurement accuracy is slightly improved when we consider the joint search of DECIGO/BBO and the third generation Einstein Telescope. This test can be performed with GW observations alone without any reference to electromagnetic observations.

DECIGO Path Finder (DPF) is a space-borne gravitational wave (GW) detector with sensitivity in the frequency band 0.1--100Hz. As a first step mission to DECIGO, it is aiming for launching in 2016--2017. Although its main objective is to demonstrate technology for GW observation in space, DPF still has a chance of detecting GW signals and performing astrophysical observations. With an observable range up to 50 kpc, its main targets are GW signals from galactic intermediate mass black hole (IMBH) binaries. By using inspiral-merger-ringdown phenomenological waveforms, we perform both pattern-averaged analysis and Monte Carlo simulations including the effect of detector motion to find that the masses and (effective) spins of the IMBHs could be determined with errors of a few percent, should the signals be detected. Since GW signals from IMBH binaries with masses above $10^4 M_\odot$ cannot be detected by ground-based detectors, these objects can be unique sources for DPF. If the inspiral signal of a $10^3M_\odot$ IMBH binary is detected with DPF, it can give alert to the ringdown signal for the ground-based detectors $10^2$--$10^3$s before coalescence. We also estimate the possible bound on the graviton Compton wavelength from a possible IMBH binary in $\omega$ Centauri. We obtain a slightly weaker constraint than the solar system experiment and an about 2 orders of magnitude stronger constraint than the one from binary pulsar tests. Unfortunately, the detection rate of IMBH binaries is rather small.

Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect ~10^6 neutron-star (NS) binaries and determine the luminosity distances to the binaries with high precision. Combining the luminosity distances with cosmologically-induced phase corrections on the GWs, cosmological expansion out to high redshift can be measured without the redshift determinations of host galaxies by electromagnetic observation and be a unique probe for dark energy. On the other hand, such a NS-binary foreground should be subtracted to detect primordial GWs produced during inflation. Thus, the constraining power on dark energy and the detectability of the primordial gravitational waves strongly depend on the detector sensitivity and are in close relation with one another. In this paper, we investigate the constraints on the equation of state of dark energy with future space-based GW detectors with/without identifying the redshifts of host galaxies. We also study the sensitivity to the primordial GWs, properly dealing with the residual of the NS binary foreground. Based on the results, we discuss the detector sensitivity required to achieve the forementioned targeted study of cosmology.

When binary black holes are embedded in a gaseous environment, a rotating disk surrounding them, the so-called circumbinary disk, will be formed. The binary exerts a gravitational torque on the circumbinary disk and thereby the orbital angular momentum is transferred to it, while the angular momentum of the circumbinary disk is transferred to the binary through the mass accretion. The binary undergoes an orbital decay due to both the gravitational wave emission and the binary-disk interaction. This causes the phase evolution of the gravitational wave signal. The precise measurement of the gravitational wave phase thus may provide information regarding the circumbinary disk. In this paper, we assess the detectability of the signature of the binary-disk interaction using the future space-borne gravitational wave detectors such as DECIGO and BBO by the standard matched filtering analysis. We find that the effect of the circumbinary disk around binary black holes in the mass range $6M_sun\le{M}\lesssim3\times10^3M_sun$ is detectable at a statistically significant level in five year observation, provided that gas accretes onto the binary at a rate greater than $\dot{M}\sim1.4\times10^{17} [gs^{-1}] j^{-1}(M/10M_sun)^{33/23}$ with 10% mass-to-energy conversion efficiency, where j represents the efficiency of the angular momentum transfer from the binary to the circumbinary disk. We show that $O(0.1)$ coalescence events are expected to occur in sufficiently dense molecular clouds in five year observation. We also point out that the circumbinary disk is detectable, even if its mass at around the inner edge is by over 10 orders of magnitude less than the binary mass.

In Randall-Sundrum II (RS-II) braneworld model, it has been conjectured according to the AdS/CFT correspondence that brane-localized black hole (BH) larger than the bulk AdS curvature scale $\ell$ cannot be static, and it is dual to a four dimensional BH emitting the Hawking radiation through some quantum fields. In this scenario, the number of the quantum field species is so large that this radiation changes the orbital evolution of a BH binary. We derived the correction to the gravitational waveform phase due to this effect and estimated the upper bounds on $\ell$ by performing Fisher analyses. We found that DECIGO/BBO can put a stronger constraint than the current table-top result by detecting gravitational waves from small mass BH/BH and BH/neutron star (NS) binaries. Furthermore, DECIGO/BBO is expected to detect 10$^5$ BH/NS binaries per year. Taking this advantage, we found that DECIGO/BBO can actually measure $\ell$ down to $\ell=0.33 \mu$m for 5 year observation if we know that binaries are circular a priori. This is about 40 times smaller than the upper bound obtained from the table-top experiment. On the other hand, when we take eccentricities into binary parameters, the detection limit weakens to $\ell=1.5 \mu$m due to strong degeneracies between $\ell$ and eccentricities. We also derived the upper bound on $\ell$ from the expected detection number of extreme mass ratio inspirals (EMRIs) with LISA and BH/NS binaries with DECIGO/BBO, extending the discussion made recently by McWilliams. We found that these less robust constraints are weaker than the ones from phase differences.

One possibility for explaining the apparent accelerating expansion of the universe is that we live in the center of a spherically inhomogeneous universe. Although current observations cannot fully distinguish $\Lambda$CDM and these inhomogeneous models, direct measurement of the acceleration of the universe can be a powerful tool in probing them. We have shown that, if $\Lambda$CDM is the correct model, DECIGO/BBO would be able to detect the positive redshift drift (which is the time evolution of the source redshift $z$) in 3--5 year gravitational wave (GW) observations from neutron-star binaries, which enables us to rule out any Lema\^itre-Tolman-Bondi (LTB) void model with monotonically increasing density profile. We may even be able to rule out any LTB model unless we allow unrealistically steep density profile at $z\sim 0$. This test can be performed with GW observations alone, without any reference to electromagnetic observations, and is more powerful than the redshift drift measurement using Lyman $\alpha$ forest.

We calculate how strongly one can constrain the alternative theories of gravity with deci-Hz gravitational wave interferometers such as DECIGO and BBO. Here we discuss Brans-Dicke theory and massive graviton theories as typical examples. We consider the inspiral of compact binaries composed of a neutron star (NS) and an intermediate mass black hole (IMBH) for Brans-Dicke (BD) theory and those composed of a super massive black hole (SMBH) and a black hole (SMBH) for massive graviton theories. Using the restricted 2PN waveforms including spin effects and taking the spin precession into account, we perform the Monte Carlo simulations of $10^4$ binaries to estimate the determination accuracy of binary parameters including the Brans-Dicke parameter $\omega_{\mathrm{BD}}$ and the graviton Compton length $\lambda_g$. Assuming a $(1.4, 10)M_{\odot}$ NS/BH binary of SNR=$\sqrt{200}$, the constraint on $\omega_{\mathrm{BD}}$ is obtained as $\omega_{\mathrm{BD}}>2.32\times 10^6$, which is 300 times stronger than the estimated constraint from LISA observation. Furthermore, we find that, due to the expected large merger rate of NS/BH binaries of $O(10^4)$ yr$^{-1}$, a statistical analysis yields $\omega_{\mathrm{BD}}>3.77\times10^8$, which is 4 orders of magnitude stronger than the current strongest bound obtained from the solar system experiment. For massive graviton theories, assuming a $(10^6, 10^5)M_{\odot}$ BH/BH binary at 3Gpc, one can put a constraint $\lambda_g>3.35\times10^{20}$cm, on average. This is three orders of magnitude stronger than the one obtained from the solar system experiment. From these results, it is understood that DECIGO/BBO is a very powerful tool for constraining alternative theories of gravity, too.

We calculate how strongly one can put constraints on alternative theories of gravity such as Brans-Dicke and massive graviton theories with LISA. We consider inspiral gravitational waves from a compact binary composed of a neutron star (NS) and an intermediate mass black hole (IMBH) in Brans-Dicke (BD) theory and that composed of 2 super massive black holes (SMBHs) in massive graviton theories. We use the restricted 2PN waveforms including the effects of spins. We also take both precession and eccentricity of the orbit into account. For simplicity, we set the fiducial value for the spin of one of the binary constituents to zero so that we can apply the approximation called \textit{simple precession}. We perform the Monte Carlo simulations of $10^4$ binaries, estimating the determination accuracy of binary parameters including the BD parameter $\omega_{\mathrm{BD}}$ and the Compton wavelength of graviton $\lambda_g$ for each binary using the Fisher matrix method. We find that including both the spin-spin coupling $\sigma$ and the eccentricity $e$ into the binary parameters reduces the determination accuracy by an order of magnitude for the Brans-Dicke case, whilst it has less influence on massive graviton theories. On the other hand, including precession enhances the constraint on $\omega_{\mathrm{BD}}$ only 20$%$ but it increases the constraint on $\lambda_g$ by an order of magnitude. Using a $(1.4+1000)M_{\odot}$ NS/BH binary of SNR=$\sqrt{200}$, one can put a constraint $\omega_{\mathrm{BD}}>6944$, whilst using a $(10^7+10^6)M_{\odot}$ BH/BH binary at 3Gpc, one can put $\lambda_g>3.06\times10^{21}$cm, on average. The latter is 4 orders of magnitude stronger than the one obtained from the solar system experiment. These results indicate that the effects of precession and eccentricity must be taken carefully in the parameter estimation analysis.