The inspiral of a stellar mass ($1 - 100\,M_\odot$) compact body into a massive ($10^5 - 10^7\,M_\odot$) black hole has been a focus of much effort, both for the promise of such systems as astrophysical sources of gravitational waves, and because they are a clean limit of the general relativistic two-body problem. Our understanding of this problem has advanced significantly in recent years, with much progress in modeling the "self force" arising from the small body's interaction with its own spacetime deformation. Recent work has shown that this self interaction is especially interesting when the frequencies associated with the orbit's $\theta$ and $r$ motions are in an integer ratio: $\Omega_\theta/\Omega_r = \beta_\theta/\beta_r$, with $\beta_\theta$ and $\beta_r$ both integers. In this paper, we show that key aspects of the self interaction for such "resonant" orbits can be understood with a relatively simple Teukolsky-equation-based calculation of gravitational-wave fluxes. We show that fluxes from resonant orbits depend on the relative phase of radial and angular motions. The purpose of this paper is to illustrate in simple terms how this phase dependence arises using tools that are good for strong-field orbits, and to present a first study of how strongly the fluxes vary as a function of this phase and other orbital parameters. Future work will use the full dissipative self force to examine resonant and near resonant strong-field effects in greater depth, which will be needed to characterize how a binary evolves through orbital resonances.

We explore perturbations about a Friedmann-Robertson-Walker background in Chern-Simons gravity. At large momenta one of the two circularly polarized tensor modes becomes ghostlike. We argue that nevertheless the theory does not exhibit classical runaway solutions, except possibly in the relativistic nonlinear regime. However, the ghost modes cause the vacuum state to be quantum mechanically unstable, with a decay rate that is naively infinite. The decay rate can be made finite only if one interprets the theory as an effective quantum field theory valid up to some momentum cutoff, which violates Lorentz invariance. By demanding that the energy density in photons created by vacuum decay over the lifetime of the Universe not violate observational bounds, we derive strong constraints on the two dimensional parameter space of the theory, consisting of the cutoff and the Chern-Simons mass.

We explore a class of effective field theory models of cosmic acceleration involving a metric and a single scalar field. These models can be obtained by starting with a set of ultralight pseudo-Nambu-Goldstone bosons whose couplings to matter satisfy the weak equivalence principle, assuming that one boson is lighter than all the others, and integrating out the heavier fields. The result is a quintessence model with matter coupling, together with a series of correction terms in the action in a covariant derivative expansion, with specific scalings for the coefficients. After eliminating higher derivative terms and exploiting the field redefinition freedom, we show that the resulting theory contains nine independent free functions of the scalar field when truncated at four derivatives. This is in contrast to the four free functions found in similar theories of single-field inflation, where matter is not present. We discuss several different representations of the theory that can be obtained using the field redefinition freedom. For perturbations to the quintessence field today on subhorizon lengthscales larger than the Compton wavelength of the heavy fields, the theory is weakly coupled and natural in the sense of t'Hooft. The theory admits a regime where the perturbations become modestly nonlinear, but very strong nonlinearities lie outside its domain of validity.

We show that transient resonances occur in the two body problem in general relativity, in the highly relativistic, extreme mass-ratio regime for spinning black holes. These resonances occur when the ratio of polar and radial orbital frequencies, which is slowly evolving under the influence of gravitational radiation reaction, passes through a low order rational number. At such points, the adiabatic approximation to the orbital evolution breaks down, and there is a brief but order unity correction to the inspiral rate. Corrections to the gravitational wave signal's phase due to resonance effects scale as the square root of the inverse of mass of the small body, and thus become large in the extreme-mass-ratio limit, dominating over all other post-adiabatic effects. The resonances make orbits more sensitive to changes in initial data (though not quite chaotic), and are genuine non-perturbative effects that are not seen at any order in a standard post-Newtonian expansion. Our results apply to an important potential source of gravitational waves, the gradual inspiral of white dwarfs, neutron stars, or black holes into much more massive black holes. It is hoped to exploit observations of these sources to map the spacetime geometry of black holes. However, such mapping will require accurate models of binary dynamics, which is a computational challenge whose difficulty is significantly increased by resonance effects. We estimate that the resonance phase shifts will be of order a few tens of cycles for mass ratios $\sim 10^{-6}$, by numerically evolving fully relativistic orbital dynamics supplemented with an approximate, post-Newtonian self-force.

