The ability to make independent detections of the signatures of exoplanets with complementary telescopes and instruments brings a new potential for robust identification of exoplanets and precision characterization. We introduce PEXO, a package for Precise EXOplanetology to facilitate the efficient modeling of timing, astrometry, and radial velocity data, which will benefit not only exoplanet science but also various astrophysical studies in general. PEXO is general enough to account for binary motion and stellar reflex motions induced by planetary companions and is precise enough to treat various relativistic effects both in the solar system and in the target system. We also model the post-Newtonian barycentric motion for future tests of general relativity in extrasolar systems. We benchmark PEXO with the pulsar timing package TEMPO2 and find that PEXO produces numerically similar results with timing precision of about 1 ns, space-based astrometry to a precision of 1 $\mu$as, and radial velocity of 1 $\mu$m/s and improves on TEMPO2 for decade-long timing data of nearby targets, due to its consideration of third-order terms of Roemer delay. PEXO is able to avoid the bias introduced by decoupling the target system and the solar system and to account for the atmospheric effects which set a practical limit for ground-based radial velocities close to 1 cm/s. Considering the various caveats in barycentric correction and ancillary data required to realize cm/s modeling, we recommend the preservation of original observational data. The PEXO modeling package is available at GitHub (https://github.com/phillippro/pexo).

Our present relativistic part of the geodetic VLBI model for Earthbound antennas is a consensus model which is considered as a standard for processing high-precision VLBI observations. It was created as a compromise between a variety of relativistic VLBI models proposed by different authors as documented in the IERS Conventions 2010. The accuracy of the consensus model is in the picosecond range for the group delay but this is not sufficient for current geodetic pur- poses. This paper provides a fully documented derivation of a new relativistic model having an accuracy substantially higher than one picosecond and based upon a well accepted formalism of relativistic celestial mechanics, astrometry and geodesy. Our new model fully confirms the consensus model at the picosecond level and in several respects goes to a great extent beyond it. More specifically, terms related to the acceleration of the geocenter are considered and kept in the model, the gravitational time-delay due to a massive body (planet, Sun, etc.) with arbitrary mass and spin-multipole moments is derived taking into account the motion of the body, and a new formalism for the time-delay problem of radio sources located at finite distance from VLBI stations is presented. Thus, the paper presents a substantially elaborated theoretical justification of the consensus model and its significant extension that allows researchers to make concrete estimates of the magnitude of residual terms of this model for any conceivable configuration of the source of light, massive bodies, and VLBI stations. The largest terms in the relativistic time delay which can affect the current VLBI observations are from the quadrupole and the angular momentum of the gravitating bodies that are known from the literature. These terms should be included in the new geodetic VLBI model for improving its consistency.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We study physical consequences of the Hubble expansion of FLRW manifold on measurement of space, time and light propagation in the local inertial frame. We analyse the solar system radar ranging and Doppler tracking experiments, and time synchronization. FLRW manifold is covered by global coordinates (t,y^i), where t is the cosmic time coinciding with the proper time of the Hubble observers. We introduce local inertial coordinates x^a=(x^0,x^i) in the vicinity of a world line of a Hubble observer with the help of a special conformal transformation. The local inertial metric is Minkowski flat and is materialized by the congruence of time-like geodesics of static observers being at rest with respect to the local spatial coordinates x^i. We consider geodesic motion of test particles and notice that the local coordinate time x^0=x^0(t) taken as a parameter along the world line of particle, is a function of the Hubble's observer time t. This function changes smoothly from x^0=t for a particle at rest (observer's clock), to x^0=t+1/2 Ht^2 for photons, where H is the Hubble constant. Thus, motion of a test particle is non-uniform when its world line is parametrized by time t. NASA JPL Orbit Determination Program presumes that motion of light (after the Shapiro delay is excluded) is uniform with respect to the time t but it does not comply with the non-uniform motion of light on cosmological manifold. For this reason, the motion of light in the solar system analysed with the Orbit Determination Program appears as having a systematic blue shift of frequency, of radio waves circulating in the Earth-spacecraft radio link. The magnitude of the anomalous blue shift of frequency is proportional to the Hubble constant H that may open an access to the measurement of this fundamental cosmological parameter in the solar system radiowave experiments.

The extremely high precision of current astronomical observations demands a much better theoretical treatment of relativistic effects in the propagation of electromagnetic signals through variable gravitational fields of isolated astronomical systems emitting gravitational waves. This paper presents a consistent approach giving a complete and exhaustive solution of this problem in the first post-Minkowskian approximation of general relativity.

We study physical consequences of the Einstein equivalence principle (EEP) for a Hubble observer in FLRW universe. We introduce the local inertial coordinates with the help of a special conformal transformation. The local inertial metric is Minkowski-flat and materialized by a congruence of time-like geodesics of static observers. The static observers are equipped with the ideal clocks measuring the proper time that is synchronized with the clocks of the Hubble observer. The local inertial metric is used for physical measurements of spacetime intervals with the ideal clocks and rulers. The special conformal transformation preserves null geodesics but does not keep invariant time-like geodesics. Moreover, it makes the rate of the local time coordinate dependent on velocity of the particle which makes impossible to rich the uniform parameterization of the world lines of static observers and light geodesics with a single parameter - they differ by the conformal factor of FLRW metric. It tells us that the metric on the light cone is not Minkowski-flat but depends on the scale factor of FLRW universe and it can be interpreted as a weak violation of EEP for photons. The importance of this violation for gravitational physics is that some of local experiments conducted with freely-propagating electromagnetic waves may be sensitive to the Hubble expansion. We show that the Hubble constant H can be measured within the solar system by means of high-precision spacecraft Doppler tracking as a blue shift of frequency of radio waves circulating in the Earth-spacecraft radio link. We also analyze the behavior of the standing wave in a microwave resonator and show that the standing wave is insensitive to the Hubble expansion.

