Recently, the secular pericentre precession was analytically computed to the second post-Newtonian (2PN) order by the present author with the Gauss equations in terms of the osculating Keplerian orbital elements in order to obtain closer contact with the observations in astronomical scenarios of potential interest. A discrepancy with previous results by other authors was found. Moreover, some of such findings by the same authors were deemed as mutually inconsistent. In this paper, it is demonstrated that, in fact, two calculational errors plagued the most recent calculation. They are explicitly disclosed and corrected. As a result, all the examined approaches mutually agree yielding the same analytical expression for the total 2PN pericentre precession once the appropriate conversions from the adopted parameterizations are made.

We perform a full analytical and numerical treatment, to the first post-Newtonian (1pN) order, of the general relativistic long-term spin precession of an orbiting gyroscope due to the mass quadrupole moment $J_2$ of its primary without any restriction on either the gyro's orbital configuration and the orientation in space of the symmetry axis $\boldsymbol{\hat{k}}$ of the central body. We apply our results to the past spaceborne Gravity Probe B (GP-B) mission by finding a secular rate of its spin's declination $\delta$ which may be as large as $\lesssim 30-40\,\mathrm{milliarcseconds\,per\,year\,(\mathrm{mas\,yr}^{-1}})$, depending on the initial orbital phase $f_0$. Both our analytical calculation and our simultaneous integration of the equations for the parallel transport of the spin 4-vector and of the geodesic equations of motion of the gyroscope confirm such a finding. For GP-B, the reported mean error in measuring the spin's declination rate amounts to $\sigma^\mathrm{GP-B}_{\dot\delta}=18.3\,\mathrm{mas\,yr}^{-1}$. We also calculate the general analytical expressions of the gravitomagnetic spin precession induced by the primary's angular momentum $\boldsymbol J$. In view of their generality, our results can be extended also to other astronomical and astrophysical scenarios of interest like, e.g., stars orbiting galactic supermassive black holes, exoplanets close to their parent stars, tight binaries hosting compact stellar corpses.

Recently, the possibility that several starless telluric planets may form around supermassive black holes (SMBHs) and receive an energy input from the hole's accretion disk, which, under certain plausible circumstances, may make them habitable in a terrestrial sense, has gained increasing attention. In particular, an observer on a planet orbiting at distance $r=100$ Schwarzschild radii from a maximally rotating Kerr SMBH with mass $M_\bullet = 1\times 10^8\,M_\odot$ in a plane slightly outside the equator of the latter, would see the gravitationally lensed accretion disk the same size as the Sun as seen from the Earth. Moreover, the accretion rate might be imagined to be set in such a way that the apparent disk's temperature would be identical to that of the solar surface. We demonstrate that the post-Newtonian (pN) de Sitter and Lense--Thirring precessions of the spin axis of such a world would rapidly change, among other things, its tilt, $\varepsilon$, to its orbital plane by tens to hundreds of degrees over a time span of, say, just $\Delta t =400\,\mathrm{yr}$, strongly depending on the obliquity $\eta_\bullet$ of the SMBH's spin to the orbital plane. Thus, such relativistic features would have per se a relevant impact on the long-term habitability of the considered planet. Other scenarios are examined as well.

We comment on a recent study reporting evidence for the general relativistic Lense-Thirring secular precession of the inclination $I$ of the orbital plane to the plane of the sky of the tight binary system PSR J1141-6545 made of a white dwarf and an emitting radiopulsar of comparable masses. The quadrupole mass moment $Q_2^\mathrm{c}$ and the angular momentum ${\boldsymbol S}^\mathrm{c}$ of the white dwarf cause the detectable effects on $I$ with respect to the present-day accuracy in the pulsar's timing. The history-dependent and model-dependent assumptions to be made on $Q_2^\mathrm{c}$ and ${\boldsymbol S}^\mathrm{c}$, required even just to calculate the analytical expressions for the resulting post-Keplerian precessions, may be deemed as too wide in order to claim a successful test of the Einsteinian gravitomagnetic effect. Moreover, depending on how $Q_2^\mathrm{c}$ is calculated, the competing quadrupole-induced rate of change, which is a major source of systematic uncertainty, may be up to $\lesssim 30-50\%$ of the Lense-Thirring effect for most of the allowed values in the 3D parameter space spanned by the white dwarf's spin period $P_\mathrm{s}$, and the polar angles $i_\mathrm{c},\,\zeta_\mathrm{c}$ of its spin axis. The possible use of the longitude of periastron $\dot\varpi$ is investigated as well. It turns out that a measurement of its secular precession, caused, among other things, also by $Q_2^\mathrm{c},\,{\boldsymbol{S}}^\mathrm{c}$, could help in further restricting the permitted regions in the white dwarf's parameter space.

At the second post-Newtonian (2PN) order, the secular pericentre precession $\dot\omega^\mathrm{2PN}$ of either a full two-body system made of well detached non-rotating monopole masses of comparable size and a restricted two-body system composed of a point particle orbiting a fixed central mass have been analytically computed so far with a variety of approaches. We offer our contribution by analytically computing $\dot\omega^\mathrm{2PN}$ in a perturbative way with the method of variation of elliptical elements by explicitly calculating both the direct contribution due to the 2PN acceleration ${\boldsymbol A}^\mathrm{2PN}$, and also an indirect part arising from the self-interaction of the 1PN acceleration ${\boldsymbol A}^\mathrm{1PN}$ in the orbital average accounting for the instantaneous shifts induced by ${\boldsymbol A}^\mathrm{1PN}$ itself. Explicit formulas are straightforwardly obtained for both the point particle and full two-body cases without recurring to simplifying assumptions on the eccentricity $e$. Two different numerical integrations of the equations of motion confirm our analytical results for both the direct and indirect precessions. The values of the resulting effects for Mercury and some binary pulsars are confronted with the present-day level of experimental accuracies in measuring/constraining their pericentre precessions. The supermassive binary black hole in the BL Lac object OJ 287 is considered as well. A comparison with some of the results appeared in the literature is made.

We numerically and analytically work out the first-order post-Newtonian (1pN) orbital effects induced on the semimajor axis $a$, the eccentricity $e$, the inclination $I$, the longitude of the ascending node $\Omega$, the longitude of perihelion $\varpi$, and the mean longitude at epoch $\epsilon$ of a test particle orbiting its primary, assumed static and spherically symmetric, by a distant massive third body X. For Mercury, the rates of change of the linear trends found are $\dot I_\mathrm{1pN}^\mathrm{X} = -4.3\,\mathrm{microarcseconds\,per\,century}\,\left(\mu\mathrm{as\,cty}^{-1}\right)$, $\dot\Omega_\mathrm{1pN}^\mathrm{X} = 18.2\,\mu\mathrm{as\,cty}^{-1}$, $\dot\varpi_\mathrm{1pN}^\mathrm{X} = 30.4\,\mu\mathrm{as\,cty}^{-1}$, $\dot\epsilon_\mathrm{1pN}^\mathrm{X} = 271.4\,\mu\mathrm{as\,cty}^{-1}$, respectively. Such values, which are due to the added actions of the other planets from Venus to Saturn, are essentially at the same level of, or larger by one order of magnitude than, the latest formal errors in the Hermean orbital precessions calculated with the EPM2017 ephemerides. The perihelion precession $\dot\varpi_\mathrm{1pN}^\mathrm{X}$ turns out to be smaller than some values recently appeared in the literature in view of a possible measurement with the ongoing BepiColombo mission. Linear combinations of the supplementary advances of the Keplerian orbital elements for several planets, if determined experimentally by the astronomers, could be set up in order to disentangle the 1pN $N$-body effects of interest from the competing larger precessions like those due to the Sun's quadrupole moment $J_2$ and angular momentum $\boldsymbol{S}$.

We consider a hypothetical planet with the same mass $m$, radius $R$, angular momentum $\boldsymbol{S}$, oblateness $J_2$, semimajor axis $a$, eccentricity $e$, inclination $I$, and obliquity $\varepsilon$ of the Earth orbiting a main-sequence star with the same mass $M_\star$ and radius $R_\star$ of the Sun at a distance $r_\bullet \simeq 1\,\mathrm{parsec}\,\left(\mathrm{pc}\right)$ from a supermassive black hole in the center of the hosting galaxy with the same mass $M_\bullet$ of, say, $\mathrm{M87}^\ast$. We preliminarily investigate some dynamical consequences of its presence in the neighborhood of such a stellar system on the planet's possibility of sustaining complex life over time. In particular, we obtain general analytic expressions for the long-term rates of change, doubly averaged over both the planetary and the galactocentric orbital periods $P_\mathrm{b}$ and $P_\bullet$, of $e,\,I,\,\varepsilon$, which are the main quantities directly linked to the stellar insolation. We find that, for certain orbital configurations, the planet's perihelion distance $q=a\left(1-e\right)$ may greatly shrink and lead to, in some cases, an impact with the star. $I$ may also notably change, with variations even of the order of tens of degrees. On the other hand, $\varepsilon$ does not seem to be particularly affected, being shifted, at most, by $\simeq 0^\circ.02$ over 1 Myr. Our results strongly depend on the eccentricity $e_\bullet$ of the galactocentric motion.

Recently, a discussion about the effects of the anisotropy in the spatial density of Dark Matter in the Solar neighbourhood due to the motion of the Sun through the Galactic halo on the orbital motion of the solar system's planets and their ability to be effectively constrained by the radiotechnical observations collected by the Cassini spacecraft appeared in the literature. We show that the semilatus rectum $p$, the eccentricity $e$, the inclination $I$, the longitude of the ascending node $\Omega$, the longitude of perihelion $\varpi$, and the mean anomaly at epoch $\eta$ of a test particle of a restricted two-body system affected by the gravity of a Dark Matter wake undergo secular rates of change. In the case of Saturn, they are completely negligible, being at the $\simeq 0.1$ millimeter per century and $\simeq 0.05-2$ nanoarcseconds per century level; the current (formal) accuracy level in constraining any anomalous orbital precessions is of the order of $\simeq 0.002-2$ milliarcseconds per century for Saturn. We also numerically simulate the Earth-Saturn range signature $\Delta\rho(t)$ due to the Dark Matter wake over the same time span (2004-2017) covered by the Cassini data record. We find that it is as little as $\simeq 0.1-0.2\,\mathrm{m}$, while the existing range residuals, computed by the astronomers without modeling any Dark Matter wake effect, are at the $\simeq 30\,\mathrm{m}$ level. The local Dark Matter density $\varrho_\mathrm{DM}$ should be larger than the currently accepted value of $\varrho_\mathrm{DM}=0.018\,\mathrm{M}_\odot\,\mathrm{pc}^{-3}$ by a factor of $2.5\times 10^6$ in order to induce a geocentric Kronian range signature so large as to make it discernible in the present-day residuals.

The distinction between the mean anomaly $\mathcal{M}(t)$ and the mean anomaly at epoch $\eta$, and the mean longitude $l(t)$ and the mean longitude at epoch $\epsilon$ is clarified in the context of a possible use of such orbital elements in post-Keplerian tests of gravity, both Newtonian and post-Newtonian. In particular, the perturbations induced on $\mathcal{M}(t),\,\eta,\,l(t),\,\epsilon$ by the post-Newtonian Schwarzschild and Lense-Thirring fields, and the classical accelerations due to the atmospheric drag and the oblateness $J_2$ of the central body are calculated for an arbitrary orbital configuration of the test particle and a general orientation of the primary's spin axis $\boldsymbol{\hat{S}}$. They provide us with further observables which could be fruitfully used, e.g., in better characterizing astrophysical binary systems and in more accurate satellite-based tests around major bodies of the Solar System. Some erroneous and misleading claims by Ciufolini and Pavlis appeared in the literature are confuted. In particular, it is shown that there are no net perturbations of the Lense-Thirring acceleration on either the semimajor axis $a$ and the mean motion $n_\mathrm{b}$. Furthermore, the quadratic signatures on $\mathcal{M}(t)$ and $l(t)$ due to certain disturbing non-gravitational accelerations like the atmospheric drag can be effectively disentangled from the post-Newtonian linear trends of interest provided that a sufficiently extended temporal interval for the data analysis is assumed. A possible use of $\eta$ along with the longitudes of the ascending node $\Omega$ in tests of general relativity with the existing LAGEOS and LAGEOS II satellites is suggested.

We consider a binary system composed of a pulsar and a massive, fast rotating, highly distorted main sequence star as a potential scenario to dynamically put to the test certain post-Keplerian effects of both Newtonian and post-Newtonian nature. We numerically produce time series of the perturbations $\Delta\left(\delta\tau\right)$ of the R{\o}mer-like, orbital component of the pulsar's time delay $\delta\tau$ induced over 10 years by the pN gravitoelectric mass monopole, quadrupole, gravitomagnetic spin dipole and octupole accelerations along with the Newtonian quadrupolar one. We do not deal with the various propagation time delays due to the travelling electromagnetic waves. It turns out that, for a Be-type star with $M = 15\ \textrm{M}_\odot$, $R_\textrm{e} = 5.96\ \textrm{R}_\odot$, $\nu = 0.203$, $S = 3.41\times 10^{45}\ \textrm{J}\ \textrm{s}$, $J_2 = 1.92\times 10^{-3}$ orbited by a pulsar with an orbital period $P_\textrm{b}\simeq 40-70\ \textrm{d}$, the classical oblateness-driven effects are at the $\lesssim 4-150\ \textrm{s}$ level, while the pN shifts are of the order of $\lesssim 1.5-20\ \textrm{s}\ \left(GMc^{-2}\right)$, $\lesssim 10-40\ \textrm{ms}\ \left(GMR^2_\textrm{e} J_2 c^{-2}\right)$, $\lesssim 0.5 - 6\ \textrm{ms}\ \left(GSc^{-2}\right)$, $\lesssim 5 - 20\ \mu\textrm{s}\ \left(GSR^2_\textrm{e} \varepsilon^2 c^{-2}\right)$, depending on their orbital configuration. The root-mean-square (rms) timing residuals $\sigma_{\tau}$ of almost all the existing non-recycled, non-millisecond pulsars orbiting massive, fast rotating main sequence stars are $\lesssim\textrm{ms}$. Thus, such kind of binaries have the potential to become interesting laboratories to measure, or, at least, constrain, some Newtonian and post-Newtonian key features of the distorted gravitational fields of the fast rotating stars hosted by them [Abridged].

HERO (Highly Eccentric Relativity Orbiter) is a space-based mission concept aimed to perform several tests of post-Newtonian gravity around the Earth with a preferably drag-free spacecraft moving along a highly elliptical path fixed in its plane undergoing a relatively fast secular precession. We considered two possible scenarios: a fast, 4-hr orbit with high perigee height of $1,047\,\mathrm{km}$, and a slow, 21-hr path with a low perigee height of $642\,\mathrm{km}$. HERO may detect, for the first time, the post-Newtonian orbital effects induced by the mass quadrupole moment $J_2$ of the Earth which affects the semimajor axis $a$ via a secular trend of $\simeq 4-12\,\mathrm{cm\,yr}^{-1}$, depending on the orbital configuration. Recently, the secular decay of the semimajor axis of the passive satellite LARES was measured with an error as little as $0.7\,\mathrm{cm\,yr}^{-1}$. Also the post-Newtonian spin dipole (Lense-Thirring) and mass monopole (Schwarzschild) effects could be tested to a high accuracy depending on the level of compensation of the non-gravitational perturbations. Moreover, the large eccentricity of the orbit would allow to constrain several long-range modified models of gravity and to accurately measure the gravitational red-shift as well. Each of the six Keplerian orbital elements could be individually monitored to extract the $GJ_2/c^2$ signature, or they could be suitably combined in order to disentangle the post-Newtonian effect(s) of interest from the competing mismodeled Newtonian secular precessions induced by the zonal harmonic multipoles $J_\ell$ of the geopotential. In the latter case, the systematic uncertainty due to the current formal errors $\sigma_{J_\ell}$ of a recent global Earth's gravity field model are better than $1\%$ for all the post-Newtonian effects considered, with a peak of $\simeq 10^{-7}$ for the Schwarzschild-like shifts. [Abridged]

We study a general relativistic gravitomagnetic 3-body effect induced by the spin angular momentum ${\boldsymbol S}_\textrm{X}$ of a rotating mass $M_\textrm{X}$ orbited at distance $r_\textrm{X}$ by a local gravitationally bound restricted two-body system $\mathcal{S}$ of size $r\ll r_\textrm{X}$ consisting of a test particle revolving around a massive body $M$. At the lowest post-Newtonian order, we analytically work out the doubly averaged rates of change of the Keplerian orbital elements of the test particle by finding non-vanishing long-term effects for the inclination $I$, the node $\Omega$ and the pericenter $\omega$. Such theoretical results are confirmed by a numerical integration of the equations of motion for a fictitious 3-body system. We numerically calculate the magnitudes of the post-Newtonian gravitomagnetic 3-body precessions for some astronomical scenarios in our solar system. For putative man-made orbiters of the natural moons Enceladus and Europa in the external fields of Saturn and Jupiter, the relativistic precessions due to the angular momenta of the gaseous giant planets can be as large as $\simeq 10-50~\textrm{milliarcseconds~per~year}~\left(\textrm{mas~yr}^{-1}\right)$. A preliminary numerical simulation shows that, for certain orbital configurations of a hypothetical Europa orbiter, its range-rate signal $\Delta\dot\rho$ can become larger than the current Doppler accuracy of the existing spacecraft Juno at Jupiter, i.e. $\sigma_{\dot\rho}=0.015~\textrm{mm~s}^{-1}$, after 1 d. The effects induced by the Sun's angular momentum on artificial probes of Mercury and the Earth are at the level of $\simeq 1-0.1~\textrm{microarcseconds~per~year}~\left(\mu\textrm{as~yr}^{-1}\right)$.

The inclination $I$ of an Earth's satellite in polar orbit undergoes a secular De Sitter precession of $-7.6$ milliarcseconds per year for a suitable choice of the initial value of its non-circulating node $\Omega$. The competing long-periodic harmonic rates of change of $I$ due to the even and odd zonal harmonics of the geopotential vanish for either a circular or polar orbit, while no secular rates occur at all. This may open up, in principle, the possibility of measuring the geodesic precession in the weak-field limit with an accurately tracked satellite by improving the current bound of $9\times 10^{-4}$ from Lunar Laser Ranging, which, on the other hand, may be even rather optimistic, by one order of magnitude, or, perhaps, even better. The most insidious competing effects are due to the solid and ocean components of the $K_1$ tide since their perturbations have nominal huge amplitudes and the same temporal pattern of the De Sitter signature. They vanish for polar orbits. Departures of $\simeq 10^{-5}-10^{-3}~\textrm{deg}$ from the ideal polar geometry allow to keep the $K_1$ tidal perturbations to a sufficiently small level. Most of the other gravitational and non-gravitational perturbations vanish for the proposed orbital configuration, while the non-vanishing ones either have different temporal signatures with respect to the De Sitter effect or can be modeled with sufficient accuracy. In order to meet the proposed goal, the measurement accuracy of $I$ should be better than $\simeq 35~\textrm{microarcseconds}=0.034~\textrm{milliarcseconds}$ over, say, 5 yr.

In a geocentric kinematically rotating ecliptical coordinate system in geodesic motion through the deformed spacetime of the Sun, both the longitude of the ascending node $\Omega$ and the inclination $I$ of an artificial satellite of the spinning Earth are affected by the post-Newtonian gravitoelectric De Sitter and gravitomagnetic Lense-Thirring effects. By choosing a circular orbit with $I = \Omega = 90\deg$ for a potential new spacecraft, which we propose to name ELXIS, it would be possible to measure each of the gravitomagnetic precessions separately at a percent level, or, perhaps, even better depending on the level of accuracy of the current and future global ocean tide models since the competing classical long-term perturbations on $I,~\Omega$ due to the even and odd zonal harmonics $J_\ell,~\ell=2,~3,~4,\ldots$ of the geopotential vanish. Moreover, a suitable linear combination of $I,~\Omega$ would be able to cancel out the solid and ocean tidal perturbations induced by the $K_1$ tide and, at the same time, enforce the geodetic precessions yielding a secular trend of $-8.3~\textrm{milliarcseconds~per~year}$, thus strengthening the goal of a $\simeq 10^{-5}$ test of the De Sitter effect recently proposed in the literature in the case of an equatorial coordinate system. Relatively mild departures $\Delta I = \Delta\Omega\simeq 0.01-0.1\deg$ from the ideal orbital configuration with $I = \Omega = 90\deg$ are allowed. [Abridged]

Recently, Will calculated an additional contribution to the Mercury's precession of the longitude of perihelion $\varpi$ of the order of $\dot\varpi_\textrm{W}\simeq 0.22$ $\textrm{milliarcseconds per century}$ ($\textrm{mas cty}^{-1}$). It is partly a direct consequence of certain 1pN third-body accelerations entering the planetary equations of motion, and partly an indirect, mixed effect due to the simultaneous interplay of the standard 1pN pointlike acceleration of the primary with the Newtonian $N$-body acceleration, to the quadrupole order, in the analytical calculation of the secular perihelion precession with the Gauss equations. We critically discuss the actual measurability of the mixed effects with respect to direct ones. The current uncertainties in either the magnitude of the Sun's angular momentum $S_\odot$ and the orientation of its spin axis ${\boldsymbol{\hat{S}}}_\odot$ impact the precessions $\dot\varpi_{J_2^\odot},~\dot\varpi_\textrm{LT}$ induced by the Sun's quadrupole mass moment and angular momentum via the Lense-Thirring effect to a level which makes almost impossible to measure $\dot\varpi_\textrm{W}$ even in the hypothesis that it comes entirely from the aforementioned 1pN third-body accelerations. On the other hand, from the point of view of the Lense-Thirring effect itself, the mismodeled quadrupolar precession $\delta\dot\varpi_{J_2^\odot}$ due to the uncertainties in ${\boldsymbol{\hat{S}}}_\odot$ corresponds to a bias of $\simeq 9\%$ of the relativistic one. The resulting simulated mismodeled range and range-rate times series of BepiColombo are at about the per cent level of the nominal gravitomagnetic ones.

Independent tests aiming to constrain the value of the cosmological constant $\Lambda$ are usually difficult because of its extreme smallness $\left(\Lambda \simeq 1\times 10^{-52}~\textrm{m}^{-2},~\textrm{or}~2.89\times 10^{-122}~\textrm{in Planck units}\right)$. Bounds on it from Solar System orbital motions determined with spacecraft tracking are currently at the $\simeq 10^{-43}-10^{-44}~\textrm{m}^{-2}~\left(5-1\times 10^{-113}~\textrm{in Planck units}\right)$ level, but they may turn out to be somewhat optimistic since $\Lambda$ has not yet been explicitly modeled in the planetary data reductions. Accurate $\left(\sigma_{\tau_\textrm{p}}\simeq 1-10~\mu\textrm{s}\right)$ timing of expected pulsars orbiting the Black Hole at the Galactic Center, preferably along highly eccentric and wide orbits, might, at least in principle, improve the planetary constraints by several orders of magnitude. By looking at the average time shift per orbit $\overline{\Delta\delta\tau}^\Lambda_\textrm{p}$, a S2-like orbital configuration with $e=0.8839,~P_\textrm{b}=16~\textrm{yr}$ would allow to obtain preliminarily an upper bound of the order of $\left|\Lambda\right|\lesssim 9\times 10^{-47}~\textrm{m}^{-2}~\left(\lesssim 2\times 10^{-116}~\textrm{in Planck units}\right)$ if only $\sigma_{\tau_\textrm{p}}$ were to be considered. Our results can be easily extended to modified models of gravity using $\Lambda-$type parameters.

The orbital dynamics of a test particle moving in the non-spherically symmetric field of a rotating oblate primary is impacted also by certain indirect, mixed effects arising from the interplay of the different Newtonian and post-Newtonian accelerations which induce known direct perturbations. We systematically calculate the indirect gravitoelectromagnetic shifts per orbit of the Keplerian orbital elements of the test particle arising from the crossing among the first even zonal harmonic $J_2$ of the central body and the post-Newtonian static and stationary components of its gravitational field. We also work out the Newtonian shifts per orbit of order $J_2^2$, and the direct post-Newtonian gravitoelectric effects of order $J_2 c^{-2}$ arising from the equations of motion. In the case of both the indirect and direct gravitoelectric $J_2 c^{-2}$ shifts, our calculation holds for an arbitrary orientation of the symmetry axis of the central body. We yield numerical estimates of their relative magnitudes for systems ranging from Earth artificial satellites to stars orbiting supermassive black holes.

Recently, the existence of a pointlike pertuber PX with $1~m_\textrm{Mars}\lesssim m_\textrm{X}\lesssim 2.4~m_\oplus$ supposedly moving at $65-80~\textrm{AU}$ along a moderately inclined orbit has been hypothesized in order to explain certain features of the midplane of the Kuiper Belt Objects (KBOs). We preliminarily selected two possible scenarios for such a PX, and numerically simulated its effect on the Earth-Saturn range $\rho(t)$ by varying some of its orbital parameters over a certain time span; then, we compared our results with some existing actual range residuals. By assuming $m_\textrm{X} = 1~m_\textrm{Mars}$ and a circular orbit, such a putative new member of our Solar System would nominally perturb $\rho(t)$ by a few km over $\Delta t = 12~\textrm{yr}~(2004-2016)$. However, the Cassini spaceraft accurately measured $\rho(t)$ to the level of $\sigma_\rho\simeq 100~\textrm{m}$. Nonetheless, such a scenario should not be considered as necessarily ruled out since the Cassini data were reduced so far without explicitly modeling any PX. Indeed, a NASA JPL team recently demonstrated that an extra-signature as large as 4 km affecting the Kronian range would be almost completely absorbed in fitting incomplete dynamical models, i.e. without PX itself, to such simulated data, thus not showing up in the standard post-fit range residuals. Larger anomalous signatures would instead occur for $m_\textrm{X} > 1~m_\textrm{Mars}$. Their nominal} amplitude could be as large as $50-150~\textrm{km}$ for $m_\textrm{X} = 2.4~m_\oplus$, thus making less plausible their existence.

Empirically determining the averaged variations of the orbital parameters of the stars orbiting the Supermassive Black Hole (SBH) hosted by the Galactic Centre (GC) in Sgr A* is, in principle, a valuable tool to put on the test the General Theory of Relativity (GTR), in regimes far stronger than those tested so far, and certain key predictions of it like the no-hair theorems. We analytically work out the long-term variations of all the six osculating Keplerian orbital elements of a test particle orbiting a non-spherical, rotating body with quadrupole moment Q_2 and angular momentum S for a generic spatial orientation of its spin axis k. This choice is motivated by the fact that, basically, we do not know the position in the sky of the spin axis of the SBH in Sgr A* with sufficient accuracy. We apply our results to S2, which is the closest star discovered so far having an orbital period P_b = 15.98 yr, and to a hypothetical closer star X with P_b = 0.5 yr. Our calculations are quite general, not being related to any specific parameterization of k, and can be applied also to astrophysical binary systems, stellar planetary systems, and planetary satellite geodesy in which different reference frames, generally not aligned with the primary's rotational axis, are routinely used.

