We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$.

The hierarchical nature of galaxy formation suggests that a supermassive black hole binary could exist in our galactic center. We propose a new approach to constraining the possible orbital configuration of such a binary companion to the galactic center black hole Sgr A* through the measurement of stellar orbits. Focusing on the star S0-2, we show that requiring its orbital stability in the presence of a companion to Sgr A* yields stringent constraints on the possible configurations of such a companion. Furthermore, we show that precise measurements of {\it time variations} in the orbital parameters of S0-2 could yield stronger constraints. Using existing data on S0-2 we derive upper limits on the binary black hole separation as a function of the companion mass. For the case of a circular orbit, we can rule out a 10^5 M_sun companion with a semimajor axis greater than 170 astronomical units or 0.8 mpc. This is already more stringent than bounds obtained from studies of the proper motion of Sgr A*. Including other stars orbiting the galactic center should yield stronger constraints that could help uncover the presence of a companion to Sgr A*. We show that a companion can also affect the accretion process, resulting in a variability which may be consistent with the measured infrared flaring timescales and amplitudes. Finally, if such a companion exists, it will emit gravitational wave radiation, potentially detectable with LISA.

The growth of a massive black hole will steepen the cold dark matter density at the center of a galaxy into a dense spike, enhancing the prospects for indirect detection. We study the impact of black hole spin on the density profile using the exact Kerr geometry of the black whole in a fully relativistic adiabatic growth framework. We find that, despite the transfer of angular momentum from the hole to the halo, rotation increases significantly the dark matter density close to the black hole. The gravitational effects are still dominated by the black hole within its influence radius, but the larger dark matter annihilation fluxes might be relevant for indirect detection estimates.

We analyze the secular evolution of hierarchical triple systems in the post-Newtonian approximation to general relativity. We expand the Newtonian three-body equations of motion in powers of the ratio $a/A$, where $a$ and $A$ are the semimajor axis of the inner binary's orbit and of the orbit of the third body relative to the center of mass of the inner binary, respectively. The leading order "quadrupole" terms, of order $(a/A)^3$ relative to the $1/a^2$ acceleration within the inner binary, are responsible for the well-known Kozai-Lidov oscillations of orbital inclination and eccentricity. The octupole terms, of order $(a/A)^4$ have been shown to allow the inner orbit to "flip" from prograde relative to the outer orbit to retrograde and back, and to permit excursions to very large eccentricities. We carry the expansion of the equations of motion to hexadecapole order, corresponding to contributions of order $(a/A)^5$. We also include the leading orbital effects of post-Newtonian theory, namely the pericenter precessions of the inner and outer orbits. Using the Lagrange planetary equations for the orbit elements of both binaries, we average over orbital timescales, obtain the equations for the secular evolution of the elements through hexadecapole order, and employ them to analyze cases of astrophysical interest. We find that, for the most part, the orbital flips found at octupole order are robust against both relativistic and hexadecapole perturbations. We show that, for equal-mass inner binaries, where the octupole terms vanish, the hexadecapole contributions can alone generate orbital flips and excursions to very large eccentricities.

We derive the secular evolution of the orbital elements of a stellar-mass object orbiting a spinning massive black hole. We use the post-Newtonian approximation in harmonic coordinates, with test-body equations of motion for the conservative dynamics that are valid through 3PN order, including spin-orbit, quadrupole and (spin)$^2$ effects, and with radiation-reaction contributions linear in the mass of the body that are valid through 4.5PN order, including the 4PN damping effects of spin-orbit coupling. The evolution equations for the osculating orbit elements are iterated to high PN orders using a two-timescale approach and averaging over orbital timescales. We derive a criterion for terminating the orbit when its Carter constant drops below a critical value, whereupon the body plunges across the event horizon at the next closest approach. The results are valid for arbitrary eccentricities and arbitrary inclinations. We then analyze numerically the orbits of objects injected into high-eccentricity orbits via interactions within a surrounding star cluster, obtaining the number of orbits and the elapsed time between injection and plunge, and the residual orbital eccentricity at plunge as a function of inclination. We derive an analytic approximation for the time to plunge in terms of initial orbital variables. We show that, if the black hole is spinning rapidly, the flux of gravitational radiation during the final orbit before plunge may be suppressed by as much as three orders of magnitude if the orbit is retrograde on the equatorial plane compared to its prograde counterpart.

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

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