We analyse the Einstein equivalence principle (EEP) for a Hubble observer in Friedmann-Lemaitre-Robertson-Walker spacetime. We show that the affine structure of light cone in the FLRW spacetime should be treated locally in terms of the optical metric which is not reduced to the Minkowski metric due to the non-uniform parametrization of the local equations of light propagation with the proper time of the observer's clock. The physical consequence of this difference is that the Doppler shift of radio waves measured locally, is affected by the Hubble expansion.

The present paper outlines theoretical principles of the post-Newtonian mechanics in the expanding universe. It is based upon the gauge-invariant theory of the Lagrangian perturbations of cosmological manifold caused by an isolated astronomical N-body system. We postulate that the background manifold is described by Friedman-Lemaitre-Robertson-Walker (FLRW) metric governed by two primary components - the dark matter and the dark energy. The dark matter is treated as an ideal fluid. The dark energy is described by a single scalar field with a potential which is hold unspecified as long as the theory permits. The Lagrangian of the dark matter and that of the scalar field are formulated in terms of the field variables. We use variational calculus to derive the gauge-invariant field equations of the post-Newtonian celestial mechanics of an isolated astronomical system in an expanding universe. These equations generalize the field equations of the post-Newtonian theory in asymptotically-flat spacetime by taking into account the cosmological effects explicitly. We introduce a new cosmological gauge which generalizes the harmonic gauge of the post-Newtonian theory in asymptotically-flat spacetime. This gauge significantly simplifies the gravitational field equations and allows finding out the approximations where the field equations can be fully decoupled and solved analytically. The residual gauge freedom is explored. The results of the present paper can be useful in the solar system for calculating more precise ephemerides of the solar system bodies on extremely long time intervals, in galactic astronomy to study the dynamics of clusters of galaxies, and in gravitational wave astronomy for discussing the impact of cosmology on generation and propagation of gravitational waves emitted by coalescing binaries and/or merging galactic nuclei.

The mathematical concept of the Newtonian limit of Einstein's field equations in the expanding Friedmann universe is formulated. The geodesic equations of motion of planets and light are derived and compared.

Post-Newtonian theory was instrumental in conducting the critical experimental tests of general relativity and in building the astronomical ephemerides of celestial bodies in the solar system with an unparalleled precision. The cornerstone of the theory is the postulate that the solar system is gravitationally isolated from the rest of the universe and the background spacetime is asymptotically flat. The present article extends this theoretical concept and formulates the principles of celestial dynamics of particles and light moving in gravitational field of a localized astronomical system embedded to the expanding Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. We formulate the precise mathematical concept of the Newtonian limit of Einstein's field equations in the conformally-flat FLRW spacetime and analyze the geodesic motion of massive particles and light in this limit. We prove that by doing conformal spacetime transformations, one can reduce the equations of motion of particles and light to the classical form of the Newtonian theory. However, the time arguments in the equations of motion of particles and light differ from each other in terms being proportional to the Hubble constant, H. This leads to the important conclusion that the equations of light propagation used currently by Space Navigation Centers for fitting range and Doppler-tracking observations of celestial bodies are missing some terms of the cosmological origin that are proportional to the Hubble constant, H. We also prove that the Hubble expansion does not affect the atomic time scale used in creation of astronomical ephemerides. We derive the cosmological correction to the light travel time equation and argue that their measurement opens an exciting opportunity to determine the local value of the Hubble constant, H, in the solar system independently of cosmological observations.

IAU 2000 resolutions on the reference frames set up a solid theoretical foundation for implementing general relativity in astronomical data processing algorithms and for unambiguous interpretation of measured relativistic effects. We discuss possible directions for further theoretical development of the IAU resolutions aimed to take into account the decadal progress in observational techniques and computer-based technologies. We address the following subjects: 1) space-time transformations and the structure of the metric tensor; -2) PPN parameters and gauge invariance of equations of motion; -3) astronomical reference frames for cosmological applications.

Since VLBI techniques give microarcsecond position accuracy of celestial objects, tests of GR using radio sources as probes of a gravitational field have been made. We present the results from two recent tests using the VLBA: In 2005, the measurement of the classical solar deflection; and in 2002, the measurement of the retarded gravitational deflection associated with Jupiter. The deflection experiment measured PPN-gamma to an accuracy of 0.0003; the Jupiter experiment measured the retarded term to 20% accuracy. The controversy over the interpretation of the retarded term is summarized.

