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 theories with multiple scalar fields coupled to the metric and each other --- a natural extension of the well studied single-scalar-tensor theories --- are interesting phenomenological frameworks to describe deviations from general relativity in the strong-field regime. In these theories, the $N$-tuple of scalar fields takes values in a coordinate patch of an $N$-dimensional Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric target-space manifold. Within this model we present a preliminary investigation of spontaneous scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, and that the properties of strongly scalarized stellar configurations additionally depend on the target-space curvature radius. In preparation for numerical relativity simulations, we also write down the $3+1$ decomposition of the field equations for generic tensor-multi-scalar theories.

We study torsional oscillations of neutron stars in the scalar-tensor theory of gravity using the relativistic Cowling approximation. We compute unperturbed neutron star models adopting realistic equations of state for the neutron star's core and crust. For scalar-tensor theories that allow for spontaneous scalarization, the crust thickness can be significantly smaller than in general relativity. We derive the perturbation equation describing torsional oscillations in scalar-tensor theory, and we solve the corresponding eigenvalue problem to find the oscillation frequencies. The fundamental mode (overtone) frequencies become smaller (larger) than in general relativity for scalarized stellar models. Torsional oscillation frequencies may yield information on the crust microphysics if microphysics effects are not degenerate with strong-gravity effects, such as those due to scalarization. To address this issue, we consider two different models for the equation of state of the crust and we look at the effects of electron screening. The effect of scalarization on torsional oscillation frequencies turns out to be smaller than uncertainties in the microphysics for all spontaneous scalarization models allowed by binary pulsar observations. Our study shows that the observation of quasi-periodic oscillations (QPOs) following giant flares can be used to constrain the microphysics of neutron star crusts, whether spontaneous scalarization occurs or not.

Scalar-tensor theories are a compelling alternative to general relativity and one of the most accepted extensions of Einstein's theory. Black holes in these theories have no hair, but could grow "wigs" supported by time-dependent boundary conditions or spatial gradients. Time-dependent or spatially varying fields lead in general to nontrivial black hole dynamics, with potentially interesting experimental consequences. We carry out a numerical investigation of the dynamics of single and binary black holes in the presence of scalar fields. In particular we study gravitational and scalar radiation from black-hole binaries in a constant scalar-field gradient, and we compare our numerical findings to analytical models. In the single black hole case we find that, after a short transient, the scalar field relaxes to static configurations, in agreement with perturbative calculations. Furthermore we predict analytically (and verify numerically) that accelerated black holes in a scalar-field gradient emit scalar radiation. For a quasicircular black-hole binary, our analytical and numerical calculations show that the dominant component of the scalar radiation is emitted at twice the binary's orbital frequency.

We examine the motion of light fields near the bottom of a potential valley in a multi-dimensional field space. In the case of two fields we identify three general scales, all of which must be large in order to justify an effective low-energy approximation involving only the light field, $\ell$. (Typically only one of these -- the mass of the heavy field transverse to the trough -- is used in the literature when justifying the truncation of heavy fields.) We explicitly compute the resulting effective field theory, which has the form of a $P(\ell,X)$ model, with $X = - 1/2(\partial \ell)^2$, as a function of these scales. This gives the leading ways each scale contributes to any low-energy dynamics, including (but not restricted to) those relevant for cosmology. We check our results with the special case of a homogeneous roll near the valley floor, placing into a broader context recent cosmological calculations that show how the truncation approximation can fail. By casting our results covariantly in field space, we provide a geometrical criterion for model-builders to decide whether or not the single-field and/or the truncation approximation is justified, identify its leading deviations, and to efficiently extract cosmological predictions.

Scalar-tensor theories are among the simplest extensions of general relativity. In theories with light scalars, deviations from Einstein's theory of gravity are determined by the scalar mass m_s and by a Brans-Dicke-like coupling parameter \omega_{BD}. We show that gravitational-wave observations of nonspinning neutron star-black hole binary inspirals can be used to set lower bounds on \omega_{BD} and upper bounds on the combination m_s/\sqrt{\omega_{BD}}$. We estimate via a Fisher matrix analysis that individual observations with signal-to-noise ratio \rho would yield (m_s/\sqrt{\omega_{BD}})(\rho/10)<10^{-15}, 10^{-16} and 10^{-19} eV for Advanced LIGO, ET and eLISA, respectively. A statistical combination of multiple observations may further improve these bounds.