Conventional approaches to incorporating general relativistic effects into the dynamics of N-body systems containing central black holes, or of hierarchical triple systems with a relativistic inner binary, may not be adequate when the goal is to study the evolution of the system over a timescale related to relativistic secular effects, such as the precession of the pericenter. For such problems, it may necessary to include post-Newtonian "cross terms" in the equations of motion in order to capture relativistic effects consistently over the long timescales. Cross terms are post-Newtonian (PN) terms that explicitly couple the two-body relativistic perturbations with the Newtonian perturbations due to other bodies in the system. In this paper, we show that the total energy and the normal component of total angular momentum of a hierarchical triple system is manifestly conserved to Newtonian order over the relativistic pericenter precession timescale of the inner binary if and only if PN cross-term effects in the equations of motion are taken carefully into account.

The increasing role of general relativity in the dynamics of stellar systems with central massive black holes and in the evolution of hierarchical triple systems inspires a close examination of how post-Newtonian effects are incorporated into N-body dynamics. The majority of approaches incorporate relativity by adding to the Newtonian N-body equations the standard two-body post-Newtonian terms for a given star around the black hole or for the close binary in a triple system. We argue that, for calculating the evolution of such systems over timescales comparable to the relativistic pericenter advance timescale, it is essential to include "cross terms" in the equations of motion. These are post-Newtonian terms that represent a coupling between the potential of the central black hole and the potential due to other stars in the system. For hierarchical triple systems, these are couplings between the potential of the inner binary and that of the distant third body. Over pericenter precession timescales, the effects of such terms can actually be "boosted" to amplitudes of Newtonian order. We write down the post-Newtonian N-body equations of motion including a central black hole in a truncated form that includes all the relevant cross terms, in a format ready to use for numerical implementation. We do the same for hierarchical triple systems, and illustrate explicitly the effects of cross terms on the orbit-averaged equations of evolution for the orbit elements of the inner binary for the special case where the third body is on a circular orbit. We also describe the inspiration for this investigation: the motion of a test body about a central body with a Newtonian quadrupole moment, including the relativistic pericenter advance, whose correct solution for the conserved total Newtonian energy requires including PN cross terms between the mass monopole and quadrupole potentials.

The status of experimental tests of general relativity and of theoretical frameworks for analyzing them are reviewed and updated. Einstein's equivalence principle (EEP) is well supported by experiments such as the Eotvos experiment, tests of local Lorentz invariance and clock experiments. Ongoing tests of EEP and of the inverse square law are searching for new interactions arising from unification or quantum gravity. Tests of general relativity at the post-Newtonian level have reached high precision, including the light deflection, the Shapiro time delay, the perihelion advance of Mercury, the Nordtvedt effect in lunar motion, and frame-dragging. Gravitational-wave damping has been detected in an amount that agrees with general relativity to better than half a percent using the Hulse-Taylor binary pulsar, and a growing family of other binary pulsar systems is yielding new tests, especially of strong-field effects. Current and future tests of relativity will center on strong gravity and gravitational waves.

The cold dark matter at the center of a galaxy will be redistributed by the presence of a massive black hole. The redistribution may be determined using an approach pioneered by Gondolo and Silk: begin with a model distribution function for the dark matter, and ``grow'' the black hole adiabatically, holding the adiabatic invariants of the motion constant. Unlike the approach of Gondolo and Silk, which adopted Newtonian theory together with ad hoc correction factors to mimic general relativistic effects, we carry out the calculation fully relativistically, using the exact Schwarzschild geometry of the black hole. We find that the density of dark matter generically vanishes at r=2R_S, not 4R_S as found by Gondolo and Silk, where R_S is the Schwarzschild radius, and that the spike very close to the black hole reaches significantly higher densities. We apply the relativistic adiabatic growth framework to obtain the final dark matter density for both cored and cusped initial distributions. Besides the implications of these results for indirect detection estimates, we show that the gravitational effects of such a dark matter spike are significantly smaller than the relativistic effects of the black hole, including frame dragging and quadrupolar effects, for stars orbiting close to the black hole that might be candidates for testing the black hole no-hair theorems.