We report on recent VLBA/VERA/IVS observational tests of General Relativity. First, we will summarize the results from the 2005 VLBA experiment that determined gamma with an accuracy of 0.0003 by measuring the deflection of four compact radio sources by the solar gravitational field. We discuss the limits of precision that can be obtained with VLBA experiments in the future. We describe recent experiments using the three global arrays to measure the aberration of gravity when Jupiter and Saturn passed within a few arcmin of bright radio sources. These reductions are still in progress, but the anticipated positional accuracy of the VLBA experiment may be about 0.01 mas.

Millisecond and binary pulsars are the most stable astronomical standards of frequency. They can be applied to solving a number of problems in astronomy and time-keeping metrology including the search for a stochastic gravitational wave background in the early universe, testing general relativity, and establishing a new time-scale. The full exploration of pulsar properties requires that proper unbiased estimates of spin and orbital parameters of the pulsar be obtained. These estimates depend essentially on the random noise components present in pulsar timing residuals. The instrumental white noise has predictable statistical properties and makes no harm for interpretation of timing observations, while the astrophysical/geophysical low-frequency noise corrupts them, thus, reducing the quality of tests of general relativity and decreasing the stability of the pulsar time scale.

This note suggests nomenclature for dealing with the units of various astronomical quantities that are used with the relativistic time scales TT, TDB, TCB and TCG. It is suggested to avoid wordings like "TDB units" and "TT units" and avoid contrasting them to "SI units". The quantities intended for use with TCG, TCB, TT or TDB should be called "TCG-compatible", "TCB-compatible", "TT-compatible" or "TDB-compatible", respectively. The names of the units second and meter for numerical values of all these quantities should be used with out any adjectives. This suggestion comes from a special discussion forum created within IAU Commission 52 "Relativity in Fundamental Astronomy".

We have used the Very Long Baseline Array (VLBA) at 43, 23 and 15 GHz to measure the solar gravitational deflection of radio waves among four radio sources during an 18-day period in October 2005. Using phase-referenced radio interferometry to fit the measured phase delay to the propagation equation of the parameterized post-Newtonian (PPN) formalism, we have determined the deflection parameter gamma = 0.9998 +/- 0.0003$ (68% confidence level), in agreement with General Relativity. The results come mainly from 43 GHz observations where the refraction effects of the solar corona were negligible beyond 3 degrees from the sun. The purpose of this experiment is three-fold: to improve on the previous results in the gravitational bending experiments near the solar limb; to examine and evaluate the accuracy limits of terrestrial VLBI techniques; and to determine the prospects and outcomes of future experiments. Our conclusion is that a series of improved designed experiments with the VLBA could increase the presented accuracy by at least a factor of 4.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is an international project, with major contributions from Europe and China and is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals. A second mission, ASTROD II is envisaged as a three-spacecraft mission which would test General Relativity to one part per billion, enable detection of solar g-modes, measure the solar Lense-Thirring effect to 10 parts per million, and probe gravitational waves at frequencies below the LISA bandwidth. In the third phase (ASTROD III or Super-ASTROD), larger orbits could be implemented to map the outer solar system and to probe primordial gravitational-waves at frequencies below the ASTROD II bandwidth.