We present two methods for integrating forced geodesic equations in the Kerr spacetime, which can accommodate arbitrary forces. As a test case, we compute inspirals under a simple drag force, mimicking the presence of gas. We verify that both methods give the same results for this simple force. We find that drag generally causes eccentricity to increase throughout the inspiral. This is a relativistic effect qualitatively opposite to what is seen in gravitational-radiation-driven inspirals, and similar to what is observed in hydrodynamic simulations of gaseous binaries. We provide an analytic explanation by deriving the leading order relativistic correction to the Newtonian dynamics. If observed, an increasing eccentricity would provide clear evidence that the inspiral was occurring in a non-vacuum environment. Our two methods are especially useful for evolving orbits in the adiabatic regime. Both use the method of osculating orbits, in which each point on the orbit is characterized by the parameters of the geodesic with the same instantaneous position and velocity. Both methods describe the orbit in terms of the geodesic energy, axial angular momentum, Carter constant, azimuthal phase, and two angular variables that increase monotonically and are relativistic generalizations of the eccentric anomaly. The two methods differ in their treatment of the orbital phases and the representation of the force. In one method the geodesic phase and phase constant are evolved together as a single orbital phase parameter, and the force is expressed in terms of its components on the Kinnersley orthonormal tetrad. In the second method, the phase constants of the geodesic motion are evolved separately and the force is expressed in terms of its Boyer-Lindquist components. This second approach is a generalization of earlier work by Pound and Poisson for planar forces in a Schwarzschild background.

We argued in a previous paper [R. A. Vanderveld et al. 2006, arXiv:astro-ph/0602476] that negative deceleration parameters at the center of symmetry in Lemaitre-Tolman-Bondi cosmological models can only occur if the model is not smooth at the origin. Here we demonstrate explicitly the connection between non-smoothness and the failure of positivity theorems for deceleration. We also address some confusion that has arisen in the literature and respond to some recent criticisms of our arguments.

Recently there have been suggestions that the Type Ia supernova data can be explained using only general relativity and cold dark matter with no dark energy. In "Swiss cheese" models of the Universe, the standard Friedmann-Robertson-Walker picture is modified by the introduction of mass compensating spherical inhomogeneities, typically described by the Lemaitre-Tolman-Bondi metric. If these inhomogeneities correspond to underdense cores surrounded by mass-compensating overdense shells, then they can modify the luminosity distance-redshift relation in a way that can mimic accelerated expansion. It has been argued that this effect could be large enough to explain the supernova data without introducing dark energy or modified gravity. We show that the large apparent acceleration seen in some models can be explained in terms of standard weak field gravitational lensing together with insufficient randomization of void locations. The underdense regions focus the light less than the homogeneous background, thus dimming supernovae in a way that can mimic the effects of acceleration. With insufficient randomization of the spatial location of the voids and of the lines of sight, coherent defocusing can lead to anomalously large demagnification effects. We show that a proper randomization of the voids and lines of sight reduces the effect to the point that it can no longer explain the supernova data.

We consider theories in which there exists a nontrivial coupling between the dark matter sector and the sector responsible for the acceleration of the universe. Such theories can possess an adiabatic regime in which the quintessence field always sits at the minimum of its effective potential, which is set by the local dark matter density. We show that if the coupling strength is much larger than gravitational, then the adiabatic regime is always subject to an instability. The instability, which can also be thought of as a type of Jeans instability, is characterized by a negative sound speed squared of an effective coupled dark matter/dark energy fluid, and results in the exponential growth of small scale modes. We discuss the role of the instability in specific coupled CDM and Mass Varying Neutrino (MaVaN) models of dark energy, and clarify for these theories the regimes in which the instability can be evaded due to non-adiabaticity or weak coupling.

It has been suggested that the acceleration of the Universe may be due to the backreaction of perturbations to the Friedmann-Robertson-Walker background. For a Universe dominated by cold dark matter, it is known that the backreaction of superhorizon perturbations can not drive acceleration. We extend this result to models with cold dark matter together with a scalar field. We show that the scalar field can drive acceleration only via the standard mechanism of a constant or nearly constant piece of its potential (i.e., a cosmological constant); there is no separate mechanism involving superhorizon backreaction. This rules out some models which have been proposed in the literature.

We consider the cosmological constraints on theories in which there exists a nontrivial coupling between the dark matter sector and the sector responsible for the acceleration of the universe, in light of the most recent supernovae, large scale structure and cosmic microwave background data. For a variety of models, we show that the strength of the coupling of dark matter to a quintessence field is constrained to be less than 7% of the coupling to gravity. We also show that long range interactions between fermionic dark matter particles mediated by a light scalar with a Yukawa coupling are constrained to be less than 5% of the strength of gravity at a distance scale of 10 Mpc. We show that all of the models we consider are quantum mechanically weakly coupled, and argue that some other models in the literature are ruled out by quantum mechanical strong coupling.