One of the directly measured quantities which are used in monitoring the orbital motions of many of the S stars revolving around the Supermassive Black Hole (SMBH) in the Galactic Center (GC) is their radial velocity (RV) $V$ obtained with near-infrared spectroscopy. Here, we devise a general approach to calculate both the instantaneous variations $\Delta V\left(t\right)$ and the net shifts per revolution $\left\langle\Delta V\right\rangle$ induced on such an observable by some post-Keplerian (pK) accelerations. In particular, we look at the general relativistic Schwarzschild (gravitoelectric) and Lense-Thirring (gravitomagnetic frame-dragging) effects, and the mass quadrupole. It turns out that we may be on the verge of measuring the Schwarzschild-type 1pN static component of the SMBH's field with the S2 star for which RV measurements accurate to about $\simeq 30-50~\textrm{km s}^{-1}$ dating back to $t_0 = 2003.271$ are currently available, and whose orbital period amounts to $P_\textrm{b} = 16$ yr. Indeed, while its expected general relativistic RV net shift per orbit amounts to just $\left\langle\Delta V^\textrm{GE}\right\rangle = -11.6~\textrm{km s}^{-1}$, it should reach a peak value as large as $\Delta V_\textrm{max}^\textrm{GE}\left(t_\textrm{max}\right) = 551~\textrm{km s}^{-1}$ at $t_\textrm{max} = 2018.35$. The periastron shift $\Delta\omega^\textrm{GE}$ of S2 over the same time span will not be larger than $0.2$ deg, while the current accuracy in estimating such an orbital element for this star is of the order of $0.6$ deg. The frame-dragging and quadrupole-induced RV shifts are far smaller for S2, amounting to, at most, $0.19~\textrm{km s}^{-1},0.0039~\textrm{km s}^{-1}$, respectively. Further studies should be dedicated to the impact on the RV of possible diffused mass distribution in the GC and of other individual stars inside and outside the orbit of S2.

We develop a general approach to analytically calculate the perturbations $\Delta\delta\tau_\textrm{p}$ of the orbital component of the change $\delta\tau_\textrm{p}$ of the times of arrival of the pulses emitted by a binary pulsar p induced by the post-Keplerian accelerations due to the mass quadrupole $Q_2$, and the post-Newtonian gravitoelectric (GE) and Lense-Thirring (LT) fields. We apply our results to the so-far still hypothetical scenario involving a pulsar orbiting the Supermassive Black Hole in in the Galactic Center at Sgr A$^\ast$. We also evaluate the gravitomagnetic and quadrupolar Shapiro-like propagation delays $\delta\tau_\textrm{prop}$. By assuming the orbit of the existing S2 main sequence star and a time span as long as its orbital period $P_\textrm{b}$, we obtain $\left|\Delta\delta\tau_\textrm{p}^\textrm{GE}\right|\lesssim 10^3~\textrm{s},~\left|\Delta\delta\tau_\textrm{p}^\textrm{LT}\right|\lesssim 0.6~\textrm{s},\left|\Delta\delta\tau_\textrm{p}^{Q_2}\right|\lesssim 0.04~\textrm{s}$. Faster $\left(P_\textrm{b} = 5~\textrm{yr}\right)$ and more eccentric $\left(e=0.97\right)$ orbits would imply net shifts per revolution as large as $\left|\left\langle\Delta\delta\tau_\textrm{p}^\textrm{GE}\right\rangle\right|\lesssim 10~\textrm{Ms},~\left|\left\langle\Delta\delta\tau_\textrm{p}^\textrm{LT}\right\rangle\right|\lesssim 400~\textrm{s},\left|\left\langle\Delta\delta\tau_\textrm{p}^{Q_2}\right\rangle\right|\lesssim 10^3~\textrm{s}$, depending on the other orbital parameters and the initial epoch. For the propagation delays, we have $\left|\delta\tau_\textrm{prop}^\textrm{LT}\right|\lesssim 0.02~\textrm{s},~\left|\delta\tau_\textrm{prop}^{Q_2}\right|\lesssim 1~\mu\textrm{s}$. The expected precision in pulsar timing in Sgr A$^\ast$ is of the order of $100~\mu\textrm{s}$, or, perhaps, even $1-10~\mu\textrm{s}$.

Detailed numerical analyses of the orbital motion of a test particle around a spinning primary are performed. They aim to investigate the possibility of using the post-Keplerian (pK) corrections to the orbiter's periods (draconitic, anomalistic and sidereal) as a further opportunity to perform new tests of post-Newtonian (pN) gravity. As a specific scenario, the S-stars orbiting the Massive Black Hole (MBH) supposedly lurking in Sgr A$^\ast$ at the center of the Galaxy is adopted. We, first, study the effects of the pK Schwarzchild, Lense-Thirring and quadrupole moment accelerations experienced by a target star for various possible initial orbital configurations. It turns out that the results of the numerical simulations are consistent with the analytical ones in the small eccentricity approximation for which almost all the latter ones were derived. For highly elliptical orbits, the size of all the three pK corrections considered turn out to increase remarkably. The periods of the observed S2 and S0-102 stars as functions of the MBH's spin axis orientation are considered as well. The pK accelerations considered lead to corrections of the orbital periods of the order of 1-100d (Schwarzschild), 0.1-10h (Lense-Thirring) and 1-10^3s (quadrupole) for a target star with a=300-800~AU and e ~ 0.8, which could be possibly measurable by the future facilities.

In the framework of the emergent gravity scenario by Verlinde, it was recently observed by Liu and Prokopec that, among other things, an anomalous pericenter precession would affect the orbital motion of a test particle orbiting an isolated central body. Here, it is shown that, if it were real, its expected magnitude for the inner planets of the Solar System would be at the same level of the present-day accuracy in constraining any possible deviations from their standard perihelion precessions as inferred from long data records spanning about the last century. The most favorable situation for testing the Verlinde-type precession seems to occur for Mars. Indeed, according to recent versions of the EPM and INPOP planetary ephemerides, non-standard perihelion precessions, of whatsoever physical origin, which are larger than some $\approx 0.02-0.11$ milliarcseconds per century are not admissible, while the putative precession predicted by Liu and Prokopec amounts to $0.09$ milliarcseconds per century. Other potentially interesting astronomical and astrophysical scenarios like, e.g., the Earth's LAGEOS II artificial satellite, the double pulsar system PSR J0737-3039A/B and the S-stars orbiting the Supermassive Black Hole in Sgr A$^\ast$ are, instead, not viable because of the excessive smallness of the predicted effects for them.

(Abridged) The trajectory of Saturn is nowadays known at essentially the same accuracy level of the inner planets due to the telemetry of the Cassini spacecraft. Thus, the expected perturbations $\dot\varpi,~\dot\Omega$ due to PX, for which we suggest the name Telisto in view of its remarkable distance, on the Kronian apsidal and draconitic orbital motions are theoretically investigated to tentatively constrain the configuration space of Telisto/Planet Nine itself. To this aim, we compare our predictions $\dot\varpi_\textrm{theo},~\dot\Omega_\textrm{theo}$ to the currently available experimental intervals of values $\Delta\dot\Omega_\textrm{obs},~\Delta\dot\varpi_\textrm{obs}$ determined by astronomers in the recent past without explicitly modeling and solving for PX itself. As such, our results, despite being plausible and in agreement to a large extent with other constraints released in the literature, should be regarded as proof-of-principle investigations aimed to encourage more accurate analyses in future. It turns out that the admissible region in its configuration space is moderately narrow as far as its position along its orbit, reckoned by the true anomaly $f_\textrm{X}$, is concerned, being concentrated around approximately $130~\textrm{deg}\lesssim f_\textrm{X}\lesssim 240~\textrm{deg}$. Telisto/Planet Nine is certainly far from its perihelion ($f_\textrm{X}=0~\textrm{deg}$), in agreement with other recent studies. The future analysis of the data from the ongoing New Horizons mission might be helpful in further constraining the scenario considered here for Telisto/Planet Nine. Its impact on the spaceraft's range over a multi-year span is investigated with a preliminary sensitivity analysis.

The S-stars discovered in the Galactic center (GC) are expected to provide unique dynamical tests of the Kerr metric of the massive black hole (MBH) orbited by them. In order to obtain unbiased measurements of its spin and the related relativistic effects, a comprehensive understanding of the gravitational perturbations of the stars and stellar remnants around the MBH is quite essential. Here, we study the perturbations on the observables of a typical target star, i.e., the apparent orbital motion and the redshift, due to both the spin-induced relativistic effects and the Newtonian attractions of a single or a cluster of disturbing object(s). We find that, in most cases, the Newtonian perturbations on the observables are mainly attributed to the perturbed orbital period of the target star, rather than the Newtonian orbital precessions. The Newtonian perturbations have their unique features when they peak around the pericenter passage in each revolution, which is quite different from those of the spin-induced effects. Looking at the currently detected star S2/S0-2, we find that its spin-induced effects on both the image position and redshift are very likely obscured by the gravitational perturbations from the star S0-102 alone. We also investigate and discuss the Newtonian perturbations on a hypothetical S-star located inside the orbits of the currently detected ones. By considering a number of possible stellar distributions near the central MBH, we find that the spin-induced effects on the apparent position and the redshift dominate over the stellar perturbations for target stars with orbital semimajor axis smaller than $100-400$AU if the MBH is maximally spinning. Our results suggest that, in principle, the stellar perturbations can be removed as they have distinctive morphologies comparing to those of the relativistic Kerr-type signatures.

The orbital motion of a binary system is characterized by various characteristic temporal intervals which, by definition, are different from each other: the draconitic, anomalistic and sidereal periods. They all coincide in the Keplerian case. Such a degeneracy is removed, in general, when a post-Keplerian acceleration is present. We analytically work out the corrections to such otherwise Keplerian periods which are induced by general relativity (Schwarzschild and Lense-Thirring) and, at the Newtonian level, by the quadrupole of the primary. In many astronomical and astrophysical systems, like exoplanets, one of the most accurately determined quantities is just the time span characterizing the orbital revolution, which is often measured independently with different techniques like the transit photometry and the radial velocities. Thus, our results could be useful, in principle, to either constrain the physical properties of the central body and/or perform new tests of general relativity, especially when no other standard observables like, e.g., the orbital precessions are accessible to observations. The difference of two independently measured periods would cancel out the common Keplerian term leaving just a post-Keplerian correction. Furthermore, by comparing the theoretically predicted post-Keplerian expressions $T^{(\textrm{pK})}$ with the experimental accuracy $\sigma_{T_\textrm{exp}}$ in measuring the orbital period(s) it is possible to identify those systems whose observations should be re-processed with genuine post-Keplerian models if $T^{(\textrm{pK})}>\sigma_{T_\textrm{exp}}$. It seems just the case for WASP-33 b since $\sigma_{T_\textrm{exp}}=0.04$ s, while $3$ s$\leq T_\textrm{dra}^{(J_2)}\leq 9.5$ s, $T_\textrm{dra}^{(\textrm{GE)}}=0.36$ s.

The Schwarzschild-de Sitter (SdS) solution exists in the large majority of modified gravity theories, as expected, and in particular the effective cosmological constant is determined by the specific parameters of the given theory. We explore the possibility to use future extended radio-tracking data from the currently ongoing New Horizons mission in the outskirts peripheries of the Solar System, at about 40 au, in order to constrain this effective cosmological constant, and thus to impose constrain on each scenario's parameters. We investigate some of the recently most studied modified gravities, namely $f(R)$ and $f(T)$ theories, dRGT massive gravity, and Ho\v{r}ava-Lifshitz gravity, and we show that New Horizons mission may bring an improvement of one-two orders of magnitude with respect to the present bounds from planetary orbital dynamics.

The radio tracking apparatus of the New Horizons spacecraft, currently traveling to the Pluto system where its arrival is scheduled for July 2015, should be able to reach an accuracy of 10 m (range) and 0.1 mm s^-1 (range-rate) over distances up to 50 au. This should allow to effectively constrain the location of a putative trans-Plutonian massive object, dubbed Planet X (PX) hereafter, whose existence has recently been postulated for a variety of reasons connected with, e.g., the architecture of the Kuiper belt and the cometary flux from the Oort cloud. Traditional scenarios involve a rock-ice planetoid with mX = 0.7mE at some 100 - 200 au, or a Jovian body with mX = 5mJ at about 10,000 - 20,000 au; as a result of our preliminary sensitivity analysis, they should be detectable by New Horizons since they would impact its range at a km level or so over a time span six years long. Conversely, range residuals statistically compatible with zero having an amplitude of 10 m would imply that PX, if it exists, could not be located at less than about 4,500 au (mX = 0.7mE) or 60,000 au (mX = 5mJ), thus making a direct detection quite demanding with the present-day technologies. As a consequence, it would be appropriate to rename such a remote body as Telisto. Also fundamental physics would benefit from this analysis since certain subtle effects predicted by MOND for the deep Newtonian regions of our Solar System are just equivalent to those of a distant pointlike mass.

The predicted Lense-Thirring perihelion precession of Mercury induced by the Sun's angular momentum through its general relativistic gravitomagnetic field amounts to 2 milliarcseconds per century. It turned out to be compatible with the latest experimental determinations of the supplementary perihelion precession of Mercury with the INPOP15a ephemerides, whose accuracy level has nowadays reached the magnitude of the predicted relativistic effect itself thanks to the analysis of some years of tracking data of the MESSENGER spacecraft, which orbited Mercury from 2011 to 2015. A dedicated analysis of three years of MESSENGER data with the DE ephemerides allowed for a $25\%$ determination of the Sun's angular momentum by means of the Lense-Thirring effect, which turned out to be highly correlated with the signature due to the Solar quadrupole mass moment $J_2^{\odot}$.

A sinusoidally time-varying pattern of the values of the Newton's constant of gravitation $G$ measured in Earth-based laboratories over the latest decades has been recently reported in the literature. We put to the test the hypothesis that the aforementioned harmonic variation may pertain $G$ itself in a direct and independent way. We numerically integrated the ad-hoc modified equations of motion of the major bodies of the Solar System by finding that the orbits of the planets would be altered by an unacceptably larger amount in view of the present-day high accuracy astrometric measurements. In the case of Saturn, its geocentric right ascension $\alpha$, declination $\delta$ and range $\rho$ would be affected up to $10^4-10^5$ milliarcseconds and $10^5$ km, respectively; the present-day residuals of such observables are as little as about $4$ milliarcseconds and $10^{-1}$ km, respectively. \textcolor{black}{We analytically calculated the long-term orbital effects induced by the putative harmonic variation of $G$ at hand finding non-zero rates of change for the semimajor axis $a$, the eccentricity $e$ and the argument of pericenter $\omega$ of a test particle. For the LAGEOS satellite, an orbital increase as large as $3.9$ m yr$^{-1}$ is predicted, in contrast with the observed decay of $-0.203\pm 0.035$ m yr$^{-1}$. An anomalous perihelion precession as large as 14 arcseconds per century is implied for Saturn, while latest observations constrain it down to the $10^{-4}$ arcseconds per century level. The rejection level provided by Mercury is of the same order of magnitude.

The laser-tracked geodetic satellites LAGEOS, LAGEOS II and LARES are currently employed, among other things, to measure the general relativistic Lense-Thirring effect in the gravitomagnetic field of the spinning Earth with the hope of providing a more accurate test of such a prediction of the Einstein's theory of gravitation than the existing ones. The secular decay $\dot a$ of the semimajor axes $a$ of such spacecrafts, recently measured in an independent way to a $\sigma_{\dot a}\approx 0.1-0.01$ m yr$^{-1}$ accuracy level, may indirectly impact the proposed relativistic experiment through its connection with the classical orbital precessions induced by the Earth's oblateness $J_2$. \textcolor{black}{Indeed,} the systematic bias due to the current measurement errors $\sigma_{\dot a}$ is of the same order of magnitude of, or even larger than, the expected relativistic signal itself; moreover, it grows linearly with the time span $T$ of the analysis. \textcolor{black}{Therefore, the parameter-fitting algorithms must be properly updated in order to suitably cope with such a new source of systematic uncertainty. Otherwise,} an improvement of one-two orders of magnitude in measuring the orbital decay of the satellites of the LAGEOS family would be required to reduce this source of systematic uncertainty to a percent fraction of the Lense-Thirring signature.

We study the motion of test particles in the metric of a localized and slowly rotating astronomical source, within the framework of linear gravitoelectromagnetism, grounded on a Post-Minkowskian approximation of general relativity. Special attention is paid to gravitational inductive effects due to time-varying gravitomagnetic fields. We show that, within the limits of the approximation mentioned above, there are cumulative effects on the orbit of the particles either for planetary sources or for binary systems. They turn out to be negligible.

By using the most recently published Doppler tomography measurements and accurate theoretical modeling of the oblateness-driven orbital precessions, we tightly constrain some of the physical and orbital parameters of the planetary system hosted by the fast rotating star WASP-33. In particular, the measurements of the orbital inclination $i_{\rm p}$ to the plane of the sky and of the sky-projected spin-orbit misalignment $\lambda$ at two epochs about six years apart allowed for the determination of the longitude of the ascending node $\Omega$ and of the orbital inclination $I$ to the apparent equatorial plane at the same epochs. As a consequence, average rates of change $\dot\Omega_{\rm exp},~\dot I_{\rm exp}$ of this two orbital elements, accurate to a $\approx 10^{-2}~{\rm deg}~{\rm yr}^{-1}$ level, were calculated as well. By comparing them to general theoretical expressions $\dot\Omega_{J_2},~\dot I_{J_2}$ for their precessions induced by an oblate star whose symmetry axis is arbitrarily oriented, we were able to determine the angle $i^{\star}$ between the line of sight the star's spin $S^{\star}$ and its first even zonal harmonic $J_2^{\star}$ obtaining $i^{\star} = 142^{+10}_{-11}~{\rm deg},~J_2^{\star} = (2.1^{+0.8}_{-0.5})\times 10^{-4}.$ As a by-product, the angle between $S^{\star}$ and the orbital angular momentum $L$ is as large as about $\psi \approx 100$ deg $(\psi^{2008} = 99^{+5}_{-4}~{\rm deg},~\psi^{2014} = 103^{+5}_{-4}~{\rm deg})$, and changes at a rate $\dot\psi = 0.7^{+1.5}_{-1.6}~{\rm deg}~{\rm yr}^{-1}$. The predicted general relativistic Lense-Thirring precessions, or the order of $\approx 10^{-3}~{\rm deg}~{\rm yr}^{-1}$, are, at present, about one order of magnitude below the measurability threshold.

In the framework of $f(T)$ theories of gravity, we solve the field equations for $f(T)=T+\alpha T^{n}$, in the weak-field approximation and for spherical symmetry spacetime. Since $f(T)=T$ corresponds to Teleparallel Gravity, which is equivalent to General Relativity, the non linearity of the Lagrangian are expected to produce perturbations of the general relativistic solutions, parameterized by $\alpha$. Hence, we use the $f(T)$ solutions to model the gravitational field of the Sun, and exploit data from accurate tracking of spacecrafts orbiting Mercury and Saturn to infer preliminary insights on what could be obtained about the model parameter $\alpha$ and the cosmological constant $\Lambda$. It turns out that improvements of about one-three orders with respect to the present-day constraints in the literature of magnitude seem possible.

The present Editorial introduces the Special Issue dedicated by the journal Universe to the General Theory of Relativity, the beautiful theory of gravitation of Einstein, a century after its birth. It reviews some of its key features in a historical perspective, and, in welcoming distinguished researchers from all over the world to contribute it, some of the main topics at the forefront of the current research are outlined.

Mindful of the anomalous perihelion precession of Mercury discovered by U. Le Verrier in the second half of the nineteenth century and its successful explanation by A. Einstein with his General Theory of Relativity in the early years of the twentieth century, discrepancies among observed effects in our Solar system and their theoretical predictions on the basis of the currently accepted laws of gravitation applied to known bodies have the potential of paving the way for remarkable advances in fundamental physics. This is particularly important now more than ever, given that most of the Universe seems to be made of unknown substances dubbed Dark Matter and Dark Energy. Should this not be directly the case, Solar system's anomalies could anyhow lead to advancements in cumulative science, as shown to us by the discovery of Neptune in the first half of the nineteenth century. Moreover, investigations in one of such directions can serendipitously enrich the other one as well. The current status of some alleged gravitational anomalies in the Solar system is critically reviewed. They are: a) Possible anomalous advances of planetary perihelia; b) Unexplained orbital residuals of a recently discovered moon of Uranus (Mab); c) The lingering unexplained secular increase of the eccentricity of the orbit of the Moon; d) The so-called Faint Young Sun Paradox; e) The secular decrease of the mass parameter of the Sun; f) The Flyby Anomaly; g) The Pioneer Anomaly; and h) The anomalous secular increase of the astronomical unit

We analytically calculate some orbital effects induced by the Lorentz-invariance/momentum-conservation PPN parameter $\alpha_3$ in a gravitationally bound binary system made of a compact primary orbited by a test particle. We neither restrict ourselves to any particular orbital configuration nor to specific orientations of the primary's spin axis $\boldsymbol{\hat{\psi}}$. We use our results to put preliminary upper bounds on $\alpha_3$ in the weak-field regime by using the latest data from Solar System's planetary dynamics. By linearly combining the supplementary perihelion precessions $\Delta\dot\varpi$ of the Earth, Mars and Saturn, determined with the EPM2011 ephemerides, we infer $|\alpha_3|\lesssim 6\times 10^{-10}$. Our result is about 3 orders of magnitude better than the previous weak-field constraints existing in the literature, and of the same order of magnitude of the bound expected from the future BepiColombo mission to Mercury. It is, by construction, independent of the other preferred-frame PPN parameters $\alpha_1,\alpha_2$, both preliminarily constrained down to a $\approx 10^{-6}$ level. The wide pulsar-white dwarf binary PSR J0407+1607 yields a preliminary upper bound on the strong-field version $\hat{\alpha}_3$ of the Lorentz-invariance/momentum-conservation PPN parameter of the order of $3\times 10^{-17}$. It relies upon certain assumptions on the unknown values of the pulsar's spin axis orientation $\boldsymbol{\hat{\psi}}$, the orbital node $\Omega$ and the inclination $I$. Neither the pulsar's proper motion, still undetected, nor a possible value of the pulsar's mass $m_{\rm p}$ up to two solar masses substantially affect our result.

We analytically compute the long-term orbital variations of a test particle orbiting a central body acted upon by an incident monochromatic plane gravitational wave. We assume that the characteristic size of the perturbed two-body system is much smaller than the wavelength of the wave. Moreover, we also suppose that the wave's frequency is much smaller than the particle's orbital one. We make neither a priori assumptions about the direction of the wavevector nor on the orbital geometry of the planet. We find that, while the semi-major axis is left unaffected, the eccentricity, the inclination, the longitude of the ascending node, the longitude of pericenter and the mean anomaly undergo non-vanishing long-term changes. They are not secular trends because of the slow modulation introduced by the tidal matrix coefficients and by the orbital elements themselves. They could be useful to indepenedently constrain the ultra-low frequency waves which may have been indirectly detected in the BICEP2 experiment. Our calculation holds, in general, for any gravitationally bound two-body system whose characteristic frequency is much larger than the frequency of the external wave. It is also valid for a generic perturbation of tidal type with constant coefficients over timescales of the order of the orbital period of the perturbed particle.

The direct long-term changes occurring in the orbital dynamics of a local gravitationally bound binary system $S$ due to the post-Newtonian tidal acceleration caused by an external massive source are investigated. A class of systems made of a test particle $m$ rapidly orbiting with orbital frequency $n_{\rm b}$ an astronomical body of mass $M$ which, in turn, slowly revolves around a distant object of mass $M^{'}$ with orbital frequency $n_{\rm b}^{'}\ll n_{\rm b}$ is considered. The characteristic frequencies of the non-Keplerian orbital variations of $m$ and of $M$ itself are assumed to be negligible with respect to both $n_{\rm b}$ and $n_{\rm b}^{'}$. General expressions for the resulting Newtonian and post-Newtonian tidal orbital shifts of $m$ are obtained. The future missions BepiColombo and JUICE to Mercury and Ganymede, respectively, are considered in view of a possible detection. The largest effects, of the order of $\approx 0.1-0.5$ milliarcseconds per year (mas yr$^{-1}$), occur for the Ganymede orbiter of the JUICE mission. Although future improvements in spacecraft tracking and orbit determination might, perhaps, reach the required sensitivity, the systematic bias represented by the other known orbital perturbations of both Newtonian and post-Newtonian origin would be overwhelming. The realization of a dedicated artificial mini-planetary system to be carried onboard and Earth-orbiting spacecraft is considered as well. Post-Newtonian tidal precessions as large as $\approx 1-10^2$ mas yr$^{-1}$ could be obtained, but the quite larger Newtonian tidal effects would be a major source of systematic bias because of the present-day percent uncertainty in the product of the Earth's mass times the Newtonian gravitational parameter.

It is known that, in the weak-field and slow-motion approximation of general relativity, the stationary component of the gravitational field of a rotating body, proportional to its angular momentum, is able to discriminate between the opposite directions of motion of a pair of counter-orbiting tests particles moving along geometrically identical paths: it is the so-called gravitomagnetic clock effect. In this paper, we show that, to the first post-Newtonian level, there is also a gravitolectric clock effect induced by the static component of the gravitational field of a non-rotating body. Indeed, the difference of the draconitic periods of two test particles moving in opposite directions along identical quasi-circular orbits turns out to be non-vanishing. It depends on the shape, the size and the pericenter location of the orbit, but it is independent of the initial position along it. On the contrary, the gravitoelectric corrections to the anomalistic periods turn out to be the same for both prograde and retrograde orbital motions. For a hypothetical pair of Earth satellites in Supertundra-like orbits, the gravitoelectric draconitic clock effect could be as large as 46 $\mu$s, while a Juno-like orbital configuration around Jupiter would allow to reach about $0.01$ s.

The recent discovery of the Sedna-like dwarf planet 2012 VP$_{\rm 113}$ by Trujillo and Sheppard has revamped the old-fashioned hypothesis that a still unseen trans-Plutonian object of planetary size, variously dubbed over the years as Planet X, Tyche, Thelisto, may lurk in the distant peripheries of the Solar System. This time, the presence of a super-Earth with mass $m_{\rm X} = 2-15m_{\oplus}$ at a distance $d_{\rm X}\approx 200-300$ astronomical units (AU) was proposed to explain the observed clustering of the arguments of perihelion $\omega$ near $\omega \approx 0^{\circ}$ but not $\omega\approx 180^{\circ}$ for Sedna, 2012 VP$_{\rm 113}$ and other minor bodies of the Solar System with perihelion distances $q>30$ AU and semimajor axes $a>150$ AU. Actually, such a scenario is strongly disfavored by the latest constraints $\Delta\dot\varpi$ on the anomalous perihelion precessions of some Solar System's planets obtained with the INPOP and EPM ephemerides. Indeed, they yield $d_{\rm X}\gtrsim 496-570$ AU ($m_{\rm X}=2m_{\oplus}$), and $d_{\rm X}\gtrsim 970-1111$ AU ($m_{\rm X} = 15 m_{\oplus}$). Much tighter constraints could be obtained in the near future from the New Horizons mission to Pluto.

An empirical formula recently appeared in the literature to explain the observed anomalies of about $\Delta\dot\rho\approx 1-10$ mm s$^{-1}$ in the geocentric range-rates $\dot\rho$ of the Galileo, NEAR and Rosetta spacecraft at some of their past perigee passages along unbound, hyperbolic trajectories. It predicts an anomaly of the order of $6$ mm s$^{-1}$ for the recent flyby of Juno, occurred on 9 October 2013. Data analyses to confirm or disproof it are currently ongoing. We numerically calculate the impact on the geocentric Juno's range rate of some classical and general relativistic dynamical effects which are either unmodelled or mismodelled to a certain level in the software used to process the data. They are: a) The first even zonal harmonic coefficient $J_2$ of the multipolar expansion of the terrestrial gravitational potential causing orbital perturbations both at the $\left.{\rm a}^{'}\right)$ Newtonian ($J_2$) and at the $\left. {\rm a}^{''}\right)$ first post-Newtonian level ($J_2 c^{-2}$) b) The post-Newtonian gravitoelectric (GE) Schwarschild-like component of the Earth's gravitational field c) The post-Newtonian gravitomagnetic (GM) Lense-Thirring effect. The magnitudes of their mismodeled and nominal range-rate signatures are: $\left. {\rm a}^{'}\right)$ $\Delta\dot\rho_{\sigma_{J_2}} \approx 1$ $\mu$m s$^{-1}$ $\left. {\rm a}^{''}\right)$ $\Delta\dot\rho_{J_2 c^{-2}} \approx 0.015$ $\mu$m s$^{-1}$ b) $\Delta\dot\rho_{\rm GE} \approx 25$ $\mu$m s$^{-1}$ c) $\Delta\dot\rho_{\rm GM} \approx 0.05$ $\mu$m s$^{-1}$. If a flyby anomaly as large as a few mm s$^{-1}$ will be finally found also for Juno, it will not be due to any of these standard gravitational effects. (Abridged)