Experimental verification of the existence of gravimagnetic fields generated by currents of matter is important for a complete understanding and formulation of gravitational physics. Although the rotational `intrinsic' gravimagnetic field has been extensively studied and is now being measured by the Gravity Probe B, the `extrinsic' gravimagnetic field generated by the translational current of matter is less well studied. The present paper uses the post-Newtonian parametrized Einstein and light geodesics equations to show that the `extrinsic' gravimagnetic field generated by the translational current of matter can be measured by observing the relativistic time delay and/or light deflection caused by the moving mass. We prove that the `extrinsic' gravimagnetic field is generated by the relativistic effect of the aberration of the gravity force caused by the Lorentz transformation of the metric tensor and the Levi-Civita connection. We show that the Lorentz transformation of the gravity field variables is equivalent to the technique of the retarded Lienard-Wiechert gravitational potentials predicting that a light particle is deflected by gravitational field of a moving body from its retarded position so that both general-relativistic phenomena -- the aberration and the retardation of gravity -- are tightly connected and observing the aberration of gravity proves that gravity has a causal nature. We explain in this framework the 2002 deflection experiment of a quasar by Jupiter where the aberration of gravity from its orbital motion was measured with accuracy 20%. We describe a theory of VLBI experiment to measure the gravitational deflection of radio waves from a quasar by the Sun, as viewed by a moving observer from the geocentric frame, to improve the measurement accuracy of the aberration of gravity to a few percent.

General relativistic deflection of light by mass, dipole, and quadrupole moments of gravitational field of a moving massive planet in the Solar system is derived. All terms of order 1 microarcsecond are taken into account, parametrized, and classified in accordance with their physical origin. We calculate the instantaneous patterns of the light-ray deflections caused by the monopole, the dipole and the quadrupole moments, and derive equations describing apparent motion of the deflected position of the star in the sky plane as the impact parameter of the light ray with respect to the planet changes due to its orbital motion. The present paper gives the physical interpretation of the observed light-ray deflections and discusses the observational capabilities of the near-future optical (SIM) and radio (SKA) interferometers for detecting the Doppler modulation of the radial deflection, and the dipolar and quadrupolar light-ray bendings by the Jupiter and the Saturn.

Astrophysical space missions deliver invaluable information about our universe, stellar dynamics of our galaxy, and motion of celestial bodies in the solar system. Astrometric space missions SIM and Gaia will determine distances to stars and cosmological objects as well as their physical characteristics and positions on the celestial sphere with microarcsecond precision. These and other space missions dedicated to exploration of the solar system are invaluable for experimental testing of general relativity. Permanently growing accuracy of space and ground-based astronomical observations require corresponding development of relativistic theory of reference frames, motion of celestial bodies, and propagation of light/radio signals from a source of light/radio to observer. Such theory must be based on Einstein's general relativity and account for various relativistic effects both in the solar system and outside of its boundary. We describe a hierarchy of the relativistic frames adopted by the International Astronomical Union in 2000, and outline directions for its theoretical and practical extentions by matching the IAU 2000 reference frames in the solar system to the cosmological Friedman-Robertson-Walker reference frame and to the frames used in the parametrized post-Newtonian formalism.

We discuss a covariant generalization of the parametrized post-Newtonian (PPN) formalism in a class of scalar-tensor theories of gravity. It includes an exact construction of a set of global and local (Fermi-like) references frames for an isolated N-body astronomical system as well as PPN multipolar decomposition of gravitational field in these frames. We derive PPN equations of translational and rotational motion of extended bodies taking into account all gravitational multipoles and analyze the body finite-size effects in relativistic dynamics that can be important at the latest stages of orbital evolution of coalescing binary systems. We also reconcile the IAU 2000 resolutions on the general relativistic reference frames in the solar system with the PPN equations of motion of the solar system bodies used in JPL ephemerides.

We describe our explicit Lorentz-invariant solution of the Einstein and null geodesic equations for the deflection experiment of 2002 September 8 when a massive moving body, Jupiter, passed within 3.7' of a line-of-sight to a distant quasar. We develop a general relativistic framework which shows that our measurement of the retarded position of a moving light-ray deflecting body (Jupiter) by making use of the gravitational time delay of quasar's radio wave is equivalent to comparison of the relativistic laws of the Lorentz transformation for gravity and light. Because, according to Einstein, the Lorentz transformation of gravity field variables must depend on a fundamental speed $c$, its measurement through the retarded position of Jupiter in the gravitational time delay allows us to study the causal nature of gravity and to set an upper limit on the speed of propagation of gravity in the near zone of the solar system as contrasted to the speed of the radio waves. We discuss the misconceptions which have inhibited the acceptance of this interpretation of the experiment. We also comment on other interpretations of this experiment by Asada, Will, Samuel, Pascual-Sanchez, and Carlip and show that their `speed of light' interpretations confuse the Lorentz transformation for gravity with that for light, and the fundamental speed of gravity with the physical speed of light from the quasar.

Recent papers by Samuel declared that the linearized post-Newtonian v/c effects are too small to have been measured in the recent experiment involving Jupiter and quasar J0842+1845 that was used to measure the ultimate speed of gravity defined as a fundamental constant entering in front of each time derivative of the metric tensor in the Einstein gravity field equations. We describe our Lorentz-invariant formulation of the Jovian deflection experiment and confirm that v/c effects are do observed, as contrasted to the erroneous claim by Samuel, and that they vanish if and only if the speed of gravity is infinite.

A prolonged timing of millisecond pulsars has revealed low-frequency uncorrelated noise, presumably of astrophysical origin, in the pulse arrival time (PAT) residuals for some of them. In most cases, pulsars in globular clusters show a low-frequency modulation of their rotational phase and spin rate. The relativistic time delay of the pulsar signal in the curved space time of randomly distributed and moving globular cluster stars (the Shapiro effect) is suggested as a possible cause of this modulation. Given the smallness of the aberration corrections that arise from the nonstationarity of the gravitational field of the randomly distributed ensemble of stars under consideration, a formula is derived for the Shapiro effect for a pulsar in a globular cluster. The derived formula is used to calculate the autocorrelation function of the low-frequency pulsar noise, the slope of its power spectrum, and the behavior of the $\sigma_z$ statistic that characterizes the spectral properties of this noise in the form of a time function. The Shapiro effect under discussion is shown to manifest itself for large impact parameters as a low-frequency noise of the pulsar spin rate with a spectral index of n=-1.8 that depends weakly on the specific model distribution of stars in the globular cluster. For small impact parameters, the spectral index of the noise is n=-1.5.

One of the main endeavors of fundamental astrometry is to establish a practical realization of an inertial reference frame anchored to celestial objects whose positions are defined in the barycentric coordinates of the solar system. The development of astrometric facilities operating from space at a microarcsecond level of precision makes the non-uniformity of the galactic motion of the barycenter an observable effect that violates the inertiality of the barycentric frame. Most of the observable effect is caused by secular acceleration of the barycenter with respect to the center of the Galaxy. The acceleration results in a pattern of secular aberration which is observable astrometrically as a systematic vector field of the proper motions of distant quasars. We employ the vector spherical harmonics to describe the predicted field of the proper motions and evaluate its amplitude at each point on the celestial sphere. It is shown that the pattern of secular aberration is fully represented by three low-order electric-type vector harmonics, and hence, it is easily distinguishable from the residual rotations of the reference frame and other possible effects, such as the hypothetical long-period gravitational waves. Comprehensive numerical simulations of the grid astrometry with SIM PlanetQuest are conducted. The full covariance matrix of the simulated grid solution is used to evaluate the covariances of the three electric harmonic coefficients, representing the secular aberration pattern of proper motions. We conclude that the grid astrometry with SIM PlanetQuest will be sensitive to the main galactocentric component of secular acceleration, while the peculiar acceleration of the Sun with respect to LSR is expected to be too small to be detected with this astrometric space interferometer.