Inspirals of stellar mass compact objects into massive black holes are an important source for future gravitational wave detectors such as Advanced LIGO and LISA. Detection of these sources and extracting information from the signal relies on accurate theoretical models of the binary dynamics. We cast the equations describing binary inspiral in the extreme mass ratio limit in terms of action angle variables, and derive properties of general solutions using a two-timescale expansion. This provides a rigorous derivation of the prescription for computing the leading order orbital motion. As shown by Mino, this leading order or adiabatic motion requires only knowledge of the orbit-averaged, dissipative piece of the self force. The two timescale method also gives a framework for calculating the post-adiabatic corrections. For circular and for equatorial orbits, the leading order corrections are suppressed by one power of the mass ratio, and give rise to phase errors of order unity over a complete inspiral through the relativistic regime. These post-1-adiabatic corrections are generated by the fluctuating piece of the dissipative, first order self force, by the conservative piece of the first order self force, and by the orbit-averaged, dissipative piece of the second order self force. We also sketch a two-timescale expansion of the Einstein equation, and deduce an analytic formula for the leading order, adiabatic gravitational waveforms generated by an inspiral.

Ground-based gravitational wave detectors may be able to constrain the nuclear equation of state using the early, low frequency portion of the signal of detected neutron star - neutron star inspirals. In this early adiabatic regime, the influence of a neutron star's internal structure on the phase of the waveform depends only on a single parameter lambda of the star related to its tidal Love number, namely the ratio of the induced quadrupole moment to the perturbing tidal gravitational field. We analyze the information obtainable from gravitational wave frequencies smaller than a cutoff frequency of 400 Hz, where corrections to the internal-structure signal are less than 10 percent. For an inspiral of two non-spinning 1.4 solar mass neutron stars at a distance of 50 Mpc, LIGO II detectors will be able to constrain lambda to lambda < 2.0 10^{37} g cm^2 s^2 with 90% confidence. Fully relativistic stellar models show that the corresponding constraint on radius R for 1.4 solar mass neutron stars would be R < 13.6 km (15.3 km) for a n=0.5 (n=1.0) polytrope.

We consider theories with a nontrivial coupling between the matter and dark energy sectors. We describe a small scale instability that can occur in such models when the coupling is strong compared to gravity, generalizing and correcting earlier treatments. The instability is characterized by a negative sound speed squared of an effective coupled dark matter/dark energy fluid. Our results are general, and applicable to a wide class of coupled models and provide a powerful, redshift-dependent tool, complementary to other constraints, with which to rule many of them out. A detailed analysis and applications to a range of models are presented in a longer companion paper.

We calculate the systematic inhomogeneity-induced correction to the cosmological constant that one would infer from an analysis of the luminosities and redshifts of Type Ia supernovae, assuming a homogeneous universe. The calculation entails a post-Newtonian expansion within the framework of second order perturbation theory, wherein we consider the effects of subhorizon density perturbations in a flat, dust dominated universe. Within this formalism, we calculate luminosity distances and redshifts along the past light cone of an observer. The resulting luminosity distance-redshift relation is fit to that of a homogeneous model in order to deduce the best-fit cosmological constant density Omega_Lambda. We find that the luminosity distance-redshift relation is indeed modified, by a small fraction of order 10^{-5}. When fitting this perturbed relation to that of a homogeneous universe, we find that the inferred cosmological constant can be surprisingly large, depending on the range of redshifts sampled. For a sample of supernovae extending from z=0.02 out to z=0.15, we find that Omega_Lambda=0.004. The value of Omega_Lambda has a large variance, and its magnitude tends to get larger for smaller redshifts, implying that precision measurements from nearby supernova data will require taking this effect into account. However, we find that this effect is likely too small to explain the observed value of Omega_Lambda=0.7. There have been previous claims of much larger backreaction effects. By contrast to those calculations, our work is directly related to how observers deduce cosmological parameters from astronomical data.

