Preprint
Article

This version is not peer-reviewed.

Massive Photon as Solution to Anomalous Acceleration of Deep-Space Probes and Some Other Cosmological Anomalies

Submitted:

24 September 2025

Posted:

26 September 2025

You are already at the latest version

Abstract
All deep-space probes seem to experience small anomalous acceleration. Although original thermal modelling of Pioneer spacecraft effectively ruled it out, newer modelling suggests that the anomalies do arise from the anisotropic thermal radiation, mainly associated with radioisotope thermoelectric generators. However, although the results are convincing, alternative solutions are not ruled out. Here, it is proposed that the thermal anisotropy is generally only a part of the solution and may not be a dominating part in all cases. The other contribution is hypothesized to come from a non-zero photon mass and its gravitational coupling. The exact solution for the Pioneer anomalies is presented, in remarkable agreement with measurements, and with the original thermal modelling. In addition, the solution provides an alternative explanation for the Schwarzschild radius of the observable universe and could potentially explain different cosmological anomalies, not limited to the anomalous acceleration of deep-space probes.
Keywords: 
;  ;  ;  

1. Introduction

The analysis of motion of deep-space probes shows anomalous constant weak long-range acceleration [1] on the order of -10 , 10 m/s , 2 (deceleration relative to the Sun). A solution has been proposed in the form of anisotropic thermal radiation [2] and it has even been claimed eventually that that this could completely resolve the Pioneer [3] anomaly [4]. However, although convincing, the result is not absolutely definitive. A large uncertainty exists in RTG surface degradation, for example. Although explained within uncertainties, thermal estimate at about 80% of the Doppler-inferred acceleration is not a perfect match. The apparent onset of the anomaly [5] only once the giant planets have been reached [6] (unexplainable by thermal anisotropy) is assumed to be a result of mismodelling of the solar thermal contribution, but this is not certain either. The observed diurnal/annual oscillation [7] in the anomalous acceleration (unexplainable by thermal anisotropy) is also assumed to be the artefact of mismodelling (including that of the spacecraft orientation). The best fit to the original data is a temporally constant anomaly [8], while the thermal anisotropy implies temporal decay of the anomaly. Although the decay cannot be ruled out - especially with the extended Pioneer data sets [9], the obtained marginal signs of decay likely are a result of mismodelled [10] solar radiation pressure [11].
The newer thermal modelling is significantly different from the original modelling (which effectively ruled out thermal anisotropy [7]), and others have criticized [12] it [13] (e.g., for the extensive parametrisation, potential bias, reliance upon an inferior quality extension of the highest quality data set to overturn the consensus view of a [mean] constant anomaly, the fact that only the "acceleration" value is considered, etc.).
Considering the involved assumptions, uncertainties, criticism, and the potential scientific value of a different explanation, alternative solutions [14] should not be easily dismissed. Here, an exact solution in the form of a massive photon is proposed, a solution that could, in addition, explain some other cosmological anomalies.

2. The Hypothesis