Recent review article by S. Samuel "On the speed of gravity and the Jupiter/Quasar measurement" published in the International Journal of Modern Physics D13, 1753 (2004) provides the reader with a misleading "theory" of the relativistic time delay in general theory of relativity. Furthermore, it misquotes original publications by Kopeikin and Fomalont & Kopeikin related to the measurement of the speed of gravity by VLBI. We summarize the general relativistic principles of the Lorentz-invariant theory of propagation of light in time-dependent gravitational field, derive Lorentz-invariant expression for the relativistic time delay, and finally explain why Samuel's "theory" is conceptually incorrect and confuses the speed of gravity with the speed of light.

This paper discusses spectral sensitivity of solitary and binary pulsars to a low-frequency (colored) noise. It is a third paper in a series of papers devoted to analysis of the low-frequency noise in pulsar timing observations and its impact on observable fitting parameters of timing model.

Recent measurements of the propagation of the quasar's radio signal past Jupiter are directly sensitive to the time-dependent effect from the geometric sector of general relativity which is proportional to the speed of propagation of gravity but not the speed of light. It provides a first confirmative measurement of the fundamental speed of the Einstein general principle of relativity for gravitational field.

We have determined the relativistic light deflection of the quasar J0842+1835 as Jupiter passed within 3.7' on 2002 September 8, by measuring the time delay using the VLBA and Effelsberg radio telescopes at 8.4 GHz. At closest approach, General Relativity (GR) predicts a radial (static) deflection of 1190 microarcsec, and tangential (retarded) deflection in the direction of Jupiter's motion of 51 microarcsec. Our experiment achieved an rms position error of <10 microarcsec, and measured this retarded deflection to be 0.98 +/- 0.19 (rms error) times that predicted by GR. The increase positional accuracy for this VLBI phase referencing experiment was achieved by using two calibrator sources. Comments on the interpretation of this experiment are given.

Gravitational physics of VLBI experiment conducted on September 8, 2002 and dedicated to measure the speed of gravity (a fundamental constant in the Einstein equations) is treated in the first post-Newtonian approximation. Explicit speed-of-gravity parameterization is introduced to the Einstein equations to single out the retardation effect caused by the finite speed of gravity in the relativistic time delay of light, passing through the variable gravitational field of the solar system. The speed-of-gravity 1.5 post-Newtonian correction to the Shapiro time delay is derived and compared with our previous result obtained by making use of the post-Minkowskian approximation. We confirm that the 1.5 post-Newtonian correction to the Shapiro delay depends on the speed of gravity $c_g$ that is a directly measurable parameter in the VLBI experiment.

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.

Equations of light, propagating from quasar to observer on earth, are integrated in the time-dependent gravitational field of the solar system by making use of either retarded or advanced solutions of the Einstein field equations. This technique allows to separate explicitly the effects associated with the propagation of gravity from those associated with light in the integral expression for the relativistic time delay of light. We prove that the relativistic correction to the Shapiro time delay discovered by Kopeikin (ApJ, 556, L1, 2001) changes sign if one retains direction of the light propagation but replaces retarded with advanced solution of the gravitational field equations. Hence, this correction is due to the propagation of gravity and allows to measure its speed. Relativistic VLBI experiment conducted by Fomalont and Kopeikin in September 2002 proves that the propagation of gravitational field is characterized by the retarded potentials.

We study the gravitational time delay in ray propagation due to rotating masses in the linear approximation of general relativity. Simple expressions are given for the gravitomagnetic time delay that occurs when rays of radiation cross a slowly rotating shell and propagate in the field of a distant rotating source. Moreover, we calculate the local gravitational time delay in the Goedel universe. The observational consequences of these results in the case of weak gravitational lensing are discussed.

We study the effect of retardation of gravity in binary pulsars. It appears in pulsar timing formula as a periodic excess time delay to the Shapiro effect. The retardation of gravity effect can be large enough for observation in binary pulsars with the nearly edgewise orbits and relatively large ratio of the projected semimajor axis to the orbital period of the pulsar. If one succeeds in measuring the retardation of gravity it will give further experimental evidence in favor of General Relativity.