In most theories of gravity involving torsion, the source for torsion is the intrinsic spin of matter. Since the spins of fermions are normally randomly oriented in macroscopic bodies, the torsion generated is normally negligible. However, in a recent paper, Mao et al. point out that there is a class of theories in which the angular momentum of macroscopic spinning bodies generates a significant amount of torsion. They argue that by the principle of action equals reaction, one would expect the angular momentum of test bodies to couple to a background torsion field, and therefore the precession of the GPB gyroscopes should be affected in these theories by the torsion generated by the Earth. We show that in fact the principle of action equals reaction does not apply to these theories. We examine in detail a generalization of the Hayashi-Shirafuji theory suggested by Mao et al. called Einstein-Hayashi-Shirafuji theory. There are a variety of different versions of this theory, depending on the precise form of the coupling to matter chosen for the torsion. We show that for any coupling to matter that is compatible with the spin transport equation postulated by Mao et al., the theory has either ghosts or an ill-posed initial value formulation. These theoretical problems can be avoided by specializing the parameters of the theory and in addition choosing the standard minimal coupling to matter of the torsion tensor. This yields a consistent theory, but one in which the action equals reaction principle is violated, and in which the angular momentum of the gyroscopes does not couple to the Earth's torsion field. Thus, the Einstein-Hayashi-Shirafuji theory does not predict a detectable torsion signal for Gravity Probe B. There may be other torsion theories which do.

There has been much debate over whether or not one could explain the observed acceleration of the Universe with inhomogeneous cosmological models, such as the spherically-symmetric Lemaitre-Tolman-Bondi (LTB) models. It has been claimed that the central observer in these models can observe a local acceleration, which would contradict general theorems. We resolve the contradiction by noting that many of the models that have been explored contain a weak singularity at the location of the observer which makes them unphysical. In the absence of this singularity, we show that LTB models must have a positive central deceleration parameter $q_{0}$, in agreement with the general theorems. We also show that it is possible to achieve a negative apparent deceleration parameter at nonzero redshifts in LTB models that do not contain this singularity. However, we find other singularities that tend to arise in LTB models when attempting to match luminosity distance data, and these generally limit the range of redshifts for which these models can mimic observations of an accelerating Universe. Exceptional models do exist that can extend to arbitrarily large redshift without encountering these pathologies, and we show how these may be constructed. These special models exhibit regions with negative effective equation of state parameter, which may fall below negative one, but we have failed to find any singularity-free models that agree with observations. Moreover, models based on dust-filled LTB metrics probably fail to reproduce observed properties of large scale structure.

The fourth order Weyl gravity theory of Mannheim and Kazanas is based on replacing the Einstein-Hilbert action with the square of the Weyl tensor, and on modifying the matter action of the standard model of particle physics to make it conformally invariant. This theory has been suggested as a model of both dark matter and dark energy. We argue that the conformal invariance is not a fundamental property of the theory, and instead is an artifact of the choice of variables used in its description. We deduce that in the limit of weak fields and slow motions the theory does not agree with the predictions of general relativity, and is therefore ruled out by Solar System observations.

Einstein's special theory of relativity revolutionized physics by teaching us that space and time are not separate entities, but join as ``spacetime''. His general theory of relativity further taught us that spacetime is not just a stage on which dynamics takes place, but is a participant: The field equation of general relativity connects matter dynamics to the curvature of spacetime. Curvature is responsible for gravity, carrying us beyond the Newtonian conception of gravity that had been in place for the previous two and a half centuries. Much research in gravitation since then has explored and clarified the consequences of this revolution; the notion of dynamical spacetime is now firmly established in the toolkit of modern physics. Indeed, this notion is so well established that we may now contemplate using spacetime as a tool for other science. One aspect of dynamical spacetime -- its radiative character, ``gravitational radiation'' -- will inaugurate entirely new techniques for observing violent astrophysical processes. Over the next one hundred years, much of this subject's excitement will come from learning how to exploit spacetime as a tool for astronomy. This article is intended as a tutorial in the basics of gravitational radiation physics.

It has recently been suggested that the acceleration of the Universe can be explained as the backreaction effect of superhorizon perturbations using second order perturbation theory. If this mechanism is correct, it should also apply to a hypothetical, gedanken universe in which the subhorizon perturbations are absent. In such a gedanken universe it is possible to compute the deceleration parameter $q_0$ measured by comoving observers using local covariant Taylor expansions rather than using second order perturbation theory. The result indicates that second order corrections to $q_0$ are present, but shows that if $q_0$ is negative then its magnitude is constrained to be less than or of the order of the square of the peculiar velocity on Hubble scales today. We argue that since this quantity is constrained by observations to be small compared to unity, superhorizon perturbations cannot be responsible for the acceleration of the Universe.