Assume that the photon is propagating as an expanding spherical wave, with mass concentrated at the expanding sphere surface (spherical shell). As such, the photon is affected by the gravitational potential of the enclosed mass (with the assumption of homogeneous/isotropic distribution, per the shell theorem, the outer mass would have no net effect) - with the gravitational coupling proportional to photon mass. Therefore, its linear momentum is converting to angular momentum, and, upon reaching its range, the photon’s linear (radial) velocity becomes equal to zero. Effectively, thus, the photon is experiencing relatively constant negative radial acceleration (or, increasing curvature in the context of General Relativity). This acceleration is the equivalent of inverse free fall, and on a return trip (assuming the photon is not absorbed and falls back upon reaching the range) it would be equivalent to free fall. The negative radial acceleration of photons (or, more precisely, the blueshift due to increasing curvature) will then be misinterpreted as a decrease in velocity (Doppler shift) of the moving source of emission if the mass of the photon is assumed to be zero. But in reality, there is no change in velocity of the source, the frequency shift occurs because the photon paths are curved and curvature changes for each new photon emitted. Note that the reference frame of the observer of an photon with the range and radius on the order of the radius of the observable universe would effectively be near an event horizon, as the enclosed mass-energy/density forms a black hole for the photon (although the photon is localized with the absorption).
Thus, the source of anomalous acceleration detected in deep-space probes is here proposed to be the massive photon. With such photon, the onset (proper term here would be "inflation") of the anomaly at particular distance can be explained with mass oscillation, settling of photons in higher mass eigenstates.
The transition to higher eigenstates should be correlated with the transition in coupling - e.g., from the density/pressure of the Solar System to that of the galaxy, or the cosmological vacuum. Interestingly, the anomaly onset seems to occur as the spacecraft themselves decouple from the Solar System [15] (pass from an elliptic orbit into an hyperbolic escape trajectory, following a gravity-assist Jupiter/Saturn flyby). In fact, since the hypothesis involves coupling of photons to gravity, the anomalies could occur with changes in gravitational coupling in general. Note that anomalous kinetic energy has been detected in several spacecraft even during gravity-assist Earth flyby [15].
The detected annual periodicity could then be explained by the motion of Earth about the Sun, as this is the place of photon absorption. Since the distance to the spacecraft is oscillating annually due to Earth’s motion, the maximal radius the photon can have before absorption is oscillating as well - hence the annual oscillation of the anomaly (this correlation, however, is not further analysed in this paper).
Note that the solution provided here also goes in favour of physical interpretation of the wave-function - the collapse of the wave-function is the collapse (localization) of the mass forming the propagating spherical waveform to the point of photon absorption.
It can be shown easily that, for a spatially flat expanding universe (where the observable radius is a Hubble radius), the observable radius is equal to the Schwarzschild radius of the total mass-energy inside the observable universe, which seems to be the case with our universe [16]. The solution presented here provides an alternative interpretation for the Schwarzschild radius, one that doesn’t require expanding universe (note, however, that the two interpretations are not mutually exclusive). If the enclosed mass grows exponentially with the expanding photon radius (dM/dr ≠ const.) - which is the case for a spherical photon and a large scale homogeneous/isotropic universe, at some point of the photon expansion the enclosed mass will effectively become a black hole for the photon, and as the photon radius becomes equal to the event horizon further expansion becomes impossible - this is the point where its radial velocity is completely converted to angular velocity.
Of course, on large distances (scales), the blueshift will be cancelled by the redshift due to universe’s expansion. The blueshift may also be negligible for the lowest mass eigenstates and potentially may also be limited with the partial localization of photons, constraining their radii - which may be more likely for higher frequencies.

3. The Solution and Discussion

