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

The 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.

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.

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.

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.

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 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.