A key source for LISA will be the inspiral of compact objects into supermassive black holes. Recently Mino has shown that in the adiabatic limit, gravitational waveforms for these sources can be computed using for the radiation reaction force the gradient of one half the difference between the retarded and advanced metric perturbations. Using post-Newtonian expansions, we argue that the resulting waveforms should be sufficiently accurate for signal detection with LISA. Data-analysis templates will require higher accuracy, going beyond adiabaticity; this remains a significant challenge. We describe an explicit computational procedure for obtaining waveforms based on Mino's result, for the case of a point particle coupled to a scalar field. We derive an expression for the time-averaged time derivative of the Carter constant, and verify that the expression correctly predicts that circular orbits remain circular while evolving via radiation reaction. The derivation uses detailed properties of mode expansions, Green's functions and bound geodesic orbits in the Kerr spacetime, which we review in detail. This paper is about three quarters review and one quarter new material. The intent is to give a complete and self-contained treatment of scalar radiation reaction in the Kerr spacetime, in a single unified notation, starting with the Kerr metric and ending with formulae for the time evolution of all three constants of the motion that are sufficiently explicit to be used immediately in a numerical code.

We describe progress evolving an important limit of binary orbits in general relativity, that of a stellar mass compact object gradually spiraling into a much larger, massive black hole. These systems are of great interest for gravitational wave observations. We have developed tools to compute for the first time the radiated fluxes of energy and angular momentum, as well as instantaneous snapshot waveforms, for generic geodesic orbits. For special classes of orbits, we compute the orbital evolution and waveforms for the complete inspiral by imposing global conservation of energy and angular momentum. For fully generic orbits, inspirals and waveforms can be obtained by augmenting our approach with a prescription for the self force in the adiabatic limit derived by Mino. The resulting waveforms should be sufficiently accurate to be used in future gravitational-wave searches.

It has frequently been claimed in the literature that the classical physical predictions of scalar tensor theories of gravity depend on the conformal frame in which the theory is formulated. We argue that this claim is false, and that all classical physical predictions are conformal-frame invariants. We also respond to criticisms by Vollick [gr-qc/0312041], in which this issue arises, of our recent analysis of the Palatini form of 1/R gravity.

We generalize the known equivalence between higher order gravity theories and scalar tensor theories to a new class of theories. Specifically, in the context of a first order or Palatini variational principle where the metric and connection are treated as independent variables, we consider theories for which the Lagrangian density is a function f of (i) the Ricci scalar computed from the metric, and (ii) a second Ricci scalar computed from the connection. We show that such theories can be written as tensor-multi-scalar theories with two scalar fields with the following features: (i) the two dimensional sigma-model metric that defines the kinetic energy terms for the scalar fields has constant, negative curvature; (ii) the coupling function determining the coupling to matter of the scalar fields is universal, independent of the choice of function f; and (iii) if both mass eigenstates are long ranged, then the Eddington post-Newtonian parameter has value 1/2. Therefore in order to be compatible with solar system experiments at least one of the mass eigenstates must be short ranged.

It has been suggested that the recent acceleration of the expansion of the Universe is due a modified gravitational action consisting of the Einstein-Hilbert term plus a term proportional to the reciprocal of the Ricci scalar. Although the original version of this theory has been shown to be in conflict with solar system observations, a modified Palatini version of the theory, in which the metric and connection are treated as independent variables in the variational principle, has been suggested as a viable model of the cosmic acceleration. We show that this theory is equivalent to a type of scalar-tensor theory in which the scalar field kinetic energy term in absent from the action, and in which the scalar field is therefore not an independent dynamical degree of freedom. Integrating out the scalar field gives rise to additional interactions among the matter fields of the standard model of particle physics at a energy scale of order 10^(-3) eV (the geometric mean of the Hubble scale and the Planck scale), and so the theory is excluded by, for example, electron-electron scattering experiments.

Rossby waves (r-modes) in rapidly rotating neutron stars are unstable because of the emission of gravitational radiation. We study saturation of this instability by nonlinear transfer of energy to stellar "inertial" oscillation modes. We present detailed calculations of stellar inertial modes in the WKB limit, their linear damping by bulk and shear viscosity, and the nonlinear coupling forces among these modes. The saturation amplitude is derived in the extreme limits of strong or weak driving by radiation reaction, as compared to the damping rate of low order inertial modes. We find the saturation energy is {\it extremely small}, at least four orders of magnitude smaller than that found by previous investigators. We discuss the consequences of this result for spin evolution of young neutron stars, and neutron stars being spun up by accretion in Low Mass X-ray Binaries.We also discuss the detection of these gravitational waves by LIGO.