A new analysis of extended data records collected with the Lunar Laser Ranging (LLR) technique performed with improved tidal models was not able to resolve the issue of the anomalous rate $\dot e$ of the eccentricity $e$ of the orbit of the Moon, which is still in place with a magnitude of $\dot e=(5\pm 2)\times 10^{-12}$ yr$^{-1}$. Some possible cosmological explanations are offered in terms of the post-Newtonian effects of the cosmological expansion, and of the slow temporal variation of the relative acceleration rate $\ddot{S} S^{-1}$ of the cosmic scale factor $S$. None of them is successful since their predicted secular rates of the lunar eccentricity are too small by several orders of magnitude.

We consider a non-rotating, massive test particle acted upon by a "pressure"-type, non-geodesic acceleration arising from a certain general class of gravitational theories with nonminimal coupling between the matter and the metric. The resulting orbital perturbations for a two-body system are investigated both analytically and numerically. Among the other long-term effects, a secular increase of the two-body relative distance occurs. In principle, it may yield a physical mechanism for the steady recession of the Earth from the Sun recently proposed to explain the Faint Young Sun Paradox in the Archean eon. At present, the theorists have not yet derived explicit expressions for some of the key parameters of the model, such as the integrated "charge" $\xi$, depending on the matter distribution of the system, and the 4-vector $K^{\mu}=\{K^0,\boldsymbol{K}\}$ connected with the nonminimal function $F$. Thus, we phenomenologically treat them as free parameters, and preliminarily infer some indications on their admissible values according to the most recent Solar System's planetary ephemerides. From the latest determinations of the corrections $\Delta\dot\varpi$ to the standard perihelion precessions, estimated by the astronomers who produced the EPM2011 ephemerides without modeling the theory considered here, we preliminarily obtain $|\xi K|\lesssim 0.1$ kg s$^{-1}$ for Mars. From guesses on what could be the current bounds on the secular rates of change of the planetary semimajor axes, we get $|\xi K^0|\lesssim 1249$ kg s$^{-1}$ for Mars. More effective constraints could be posed by reprocessing the same planetary data sets with dedicated dynamical models including the effects studied here, and explicitly estimating the associated parameters. COBE and GP-B terrestrial satellites yield $|\xi K|\lesssim 2\times 10^{-4}$ kg s$^{-1}$ and $|\xi K_0|\lesssim 2\times 10^{-10}$ kg s$^{-1}$.

Analytical expressions for the orbital precessions affecting the relative motion of the components of a local binary system induced by Lorentz-violating Preferred Frame Effects (PFE) are explicitly computed in terms of the PPN parameters $\alpha_1$, $\alpha_2$. A linear combination of the supplementary perihelion precessions of all the inner planets of the Solar System, able to remove the a-priori bias of unmodelled/mismodelled standard effects such as the general relativistic Lense-Thirring precessions and the classical rates due to the Sun's oblateness $J_2$, allows to infer $|\alpha_1| \leq 6\times 10^{-6}, |\alpha_2| \leq 3.5\times 10^{-5}$. Such bounds should be improved in the near future after processing the data that are being collected by the MESSENGER spacecraft, currently orbiting Mercury. Further improvements may come in the mid-future from the approved BepiColombo mission to Mercury. The constraint $|\alpha_2|\leq 10^{-7}$ existing in the literature is critically discussed (Abridged).

By phenomenologically assuming a slow temporal variation of the percent acceleration rate $\ddot S S^{-1}$ of the cosmic scale factor $S(t)$, it is shown that the orbit of a local binary undergoes a secular expansion. To first order in the power expansion of $\ddot S S^{-1}$ around the present epoch $t_0$, a non-vanishing shift per orbit $\left\langle\Delta r\right\rangle$ of the two-body relative distance $r$ occurs for eccentric trajectories. A general relativistic expression, which turns out to be cubic in the Hubble parameter $H_0$ at the present epoch, is explicitly calculated for it in the case of matter-dominated epochs with Dark Energy. For a highly eccentric Oort comet orbit with period $P_{\rm b}\approx 31$ Myr, the general relativistic distance shift per orbit turns out to be of the order of $\left\langle\Delta r\right\rangle\approx 70$ km. For the Large Magellanic Cloud, assumed on a bound elliptic orbit around the Milky Way, the shift per orbit is of the order of $\left\langle\Delta r\right\rangle\approx 2-4$ pc. Our result has a general validity since it holds in any cosmological model admitting the Hubble law and a slowly varying $\ddot S S^{-1}(t)$. More generally, it is valid for an arbitrary Hooke-like extra-acceleration whose \virg{elastic} parameter $\mathcal{K}$ is slowly time-dependent, irrespectively of the physical mechanism which may lead to it. The coefficient $\mathcal{K}_1$ of the first-order term of the power expansion of $\mathcal{K}(t)$ can be preliminarily constrained in a model-independent way down to a $\mathcal{K}_1\lesssim 2\times 10^{-13}$ yr$^{-3}$ level from latest Solar System's planetary observations. The radial velocities of the double lined spectroscopic binary $\alpha$ Cen AB yield $\mathcal{K}_1\lesssim 10^{-8}$ yr$^{-3}$.

It has been recently proposed to combine the node drifts of the future constellation of 27 Galileo spacecraft together with those of the existing LAGEOS-type satellites to improve the accuracy of the past and ongoing tests of the Lense-Thirring (LT) effect by removing the bias of a larger number of even zonal harmonics $J_{\ell}$ than either done or planned so far. Actually, it seems a difficult goal to be realistically achieved for a number of reasons. First, the LT range signature of a Galileo-type satellite is as little as $0.5$ mm over 3-days arcs, corresponding to a node rate of just $\dot\Omega_{\rm LT}=2$ milliarcseconds per year (mas yr$^{-1}$). Some tesseral and sectorial ocean tides such as $K_1,K_2$ induce long-period harmonic node perturbations with frequencies which are integer multiples of the extremely slow Galileo's node rate $\dot\Omega$ completing a full cycle in about 40 yr. Thus, over time spans $T$ of some years they would act as superimposed semi-secular aliasing trends. Since the coefficients of the $J_{\ell}$-free multisatellite linear combinations are determined only by the semimajor axis $a$, the eccentricity $e$ and the inclination $I$, which are nominally equal for all the Galileo satellites, it is not possible to include all of them. Even by using only one Galileo spacecraft together with the LAGEOS family would be unfeasible because the resulting Galileo coefficient would be $\gtrsim 1$, thus enhancing the aliasing impact of the uncancelled non-conservative and tidal perturbations.

We comment on some statements in a recent paper by Paolozzi and Ciufolini concerning certain remarks raised by us on the realistic accuracy obtainable in testing the general relativistic Lense-Thirring effect in the gravitational field of the Earth with the newly launched LARES satellite together with the LAGEOS and LAGEOS II spacecraft in orbit for a long time. The orbital configuration of LARES is different from that of the originally proposed LAGEOS-3. Indeed, while the latter one should have been launched to the same altitude of LAGEOS (i.e. about $h_{\rm L}=5890$ km) in an orbital plane displaced by $180$ deg with respect to that of LAGEOS ($I_{\rm L}=110$ deg, $I_{\rm L3}=70$ deg), LARES currently moves at a much smaller altitude (about $h_{\rm LR}=1440$ km) and at a slightly different inclination ($I_{\rm LR} = 69.5$ deg). As independently pointed out in the literature by different authors, the overall accuracy of a LARES-LAGEOS-LAGEOS II Lense-Thirring test may be unfavorably \textcolor{black}{impacted} by the lower altitude of LARES with respect to the expected $\approx 1%$ level claimed by Ciufolini \textit{et al.} because of an enhanced sensitivity to the low-degree even zonal geopotential coefficients inducing orbital precessions competing with the relativistic ones. Concerning the previous tests performed with the combined nodes of only LAGEOS and LAGEOS II, an independent analysis recently appeared in the literature indirectly confirms that the total uncertainty in them is likely far from being as little as $10%$.

In the framework of the Parameterized Post-Newtonian (PPN) formalism, we calculate the long-term Preferred Location (PL) effects, proportional to the Whitehead parameter $\xi$, affecting all the Keplerian orbital elements of a localized two-body system, apart from the semimajor axis $a$. They violate the gravitational Local Position Invariance (LPI), fulfilled by General Relativity (GR). We obtain preliminary bounds on $\xi$ by using the latest results in the field of the Solar System planetary ephemerides. The non-detection of any anomalous perihelion precession for Mars allows us to indirectly infer $|\xi|\leq 5.8\times 10^{-6}$. \textcolor{black}{Such a bound is close to the constraint, of the order of $10^{-6}$, expected from the future BepiColombo mission to Mercury. As a complementary approach, the PL effects should be explicitly included in the dynamical models fitted to planetary data sets to estimate $\xi$ in a least-square fashion in a dedicated ephemerides orbit solution.} The ratio of the anomalous perihelion precessions for Venus and Jupiter, determined with the EPM2011 ephemerides at the $<3\sigma$ level, if confirmed as genuine physical effects needing explanation by future studies, rules out the hypothesis $\xi \neq 0$. A critical discussion of the $|\xi| \lesssim 10^{-6}-10^{-7}$ upper bounds obtained in the literature from the close alignment of the Sun's spin axis and the total angular momentum of the Solar System is presented.

Given a solar luminosity L_Ar = 0.75 L_0 at the beginning of the Archean 3.8 Gyr ago, where L_0 is the present-day one, if the heliocentric distance r of the Earth was r_Ar = 0.956 r_0, the solar irradiance would have been as large as I_Ar = 0.82 I_0. It would allowed for a liquid ocean on the terrestrial surface which, otherwise, would have been frozen, contrary to the empirical evidence. By further assuming that some physical mechanism subsequently displaced the Earth towards its current distance in such a way that the irradiance stayed substantially constant over the entire Archean from 3.8 Gyr to 2.5 Gyr ago, a relative recession rate as large as \dot r/r \simeq 3.4 x 10^-11 yr^-1 would have been required. Although such a figure is roughly of the same order of magnitude of the value of the Hubble parameter 3.8 Gyr ago H_Ar = 1.192 H_0 = 8.2 x 10^-11 yr^-1, standard general relativity rules out cosmological explanations for the hypothesized Earth' s recession rate. Instead, a class of modified theories of gravitation with nonminimal coupling between the matter and the metric naturally predicts a secular variation of the relative distance of a localized two-body system, thus yielding a potentially viable candidate to explain the putative recession of the Earth' s orbit. Another competing mechanism of classical origin which could, in principle, allow for the desired effect is the mass loss which either the Sun or the Earth itself may have experienced during the Archean. On the one hand, this implies that our planet should have lost 2% of its present mass in the form of eroded/evaporated hydrosphere which, thus, should have been two orders of magnitude larger than now. On the other hand, it is widely believed that the Sun could have lost mass at an enhanced rate due to a stronger solar wind in the past for not more than \sim 0.2-0.3 Gyr.

In this paper we directly constrain possible spatial variations of the Newtonian gravitational constant G over ranges 0.01-5 AU in various extrasolar multi-planet systems. By means of the third Kepler's law we determine the quantity \Gamma_ XY=G_X/G_Y for each couple of planets X and Y located at different distances from their parent star: deviations of the measured values of \Gamma from unity would signal variations of G. The obtained results for \eta=1-\Gamma are found to be well compatible with zero within the experimental errors (\eta/\delta\eta= 0.2-0.3). A comparison with an analogous test previously performed in our Solar System is made.

The expansion in multipoles of the gravitational potential of a rotating body affects the orbital motion of a test particle orbiting it with long-term perturbations both at a classical and at a relativistic level. In this preliminary sensitivity analysis, we show that, for the first time, the J2 c^-2 effects could be measured by the ongoing Juno mission in the gravitational field of Jupiter during its yearlong science phase (10 November 2016-5 October 2017) thanks to its high eccentricity (e=0.947) and to the huge oblateness of Jupiter (J2=1.47 10^-2). The semi-major axis a and the perijove \omega\ of Juno are expected to be shifted by \Delta a =700-900 m and \Delta\omega = 50-60 milliarcseconds, respectively, over 1-2 yr. A numerical analysis shows also that the expected J2c^-2 range-rate signal for Juno should be as large as 280 microns per second during a typical 6 h pass at its closest approach. Independent analyses previously performed by other researchers about the measurability of the Lense-Thirring effect showed that the radio science apparatus of Juno should reach an accuracy in Doppler range-rate measurements of 1-5 microns per second over such passes. The range-rate signature of the classical even zonal perturbations is different from the 1PN one. Thus, further investigations, based on covariance analyses of simulated Doppler data and dedicated parameters estimation, are worth of further consideration. It turns out that the J2 c^-2 effects cannot be responsible of the flyby anomaly in the gravitational field of the Earth. A dedicated spacecraft in a 6678 km X 57103 km polar orbit would experience a geocentric J2 c^-2 range-rate shift of 0.4 mm s^-1.

Non-spherical systems described by MOND theories of modified gravity arising from generalizations of the Poisson equations are affected by a MONDian extra-quadrupolar potential \phi_M even if they are isolated and they are in deep Newtonian regime. In general MOND theories quickly approaching Newtonian dynamics for accelerations beyond A_0, \phi_M is proportional to a multiplicative scaling coefficient \alpha \sim 1, while in MOND models becoming Newtonian beyond \kappa A_0, \kappa >> 1, it is enhanced by \kappa^2. We analytically work out some orbital effects due to \phi_M in the framework of QUMOND, and compare them with the latest observational determinations of Solar System's planetary dynamics, exoplanets and double lined spectroscopic binary stars. The current admissible range for the anomalous perihelion precession of Saturn yields |\kappa| <= 3.5 x 10^3, while the radial velocity of \alpha Cen AB allows to infer |\kappa|<= 6.2 x 10^4 (A) and |\kappa|<= 4.2 x 10^4 (B). In evaluating such preliminary constraints it must be recalled that QUMOND is not the nonrelativistic limit of TeVeS.

A realistic assessment of the uncertainties in the even zonals of a given geopotential model must be made by directly comparing its coefficients with those of a wholly independent solution of superior formal accuracy. Otherwise, a favorable selective bias is introduced in the evaluation of the total error budget of the LAGEOS-based Lense-Thirring tests yielding likely too optimistic figures for it. By applying a novel approach which recently appeared in the literature, the second (L = 4) and the third (L = 6) even zonals turn out to be uncertain at a 2-3 10^-11 (L = 4) and 3-4 10^-11 (L = 6) level, respectively, yielding a total gravitational error of about 27-28%, with an upper bound of 37-39%. The results by Ries et al. themselves yield an upper bound for it of about 33%. The low-degree even zonals are not exclusively determined from the GRACE Satellite-to-Satellite Tracking (SST) range since they affect it with long-period, secular-like signatures over orbital arcs longer than one orbital period: GRACE SST is not accurately sensitive to such signals. Conversely, general relativity affects it with short-period effects as well. Thus, the issue of the a priori "imprinting" of general relativity itself in the GRACE-based models used so far remains open.

We analytically calculate the secular precession of the pericenter of a test particle orbiting a central body surrounded by a continuous distribution of Dark Matter (DM) by using some commonly adopted density profiles for it. We obtain exact expressions without resorting to a-priori simplifying assumptions on the orbital geometry of the test particle. Our formulas allow us to put constraints on the parameters of the DM distributions considered in several local astronomical and astrophysical scenarios, such as the Sun's planetary system, the double pulsar, and the stellar system around the supermassive black hole in Sgr A$^{\ast}$, all characterized by a wide variety of orbital configurations.

LETSGO (LEnse-Thirring Sun-Geo Orbiter) is a proposed space-based mission involving the use of a spacecraft moving along a highly eccentric heliocentric orbit perpendicular to the ecliptic. It aims to accurately measure some important physical properties of the Sun and to test some post-Newtonian features of its gravitational field by continuously monitoring the Earth-probe range. Preliminary sensitivity analyses show that, by assuming a cm-level accuracy in ranging to the spacecraft, it would be possible to detect, in principle, the Lense-Thirring effect on it at a 10^-3-10^-4 level over a timescale of 2 yr, while the larger Schwarzschild component of the solar gravitational field may be sensed with a relative accuracy of about 10^-8-10^-9 during the same temporal interval. The competing range perturbation due to the non-sphericity of the Sun would be a source of systematic error, but it turns out that all the three dynamical features of motion examined affect the Earth-probe range in different ways, allowing for a separation in data analyses. The high eccentricity would help in reducing the impact of the non-gravitational perturbations whose impact would certainly be severe when LETSGO would approach the Sun at just a few solar radii. Further studies should be devoted to investigate both the consequences of the non-conservative forces and the actual measurability of the effects of interest by means of extensive numerical data simulations, parameter estimations and covariance analyses. Also an alternative, fly-by configuration is worth of consideration.

We analytically work out several effects which a violation of the Strong Equivalence Principle (SEP) induces on the orbital motion of a binary system constituted of self-gravitating bodies immersed in a constant and uniform external field. We do not restrict to the small eccentricity limit. Moreover, we do not select any specific spatial orientation of the external polarizing field. We explicitly calculate the SEP-induced mean rates of change of all the osculating Keplerian orbital elements of the binary, the perturbation of the projection of the binary orbit onto the line-of-sight, the shift of the radial velocity, and the range and range-rate signatures and as well. We find that the ratio of the SEP precessions of the node and the inclination of the binary depends only on and the pericenter of the binary itself, being independent on both the magnitude and the orientation of the polarizing field, and on the semimajor axis, the eccentricity and the node of the binary. Our results, which do not depend on any particular SEP-violating theoretical scheme, can be applied to quite general astronomical and astrophysical scenarios. They can be used to better interpret present and future SEP experiments, especially when several theoretical SEP mechanisms may be involved, and to suitably design new dedicated tests.

A two-body system hypothetically affected by an additional radial acceleration H v_r, where v_r is the radial velocity of the binary's proper orbital motion, would experience long-term temporal changes of both its semimajor axis a and the eccentricity e qualitatively different from any other standard competing effect for them. Contrary to what one might reasonably expect, the analytical expressions of such rates do not vanish in the limit M--> 0, where M is the mass of the primary, being independent of it. This is a general requirement that any potentially viable physical mechanism able to provide such a putative acceleration should meet. Nonetheless, if H had the same value H_0 of the Hubble parameter at present epoch, such rates of change would have magnitude close to the present-day level of accuracy in determining planetary orbital motions in our Solar System. A tension with recent observations may even be present for Mercury and Mars. However, general relativity, applied to a localized gravitationally bound binary system immersed in an expanding Friedmann-Lemaitre-Robertson-Walker, does not predict the existence of such a putative radial acceleration at Newtonian level. Instead, it was recently shown in literature that an acceleration of order H and directed along the velocity v of the test particle occurs at post-Newtonian level. We worked out its orbital effects finding well-behaved secular rates of change for both a and e proportional to the Schwarzschild radius r_s of the primary. Their magnitude is quite small: the rate of change of a amounts to just 20 microns per century in our Solar System. Finally, we discussed certain basic criteria of viability that modified models of gravity should generally meet when their observable effects are calculated.

We use recent observations from solar system orbital motions in order to constrain f(T) gravity. In particular, imposing a quadratic f(T) correction to the linear-in-T form, which is a good approximation for every realistic case, we extract the spherical solutions of the theory. Using these spherical solutions to describe the Sun's gravitational field, we use recently determined supplementary advances of planetary perihelia, to infer upper bounds on the allowed f(T) corrections. We find that the maximal allowed divergence of the gravitational potential in f(T) gravity from that in the teleparallel equivalent of General Relativity is of the order of 6.2 \times 10^{-10}, in the applicability region of our analysis. This is much smaller than the corresponding (significantly small too) divergence that is predicted from cosmological observations, as expected. Such a tiny allowed divergence from the linear form should be taken into account in f(T) model building.

In this paper we critically analyze the so far performed and proposed tests for measuring the general relativistic Lense-Thirring effect in the gravitational field of the Earth with some of the existing accurately tracked artificial satellites. The impact of the 2nd generation GRACE-only EIGEN-GRACE02S Earth gravity model and of the 1st CHAMP+GRACE+terrestrial gravity combined EIGEN-CG01C Earth gravity model is discussed. The role of the proposed LARES is discussed as well.

The angular momentum of a star is an important astrophysical quantity related to its internal structure, formation and evolution. On average, helioseismology yields S = 1.92 10^41 kg m^2 s^-1 for the angular momentum of the Sun. We show how it should be possible to measure or, at least, constrain it in a near future by using the gravitomagnetic Lense-Thirring effect predicted by general relativity for the orbit of a test particle moving around a central rotating body. We also discuss the present-day situation in view of the latest determinations of the supplementary perihelion precession of Mercury.

We use latest data from solar system planetary orbital motions to put constraints on some Galileon-induced precessional effects. Due to the Vainshtein mechanism, the Galileon-type spherically symmetric field of a monopole induces a small, screened correction proprtional to \sqrt{r} to its usual r^-1 Newtonian potential which causes a secular precession of the pericenter of a test particle. In the case of our solar system, latest data from Mars allow to constrain the magnitude of such an interaction down to \alpha <= 0.3 level. Another Galileon-type effect which might impact solar system dynamics is due to an unscreened constant gradient induced by the peculiar motion of the Galaxy. The magnitude of such an effect, depending on the different gravitational binding energies of the Sun and the planets, is \xi <= 0.004 from the latest bounds on the supplementary perihelion precession of Saturn.

The Lense-Thirring precession of the longitude of perihelion of Mercury, as predicted by general relativity by using the value of the Sun's angular momentum S = 190 x 10^39 kg m^2 s^-1 from helioseismology, is -2.0 milliarcseconds per century, computed in a celestial equatorial reference frame. It disagrees at 4-{\sigma} level with the correction 0.4 +/- 0.6 milliarcseconds per century to the standard Newtonian/Einsteinian precession, provided that the latter is to be entirely attributed to frame-dragging. The supplementary precession was recently determined in a global fit with the INPOP10a ephemerides to a long planetary data record (1914-2010) including also 3 data points collected in 2008-2009 from the MESSENGER spacecraft. The INPOP10a models did not include the solar gravitomagnetic field at all, so that its signature might have partly been removed in the data reduction process. On the other hand, the Lense-Thirring precession may have been canceled to a certain extent by the competing precessions caused by small mismodeling in the quadrupole mass moment of the Sun and in the PPN parameter beta entering the Schwarzschild-like 1PN precession, both modeled in INPOP10a. On the contrary, the oblateness of Mercury itself has a negligible impact on its perihelion. The same holds for the mismodelled actions of both the largest individual asteroids and the ring of the minor asteroids. Future analysis of more observations from the currently ongoing MESSENGER mission will shed further light on such an issue which, if confirmed, might potentially challenge our present-day picture of the currently accepted laws of gravitation and/or of the physical properties of the Sun.

We analytically work out the long-term orbital perturbations induced by a homogeneous circular ring of radius Rr and mass mr on the motion of a test particle in the cases (I): r > R_r and (II): r < R_r. In order to extend the validity of our analysis to the orbital configurations of, e.g., some proposed spacecraftbased mission for fundamental physics like LISA and ASTROD, of possible annuli around the supermassive black hole in Sgr A* coming from tidal disruptions of incoming gas clouds, and to the effect of artificial space debris belts around the Earth, we do not restrict ourselves to the case in which the ring and the orbit of the perturbed particle lie just in the same plane. From the corrections to the standard secular perihelion precessions, recently determined by a team of astronomers for some planets of the Solar System, we infer upper bounds on mr for various putative and known annular matter distributions of natural origin (close circumsolar ring with R_r = 0.02-0.13 au, dust ring with R_r = 1 au, minor asteroids, Trans-Neptunian Objects). We find m_r <= 1.4 10^-4 m_E (circumsolar ring with R_r = 0.02 au), m_r <= 2.6 10^-6 m_E (circumsolar ring with R_r = 0.13 au), m_r <= 8.8 10^-7 m_E (ring with R_r = 1 au), m_r <= 7.3 10^-12 M_S (asteroidal ring with R_r = 2.80 au), m_r <= 1.1 <= 10^-11 M_S (asteroidal ring with R_r = 3.14 au), m_r <= 2.0 10^-8 M_S (TNOs ring with R_r = 43 au). In principle, our analysis is valid both for baryonic and non-baryonic Dark Matter distributions.

We numerically investigate the impact of GTR on the orbital part of the satellite-to-satellite range \rho and range-rate \dot\rho of the twin GRACE A/B spacecrafts through their dynamical equations of motion integrated in an Earth-centered frame over a time span \Delta t=1 d. Instead, the GTR effects connected with the propagation of the electromagnetic waves linking the spacecrafts are neglected. The present-day accuracies in measuring the GRACE biased range and range-rate are \sigma_\rho\sim 1-10 \mum, \sigma_\dot\rho\sim 0.1-1 \mum s^-1; studies for a follow-on of such a mission points toward a range-rate accuracy of the order of \sigma_\dot\rho\sim 1 nm s^-1 or better. The GTR range and range-rate effects turn out to be \Delta\rho=80 \mum and \Delta\dot\rho=0.012 \mum s^-1 (Lense-Thirring), and \Delta\rho=6000 \mum and \Delta\dot\rho=10 \mum s^-1 (Schwarzschild). We also compute the dynamical range and range-rate perturbations caused by the first six zonal harmonic coefficients J_L, L=2,3,4,5,6,7 of the classical multipolar expansion of the terrestrial gravitational potential in order to evaluate their aliasing impact on the relativistic effects. Conversely, we also quantitatively, and preliminarily, assess the possible a-priori \virg{imprinting} of GTR itself, not solved-for in all the GRACE-based Earth's gravity models produced so far, on the estimated values of the low degree zonals of the geopotential. The present sensitivity analysis can also be extended, in principle, to different orbital configurations in order to design a suitable dedicated mission able to accurately measure the relativistic effects considered.

We analytically work out the long-term orbital perturbations induced by the first term of the expansion of the perturbing potential arising from the local modification of the Newton's inverse square law due to a topology R^2 x S^1 with a compactified dimension of radius R recently proposed by Floratos and Leontaris. We neither restrict to any specific spatial direction for the asymmetry axis nor to particular orbital configurations of the test particle. Thus, our results are quite general. Non-vanishing long-term variations occur for all the usual osculating Keplerian orbital elements, apart from the semimajor axis which is left unaffected. By using recent improvements in the determination of the orbital motion of Saturn from Cassini data, we preliminarily inferred R >= 4-6 kau. As a complementary approach, the putative topological effects should be explicitly modeled and solved-for with a modified version of the ephemerides dynamical models with which the same data sets should be reprocessed.