We obtain approximate analytic expressions for the critical value of the total angular momentum of a non-relativistic test particle moving in the Kerr geometry, such that it will be captured by the black hole. The expressions apply to arbitrary orbital inclinations, and are accurate over the entire range of angular momentum for the Kerr black hole. The expressions can be easily implemented in N-body simulations of the evolution of star clusters around massive galactic black holes, where such captures play an important role.

We derive the equations of motion, the periastron shift, and the gravitational radiation damping for quasicircular compact binaries in a massive variant of the Brans-Dicke theory of gravity. We also study the Shapiro time delay and the Nordtvedt effect in this theory. By comparing with recent observational data, we put bounds on the two parameters of the theory: the Brans-Dicke coupling parameter \omega_{BD} and the scalar mass m_s. We find that the most stringent bounds come from Cassini measurements of the Shapiro time delay in the Solar System, that yield a lower bound \omega_{BD}>40000 for scalar masses m_s<2.5x10^{-20} eV, to 95% confidence. In comparison, observations of the Nordtvedt effect using Lunar Laser Ranging (LLR) experiments yield \omega_{BD}>1000 for m_s<2.5x10^{-20} eV. Observations of the orbital period derivative of the quasicircular white dwarf-neutron star binary PSR J1012+5307 yield \omega_{BD}>1250 for m_s<10^{-20} eV. A first estimate suggests that bounds comparable to the Shapiro time delay may come from observations of radiation damping in the eccentric white dwarf-neutron star binary PSR J1141-6545, but a quantitative prediction requires the extension of our work to eccentric orbits.

Observations of the precessing orbits of stars very near the massive black hole in the galactic center could provide measurements of the spin and quadrupole moment of the hole and thereby test the no-hair theorem of general relativity. Since the galactic center is likely to be populated by a distribution of stars and small black holes, their gravitational interactions will perturb the orbit of any given star. We estimate the effects of such perturbations using analytic orbital perturbation theory, and show that for a range of possible stellar distributions, and for an observed star sufficiently close to the black hole, the relativistic spin and quadrupole effects will be larger than the effects of stellar cluster perturbations. Our results are consistent those from recent numerical N-body simulations by Merritt et al.

Inspiral of compact stellar remnants into massive black holes (MBHs) is accompanied by the emission of gravitational waves at frequencies that are potentially detectable by space-based interferometers. Event rates computed from statistical (Fokker-Planck, Monte-Carlo) approaches span a wide range due to uncertaintities about the rate coefficients. Here we present results from direct integration of the post-Newtonian N-body equations of motion descrbing dense clusters of compact stars around Schwarzschild MBHs. These simulations embody an essentially exact (at the post-Newtonian level) treatment of the interplay between stellar dynamical relaxation, relativistic precession, and gravitational-wave energy loss. The rate of capture of stars by the MBH is found to be greatly reduced by relativistic precession, which limits the ability of torques from the stellar potential to change orbital angular momenta. Penetration of this "Schwarzschild barrier" does occasionally occur, resulting in capture of stars onto orbits that gradually inspiral due to gravitational wave emission; we discuss two mechanisms for barrier penetration and find evidence for both in the simulations. We derive an approximate formula for the capture rate, which predicts that captures would be strongly disfavored from orbits with semi-major axes below a certain value; this prediction, as well as the predicted rate, are verified in the N-body integrations. We discuss the implications of our results for the detection of extreme-mass-ratio inspirals from galactic nuclei with a range of physical properties.