A relativistic sub-picosecond model of gravitational time delay in radio astronomical observations is worked out and a new experimental test of general relativity is discussed in which the effect of retardation of gravity associated with its finite speed can be observed. As a consequence, the speed of gravity can be measured by differential VLBI observations. Retardation in propagation of gravity is a central part of the Einstein theory of general relativity which has not been tested directly so far. The idea of the proposed gravitational experiment is based on the fact that gravity in general relativity propagates with finite speed so that the deflection of light caused by the body must be sensitive to the ratio of the body's velocity to the speed of gravity. The interferometric experiment can be performed, for example, during the very close angular passage of a quasar by Jupiter. Due to the finite speed of gravity and orbital motion of Jupiter, the variation in its gravitational field reaches observer on Earth not instantaneously but at the retarded instant of time and should appear as a velocity-dependent excess time delay in addition to the well-known Shapiro delay, caused by the static part of the Jupiter's gravitational field. Such Jupiter-QSO encounter events take place once in a decade. The next such event will occur on September 8, 2002 when Jupiter will pass by quasar J0842+1835 at the angular distance 3.7 arcminutes. If radio interferometric measurement of the quasar coordinates in the sky are done with the precision of a few picoseconds (about 5 microarcseconds) the effect of retardation of gravity and its speed of propagation may be measured with an accuracy about 10%.

A new gauge-invariant approach for describing cosmological perturbations is developed. It is based on a physically motivated splitting of the stress-energy tensor of the perturbation into two parts - the bare perturbation and the complementary perturbation associated with stresses in the background gravitational field induced by the introduction of the bare perturbation. The complementary perturbation of the stress-energy tensor is explicitly singled out and taken to the left side of the perturbed Einstein equations so that the bare stress-energy tensor is the sole source for the perturbation of the metric tensor and both sides of these equations are gauge invariant with respect to infinitesimal coordinate transformations. For simplicity we analyze the perturbations of the spatially-flat Friedmann-Lemaitre-Robertson-Walker dust model. A cosmological gauge can be chosen such that the equations for the perturbations of the metric tensor are completely decoupled for the h_{00}, h_{0i}, and h_{ij} metric components and explicitly solvable in terms of retarded integrals.

Propagation of light in the gravitational field of self-gravitating spinning bodies moving with arbitrary velocities is discussed. The gravitational field is assumed to be "weak" everywhere. Equations of motion of a light ray are solved in the first post-Minkowskian approximation that is linear with respect to the universal gravitational constant $G$. We do not restrict ourselves with the approximation of gravitational lens so that the solution of light geodesics is applicable for arbitrary locations of source of light and observer. This formalism is applied for studying corrections to the Shapiro time delay in binary pulsars caused by the rotation of pulsar and its companion. We also derive the correction to the light deflection angle caused by rotation of gravitating bodies in the solar system (Sun, planets) or a gravitational lens. The gravitational shift of frequency due to the combined translational and rotational motions of light-ray-deflecting bodies is analyzed as well. We give a general derivation of the formula describing the relativistic rotation of the plane of polarization of electromagnetic waves (Skrotskii effect). This formula is valid for arbitrary translational and rotational motion of gravitating bodies and greatly extends the results of previous researchers. Finally, we discuss the Skrotskii effect for gravitational waves emitted by localized sources such as a binary system. The theoretical results of this paper can be applied for studying various relativistic effects in microarcsecond space astrometry and developing corresponding algorithms for data processing in space astrometric missions such as FAME, SIM, and GAIA.

It is shown that the finite speed of gravity affects very-long baseline interferometric observations of quasars during the time of their line-of-sight close angular encounter with Jupiter. The next such event will take place in 2002, September 8. The present Letter suggests a new experimental test of general relativity in which the effect of propagation of gravity can be directly measured by very-long baseline interferometry as an excess time delay in addition to the logarithmic Shapiro time delay (Shapiro, I. I., 1964, Phys. Rev. Lett., 13, 789).

Recent modern space missions deliver invaluable information about origin of our universe, physical processes in the vicinity of black holes and other exotic astrophysical objects, stellar dynamics of our galaxy, etc. On the other hand, space astrometric missions make it possible to determine with unparalleled precision distances to stars and cosmological objects as well as their physical characteristics and positions on the celestial sphere. Permanently growing accuracy of space astronomical observations and the urgent need for adequate data processing algorithms require corresponding development of an adequate theory of reference frames along with unambiguous description of propagation of light rays from a source of light to observer. Such a theory must be based on the Einstein's general relativity and account for numerous relativistic effects both in the solar system and outside of its boundary. The main features of the relativistic theory of reference frames are presented in this work. A hierarchy of the frames is described starting from the perturbed cosmological Friedmann-Robertson-Walker metric and going to the observer's frame through the intermediate barycentric and geocentric frames in the solar system. Microarcsecond astrometry and effects of propagation of light rays in time-dependent gravitational fields are discussed as well.

The accuracy of astrometric observations conducted via a space-borne optical interferometer orbiting the Earth is expected to approach a few microarcseconds. Data processing of such extremely high-precision measurements requires access to a rigorous relativistic model of light ray propagation developed in the framework of General Relativity. The data-processing of the space interferometric observations must rely upon the theory of general-relativistic transformations between the spacecraft, geocentric, and solar barycentric reference systems allowing unique and unambiguous interpretation of the stellar aberration and parallax effects. On the other hand, the algorithm must also include physically adequate treatment of the relativistic effect of light deflection caused by the spherically-symmetric (monopole-dependent) part of the gravitational field of the Sun and planets as well as the quadrupole- and spin-dependent counterparts of it. In some particular cases the gravitomagnetic field induced by the translational motion of the Sun and planets should be also taken into account for unambigious prediction of the light-ray deflection angle. In the present paper we describe the corresponding software program for taking into account all classical (proper motion, parallax, etc.) and relativistic (aberration, deflection of light) effects up to the microarcsecond threshold and demonstrate, using numerical simulations, how observations of stars and/or quasars conducted on board a space optical interferometer orbiting the Earth can be processed and disentangled.