We address the issue of finding an optimal detection method for a discontinuous or intermittent gravitational wave stochastic background. Such a signal might sound something like popcorn popping. We derive an appropriate version of the maximum likelihood detection statistic, and compare its performance to that of the standard cross-correlation statistic both analytically and with Monte Carlo simulations. The maximum likelihood statistic performs better than the cross-correlation statistic when the background is sufficiently non-Gaussian. For both ground and space based detectors, this results in a gain factor, ranging roughly from 1 to 3, in the minimum gravitational-wave energy density necessary for detection, depending on the duty cycle of the background. Our analysis is exploratory, as we assume that the time structure of the events cannot be resolved, and we assume white, Gaussian noise in two collocated, aligned detectors. Before this detection method can be used in practice with real detector data, further work is required to generalize our analysis to accommodate separated, misaligned detectors with realistic, colored, non-Gaussian noise.

We develop the formalism required to study the nonlinear interaction of modes in rotating Newtonian stars in the weakly nonlinear regime. The formalism simplifies and extends previous treatments. At linear order, we elucidate and extend slightly a formalism due to Schutz, show how to decompose a general motion of a rotating star into a sum over modes, and obtain uncoupled equations of motion for the mode amplitudes under the influence of an external force. Nonlinear effects are added perturbatively via three-mode couplings. We describe a new, efficient way to compute the coupling coefficients, to zeroth order in the stellar rotation rate, using spin-weighted spherical harmonics. We apply this formalism to derive some properties of the coupling coefficients relevant to the nonlinear interactions of unstable r-modes in neutron stars, postponing numerical integrations of the coupled equations of motion to a later paper. From an astrophysical viewpoint, the most interesting result of this paper is that many couplings of r-modes to other rotational modes (modes with zero frequencies in the non-rotating limit) are small: either they vanish altogether because of various selection rules, or they vanish to lowest order in the angular velocity. In zero-buoyancy stars, the coupling of three r-modes is forbidden entirely and the coupling of two r-modes to one hybrid rotational mode vanishes to zeroth order in rotation frequency. In incompressible stars, the coupling of any three rotational modes vanishes to zeroth order in rotation frequency.

We derive a detection method for a stochastic background of gravitational waves produced by events where the ratio of the average time between events to the average duration of an event is large. Such a signal would sound something like popcorn popping. Our derivation is based on the somewhat unrealistic assumption that the duration of an event is smaller than the detector time resolution.

We elaborate on the recently proposed static brane world scenario, where the effective 4-D cosmological constant is exponentially small when parallel 3-branes are far apart. We extend this result to a compactified model with two positive tension branes. Besides an exponentially small effective 4-D cosmological constant, this model incorporates a Randall-Sundrum-like solution to the hierarchy problem. Furthermore, the exponential factors for the hierarchy problem and the cosmological constant problem obey an inequality that is satisfied in nature. This inequality implies that the cosmological constant problem can be explained if the hierarchy problem is understood. The basic idea generalizes to the multibrane world scenario. We discuss models with piecewise adjustable bulk cosmological constants (to be determined by the 5-dimensional Einstein equation), a key element of the scenario. We also discuss the global structure of this scenario and clarify the physical properties of the particle (Rindler) horizons that are present. Finally, we derive a 4-D effective theory in which all observers on all branes not separated by particle horizons measure the same Newton's constant and 4-D cosmological constant.

We develop a possible cosmology for a Universe with n additional spatial dimensions of variable scale, and an associated scalar field, the radion, which is distinct from the field responsible for inflation, the inflaton. Based on a particular ansatz for the effective potential for the inflaton and radion (which may emerge in string theory), we show that the early expansion of the Universe may proceed in three stages. First, the radion becomes trapped at a value much smaller than the size of the extra dimensions today. Second, the Universe expands exponentially, but with a Planck mass smaller than its present value. Because the Planck mass during inflation is small, density fluctuations in agreement with observations can arise naturally. Third, when inflation ends, the Universe reheats, and the radion becomes free to expand once more. During the third phase the Universe is ``radiation-dominated'' and tends toward a fixed-point evolutionary model in which the radion grows but the temperature remains unchanged. Ultimately, the radion becomes trapped once again at its present value, and a short period of exponential expansion, which we identify with the electroweak phase transition, ensues. Once this epoch is over, the Universe reheats to a temperature of order the electroweak scale, and the mature Universe evolves as in standard cosmological models. We show that the present day energy density in radions can be smaller than the closure density if the second inflationary epoch lasts about 8 e-foldings or more; the present-day radion mass turns out to be small (less than or of the order of an eV, depending on parameters). We argue that although our model envisages considerable time evolution in the Planck mass, substantial spatial fluctuations in Newton's constant are not produced.