In established quantum physics photon is assumed to have zero rest mass, however, this doesn’t have to be the case in reality, and there are reasons to believe it is not. Localized photon mass in experiments is on the order of 10 , 50 to 10 , 54 kg [17]. This is usually interpreted as the upper limit on its rest mass. However, it can also be interpreted as the actual photon mass associated with particular photon scale and correlated density/pressure. If the photon has mass the electro-magnetic potential is a Yukawa potential and the photon also has a range, equal to the [reduced] Compton wavelength, i.e.,:
r = m p c
where is the reduced Planck constant, c is the standard vacuum speed of light, and m , p is photon mass. This range should be, obviously, equal or larger than the observable universe for cosmological photons, which translates to the upper limit on the order of ∼10 , 69 kg. Others have calculated this mass, with the assumption of dS vacuum and a Ricci scalar of 4 Λ (where Λ is a positive cosmological constant), to be ≈2 × 10 , 69 kg [18]. Using matter density and pressure of the Solar System (Sun magnetosphere) in the Ricci scalar instead, and a zero cosmological constant, one obtains photon mass of ≈2 × 10 , 72 kg [18]. The author has obtained similar values in a different approach [19], where 3 mass eigenstates of the photon are hypothesized as well. In example, the mass of 6.335 × 10 , 69 kg was obtained as the cosmological photon tau equivalent eigenstate, and 1.822 × 10 , 72 kg as the lowest mass eigenstate. Using the higher mass, one obtains the radial acceleration:
a = 1 2 c 2 r = 1 2 c 2 c M γ τ = 8.093 × 10 10 m s 2
c = standard vacuum speed of light = 2.99792458 × 10 , 8 m/s
= reduced Planck constant = 1.054573 × 10 , 34 Js
M , γ τ = tau photon mass = 6.335 × 10 , 69 kg
which is in remarkable agreement with the observed Pioneer 10 anomaly of -8.09±0.20 × 10 , 10 m/s , 2 [1]. Similarly, if one assumes a superposition (sum) of tau and muon photon eigenstates, one obtains:
a = 1 2 c 2 c M γ τ + M γ μ = 8.574 × 10 10 m s 2
M , γ τ = tau photon mass = 6.335 × 10 , 69 kg
M , γ μ = muon photon mass = 3.767 × 10 , 70 kg
which is in remarkable agreement with the observed Pioneer 11 anomaly of -8.56±0.15 × 10 , 10 m/s , 2 [1]. The obtained results are also in agreement with the anomalous Galileo acceleration of -8±3 × 10 , 10 m/s , 2 [1].
Note that the ratio of mass between tau photon and muon photon is the same as the ratio of mass between tau electron and muon electron, as the author hypothesises that muon/tau eigenstates are not limited to electrons, rather represent a more general oscillation of particles. Note also that, assuming photons here generally couple to mass/gravity in pairs, an additional eigenstate should be added to the first equation. However, if this is the lowest mass eigenstate, the result wouldn’t change significantly. For both photons in lowest mass eigenstates, the acceleration becomes relatively negligible (on the order of 10 , 13 m/s , 2 ). The onset of the anomaly at particular distance is then explained as transition to higher mass eigenstates (which, however, may be more correlated with the decoupling from the Solar System, rather than with the distance itself).
However, even though the results above suggest the thermal anisotropy in Pioneer spacecraft is likely to be negligible (in agreement with the original modelling, where the value of 0.55±0.55 × 10 , 10 m/s , 2 has been obtained for the thermal anisotropy, with reasons given to consider this an upper bound [7]), higher anisotropy of thermal dissipation in the probes generally cannot be ruled out. In fact, thermal anisotropy could dominate over the proposed effect, especially for smaller photon mass eigenstates.
Coupling thermal modelling with the above, one could then obtain solutions in agreement with anomalies of other probes (e.g., detected Ulysses anomaly of 12±3 × 10 , 10 m/s , 2 [1], New Horizons anomaly of 13.2±0.6 × 10 , 10 m/s , 2 [20]).
In any case, at smaller distances gravitational coupling prevails, however, in an expanding universe, at larger distances the deceleration (blueshift) will be counteracted by the acceleration (redshift) due to the expansion. Therefore, if photons are generally propagating by the proposed mechanism, expansion of the universe may be underestimated, however, blueshift can also be limited by smaller mass eigenstates and possibly by partial localization of photons during propagation.
Note that the proposed photon nature and mechanism of propagation can potentially explain some other anomalies as well. By the hypothesis, a photon emitted from an celestial object may be reflected back towards the original point of emission upon reaching the range. Since celestial objects are generally in motion, an observer receiving both direct light from the object and reflected light will observe two images of the same object at different points in time, which could then be interpreted as two different but highly correlated objects even if they appear far away from each other. This could explain, for example, the alignment of many quasar polarization vectors over extremely large regions of the sky - billions of light-years apart, even though the quasars are not gravitationally bound [21]. Another potential example is the Huge Large Quasar Group [22]. Here, part of the group may be formed by reflections, so the actually physical group is smaller. However, the examination of the plausibility of the proposed explanation for these particular cases is beyond the scope of this paper.

3.1. Is the Universe Contracting?

It has been hypothesized that the observed effect depends on the enclosed mass. If the enclosed mass is increasing with each photon emitted, a blueshift is produced. With constant energy density of space, the enclosed mass is increasing with each photon emitted if the source and absorber are moving apart. Thus, a blueshift is expected (although it will be counteracted with the redshift produced with motion of objects, or universe’s expansion).
Conversely, if the universe is contracting, without increasing energy density, the enclosed mass should be decreasing with each photon emitted. Thus, a redshift would be produced instead of blueshift, even for otherwise relatively stationary sources. If, however, the contraction of the universe dominates on larger scales (like dark energy), on smaller scales, for two objects moving apart, a blueshift would be produced instead (the effect depends solely on eclosed mass and photon mass/range, so as long as the enclosed mass is increasing with each new photon a blueshift is produced).
Note that, since this is effectively a gravitational frequency shift, it includes other relativistic effects, such as time dilation. If one now assumes there is no partial localization of photons during propagation (affecting enclosed mass), the observed increasing redshift with distance in the observable universe may have been misinterpreted. Instead of expanding, the universe may be actually contracting.
Note that the lowest mass eigenstate of the photon, per the hypothesis, is equal to 1.822 × 10 , 72 kg. Per Equation (2), this produces the acceleration of -2.328 × 10 , 13 m/s , 2 , or double that value for the sum of two eigenstates. Expressing this in a value analogous to the Hubble constant (which here, being negative, may be interpreted as contraction), one obtains:
A 0 = 2.328 × 10 13 m / s 2 2 × 3.086 × 10 22 m 2.328 × 10 13 m / s 2 1 M p c = 119.86 ( k m / s ) / M p c
And, for the sum of two eigenstates:
B 0 = 2 × 2.328 × 10 13 m / s 2 2 × 3.086 × 10 22 m 2 × 2.328 × 10 13 m / s 2 1 M p c = 2 2 A 0 = 169.50 ( k m / s ) / M p c
It may be a coincidence, but these values are suspiciously close to the inferred values of the Hubble constant (H , 0 ), suggesting that the dark energy may be coupled to a particle of similar mass to photon. Assuming now that the expansion of the universe is actually equal to one of these values (albeit positively signed), one can obtain the Hubble constant as the superposition of expansion and the effective contraction (a product of the proposed nature of the photon [propagation]). Simple addition gives, for the expansion equal to B , 0 in value:
H 0 = B 0 + A 0 = 49.64 ( k m / s ) / M p c
Associating the sum of 4 particles with dark energy and 2 with photon coupling, one obtains:
H 0 = 4 4 2 2 A 0 = 70.21 ( k m / s ) / M p c
A value that is about the average of the measured values, and is in agreement with many measurements/models [23]. Based on the above, one possible solution to the Hubble tension may be in the changing, or the running, gravitational coupling.