We work out the impact that the recently determined time-dependent component of the Pioneer Anomaly (PA), interpreted as an additional exotic acceleration of gravitational origin with respect to the well known PA-like constant one, may have on the orbital motions of some planets of the solar system. By assuming that it points towards the Sun, it turns out that both the semi-major axis a and the eccentricity e of the orbit of a test particle experience secular variations. For Saturn and Uranus, for which modern data records cover at least one full orbital revolution, such predicted anomalies are up to 2-3 orders of magnitude larger than the present-day accuracies in empirical determinations their orbital parameters from the usual orbit determination procedures in which the PA was not modeled. Given the predicted huge sizes of such hypothetical signatures, it is unlikely that their absence from the presently available processed data can be attributable to an "absorption" for them in the estimated parameters caused by the fact that they were not explicitly modeled. The magnitude of a constant PA-type acceleration at 9.5 au cannot be larger than 9 10^-15 m s^-2 according to the latest observational results for the perihelion precession of Saturn.

It has been recently suggested by Dvali and Vikman that the superluminal neutrino phenomenology of the OPERA experiment may be due to an environmental feature of the Earth, naturally yielding a long-range fifth force of gravitational origin whose coupling with the neutrino is set by the scale M_*, in units of reduced Planck mass. Its characteristic length lambda should not be smaller than one Earth's radius R_e, while its upper bound is expected to be slightly smaller than the Earth-Moon distance (60 R_e). We analytically work out some orbital effects of a Yukawa-type fifth force for a test particle moving in the modified field of a central body. Our results are quite general since they are not restricted to any particular size of lambda; moreover, they are valid for an arbitrary orbital configuration of the particle, i.e. for any value of its eccentricity $e$. We find that the dimensionless strength coupling parameter alpha is constrained to |alpha| <= 1 10^-10-4 10^-9 for 1 R_e <= lambda <= 10 R_e by the laser data of the Earth's artificial satellite LAGEOS II, corresponding to M_* >= 4 10^9 -1.6 10^10. The Moon perigee allows to obtain |alpha| <= 3 10^-11 for the Earth-Moon pair in the range 15 R_e <= lambda <= 60 R_e, which translates as M_* >= 3 10^10 - 4.5 10^10. Our results are neither necessarily limited to the superluminal OPERA scenario nor to the Dvali-Vikman model, in which it is M_* = 10^-6 at lambda = 1 R_e, in contrast with our bounds: they generally extend to any theoretical scenario implying a fifth-force of Yukawa-type.

We analytically work out the long-term rates of change of the six osculating Keplerian orbital elements of a test particle acted upon by the Lorentz-violating gravitomagnetic acceleration due to a static body, as predicted by the Standard Model Extension (SME). We neither restrict to any specific spatial orientation for the symmetry-violating vector s nor make a priori simplifying assumptions concerning the orbital configuration of the perturbed test particle. Thus, our results are quite general, and can be applied for sensitivity analyses to a variety of specific astronomical and astrophysical scenarios. We find that, apart from the semimajor axis a, all the other orbital elements undergo non-vanishing secular variations. By comparing our results to the latest determinations of the supplementary advances of the perihelia of some planets of the solar system we preliminarily obtain s_x = (0.9 +/- 1.5) 10^-8, s_y = (-4 +/- 6) 10^-9, s_z = (0.3 +/- 1) 10^-9. Bounds from the terrestrial LAGEOS and LAGEOS II satellites are of the order of s\sim 10^-3-10^-4.

We analytically work out the long-term variations caused on the motion of a planet orbiting a star by a very distant, pointlike massive object X. Apart from the semi-major axis a, all the other Keplerian osculating orbital elements experience long-term variations which are complicated functions of the orbital configurations of both the planet itself and of X. We infer constraints on the minimum distance d_X at which X may exist by comparing our prediction of the long-term variation of the longitude of the perihelion \varpi to the latest empirical determinations of the corrections \Delta\dot\varpi to the standard Newtonian/Einsteinian secular precessions of several solar system planets recently estimated by independent teams of astronomers. We obtain the following approximate lower bounds on dX for the assumed masses of X quoted in brackets: 150 - 200 au (m_Mars), 250 - 450 au (0.7 m_Earth), 3500 - 4500 au (4 m_Jup).

We analytically compute the orbital effects induced on the motion of a spinning particle geodesically traveling around a central rotating body by the general relativistic two-body spin-spin and spin-orbit interactions. We use the weak-field and slow-motion leading-order approximation. Neither assumptions on the eccentricity nor the inclination of the orbit of the particle have been assumed. No specific orientations of the spins involved have been chosen. (Abridged)

A phenomenological anisotropic variation \Delta G/G of the Newtonian gravitational coupling parameter G, if real, would affect the orbital dynamics of a two-body gravitationally bound system in a specific way. We analytically work out the long-term effects that such a putative modification of the usual Newtonian inverse-square law would induce on the trajectory of a test particle orbiting a central mass. Without making any a-priori simplifying assumptions concerning the orbital configuration of the test particle, it turns out that its osculating semi-major axis a, eccentricity e, pericenter \varpi and mean anomaly M undergo long-term temporal variations, while the inclination I and the node \Omega are left unaffected. Moreover, the radial and the transverse components of the position and the velocity vectors r and v of the test particle experience non-vanishing changes per orbit, contrary to the out-of-plane ones. Then, we compute our theoretical predictions for some of the major bodies of the solar system by orienting the gradient of G(r) towards the Galactic Center and keeping it fixed over the characteristic timescales involved. By comparing our calculation to the latest observational determinations for the same bodies, we infer \Delta G/G <= 10^-17 over about 1 au. Finally, we consider also the Supermassive Black Hole hosted by the Galactic Center in Sgr A^* and the main sequence star S2 orbiting it in about 16 yr, obtaining just \Delta G/G <= 10^-2 over 1 kau.

We analytically and numerically investigate the long-term, i.e. averaged over one full revolution, orbital effects of the non-isotropic percent mass loss \dot m/m experienced by several transiting hot Jupiters whose atmospheres are hit by severe radiations flows coming from their close parent stars. The semi-major axis a, the argument of pericenter \omega and the mean anomaly M experience net variations, while the eccentricity e, the inclination I and the longitude of the ascending node remain unchanged, on average. In particular, a increases independently of e and of the speed Vesc of the ejected mass. By assuming |\dot m| <= 10^17 kg yr-1, corresponding to |\dot m/m| <= 10^-10 yr-1 for a Jupiter-like planet, it turns out \dot a = 2.5 m yr^-1 for orbits with a = 0.05 au. Such an effect may play a role in the dynamical history of the hot Jupiters, especially in connection with the still unresolved issue of the arrest of the planetary inward migrations after a distance a >= 0.01 au is reached. The retrograde pericenter variation depends, instead, on e and V_esc. It may, in principle, act as a source of systematic uncertainty in some proposed measurements of the general relativistic pericenter precession; however, it turns out to be smaller than it by several orders of magnitude.

The cancelation of the first even zonal harmonic coefficient J2 of the multipolar expansion of the Newtonian part of the Earth's gravitational potential from the linear combination f(2L) of the nodes of LAGEOS and LAGEOS II used in the latest tests of the general relativistic Lense-Thirring effect cannot be perfect, as assumed so far. It is so, among other things, because of the uncertainties in the spatial orientation of the terrestrial spin axis as well. As a consequence, the coefficient c1 entering f(2L), which is not a solve-for parameter being theoretically computed from the analytical expressions of the classical node precessions due to J2, is, on average, uncertain at a 10-8 level over multi-decadal time spans DT comparable to those used in the data analyses performed so far. A further \simeq 20% systematic uncertainty in the theoretically predicted gravitomagnetic signal, thus, occurs. The shift due to the gravitomagnetic frame-dragging on the station-spacecraft range is numerically computed over DT = 15 d and DT = 1 yr. The need of looking at such a directly observable quantity is pointed out, along with some critical remarks concerning the methodology used so far to measure the Lense-Thirring effect with the LAGEOS satellites. Suggestions for a different, more trustable and reliable approach are offered.

According to a recent modified model of gravity at large distances, a radial constant and uniform extra-acceleration A_Rin of Rindler type acts upon a test particle p in the static field of a central mass M if certain conditions are satisfied. Among other things, it was proposed as a potentially viable explanation of a part of the Pioneer anomaly. We study the impact that an anomalous Rindler-type term as large as |A_Rin| \sim 10^-10 m s^-2 may have on the the orbital dynamics of a typical object of the Oort cloud whose self-energy is quite smaller than its putative Rindler energy. By taking a typical comet moving along a highly eccentric and inclined orbit throughout the expected entire extension of the Oort cloud (0.02 pc-1 pc), it turns out that the addition of an outward Rindler-like acceleration, i.e. for A_Rin > 0, does not allow bound orbits. Instead, if A_Rin < 0, the resulting numerically integrated trajectory is limited in space, but it radically differs from the standard Keplerian ellipse. In particular, the heliocentric distance of the comet gets markedly reduced and experiences high frequency oscillations, its speed is increased, and the overall pattern of the trajectory is quite isotropic. As a consequence, the standard picture of the Oort cloud is radically altered since its modified orbits are much less sensitive to the disturbing actions of the Galactic tide and nearby passing stars whose effects, in the standard scenario, are responsible for the phenomenology on which our confidence in the existence of the cloud itself is based. The present analysis may be supplemented in future by further statistical Monte Carlo-type investigations by randomly varying the initial conditions of the comets.