The Lorentz covariant theory of precise Doppler measurements (PDM) based on the retarded Li\'{e}nard-Wiechert solution of the Einstein equations is described. An exact solution of equations of light propagation in the field of arbitrary moving bodies, which drastically extends the range of applicability of the new theory of PDM, is obtained. An explicit formula for the gravitational shift of frequency is given in analytic form. The limiting cases of the Doppler observations in gravitational lensing and of the spacecraft's Doppler tracking are described in more detail. We also present the post-Newtonian theory of the PDM developed for searching relativistic effects in close optical binaries and massive planetary systems.

The influence of the low-frequency timing noise on the precision of measurements of the Keplerian and post-Keplerian orbital parameters in binary pulsars is studied. Fundamental limits on the accuracy of tests of alternative theories of gravity in the strong-field regime are established. The gravitational low-frequency timing noise formed by an ensemble of binary stars is briefly discussed.

The Lorentz covariant theory of propagation of light in the (weak) gravitational fields of N-body systems consisting of arbitrarily moving point-like bodies with constant masses is constructed. The theory is based on the Lienard-Wiechert presentation of the metric tensor. A new approach for integrating the equations of motion of light particles depending on the retarded time argument is applied. In an approximation which is linear with respect to the universal gravitational constant, G, the equations of light propagation are integrated by quadratures and, moreover, an expression for the tangent vector to the perturbed trajectory of light ray is found in terms of instanteneous functions of the retarded time. General expressions for the relativistic time delay, the angle of light deflection, and gravitational red shift are derived. They generalize previously known results for the case of static or uniformly moving bodies. The most important applications of the theory are given. They include a discussion of the velocity dependent terms in the gravitational lens equation, the Shapiro time delay in binary pulsars, and a precise theoretical formulation of the general relativistic algorithm of data processing of radio and optical astrometric measurements in the non-stationary gravitational field of the solar system. Finally, proposals for future theoretical work being important for astrophysical applications are formulated.

Localized astronomical sources like a double stellar system, rotating neutron star, or a massive black hole at the center of the Milky Way emit periodic gravitational waves. For a long time only a far-zone contribution of gravitational fields of the localized sources (plane-wave-front approximation) were a matter of theoretical analysis. We demonstrate how this analysis can be extended to take into account near-zone and intermediate-zone contributions as well. The formalism is used to calculate gravitational-wave corrections to the Shapiro time delay in binary pulsars and low-frequency (LF) pulsar timing noise produced by an ensemble of double stars in our galaxy.

A consistent approach for an exhaustive solution of the problem of propagation of light rays in the field of gravitational waves emitted by a localized source of gravitational radiation is developed in the first post-Minkowskian and quadrupole approximation of General Relativity. We demonstrate that the equations of light propagation in the retarded gravitational field of an arbitrary localized source emitting quadrupolar gravitational waves can be integrated exactly. The influence of the gravitational field on the light propagation is examined not only in the wave zone but also in cases when light passes through the intermediate and near zones of the source. Explicit analytic expressions for light deflection and integrated time delay (Shapiro effect) are obtained accounting for all possible retardation effects and arbitrary relative locations of the source of gravitational waves, that of light rays, and the observer. It is shown that the ADM and harmonic gauge conditions can both be satisfied simultaneously outside the source of gravitational waves. Their use drastically simplifies the integration of light propagation equations and those for the motion of light source and observer in the field of the source of gravitational waves, leading to the unique interpretation of observable effects. The two limiting cases of small and large values of impact parameter are elaborated in more detail. Explicit expressions for Shapiro effect and deflection angle are obtained in terms of the transverse-traceless part of the space-space components of the metric tensor. We also discuss the relevance of the developed formalism for interpretation of radio interferometric and timing observations, as well as for data processing algorithms for future gravitational wave detectors.

The determination of velocities of stars from precise Doppler measurements is described here using relativistic theory of astronomical reference frames so as to determine the Keplerian and post-Keplerian parameters of binary systems. We apply successive Lorentz transformations and the relativistic equation of light propagation to establish the exact treatment of Doppler effect in binary systems both in special and general relativity theories. As a result, the Doppler shift is a sum of (1) linear in $c^{-1}$ terms, which include the ordinary Doppler effect and its variation due to the secular radial acceleration of the binary with respect to observer; (2) terms proportional to $c^{-2}$, which include the contributions from the quadratic Doppler effect caused by the relative motion of binary star with respect to the Solar system, motion of the particle emitting light and diurnal rotational motion of observer, orbital motion of the star around the binary's barycenter, and orbital motion of the Earth; and (3) terms proportional to $c^{-2}$, which include the contributions from redshifts due to gravitational fields of the star, star's companion, Galaxy, Solar system, and the Earth. After parameterization of the binary's orbit we find that the presence of periodically changing terms in the Doppler schift enables us disentangling different terms and measuring, along with the well known Keplerian parameters of the binary, four additional post-Keplerian parameters, including the inclination angle of the binary's orbit, $i$. We briefly discuss feasibility of practical implementation of these theoretical results, which crucially depends on further progress in the technique of precision Doppler measurements.