3.2. Derivation of Photon Mass

In the theory of Complete Relativity [19] the author postulates discrete vertical energy levels. These energy levels are somewhat analogous to the energy levels of particles bound to atomic nuclei. One major difference is that the progression is logarithmic so the difference between energies is in the order of magnitude (vertical), rather than in the value (horizontal). Another salient difference is that the equivalence between the particle on one energy level and the other is very relative. And this is due to the associated running coupling and energy transformation between scales (vertical energy levels). On one energy level, for example, the particle may be electrically charged with its gravitational mass being negligible, on the other energy level gravity may dominate (one type of energy is always exchanged for the other with the transition between levels).
The vertical energy levels can explain many unexplained features of the observable universe. Here, however, due to limited scope, only the equation for their calculation will be provided (additional details, including the evidence for their existence can be found in other [19] works [24] of the author):
log ( E n ) = log ( E n 1 ) + ( n 1 ) n = log ( 10 ( n 1 ) n E n 1 )
n = vertical energy level (integer)
Assuming now that standard particles occupy level n = 7, electron’s [equivalent] mass on level n = 8 is:
M N = M e 8 = 10 log ( M e ) + 7 × 8 = 0.910938356 × 10 26 k g
M , e = standard electron mass = 9.10938356 × 10 , 31 kg
Note that this is on the order of Neptunian planets, and very close to the mass of Neptune (which, author argues is not a coincidence [24]). Similarly, using masses of atoms and other standard particles, masses on the order of stars and other celestial objects can be obtained.
Calculating now the electron mass on level n = 6, one obtains:
M e 6 = 10 log ( M e ) 6 × 7 = 9.10938356 × 10 73 k g
How to interpret this particle? Well, there are not many known candidates for such mass, but it is within the orders of mass expected for cosmological photons, and since photons mirror the source of emission (charge), it shouldn’t be surprising that the electron equivalent on level n = 6 is a component of photons. The simplest composite structure of a photon involves 2 fermions. Assuming now that one component is electron on n = 6, and the other is a positron on the same level, and assuming negligible binding energy compared to component masses, one obtains a photon mass of:
M γ = 2 × M e 6 = 1.821876712 × 10 72 k g
Now, since the equation for vertical energy levels is the same for particles of different mass, obviously, the photon should have 3 mass eigenstates and the ratios between these eigenstates should be equal to the ratios between standard tau, muon and electron (e) eigenstates. Therefore, assuming the obtained photon mass is the mass of the e eigenstate, masses of other eigenstates can be easily obtained:
M γ μ = M μ M e M γ = 3.767 × 10 70 k g
M γ τ = M τ M e M γ = 6.335 × 10 69 k g
M , μ = standard muon mass = 105.6584 MeV/c , 2 [25]
M , τ = standard tau mass = 1776.86±0.12 MeV/c , 2 [25]
M , e = standard electron mass = 0.510999 MeV/c , 2 [25]

4. Conclusion

It has been shown that anomalous acceleration of Pioneer probes can be fully explained with massive photons propagating by the proposed mechanism, going in favour of the original thermal modelling of Pioneer spacecraft. Although newer thermal modelling has produced convincing results, it remains questionable whether it is the proper solution of the anomalies. The solution provided here can explain multiple aspects of the anomaly, not limited to the value of "acceleration". It can also explain anomalies in other deep-space probes, but may in some cases represent only a partial solution (the other being thermal recoil). The solution provides an alternative interpretation for the Schwarzschild radius of the observable universe - without excluding the established one. Additionally, it could potentially explain unexpected cosmological correlation and some other cosmological anomalies. In any case, experimental verification of the proposed solution is desirable, as any deviation from conventional physics could have large implications for cosmology, astrophysics and particle physics. A deep-space probe, for example, designed in such way to minimize or rule out anisotropic thermal dissipation could provide an unambiguous confirmation or refutation of alternative solutions, so it would be useful even in the case the solution proposed here doesn’t hold up to scrutiny.