We numerically work out the impact of the general relativistic Lense-Thirring effect on the Earth-Mercury range caused by the gravitomagnetic field of the rotating Sun. The peak-to peak nominal amplitude of the resulting time-varying signal amounts to 1.75 10^1 m over a temporal interval 2 yr. Future interplanetary laser ranging facilities should reach a cm-level in ranging to Mercury over comparable timescales; for example, the BepiColombo mission, to be launched in 2014, should reach a 4.5 - 10 cm level over 1 - 8 yr. We looked also at other Newtonian (solar quadrupole mass moment, ring of the minor asteroids, Ceres, Pallas, Vesta, Trans-Neptunian Objects) and post-Newtonian (gravitoelectric Schwarzschild solar field) dynamical effects on the Earth-Mercury range. They act as sources of systematic errors for the Lense-Thirring signal which, in turn, if not properly modeled, may bias the recovery of some key parameters of such other dynamical features of motion. Their nominal peak-to-peak amplitudes are as large as 4 10^5 m (Schwarzschild), 3 10^2 m (Sun's quadrupole), 8 10^1 m (Ceres, Pallas, Vesta), 4 m (ring of minor asteroids), 8 10^-1 m (Trans-Neptunian Objects). Their temporal patterns are different with respect to that of the gravitomagnetic signal.

We deal with the effects induced on the orbit of a test particle revolving around a central body by putative spatial variations of fundamental coupling constants $\zeta$. In particular, we assume a dipole gradient for $\zeta(\bds r)/\bar{\zeta}$ along a generic direction $\bds{\hat{k}}$ in space. We analytically work out the long-term variations of all the six standard Keplerian orbital elements parameterizing the orbit of a test particle in a gravitationally bound two-body system. It turns out that, apart from the semi-major axis $a$, the eccentricity $e$, the inclination $I$, the longitude of the ascending node $\Omega$, the longitude of pericenter $\pi$ and the mean anomaly $\mathcal{M}$ undergo non-zero long-term changes. By using the usual decomposition along the radial ($R$), transverse ($T$) and normal ($N$) directions, we also analytically work out the long-term changes $\Delta R,\Delta T,\Delta N$ and $\Delta v_R,\Delta v_T,\Delta v_N$ experienced by the position and the velocity vectors $\bds r$ and $\bds v$ of the test particle. It turns out that, apart from $\Delta N$, all the other five shifts do not vanish over one full orbital revolution. In the calculation we do not use \textit{a-priori} simplifying assumptions concerning $e$ and $I$. Thus, our results are valid for a generic orbital geometry; moreover, they hold for any gradient direction (abridged).

Both the recently reported anomalous secular increase of the astronomical unit, of the order of a few cm yr^-1, and of the eccentricity of the lunar orbit e_ = (9+/-3) 10^-12 yr^-1 can be phenomenologically explained by postulating that the acceleration of a test particle orbiting a central body, in addition to usual Newtonian component, contains a small additional radial term proportional to the radial projection vr of the velocity of the particle's orbital motion. Indeed, it induces secular variations of both the semi-major axis a and the eccentricity e of the test particle's orbit. In the case of the Earth and the Moon, they numerically agree rather well with the measured anomalies if one takes the numerical value of the coefficient of proportionality of the extra-acceleration approximately equal to that of the Hubble parameter H0 = 7.3 10^-11 yr^-1.

Possible explanations of the recently reported anomalous increase of the eccentricity of the lunar orbit are sought in terms of classical Newtonian mechanics, general relativity, and long-range modifications of gravity.

We analytically work out the orbital effects caused by a Rindlertype extra-acceleration ARin which naturally arises in some recent models of modified gravity at large distances. In particular, we focus on the perturbations induced by it on the two-body range {\rho} and range-rate {\rho}\cdot which are commonly used in satellite and planetary investigations as primary observable quantities. The constraints obtained for ARin by comparing our calculations with the currently available range and range-rate residuals for some of the major bodies of the solar system, obtained without explicitly modeling ARin, are 1 - 2 \times 10-13 m s-2 (Mercury and Venus), 1 \times 10-14 m s-2 (Saturn), 1 \times 10-15 m s-2 (Mars), while for a terrestrial Rindler acceleration we have 5 \times 10-16 m s-2 (Moon). Another approach which could be followed consists of taking into account ARin in re-processing all the available data sets with accordingly modified dynamical models, and estimating a dedicated solve-for parameter explicitly accounting for it. Anyway, such a method is time-consuming. A preliminary analysis likely performed in such a way by a different author yields A <= 8\times10-14 m s-2 at Mars' distance and A < = 1\times10-14 m s-2 at Saturn's distance. The method adopted here can be easily and straightforwardly extended to other long-range modified models of gravity as well.

We summarize some critical issues pertaining the tests of the general relativistic Lense-Thirring effect performed by I. Ciufolini and coworkers in the gravitational field of the Earth with the geodetic satellites LAGEOS and LAGEOS II tracked with the Satellite Laser Ranging technique.

A recent analysis of a Lunar Laser Ranging (LLR) data record spanning 38.7 yr revealed an anomalous increase of the eccentricity of the lunar orbit amounting to de/dt_meas = (9 +/- 3) 10^-12 yr^-1. The present-day models of the dissipative phenomena occurring in the interiors of both the Earth and the Moon are not able to explain it. We examine several dynamical effects, not modeled in the data analysis, in the framework of long-range modified models of gravity and of the standard Newtonian/Einsteinian paradigm. It turns out that none of them can accommodate de/dt_meas. Many of them do not even induce long-term changes in e; other models do, instead, yield such an effect, but the resulting magnitudes are in disagreement with de/dt_meas. In particular, the general relativistic gravitomagnetic acceleration of the Moon due to the Earth's angular momentum has the right order of magnitude, but the resulting Lense-Thirring secular effect for the eccentricity vanishes. A potentially viable Newtonian candidate would be a trans-Plutonian massive object (Planet X/Nemesis/Tyche) since it, actually, would affect e with a non-vanishing long-term variation. On the other hand, the values for the physical and orbital parameters of such a hypothetical body required to obtain the right order of magnitude for de/dt are completely unrealistic. Moreover, they are in neat disagreement with both the most recent theoretical scenarios envisaging the existence of a distant, planetary-sized body and with the model-independent constraints on them dynamically inferred from planetary motions. Thus, the issue of finding a satisfactorily explanation for the anomalous behavior of the Moon's eccentricity remains open.

We analytically work out the perturbation induced by the Kehagias-Sfetsos (KS) space-time solution of the Horava-Lifshitz (HL) modified gravity at long distances on the two-body range for a pair of test particles A and B orbiting the same mass M. We apply our results to the most recently obtained range-residuals \delta\rho for some planets of the solar system (Mercury, Mars, Saturn) ranged from the Earth to effectively constrain the dimensionsless KS parameter \psi_0 for the Sun. We obtain \psi_0 >= 7.2 x 10^-10 (Mercury), \psi_0 >= 9 x 10^-12 (Mars), \psi_0 >= 1.7 x 10^-12 (Saturn). Such lower bounds are tighter than other ones existing in literature by several orders of magnitude. We also preliminarily obtain \psi_0 >= 8 x 10^-10 for the system constituted by the S2 star orbiting the Supermassive Black Hole (SBH) in the center of the Galaxy.

The strategy followed so far in the performed or proposed tests of the general relativistic Lense-Thirring effect in the gravitational field of the Earth with laser-ranged satellites of LAGEOS type relies upon the cancelation of the disturbing huge precessions induced by the first even zonal harmonic coefficient J_2 of the multipolar expansion of the Newtonian part of the terrestrial gravitational potential by means of suitably designed linear combinations of the nodes \Omega of more than one spacecraft. Actually, such a removal does depend on the accuracy with which the coefficients of the combinations adopted can be realistically known. Uncertainties of the order of 2 cm in the semimajor axes a and 0.5 milliarcseconds in the inclinations I of LAGEOS and LAGEOS II, entering the expression of the coefficient c_1 of the combination of their nodes used so far, yield an uncertainty \delta c_1 = 1.30 10^-8. It gives an imperfectly canceled J_2 signal of 10.8 milliarcseconds per year corresponding to 23% of the Lense-Thirring signature. Uncertainties of the order of 10-30 microarcseconds in the inclinations yield \delta c_1=7.9 10^-9 which corresponds to an uncanceled J_2 signature of 6.5 milliarcseconds per year, i.e. 14% of the Lense-Thirring signal. Concerning a future LAGEOS-LAGEOS II-LARES combination with coefficients k_1 and k_2, the same uncertainties in a and the less accurate uncertainties in I as before yield \delta k_1=1.1 10^-8, \delta k_2=2 10^-9; they imply a residual J_2 combined precession of 14.7 milliarcseconds per year corresponding to 29% of the Lense-Thirring trend. Uncertainties in the inclinations at about 10 microarcseconds level give \delta k_1=5 10^-9, \delta k_2 = 2 10^-9; the uncanceled J_2 effect is 7.9 milliarcseconds per year, i.e. 16% of the relativistic effect.

We numerically investigate the impact on the two-body range by several Newtonian and non-Newtonian dynamical effects for some Earth-planet (Mercury, Venus, Mars, Jupiter, Saturn) pairs in view of the expected cm-level accuracy in some future planned or proposed interplanetary ranging operations (abridged).

We show that the retrograde perihelion precession of Saturn \Delta\dot\varpi, recently estimated by different teams of astronomers by processing ranging data from the Cassini spacecraft and amounting to some milliarcseconds per century, can be explained in terms of a localized, distant body X, not yet directly discovered. From the determination of its tidal parameter K = GM_X/r_X^3 as a function of its ecliptic longitude \lambda_X and latitude \beta_X, we calculate the distance at which X may exist for different values of its mass, ranging from the size of Mars to that of the Sun. The minimum distance would occur for X located perpendicularly to the ecliptic, while the maximum distance is for X lying in the ecliptic. We find for rock-ice planets of the size of Mars and the Earth that they would be at about 80-150 au, respectively, while a Jupiter-sized gaseous giant would be at approximately 1 kau. A typical brown dwarf would be located at about 4 kau, while an object with the mass of the Sun would be at approximately 10 kau, so that it could not be Nemesis for which a solar mass and a heliocentric distance of about 88 kau are predicted. If X was directed towards a specific direction, i.e. that of the Galactic Center, it would mimick the action of a recently proposed form of the External Field Effect (EFE) in the framework of the MOdified Newtonian Dynamics (MOND).

We analytically compute the effects that a mass variation rate \dot M/M of a pulsar may have on the changes \Delta\tau in the times of arrival of its pulses due to test particle companions, and on their orbital dynamics. We apply our results to the planetary system of the PSR B1257+12 pulsar, located in the Galaxy at 600 pc from us, to phenomenologically constrain a putative accretion of non-annihilating dark matter on the hosting neutron star. (Abridged)

We focus on Horava-Lifshitz (HL) theory of gravity, and, in particular, on the Kehagias and Sfetsos s solution that is the analog of Schwarzschild black hole of General Relativity. In the weak-field and slow-motion approximation we analytically work out the secular precession of the longitude of the pericentre of a test particle induced by this solution. Its analytical form is different from that of the general relativistic Einstein's pericentre precession. Then, we compare it to the latest determinations of the corrections to the standard Newtonian/ Einsteinian planetary perihelion precessions recently estimated by E.V. Pitjeva with the EPM2008 ephemerides. It turns out that the planets of the solar system, taken singularly one at a time, allow to put lower bounds on the adimensional HL parameter \psi_0 of the order of 10^-12 (Mercury) 10^-24 (Pluto). They are not able to account for the Pioneer anomalous acceleration for r > 20 AU.

We analytically work out the long-term, i.e. averaged over one orbital revolution, time variations of some direct observable quantities Y induced by classical and general relativistic dynamical perturbations of the two-body pointlike Newtonian acceleration in the case of transiting exoplanets moving along elliptic orbits. More specifically, the observables $Y$ with which we deal are the transit duration, the radial velocity and the time interval between primary and secondary eclipses. The dynamical effects considered are the centrifugal oblateness of both the star and the planet, their tidal bulges mutually raised on each other, a distant third body X, and general relativity (both Schwarzschild and Lense-Thirring). We take into account the effects due to the perturbations of all the Keplerian orbital elements involved in a consistent and uniform way. First, we explicitly compute their instantaneous time variations due to the dynamical effects considered and plug them in the general expression for the instantaneous change of Y; then, we take the overall average over one orbital revolution of the so-obtained instantaneous rate $\dot Y(t)$ specialized to the perturbations considered. Instead, somewhat hybrid expressions can be often found in literature: in them, the secular precession of, typically, the periastron only is straightforwardly inserted into instantaneous formulas. Numerical evaluations of the obtained results are given for a typical star-planet scenario and compared with the expected observational accuracies over a time span 10 yr long. Our results are, in principle, valid also for other astronomical scenarios. They may allow, e.g., for designing various tests of gravitational theories with natural and artificial bodies in our solar system. (Abridged)

We consider a spherically symmetric and asymptotically flat vacuum solution of the Horava-Lifshitz (HL) gravity that is the analog of the general relativistic Schwarzschild black hole. In the weak-field and slow-motion approximation, we work out the correction to the third Kepler law of a test particle induced by such a solution and compare it to the phenomenologically determined orbital period of the transiting extrasolar planet HD209458b Osiris to preliminarily obtain an order-of-magnitude lower bound on the KS dimensionless parameter \omega_0 >= 1.4\times 10^-18. As suggestions for further analyses, the entire data set of HD209458b should be re-processed by explicitly modeling KS gravity as well, and one or more dedicated solve-for parameter(s) should be estimated.

Recent years have seen increasing efforts to directly measure some aspects of the general relativistic gravitomagnetic interaction in several astronomical scenarios in the solar system. After briefly overviewing the concept of gravitomagnetism from a theoretical point of view, we review the performed or proposed attempts to detect the Lense-Thirring effect affecting the orbital motions of natural and artificial bodies in the gravitational fields of the Sun, Earth, Mars and Jupiter. In particular, we will focus on the evaluation of the impact of several sources of systematic uncertainties of dynamical origin to realistically elucidate the present and future perspectives in directly measuring such an elusive relativistic effect.

We analytically work out the cumulative, i.e. averaged over one orbital revolution, time variations of the radial velocity v_r of a typical S star orbiting the Supermassive Black Hole (SBH) hosted by the Galactic Center (GC) in Sgr A* caused by several dynamical effects. They are the general relativistic gravitoelectromagnetic (GEM) fields of the SBH, its quadrupole mass moment Q2, and a diffuse dark matter distribution around the SBH. All of them induce non-zero secular radial accelerations proportional to the eccentricity e of the orbit. By taking the S2 star, orbiting the SBH along a highly eccentric (e = 0.8831) ellipse with a period Pb = 15.9 yr, we numerically compute the magnitudes of its radial accelerations (Abridged).

We analytically and numerically investigate the possibility that a still undiscovered body X, moving along an unbound hyperbolic path from outside the solar system, may penetrate its inner regions in the next few years posing a threat to the Earth. By conservatively using as initial position the lower bounds on the present-day distance dX of X dynamically inferred from the gravitational perturbations induced by it on the orbital motions of the planets of the solar system, both the analyses show that, in order to reach the Earth's orbit in the next 2 yr, X should move at a highly unrealistic speed v, whatever its mass MX is. For example, by assuming for it a solar (MX =M_Sun) or brown dwarf mass (MX = 80mJup), now at not less than dX = 11-6 kau (1 kau=1000 astronomical units), v would be of the order of 6-10% and 3-5% of the speed of light c, respectively. By assuming larger present-day distances for X, on the basis of the lacking of direct observational evidences of electromagnetic origin for it, its speed would be even higher. Instead, the fastest solitary massive objects known so far, like hypervelocity stars (HVSs) and supernova remnants (SRs), travel at v = 0.002-0.005c, having acquired so huge velocities in some of the most violent astrophysical phenomena like interactions with supermassive galactic black holes and supernova explosions. It turns out that the orbit of the Earth would not be macroscopically altered by a close (0.2 au) passage of such an ultrafast body X in the next 2 yr. On the contrary, our planet would be hurled into the space if a Sun-sized body X would encounter it by moving at v/c = 10^-4. On the other hand, this would imply that such a X should be now at just 20-30 au, contrary to all direct observational and indirect dynamical evidences.

We numerically investigate the features of typical orbits occurring in the Oort cloud (r\approx 50-150 kAU) in the low-acceleration regime of the MOdified Newtonian Dynamics (MOND). We take into account the so-called External Field Effect (EFE) because the solar system is embedded in the Milky Way. In the framework of MOND this does matter since the gravitational acceleration of Galactic origin felt by the solar system is of the same order of magnitude of the characteristic MOND acceleration scale A_0\approx 10^-10 m s^-2. We use three different forms of the MOND interpolating function \mu(x) and two different values for the Galactic field at the Sun's location. We find that MOND produces highly distorted trajectories with respect to the Newtonian case, especially for very eccentric orbits. It turns out that the shape of the MOND orbits strongly depend on the initial conditions. For particular initial state vectors, the MOND paths in the ecliptic plane get shrunk extending over much smaller spatial regions than in the Newtonian case, and experience high frequency variations over one Keplerian orbital period. Ecliptic orbits with different initial conditions and nearly polar orbits are quite different getting distorted as well, but they occupy more extended spatial regions. This facts may have consequences on the composition and the dynamical history of the Oort cloud which are difficult to predict.

WASP-33 is a fast rotating, main sequence star which hosts a hot Jupiter moving along a retrograde and almost polar orbit with semi-major axis a = 0.02 au and eccentricity provisionally set to e = 0. The quadrupole mass moment J_2 and the proper angular momentum S of the star are 1900 and 400 times, respectively, larger than those of the Sun. Thus, huge classical and relativistic non-Keplerian orbital effects should take place in such a system. In this paper we investigate the perspectives in detecting them (Abridged)

We investigate on the relative inclination of the planets B and C orbiting the pulsar PSR B1257+12 in connection with potential violations of the equivalence principle (Abridged).

We phenomenologically put local constraints on the rotation of distant masses by using the planets of the solar system. First, we analytically compute the orbital secular precessions induced on the motion of a test particle about a massive primary by a Coriolis-like force, treated as a small perturbation of first order in the rotation, in the case of a constant angular velocity vector \Psi directed along a generic direction in space. The semimajor axis a and the eccentricity e of the test particle do not secularly precess, contrary to the inclination I, the longitude of the ascending node \Omega, the longitude of the pericenter \varpi and the mean anomaly M. Then, we compare our prediction for <\dot\varpi> with the corrections \Delta\dot\varpi to the usual perihelion precessions of the inner planets recently estimated by fitting long data sets with different versions of the EPM ephemerides. We obtain |\Psi_z| <= 0.0006-0.013 arcsec cty^-1, |\Psi_x| <= 0.1-2.7 arcsec cty-1, |\Psi_y| <= 0.3-2.3 arcsec cty^-1. Interpreted in terms of models of space-time involving cosmic rotation, our results are able to yield constraints on cosmological parameters like the cosmological constant \Lambda and the Hubble parameter H_0 not too far from their values determined with cosmological observations and, in some cases, several orders of magnitude better than the constraints usually obtained so far from space-time models not involving rotation. In the case of the rotation of the solar system throughout the Galaxy, occurring clockwise about the North Galactic Pole, our results for \Psi_z are in disagreement with the expected value of it at more than 3-\sigma level.

We, first, analytically work out the long-term, i.e. averaged over one orbital revolution, perturbations on the orbit of a test particle moving in a local Fermi frame induced therein by the cosmological tidal effects of the inhomogeneous Lemaitre-Tolman-Bondi (LTB) model. The LTB solution has recently attracted attention, among other things, as a possible explanation of the observed cosmic acceleration without resorting to dark energy. Then, we phenomenologically constrain both the parameters K_1 = -\ddot R/R and K_2 = -\ddot R^'/R^' of the LTB metric in the Fermi frame by using different kinds of solar system data. The corrections $\Delta\dot\varpi$ to the standard Newtonian/Einsteinian precessions of the perihelia of the inner planets recently estimated with the EPM ephemerides, compared to our predictions for them, yield K_1 = (4+8) 10^-26 s^-2, K_2 = (3+7) 10^-23 s^-2. The residuals of the Cassini-based Earth-Saturn range, compared with the numerically integrated LTB range signature, allow to obtain K_1/2 = 10^-27 s^-2. The LTB-induced distortions of the orbit of a typical object of the Oort cloud with respect to the commonly accepted Newtonian picture, based on the observations of the comet showers from that remote region of the solar system, point towards K_1/2 <= 10^-30-10^-32 s^-2. Such figures have to be compared with those inferred from cosmological data which are of the order of K1 \approx K2 = -4 10^-36 s^-2.

We apply our recent results on orbital dynamics around a mass-varying central body to the phenomenon of accretion of Dark Matter-assumed not self-annihilating-on the Sun and the major bodies of the solar system due to its motion throughout the Milky Way halo. We inspect its consequences on the orbits of the planets and their satellites over timescales of the order of the age of the solar system. It turns out that a solar Dark Matter accretion rate of \approx 10^-12 yr^-1, inferred from the upper limit \Delta M/M= 0.02-0.05 on the Sun's Dark Matter content, assumed somehow accumulated during last 4.5 Gyr, would have displaced the planets faraway by about 10^-2-10^1 au 4.5 Gyr ago. Another consequence is that the semimajor axis of the Earth's orbit, approximately equal to the Astronomical Unit, would undergo a secular increase of 0.02-0.05 m yr^-1, in agreement with the latest observational determinations of the Astronomical Unit secular increase of 0.07 +/- 0.02 m yr^-1 and 0.05 m yr^-1. By assuming that the Sun will continue to accrete Dark Matter in the next billions year at the same rate as in the past, the orbits of its planets will shrink by about 10^-1-10^1 au (\approx 0.2-0.5 au for the Earth), with consequences for their fate, especially of the inner planets. On the other hand, lunar and planetary ephemerides set upper bounds on the secular variation of the Sun's gravitational parameter GM which are one one order of magnitude smaller than 10^-12 yr^-1. Dark Matter accretion on planets has, instead, less relevant consequences for their satellites. Indeed, 4.5 Gyr ago their orbits would have been just 10^-2-10^1 km wider than now. (Abridged)

Binary systems containing at least one radiopulsar are excellent laboratories to test several aspects of fundamental physics like matter properties in conditions of extreme density and theories of gravitation like the Einstein's General Theory of Gravitation (GTR) along with modifications/extensions of it. In this Chapter we focus on the perspectives on measuring the moment of inertia of the double pulsar, its usefulness in testing some modified models of gravity, and the possibility of using the mean anomaly as a further post-Keplerian orbital parameter to probe GTR.

We study the effects of a time-varying gravitomagnetic field on the motion of test particles. Starting from recent results, we consider the gravitomagnetic field of a source whose spin angular momentum has a linearly time-varying magnitude. The acceleration due to such a time-varying gravitomagnetic field is considered as a perturbation of the Newtonian motion, and we explicitly evaluate the effects of this perturbation on the Keplerian elements of a closed orbit. The theoretical predictions are compared with actual astronomical and astrophysical scenarios, both in the solar system and in binary pulsars systems, in order to evaluate the impact of these effects on real systems.

If the Pioneer Anomaly (PA) was a genuine dynamical effect of gravitational origin, it should also affect the orbital motions of the solar system's bodies moving in the space regions in which the PA manifested itself in its presently known form, i.e. as a constant and uniform acceleration approximately directed towards the Sun with a non-zero magnitude APio = (8.74 +/- 1.33) x 10^-10 m s^-2 after 20 au from the Sun. In this paper we preliminarily investigate its effects on the orbital motions of the Neptunian satellites Triton, Nereid and Proteus, located at about 30 au from the Sun, both analytically and numerically. Extensive observational records covering sev- eral orbital revolutions have recently been analyzed for them, notably improving the knowledge of their orbits. Both analytical and numerical calculations, limited to the direct, Neptune-satellite interaction, show that the peak-to-peak amplitudes of the PA-induced radial, transverse and out-of-plane perturbations over one century are up to 300 km, 600 km, 8 m for Triton, 17,500 km, 35,000 km, 800 km for Nereid, and 60 km, 120 km, 30 m for Proteus. The corresponding orbital uncertainties ob- tained from a recent analysis of all the data available for the satellites considered are, in general, smaller by one-two orders of magnitude, although obtained without modeling a Pioneer-like extra-force. Further investigations based on a re-processing of the satellites' real or simulated data with modified equations of motions including an additional Pioneer-type force as well are worth being implemented and may shed further light on this important issue.

The goal of the recently approved space-based LARES mission is to measure the general relativistic Lense-Thirring effect in the gravitational field of the spinning Earth at a repeatedly claimed 1% accuracy by combining its node Omega with those of the existing LAGEOS and LAGEOS II laser-ranged satellites. In this paper we show that, in view of the lower altitude of LARES (h=1450 km) with respect to LAGEOS and LAGEOS II (h\approx 6000 km), the cross-coupling between the effect of the atmospheric drag, both neutral and charged, on the inclination of LARES and its classical node precession due to the Earth's oblateness may induce a 3-9% year^-1 systematic bias on the total relativistic precession. Since its extraction from the data will take about 5-10 years, such a perturbing effect may degrade the total accuracy of the test, especially in view of the large uncertainties in modeling the drag force.

We investigate the effect of possible a-priori "imprinting" effects of general relativity itself on satellite/spaceraft-based tests of it. We deal with some performed or proposed time-delay ranging experiments in the Sun's gravitational field. It turns out that the "imprint" of general relativity on the Astronomical Unit and the solar gravitational constant GM_{\odot}, not solved for in the so far performed spacecraft-based time-delay tests, induces an a-priori bias of the order of 10^-6 in typical solar system ranging experiments aimed to measuring the space curvature PPN parameter gamma. It is too small by one order of magnitude to be of concern for the performed Cassini experiment, but it would affect future planned or proposed tests aiming to reach a 10^-7-10^-9 accuracy in determining gamma.

After extensively reviewing general relativistic gravitomagnetism, both historically and phenomenologically, we review in detail the so-called magnetic components of gravitational waves (GWs), which have to be taken into account in the context of the total response functions of interferometers for GWs propagating from arbitrary directions. Following the more recent approaches of this important issue, the analysis of such magnetic components will be reviewed in both of standard General Theory of Relativity (GTR) and Scalar Tensor Gravity. Thus, we show in detail that such a magnetic component becomes particularly important in the high-frequency portion of the range of ground based interferometers for GWs which arises from the two different theories of gravity. Our reviewed results show that if one neglects the magnetic contribution to the gravitational field of a GW, approximately 15% of the potential observable signal could, in principle, be lost.

The impact of possible a-priori "imprinting" effects of general relativity itself on recent attempts to measure the Lense-Thirring precessions with the LAGEOS satellites orbiting the Earth and the terrestrial geopotential models by the dedicated mission GRACE is investigated. It is analytically shown that general relativity, not explicitly solved for in the GRACE-based models, may "imprint" their even zonal harmonic coefficients J_L at a non-negligible level, given the present-day accuracy in recovering them. This translates into a bias of the LAGEOS-based relativistic tests as large as the Lense-Thirring effect itself. Further analyses should include general relativity itself in the GRACE data processing by explicitly solving for it.

The recently approved Juno mission will orbit Jupiter for one year in a highly eccentric (r_min=1.06R_Jup, r_max=39R_Jup) polar orbit (i=90 deg) to accurately map, among other things, the jovian magnetic and gravitational fields. Such an orbital configuration yields an ideal situation, in principle, to attempt a measurement of the general relativistic Lense-Thirring effect through the Juno's node Omega which would be displaced by about 570 m over the mission's duration. Conversely, by assuming the validity of general relativity, the proposed test can be viewed as a direct, dynamical measurement of the Jupiter's angular momentum S which would give important information concerning the internal structure and formation of the giant planet. The long-period orbital perturbations due to the zonal harmonic coefficients J_L, L=2,3,4,6 of the multipolar expansion of the jovian gravitational potential accounting for its departures from spherical symmetry are a major source of systematic bias. While the Lense-Thirring node rate is independent of the inclination i, the node zonal perturbations vanish for i=90. In reality, the orbit injection errors will induce departures \delta i from the ideal polar geometry, so that the zonal perturbations will come into play at an unacceptably high level, in spite of the expected improvements in the low-degree zonals by Juno. A linear combination of Omega, the periJove omega and the mean anomaly M cancels out the impact of J_2 and J_6. A two orders of magnitude improvement in the uncanceled J_3 and J_4 would be needed to reduce their bias on the relativistic signal to the percent level; it does not seem unrealistic because the expected level of improvement in such zonals is three orders of magnitude.

I calculate the classical effects induced by an isotropic mass loss of a body on the orbital motion of a test particle around it. By applying my results to the phase in which the radius of the Sun, already moved to the Red Giant Branch of the Hertzsprung-Russell Diagram, will become as large as 1.20 AU in about 1 Myr, I find that the Earth's perihelion position on the fixed line of the apsides will increase by about 0.22-0.25 AU (for \dot M/M = 2 x 10^-7 yr^-1); other researchers point towards an increase of 0.37-0.63 AU. Mercury will be destroyed already at the end of the Main Sequence, while Venus should be engulfed in the initial phase of the Red Giant Branch phase; the orbits of the outer planets will increase by 1.2-7.5 AU. Simultane- ous long-term numerical integrations of the equations of motion of all the major bodies of the solar system, with the inclusion of a mass-loss term in the dynamical force models as well, are required to check if the mutual N-body interactions may substantially change the picture analytically outlined here, especially in the Red Giant Branch phase in which Mercury and Venus may be removed from the integration.

I calculate the classical effects induced by an isotropic mass loss of a body on the orbital motion of a test particle around it; the present analysis is also valid for a variation of the Newtonian constant of gravitation. I perturbatively obtain negative secular rates for the osculating semimajor axis, the eccentricity and the mean anomaly, while the argument of pericenter does not undergo secular precession; the node and the inclination i remain unchanged. The anomalistic period is different from the Keplerian one, being larger than it. The true orbit, instead, expands, as shown by a numerical integration of the equations of motion in Cartesian coordinates; in fact, this is in agreement with the seemingly counter-intuitive decreasing of the semimajor axis and the eccentricity because they refer to the osculating Keplerian ellipses which approximate the trajectory at each instant. A comparison with the results obtained with different approaches by other researchers is made. General relativity induces positive secular rates of the semimajor axis and the eccentricity completely negligible in the present and future evolution of the solar system.

We investigate possible a-priori "imprinting" of general relativity itself on spaceraft-based tests of it. We deal with some performed or proposed time-delay ranging experiments in the Sun's gravitational field. The "imprint" of general relativity on the Astronomical Unit and the solar gravitational constant GM_\odot, not solved for in the so far performed spacecraft-based time-delay tests, may induce an a-priori bias of the order of 10^-6 in typical solar system ranging experiments aimed to measuring the space curvature PPN parameter \gamma. It is too small by one order of magnitude to be of concern for the performed Cassini experiment, but it would affect future planned or proposed tests aiming to reach a 10^-7-10^-9 accuracy in determining \gamma.

We tentatively look at anthropic constraints on the Cosmological Constant (CC) \Lambda at galactic scales by investigating its influence on the motion of the Sun throughout the Milky Way (MW) for -4.5 <= t <=0 Gyr. In particular, we look at the Galactocentric distance at which the Sun is displaced at the end of the numerical integration of its equations of motion modified in order to include the effect of \Lambda as well. Values of it placing our star at its birth at more than 10 kpc from the Galactic center (GC) are to be considered implausible, according to the current views on the Galactic Habitable Zone (GHZ) on the metallicity level needed for stars' formation. Also values yielding too close approaches to GC should be excluded because of the risks to life's evolution coming from too much nearby supernovae (SN) explosions and Gamma Ray Bursts (GRB). We investigate the impact on our results of the uncertainties on both the MW model's parameters and the Sun's initial conditions, in particular the Hubble parameter H_0 and the Local Standard Rest (LSR) speed \Theta_0 accurate at 2% and 6.2% level, respectively. While H_0=70.1 km s^-1 Mpc^-1, \Theta_0=254 km s^-1 and \Lambda <= 10^-55 cm^-2 locates the place of birth of the Sun at 19.6 kpc from GC, the same values for H_0 and \Lambda, and \Theta_0^max=270 km s^-1, places it at the plausible Galactocentric distance of 8.5 kpc. \Lambda = 10^-54 cm^-2 and \Lambda = 10^-53 cm^-2 place the Sun at 10.6 kpc and 18.7 kpc, respectively.

According to general relativity, a spinning body of mass M and angular momentum S, like a star or a planet, generates a gravitomagnetic field which induces, among other phenomena, also the Lense-Thirring effect, i.e. secular precessions of the path of a test particle orbiting it. Direct and indisputable tests of such a relativistic prediction are still missing. We discuss some performed attempts to measure it in the gravitational fields of several bodies in the Solar System with natural and artificial objects. The focus is on the realistic evaluation of the impact of some competing classical forces regarded as sources of systematic uncertainties degrading the total accuracy obtainable.

We simultaneously integrate in a numerical way the equations of motion of both the Magellanic Clouds (MCs) in MOND, MOG and CDM for -1 <= t <= +1 Gyr to see if, at least in principle, it is possible to discriminate between them (Abridged version).

We looked at the orbital motions of test particles according to the External Field Effect (EFE) predicted by the MOdified Newtonian Dynamics (MOND) in the Oort cloud which falls in the deep MONDian regime (r\approx 50-150 kAU). Concerning the interpolating function \mu(x), we extensively used the forms \mu_1=1/(1+x),\mu_2=x/(1+x^2)^1/2,\mu_3/2=x/(1+x^3/2)^2/3. We integrated both the MOND and the Newtonian equations of motion in Cartesian coordinates sharing the same initial conditions. We considered both ecliptic and nearly polar trajectories, all with high eccentricities (e>0.1). In order to evaluate the characteristic MOND parameters \mu_g and L_g entering the problem, we used two different values (V=220 km s^-1 and V=254 km s^-1) of the circular speed of the solar system's motion through the Galaxy; $V$ allows to evaluate the Milky Way's gravitational field at the Sun's location. It turns out that EFE induces strong distortions of the Newtonian ellipses, especially in the ecliptic plane yielding more involved paths which span less extended spatial regions.

We numerically integrate the equations of motion of the Sun in Galactocentric Cartesian rectangular coordinates for -4.5 Gyr <= t <= 0 in Newtonian mechanics with two different models for the Cold Dark Matter (CDM) halo, in MOdified Newtonian Dynamics (MOND) and in MOdified Gravity (MOG) without resorting to CDM. The initial conditions used come from the latest kinematical determination of the 3D Sun's motion in the Milky Way (MW) by assuming for the rotation speed of the Local Standard of Rest (LSR) the recent value \Theta_0=268 km s^-1 and the IAU recommended value \Theta_0=220 km s^-1; the Sun is assumed located at 8.5 kpc from the Galactic Center (GC). For \Theta_0=268 km s^-1 the birth of the Sun, 4.5 Gyr ago, would have occurred at large Galactocentric distances (12-27 kpc depending on the model used), while for \Theta_0=220 km s^-1 it would have occurred at about 8.8-9.3 kpc for almost all the models used. The integrated trajectories are far from being circular, especially for \Theta_0=268 km s^-1, and differ each other with the CDM models yielding the widest spatial extensions for the Sun's orbital path.

We put full 3D constraints on a putative planet X by using the dynamics of the inner planets of the solar system. In particular, we compute the mimium distance of X as a function of its heliocentric latitude and longitude for different values of its mass.

After an introduction emphasizing the importance of the gravitomag- netic effect in general relativity, with a resume of some space-based appli- cations, we discuss the so-called magnetic components of gravitational waves (GWs), which have to be taken into account in the context of the total response functions of interferometers for GWs propagating from ar- bitrary directions.

I discuss some aspects of the recent test of frame-dragging performed by me by exploiting the Root-Mean-Square (RMS) orbit overlap differences of the out-of-plane component N of the orbit of the Mars Global Surveyor (MGS) spacecraft in the gravitational field of Mars. A linear fit of the full time series of the entire MGS data (4 February 1999-14 January 2005) yields a normalized slope 1.03 +/- 0.41 (with 95% confidence bounds). Other linear fits to different data sets confirm the agreement with general relativity. The huge systematic effects induced by the mismodeling in the martian gravitational field claimed by some authors are absent in the MGS out-of-plane record. The non-gravitational forces affect at the same level of the gravitomagnetic one the in-plane orbital components of MGS, not the out-of-plane one. Moreover, they experience high-frequency variations which does not matter in the present case in which secular effects are relevant.

In the framework of the MOdified Newtonian Dynamics (MOND), the internal dynamics of a gravitating system s embedded in a larger one S is affected by the external background field E of S even if it is constant and uniform, thus implying a violation of the Strong Equivalence Principle: it is the so-called External Field Effect (EFE). In the case of the solar system, E would be A_cen\approx 10^-10 m s^-2 because of its motion through the Milky Way: it is orders of magnitude smaller than the main Newtonian monopole terms for the planets. We address here the following questions in a purely phenomenological manner: are the Sun's planets affected by an EFE as large as 10^-10 m s^-2? Can it be assumed that its effect is negligible for them because of its relatively small size? Does $\bds E$ induce vanishing net orbital effects because of its constancy over typical solar system's planetary orbital periods? It turns out that a constant and uniform acceleration, treated perturbatively, does induce non-vanishing long-period orbital effects on the longitude of the pericenter of a test particle. In the case of the inner planets of the solar system and with E\approx 10^-10 m s^-2, they are 4-6 orders of magnitude larger than the present-day upper bounds on the non-standard perihelion precessions recently obtained with by E.V. Pitjeva with the EPM ephemerides in the Solar System Barycentric frame. The upper limits on the components of E are E_x <= 1 x 10^-15 m s^-2, E_y <= 2 x 10^-16 m s^-2, E_z <= 3 x 10^-14 m s^-2. This result is in agreement with the violation of the Strong Equivalence Principle by MOND.

Many long-range modifications of the Newtonian/Einsteinian standard laws of gravity have been proposed in the recent past to explain various celestial phenomena occurring at different scales ranging from solar system to the entire universe. The most famous ones are the so-called Pioneer anomaly, {i.e.} a still unexplained acceleration detected in the telemetry of the Pioneer 10/11 spacecraft after they passed the 20 AU threshold in the solar system, the non-Keplerian profiles of the velocity rotation curves of several galaxies and the cosmic acceleration. We use the latest observational determinations of the planetary motions in the solar system and in the double pulsar system to put constraints on such models independently of the phenomena for which they were originally proposed. We also deal with the recently detected anomalous perihelion precession of Saturn and discuss the possibility that it can be explained by some of the aforementioned models of modified gravity.

Analytical and numerical calculations show that a putative temporal variation of the speed of light c, with the meaning of space-time structure constant c_ST, assumed to be linear over timescales of about one century, would induce a secular precession of the longitude of the pericenter \varpi of a test particle orbiting a spherically symmetric body. By comparing such a predicted effect to the corrections \Delta\dot\varpi to the usual Newtonian/Einsteinian perihelion precessions of the inner planets of the Solar System, recently estimated by E.V. Pitjeva by fitting about one century of modern astronomical observations with the standard dynamical force models of the EPM epehemerides, we obtained \dot c/c =(0.5 +/- 2)\times 10^-7 yr^-1. Moreover, the possibility that \dot c/c\neq 0 over the last century is ruled out at 3-12\sigma level by taking the ratios of the perihelia for different pairs of planets. Our results are independent of any measurement of the variations of other fundamental constants which may be explained by a variation of $c$ itself (with the meaning of electromagnetic constant c_EM). It will be important to repeat such tests if and when other teams of astronomers will estimate their own corrections to the standard Newtonian/Einsteinian planetary perihelion precessions.

We deal with the test of the general relativistic gravitomagnetic Lense-Thirring effect currently ongoing in the Earth's gravitational field with the combined nodes \Omega of the laser-ranged geodetic satellites LAGEOS and LAGEOS II. One of the most important source of systematic uncertainty on the orbits of the LAGEOS satellites, with respect to the Lense-Thirring signature, is the bias due to the even zonal harmonic coefficients J_L of the multipolar expansion of the Earth's geopotential which account for the departures from sphericity of the terrestrial gravitational potential induced by the centrifugal effects of its diurnal rotation. The issue addressed here is: are the so far published evaluations of such a systematic error reliable and realistic? The answer is negative. Indeed, if the difference \Delta J_L among the even zonals estimated in different global solutions (EIGEN-GRACE02S, EIGEN-CG03C, GGM02S, GGM03S, ITG-Grace02, ITG-Grace03s, JEM01-RL03B, EGM2008, AIUB-GRACE01S) is assumed for the uncertainties \delta J_L instead of using their more or less calibrated covariance sigmas \sigma_{J_L}, it turns out that the systematic error \delta\mu in the Lense-Thirring measurement is about 3 to 4 times larger than in the evaluations so far published based on the use of the sigmas of one model at a time separately, amounting up to 37% for the pair EIGEN-GRACE02S/ITG-Grace03s. The comparison among the other recent GRACE-based models yields bias as large as about 25-30%. The major discrepancies still occur for J_4, J_6 and J_8, which are just the zonals the combined LAGEOS/LAGOES II nodes are most sensitive to.

We explicitly worked out the orbital effects induced on the trajectory of a test particle by the the weak-field approximation of the Kerr-de Sitter metric. It results that the node, the pericentre and the mean anomaly undergo secular precessions proportional to k, which is a measure of the non linearity of the theory. We used such theoretical predictions and the latest observational determinations of the non-standard precessions of the perihelia of the inner planets of the Solar System to put a bound on k getting k <= 10^-29 m^-2. The node rate of the LAGEOS Earth's satellite yields k <= 10^-26 m^-2. The periastron precession of the double pulsar PSR J0737-3039A/B allows to obtain k <= 3 10^-21 m^-2. Interpreting k as a cosmological constant \Lambda, it turns out that such constraints are weaker than those obtained from the Schwarzschild-de Sitter metric.

In this note we mainly explore the possibility of measuring the action of the intrinsic gravitomagnetic field of the rotating Earth on the orbital motion of the Moon with the Lunar Laser Ranging (LLR) technique. Expected improvements in it should push the precision in measuring the Earth-Moon range to the mm level; the present-day Root-Mean-Square (RMS) accuracy in reconstructing the radial component of the lunar orbit is about 2 cm; its harmonic terms can be determined at the mm level. The current uncertainty in measuring the lunar precession rates is about 10^-1 milliarcseconds per year. The Lense-Thirring secular, i.e. averaged over one orbital period, precessions of the node and the perigee of the Moon induced by the Earth's spin angular momentum amount to 10^-3 milliarcseconds per year yielding transverse and normal shifts of 10^-1-10^-2 cm yr^-1. In the radial direction there is only a short-period, i.e. non-averaged over one orbital revolution, oscillation with an amplitude of 10^-5 m. Major limitations come also from some systematic errors induced by orbital perturbations of classical origin like, e.g., the secular precessions induced by the Sun and the oblateness of the Moon whose mismodelled parts are several times larger than the Lense-Thirring signal. The present analysis holds also for the Lue-Starkman perigee precession due to the multidimensional braneworld model by Dvali, Gabadadze and Porrati (DGP); indeed, it amounts to about 5 X 10^-3 milliarcseconds per year.

We investigate qualitatively and quantitatively the impact of the general relativistic gravito-electromagnetic forces on hyperbolic orbits around a massive spinning body. The gravito-magnetic field, which is the cause of the well known Lense-Thirring precessions of elliptic orbits, is generated by the spin S of the central body. It deflects and displaces the trajectories differently according to the mutual orientation of S and the orbital angular momentum L of the test particle. The gravito-electric force, which induces the Einstein precession of the perihelion of the orbit of Mercury, always deflects the trajectories inward irrespective of the L-S orientation. We numerically compute their effect on the range r, radial and transverse components v_r and v_\tau of the velocity and speed v of the NEAR spacecraft at its closest approach with the Earth in January 1998 when it experienced an anomalous increase of its asymptotic outgoing velocity v_\infty o of 13.46 +/- 0.01 mm sec^-1; while the gravito-electric force was modeled in the software used to process the NEAR data, this was not done for the gravito-magnetic one. The range-rate and the speed are affected by general relativistic gravito-electromagnetism at 10^-2 (gravito-electric)-10^-5 (gravito-magnetic) mm sec^-1 level. The changes in the range are of the order of 10^-2 (gravito-magnetic)-10^1 (gravito-electric) mm.

The astronomer E.V. Pitjeva, by analyzing with the EPM2008 ephemerides a large number of planetary observations including also two years (2004-2006) of normal points from the Cassini spacecraft, phenomenologically estimated a statistically significant non-zero correction to the usual Newtonian/Einsteinian secular precession of the longitude of the perihelion of Saturn, i.e. \Delta\dot\varpi_Sat = -0.006 +/- 0.002 arcsec/cy; the formal, statistical error is 0.0007 arcsec/cy. It can be explained neither by any of the standard classical and general relativistic dynamical effects mismodelled/unmodelled in the force models of the EPM2008 ephemerides nor by several exotic modifications of gravity recently put forth to accommodate certain cosmological/astrophysical observations without resorting to dark energy/dark matter. Both independent analyses by other teams of astronomers and further processing of larger data sets from Cassini will be helpful in clarifying the nature and the true existence of the anomalous precession of the perihelion of Saturn.

After the approval by the Italian Space Agency of the LARES satellite, which should be launched at the end of 2009 with a VEGA rocket and whose claimed goal is a about 1% measurement of the general relativistic gravitomagnetic Lense-Thirring effect in the gravitational field of the spinning Earth, it is of the utmost importance to reliably assess the total realistic accuracy that can be reached by such a mission. The observable is a linear combination of the nodes of the existing LAGEOS and LAGEOS II satellites and of LARES able to cancel out the impact of the first two even zonal harmonic coefficients of the multipolar expansion of the classical part of the terrestrial gravitational potential representing a major source of systematic error. While LAGEOS and LAGEOS II fly at altitudes of about 6000 km, LARES will be placed at an altitude of 1450 km. Thus, it will be sensitive to much more even zonals than LAGEOS and LAGEOS II. Their corrupting impact \delta\mu has been evaluated by using the standard Kaula's approach up to degree L=70 along with the sigmas of the covariance matrices of eight different global gravity solutions (EIGEN-GRACE02S, EIGEN-CG03C, GGM02S, GGM03S, JEM01-RL03B, ITG-Grace02s, ITG-Grace03, EGM2008) obtained by five institutions (GFZ, CSR, JPL, IGG, NGA) with different techniques from long data sets of the dedicated GRACE mission. It turns out \delta\mu about 100-1000% of the Lense-Thirring effect. An improvement of 2-3 orders of magnitude in the determination of the high degree even zonals would be required to constrain the bias to about 1-10%.

We deal with the attempts to measure the Lense-Thirring effect with the Satellite Laser Ranging (SLR) technique applied to the existing LAGEOS and LAGEOS II terrestrial satellites and to the recently approved LARES spacecraft.The first issue addressed here is: are the so far published evaluations of the systematic uncertainty induced by the bad knowledge of the even zonal harmonic coefficients J_L of the multipolar expansion of the Earth's geopotential reliable and realistic? Our answer is negative. Indeed, if the differences Delta J_L among the even zonals estimated in different Earth's gravity field global solutions from the dedicated GRACE mission are assumed for the uncertainties delta J_L instead of using their covariance sigmas sigma_JL, it turns out that the systematic uncertainty \delta\mu in the Lense-Thirring test with the nodes Omega of LAGEOS and LAGEOS II may be up to 3 to 4 times larger than in the evaluations so far published ($5-10%$) based on the use of the sigmas of one model at a time separately. The second issue consists of the possibility of using a different approach in extracting the relativistic signature of interest from the LAGEOS-type data. The third issue is the possibility of reaching a realistic total accuracy of 1% with LAGEOS, LAGEOS II and LARES, which should be launched in November 2009 with a VEGA rocket. While LAGEOS and LAGEOS II fly at altitudes of about 6000 km, LARES will be likely placed at an altitude of 1450 km. Thus, it will be sensitive to much more even zonals than LAGEOS and LAGEOS II. Their corrupting impact has been evaluated with the standard Kaula's approach up to degree L=60 by using Delta J_L and sigma_JL; it turns out that it may be as large as some tens percent.

We deal with a Yukawa-like long-range modified model of gravity (MOG) which recently allowed to successfully accommodate many astrophysical and cosmological features without resorting to dark matter. On Solar System scales MOG predicts retrograde secular precessions of the planetary longitudes of the perihelia \varpi whose existence has been put on the test here by taking the ratios of the observationally estimated Pitjeva's corrections to the standard Newtonian/Einsteinian perihelion precessions for different pairs of planets. It turns out that MOG, in the present form which turned out to be phenomenologically successful on astrophysical scales, is ruled out at more than 3sigma level in the Solar System. If and when other teams of astronomers will independently estimate their own extra-precessions of the perihelia it will be possible to repeat such a test.

We consider the recently estimated corrections \Delta\dot\varpi to the Newtonian/Einsteinian secular precessions of the longitudes of perihelia \varpi of several planets of the Solar System in order to evaluate whether they are compatible with the predicted precessions due to models of long-range modified gravity put forth to account for certain features of the rotation curves of galaxies without resorting to dark matter. In particular, we consider a logarithmic-type correction and a f(R) inspired power-law modification of the Newtonian gravitational potential. The results obtained by taking the ratio of the apsidal rates for different pairs of planets show that the modifications of the Newtonian potentials examined in this paper are not compatible with the secular extra-precessions of the perihelia of the Solar System's planets estimated by E.V. Pitjeva as solve-for parameters processing almost one century of data with the latest EPM ephemerides.

After the recent approval by the Italian Space Agency (ASI) of the LARES mission, which will be launched at the end of 2008 by a VEGA rocket to measure the general relativistic gravitomagnetic Lense-Thirring effect by combining LARES data with those of the existing LAGEOS and LAGEOS II satellites, it is of the utmost importance to assess if the claimed accuracy \lesssim 1% will be realistically obtainable. A major source of systematic error is the mismodelling \delta J_L in the static part of the even zonal harmonic coefficients J_L, L=2,4,6,.. of the multipolar expansion of the classical part of the terrestrial gravitational potential; such a bias crucially depends on the orbital configuration of LARES. If for \delta J_L the difference between the best estimates of different Earth's gravity solutions from the dedicated GRACE mission is conservatively taken instead of optimistically considering the statistical covariance sigmas of each model separately, as done so far in literature, it turns out that, since LARES will be likely launched in a low-orbit (semimajor axis a\lesssim 7600 km), the bias due to the geopotential may be up to ten times larger than what claimed, according to a calculation up to degree L=20. Taking into account also the even zonal harmonics with L>20, as required by the relatively low altitude of LARES, may further degrade the total accuracy. Should a nearly polar configuration (inclination to the Earth's equator i\approx 90 deg) be finally implemented, also other perturbations would come into play, further corrupting the measurement of the Lense-Thirring effect. The orbital configuration of LARES may also have some consequences in terms of non-gravitational perturbations and measurement errors.

E.V. Pitjeva, by processing more than 400,000 planetary observations of various types with the dynamical models of the EPM2006 ephemerides, recently estimated a correction to the canonical Newtonian-Einsteinian Venus' perihelion precession of -0.0004 +/- 0.0001 arcseconds per century. The prediction of general relativity for the Lense-Thirring precession of the perihelion of Venus is -0.0003 arcseconds per century. It turns out that neither other mismodelld/unmodelled standard Newtonian/Einsteinian effects nor exotic ones, postulated to, e.g., explain the Pioneer anomaly, may have caused the determined extra-precession of the Venus orbit which, thus, can be reasonably attributed to the gravitomagnetic field of the Sun, not modelled in the routines of the EPM2006 ephemerides. However, it must be noted that the quoted error is the formal, statistical one; the realistic uncertainty might be larger. Future improvements of the inner planets' ephemerides, with the inclusion of the Messenger and Venus-Express tracking data, should further improve the accuracy and the consistency of such a test of general relativity which would also benefit of the independent estimation of the extra-precessions of the perihelia (and the nodes) by other teams of astronomers.

In this paper we preliminarily explore the possibility of designing a dedicated satellite-based mission to measure the general relativistic gravitomagnetic Lense-Thirring effect in the gravitational field of Mars. The focus is on the systematic error induced by the multipolar expansion of the areopotential and on possible strategies to reduce it. It turns out that the major sources of bias are the Mars'equatorial radius R and the even zonal harmonics J_L, L = 2,4,6... of the areopotential. An optimal solution, in principle, consists of using two probes at high-altitudes (a\approx 9500-9600 km) and different inclinations, and suitably combining their nodes in order to entirely cancel out the bias due to \delta R. The remaining uncancelled mismodelled terms due to \delta J_L, L = 2,4,6,... would induce a bias \lesssim 1%, according to the present-day MGS95J gravity model, over a wide range of admissible values of the inclinations. The Lense-Thirring out-of-plane shifts of the two probes would amount to about 10 cm yr^-1.

In this paper we analyze the impact on the orbital motions of the outer planets of the solar system from Jupiter to Pluto of some velocity-dependent forces recently proposed to phenomenologically explain the Pioneer anomaly, and compare their predictions (secular variations of the longitude of perihelion \varpi or of the semimajor axis a and the eccentricity e) with the latest observational determinations by E.V. Pitjeva with the EPM2006 ephemerides. It turns out that while the predicted centennial shifts of a are so huge that they would have been easily detected for all planets with the exception of Neptune, the predicted anomalous precessions of \varpi are too small, with the exception of Jupiter, so that they are still compatible with the estimated corrections to the standard Newton-Einstein perihelion precessions. As a consequence, we incline to discard those extra-forces predicting secular variations of a and e, also for some other reasons, and to give a chance, at least observationally, to those models predicting still undetectable perihelion precessions. Of course, adequate theoretical foundations for them should be found.

The fast spinning B-star Regulus has recently been found to be orbited by a fainter companion in a close circular path with orbital period P_b = 40.11(2) d. Being its equatorial radius R_e 32% larger than the polar one R_p, Regulus possesses a remarkable quadrupole mass moment Q. We investigate the effects of Q on the orbital period P_b of its companion in order to see if they are measurable, given the present-day level of accuracy in measuring P_b. Conversely, we will look for deviations from the third Kepler law, attributed to the quadrupole mass moment Q of Regulus, to constrain the ratio \gamma=m/M of the system's masses.

We consider the double pulsar binary system as a laboratory to locally test the orbital effects induced by an uniform cosmological constant $\Lambda$ in the framework of the known general relativistic laws of gravity, and the DGP braneworld model of gravity independently of the cosmological acceleration itself for which they were introduced. We, first, construct the ratio R=Delta\dot\omega/Delta P of the discrepancies between the phenomenologically determined periastron rate \dot\omega and orbital period and their predicted values from the 1PN approximation and the third Kepler law. Then, we compare its value R = (0.3 +/- 4) \times 10^-11 s^-2, compatible with zero within the errors, to the ratios R_Lambda and R_DGP of the effects induced on the apsidal rate and the orbital period by Lambda and the DGP gravity; we find them neatly incompatible with R being R_Lambda = (3.4 +/- 0.3) \times 10^-8 s^-2 and R_DGP = (1.4 +/- 0.1) \times 10^-7 s^-2, respectively. Such a result, which for the case of Lambda is valid also for any other Hooke-like exotic force proportional to r, is in agreement with other negative local tests recently performed in the Solar System with the ratios of the non-Newtonian/Einsteinian perihelion precessions for several pairs of planets.

In this Chapter we deal with the attempts to measure the general relativistic gravitomagnetic Lense-Thirring effect with the Satellite Laser Ranging (SLR) technique applied to the existing LAGEOS and LAGEOS II terrestrial satellites and to the recently approved LARES.

In this paper we explore the consequences of the recent determination of the mass m=(8.7 +/- 0.8)M_Sun of Cygnus X-1, obtained from the Quasi-Periodic Oscillation (QPO)-photon index correlation scaling, on the orbital and physical properties of the binary system HDE 226868/Cygnus X-1. By using such a result and the latest spectroscopic optical data of the HDE 226868 supergiant star we get M=(24 +/- 5)M_Sun for its mass. It turns out that deviations from the third Kepler law significant at more than 1-sigma level would occur if the inclination i of the system's orbital plane to the plane of the sky falls outside the range 41-56 deg: such deviations cannot be due to the first post-Newtonian (1PN) correction to the orbital period because of its smallness; interpreted in the framework of the Newtonian theory of gravitation as due to the stellar quadrupole mass moment Q, they are unphysical because Q would take unreasonably large values. By conservatively assuming that the third Kepler law is an adequate model for the orbital period we obtain i=(48 +/- 7) deg which yields for the relative semimajor axis a=(42 +/- 9)R_Sun. Our estimate for the Roche's lobe of HDE 226868 is r_M = (21 +/- 6)R_Sun.

Here we consider the possibility of measuring at 10% or better the moment of inertia I of the J0737-3039A via the gravitomagnetic spin-orbit periastron precession. It turns out that the prospect of measuring the moment of inertia of PSR J0737-3039A at 10% accuracy or better seems unlikely given the limitations to the precision with which the system's basic binary and post-Keplerian parameters can be measured via radio timing.

We use the corrections to the Newton-Einstein secular precessions of the longitudes of perihelia of some planets (Mercury, Earth, Mars, Jupiter, Saturn) of the Solar System, phenomenologically estimated as solve-for parameters by the Russian astronomer E.V. Pitjeva in a global fit of almost one century of data with the EPM2004 ephemerides, in order to put on the test the expression for the perihelion precession induced by an uniform cosmological constant Lambda in the framework of the Schwarzschild-de Sitter (or Kottler) space-time. We compare such an extra-rate to the estimated corrections to the planetary perihelion precessions by taking their ratio for different pairs of planets instead of using one perihelion at a time for each planet separately, as done so far in literature. The answer is neatly negative, even by further re-scaling by a factor 10 (and even 100 for Saturn) the errors in the estimated extra-precessions of the perihelia released by Pitjeva. Our conclusions hold also for any other metric perturbation having the same dependence on the spatial coordinates, as those induced by other general relativistic cosmological scenarios and by many modified models of gravity. Currently ongoing and planned interplanetary spacecraft-based missions should improve our knowledge of the planets' orbits allowing for more stringent constraints.

In this letter we investigate the deep Newtonian regime of the MOND paradigm from a purely phenomenological point of view by exploiting the least-square estimated corrections to the secular rates of the perihelia of the inner and of some of the outer planets of the Solar System by E.V. Pitjeva with the EPM2004 ephemerides. By using $\mu(x)\approx 1-k_0(1/x)^n$ for the interpolating MONDian function, and by assuming that $k_0$, considered body-independent so to avoid violations of the equivalence principle, experiences no spatial variations throughout the Solar System we tightly constrain $n$ with the ratios of the perihelion precessions for different pairs of planets. We find that the range $1\leq n\leq 2$ is neatly excluded at much more than $3-\sigma$ level. Such a test would greatly benefit from the use of extra-precessions of perihelia independently estimated by other groups as well.

In this paper we determine the tidal distortion parameter k_m of the secondary partner (mass loser) of the semi-detached eclipsing binary system V621 Cen by comparing the phenomenologically determined orbital period P_b=3.683549(11) d to the Keplerian one P^Kep computed with the values of the relevant system's parameters determined independently of the third Kepler law itself. Our result is k_m = -1.5 +/- 0.6. Using the periastron precession, as traditionally done with other eclipsing binaries in eccentric orbits, would have not been possible because of the circularity of the V621 Cen path.

In this paper we use the ratio of the corrections to the standard Newtonian/Einsteinian secular precessions of the longitudes of perihelia of Jupiter and Saturn, recently estimated by the Russian astronomer E.V. Pitjeva by fitting almost one century of data with the EPM ephemerides, to make an independent, planetary-based test of the hypothesis that the Pioneer anomaly (PA), as it is presently known in the 5-10 AU region, is of gravitational origin. Accounting for the errors in the determined apsidal extra-rates and in the values of the PA acceleration at the orbits of Jupiter and Saturn the answer is negative. If and when the re-analysis of the entire Pioneer 10/11 will be completed more firm conclusions could be reached. Moreover, it would also be important that other teams of astronomers estimate independently their own corrections to the perihelion precessions.

The current bounds on the PPN parameters gamma and beta are of the order of 10^-4-10^-5. Various missions aimed at improving such limits by several orders of magnitude have more or less recently been proposed like LATOR, ASTROD, BepiColombo and GAIA. They involve the use of various spacecraft, to be launched along interplanetary trajectories, for measuring the effects of the solar gravity on the propagation of electromagnetic waves. In this paper we investigate what is needed to measure the combination nu=(2+2gamma-beta)/3 of the post-Newtonian gravitoelectric Einstein perigee precession of a test particle to an accuracy of about 10^-5 with a pair of drag-free spacecraft in the Earth's gravitational field. It turns out that the latest gravity models from the dedicated CHAMP and GRACE missions would allow to reduce the systematic error of gravitational origin just to this demanding level of accuracy. In regard to the non-gravitational errors, the spectral noise density of the drag-free sensors required to reach such level of accuracy would amounts to 10^-8-10^-9 cm s^-2 Hz^-1/2 over very low frequencies. Although not yet obtainable with the present technologies, such level of compensation is much less demanding than those required for, e.g., LISA. As a by-product, an independent measurement of the post-Newtonian gravitomagnetic Lense-Thirring effect with a 0.9% accuracy would be possible as well. The forthcoming Earth gravity models from CHAMP and GRACE will further reduce the systematic gravitational errors in both of such tests.

The post-Newtonian gravitoelectric secular rate of the mean anomaly M is worked out for a two-body system in the framework of the General Theory of Relativity. The possibility of using such an effect, which is different from the well known decrease of the orbital period due to gravitational wave emission, as a further post-Keplerian parameter in binary systems including one pulsar is examined. The resulting effect is almost three times larger than the periastron advance \dot\omega. E.g., for the recently discovered double pulsar system PSR J0737-3039 A+B it would amount to -47.79 deg yr^-1. This implies that it could be extracted from the linear part of a quadratic fit of the orbital phase because the uncertainties both in the linear drift due to the mean motion and in the quadratic shift due to the gravitational wave are smaller. The availability of such additional post-Keplerian parameter would be helpful in further constraining the General Theory of Relativity, especially for such systems in which some of the other post-Keplerian parameters can be measured with limited accuracy. Moreover, also certain pulsar-white dwarf binary systems, characterized by circular orbits like PSR B1855+09 and a limited number of measured post-Keplerian parameters, could be used for constraining competing theories of gravity.

In this paper we use the corrections to the usual Newton-Einstein secular precessions of the perihelia of the inner planets of the Solar System, phenomenologically estimated as solve-for parameters by the Russian astronomer E.V. Pitjeva by fitting almost one century of data with the EPM2004 ephemerides, in order to constrain some long-range models of modified gravity recently put forth to address the dark energy and dark matter problems. The models examined here are the four-dimensional ones obtained with the addition of inverse powers and logarithm of some curvature invariants, and the multidimensional braneworld model by Dvali, Gabadadze and Porrati (DGP). After working out the analytical expressions of the secular perihelion precessions induced by the corrections to the Newtonian potential of such models, we compare them to the estimated corrections to the rates of perihelia by taking their ratio for different pairs of planets instead of using one perihelion at a time for each planet separately, as done so far in literature. As a result, the curvature invariants-based models are ruled out, even by re-scaling by a factor 10 the errors in the planetary orbital parameters estimated by Pitjeva. Less neat is the situation for the DGP model. Only the general relativistic Lense-Thirring effect, not included, as the other exotic models considered here, by Pitjeva in the dynamical force models used in the estimation process, passes such a test. It would be important to repeat the present analysis by using corrections to the precessions of perihelia independently estimated by other teams of astronomers as well, but, at present, such rates are not yet available.

It has often been claimed that the proposed Earth artificial satellite LARES/WEBER-SAT-whose primary goal is, in fact, the measurement of the general relativistic Lense-Thirring effect at a some percent level-would allow to greatly improve, among (many) other things, the present-day (10^-13) level of accuracy in testing the equivalence principle as well. Recent claims point towards even two orders of magnitude better, i.e. 10^-15. In this note we show that such a goal is, in fact, unattainable by many orders of magnitude being, instead, the achievable level of the order of 10^-9.

In this paper we use the latest corrections to the Newton-Einstein secular perihelion rates of some planets of the Solar System, phenomenologically estimated with the EPM2004 ephemerides by the Russian astronomer E.V. Pitjeva, to put severe constraints on the range parameter lambda characterizing the Yukawa-like modifications of the Newtonian inverse-square law of gravitation. It turns out that the range cannot exceed about one tenth of an Astronomical Unit. We assumed neither equivalence principle violating effects nor spatial variations of $\alpha$ and $\lambda$. This finding may have important consequences on all the modified theories of gravity involving Yukawa-type terms with range parameters much larger than the Solar System size. However, caution is advised since we, currently have at our disposal only the periehlion extra-rates estimated by Pitjeva: if and when other groups will estimate their own corrections to the secular motion of perihelia, more robust and firm tests may be conducted.

In this paper I deal with the WD0137-349 binary system consisting of a white dwarf (WD) and a brown dwarf (BD) in a close circular orbit of about 116 min. I, first, constrain the admissible range of values for the inclination i by noting that, from looking for deviations from the third Kepler law, the quadrupole mass moment Q would assume unlikely large values, incompatible with zero at more than 1-sigma level for i < 35 deg and i > 43 deg. Then, by conservatively assuming that the most likely values for i are those that prevent such an anomalous behavior of Q, i.e. those for which the third Kepler law is an adequate modeling of the orbital period, I obtain i=39 +/- 2 deg. Such a result is incompatible with the value i=35 deg quoted in literature by more than 2 sigma. Conversely, it is shown that the white dwarf's mass range obtained from spectroscopic measurements is compatible with my experimental range, but not for i=35 deg. As a consequence, my estimate of $i$ yields an orbital separation of a=(0.59 +/- 0.05)R_Sun and an equilibrium temperature of BD of T_eq=(2087 +/- 154)K which differ by 10% and 4%, respectively, from the corresponding values for i=35 deg.

In this paper we use Delta P = -1.772341 +/- 13.153788 s between the phenomenologically determined orbital period P_b of the PSR J0737-3039 double pulsar system and the purely Keplerian period P^(0)=2\pi\sqrt{a^3/G(m_A+m_B)} calculated with the system's parameters, determined independently of the third Kepler law itself, in order to put constraints on some models of modified gravity (f(R), Yukawa-like fifth force, MOND). The major source of error affecting Delta P is not the one in the phenomenologically measured period (\delta P_b=4 10^-6 s), but the systematic uncertainty \delta P^(0) in the computed Keplerian one due to the relative semimajor axis a mainly caused, in turn, by the errors in the ratio R of the pulsars' masses and in sin i. We get |\kappa|< 0.8 10^-26 m^-2 for the parameter that in the f(R) framework is a measure of the non linearity of the theory, |\alpha|< 5.5 10^-4 for the fifth-force strength parameter (for \lambda\approx a=0.006 AU). The effects predicted by the strong-acceleration regime of MOND are far too small to be constrained with some effectiveness today and in the future as well. In view of the continuous timing of such an important system, it might happen that in the near future it will be possible to obtain somewhat tighter constraints.

We preliminarily investigate the impact of the Kuiper Belt Objects (KBOs) and of the asteroid ring on some proposed high-precision tests of Newtonian and post-Newtonian gravity to be performed in the Solar System by means of spacecraft in heliocentric \approx 1 AU orbits and accurate orbit determination of some of the inner planets. It turns out that the Classical KBOSs (CKBOS), which amount to \approx 70% of the observed population of Trans-Neptunian bodies, induce a systematic secular error of about 1 m after one year in the transverse direction T of the orbit of a test particle orbiting at 1 AU from the Sun. For Mercury the ratios of the secular perihelion precessions induced by CKBOs to the ones induced by the general relativity and the solar oblateness J_2 amount to 6 10^-7 and 8 10^-4, respectively. The secular transverse perturbation induced on a \approx 1 AU orbit by the asteroid ring, which globally accounts for the action of the minor asteroids whose mass is about 5 10^-10 solar masses, is 10 m yr^-1; the bias on the relativistic and J_2 Mercury perihelion precessions is 6.1 10^-6 and 1 10^-2, respectively. Given the very ambitious goals of many expensive and complex missions aimed to testing gravitational theories to unprecedented levels of accuracy, these notes may suggest further and more accurate investigations of such sources of potentially insidious systematic bias.

In this paper we propose an interpretation of the time series of the RMS overlap differences of the out-of-plane part $N$ of the orbit of the Martian polar orbiter Mars Global Surveyor (MGS) over a time span Delta P of about 5 years (14 November 1999-14 January 2005) in terms of the general relativistic gravitomagnetic Lense-Thirring effect. It turns out that the average of such a time series over Delta P, normalized to the predicted Lense-Thirring out-of-plane mean shift over the same time span, is \mu= 1.0018 +/- 0.0053. A linear fit of the considered time series over Delta P yields a normalized slope 1.1 +/- 0.3 strengthening our interpretation of the determined MGS out-of-plane features of motion as due to the gravitomagnetic force of Mars.

The general relativistic Lense-Thirring precessions of the perihelia of the inner planets of the Solar System are about 10^-3 arcseconds per century. Recent improvements in planetary orbit determination may yield the first observational evidence of such a tiny effect. Indeed, corrections to the known perihelion rates of -0.0036 +/- 0.0050, -0.0002 +/- 0.0004 and 0.0001 +/- 0.0005 arcseconds per century were recently estimated by E.V. Pitjeva for Mercury, the Earth and Mars, respectively, on the basis of the EPM2004 ephemerides and a set of more than 317,000 observations of various kinds. The predicted relativistic Lense-Thirring precessions for these planets are -0.0020, -0.0001 and -3 10^-5 arcseconds per century, respectively and are compatible with the determined perihelia corrections. The relativistic predictions fit better than the zero-effect hypothesis, especially if a suitable linear combination of the perihelia of Mercury and the Earth, which a priori cancels out any possible bias due to the solar quadrupole mass moment, is considered. However, the experimental errors are still large. Also the latest data for Mercury processed independently by Fienga et al. with the INPOP ephemerides yield preliminary insights about the existence of the solar Lense-Thirring effect. The data from the forthcoming planetary mission BepiColombo will improve our knowledge of the orbital motion of this planet and, consequently, the precision of the measurement of the Lense-Thirring effect. As a by-product of the present analysis, it is also possible to constrain the strength of a Yukawa-like fifth force to a 10^-12-10^-13 level at scales of about one Astronomical Unit (10^11 m).

In this note we comment on a recent paper by I.Ciufolini about the possibility of placing the proposed terrestrial satellite LARES in a low-altitude, nearly polar orbit in order to measure the general relativistic Lense-Thirring effect with its node. Ciufolini claims that, for a departure of 4 deg in the satellite's inclination $i$ from the ideal polar configuration (i=90 deg), the impact of the errors in the even zonal harmonics of the geopotential, modelled with EIGEN-GRACE02S, would be nearly zero allowing for a few-percent measurement of the Lense-Thirring effect. Instead, we find that, with the same Earth gravity model and for the same values of the inclination, the upper bound of the systematic error due to the mismodelling in the even zonals amounts to 64% of the relativistic effect investigated.

In this letter we dynamically determine the quadrupole mass moment Q of the magnetic white dwarf WD 0137-349 by analyzing the period of the recently discovered brown dwarf moving around it in a close 2-hr orbit. It turns out that a purely Newtonian model for the orbit of WD 0137-349B, assumed circular and equatorial, is adequate, given the present-day accuracy in knowing the orbital parameters of such a binary system. Our result is Q=(-1.4615 +/- 0.9004) 10^47 kg m^2 for i=35 deg. It is able to accommodate the 3-sigma significant discrepancy of (1.0 +/- 0.3) 10^-8 s^-2 between the inverse square of the phenomenologically determined orbital period and the inverse square of the calculated Keplerian one. The impact of i, for which a range Delta i of possible values close to 35 deg is considered, is investigated as well.

In this note we reply to the criticisms by Felici concerning some aspects of the recent frame-dragging test performed by Iorio with the Mars Global Surveyor (MGS) spacecraft in the gravitational field of Mars.

In order to satisfy the equivalence principle, any non-conventional mechanism proposed to gravitationally explain the Pioneer anomaly, in the form in which it is presently known from the so-far analyzed Pioneer 10/11 data, cannot leave out of consideration its impact on the motion of the planets of the Solar System as well, especially those orbiting in the regions in which the anomalous behavior of the Pioneer probes manifested itself. In this paper we, first, discuss the residuals of the right ascension \alpha and declination \delta of Uranus, Neptune and Pluto obtained by processing various data sets with different, well established dynamical theories (JPL DE, IAA EPM, VSOP). Second, we use the latest determinations of the perihelion secular advances of some planets in order to put on the test two gravitational mechanisms recently proposed to accommodate the Pioneer anomaly based on two models of modified gravity. Finally, we adopt the ranging data to Voyager 2 when it encountered Uranus and Neptune to perform a further, independent test of the hypothesis that a Pioneer-like acceleration can also affect the motion of the outer planets of the Solar System. The obtained answers are negative.

We report on a \lessim 1% test of the Lense-Thirring effect with the Mars Global Surveyor (MGS) spacecraft and on certain features of motion of Uranus, Neptune and Pluto which contradict the hypothesis that the Pioneer anomaly can be caused by some gravitational mechanism.

In this note we reply to the criticisms by Sindoni, Paris and Ialongo concerning some aspects of the recent frame-dragging test performed by Iorio with the Mars Global Surveyor (MGS) spacecraft in the gravitational field of Mars.

Recently, Brownstein and Moffat proposed a gravitational mechanism to explain the Pioneer anomaly based on their scalar-tensor-vector (STVG) metric theory of gravity. In this paper we show that their model, fitted to the presently available data for the anomalous Pioneer 10/11 acceleration, is in contrast with the latest determinations of the perihelion extra-rates of Jupiter, Saturn and Uranus.

In this talk we review the perspectives of testing the multidimensional Dvali-Gabadadze-Porrati (DGP) model of modified gravity in the Solar System. The inner planets, contrary to the giant gaseous ones, yield the most promising scenario for the near future.

In this paper we study the effects of $f(R)$ Theories of Gravity on Solar System gravitational tests. In particular, starting from an exact solution of the field equation in vacuum, in the Palatini formalism, we work out the effects that the modifications to the Newtonian potential would induce on the Keplerian orbital elements of the Solar System planets, and compare them with the latest results in planetary orbit determination from the EPM2004 ephemerides. It turns out that the longitudes of perihelia and the mean longitudes are affected by secular precessions. We obtain the resulting best estimate of the parameter $k$ which, being simply related to the scalar curvature, measures the non linearity of the gravitational theory. We use our results to constrain the cosmological constant and show how $f(R)$ functions can be constrained, in principle. What we obtain suggests that, in agreement with other recent papers, the Solar System experiments are not effective to set such constraints, if compared to the cosmologically relevant values.

The general relativistic gravitomagnetic clock effect involves a coupling between the orbital motion of a test particle and the rotation of the central mass and results in a difference in the proper periods of two counter-revolving satellites. It is shown that at O(c^-2) this effect has a simple analogue in the electromagnetic case. Moreover, in view of a possible measurement of the clock effect in the gravitational field of the Earth, we investigate the influence of some classical perturbing forces of the terrestrial space environment on the orbital motion of test bodies along opposite trajectories.

In this paper we put on the test the new mechanism of gravitational origin recently put forth by Jaekel and Reynaud in order to explain the Pioneer anomaly in the framework of their post-Einsteinian metric extension of general relativity. According to such a proposal, the secular part of the anomalous acceleration experienced by the twin spacecraft of about 1 nm s^-2 could be caused by an extra-potential \delta\Phi_P=c^2\chi r^2, with \chi=4 10^-8 AU^-2, coming from the second sector of the considered model. When applied to the motion of the planets of the Solar System, it would induce anomalous secular perihelion advances which amount to tens-hundreds of arcseconds per century for the outer planets. As for other previously proposed non-conventional gravitational explanations of the Pioneer anomaly, the answer of the latest observational determinations of the residual perihelion rates by RAS IAA is neatly and unambiguously negative. The presence of another possible candidate to explain the Pioneer anomaly, i.e. the extra-potential \delta\Phi_N, linear in distance, from the first sector of the Jaekel and Reynaud model, is ruled out not only by the residuals of the optical data of the outer planets processed with the recent RAS IAA EPM2004 ephemerides but also by those produced with other, older dynamical theories like, e.g., the well known NASA JPL DE200.

The aim of this paper is to dynamically constrain the quadrupole mass moment J_2 of the main-sequence HD 209458 star. The adopted method is the confrontation between the measured orbital period of its transiting planet Osiris and a model of it. Osiris is assumed to move along an equatorial and circular orbit. Our estimate, for given values of the stellar mass M and radius R and by assuming the validity of general relativity, is J_2=(3.5 +/- 385.1) 10^-5. Previous fiducial evaluations based on indirect, spectroscopic measurements yielded J_2\sim 10^-6: such a value is compatible with our result.

In this paper we dynamically determine the mass of the Kuiper Belt Objects by exploiting the latest observational determinations of the orbital motions of the inner planets of the Solar System. Our result, in units of terrestrial masses, is 0.033 +/- 0.115 by modelling the Classical Kuiper Belt Objects as an ecliptic ring of finite thickness. A two-rings model yields for the Resonant Kuiper Belt Objects a value of 0.018 +/- 0.063. Such figures are consistent with recent determinations obtained with ground and space-based optical techniques. Some implications for precise tests of Einsteinian and post-Einsteinian gravity are briefly discussed.

Recently, Jaekel and Reynaud put forth a metric linear extension of general relativity which, in the intentions of its proponents, would be able, among other things, to provide a gravitational mechanism for explaining the Pioneer anomaly without contradicting either the equivalence principle or what we know about the planetary motions. In this paper we perform an independent test of such an hypothesis by showing that the planets' orbits are, in fact, affected by the suggested mechanism as well and comparing the resulting effects with the latest observational determinations. It turns out that the predicted perihelion precessions, expressed in terms of an adjustable free parameter zeta_P set equal to the value used to reproduce the magnitude of the Pioneer anomalous acceleration, are quite different from the observationally determined extra-advances of such Keplerian element for the inner planets. Conversely, the values obtained for zeta_P from the determined perihelion extra-rates of the inner planets turn out to be in disagreement with the value which would be required to accommodate the Pioneer anomaly. As a consequence, the suggested explanation of gravitational origin for the Pioneer anomaly, based on the assumption that zeta_P is constant throughout the Solar System, should be rejected, at least in its present form.

Here we analyze in detail some aspects of the proposed use of Ajisai and Jason-1, together with the LAGEOS satellites, to measure the general relativistic Lense-Thirring effect in the gravitational field of the Earth. A linear combination of the nodes of such satellites is the proposed observable. The systematic error due to the mismodelling in the uncancelled even zonal harmonics would be \sim 1% according to the latest present-day CHAMP/GRACE-based Earth gravity models. In regard to the non-gravitational perturbations especially affecting Jason-1, only relatively high-frequency harmonic perturbations should occur: neither semisecular nor secular bias of non-gravitational origin should affect the proposed combination: their maximum impact is evaluated to \sim 4% over 2 years. Our estimation of the root-sum-square total error is about 4-5% over at least 3 years of data analysis required to average out the uncancelled tidal perturbations.

In this paper we address the following question: do the recent advances in the orbit determination of the major natural satellites of Saturn obtained with the analysis of the first data sets from the Cassini mission allow to detect the general relativistic gravitoelectric orbital precessions of such moons? The answer is still negative. The present-day down-track accuracy would be adequate for Mimas, Enceladus, Thetys, Dione, Rhea and Titan and inadequate for Hyperion, Iapetus and Phoebe. Instead, the size of the systematic errors induced by the mismodelling in the key parameters of the Saturnian gravitational field like the even zonal harmonics Jl are larger than the relativistic down-track shifts by about one order of magnitude, mainly for the inner satellites like Mimas, Enceladus, Thetys, Dione, Rhea, Titan and Hyperion. Iapetus and Phoebe are not sensibly affected by such kind of perturbations. Moreover, the bias due to the uncertainty in Saturn's GM is larger than the relativistic down-track effects for all such moons. Proposed linear combinations of the satellites' orbital elements would allow to cancel out the impact of the mismodelling in the low-degree even zonal harmonics and GM, but the combined down-track errors would be larger than the combined relativistic signature.

In this paper we reply to recent claims by Ciufolini and Pavlis about certain aspects of the measurement of the general relativistic Lense-Thirring effect in the gravitational field of the Earth. I) The proposal by such authors of using the existing satellites endowed with some active mechanism of compensation of the non-gravitational perturbations as an alternative strategy to improve the currently ongoing Lense-Thirring tests is unfeasible because of the impact of the uncancelled even zonal harmonics of the geopotential and of some time-dependent tidal perturbations. II) It is shown that their criticisms about the possibility of using the existing altimeter Jason-1 and laser-ranged Ajisai satellites are groundless.III) Ciufolini and Pavlis also claimed that we would have explicitly proposed to use the mean anomaly of the LAGEOS satellites in order to improve the accuracy of the Lense-Thirrring tests. We prove that it is false. In regard to the mean anomaly of the LAGEOS satellites, Ciufolini himself did use such an orbital element in some previously published tests. About the latest test performed with the LAGEOS satellites, IV) we discuss the cross-coupling between the inclination errors and the first even zonal harmonic as another possible source of systematic error affecting it with an additional 9% bias. V) Finally, we stress the weak points of the claims about the origin of the two-nodes LAGEOS-LAGEOS II combination used in that test.

In this paper we test the hypothesis that the Pioneer anomaly can be of gravitational origin by comparing the predicted model-independent shifts Delta a/a for the semimajor axis of Uranus and Neptune with the Voyager 2 radio-technical distance measurements performed at JPL-NASA. As in the case of other tests based on different methods and data sets (secular perihelion advance, right ascension/declination residuals over about one century), the orbits of the investigated planets are not affected by any anomalous acceleration like that experienced by the Pioneer 10/11 spacecraft.

The multidimensional braneworld gravity model by Dvali, Gabadadze and Porrati was primarily put forth to explain the observed acceleration of the expansion of the Universe without resorting to dark energy. One of the most intriguing features of such a model is that it also predicts small effects on the orbital motion of test particles which could be tested in such a way that local measurements at Solar System scales would allow to get information on the global properties of the Universe. Lue and Starkman derived a secular extra-perihelion \omega precession of 5\times 10^-4 arcseconds per century, while Iorio showed that the mean longitude \lambda is affected by a secular precession of about 10^-3 arcseconds per century. Such effects depend only on the eccentricities e of the orbits via second-order terms: they are, instead, independent of their semimajor axes a. Up to now, the observational efforts focused on the dynamics of the inner planets of the Solar System whose orbits are the best known via radar ranging. Since the competing Newtonian and Einsteinian effects like the precessions due to the solar quadrupole mass moment J2, the gravitoelectric and gravitomagnetic part of the equations of motion reduce with increasing distances, it would be possible to argue that an analysis of the orbital dynamics of the outer planets of the Solar System, with particular emphasis on Saturn because of the ongoing Cassini mission with its precision ranging instrumentation, could be helpful in evidencing the predicted new features of motion. In this note we investigate this possibility in view of the latest results in the planetary ephemeris field. Unfortunately, the current level of accuracy rules out this appealing possibility and it appears unlikely that Cassini and GAIA will ameliorate the situation.

We discuss the impact of the present-day uncertainties in the recently released CHAMP and/or GRACE Earth gravity models on the measurement of the Lense-Thirring effect with the nodes of the LAGEOS satellites. Also the role of the secular variations Jdots of the even zonal harmonics is quantitatively assessed via numerical simulations and tests. While the systematic error due to the static part of the geopotential ranges from 4% (EIGEN-GRACE02S) to 9% (GGM02S), the impact of the Jdots amounts to 13% over 11 years. This yields a 19% 1-sigma total error in the test recently performed with EIGEN-GRACE02S.

According to general relativity, the gravitomagnetic Lense-Thirring force of Mars would secularly shift the orbital plane of the Mars Global Surveyor (MGS) spacecraft by an amount of 1.5 m, on average, in the cross-track direction over 5 years. The determined cross-track post-fit residuals of MGS, built up by neglecting just the gravitomagnetic force in the dynamical force models and without fitting any empirical cross-track acceleration which could remove the relativistic signal, amount to 1.6 m, on average, over a 5-years time interval spanning from 10 February 2000 to 14 January 2005. The discrepancy with the predictions of general relativity is 6%.

We study Doppler effects in curved space-time, i.e. the frequency shifts induced on electromagnetic signals propagating in the gravitational field. In particular, we focus on the frequency shift due to the bending of light rays in weak gravitational fields. We consider, using the PPN formalism, the gravitational field of an axially symmetric distribution of mass. The zeroth order, i.e. the sphere, is studied then passing to the contribution of the quadrupole moment, and finally to the case of a rotating source. We give numerical estimates for situations of physical interest, and by a very preliminary analysis, we argue that analyzing the Doppler effect could lead, in principle, in the foreseeable future, to the measurement of the quadrupole moment of the giant planets of the Solar System.

In this paper we very preliminarily investigate the possibility of measuring the post-Newtonian general relativistic gravitoelectric and gravitomagnetic components of the acceleration of gravity on the Earth, in continuous regime, with two absolute measurements at the equator and the south pole with superconducting gravimeters. The magnitudes of such relativistic effects are 10^-10 m s^-2 and 10^-11 m s^-2, respectively. Unfortunately, the present-day uncertainties in the Earth's geodetic parameters which enter the classical Newtonian terms induce systematic errors 1-2 orders of magnitude larger than the relativistic ones. Moreover, a \sim 1 ngal sensitivity can be reached by the currently available superconducting gravimeters, but only for relative measurements.

In this paper we investigate the effects that an anomalous acceleration as that experienced by the Pioneer spacecraft after they passed the 20 AU threshold would induce on the orbital motions of the Solar System planets placed at heliocentric distances of 20 AU or larger as Uranus, Neptune and Pluto. It turns out that such an acceleration, with a magnitude of 8.74\times 10^-10 m s^-2, would affect their orbits with secular and short-period signals large enough to be detected according to the latest published results by E.V. Pitjeva, even by considering errors up to 30 times larger than those released. The absence of such anomalous signatures in the latest data rules out the possibility that in the region 20-40 AU of the Solar System an anomalous force field inducing a constant and radial acceleration with those characteristics affects the motion of the major planets.

In this paper we explicitly work out the effects that a spherically symmetric distribution of dark matter with constant density would induce on the Keplerian orbital elements of the Solar System planets and compare them with the latest results in planetary orbit determination from the EPM2004 ephemerides. It turns out that the longitudes of perihelia and the mean longitudes are affected by secular precessions. The resulting upper bounds on dark matter density, obtained from the EPM2004 formal errors in the determined mean longitude shifts over 90 years, lie in the range 10^-19-10^-20 g cm^-3 with a peak of 10^-22 g cm^-3 for Mars. Suitable combinations of the planetary mean longitudes and perihelia, which cancel out the aliasing impact of some of the unmodelled or mismodelled forces of the dynamical models of EPM2004, yield a global upper bound of 7 10^-20 g cm^-3 and 4 10^-19 g cm^-3, respectively.

In this paper we investigate the possibility of measuring the general relativistic gravitoelectric contribution P^(GE) to the orbital period P of the transiting exoplanet HD 209458b 'Osiris'. It turns out that the predicted magnitude of such an effect is \sim 0.1 s, while the most recent determinations of the orbital period of HD 209458b with the photometric transit method are accurate to \sim 0.01 s. The present analysis shows that the major limiting factor is the \sim 1 m s^-1 sensitivity in the measurement of the projected semiamplitude of the star's radial velocity K. Indeed, it affects the determination of the mass m of the planet which, in turn, induces a systematic error in the Keplerian period P^(0) of \sim 8 s. It is of crucial importance because P^(0) should be subtracted from the measured period in order to extract the relativistic correction. The present-day uncertainty in $m$ does not yet make necessary the inclusion of relativistic corrections in the data-reduction process of the determination of the system's parameters. The present situation could change only if improvements of one-two orders of magnitude in the ground-based Doppler spectroscopy technique occurred.

In addition to the pericentre \omega, the mean anomaly M and, thus, the mean longitude \lambda, also the orbital period Pb and the mean motion $n$ of a test particle are modified by the Dvali-Gabadadze-Porrati gravity. While the correction to Pb depends on the mass of the central body and on the geometrical features of the orbital motion around it, the correction to $n$ is independent of them, up to terms of second order in the eccentricity $e$. The latter one amounts to about 2\times 10^-3 arcseconds per century. The present-day accuracy in determining the mean motions of the inner planets of the Solar System from radar ranging and differential Very Long Baseline Interferometry is 10^-2-5\times 10^-3 arcseconds per century, but it should be improved in the near future when the data from the spacecraft to Mercury and Venus will be available.

In this note we show that the latest determinations of the residual Mercury's perihelion advance, obtained by accounting for almost all known Newtonian and post-Newtonian orbital effects, yields only very broad constraints on the cosmological constant. Indeed, from \delta\dot\omega=-0.0036 + - 0.0050 arcseconds per century one gets -2 10^-34 km^-2 < Lambda < 4 10^-35 km^-2. The currently accepted value for Lambda, obtained from many independent cosmological and large-scale measurements, amounts to almost 10^-46 km^-2.

In this letter we work out the secular precession of the spin of a gyroscope in geodesic motion around a central mass in the framework of the Dvali-Gabadadze-Porrati multidimensional gravity model. Such an effect, which depends on the mass of the central body and on the orbit radius of the gyroscope, contrary to the precessions of the orbital elements of the orbit of a test body, is far too small to be detected.

The impact of the latest combined CHAMP/GRACE/terrestrial measurements Earth gravity model EIGEN-CG03C on the measurement of the Lense-Thirring effect with some linear combinations of the nodes of some of the existing Earth's artificial satellites is presented. The 1-sigma upper bound of the systematic error in the node-node LAGEOS-LAGEOS II combination is 3.9% (4% with EIGEN-GRACE02S, \sim 6% with EIGEN-CG01C and \sim 9% with GGM02S), while it is 1$% for the node-only LAGEOS-LAGEOS II-Ajisai-Jason-1 combination (2% with EIGEN-GRACE02S, 1.6% with EIGEN-CG01C and 2.7% with GGM02S).

In this paper we investigate the possibility of measuring the post-Newtonian general relativistic gravitomagnetic Lense-Thirring effect in the Jovian system of its Galilean satellites Io, Europa, Ganymede and Callisto in view of recent developments in processing and modelling their optical observations spanning a large time interval (125 years). The present day best observations have an accuracy between several kilometers to few tens of kilometers, which is just the order of magnitude of the Lense-Thirring shifts of the orbits of the Galilean satellites over almost a century. From a comparison between analytical development and numerical integration it turns out that, unfortunately, most of the secular component of the gravitomagnetic signature is removed in the process of fitting the initial conditions. Indeed, an estimation of the magnitude of the Lense-Thirring effect in the ephemerides residuals is given; the resulting residuals have a maximum magnitude of 20 meters only (over 125 years).

An unexpected secular increase of the Astronomical Unit, the length scale of the Solar System, has recently been reported by three different research groups (Krasinsky and Brumberg, Pitjeva, Standish). The latest JPL measurements amount to 7+-2 m cy^-1. At present, there are no explanations able to accommodate such an observed phenomenon, neither in the realm of classical physics nor in the usual four-dimensional framework of the Einsteinian General Relativity. The Dvali-Gabadadze-Porrati braneworld scenario, which is a multi-dimensional model of gravity aimed to the explanation of the observed cosmic acceleration without dark energy, predicts, among other things, a perihelion secular shift, due to Lue and Starkman, of 5 10^-4 arcsec cy^-1 for all the planets of the Solar System. It yields a variation of about 6 m cy^-1 for the Earth-Sun distance which is compatible at 1-sigma level with the observed rate of the Astronomical Unit. The recently measured corrections to the secular motions of the perihelia of the inner planets of the Solar System are in agreement, at 1-sigma level, with the predicted value of the Lue-Starkman effect for Mercury and Mars and at 2-sigma level for the Earth.

In this paper we explicitly work out the secular perturbations induced on all the Keplerian orbital elements of a test body to order O(e^2) in the eccentricity e by the weak-field long-range modifications of the usual Newton-Einstein gravity due to the Dvali-Gabadadze-Porrati (DGP) braneworld model. The Gauss perturbative scheme is used. It turns out that the argument of pericentre and the mean anomaly are affected by secular rates which are independent of the semimajor axis of the orbit of the test particle. The first nonvaishing eccentricity-dependent corrections are of order O(e^2). For circular orbits the Lue-Starkman (LS) effect on the pericentre is obtained. Some observational consequences are discussed for the Solar System planetary mean longitudes lambda which would undergo a 1.2\cdot 10^-3 arcseconds per century braneworld secular precession. According to recent data analysis over 92 years for the EPM2004 ephemerides, the 1-sigma formal accuracy in determining the Martian mean longitude amounts to 3\cdot 10^-3 milliarcseconds, while the braneworld effect over the same time span would be 1.159 milliarcseconds. The major limiting factor is the 2.6\cdot 10^-3 arcseconds per century systematic error due to the mismodelling in the Keplerian mean motion of Mars. A suitable linear combination of the mean longitudes of Mars and Venus may overcome this problem. The formal, 1-sigma obtainable observational accuracy would be \sim 7%. The systematic error due to the present-day uncertainties in the solar quadrupole mass moment, the Keplerian mean motions, the general relativistic Schwarzschild field and the asteroid ring would amount to some tens of percent.

In the context of the currently ongoing efforts to improve the accuracy and reliability of the measurement of the Lense-Thirring effect in the gravitational field of the Earth it has recently been proposed to use the data from the existing spacecraft endowed with some active mechanisms of compensation of the non-gravitational accelerations like GRACE. In this paper we critically discuss this interesting possibility. Unfortunately, it turns out to be unpracticable because of the impact of the uncancelled even zonal harmonic coefficients of the multipolar expansion of the terrestrial gravitational potential and of some time-dependent tidal perturbations which would resemble as superimposed linear trends over the necessarily limited observational time span of the analysis.

In this paper we try to suggest a possible novel method to determine some selected even zonal harmonics J_l of the Earth's geopotential. Time series many years long of suitably linearly combined residuals of some Keplerian orbital elements of certain existing geodetic SLR satellites would be examined. A CHAMP/GRACE-only background reference model should be used for the part of the geopotential which we are not interested in. The retrieved values for the even zonal harmonics of interest would be, by construction, independent of each other and of any post-Newtonian features. The so obtained mini-model could, subsequently, be used in order to enhance the accuracy and the reliability of some tests of post-Newtonian gravity, with particular emphasis to the measurement of the Lense-Thirring effect by means of LAGEOS and LAGEOS II.

With quantitative numerical simulations we show that the secular variations of the even zonal harmonics of the geopotential do affect the measurement of the Lense-Thirring effect with the nodes of LAGEOS and LAGEOS II at a \sim 14% level. Thus, the 1-sigma total error in the performed test with the 2nd generation CHAMP-only EIGEN2 Earth gravity model amounts to 51%, contrary to 18% claimed by Ciufolini.

Recently it has been suggested to use the perigee of the proposed LARES/WEBER-SAT satellite in order to measure the secular precession which would be induced on such Keplerian orbital element by a weak-field modification of gravity occurring in some Brane-World scenarios put forth by Dvali, Gabadadze and Porrati. This precession, derived for the first time by Lue and Starkman, amounts to 4\times 10-3 milliarcseconds per year for an Earth orbiting satellite. In this paper we show that, according to the recently released EIGEN-CG01C Earth gravity model the quite larger systematic errors due to the Newtonian part of the terrestrial gravitational potential would vanish any attempts to detect a so small effect. The situation is much more favorable in the Solar System scenario. Indeed, the non-Newtonian perihelion advance of Mars, which is the currently best tracked planet, has recently been measured with an accuracy of 1\times 10^-4 arcseconds per century; the Dvali precession is 4\times 10^-4 arcseconds per century. A suitable combination of the perihelia and the nodes of Mars and Mercury, which disentangles the Dvali effect from the competing larger Newtonian and general relativistic precessions, might allow to reach a 1-sigma 47% level of accuracy.

In this paper we will show in detail that the performed attempts aimed at the detection of the general relativistic Lense-Thirring effect in the gravitational field of the Earth with the existing LAGEOS satellites are often presented in an optimistic and misleading way which is inadequate for such an important test of fundamental physics. E.g., in the latest reported measurement of the gravitomagnetic shift with the nodes of the LAGEOS satellites and the 2nd generation GRACE-only EIGEN-GRACE02S Earth gravity model over an observational time span of 11 years a 5-10% total accuracy is claimed at 1-3sigma, respectively. We will show that, instead, it might be 15-45% (1-3sigma) if the impact of the secular variations of the even zonal harmonics is considered as well.

In this paper we investigate the opportunities offered by the new Earth gravity models from the dedicated CHAMP and, especially, GRACE missions to the project of measuring the general relativistic Lense-Thirring effect with a new Earth's artificial satellite. It turns out that it would be possible to abandon the stringent, and expensive, requirements on the orbital geometry of the originally prosed LARES mission (same semimajor axis a=12270 km of the existing LAGEOS and inclination i=70 deg) by inserting the new spacecraft in a relatively low, and cheaper, orbit (a=7500-8000 km, i\sim 70 deg) and suitably combining its node Omega with those of LAGEOS and LAGEOS II in order to cancel out the first even zonal harmonic coefficients of the multipolar expansion of the terrestrial gravitational potential J_2, J_4 along with their temporal variations. The total systematic error due to the mismodelling in the remaining even zonal harmonics would amount to \sim 1% and would be insensitive to departures of the inclination from the originally proposed value of many degrees. No semisecular long-period perturbations would be introduced because the period of the node, which is also the period of the solar K_1 tidal perturbation, would amount to \sim 10^2 days. Since the coefficient of the node of the new satellite would be smaller than 0.1 for such low altitudes, the impact of the non-gravitational perturbations of it on the proposed combination would be negligible. Then, a particular financial and technological effort for suitably building the satellite in order to minimize the non-conservative accelerations would be unnecessary.

The possibility of measuring the post-Newtonian gravitoelectric correction to the orbital period of a test particle freely orbiting a spherically symmetric mass in the Solar System is analyzed. It should be possible, in principle, to detect it for Mercury at a precision level of 10^-4. This level is mainly set by the unavoidable systematic errors due to the mismodelling in the Keplerian period which could not be reduced by accumulating a large number of orbital revolutions. Future missions like Messenger and BepiColombo should allow to improve it by increasing our knowledge of the Mercury's orbital parameters. The observational accuracy is estimated to be 10^-4 from the knowledge of the International Celestial Reference Frame (ICRF) axes. It could be improved by observing as many planetary transits as possible. It is not possible to measure such an effect in the gravitational field of the Earth by analyzing the motion of artificial satellites or the Moon because of the unavoidable systematic errors related to the uncertainties in the Keplerian periods. In the case of some recently discovered exoplanets the problems come from the observational errors which are larger than the relativistic effect.

In this paper we discuss and compare a node-only LAGEOS-LAGEOS II combination and a node-only LAGEOS-LAGEOS II-Ajisai-Jason1 combination for the determination of the Lense-Thirring effect. The new combined EIGEN-CG01C Earth gravity model has been adopted. The second combination cancels the first three even zonal harmonics along with their secular variations but introduces the non-gravitational perturbations of Jason1. The first combination is less sensitive to the non-conservative forces but is sensitive to the secular variations of the uncancelled even zonal harmonics of low degree J4 and J6 whose impact grows linearly in time.

In this paper the effect of the post-Newtonian gravitomagnetic force on the mean longitudes $l$ of a pair of counter-rotating Earth artificial satellites following almost identical circular equatorial orbits is investigated. The possibility of measuring it is examined. The observable is the difference of the times required to $l$ in passing from 0 to 2$\pi$ for both senses of motion. Such gravitomagnetic time shift, which is independent of the orbital parameters of the satellites, amounts to 5$\times 10^{-7}$ s for Earth; it is cumulative and should be measured after a sufficiently high number of revolutions. The major limiting factors are the unavoidable imperfect cancellation of the Keplerian periods, which yields a constraint of 10$^{-2}$ cm in knowing the difference between the semimajor axes $a$ of the satellites, and the difference $I$ of the inclinations $i$ of the orbital planes which, for $i\sim 0.01^\circ$, should be less than $0.006^\circ$. A pair of spacecrafts endowed with a sophisticated intersatellite tracking apparatus and drag-free control down to 10$^{-9}$ cm s$^{-2}$ Hz$^{-{1/2}}$ level might allow to meet the stringent requirements posed by such a mission.

In this brief note some comments about the observable used in a recently published paper on the measurement of the general relativistic Lense-Thirring in the gravitational field of the Earth are presented. It turns out that, among other things, the authors might have yielded an optimistic evaluation of the error budget because of an underestimation of the impact of the secular variations of the even zonal harmonics of the geopotential. More realistic evaluations point towards a 15-45% error at 1-3sigma level, respectively.

In this paper we first calculate the post-Newtonian gravitoelectric secular rate of the mean anomaly of a test particle freely orbiting a spherically symmetric central mass. Then, we propose a novel approach to suitably combine the presently available planetary ranging data to Mercury, Venus and Mars in order to determine, simultaneously and independently of each other, the Sun's quadrupole mass moment J_2 and the secular advances of the perihelion and the mean anomaly. This would also allow to obtain the PPN parameters gamma and beta independently. We propose to analyze the time series of three linear combinations of the experimental residuals of the rates of the nodes, the longitudes of perihelia and mean anomalies of Mercury, Venus and Mars built up in order to absorb the secular precessions induced by the solar oblateness and the post-Newtonian gravitoelectric forces. The values of the three investigated parameters can be obtained by fitting the expected linear trends with straight lines, determining their slopes in arcseconds per century and suitably normalizing them. According to the present-day EPM2000 ephemerides accuracy, the obtainable precision would be of the order of 10^-4-10^-5 for the PPN parameters and, more interestingly, of 10^-9 for J_2. The future BepiColombo mission should improve the Mercury's orbit by one order of magnitude.

Here we explore a novel approach in order to try to measure the post-Newtonian 1/c^2 Lense-Thirring secular effect induced by the gravitomagnetic field of the Sun on the planetary orbital motion. Due to the relative smallness of the solar angular momentum J and the large values of the planetary semimajor axes a, the gravitomagnetic precessions, which affect the nodes Omega and the perihelia omega and are proportional to J/a^3, are of the order of 10^-3 arcseconds per century only for, e.g., Mercury. This value lies just at the edge of the present-day observational sensitivity in reconstructing the planetary orbits, although future missions to Mercury like Messenger and BepiColombo could allow to increase it. The major problems come from the main sources of systematic errors. They are the aliasing classical precessions induced by the multipolar expansion of the Sun's gravitational potential and the classical secular N-body precessions which are of the same order of magnitude or much larger than the Lense-Thirring precessions of interest. This definitely rules out the possibility of analyzing only one orbital element of, e.g., Mercury. In order to circumvent these problems, we propose a suitable linear combination of the orbital residuals of the nodes of Mercury, Venus and Mars which is, by construction, independent of such classical secular precessions. A 1-sigma reasonable estimate of the obtainable accuracy yields a 36% error. Since the major role in the proposed combination is played by the Mercury's node, it could happen that the new, more accurate ephemerides available in future thanks to the Messenger and BepiColombo missions will offer an opportunity to improve the present unfavorable situation.

The possibility of analyzing the node Omega of the GP-B satellite in order to measure also the Lense-Thirring effect on its orbit is examined. This feature is induced by the general relativistic gravitomagnetic component of the Earth gravitational field. The GP-B mission has been launched in April 2004 and is aimed mainly to the measurement of the gravitomagnetic precession of four gyroscopes carried onboard at a claimed accuracy of 1% or better. The aliasing effect of the solid Earth and ocean components of the solar K_1 tidal perturbations would make the measurement of the Lense-Thirring effect on the orbit unfeasible. Indeed, the science period of the GP-B mission amounts to almost one year. During this time span the Lense-Thirring shift on the GP-B node would be 164 milliarcseconds (mas), while the tidal perturbations on its node would have a period of the order of 10^3 years and amplitudes of the order of 10^5 mas.

Here we explore the possibility of extending to the Solar System scenario the combined residual approach employed for the deterination of the Lense-Thirring effect in the Earth-LAGEOS satellites system. After calculating the secular gravitoelectric advance of the mean anomaly, we propose several linear combinations of the residuals of the nodes, the perihelia and the mean anomalies of some inner planets in order to measure the gravitoelectric advance, the Lense-Thirring effect and the solar quadrupole mass moment. The obtainable accuracies are comparable to those of new complex and expensive missions involving the launch of one or more spacecrafts. They are 10^-4 for the PPN parameter beta, 10^-5 for gamma, 10^-9 for the solar J2, 20-30% for the Lense-Thirring effect. In the case of the solar gravitomagnetic field, it would be the first attempt to measure its effect on the planetary motions. However, the improvements in the ephemerides of Mercury thanks to the Messenger and BepiColombo missions will map into an increasing accuracy of the proposed tests.

In this paper we treat the so called clock paradox in an analytical way by assuming that a constant and uniform force F of finite magnitude acts continuously on the moving clock along the direction of its motion assumed to be rectilinear. No inertial motion steps are considered. The rest clock is denoted as (1), the to-and-fro moving clock is (2), the inertial frame in which (1) is at rest in its origin and (2) is seen moving is I and, finally, the accelerated frame in which (2) is at rest in its origin and (1) moves forward and backward is A. We deal with the following questions: I) What is the effect of the finite force acting on (2) on the proper time intervals measured by the two clocks when they reunite? Does a differential aging between the two clocks occur, as it happens when inertial motion and infinite values of the accelerating force is considered? The Special Theory of Relativity is used in order to describe the hyperbolic motion of (2) in the frame I II) Is this effect an absolute one, i.e. does the accelerated observer A comoving with (2) obtain the same results as that in I, both qualitatively and quantitatively, as it is expected? We use the General Theory of Relativity in order to answer this question.

In this paper we deal analytically with a version of the so called clock paradox in which the moving clock performs a circular motion of constant radius. The rest clock is denoted as (1), the rotating clock is (2), the inertial frame in which (1) is at rest and (2) moves is I and, finally, the accelerated frame in which (2) is at rest and (1) rotates is A. By using the General Theory of Relativity in order to describe the motion of (1) as seen in A we will show the following features. I) A differential aging between (1) and (2) occurs at their reunion and it has an absolute character, i.e. the proper time interval measured by a given clock is the same both in I and in A. II) From a quantitative point of view, the magnitude of the differential aging between (1) and (2) does depend on the kind of rotational motion performed by A. Indeed, if it is uniform there is no any tangential force in the direction of motion of (2) but only normal to it. In this case, the proper time interval reckoned by (2) does depend only on its constant velocity v=romega. On the contrary, if the rotational motion is uniformly accelerated, i.e. a constant force acts tangentially along the direction of motion, the proper time intervals $do depend$ on the angular acceleration alpha. III) Finally, in regard to the sign of the aging, the moving clock (2) measures always a $shorter$ interval of proper time with respect to (1).

We examine how the new forthcoming Earth gravity models from the CHAMP and, especially, GRACE missions could improve the measurement of the general relativistic Lense-Thirring effect according to the various kinds of observables which could be adopted. In a very preliminary way, we use the recently released EIGEN2 CHAMP-only and GRACE01S GRACE-only Earth gravity models in order to assess the impact of the mismodelling in the even zonal harmonic coefficients of geopotential which represents one of the major sources of systematic errors in this kind of measurement.

In this paper we wish to preliminary investigate if it would be possible to use the orbital data from the proposed OPTIS mission together with those from the existing geodetic passive SLR LAGEOS and LAGEOS II satellites in order to perform precise measurements of some general relativistic gravitoelectromagnetic effects, with particular emphasis on the Lense-Thirring effect.

In this paper we use, in a preliminary way, the recently released EIGEN2 Earth gravity model, which is based on six months of data of CHAMP only, in order to reassess the systematic error due to the mismodelling in the even zonal harmonics of geopotential in the LAGEOS-LAGEOS II Lense-Thirring experiment involving the nodes of both the LAGEOS satellites and the perigee of LAGEOS II. The first results from the GGM01C Earth gravity model including the first GRACE data are very promising.

The possibility of performing post-Newtonian gravitoelectromagnetic measurements with a joint LAGEOS-LAGEOS II-OPTIS space-based mission is investigated

In this contribution we present two new proposals for measuring the general relativistic gravitomagnetic component of the gravitational field of the Earth. One proposal consists of the measurement of the difference of the rates of the perigee $\psi$ from the analysis of the laser--ranged data of two identical Earth'artificial satellites placed in equal orbits with supplementary inclinations. In this way the impact of the aliasing classical secular precessions due to the even zonal harmonics of the geopotential would be canceled out, although the non--gravitational perturbations, to which the perigees of LAGEOS--type satellites are particularly sensitive, should be a limiting factor in the obtainable accuracy. With a suitable choice of the inclinations of the orbital planes it would be possible to reduce the periods of such insidious perturbations so to use not too long observational time spans. However, the use of a pair of drag--free satellites would greatly reduce this problem, provided that the time span of the data analysis does not excess the lifetime of the drag--free apparatus. In the other proposal the difference of the rotational periods of two counter-revolving particles placed on a friction-free plane in a vacuum chamber at the South Pole should be measured in order to extract the relativistic gravitomagnetic signal. Among other very challenging practical implications, the Earth's angular velocity $\omega_{\oplus}$ should be known at a $10^{-15}$ rad s$^{-1}$ level from VLBI and the friction force of the plane should be less than $2\times 10^{-9}$ dyne.

Among the effects predicted by the General Theory of Relativity for the orbital motion of a test particle, the post-Newtonian gravitomagnetic Lense-Thirring effect is very interesting and, up to now, there is not yet an undisputable experimental direct test of it. To date, the data analysis of the orbits of the existing geodetic LAGEOS and LAGEOS II satellites has yielded a test of the Lense-Thirring effect with a claimed accuracy of 20%-30%. According to some scientists such estimates could be optimistic. Here we wish to discuss the improvements obtainable in this measurement, in terms of reliability of the evaluation of the systematic error and reduction of its magnitude, due to the new CHAMP and GRACE Earth gravity models.

In this paper we examine the possibility of testing the equivalence principle, in its weak form, by analyzing the orbital motion of a pair of artificial satellites of different composition moving along orbits of identical shape and size in the gravitational field of Earth. It turns out that the obtainable level of accuracy is, realistically, of the order of 10^-10 or slightly better. It is limited mainly by the fact that, due to the unavoidable orbital injection errors, it would not be possible to insert the satellites in orbits with exactly the same radius and that such difference could be known only with a finite precision. The present-day level of accuracy, obtained with torsion balance Earth-based measurements and the analysis of Earth-Moon motion in the gravitational field of Sun with the Lunar Laser Ranging technique, is of the order of 10^-13. The proposed space-based missions STEP, \muSCOPE, GG and SEE aim to reach a 10^-15-10^-18 precision level.

In this paper the most recent developments in testing General Relativity in the gravitational field of the Earth with the technique of Satellite Laser Ranging are presented. In particular, we concentrate our attention on some gravitoelectric and gravitomagnetic post--Newtonian orbital effects on the motion of a test body in the external field of a central mass.

Testing the effects predicted by the General Theory of Relativity, in its linearized weak field and slow motion approximation, in the Solar System is difficult because they are very small. Among them the post-Newtonian gravitomagnetic Lense-Thirring effect, or dragging of the inertial frames, on the orbital motion of a test particle is very interesting and, up to now, there is not yet an undisputable experimental direct test of it. Here we illustrate how it could be possible to measure it with an accuracy of the order of 1%, together with other tests of Special Relativity and post-Newtonian gravity, with a joint space based OPTIS/LARES mission in the gravitational field of Earth. Up to now, the data analysis of the orbits of the existing geodetic LAGEOS and LAGEOS II satellites has yielded a test of the Lense-Thirring effect with a claimed accuracy of 20%-30%.

In this paper we propose, in a preliminary way, a new Earth-based laboratory experiment aimed to the detection of the gravitomagnetic field of the Earth. It consists of the measurement of the difference of the circular frequencies of two rotators moving along identical circular paths, but in opposite directions, on a horizontal friction-free plane in a vacuum chamber placed at South Pole. The accuracy of our knowledge of the Earth's rotation from VLBI and the possibility of measuring the rotators'periods over many revolutions should allow for the feasibility of the proposed experiment.

In this paper we consider the possibility of measuring the corrections induced by the square of the parameter a_g of the Kerr metric to the general relativistic deflection of electromagnetic waves and time delay in an Earth based experiment. It turns out that, while at astronomical scale the well known gravitoelectric effects are far larger than the gravitomagnetic ones, at laboratory scale the situation is reversed: the gravitomagnetic effects exceed definitely the gravitoelectric ones which are totally negligible. By using a small rapidly rotating sphere as gravitating source on the Earth the deflection of a grazing light ray amounts to 10^{-13} rad and the time delay is proportional to 10^{-23} s. These figures are determined by the upper limit in the attainable values of a_g due to the need of preventing the body from exploding under the action of the centrifugal forces. Possible criticisms to the use of the Kerr metric at a_g^2 level are discussed.

The gravitomagnetic corrections to the Keplerian period for a circular, geodesic orbit of a test particle in a polar plane containing the proper angular momentum J of a central rotating body are considered.

In this paper we analyze quantitatively the concept of LAGEOS--type satellites in critical supplementary orbit configuration (CSOC) which has proven capable of yielding various observables for many tests of General Relativity in the terrestrial gravitational field, with particular emphasis on the measurement of the Lense--Thirring effect.

In this paper we present a rather extensive error budget for the difference of the perigees of a pair of supplementary SLR satellites aimed to the detection of the Lense-Thirring effect.

In this letter we reexamine the evaluation of \zone in some proposed tests of relativistic gravitomagnetism with existing and proposed laser--ranged LAGEOS--like satellites in the gravitational field of the Earth. A more conservative and realistic approach is followed by using the diagonal part only of the covariance matrix of the EGM96 Earth's gravity model up to degree l=20. It turns out that, within this context and according to the present level of knowledge of the terrestrial gravitational field, the best choice would be the use of a recently proposed combination using the nodes \Omega of LAGEOS, LAGEOS II and LARES and the perigees \omega of LAGEOS II and LARES. Indeed, it turns out to be insensitive both to the even zonal harmonics of degree higher than l=20 and to the correlation among them

In this letter we further analyze the concept of the use of a pair of identical geodetic satellites placed in identical orbits except for the inclinations which should be supplementary for measuring the Lense-Thirring effect in the terrestrial gravitational field. It turns out that not only the sum of the nodes, as already proposed for the LAGEOS-LARES mission, but also the difference of the perigees could be fruitfully adopted.

In this paper we calculate explicitly the secular classical precessions of the node \Omega and the perigee \omega of an Earth artificial satellite induced by the static, even zonal harmonics of the geopotential up to degree l=20. Subsequently, their systematic errors induced by the mismodelling in the even zonal geopotential coefficients J_l are compared to the general relativistic secular gravitomagnetic and gravitoelectric precessions of the node and the perigee of the existing laser-ranged geodetic satellites and of the proposed LARES.

In this paper we look at the gravitational spin--spin interaction between macroscopic astronomical bodies. In particular, we calculate their post--Newtonian orbital effects of order $\mathcal{O}(c^{-2})$ on the trajectory of a spinning particle with proper angular momentum ${\bf s}$ moving in the external gravitomagnetic field generated by a central spinning mass with proper angular momentum ${\bf J}$. It turns out that, at order $\mathcal{O}(e)$ in the orbiter's eccentricity, the eccentricity the pericenter and the mean anomaly rates of the moving particle are affected by long--term harmonic effects. If, on one hand, they are undetectable in the Solar System, on the other, maybe that in an astrophysical context like that of the binary millisecond pulsars there will be some hopes of measuring them in the future.

According to very recent developments of the LARES mission, which would be devoted to the measurement of the general relativistic Lense--Thirring effect in the gravitational field of the Earth with Satellite Laser Ranging, it seems that the LARES satellite might be finally launched in a polar, low--altitude orbit by means of a relatively low--cost rocket. The observable would be the node only. In this letter we critically analyze this scenario.

A new proposal for the implementation of LARES mission is presented. In particular, a new observable is proposed and alternative scenarios are discussed.

One of the most famous classical tests of General Relativity is the gravitoelectric secular advance of the pericenter of a test body in the gravitational field of a central mass. In this paper we explore the possibility of performing a measurement of the gravitoelectric pericenter advance in the gravitational field of the Earth by analyzing the laser-ranged data to some existing, or proposed, laser-ranged geodetic satellites. At the present level of knowledge of various error sources, the relative precision obtainable with the data from LAGEOS and LAGEOS II, suitably combined, is of the order of $10^{\rm -3}$. Nevertheless, these accuracies could sensibly be improved in the near future when the new data on the terrestrial gravitational field from the CHAMP and GRACE missions will be available. The use of the perigee of LARES (LAser RElativity Satellite), in the context of a suitable combination of orbital residuals including also LAGEOS II, should further raise the precision of the measurement. As a secondary outcome of the proposed experiment, with the so obtained value of $\ppn$ and with $\et=4\beta-\gamma-3$ from Lunar Laser Ranging it could be possible to obtain an estimate of the PPN parameters $\gamma$ and $\beta$ at the $10^{-2}-10^{-3}$ level.

In this work some aspects of the detection of certain general relativistic effects in the weak gravitational field of the Earth via laser-ranged data to some existing or proposed geodetic satellites are examined. The focus is on the Lense-Thirring drag of the orbit of a test body, the gravitomagnetic clock effect and the gravitoelectric perigee shift. The impact of some sources of systematic errors is investigated. An experiment whose goal is the measurement of the PPN parameters beta and gamma in the terrestrial field with LAGEOS satellites at a level of 10^(-3)is presented. A modified version of the proposed LARES mission is examined.

In this paper we investigate the possibility of constraining the hypothesis of a fifth force at the length scale of two Earth's radii by investigating the effects of a Yukawa gravitational potential on the orbits of the laser--ranged LAGEOS satellites. The existing constraints on the Yukawa coupling $\alpha$, obtained by fitting the LAGEOS orbit, are of the order of | \alpha | < 10^{-5}-10^{-8} for distances of the order of 10^9 cm. Here we show that with a suitable combination of the orbital residuals of the perigee \omega of LAGEOS II and the nodes \Omega of LAGEOS II and LAGEOS it should be possible to constrain \alpha at a level of 4 X 10^{-12} or less. Various sources of systematic errors are accounted for, as well. Their total impact amounts to 1 X 10^{-11} during an observational time span of 5 years. In the near future, when the new data on the terrestrial gravitational field will be available from the CHAMP and GRACE missions, these limits will be further improved. The use of the proposed LARES laser--ranged satellite would yield an experimental accuracy in constraining \alpha of the order of 1 X 10^{-12}.

In this paper the possibility of measuring some general relativistic effects in the gravitational field of the Moon via selenodetic missions, with particular emphasis to the future Japanese SELENE mission, is investigated. For a typical selenodetic orbital configuration the post-Newtonian Lense-Thirring gravitomagnetic and the Einstein's gravitoelectric effects on the satellites orbits are calculated and compared to the present-day orbit accuracy of lunar missions. It turns out that for SELENE's Main Orbiter, at present, the gravitoelectric periselenium shift, which is the largest general relativistic effect, is 1 or 2 orders of magnitude smaller than the experimental sensitivity. The systematic error induced by the mismodelled classical periselenium precession due to the first even zonal harmonic J2 of the Moon's non-spherical gravitational potential is 3 orders of magnitude larger than the general relativistic gravitoelectric precession. The estimates of this work could be used for future lunar missions having as their goals relativistic measurements as well.

The well known general relativistic Lense-Thirring drag of the orbit of a test particle in the stationary field of a central slowly rotating body is generated, in the weak-field and slow-motion approximation of General Relativity, by a gravitomagnetic Lorentz-like acceleration in the equations of motion of the test particle. In it the gravitomagnetic field is due to the central body's angular momentum supposed to be constant. In the context of the gravitational analogue of the Larmor theorem, such acceleration looks like a Coriolis inertial term in an accelerated frame. In this paper the effect of the variation in time of the central body's angular momentum on the orbit of a test mass is considered. It can be shown that it is analogue to the inertial acceleration due to the time derivative of the angular velocity vector of an accelerated frame. The possibility of detecting such effect in the gravitational field of the Earth with LAGEOS-like satellites is investigated. It turns out that the orbital effects are far too small to be measured.

The accuracy reached in the past few years by Satellite Laser Ranging (SLR) allows for measuring even tiny features of the Earth's gravitational field predicted by Einstein's General Relativity by means of artificial satellites. The gravitomagnetic dragging of the orbit of a test body is currently under measurement by analyzing a suitable combination of the orbital residuals of LAGEOS and LAGEOS II. The lower bound of the error in this experiment amount to 12.92%. It is due to the mismodeling in the even zonal harmonics of the geopotential which are the most important sources of systematic error. A similar approach could be used in order to measure the relativistic gravitoelectric pericenter shift in the field of the Earth with a lower bound of the systematic relative error of 0.6% due to the even zonal harmonics as well. The inclusion of the ranging data to the Japanese passive geodetic satellite Ajisai would improve such limits to 10.78% and 0.08% respectively and would allow to improve the accuracy in the determination of the PPN parameters beta and gamma.

The possibility of detecting the gravitomagnetic clock effect using artificial Earth satellites provides the incentive to develop a more intuitive approach to its derivation. We first consider two test electric charges moving on the same circular orbit but in opposite directions in orthogonal electric and magnetic fields and show that the particles take different times in describing a full orbit. The expression for the time difference is completely analogous to that of the general relativistic gravitomagnetic clock effect in the weak-field and slow-motion approximation. The latter is obtained by considering the gravitomagnetic force as a small classical non-central perturbation of the main central Newtonian monopole force. A general expression for the clock effect is given for a spherical orbit with an arbitrary inclination angle. This formula differs from the result of the general relativistic calculations by terms of order c^{-4}.

The general relativistic gravitomagnetic clock effect consists in the fact that two point particles orbiting a central spinning object along identical, circular equatorial geodesic paths, but in opposite directions, exhibit a time difference in describing a full revolution. It turns out that the particle rotating in the same sense of the central body is slower than the particle rotating in the opposite sense. In this paper it is proposed to measure such effect in an Earth laboratory experiment involving interferometry of slow neutrons. With a sphere of 2.5 cm radius and spinning at 4.3 x 10^4 rad/s as central source, and using neutrons with wavelength of 1 Angstrom it should be possible to obtain, for a given sense of rotation of the central source, a phase shift of 0.18 rad, well within the experimental sensitivity. By reversing the sense of rotation of the central body it should be possible to obtain a 0.06 fringe shift.

The general relativistic gravitomagnetic clock effect consists in the fact that two massive test bodies orbiting a central spinning mass in its equatorial plane along two identical circular trajectories, but in opposite directions, take different times in describing a full revolution with respect to an asymptotically inertial observer. In the field of the Earth such time shift amounts to 10^{-7} s. Detecting it by means of a space based mission with artificial satellites is a very demanding task because there are severe constraints on the precision with which the radial and azimuthal positions of a satellite must be known: delta r= 10^{-2} cm and delta phi= 10^{-2} milliarcseconds per revolution. In this paper we assess if the systematic errors induced by various non-gravitational perturbations allow to meet such stringent requirements. A couple of identical, passive laser-ranged satellites of LAGEOS type with their spins aligned with the Earth's one is considered. It turns out that all the non vanishing non-gravitational perturbations induce systematic errors in r and phi within the required constraints for a reasonable assumption of the mismodeling in some satellite's and Earth's parameters and/or by using dense satellites with small area-to-mass ratio. However, the error in the Earth's GM is by far the largest source of uncertainty in the azimuthal location which is affected at a level of 1.2 milliarcseconds per revolution.

The general relativistic Lense-Thirring effect can be measured by inspecting a suitable combination of the orbital residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. The solid and ocean Earth tides affect the recovery of the parameter by means of which the gravitomagnetic force is accounted for in the combined residuals. Thus an extensive analysis of the perturbations induced on these orbital elements by the solid and ocean Earth tides is carried out. It involves the l=2 terms for the solid tides and the l=2,3,4 terms for the ocean tides. The perigee of LAGEOS II turns out to be very sensitive to the l=3 part of the ocean tidal spectrum, contrary to the nodes of LAGEOS and LAGEOS II. The uncertainty in the solid tidal perturbations, mainly due to the Love number k2, ranges from 0.4% to 1.5%, while the ocean tides are uncertain at 5%-15% level. The obtained results are used in order to check in a preliminary way which tidal constituents the Lense-Thirring shift is sensitive to. In particular it is tested if the semisecular 18.6-year zonal tide really does not affect the combined residuals. It turns out that, if modeled at the level of accuracy worked out in the paper, the l=2,4 m=0 and also, to a lesser extent, the l=3 m=0 tidal perturbations cancel out.

The general relativistic Lense-Thirring effect can be detected by means of a suitable combination of orbital residuals of the laser-ranged LAGEOS and LAGEOS II satellites. While this observable is not affected by the orbital perturbation induced by the zonal Earth solid and ocean tides, it is sensitive to those generated by the tesseral and sectorial tides. The assessment of their influence on the measurement of the parameter mu, with which the gravitomagnetic effect is accounted for, is the goal of this paper. After simulating the combined residual curve by calculating accurately the mismodeling of the more effective tidal perturbations, it has been found that, while the solid tides affect the recovery of mu at a level always well below 1%, for the ocean tides and the other long-period signals Delta mu depends strongly on the observational period and the noise level: Delta mu(tides) amounts to almost 2% after 7 years. The aliasing effect of K1 l=3 p=1 tide and SRP(4241) solar radiation pressure harmonic, with periods longer than 4 years, on the perigee of LAGEOS II yield to a maximum systematic uncertainty on $\m_{LT}$ of less than 4% over different observational periods. The zonal 18.6-year tide does not affect the combined residuals.

The tiny general relativistic Lense-Thirring effect can be measured by means of a suitable combination of the orbital residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. This observable is affected, among other factors, by the Earth' s solid and ocean tides. They induce long-period orbital perturbations that, over observational periods of few years, may alias the detection of the gravitomagnetic secular trend of interest. In this paper we calculate explicitly the most relevant tidal perturbations acting upon LAGEOSs and assess their influence on the detection of the Lense-Thirring effect. The present day level of knowledge of the solid and ocean tides allow us to conclude that their influence on it ranges from almost 4% over 4 years to less than 2% over 7 years.

The detection of some tiny gravitomagnetic effects in the field of the Earth by means of artificial satellites is a very demanding task because of the various other perturbing forces of gravitational and non-gravitational origin acting upon them. Among the gravitational perturbations a relevant role is played by the Earth solid and ocean tides. In this communication I outline their effects on the detection of the Lense-Thirring drag of the orbits of LAGEOS and LAGEOS II, currently analyzed, and the proposed GP-C experiment devoted to the measurement of the clock effect.