The post-Newtonian approximation is a method for solving Einstein's field equations for physical systems in which motions are slow compared to the speed of light and where gravitational fields are weak. Yet it has proven to be remarkably effective in describing certain strong-field, fast-motion systems, including binary pulsars containing dense neutron stars and binary black hole systems inspiraling toward a final merger. The reasons for this effectiveness are largely unknown. When carried to high orders in the post-Newtonian sequence, predictions for the gravitational-wave signal from inspiraling compact binaries will play a key role in gravitational-wave detection by laser-interferometric observatories.

The gravitational recoil or "kick" of a black hole formed from the merger of two orbiting black holes, and caused by the anisotropic emission of gravitational radiation, is an astrophysically important phenomenon. We combine (i) an earlier calculation, using post-Newtonian theory, of the kick velocity accumulated up to the merger of two non-spinning black holes, (ii) a "close-limit approximation" calculation of the radiation emitted during the ringdown phase, and based on a solution of the Regge-Wheeler and Zerilli equations using initial data accurate to second post-Newtonian order. We prove that ringdown radiation produces a significant "anti-kick". Adding the contributions due to inspiral, merger and ringdown phases, our results for the net kick velocity agree with those from numerical relativity to 10-15 percent over a wide range of mass ratios, with a maximum velocity of 180 km/s at a mass ratio of 0.38.

The spin and quadrupole moment of the supermassive black hole at the Galactic center can in principle be measured via astrometric monitoring of stars orbiting at milliparsec (mpc) distances, allowing tests of general relativistic "no-hair" theorems (Will 2008). One complicating factor is the presence of perturbations from other stars, which may induce orbital precession of the same order of magnitude as that due to general relativistic effects. The expected number of stars in this region is small enough that full N-body simulations can be carried out. We present the results of a comprehensive set of such simulations, which include a post-Newtonian treatment of spin-orbit effects. A number of possible models for the distribution of stars and stellar remnants are considered. We find that stellar perturbations are likely to obscure the signal due to frame-dragging for stars beyond ~0.5 mpc from the black hole, while measurement of the quadrupole moment is likely to require observation of stars inside ~0.2 mpc. A high fraction of stellar remnants, e.g. 10-Solar-mass black holes, in this region would make tests of GR problematic at all radii. We discuss the possibility of separating the effects of stellar perturbations from those due to GR.

Spin induced precessional modulations of gravitational wave signals from supermassive black hole binaries can improve the estimation of luminosity distance to the source by space based gravitational wave missions like the Laser Interferometer Space Antenna (LISA). We study how this impacts the ablity of LISA to do cosmology, specifically, to measure the dark energy equation of state (EOS) parameter $w$. Using the $\Lambda$CDM model of cosmology, we show that observations of precessing binaries by LISA, combined with a redshift measurement, can improve the determination of $w$ up to an order of magnitude with respect to the non precessing case depending on the masses, mass ratio and the redshift.

If a class of stars orbits the central black hole in our galaxy in short period (~ 0.1 year), high eccentricity (~ 0.9) orbits, they will experience precessions of their orbital planes induced by both relativistic frame-dragging and the quadrupolar gravity of the hole, at levels that could be as large as 10 microarcseconds per year, if the black hole is rotating faster than 1/2 of its maximum rotation rate. Astrometric observations of the orbits of at least two such stars can in principle lead to a determination of the angular momentum vector J of the black hole and its quadrupole moment Q_2. This could lead to a test of the general relativistic no-hair theorems, which demand that Q_2 = - J^2/M. Future high-precision adaptive infrared optics instruments make make such a fundamental test of the black-hole paradigm possible.

We use a post-Newtonian diagnostic tool to examine numerically generated quasiequilibrium initial data sets for non-spinning double neutron star and neutron star-black hole binary systems. The PN equations include the effects of tidal interactions, parametrized by the compactness of the neutron stars and by suitable values of ``apsidal'' constants, which measure the degree of distortion of stars subjected to tidal forces. We find that the post-Newtonian diagnostic agrees well with the double neutron star initial data, typically to better than half a percent except where tidal distortions are becoming extreme. We show that the differences could be interpreted as representing small residual eccentricity in the initial orbits. In comparing the diagnostic with preliminary numerical data on neutron star-black hole binaries, we find less agreement.

Newly formed black holes are expected to emit characteristic radiation in the form of quasi-normal modes, called ringdown waves, with discrete frequencies. LISA should be able to detect the ringdown waves emitted by oscillating supermassive black holes throughout the observable Universe. We develop a multi-mode formalism, applicable to any interferometric detectors, for detecting ringdown signals, for estimating black hole parameters from those signals, and for testing the no-hair theorem of general relativity. Focusing on LISA, we use current models of its sensitivity to compute the expected signal-to-noise ratio for ringdown events, the relative parameter estimation accuracy, and the resolvability of different modes. We also discuss the extent to which uncertainties on physical parameters, such as the black hole spin and the energy emitted in each mode, will affect our ability to do black hole spectroscopy.

We provide some considerations on the excitation of black hole quasinormal modes (QNMs) in different physical scenarios. Considering a simple model in which a stream of particles accretes onto a black hole, we show that resonant QNM excitation by hyperaccretion requires a significant amount of fine-tuning, and is quite unlikely to occur in nature. Then we summarize and discuss present estimates of black hole QNM excitation from gravitational collapse, distorted black holes and head-on black hole collisions. We emphasize the areas that, in our opinion, are in urgent need of further investigation from the point of view of gravitational wave source modeling.

The loss of linear momentum by gravitational radiation and the resulting gravitational recoil of black-hole binary systems may play an important role in the growth of massive black holes in early galaxies. We calculate the gravitational recoil of non-spinning black-hole binaries at the second post-Newtonian order (2PN) beyond the dominant effect, obtaining, for the first time, the 1.5PN correction term due to tails of waves and the next 2PN term. We find that the maximum value of the net recoil experienced by the binary due to the inspiral phase up to the innermost stable circular orbit (ISCO) is of the order of 22 km/s. We then estimate the kick velocity accumulated during the plunge from the ISCO up to the horizon by integrating the momentum flux using the 2PN formula along a plunge geodesic of the Schwarzschild metric. We find that the contribution of the plunge dominates over that of the inspiral. For a mass ratio m_2/m_1=1/8, we estimate a total recoil velocity (due to both adiabatic and plunge phases) of 100 +/- 20 km/s. For a ratio 0.38, the recoil is maximum and we estimate it to be 250 +/- 50 km/s. In the limit of small mass ratio, we estimate V/c to be approximately 0.043 (1 +/- 20%)(m_2/m_1)^2. Our estimates are consistent with, but span a substantially narrower range than, those of Favata et al. (2004).

Observations of binary inspirals with LISA will allow us to place bounds on alternative theories of gravity and to study the merger history of massive black holes (MBH). These possibilities rely on LISA's parameter estimation accuracy. We update previous studies of parameter estimation including non-precessional spin effects. We work both in Einstein's theory and in alternative theories of gravity of the scalar-tensor and massive-graviton types. Inclusion of non-precessional spin terms in MBH binaries has little effect on the angular resolution or on distance determination accuracy, but it degrades the estimation of the chirp mass and reduced mass by between one and two orders of magnitude. The bound on the coupling parameter of scalar-tensor gravity is significantly reduced by the presence of spin couplings, while the reduction in the graviton-mass bound is milder. LISA will measure the luminosity distance of MBHs to better than ~10% out to z~4 for a (10^6+10^6) Msun binary, and out to z~2 for a (10^7+10^7) Msun binary. The chirp mass of a MBH binary can always be determined with excellent accuracy. Ignoring spin effects, the reduced mass can be measured within ~1% out to z=10 and beyond for a (10^6+10^6) Msun binary, but only out to z~2 for a (10^7+10^7) Msun binary. Present-day MBH coalescence rate calculations indicate that most detectable events should originate at z~2-6: at these redshifts LISA can be used to measure the two black hole masses and their luminosity distance with sufficient accuracy to probe the merger history of MBHs. If the low-frequency LISA noise can only be trusted down to 10^-4 Hz, parameter estimation for MBHs (and LISA's ability to perform reliable cosmological observations) will be significantly degraded.

We investigate the effect of spin-orbit and spin-spin couplings on the estimation of parameters for inspiralling compact binaries of massive black holes, and for neutron stars inspiralling into intermediate-mass black holes, using hypothetical data from the proposed Laser Interferometer Space Antenna (LISA). We work both in Einstein's theory and in alternative theories of gravity of the scalar-tensor and massive-graviton types. We restrict the analysis to non-precessing spinning binaries, i.e. to cases where the spins are aligned normal to the orbital plane. We find that the accuracy with which intrinsic binary parameters such as chirp mass and reduced mass can be estimated within general relativity is degraded by between one and two orders of magnitude. We find that the bound on the coupling parameter omega_BD of scalar-tensor gravity is significantly reduced by the presence of spin couplings, while the reduction in the graviton-mass bound is milder. Using fast Monte-Carlo simulations of 10^4 binaries, we show that inclusion of spin terms in massive black-hole binaries has little effect on the angular resolution or on distance determination accuracy. For stellar mass inspirals into intermediate-mass black holes, the angular resolution and the distance are determined only poorly, in all cases considered. We also show that, if LISA's low-frequency noise sensitivity can be extrapolated from 10^-4 Hz to as low as 10^-5 Hz, the accuracy of determining both extrinsic parameters (distance, sky location) and intrinsic parameters (chirp mass, reduced mass) of massive binaries may be greatly improved.

We estimate the rate at which the proposed space gravitational-wave interferometer LISA could detect intermediate-mass black-hole binaries, that is, binaries containing a black hole in the mass range 10 -- 100 Msun orbiting a black hole in the mass range 100 -- 1000 Msun. For one-year integrations leading up to the innermost stable orbit, and a signal-to-noise ratio of 10, we estimate a detection rate of only 1 per million years for 10 Msun/100 Msun binaries. The estimate uses the method of parameter estimation via matched filtering, incorporates a noise curve for LISA established by the LISA Science Team that is available online, and employs an IMBH formation rate model used by Miller (2002). We find that the detectable distance is relatively insensitive to LISA arm lengths or acceleration noise, but is roughly inversely proportional to LISA position errors, and varies roughly as T^(1/2), where T is the integration time in years. We also show that, while all IMBH systems in this mass range may be detected in the Virgo cluster up to 40 years before merger, none can be detected there earlier than 400 years before merger. An extended LISA mission that enabled 10-year integrations could detect IMBH systems at the Virgo cluster 1000 years before merger, and systems in galactic globular clusters a million years before merger. We compare and contrast these estimates with earlier estimates by Miller (2002).

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before.

We calculate the delay in the propagation of a light signal past a massive body that moves with speed v, under the assumption that the speed of propagation of the gravitational interaction c_g differs from that of light. Using the post-Newtonian approximation, we consider an expansion in powers of v/c beyond the leading ``Shapiro'' time delay effect, while working to first order only in Gm/c^2, and show that the altered propagation speed of the gravitational signal has no effect whatsoever on the time delay to first order in v/c beyond the leading term, although it will have an effect to second and higher order. We show that the only other possible effects of an altered speed c_g at this order arise from a modification of the parametrized post-Newtonian (PPN) coefficient \alpha_1 of the metric from the value zero predicted by general relativity. Current solar-system measurements already provide tight bounds on such a modification. We conclude that recent measurements of the propagation of radio signals past Jupiter are sensitive to \alpha_1, but are not directly sensitive to the speed of propagation of gravity.

We examine the possibility that an epoch of inflationary expansion induces high-frequency oscillations of Newton's constant, $G$. The effect occurs because inflation can shift the expectation value of a non-minimally coupled, Brans-Dicke-like field away from the minimum of its effective potential. At some time after inflation ends, the field begins to oscillate, resulting in periodic variations in $G$. We find conditions for which the oscillation energy would be sufficient to close the universe, consistent with all known constraints from cosmology and local tests of general relativity.