References

  1. Indication, from Pioneer 10/11, Galileo, and Ulysses Data, of an Apparent Anomalous, Weak, Long-Range Acceleration (1998), J. D. Anderson et al., Phys. Rev. Lett. 81, 2858. [CrossRef]
  2. The constancy of the Pioneer anomalous acceleration (2007), Ø. Olsen, A&A 463, 393-397. [CrossRef]
  3. High precision thermal modeling of complex systems with application to the flyby and Pioneer anomaly (2011), B. Rievers and C. Lämmerzahl, Annalen der Physik 523(6), 439-449. [CrossRef]
  4. Support for the Thermal Origin of the Pioneer Anomaly (2012), S. G. Turyshev et al., Phys. Rev. Lett. 108(24), 241101. [CrossRef]
  5. Using early data to illuminate the Pioneer anomaly (2005), M. M. Nieto and J. D. Anderson, Class. Quantum Grav. 22(24), 5343. [CrossRef]
  6. New Horizons and the onset of the Pioneer anomaly (2008), M. M. Nieto, Physics Letters B 659(3), 483-485. [CrossRef]
  7. Study of the anomalous acceleration of Pioneer 10 and 11 (2002), J. D. Anderson et al., Phys. Rev. D 65(8), 082004. [CrossRef]
  8. Pioneer 10 Doppler data analysis: Disentangling periodic and secular anomalies (2009), A. Levy et al., Advances in Space Research 43(10), 1538-1544. [CrossRef]
  9. Support for Temporally Varying Behavior of the Pioneer Anomaly from the Extended Pioneer 10 and 11 Doppler Data Sets (2011), S. G. Turyshev et al., Phys. Rev. Lett. 107(8), 081103. [CrossRef]
  10. Mystery Pioneer Anomaly Is Real But Still a Mystery (2011), R. A. Kerr, Science 333(6047), 1208. [CrossRef]
  11. Chameleon effect and the Pioneer anomaly (2012), J. D. Anderson and J. R. Morris, Phys. Rev. D 85(8), 084017. [CrossRef]
  12. The Pioneer Anomaly: an inconvenient reality or NASA’s 12 year misconception? (2013), P. G. ten Boom, preprint [arXiv]. [CrossRef]
  13. Comments on “The Pioneer Anomaly: an inconvenient reality or NASA’s 12 year misconception” (2013), J. C. Hodge, last access 2025.09.11https://intellectualarchive.com/?link=item&id=1133.
  14. Scalar potential model of the Pioneer Anomaly (2007), J. C. Hodge, preprint [arXiv]. [CrossRef]
  15. The Pioneer Anomaly (2010), S. G. Turyshev and V. T. Toth, Living Rev Relativ. 13(1), 4. [CrossRef]
  16. All objects and some questions (2023), C. H. Lineweaver and V. M. Patel, American Journal of Physics 91(10), 819-825. [CrossRef]
  17. Rest mass of photon on the surface of matter (2020), M. Goray et al., Results Phys. 16, 102866. [CrossRef]
  18. Photon mass as a probe to extra dimensions (2016), G. Alencar et al., Physics Letters B 759, 138-140. [CrossRef]
  19. Complete Relativity: Nature of Observables (2025), M. Ljubičić, [draft]. [CrossRef]
  20. When leaving the Solar system: Dark matter makes a difference (2022), E. Belbruno and J. Green, Mon. Not. R. Astron. Soc. 510(4), 5154-5163. [CrossRef]
  21. Alignment of quasar polarizations with large-scale structure (2014), D. Hutsemékers et al., Astronomy and Astrophysics 572, A18. [CrossRef]
  22. A structure in the early Universe at z ∼ 1.3 that exceeds the homogeneity scale of the R-W concordance cosmology (2013), R. G. Clowes et al., Mon. Not. R. Astron. Soc. 429(4), 2910–2916. [CrossRef]
  23. In the realm of the Hubble tension - a review of solutions (2021), E. Di Valentino et al., Class. Quantum Grav. 38(15), 153001. [CrossRef]
  24. The Solar System: Nature and Mechanics (2025), M. Ljubičić, [draft]. [CrossRef]
  25. Review of particle physics (2022), R. L. Workman et al (Particle Data Group), Prog. Theor. Exp. Phys. 2022(8), 28-1245. [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated