1. Introduction
From a practical standpoint, general relativity (GR) has been validated by numerous observations spanning various fields. However, its current status presents challenges. While supported by substantial observational evidence, GR faces significant hurdles. Its verifications in cosmology and astrophysics are under scrutiny, primarily because it relies on unproven concepts like dark matter and energy to explain phenomena on a large scale, such as galaxies and the universe. Often, these unconfirmed elements are employed while simultaneously overlooking a crucial aspect of GR: retardation, which contradicts Newtonian principles of action at a distance. This discrepancy may be connected to the problems we have understanding of gravitational interactions on cosmic scales.
The mystery surrounding dark matter has long been a topic of discussion within the astronomical community, dating back to the 1930s, and possibly even earlier in the 1920s when it was referred to as the "question of missing mass." Over time, this enigma has only grown more prominent, particularly as the need for dark matter (and the tendency to overlook retardation) has increased on larger scales under examination. Despite extensive and costly efforts, including a forty-year search conducted underground and using accelerators, dark matter’s existence remains unproven. Recent years have only added to the challenge, as the Large Hadron Collider’s failure to detect any super-symmetric particles—a favored form of dark matter among astro-particle physicists—poses further complications. These particles are not only crucial for understanding dark matter but also play a pivotal role in string theory, which is anticipated to provide insights into the quantization of gravity.
As far back as 1933, Zwicky observed anomalies in the velocities of galaxies inside the Coma Cluster that exceeded predictions based on Newtonian theory. He calculated [
1,
2] that the required amount of matter to explain these velocities could be about 400 times greater than that of visible matter, although later adjustments mitigated this discrepancy to some extent. Had Zwicky utilized the concept of retarded gravity in his calculations, this issue might have been resolved without significant complication [
3]. In 1959, Volders [
4] noted similar discrepancies on a smaller scale within the outer parts of the nearby spiral galaxy M33, where velocities did not follow the expected
pattern. This observation was corroborated in subsequent years by Rubin and Ford [
5,
6], who demonstrated that velocities at the outer rim of many spiral galaxies either plateaued or continued to increase, each at a different velocity. Previous studies have indicated that such velocity patterns can be directly inferred from general relativity (GR) if retardation effects are considered [
7,
8,
9,
10,
11,
12,
13,
14]. The mechanism underlying retardation is intricately linked to the dynamics of the matter density within galaxies, specifically to the second derivative of density. Changes in density may result from various factors, including gas depletion in surrounding intergalactic gas [
12] or dynamic processes such as star formation and supernova explosions [
10,
11]. These processes can be characterized by three different typical length scales: the density gradient length, velocity field gradient length, and dynamical length scale. The importance of retardation is determined by the shortest among these length scales [
15].
The famous Tully-Fisher relation [
16], which links the baryonic mass of a galaxy to the fourth power of its rotational velocity at the outer rim, can also be derived from the principles of retarded gravity [
17]. It was shown that the effects of retarded gravity extend beyond just slowly moving particles and also apply to photons. While there may be some differences in the mathematical analysis for each case, it is ultimately concluded that the observed "dark mass" inferred from galactic rotation curves must be identical in both the lensing and rotation curves scenarios [
15,
18].
While the prevailing notion of dark matter remains prominent, the current circumstances warrant consideration of the case that this prevailing paradigm may need to be reevaluated. Several challenges cast doubt on this common idea:
Firstly, in order to align with observed phenomena and structure formation simulations, a set of properties has been assigned to dark matter (DM) [
2]. However, despite being over 50 years since its inception, dark matter has yet to be directly observed, nor have any known particles been identified that match its purported properties.
Secondly, simulations involving dark matter often encounter what is known as the core-cusp problem. The Navarro-Frenk-White (NFW) [
19] profile, derived from Cold Dark Matter (CDM) simulations and commonly used to fit rotation curves, is a prime example. However, this profile faces challenges, particularly when applied to Low Surface Brightness galaxies (LSBs). Derivations made using the NFW profile regarding rotational velocities frequently diverge from actual observations, leading to discrepancies. Specifically, while the NFW profile anticipates a "cuspy" internal region for a dark halo (where density changes rapidly), observations tend to favor a "core-like" behavior (where density remains approximately constant). Efforts to address this issue have often relied on specific and somewhat contrived adjustments, raising doubts about whether these solutions were devised primarily to maintain the current paradigm.
Thirdly, Sancisi’s Law [
20] presents a significant and broadly applicable observation. It suggests that changes in the luminosity profile of a galaxy correspond to changes in its rotation curve, and vice versa. This phenomenon applies to various types of dark halos. However, from a dark matter perspective, this relationship is unexpected: the dark halo is typically assumed to be much more massive than the baryonic matter. Consequently, fluctuations in the distribution of baryonic matter should not significantly affect the velocity distribution, contrary to what is observed. This discrepancy is particularly pronounced in LSBs, where the dark halo is believed to dominate at every radius, yet the velocity distribution exhibits fluctuations corresponding to each "baryonic bump." This suggests that, somehow, the overall velocity distribution is influenced by small fluctuations in baryonic matter.
Hence, the current retarded gravity proposition offers a unique perspective. Unlike alternative theories that propose modifications to general relativity, such as Milgrom’s Modified Newtonian Dynamics (MOND) [
21], Mannheim’s Conformal Gravity [
22,
23], or Moffat’s Modified Gravity (MOG) [
24], our approach does not seek to alter the fundamental framework of general relativity. Instead, we adhere strictly to the principle of Occam’s razor, as advocated by both Newton and Einstein. Our objective is to replace the need for dark matter with phenomena inherent within standard General Relativity itself (retardation). It’s worth noting that recent research has shed light on the relationship between retardation and MOND [
25], demonstrating how criteria for low acceleration MOND can be deduced from retardation theory, and how the MOND interpolation function can approximate retarded gravity effectively.
It’s essential to highlight that significant retardation effects are not contingent upon high velocities of matter within galaxies, although higher velocities may enhance these effects. In reality, the majority of galactic constituents, such as stars and gas, move relatively slowly compared to the speed of light. This is indicated by the ratio of the velocity v to the speed of light in vacuum c, denoted as , which is much smaller than 1. For instance, typical velocities within galaxies are around 100 km/s, resulting in a ratio of or smaller. This will be discussed in more details in the sections that follow.
In contrast to the solar system, where retardation effects are considered negligible [
26,
27], observations of galaxies’ velocity curves suggest that these effects become significant beyond a certain distance [
10,
11,
12]. Recent research [
28] has expanded the empirical basis for the theory of retarded gravity. Building upon previous studies that analyzed eleven galaxies [
11], the latest research extends its scope to a larger sample of 143 galaxies sourced from the SPARC Galaxy collection. These galaxies vary in type, size, and luminosity. The analysis indicates that in most cases, an excellent fit to the observed data is achieved without the need to postulate dark matter or modify general relativity (see
Figure 1 for some examples).
As we shall show below this is not an accident but is rather dictated by general relativity.
2. A Discrete Model
Retarded gravity emerges from the weak field approximation of general relativity [
12]. In this framework, the metric perturbation
can be expressed in terms of a retarded potential
as follows [
12,
15]:
where
G represents the gravitational constant,
denotes the location where the potential is calculated,
signifies the whereabout of the mass element producing the potential,
, and
represents the mass distribution. The characteristic duration
for galaxies might span a few tens of thousands of years, but small relative to the timescale over which galactic density significantly changes. Therefore, we can express the density using a Taylor series expansion:
By substituting
rhotay into Equation photonequastartpoint and retaining the first three terms, we can derive:
The initial term in the series is referred to as the Newtonian potential:
The second term doesn’t influence the force acting on subluminal particles since its gradient is zero. As for the third term, it serves as a lower-order ammendment to the Newtonian potential.
The geodesic equation governing the motion of a "slow" test particle in the given space-time metric can be evaluated by utilizing the force per unit mass [
12]:
The total acceleration is thus:
Now, let’s examine a point particle with a mass
positioned at
. The particle will possess a mass density given by:
Here,
represents a three-dimensional Dirac delta distribution. This particle induces a Newtonian potential given by:
and a retardation potential in the following form:
A test particle at the vicinity of particle
j will be accelerated as follows:
in which the reader should not confuse the acceleration of the point particle
j denoted
and the acceleration caused by particle on a test particle located at point
denoted
.
3. The Uncertainty Relation of Retarded Gravity
For non relativistic matter:
hence we may write approximately:
Thus in order for retarded gravity to have a significant effect
must be of the order of unity (
) or larger, this leads to the retarded gravity "uncertainty" type relation:
Consider a point mass located on a circle which serves as border of the M33 galaxy, then we may ask what will be the amount of acceleration suffered by the point mass that will cause a retardation effect on a test particle located across the diameter of the galaxy (which is the furthest point on the imaginary circle from the point mass). Now the radius of the galaxy M33 is
light years
meters. Hence we will need an acceleration of about:
to observe the effect of retarded gravity. This does not seem to be such a huge acceleration and many point masses (atoms, molecules etc.) in the galaxy may have accelerations that need to be considered in the total galactic balance of gravitational forces. From this point of view we may partition the galactic point masses in to two classes: Newtonian gravity particles and retarded (+Newtonian) gravity particles the difference depends on how big is the Newtonian radius (which depends on the particle’s acceleration):
This reality is depicted in
Figure 2. Of course each test particle in not affected by just one massive particle but by all
massive particles this leads to the equation:
in which RG means particles that have a retarded influence in point
, that is particles in which point
is outside their Newtonian radius. We point out that the stellar component of disk galaxies is not responsible for the retardation effects. To see this look at the velocity & acceleration curves of the M33 galaxy depicted in
Figure 3:
Assuming stars (and gas) moving in circles we obtain accelerations which do not exceed
m/s
2 yielding a Newtonian Radius not smaller
kpc, hence the effects of those stars and rotating gas is completely Newtonian within the galaxy. Thus to obtain retardation corrections one must look at other processes taking place in galaxy like stellar winds, supernovae explosions [
11] and the depletion of gas in the intergalactic medium [
12] leading to the deceleration of the rate of mass accretion by the galaxy and to an attractive gravitational retardation effect.