- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Physics Research International
Volume 2012 (2012), Article ID 469095, 12 pages
Conformal Cosmology and the Pioneer Anomaly
Department of Physics, Loyola Marymount University, Los Angeles, CA 90045, USA
Received 9 August 2011; Revised 9 October 2011; Accepted 10 October 2011
Academic Editor: A. Beesham
Copyright © 2012 Gabriele U. Varieschi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
We review the fundamental results of a new cosmological model, based on conformal gravity, and apply them to the analysis of the early data of the Pioneer anomaly. We show that our conformal cosmology can naturally explain the anomalous acceleration of the Pioneer 10 and 11 spacecrafts, in terms of a local blueshift region extending around the solar system and therefore affecting the frequencies of the navigational radio signals exchanged between Earth and the spacecraft. By using our model, we explain the numerical coincidence between the value of the anomalous acceleration and the Hubble constant at the present epoch and also confirm our previous determination of the cosmological parameters and . New Pioneer data are expected to be publicly available in the near future, which might enable more precise evaluations of these parameters.
The Pioneer 10 and 11 spacecrafts were launched in the early 1970s, to conduct explorations in the region of the solar system beyond the orbit of Mars and to perform close observations of Jupiter. They were also the first spacecraft to explore the outer solar system and to send back to Earth their navigational signals for almost thirty years (for a review see  and references therein).
In recent years, the orbits of Pioneer 10 and 11 were reconstructed very accurately, by using the original radiometric Doppler tracking data, based on the signals exchanged between the spacecraft and NASA’s terrestrial tracking stations. This reconstruction yielded a persistent discrepancy between the observed and predicted data, equivalent to an unexplained small acceleration of the spacecraft in the direction of the Sun. This effect is evidenced by measuring a small frequency shift (toward higher frequencies, i.e., a “blueshift”) of the signal reaching us from the spacecraft. The nature of this anomalous acceleration or of the related blueshift remains unexplained; this effect has become known as the “Pioneer anomaly” ([2–4]).
This is not the only known gravitational anomaly in the solar system, since several others are currently under investigation (for general reviews see [5, 6]), such as the secular increase of the astronomical unit , the anomalies in planetary flybys ([8–10]), the anomalous perihelion precession of Saturn ([11, 12]), the increase in the eccentricity of the Moon’s orbit ([13, 14]), and other effects related in general with ephemerides of planets and the Moon ([15, 16]).
The importance of all these effects is not related to how they affect the spacecraft navigation, since they all produce very small corrections to the orbits, but to the possibility that these anomalies might be an indication of new gravitational physics. In particular, several nonconventional explanations of these effects have been proposed (see general discussion in [1, 2, 4]) such as modifications of the law of gravity, or a modified inertia, as proposed by the Modified Newtonian Dynamics (MOND) theory, the existence of a dark matter halo around the Earth, or in the solar system, which might slightly alter the gravitational force acting on the spacecraft, and several others ([17–21]).
In this line of reasoning, alternative gravitational theories such as conformal gravity (CG), originally proposed by Weyl in 1918 ([22–24]) and revisited by Mannheim and Kazanas ([25–27]), have provided a new framework for cosmological models, with the advantage of avoiding some of the most controversial elements of current standard cosmology, such as dark matter, dark energy, inflation, and others.
Following the original CG, we have recently studied an alternative approach to these models which was named “kinematical conformal cosmology” , but that for brevity will be called conformal cosmology (CC) in the rest of this paper. This approach was based on the direct application of the conformal symmetry to the Universe, that is, considering the possibility that a “stretching” of the spacetime fabric might be acting over cosmological scales. In a second part of this work , it was shown that this model can successfully fit type-Ia Supernovae data, without assuming the existence of dark matter or dark energy.
A preliminary analysis also performed in our second paper  indicated that CC might be able to explain the existence of the Pioneer anomaly, since the observed blueshift of the spacecraft signal could be due to a region of cosmological blueshift surrounding our solar system, which is naturally predicted by our model. A new comprehensive review of the Pioneer anomaly has been published , together with more details of the Pioneer early data , thus prompting us to reconsider and improve our previous analysis , based on the conformal cosmology approach. In addition, a revised analysis of the Pioneer anomaly, based on extended data sets, has recently appeared in the literature , confirming the existence of the anomaly and adding new insights into its temporally varying behavior.
In the next section, we will briefly review our CC solutions, showing how a local blueshift region can naturally emerge, while, in Section 3, we will fit all current Pioneer data  with our cosmological solutions and determine the values of the parameters in our model. Finally, in Section 4, we will discuss our results and compare them to the existing physical limits of standard gravity in the solar system.
2. Conformal Cosmology
In our first CC paper , we used as a starting point the line element originally derived by Mannheim and Kazanas  as an exterior solution for a static, spherically symmetric source in conformal gravity theory, that is, the analogue of the Schwarzschild exterior solution in general relativity: where in spherical coordinates and with the parameters , , , where is the mass of the (spherically symmetric) source and is the gravitational constant. Conformal gravity introduces two new parameters and which are not present in standard general relativity, while the familiar Schwarzschild solution is recovered in the limit for , in the equations above.
We then considered regions far away from matter distributions, thus ignoring the matter dependent terms, and rewrote the last equation in a simplified form: where the parameter is linked to and , through , and it is ultimately connected to the so-called trichotomy constant (in bold) of a Robertson-Walker (RW) metric, defined as . This is also related to another fundamental aspect of CG: the existence of coordinate and conformal transformations connecting the static, spherically symmetric solution represented by (1) and (3), with the classical Robertson-Walker metric (see details in ).
It was precisely this connection between the two solutions which prompted us to consider the CG static, spherically symmetric solution as an alternative description of the standard cosmological evolution, based on the RW metric. In other words, the CG static solution might also contain information about the cosmological redshift, the expansion of the Universe, and so forth, and constitute an alternative approach to cosmology. In particular, the CG expressions in (2) or (3) contain a linear and a quadratic term, in the radial coordinate , which might yield considerable effects at large distances, (by using the contribution of the linear term , the flat galactic rotation curves were in fact explained by Mannheim ([32, 33]) without the need of dark matter.) including a strong gravitational redshift which could be, at least in part, responsible for the observed cosmological redshift.
Therefore, we postulated in  that the observed redshift is due to this gravitational effect, which influences the wavelength or frequency of a light signal emitted at time and position , and observed at the origin () at the current time , in the following way:
In the previous equation, the redshift factor is related to the ratio of cosmic scale factors , which simply depend on the radial distance , in view of (3). Alternatively, to obtain a time-dependent form of the cosmic scale factor, we considered that the radial distance is associated with a look-back time , related to the time of travel of a light signal. Integrating the CG metric in (1) and (3) along the null geodesic, we obtained  for the three possible values of the parameter . In the previous equation, we preferred to use dimensionless quantities and parameters, defined as follows: so that the fundamental parameters of our conformal cosmology are now expressed by and the dimensionless ( is the speed of light in vacuum, assumed constant).
In Figure 1, we plot the results of (5) in terms of the inverse ratio which describes better the cosmic evolution. The dimensionless quantity , on the horizontal axis, represents a look-back time, so that the universal evolution of the cosmic scale factor, from the past to the future, can be seen by following our curves from right to left. The circular dot on the vertical axis represents our “current time” (). We can clearly see that the only solution which shows a redshift in the past (values below the horizontal black dashed line, representing ) is the red-solid curve, corresponding to . Therefore, the other two solutions, for , 0, are ruled out; only the solution will be considered in the following.
Our preferred solution in Figure 1 (red-solid) also shows a blueshift region in the immediate past of our current time, which in Section 3 will be related directly to the Pioneer anomaly. This blueshift region is greatly exaggerated in the figure, since the different curves were plotted for , an unrealistically high value. We will show in the next sections that is positive and close to zero, resulting in a very small-sized blueshift region, compared to the overall size of the Universe. Similar plots can be obtained for the ratio expressed in terms of the radial distance (see  for details), which also suggest the existence of a blueshift region localized around the observer’s position, that is, the Earth could be surrounded by a natural blueshift region, extending at least over the solar system region. This might be the origin of the Pioneer anomaly. (Obviously, the Earth’s observer is not located at any privileged position. The same cosmological evolution described by CC would be seen by any other observer in the Universe, provided that the local values of the cosmological parameters and are the same. In our previous work ([28, 29]), we have suggested that might play the role of a universal time, so that for a certain value of this parameter the evolution of the Universe would look the same for any observer. In this way, conformal cosmology does not violate the cosmological principle, which postulates a homogeneous and isotropic Universe.)
Before we proceed to analyze this possible explanation for the anomaly, we recall a few more results obtained in our second paper . Since we have closed-form expressions for our scale factor , in (4) and (5), it is straightforward to obtain the Hubble parameter () and the deceleration parameter () as a function of time or redshift . For the case, we obtained  and, in particular for or ,
The signs of the quantities in (7) and (8) can be explained by considering again the red-solid curve in Figure 1, which represents the ratio , or equivalently , over different cosmological epochs. This bell-shaped curve was plotted for a positive value of and shows a local blueshift area in the “past” evolution of the Universe, extending back to (represented by the square point in Figure 1) or , for the look-back time at which the redshift () starts being observed. The red curve has a maximum at or (we can also find or ), and it is evidently symmetric around this point of maximum expansion of the Universe.
Therefore, for each value of , that is, for each value of , we have two corresponding values of the Hubble parameter (except at the maximum, for , where ). The two related points on the curve, at the same redshift level, will have equal and opposite expansion rates. This yields the double sign in the previous expressions for , when given as a function of . This argument applies also to the case, corresponding to the current time , at which is negative, showing that the Universe is already in a contracting phase. (This is also a consequence of the signs of our conformal parameters, in particular the positive value of . Our estimate of will be given in Section 3, but we recall that Mannheim has independently evaluated as a small but positive quantity (), by fitting rotational velocity curves for several spiral galaxies, using conformal gravity . If was to have a negative value, we would still be in an expanding phase of the Universe.) As discussed above, the same value can also refer to the time in the past () at which we start observing the cosmological redshift, with , a positive quantity. This analysis does not contradict the current astrophysical estimates of as a positive quantity. They are based on redshift observations of light coming from galaxies at times in the past ; therefore, what is denoted by in standard cosmology should be actually indicated as , again a positive quantity related to the expanding phase of the Universe. The same analysis can be done in terms of radial distances . The blueshift region would extend from up to a distance given by where is the distance at which we start observing the cosmological redshift. In general, the slope of the red-solid curve in Figure 1 is related to the value of the Hubble parameter at that point, while its curvature is connected to the deceleration parameter, through the expressions given above.
In particular, following (8), the slope of the plot and its curvature at current time are basically connected to our two fundamental parameters and . In the next section, we will show that the slope of the red-solid plot at is closely related to the value of the Pioneer anomalous acceleration , which can therefore be used to determine . Similarly, the curvature of the plot at will be related to the rate of change of the anomalous acceleration (i.e., the “jerk” ) and will be used to determine the value of our other parameter .
We conclude this section by noting that the values of our parameters ( and ) could be derived directly from standard cosmological observations, in view of (8). Using the current best estimate of  and the positive sign in (8), we obtain The direct determination of is more difficult, since the deceleration parameter is not known explicitly. In , we based our analysis on recent luminosity data for type-Ia Supernovae, obtaining an estimate of , but this analysis needs to be confirmed by further studies.
3. The Pioneer Anomaly
In the previous section, we briefly reviewed our conformal cosmology and outlined the reasons why we consider the solution as a possible description of the evolution of the Universe. This solution can explain the observed cosmological redshift, but it requires the existence of a blueshift region in the immediate vicinity of our current spacetime position in the Universe.
This could be a serious problem for our model, since we do not observe blueshift of nearby astrophysical objects except for the one caused by the peculiar velocities of nearby galaxies, presumably due to standard Doppler shift. However, as already mentioned in Section 1, experimental evidence of a local region of blueshift might come from the analysis of the Pioneer anomaly ([1–4, 30, 35–42]).
This is a small frequency drift (blueshift), observed analyzing the navigational data of the Pioneer 10-11 spacecraft, received from distances between 20 and 70 AU (astronomical units) from the Sun, while these spacecraft were exploring the outer solar system. This anomaly is usually reported as a positive rate of change of the signal frequency, (blueshift), resulting in a frequency drift of about 1.5 Hz every 8 years, or as an almost constant sunward acceleration, , or even as a “clock acceleration” . More precisely ([1, 35]),
An attempt was made to detect such anomaly also in the radiometric data from other spacecraft traveling at the outskirts of the solar system, such as the Galileo and Ulysses missions . In the case of Galileo, the effects of solar radiation made such detection impossible, while for Ulysses a possible anomalous acceleration cm/s2 was seen in the data. Other spacecraft, such as the New Horizons mission to Pluto, launched in 2006, might provide new data in the near future. These discoveries prompted a complete reanalysis of all the historical navigational data of these space missions, which is currently underway ([1, 31, 37–39, 42, 43]) and will be completed in the near future . This new analysis will try to determine additional characteristics of the anomaly, such as its precise direction, the possible temporal and spatial variations, and its dependence on heliocentric or geocentric distance. A future dedicated mission is also being proposed ([45–48]) to test directly this puzzling phenomenon.
Currently, the origin and nature of this anomaly remains unexplained; all possible sources of systematic errors have been considered ([1, 4, 35, 36, 39, 43, 49]), but they cannot fully account for the observed effect. The current focus of conventional explanations of the anomaly seems to be the thermal recoil force, that is, anisotropically emitted thermal radiation, originating from the spacecraft four radioisotope thermoelectric generators (RTGs), which can contribute significantly to the measured acceleration. The natural decay of the radioactive material in the RTGs, the aging of the thermocouples in the system, and other effects all contribute to the decrease of the total thermal power during the spacecraft life. This might explain the decrease overtime of the measured Pioneer acceleration (in absolute value), that is, the negative “jerk” , already seen in the early Pioneer data ([1, 30, 44]).
Although the anomaly can be caused by these standard physical effects, we will try in the following to explain its origin by using the cosmological model outlined in the previous section. The phenomenology of the Pioneer anomaly is related to a complex exchange of radiometric signals between the tracking stations on Earth (of the deep space network (DSN)) and the spacecraft, using S-band Doppler frequencies (1.55–5.20 GHz). Typically, an uplink signal is sent from the DSN to the spacecraft at a frequency of 2.11 GHz, based on a very stable hydrogen maser system, then an S-band transponder on board the spacecraft applies an exact and fixed turn-around ratio of 240/221 to the uplink signal, so that the Pioneer returns a downlink signal at a slightly different frequency of about 2.29 GHz, to avoid interference with the uplink one.
This procedure is known as a two-way Doppler coherent mode and allows for very precise tracking of the spacecraft, since the returning signal is directly compared to the original one. On the contrary, a one-way Doppler signal (with a fixed signal source on the spacecraft, whose frequency cannot be monitored for accuracy) is less effective. This type of tracking system added to the propulsion and navigational characteristics of the Pioneer spaceship (especially the presence of a spin-stabilization system) resulted in a very good acceleration sensitivity of about 10−8 cm/s2, once the influence of solar radiation pressure can be neglected (for distances AU from the Sun).
The DSN station acquires the downlink signal after a time delay ranging from a few minutes to some hours, depending on the distance involved, and compares it to the reference frequency to determine the Doppler shift due to the actual motion of the spacecraft. The navigational software can also model with great precision the expected frequency of the signal returned from the Pioneer, which should coincide with the one observed on Earth. As already mentioned, a discrepancy was found, corresponding to the values in (11), whose origin cannot be traced to any systematic effect due to either the performance of the spacecraft or the theoretical modeling of its navigation.
The Pioneer anomaly was first reported ([2–4]) as an almost constant value of the anomalous acceleration, with temporal and space variation of within 10%, over a range of heliocentric distances ~20–70 AU, and possibly at even closer distances AU, so that we will concentrate first on the average value of and later on its variation with time and distance. In our view, the Pioneer phenomenology represents the most basic experiment we could perform in order to check if the cosmic evolution is really affecting the frequency of electromagnetic radiation emitted and observed at different spacetime locations, following (4) and (5).
In the standard analysis of the Pioneer anomaly, the signal coming back to Earth is affected by the relativistic Doppler effect. Following this model, will be the frequency of the expected signal and will be related to the signal reference frequency GHz (for the uplink signal in a two-way system) by the standard relativistic Doppler formula (see equation 2.2.2 in ): where is the spacecraft radial velocity and the approximation on the right-hand side holds to first order in .
Since we have a two-way system, the Doppler shift involved is actually double, so we can use the previous equation but with , where is the expected velocity of the spacecraft, according to the theoretical navigation model, at time , when the spaceship receives and immediately retransmits the signal. We use here a time variable which can be simply considered the elapsed time since the spacecraft launch ( at ), and then later we will simply identify with our cosmological look-back time in (6). With this radial velocity, (12) to first order in becomes and this frequency is expected to be observed with high precision, due to the reported excellent navigational control of the spacecraft.
On the contrary, a different frequency is observed, , involving an additional unexplained blueshift: this is the Pioneer anomaly. Following (11), the frequency difference is reported as where the factor of two in the first line of the previous equation is due to the two-way system. We also remark here that several of the cited references adopt a rather confusing “DSN sign convention” for the frequency difference in (14) (see [1, 4, 38] and of ), resulting in a change of sign in most of their equations. We prefer to use here our definition of as given in the previous equation.
The anomalous acceleration is introduced as an alternative way of describing the effect, although, in our view, it does not correspond to a real spacecraft acceleration. As in (13), we can write the observed frequency to first order in as where the “observed” velocity of the spacecraft refers to the time of interest . Combining together the last three equations, we can write the frequency difference as
These frequency differences (also called frequency residuals in the literature cited) are therefore equivalent to the corresponding velocity residuals (), and they are usually plotted as a function of the elapsed time , showing an almost linear increase with time of these residuals, which is the essence of the Pioneer anomaly (see, e.g., Figure 5.2 in ). The Pioneer anomalous acceleration can be defined as the rate of change of the velocity residuals, related to the corresponding rate of change of the frequency residuals, in view of (16). Therefore, if we define , the Pioneer acceleration can be related to the frequency differences which are more significant quantities in our analysis. We will assume that these frequency differences are intrinsically due to the different locations of the spacecraft (at position ) and of the Earth’s observer (at ). Therefore, we identify the reference frequency in (13) with and the similar quantity in (15) with . Then, we subtract (13) from (15): where the common factor of two in all the parts of the previous equation was added again because of the two-way effect, which has to be included also in our gravitational blueshift model. The velocities and from (13) and (15) are assumed to be the same, so that the common factor is close to unity and can be neglected, since the average Pioneer speed is . We also identified with the Earth reference frequency and used our fundamental equation (4) and (5), case, to first order in . (The elapsed time for the Pioneer spacecraft missions is of the order of a few years (); we can assume ; therefore, .) Using these results, (17) simplifies as follows: in view also of our evaluation of (a negative quantity) from (8).
This result immediately explains the often cited “numerical coincidence,” that is, the simple relation between the Pioneer acceleration and the standard (positive) Hubble constant, with the correct negative sign for both quantities in (19), in view of our previous discussion of the sign of . Equation (19) can also be used to determine and (as a positive quantity), using the reported value of from (11):
The value of (considered measured at the current time , even if the Pioneer data are a few years old) is close to our first direct estimate in (10), and the corresponding value of the Hubble constant is close to the value of standard cosmology. We remark here again that our model fully explains the reason of this “numerical coincidence” and provides also the correct signs for all the quantities involved. (The numerical “coincidence” between the Hubble constant and the value of the Pioneer acceleration divided by was noticed immediately after the discovery of the Pioneer effect and prompted many speculations and different explanations. This coincidence is even more striking if one uses the value cited in  as the experimental value for Pioneer 10 data before systematics, , thus obtaining and .)
Following (17)–(19) and the related discussion, we can generalize our expression of the Pioneer acceleration, as a function of time : with . (Although is the elapsed time since the spacecraft launch, it is treated here as equivalent to a look-back time because the Pioneer is moving toward increasing distances , therefore corresponding to increased look-back times in our original redshift interpretation.) In particular, by using the previous equation and taking another time derivative, it is easy to derive the “jerk” and its value in the limit for : expressed in terms of our fundamental parameters and . The current value of in the last equation is positive (for small values of ), but the Pioneer acceleration, as in (11) or (19), is considered negative in this paper so that a positive jerk means that the absolute value of will decrease for increasing times or radial distances, which is indeed shown in the early Pioneer data, as it was already mentioned at the beginning of this section.
In Figure 2, we illustrate the early Pioneer 10/11 data, as originally reported in , where the absolute value of the Pioneer acceleration is plotted as a function of the radial heliocentric distance in AU. The red-dashed horizontal line and the green-dotted lines represent, respectively, the value of and the related 1-sigma error range quoted in (11). The first three data points for Pioneer 11, at smaller distances, lie outside the considered error range probably because the anomalous acceleration was masked by solar radiation or other effects. We will not include these first three data points in our subsequent discussion. We will concentrate our analysis on either just the Pioneer 10 data points or the combination of data points for both spacecraft, but within the 1-sigma error range (“error range data” in the following).
These two sets of data clearly show a possible decrease of the Pioneer anomaly (in the absolute value ) with increasing heliocentric distance. The black (dash-dotted) line and the blue-solid line in the same figure represent linear fits for the Pioneer 10 and the error range data, respectively, both of them indicating a decrease of .
If our conformal cosmology is the origin of the Pioneer anomaly, and not the thermal recoil force mentioned at the beginning of this section, our “jerk” equation (22) will explain the decrease of and can also be used to determine our second parameter .
We computed the slopes of our two linear fits in Figure 2 and used them as (positive) values of in (22), together with the value from (20). (The radial distances of the data plotted in Figure 2 were converted into elapsed times , by using a simple approximation: , where is the average Pioneer speed. From the original data (available from the NASA website at: http://cohoweb.gsfc.nasa.gov/helios/), we estimated km/s, km/s, and used an average km/s when combining data for both spacecraft.) Solving (22) for , we obtain and these values for are very close to the one we obtained in  (), which was based solely on the analysis of type-Ia Supernovae data.
Another type of analysis is illustrated in Figure 3. The Pioneer 10/11 data, the standard value of , and the related error range are the same as in the previous figure, but this time we used the generalized expression of in (21) to fit the data within the error range. We allowed both quantities and to be free parameters in our fitting procedure, and we converted the elapsed time in (21) into the radial distance by using the approximation , where is the average Pioneer speed, as it was done also for the data in the previous figure. The radial distance should be more properly identified with the geocentric distance of the spacecraft, rather than the heliocentric one, since should be the distance from the Earth observer. We also performed fits using the geocentric distance, but the results were very similar to those obtained by using heliocentric distances, so we will not include them in the following analysis.
Again, in Figure 3, we used the expression in (21) to fit the data, although the fitting curves appear almost as straight lines in this figure. The first conformal cosmology fit, illustrated by the black (dash-dotted) curve, was obtained by using only the Pioneer 10 data and yielded the following values of the parameters: The second fit (blue-solid curve) was obtained by using all the data within the error range (again omitting the first three Pioneer 11 data points) and produced the following results:
Comparing the results in the last two equations with those for in (23), obtained with a fixed as in (20), we can see that all the values of our parameters are in agreement. In particular, from the different analyses, we consistently obtain and , where the different values depend on the Pioneer data being used. As already remarked, the values for quoted above are also close to the one we obtained in  (), based on type-Ia Supernovae data. In the next section, we will discuss our results and compare them to the current limits of standard gravity in the solar system.
As already mentioned at the beginning of this paper, a new analysis of extended Pioneer data has recently appeared , and these new results will also be discussed in the next section. However, complete new data points from this extended analysis are not yet available but will be published in the future . Due to this reason, we have based our analysis in this section only on the early Pioneer data which were available at the time of our study.
4. Discussion of Our Results and Conclusions
In the previous sections, we discussed how conformal cosmology provides a natural explanation for the Pioneer anomalous acceleration, in both magnitude and direction (i.e., the negative sign of the radial acceleration). We also explained the “numerical coincidence,” connecting with the Hubble constant, and the observed decrease with heliocentric distance of , related to the Pioneer jerk . Although the Pioneer data are still not very accurate, our analysis consistently indicated that our conformal parameters are approximately given by and (see (10), (20), and (23)–(25)). In this final section, we will discuss the implications of the values of our parameters in relation to other studies in the field.
We first remark that a new analysis of rotational velocity data for spiral galaxies, based on conformal gravity, has recently appeared ([59, 60]) improving the original work on the subject ([32, 33]). This new study uses the full line element of conformal gravity in (1)-(2), including the effects of the quadratic , which were previously neglected, thus obtaining a global gravitational potential of cosmological origin. In addition to this, a local gravitational potential is obtained by integrating over the visible galactic mass distribution a gravitational potential per unit solar mass of the form . The two potentials, global and local, are then combined together to model the rotational motion of galaxies. The fits to galactic rotation data ([59, 60]), performed without any dark matter contribution, show a remarkable success of conformal gravity, even at the largest distances from the galactic centers, where the quadratic becomes important and comparable to the linear term . Mannheim and collaborators ([59, 60]) were then able to determine the values of the global universal parameters as and . The related terms of the global gravitational potential were associated respectively to the cosmological background and to cosmological inhomogeneities. The local parameter was also evaluated as .
The values of the dimensionful parameters and obtained through this analysis of galactic rotation curves are somewhat different from our values, reported in this paper or in our previous work  (cm−1 and cm−2). This difference could be due, as we explained in , to a possible redefinition of the luminosity distance and other distance indicators, which might affect even the radial distances (from the galactic centers) which are employed in the galactic rotation analysis.
However, it is instructive to compute the dimensionless parameter, using Mannheim’s values and , because this dimensionless constant should not be affected by a revision of the cosmological distances. As explained at the beginning of Section 2, the parameters , , , and are related through and also so that we obtain Therefore, the conformal gravity analysis by Mannheim and collaborators suggests a Universe, consistent with our cosmological model and also a value of , close to our quoted values of .
Conformal gravity considers local gravitational effects as being due to the local potential or simply to the potential , where is the radial distance from the center of our solar system. Since the value of the local constant is very small, the modifications to standard dynamics of the solar system are negligible. (For example, the ratio between the conformal gravitational potential and the standard Newtonian term at a heliocentric distance of 1 AU is ~10−20, while, at a distance of 100 AU (outer solar system), the same ratio is ~10−16. The galactic potential, related to another conformal gravity term of the form ( distance from galactic center) would generate only a very small tidal force effect on the solar system. This correction would behave as the ratio of the distance of a planet or satellite from the sun divided by the distance from the solar system to the center of the galaxy . Therefore, the “conformal gravity force” is negligible, compared to the standard Newtonian one, over the whole solar system region.) Therefore, conformal gravity is not in any way in contradiction with the very stringent limits on alternative gravity theories imposed by studies of planetary ephemerides or other solar system observations ([62, 63]).
In particular, more recent studies have focused their attention on the critical issue of the influence that a gravitational Pioneer anomaly acceleration would have on the motion of bodies in the solar system, such as inner and outer planets, comets, and asteroids ([64–75]). These studies indicate that an anomalous acceleration of gravitational origin would have a significant effect on the motion of these solar system objects, but the resulting orbital anomalies were not detected in the latest observations, therefore ruling out the gravitational origin of the Pioneer anomaly. As already mentioned above, we also believe that the dynamical corrections due to conformal gravity are in fact negligible at the solar system level. Therefore, no significant changes in the orbits of bodies in the solar system should be detected, in agreement with the cited references.
To further clarify the issue, in our analysis of the Pioneer anomaly, we used the reported values of the anomalous acceleration to determine the cosmological parameters, simply because such was the way these data were reported in the literature cited. However, it should be clear from the discussion in Section 3 that we explain the Pioneer anomaly in terms of our cosmological-gravitational blueshift, based on the global values of the parameters , , and . In this view, there is no real dynamic acceleration of the Pioneer spacecraft (or of any other object in the solar system) oriented toward the Sun, due to some new gravitational force or modification of existing gravity, except for the tiny corrections coming from local conformal gravity mentioned above. In fact, in our analysis, we assume that there is no difference between the two velocities and in (13) and (15); therefore, the anomalous acceleration defined as with is actually zero.
In this way, we also overcome the original objection, reported in , that “the anomalous acceleration is too large to have gone undetected in planetary orbits, particularly for Earth and Mars,” since “NASA’s Viking mission provided radio-ranging measurements  to an accuracy of about 12 m,” which should have shown the effect of the anomalous acceleration on the orbits of these two planets.
In our view, precision ranging measurements with radio signals or lasers, based on the round-trip travel time from Earth to other bodies in the solar system, would not show any anomalous effect because the speed of light is not affected by our cosmological model and the corrections to the dynamics of the solar system due to conformal gravity are negligible.
On the contrary, we would observe an effect similar to the anomalous acceleration for a spacecraft, a planet, or any other object in the solar system, if we were to study its motion through Doppler frequency ranging, because of the intrinsic differences in frequency or wavelength for light emitted at different spacetime positions, due to our cosmological model.
A similar discussion can be done regarding possible explanations of the Pioneer anomaly of cosmological origin. Recent studies based on the standard Friedmann-Lemaitre metric ([77, 78]) have shown that cosmological effects fail to account for the anomaly by several orders of magnitude. Again, our study does not propose dynamical corrections to the orbits, due to conformal gravity or to a conformal cosmological model, since these corrections would be negligible in the solar system, as are those due to standard cosmology. Our model simply assumes the existence of a local blueshift region, which is able to affect the frequencies of the signals exchanged between the spacecraft and Earth.
The size of the local blueshift region, which in our model is responsible for the frequency differences, can be easily estimated by using (9) and the values of our parameters. For example, using the values from our conformal cosmology fits in (24) and (25), we obtain , corresponding to a distance comparable to the one between Earth and the nearest bright stars (which is about ). This blueshift region would extend well beyond the solar system but would cover a small portion of our galaxy, since .
The maximum blueshift effect would be seen at and would correspond to a , a very small value. Therefore, the blueshift region and the related effects are so small that they cannot be practically observed in the radiation spectrum of stars or other radiation-emitting objects within this region. These effects are only small corrections to the Doppler signals coming from the Pioneer or other similar spacecraft.
We also want to compare our estimates of the rate of change of the anomalous acceleration (i.e., the jerk ) with those presented by independent verifications of the Pioneer anomaly (see review in ). The first of these studies was performed by Markwardt , who reviewed data for Pioneer 10 and reported , with . (We prefer to report here, as also done in the rest of the paper, the anomalous acceleration as a negative quantity and the related jerk as a positive quantity. Some of the papers in the literature adopt the opposite sign convention, which might generate some confusion.) Using Markwardt values in (19) and (22), we obtain and , consistent with our values in (20) and (23)-(24) for Pioneer 10.
The second independent study was done by Toth  and reported results separately for the two spacecraft. From Toth’s results for Pioneer 10 (, ), we compute and . Using instead Toth’s results for Pioneer 11 (, ), we obtain and , and all these results are also consistent with those discussed in Section 3.
Finally, we wish to comment briefly on the recent analysis of extended Pioneer data by Turyshev and collaborators . This new study confirmed the anomalous acceleration and its temporal decrease in absolute value () using data spans more than twice as long as those used in previous studies. To analyze the temporal behavior, different models were used including a linear fit similar to the one which was employed in our analysis. Their results were reported separately for the two spacecrafts as follows (using our sign conventions).
Pioneer 10: cm/s2, , from which we compute and . Pioneer 11: , , from which we obtain and . Therefore, these results are also in agreement with those discussed in Section 3. A more detailed comparison of the latest Pioneer results  with those obtained using conformal cosmology will be presented in a future study on the subject.
In conclusion, the detailed analysis of the Pioneer anomaly presented in this work has indicated that our conformal cosmology might be the origin of this effect, while conformal gravity alone cannot account for the anomalous acceleration of the spacecraft. If our analysis is correct, it explains naturally the numerical coincidence between the Pioneer acceleration and the Hubble constant, including the signs of these quantities. In addition, we confirm our previous evaluations of the cosmological parameters, and , also in agreement with independent evaluations. Further studies will be needed when the reanalysis of all the historical navigational data for the Pioneer spacecraft will be completed by Turyshev and collaborators, and new data will be publicly available.
This work was supported by a grant from the Frank R. Seaver College of Science and Engineering, Loyola Marymount University. The author would like to acknowledge Dr. S. Turyshev and Dr. P. Mannheim for very useful discussions and advice on the subject. The author also thanks the anonymous reviewers for the useful comments which helped improve the final version of the paper.
- S. G. Turyshev and V. T. Toth, “The pioneer anomaly,” Living Reviews in Relativity, vol. 13, no. 4, pp. 9–175, 2010.
- J. D. Anderson, P. A. Laing, E. L. Lau, A. S. Liu, M. M. Nieto, and S. G. Turyshev, “Indication, from pioneer 10/11, Galileo, and Ulysses data, of an apparent anomalous, weak, long-range acceleration,” Physical Review Letters, vol. 81, no. 14, pp. 2858–2861, 1998.
- S. G. Turyshev, J. D. Anderson, P. A. Laing, E. L. Lau, A. S. Liu, and M. M. Nieto, “The apparent anomalous, weak, long-range acceleration of pioneer 10 and 11,” http://arxiv.org/abs/gr-qc/9903024.
- J. D. Anderson, P. A. Laing, E. L. Lau, A. S. Liu, M. M. Nieto, and S. G. Turyshev, “Study of the anomalous acceleration of Pioneer 10 and 11,” Physical Review D, vol. 65, no. 8, Article ID 082004, pp. 820041–8200450, 2002.
- C. Lämmerzahl, O. Preuss, and H. Dittus, “Is the physics within the Solar system really understood?” in Lasers, Clocks and Drag-Free Control: Exploration of Relativistic Gravity in Space, H. Dittus, C. Lämmerzahl, and S. G. Turyshev, Eds., vol. 349 of Astrophysics and Space Science Library, p. 75, 2008.
- J. D. Anderson and M. M. Nieto, “Astrometric solar-system anomalies,” in Relativity in Fundamental Astronomy: Dynamics, Reference Frames, and Data Analysis, S. A. Klioner, P. K. Seidelmann, and M. H. Soffel, Eds., Proceedings IAU Symposium no. 261, pp. 189–197, 2010.
- G. A. Krasinsky and V. A. Brumberg, “Secular increase of astronomical unit from analysis of the major planet motions, and its interpretation,” Celestial Mechanics and Dynamical Astronomy, vol. 90, no. 3-4, pp. 267–288, 2004.
- J. D. Anderson, J. K. Campbell, and M. M. Nieto, “The energy transfer process in planetary flybys,” New Astronomy, vol. 12, no. 5, pp. 383–397, 2007.
- J. D. Anderson, J. K. Campbell, J. E. Ekelund, J. Ellis, and J. F. Jordan, “Anomalous orbital-energy changes observed during spacecraft flybys of Earth,” Physical Review Letters, vol. 100, no. 9, Article ID 091102, 2008.
- S. G. Turyshev and V. T. Toth, “The puzzle of the flyby anomaly,” Space Science Reviews, vol. 148, no. 1–4, pp. 169–174, 2009.
- L. Iorio, “The recently determined anomalous perihelion precession of saturn,” Astronomical Journal, vol. 137, no. 3, pp. 3615–3618, 2009.
- L. Iorio, “The perihelion precession of Saturn, planet X/Nemesis and MOND,” The Open Astronomy Journal, vol. 3, no. 1, 14 pages, 2010.
- L. Iorio, “On the anomalous secular increase of the eccentricity of the orbit of the Moon,” Monthly Notices of the Royal Astronomical Society, vol. 415, no. 2, pp. 1266–1275, 2011.
- L. Iorio, “An empirical explanation of the anomalous increases in the astronomical unit and the lunar eccentricity,” Astronomical Journal, vol. 142, no. 3, article 68, 2011.
- E. V. Pitjeva, “EPM ephemerides and relativity,” in Relativity in Fundamental Astronomy: Dynamics, Reference Frames, and Data Analysis, S. A. Klioner, P. K. Sei-delmann, and M. H. Soffel, Eds., Proceedings IAU Symposium no. 261, pp. 170–178, Cambridge University Press, Cambridge, UK, 2010.
- E. V. Pitjeva and N. P. Pitjev, “Estimations of changes of the Sun's mass and the gravitation constant from the modern observations of planets and spacecraft,” http://arxiv.org/abs/1108.0246.
- D. Grumiller, “Model for gravity at large distances,” Physical Review Letters, vol. 105, no. 21, Article ID 211303, 2010.
- S. Carloni, D. Grumiller, and F. Preis, “Solar system constraints on Rindler acceleration,” Physical Review D, vol. 83, no. 12, Article ID 124024, 2011.
- D. Grumiller and F. Preis, “Rindler force at large distances,” http://arxiv.org/abs/1107.2373.
- L. Iorio, “Solar system constraints on a Rindler-type extra-acceleration from modified gravity at large distances,” Journal of Cosmology and Astroparticle Physics, vol. 2011, no. 5, article 019, 2011.
- L. Iorio, “Impact of a Pioneer/Rindler-type acceleration on the Oort cloud,” http://arxiv.org/abs/1108.0409.
- H. Weyl, “Reine Infinitesimalgeometrie,” Mathematische Zeitschrift, vol. 2, pp. 384–411, 1918.
- H. Weyl, “Gravitation und Elektrizität,” Sitzungsberichte der Preussischen Akademie der Wissenschaften. Physikalisch-Mathematische Klasse, vol. 1, pp. 465–480, 1918.
- H. Weyl, “Eine Neue Erweiterung der Relativitaetstheorie,” Annalen der Physik, vol. 59, pp. 101–103, 1919.
- P. D. Mannheim and D. Kazanas, “Exact vacuum solution to conformal Weyl gravity and galactic rotation curves,” Astrophysical Journal, vol. 342, p. 635, 1989.
- D. Kazanas and P. D. Mannheim, “General structure of the gravitational equations of motion in conformal weyl gravity,” Astrophysical Journal, Supplement Series, vol. 76, no. 2, pp. 431–453, 1991.
- P. D. Mannheim, “Alternatives to dark matter and dark energy,” Progress in Particle and Nuclear Physics, vol. 56, no. 2, pp. 340–445, 2006.
- G. U. Varieschi, “A kinematical approach to conformal cosmology,” General Relativity and Gravitation, vol. 42, no. 4, pp. 929–974, 2010.
- G. U. Varieschi, “Kinematical conformal cosmology: fundamental parameters from astrophysical observations,” ISRN Astronomy and Astrophysics, vol. 2011, Article ID 806549, 24 pages, 2011.
- M. M. Nieto and J. D. Anderson, “Using early data to illuminate the Pioneer anomaly,” Classical and Quantum Gravity, vol. 22, no. 24, pp. 5343–5354, 2005.
- S. G. Turyshev, V. T. Toth, J. Ellis, and C. B. Markwardt, “Support for temporally varying behavior of the Pioneer anomaly from the extended Pioneer 10 and 11 doppler data sets,” Physical Review Letters, vol. 107, no. 8, Article ID 081103, 2011.
- P. D. Mannheim, “Linear potentials and galactic rotation curves,” Astrophysical Journal, vol. 419, no. 1, pp. 150–154, 1993.
- P. D. Mannheim, “Are galactic rotation curves really flat?” Astrophysical Journal, vol. 479, no. 2, pp. 659–664, 1997.
- K. Nakamura and Particle Data Group, “Review of particle physics,” Journal of Physics A, vol. G37, Article ID 075021, 2010.
- J. D. Anderson, E. L. Lau, S. G. Turyshev, P. A. Laing, and M. M. Nieto, “Search for a standard explanation of the pioneer anomaly,” Modern Physics Letters A, vol. 17, no. 14, pp. 875–885, 2002.
- S. G. Turyshev, M. M. Nieto, and J. D. Anderson, “Study of the Pioneer anomaly: a problem set,” American Journal of Physics, vol. 73, no. 11, pp. 1033–1044, 2005.
- S. G. Turyshev, V. T. Toth, L. R. Kellogg, E. L. Lau, and K. J. Lee, “A study of the pioneer anomaly: new data and objectives for new investigation,” International Journal of Modern Physics D, vol. 15, no. 1, pp. 1–55, 2006.
- V. T. Toth and S. G. Turyshev, “The Pioneer anomaly: seeking an explanation in newly recovered data,” Canadian Journal of Physics, vol. 84, no. 12, pp. 1063–1087, 2006.
- V. T. Toth and S. G. Turyshev, “Pioneer anomaly: evaluating newly recovered data,” in Proceedings of the 3rd Mexican Meeting on Mathematical and Experimental Physics, vol. 977 of AIP Conference Proceedings, pp. 264–283, September 2007.
- M. M. Nieto and J. D. Anderson, “Search for a solution of the Pioneer anomaly,” Contemporary Physics, vol. 48, no. 1, pp. 41–54, 2007.
- V. T. Toth and S. G. Turyshev, “Thermal recoil force, telemetry, and the Pioneer anomaly,” Physical Review D, vol. 79, no. 4, Article ID 043011, 2009.
- S. G. Turyshev and V. T. Toth, “The pioneer anomaly in the light of new data,” Space Science Reviews, vol. 148, no. 1–4, pp. 149–167, 2009.
- S. G. Turyshev and V. T. Toth, “Physics engineering in the study of the pioneer anomaly,” http://arxiv.org/abs/0710.0191.
- S. Turyshev, private communication, 2010.
- J. D. Anderson, M. M. Nieto, and S. G. Turyshev, “A mission to test the Pioneer anomaly,” International Journal of Modern Physics D, vol. 11, no. 10, pp. 1545–1551, 2002.
- M. M. Nieto, S. G. Turyshev, and J. D. Anderson, “The pioneer anomaly: the data, its meaning, and a future test,” in Proceedings of the 2nd Mexican Meeting on Mathematical and Experimental Physics, vol. 758 of AIP Conference Proceedings, pp. 113–128, September 2004.
- H. Dittus, S. G. Turyshev, C. Lämmerzahl et al., “A mission to explore the Pioneer anomaly,” in Proceedings of the 39th ESLAB Symposium: Trends in Space Science and Cosmic Vision 2020, ESA Spec.Publ. 588, pp. 3–10, April 2005.
- S. G. Turyshev, M. M. Nieto, and J. D. Anderson, “A route to understanding of the pioneer anomaly,” in Proceedings of the 22nd Texas Symposium on Relativistic Astrophysic, December 2004, paper no. 0310.
- M. M. Nieto and S. G. Turyshev, “Finding the origin of the Pioneer anomaly,” Classical and Quantum Gravity, vol. 21, no. 17, pp. 4005–4023, 2004.
- O. Bertolami, F. Francisco, P. J. S. Gil, and J. Páramos, “Thermal analysis of the Pioneer anomaly: a method to estimate radiative momentum transfer,” Physical Review D, vol. 78, no. 10, Article ID 103001, 2008.
- B. Rievers, C. Lämmerzahl, M. List, S. Bremer, and H. Dittus, “New powerful thermal modelling for high-precision gravity missions with application to Pioneer 10/11,” New Journal of Physics, vol. 11, Article ID 113032, 2009.
- O. Bertolami, F. Francisco, P. J. S. Gil, and J. Páramos, “Estimating radiative momentum transfer through a thermal analysis of the pioneer anomaly,” Space Science Reviews, vol. 151, no. 1–3, pp. 75–91, 2010.
- B. Rievers, S. Bremer, M. List, C. Lämmerzahl, and H. Dittus, “Thermal dissipation force modeling with preliminary results for Pioneer 10/11,” Acta Astronautica, vol. 66, no. 3-4, pp. 467–476, 2010.
- B. Rievers, C. Lämmerzahl, and H. Dittus, “Modeling of thermal perturbations using raytracing method with preliminary results for a test case model of the pioneer 10/11 radioisotopic thermal generators,” Space Science Reviews, vol. 151, no. 1–3, pp. 123–133, 2010.
- B. Rievers and C. Lämmerzahl, “High precision thermal modeling of complex systems with application to the flyby and Pioneer anomaly,” Annalen der Physik (Leipzig), vol. 523, no. 6, pp. 439–449, 2011.
- F. Francisco, O. Bertolami, P. J. S. Gil, and J. Páramos, “Modelling the reflective thermal contribution to the acceleration of the Pioneer spacecraft,” http://arxiv.org/abs/1103.5222.
- S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity, John Wiley & Sons, New York, NY, USA, 1972.
- C. B. Markwardt, “Independent confirmation of the pioneer 10 anomalous acceleration,” http://arxiv.org/abs/gr-qc/0208046.
- P. D. Mannheim and J. G. O'Brien, “Impact of a global quadratic potential on galactic rotation curves,” Physical Review Letters, vol. 106, no. 12, Article ID 121101, 2011.
- P. D. Mannheim and J. G. O'Brien, “Fitting galactic rotation curves with conformal gravity and a global quadratic potential,” http://arxiv.org/abs/1011.3495.
- P. Mannheim, private communication, 2011.
- E. M. Standish, “Planetary and Lunar Ephemerides: testing alternate gravitational theories,” in Recent Developments in Gravitation and Cosmology, A. Macias, C. Lämmerzahl, and A. Camacho, Eds., vol. 977 of AIP Conference Proceedings, pp. 254–263, American Institute of Physics, Melville, NJ, USA, 2008.
- E. M. Standish, “Testing alternate gravitational theories,” in Fundamental Astronomy: Dynamics, Reference Frames, and Data Analysis, Proceedings of the International Astronomical Union, S. A. Klioner, P. K. Seidelmann, and M. H. Soffel, Eds., vol. 261 of IAU Symposium, pp. 179–182, 2010.
- L. Iorio and G. Giudice, “What do the orbital motions of the outer planets of the Solar System tell us about the Pioneer anomaly?” New Astronomy, vol. 11, no. 8, pp. 600–607, 2006.
- G. L. Page, D. S. Dixon, and J. F. Wallin, “Can minor planets be used to assess gravity in the outer solar system?” Astrophysical Journal, vol. 642, no. 1, pp. 606–614, 2006.
- K. Tangen, “Could the Pioneer anomaly have a gravitational origin?” Physical Review D, vol. 76, no. 4, 2007.
- L. Iorio, “Can the pioneer anomaly be of gravitational origin? A phenomenological answer,” Foundations of Physics, vol. 37, no. 6, pp. 897–918, 2007.
- L. Iorio, “Jupiter, Saturn and the Pioneer anomaly: a planetary-based independent test,” Journal of Gravitational Physics, vol. 1, no. 1, pp. 5–8, 2007.
- J. F. Wallin, D. S. Dixon, and G. L. Page, “Testing gravity in the outer solar system: results from trans-neptunian objects,” Astrophysical Journal, vol. 666, no. 2, pp. 1296–1302, 2007.
- L. Iorio, “The Lense-Thirring effect and the Pioneer anomaly: solar system tests,” in Proceedings of the the 11th Marcel Grossmann Meeting on Recent Developments in Theoretical and Experimental General Relativity, Gravitation and Relativistic Field Theorie, H. Kleinert, R.T. Jantzen, and R. Ruffini, Eds., pp. 2558–2560, World Scientific, 2008.
- A. Fienga, J. Laskar, P. Kuchynka, C. Leponcin-Lafitte, H. Manche1, and M. Gastineau, “Gravity tests with INPOP planetary ephemerides,” in Relativity in Fundamental Astronomy, S. A. Klioner, P. K. Seidelman, and M. H. Soffel, Eds., Proceedings IAU Symposium no. 261, pp. 159–169, 2010.
- G. L. Page, J. F. Wallin, and D. S. Dixon, “How well do we know the orbits of the outer planets?” Astrophysical Journal, vol. 697, no. 2, pp. 1226–1241, 2009.
- L. Iorio, “Does the Neptunian system of satellites challenge a gravitational origin for the Pioneer anomaly,” Monthly Notices of the Royal Astronomical Society, vol. 405, no. 4, pp. 2615–2622, 2010.
- G. L. Page, “Exploring the weak limit of gravity at solar system scales,” Publications of the Astronomical Society of the Pacific, vol. 122, no. 888, pp. 259–260, 2010.
- L. Iorio, “Orbital effects of the time-dependent component of the Pioneer anomaly,” http://arxiv.org/abs/1107.3445.
- R.D. Reasenberg, I. I. Shapiro, P. E. MacNeil, et al., “Viking relativity experiment: verification of signal retardation by solar gravity,” Astrophysical Journal, vol. 234, pp. L219–L221, 1979.
- M. Mizony and M. Lachièze-Rey, “Cosmological effects in the local static frame,” Astronomy and Astrophysics, vol. 434, no. 1, pp. 45–52, 2005.
- M. Lachièze-Rey, “Cosmology in the solar system: the Pioneer effect is not cosmological,” Classical and Quantum Gravity, vol. 24, no. 10, article 016, pp. 2735–2742, 2007.
- V. T. Toth, “Independent analysis of the orbits of pioneer 10 and 11,” International Journal of Modern Physics D, vol. 18, no. 5, pp. 717–741, 2009.