Pulsars are the most stable natural frequency standards. They can be applied to a number of principal problems of modern astronomy and time-keeping metrology. The full exploration of pulsar properties requires obtaining unbiased estimates of the spin and orbital parameters. These estimates depend essentially on the random noise component being revealed in the residuals of time of arrivals (TOA). In the present paper, the influence of low-frequency ("red") timing noise with spectral indices from 1 to 6 on TOA residuals, variances, and covariances of estimates of measured parameters of single and binary pulsars are studied. In order to determine their functional dependence on time, an analytic technique of processing of observational data in time domain is developed which takes into account both stationary and non-stationary components of noise. Our analysis includes a simplified timing model of a binary pulsar in a circular orbit and procedure of estimation of pulsar parameters and residuals under the influence of red noise. We reconfirm that uncorrelated white noise of errors of measurements of TOA brings on gradually decreasing residuals, variances and covariances of all parameters. On the other hand, we show that any red noise causes the residuals, variances, and covariances of certain parameters to increase with time. Hence, the low frequency noise corrupts our observations and reduces experimental possibilities for better tests of General Relativity Theory. We also treat in detail the influence of a polynomial drift of noise on the residuals and fitting parameters. Results of the analitic analysis are used for discussion of a statistic describing stabilities of kinematic and dynamic pulsar time scales.

Millisecond and binary pulsars are the most stable natural frequency standards which admits to introduce modified versions of universal and ephemeris time scales based correspondingly on the intrinsic rotation of pulsar and on its orbital motion around barycenter of a binary system. Measured stability of these time scales depends on numerous physical phenomena which affect rotational and orbital motion of the pulsar and observer on the Earth, perturb propagation of electromagnetic pulses from pulsar to the observer and bring about random fluctuations in the rate of atomic clock used as a primary time reference in timing observations. On the long time intervals the main reason for the instability of the pulsar time scales is the presence of correlated, low-frequency timing noise in residuals of times of arrivals (TOA) of pulses from the pulsar which has both astrophysical and geophysical origin. Hence, the timing noise can carry out the important physical information about interstellar medium, interior structure of the pulsar, stochastic gravitational waves coming from the early universe, etc. Each specific type of the low-frequency noise can be described in framework of power law spectrum model. Although the data processing of pulsar timing observations in time domain seems to be the most imformative it is significantly important to know to which spectral bands single and binary pulsars, considered as detectors of the low-frequency noise signal, are the most sensitive. The given problem is examined in the present article.

For radio pulsars in orbit with a compact companion, pulsar timing observations have proved to be a powerful tool for identifying the physical nature of the companion. Unfortunately, perhaps the most intriguing system where such a tool could be used, a pulsar in orbit with a black hole, has yet to be discovered. In this paper we give a detailed investigation of what one can learn about the black hole companion via timing observations of the pulsar. We present an analytic calculation for the propagation delay caused by the frame-dragging effect and show that it has the same functional behavior as the modulation of the observed rotational phase of the pulsar caused by the deflection of the radio signals in the gravitational field of the companion (bending delay). Thus, contrary to statements of other authors, the frame-dragging delay is unlikely to be separately measurable in pulsar binaries where the companion is a stellar mass black hole. We demonstrate, however, that the precession of the binary orbit caused by the relativistic spin-orbit coupling can lead to observable effects which can be used to set a lower limit to the black-hole spin, or possibly allow the determination of its magnitude and orientation. We give parameter estimates for two possible systems, a 10 solar mass black hole in orbit either with a young (~0.1 s) pulsar or with a millisecond pulsar. Finally, we discuss the measurability of the quadrupole moment of the rotating black hole companion which would test the presence of a Kerr black hole. As an interesting side result of our calculations, we can give a further argument why the companion of PSR J0045-7319 cannot be a Kerr black hole.

An improved formula for the timing of binary pulsars that accounts for the relativistic deflection of light in the gravitational field of the pulsar's companion is presented, and the measurability of this effect together with its variance estimates are discussed. The deflection of the pulsar's beam trajectory in the gravitational field of its companion leads to variation in the pulsar's rotational phase. This variation appears as the narrow sharp growth of the magnitude of the post-fit pulsar timing residuals in the vicinity of the moment of the superior conjunction of the pulsar with its companion. In contrast to the relativistic Shapiro effect the amplitude of the effect of gravitational deflection of the pulsar radio beam has two peaks with opposite signs, which become sharper as the inclination $i$ of the pulsar's orbit approaches 90$^{\circ}$. The effect under consideration influences the estimation of parameters of the relativistic Shapiro effect in the binary pulsars with nearly edgewise orbits. Its inclusion in the fitting procedure provide more careful measurement of sine of the orbital inclination $i$ as well as the masses of the pulsar $m_p$ and its companion $m_c$. This permits an improved testing of alternative theories of gravity in the strong field regime. Moreover, the measurement of this effect permits independent geometric constraints on the position of the pulsar rotational axis in space. The effect of the gravitational deflection of light has been numerically investigated for binary pulsars with nearly edgewise orbits. It is shown that the effect is observed in general only when $\cos i\leq 0.003.$ This estimate becomes less restrictive as the pulsar's spin axis approaches the line of the sight.