A possible explanation is suggested for the controversial star-crushing effect seen in numerical simulations of inspiraling neutron star binaries by Wilson, Mathews and Marronetti (WMM). An apparently incorrect definition of momentum density in the momentum constraint equation used by WMM gives rise to a post-1-Newtonian error in the approximation scheme. We show by means of an analytic, post-1-Newtonian calculation that this error causes an increase of the stars' central densities which is of the order of several percent when the stars are separated by a few stellar radii, in agreement with what is seen in the simulations.

We analyze the coupling between the internal degrees of freedom of neutron stars in a close binary, and the stars' orbital motion. Our analysis is based on the method of matched asymptotic expansions and is valid to all orders in the strength of internal gravity in each star, but is perturbative in the ``tidal expansion parameter'' (stellar radius)/(orbital separation). At first order in the tidal expansion parameter, we show that the internal structure of each star is unaffected by its companion, in agreement with post-1-Newtonian results of Wiseman (gr-qc/9704018). We also show that relativistic interactions that scale as higher powers of the tidal expansion parameter produce qualitatively similar effects to their Newtonian counterparts: there are corrections to the Newtonian tidal distortion of each star, both of which occur at third order in the tidal expansion parameter, and there are corrections to the Newtonian decrease in central density of each star (Newtonian ``tidal stabilization''), both of which are sixth order in the tidal expansion parameter. There are additional interactions with no Newtonian analogs, but these do not change the central density of each star up to sixth order in the tidal expansion parameter. These results, in combination with previous analyses of Newtonian tidal interactions, indicate that (i) there are no large general-relativistic crushing forces that could cause the stars to collapse to black holes prior to the dynamical orbital instability, and (ii) the conventional wisdom with respect to coalescing binary neutron stars as sources of gravitational-wave bursts is correct: namely, the finite-stellar-size corrections to the gravitational waveform will be unimportant for the purpose of detecting the coalescences.

In the coming decade, the LIGO/VIRGO/GEO network of ground-based kilometer-scale laser interferometer gravitational wave detectors will open up a new astronomical window on the Universe: gravitational waves in the frequency band 10 to 10^4 Hz. In addition, if the proposed, 5 million kilometer long, space based interferometer LISA flies, another window will be opened in the frequency band 10^(-4) to 1 Hz. I review the various possible sources that might be detected in these frequency bands, and the information that might be obtainable from observed sources. This article closely parallels my talk at the GR15 conference, but is updated to include discussions of some recent developments.

We estimate the expected signal-to-noise ratios (SNRs) from the three phases (inspiral,merger,ringdown) of coalescing binary black holes (BBHs) for initial and advanced ground-based interferometers (LIGO/VIRGO) and for space-based interferometers (LISA). LIGO/VIRGO can do moderate SNR (a few tens), moderate accuracy studies of BBH coalescences in the mass range of a few to about 2000 solar masses; LISA can do high SNR (of order 10^4) high accuracy studies in the mass range of about 10^5 to 10^8 solar masses. BBHs might well be the first sources detected by LIGO/VIRGO: they are visible to much larger distances (up to 500 Mpc by initial interferometers) than coalescing neutron star binaries (heretofore regarded as the "bread and butter" workhorse source for LIGO/VIRGO, visible to about 30 Mpc by initial interferometers). Low-mass BBHs (up to 50 solar masses for initial LIGO interferometers; 100 for advanced; 10^6 for LISA) are best searched for via their well-understood inspiral waves; higher mass BBHs must be searched for via their poorly understood merger waves and/or their well-understood ringdown waves. A matched filtering search for massive BBHs based on ringdown waves should be capable of finding BBHs in the mass range of about 100 to 700 solar masses out to 200 Mpc (initial LIGO interferometers), and 200 to 3000 solar masses out to about z=1 (advanced interferometers). The required number of templates is of order 6000 or less. Searches based on merger waves could increase the number of detected massive BBHs by a factor of order 10 or more over those found from inspiral and ringdown waves, without detailed knowledge of the waveform shapes, using a "noise monitoring" search algorithm. A full set of merger templates from numerical relativity could further increase the number of detected BBHs by an additional factor of up to 4.