Observations of pulsar timing provide strong constraints on scalar-tensor theories of gravity, but these constraints are traditionally quoted as limits on the microscopic parameters (like the Brans-Dicke coupling, for example) that govern the strength of scalar-matter couplings at the particle level in particular models. Here we present fits to timing data for several pulsars directly in terms of the phenomenological couplings (masses, scalar charges, moment of inertia sensitivities and so on) of the stars involved, rather than to the more microscopic parameters of a specific model. For instance, for the double pulsar PSR J0737-3039A/B we find at the 68% confidence level that the masses are bounded by 1.28 < m_A/m_sun < 1.34 and 1.19 < m_B/m_sun < 1.25, while the scalar-charge to mass ratios satisfy |a_A| < 0.21, |a_B| < 0.21 and |a_B - a_A| < 0.002$. These constraints are independent of the details of the scalar tensor model involved, and of assumptions about the stellar equations of state. Our fits can be used to constrain a broad class of scalar tensor theories by computing the fit quantities as functions of the microscopic parameters in any particular model. For the Brans-Dicke and quasi-Brans-Dicke models, the constraints obtained in this manner are consistent with those quoted in the literature.

An old result ({\tt astro-ph/9905303}) by Jacobson implies that a black hole with Schwarzschild radius $r_s$ acquires scalar hair, $Q \propto r_s^2 \mu$, when the (canonically normalized) scalar field in question is slowly time-dependent far from the black hole, $\partial_t \phi \simeq \mu M_p$ with $\mu r_s \ll 1$ time-independent. Such a time dependence could arise in scalar-tensor theories either from cosmological evolution, or due to the slow motion of the black hole within an asymptotic spatial gradient in the scalar field. Most remarkably, the amount of scalar hair so induced is independent of the strength with which the scalar couples to matter. We argue that Jacobson's Miracle Hair-Growth Formula${}^\copyright$ implies, in particular, that an orbiting pair of black holes can radiate {\em dipole} radiation, provided only that the two black holes have different masses. Quasar OJ 287, situated at redshift $z \simeq 0.306$, has been argued to be a double black-hole binary system of this type, whose orbital decay recently has been indirectly measured and found to agree with the predictions of General Relativity to within 6%. We argue that the absence of observable scalar dipole radiation in this system yields the remarkable bound $|\,\mu| < (16 \, \hbox{days})^{-1}$ on the instantaneous time derivative at this redshift (as opposed to constraining an average field difference, $\Delta \phi$, over cosmological times), provided only that the scalar is light enough to be radiated --- i.e. $m \lsim 10^{-23}$ eV --- independent of how the scalar couples to matter. This can also be interpreted as constraining (in a more model-dependent way) the binary's motion relative to any spatial variation of the scalar field within its immediate vicinity within its host galaxy.

Precision tests of gravity can be used to constrain the properties of hypothetical very light scalar fields, but these tests depend crucially on how macroscopic astrophysical objects couple to the new scalar field. We develop quasi-analytic methods for solving the equations of stellar structure using scalar-tensor gravity, with the goal of seeing how stellar properties depend on assumptions made about the scalar coupling at a microscopic level. We illustrate these methods by applying them to Brans-Dicke scalars, and their generalization in which the scalar-matter coupling is a weak function of the scalar field. The four observable parameters that characterize the fields external to a spherically symmetric star (the stellar radius, R, mass, M, scalar `charge', Q, and the scalar's asymptotic value, phi_infty) are subject to two relations because of the matching to the interior solution, generalizing the usual mass-radius, M(R), relation of General Relativity. We identify how these relations depend on the microscopic scalar couplings, agreeing with earlier workers when comparisons are possible. Explicit analytical solutions are obtained for the instructive toy model of constant-density stars, whose properties we compare to more realistic equations of state for neutron star models.