We discuss the extraction of information from detected binary black hole (BBH) coalescence gravitational waves, focusing on the merger phase that occurs after the gradual inspiral and before the ringdown. Our results are: (1) If numerical relativity simulations have not produced template merger waveforms before BBH detections by LIGO/VIRGO, one can band-pass filter the merger waves. For BBHs smaller than about 40 solar masses detected via their inspiral waves, the band pass filtering signal to noise ratio indicates that the merger waves should typically be just barely visible in the noise for initial and advanced LIGO interferometers. (2) We derive an optimized (maximum likelihood) method for extracting a best-fit merger waveform from the noisy detector output; one "perpendicularly projects" this output onto a function space (specified using wavelets) that incorporates our prior knowledge of the waveforms. An extension of the method allows one to extract the BBH's two independent waveforms from outputs of several interferometers. (3) If numerical relativists produce codes for generating merger templates but running the codes is too expensive to allow an extensive survey of the merger parameter space, then a coarse survey of this parameter space, to determine the ranges of the several key parameters and to explore several qualitative issues which we describe, would be useful for data analysis purposes. (4) A complete set of templates could be used to test the nonlinear dynamics of general relativity and to measure some of the binary parameters. We estimate the number of bits of information obtainable from the merger waves (about 10 to 60 for LIGO/VIRGO, up to 200 for LISA), estimate the information loss due to template numerical errors or sparseness in the template grid, and infer approximate requirements on template accuracy and spacing.

A recent study of the rotation of the plane of polarization of light from 160 cosmological sources claims to find significant evidence for cosmological anisotropy. We point out methodological weaknesses of that study, and reanalyze the same data using Bayesian methods that overcome these problems. We find that the data always favor isotropic models for the distribution of observed polarizations over counterparts that have a cosmological anisotropy of the type advocated in the earlier study. Although anisotropic models are not completely ruled out, the data put strong lower limits on the length scale $\lambda$ (in units of the Hubble length) associated with the anisotropy; the lower limits of 95% credible regions for $\lambda$ lie between 0.43 and 0.62 in all anisotropic models we studied, values several times larger than the best-fit value of $\lambda \approx 0.1$ found in the earlier study. The length scale is not constrained from above. The vast majority of sources in the data are at distances closer than 0.4 Hubble lengths (corresponding to a redshift of $\approx$0.8); the results are thus consistent with there being no significant anisotropy on the length scale probed by these data.

The most promising source of gravitational waves for the planned detectors LIGO and VIRGO are merging compact binaries, i.e., neutron star/neutron star (NS/NS), neutron star/black hole (NS/BH), and black hole/black-hole (BH/BH) binaries. We investigate how accurately the distance to the source and the masses and spins of the two bodies will be measured from the gravitational wave signals by the three detector LIGO/VIRGO network using ``advanced detectors'' (those present a few years after initial operation). The combination ${\cal M} \equiv (M_1 M_2)^{3/5}(M_1 +M_2)^{-1/5}$ of the masses of the two bodies is measurable with an accuracy $\approx 0.1\%-1\%$. The reduced mass is measurable to $\sim 10\%-15\%$ for NS/NS and NS/BH binaries, and $\sim 50\%$ for BH/BH binaries (assuming $10M_\odot$ BH's). Measurements of the masses and spins are strongly correlated; there is a combination of $\mu$ and the spin angular momenta that is measured to within $\sim 1\%$. We also estimate that distance measurement accuracies will be $\le 15\%$ for $\sim 8\%$ of the detected signals, and $\le 30\%$ for $\sim 60\%$ of the signals, for the LIGO/VIRGO 3-detector network.

We analyze the sensitivity of a network of interferometer gravitational-wave detectors to the gravitational-wave stochastic background, and derive the dependence of this sensitivity on the orientations of the detector arms. We build on and extend the recent work of Christensen, but our conclusion for the optimal choice of orientations of a pair of detectors differs from his. For a pair of detectors (such as LIGO) that subtends an angle at the center of the earth of $\,\alt 70^\circ$, we find that the optimal configuration is for each detector to have its arms make an angle of $45^\circ$ (modulo $90^\circ$) with the arc of the great circle that joins them. For detectors that are farther separated, each detector should instead have one arm aligned with this arc. We also describe in detail the optimal data-analysis algorithm for searching for the stochastic background with a detector network, which is implicit in earlier work of Michelson. The LIGO pair of detectors will be separated by $\sim 3000 \, {\rm km}$. The minimum detectable stochastic energy-density for these detectors with their currently planned orientations is $\sim 3\%$ greater than what it would be if the orientations were optimal.

Gravitational-wave interferometers are expected to monitor the last three minutes of inspiral and final coalescence of neutron star and black hole binaries at distances approaching cosmological, where the event rate may be many per year. Because the binary's accumulated orbital phase can be measured to a fractional accuracy $\ll 10^{-3}$ and relativistic effects are large, the waveforms will be far more complex, carry more information, and be far harder to model theoretically than has been expected. Theorists must begin now to lay a foundation for extracting the waves' information.