Lensing Effects in Retarded Gravity

Galaxies are gigantic physical systems having a typical size of many tens of thousands of light years. Thus any change at the center of the galaxy will affect the rim only tens of millennia later. Those retardation effects seems to be ignored in present day modelling used to calculate rotational velocities of matter in the outskirts of the galaxy and the surrounding gas. The significant discrepancies between the velocities predicted by Newtonian theory and observed velocities are usually handled by either assuming an unobservable type of matter denoted"dark matter"or by modifying the laws of gravity (MOND as an example). Here we will show that considering general relativistic effects without neglecting retardation one can explain the apparent excess matter leading to gravitational lensing in both galaxies and galaxy clusters.


Introduction
Einstein's general relativity (GR) is known to be invariant under general coordinate modifications.This group of general transformations has a Lorentz subgroup, which is valid even in the weak field approximation.This is seen through the field equations containing the d'Alembert (wave) operator, which can be solved using a retarded potential solution.
It is known that GR is verified by many types of observations.However, currently, Newton-Einstein gravitational theory is at a crossroads.It has much in its favor observationally, and it has some very disquieting challenges.The successes that it has achieved in both astrophysical and cosmological scales have to be considered in light of the fact that GR needs to appeal to two unconfirmed ingredients, dark matter and energy, to achieve these successes.Dark matter has not only been with us since the 1920s (when it was initially known as the missing mass problem), but it has also become severe as more and more of it had to be introduced on larger distance scales as new data have become available.Here we will be particularly concerned in the excess dark matter needed to justify observed gravitational lensing.Moreover, 40-year-underground and accelerator searches and experiments have failed to establish its existence.The dark matter situation has become even more disturbing in recent years as the Large Hadron Collider was unable to find any super symmetric particle candidates, the community's preferred form of dark matter.
While things may still take turn in favor of the dark matter hypothesis, the current situation is serious enough to consider the possibility that the popular paradigm might need to be amended in some way if not replaced altogether.The present paper seeks such a modification.Unlike other ideas such as Milgrom's MOND [1], Mannheim's Conformal Gravity [2][3][4], Moffat's MOG [5] or f (R) theories and scalar-tensor gravity [6], the present approach is, the minimalist one adhering to the razor of Occam.It suggests to replace dark matter by effects within standard GR.
Fritz Zwicky noticed in 1933 that the velocities of Galaxies within the Comma Cluster are much higher than those predicted by the virial calculation that assumes Newtonian theory [7].He calculated that the amount of matter required to account for the velocities could be 400 times greater with respect to that of visible matter, which led to suggesting dark matter throughout the cluster.In 1959, Volders, pointed out that stars in the outer rims of the nearby galaxy M33 do not move according to Newtonian theory [8].The virial theorem coupled with Newtonian Gravity implies that MG/r ∼ Mv 2 , thus the expected rotation curve should at some point decrease as 1/ √ r.During the seventies Rubin and Ford [9,10] showed that, for a large number of spiral galaxies, this behavior can be considered generic: velocities at the rim of the galaxies do not decrease-but they attain a plateau at some unique velocity, which differs for every galaxy.We have shown that such velocity curves can be deduced from GR if retardation is not neglected.The derivation of the retardation force is described in previous publications [11][12][13][14][15][16][17], see for example figure 1. PhD student under my supervision, using [40]; the full line describes the rotation curve, which is the sum of the dotted line, describing the retardation contribution, and the dashed line, which is Newtonian.
Previous work assumed a test particle moving slowly with respect to the speed of light as is appropriated for the case of galactic rotation curves, this is not adequate when the test particle is a photon moving in the speed of light as in the case relevant to gravitational lensing.Here, a different mathematical approach is needed as described in the current paper.
A gravitational lens is some form of matter (for example a cluster of galaxies) between a distant source of light and the observer, that is bending the light as it travels towards the observer.This effect is denoted gravitational lensing, the amount of bending is one of the predictions of GR [18][19][20].It should be noted that Newtonian physics also predicts light bending, but only half of that predicted by GR [21].
Einstein made unpublished work on gravitational lensing as early as 1912 [22].In 1915 Einstein showed how GR explained the anomalous perihelion advance of the planet Mercury without any arbitrary parameters [23], in 1919 an expedition led by Eddington confirmed GR's prediction for the deflection of starlight by the Sun in the total solar eclipse of May 29, 1919,[24,25] making Einstein famous [23] instantly.The reader should recall that there was a special significance to a British scientist confirming the prediction of a German scientist after the bloody battles of world war I.
The fact that distortion of space time is proportional to the amount of mass that causes the distortion led to the use of gravitational lensing as a tool for proving the existence of dark matter.
Strong lensing is the distortion of background galaxies into arcs when their light passes through a gravitational lens.This has been observed around many clusters such as Abell 1689 [26].By measuring the observed geometry, the mass of the in-between cluster can be calculated.In the cases where this was done, the mass-to-light ratios corresponded to the dynamical dark matter of clusters [27].Here we will show that is not a coincidence, and retardation dictates that this should be so.Lensing can cause multiple copies of an image.By analyzing multiple image copies, astronomers have been able to map the distribution of dark matter around the MACS J0416.1-2403galaxy cluster [28,29].
Weak gravitational lensing is concerned with small distortions of galaxies, using statistical methods from huge galaxy surveys.By examining the shear deformation of the adjacent background galaxies, the distribution of dark matter can be calculated.The mass-to-light ratios correspond to dark matter densities predicted by other large-scale structure measurements [30].
We underline that dark matter does not bend light itself; mass (in this case the alleged mass of the dark matter) bends spacetime.Light follows the curvature of spacetime, resulting in the lensing effect [31,32].Here we will show that the above described effects may be attributed to retardation and baryonic matter, with no additional form of (dark) matter.
The structure of the paper is as follows: First we describe the main results of general relativity, this is followed by the weak field approximation.Next we derive the geodesic equations for a particle moving at a speed of light ("a photon").This is followed by a comparison of the current approach to that of Weinberg [35].The classical result of Einstein and Eddington are re-derived.Finally we discuss the effect of retardation on the lensing phenomena and show in what way can retardation cause a "dark matter" phenomena.

General Relativity
The general theory of relativity is based on two fundamental equations, the first being Einstein equations [33][34][35][36]: G µν stands for the Einstein tensor (see equation ( 7)), T µν indicates the stress-energy tensor (see equation ( 3)), G 6.67 10 −11 m 3 kg −1 s −2 is the universal gravitational constant and c 3 10 8 m s −1 indicates the velocity of light in the absence of matter (Greek letters are indices in the range 0 − 3).The second fundamental equation that GR is based on is the geodesic equation: x α (p) are the coordinates of the particle in spacetime, p is a typical parameter along the trajectory that for massive particles is chosen to be the length of the trajectory (p = s), u µ = dx µ dp is the µ-th component of the 4-velocity of a massive particle moving along the geodesic trajectory (increment of x per p) and Γ α µν is the affine connection (Einstein summation convention is assumed).The stress-energy tensor of matter is usually taken in the form: In the above, pr is the pressure and ρ is the mass density.We remind the reader that lowering and raising indices is done through the metric g µν and inverse metric g µν , such that u µ = g µν u ν .The same metric serves to calculate s: for a photon ds = 0.The affine connection is connected to the metric as follows: Using the affine connection we calculate the Riemann and Ricci tensors and the curvature scalar: which, in turn, serves to calculate the Einstein tensor: Hence, the given matter distribution determines the metric through Equation ( 1) and the metric determines the geodesic trajectories through Equation (2).Those equations are well known to be symmetric under smooth coordinate transformations (diffeomorphism).

Linear Approximation of GR -The Metric
Only in extreme cases of compact objects (black holes and neutron stars) and the primordial reality or the very early universe does one need not consider the solution of the full non-linear Einstein Equation [11].In typical cases of astronomical interest (the galactic case included) one can use a linear approximation to those equations around the flat Lorentz metric η µν such that (Private communication with the late Professor Donald Lynden-Bell): One then defines the quantity: hµν = h µν for non diagonal terms.For diagonal terms: The general coordinate transformation symmetry of Equation ( 8) has a subgroup of infinitesimal transformations which are manifested in the gauge freedom of h µν in the weak field approximation.It can be shown ([33] page 75, exercise 37, see also [34][35][36]) that one can choose a gauge such that the Einstein equations are: The d'Alembert operator 2 is clearly invariant under the Lorentz symmetry group (another subgroup of the general coordinate transformation symmetry described by Equation ( 8)), of which the Newtonian Laplace operator ∇ 2 is not, but this comes with the price that "action at a distance" solutions are forbidden and only retarded solutions are allowed.The T µν stress energy tensor should be calculated at the appropriate frame and thus for the massive body causing the lensing, matter is approximately at rest.Equation (12) can always be integrated to take the form [37]: For reasons why the symmetry between space and time is broken, see [38,39].The factor before the integral is small: 4G c 4 3.3 × 10 −44 in MKS units; hence, in the above calculation one can take T µν , which is zero order in h αβ .
In the zeroth order: in which we assume that massive body causing the lensing effect is composed of massive particles (those equations will not be correct for the photon which is affected by the lensing potential).
Assuming the reasonable assumption that the said massive body is composed of particles of non relativistic velocities: Let us now look at equation (3).We assume ρc 2 pr and, taking into account equation (15), we arrive at T 00 = ρc 2 , while other tensor components are significantly smaller.Thus, h00 is significantly larger than other components of hµν which are ignored from now on.One should notice that it is possible to deduce from the gauge condition in equation ( 12) the relative order of magnitude of the relative components of Thus, the zeroth derivative of h00 (which contains a 1 c as x 0 = ct) is the same order as the spatial derivative of h0a meaning that h0a is of order v c smaller than h00 .And the zeroth derivative of h0a (which appears in Equation ( 16)) is the same order as the spatial derivative of hab .Meaning that hab is of order v c with respect to h0a and of order ( v c ) 2 with respect to h00 .In the current approximation, the following results hold: (The underline aa signifies that the Einstein summation convention is not assumed).
We can summarize the above results in a concise formula: in which δ µν is Kronecker's delta.It will be useful to introduce the gravitational potential φ which is defined below and can be calculated using Equation ( 13): from the above definition and equation ( 18) it follows that:

Linear Approximation of GR -The Lensing Trajectory
Let us start calculating the lensing trajectory by writing u µ = dx µ dp in terms of the notation introduced in equation ( 13): it thus follows that: Taking into account equation ( 2) we obtain: multiplying by dp dt 2 and using the notation of equation ( 14) we obtain: However, according to equation ( 2): Inserting equation (28) into equation ( 27) we arrive at the form: Let us now calculate the affine connection in the linear approximation: Taking into account equation (30) and equation ( 21) we obtain: The affine connection has only first order terms in h αβ ; hence, to the first order Γ α µν v µ v ν appearing in the geodesic, v µ v ν is of the zeroth order.By definition v 0 = dx 0 dt = c.Also the null interval of the photon is given to zeroth order in h 00 (that is in the absence of a gravitational field) by: It follows that in the zeroth order approximation: and also in the same approximation: We now combine the above results and write: We now turn our attention to the second term in equation ( 29) inserting equation ( 21) into equation ( 36) we arrive at the result: Combining the results from equation ( 35) and equation ( 37) into equation ( 29) we arrive at the photon's equation of motion: In terms of the gravitational potential this can be rewritten using equation (23) as follows: This equation is almost identical to equation 9.2.6 of Weinberg [35].Notice, however, that Weinberg neglects the term 2 v c 2 ∂ t φ in his post Newtonian approximation.Although this term can be neglected or shown to be small in specific circumstances such as case that φ is static or slowly varying, its removal may lead to inconsistencies as we explain below.
Let us inspect the null interval equation ( 4): taking into account equation ( 9) and equation (21) this takes the form: Hence: Or: compare to Weinberg [35] equation 9.2.5.Since by equation (22) and equation ( 23) h 00 is a small negative number, it follows that v < c (notice that this result holds for a global coordinate system, in the local flat coordinate system the velocity will of course be exactly c).Now: To first order in h 00 : Taking into account equation (39): Thus to the same first order of h 00 we have: which is dependent on not neglecting the ∂ t h 00 term.Thus this term is absolutely necessary to maintain the identity of equation ( 44) and without it we are led to a contradiction.We conclude that for a time dependent gravitational potential the ∂ t h 00 term is required.
Let us consider a photon travelling at a straight line in the direction v0 in the absence of gravity, the velocity of this photon would be v 0 = c v0 .Now let us assume that the photon passes near a weak gravitational source such that equation ( 39) is valid.We thus decompose the photon velocity field to two components one which parallel and one which is perpendicular to its original direction: Thus it follows from equation ( 39) that: Thus to first order in h 00 : in which we remember that v = v 0 + O(h 00 ) = c v0 + O(h 00 ).This can be integrated as follows: in which k is a constant.Far away from the gravitational source lim , and we may write: Now: And according to equation ( 43): It now follows from equation ( 48): The lensing angle is defined (see figure 2) as the angle α is small, thus for lensing in the linear approximation only the perpendicular component v ⊥ is important.Now to the first order in h 00 we may write equation (39) as: in which we define the perpendicular gradient as: ∇ ⊥ ≡ ∇ − v0 ( v0 • ∇).Now taking into account equation (48) and equation ( 50) we arrive at the result: in which we have used equation (23) to give the perpendicular acceleration in terms of a gravitational potential.

Other approaches to the problem of lensing
Another approach to the problem of lensing than the one given above is to start from a Schwarzschild metric as is done by Weinberg [35].This metric describe a static spherically symmetric mass distribution and thus is less general than the approach taken in this paper.It does have one advantage, however, and this is the ability to take into account strong gravitational fields and not just weak ones.This advantage is irrelevant in most astronomical cases in which gravity is weak and must be only considered for light trajectories near compact objects (black holes & neutron stars).The Schwarzschild squared interval can be written as: In which r , θ, φ are spherical coordinates and the point massive body is located at r = 0.The Schwarzschild radius is defined as: in which M is the mass of the point particle.Comparing the g 00 component of equation ( 59) and equation ( 41) it follows that we can identify: provided r s r 1 in accordance with equation ( 9), hence the two approaches should coincide for: Going back to equation (41) we have: it tempting to write d x 2 = dr 2 + r 2 dθ 2 + sin 2 θdφ 2 as is usually done for spherical coordinates.But notice that r is not strictly a radial coordinate which is defined as the circumference, divided by 2π, of a sphere centered around the massive body.In fact from equation (63) it is clear that the appropriate radial coordinate is: which is a small correction to r.Now calculating the differential dr it follows that: According to equation (61): where the last equality is correct to first order.Alternatively one can use equation ( 22) to calculate the gravitational potential for a static point mass to obtain: and then plug this into equation ( 23) to obtain again: It now follows that: Plugging equation (69) into equation (65) leads to: Hence to first order in h 00 : Using the results equation (64) and equation (71) the interval given in equation ( 63) can be rewritten as: As to first order in h 00 : and taking into account equation (61) we obtain: Thus to first order our metric is identical to Schwarzschild's for the case of a static point particle.This makes our analysis superior as it addresses the case of a general density distribution and does not ignore the possibility of time dependence which is crucial for retardation effects to take place.

Lensing in the static case
Newtonian theory dictates that a body with any mass (since inertial mass and gravitational mass are equal) moving under the influence of gravity alone must follow a trajectory which is dictated by the equation: The Newtonian potential φ N resembles φ given by equation ( 22) but neglects the retardation effect such that: In [15] we have shown that the geodesic equations reduce for slow moving test particles to: this will coincide with equation (75) for a static density distribution or for a slowly changing mass distribution as is well known.For light rays we obtained equation (58), this also resembles equation (75) but carries some major differences even for a static mass distribution.First there is a factor 2 multiplying the potential which is missing in equation ( 75).Second this equation only describe motion perpendicular to the original direction of the light ray and not the propagation in the direction of the light ray itself.Suppose the gravitational field is orthogonal to the light ray at least when the light ray is moving in proximity to the gravitating body, that is when the gravitational force is most significant.In this case we can write approximately: It is tempting to make a further step and write: but the parallel component of v satisfies equation ( 52) and thus equation ( 79) is not correct even to the first order in h 00 at it simply states that this component is fixed (remember we assume the force to be orthogonal to the light ray).Notice, however, that according to equation ( 56) the lensing phenomena is not affected by the parallel component, hence, no harm is done if we calculate this component to zero order instead of first order.It follows that for a static mass distribution it suffices for the purpose of deducing the lensing effect to solve the equation: which is just the Newtonian trajectory but with double the gravitational force.Now, the trajectory of gravitating test particles in a Newtonian point mass gravitational field is known and is given by the equation (Goldstein [41] equation 3.55) : in which r, θ are standard cylindrical coordinates.A static point gravitating body is assumed located at the origin of axis and generating a gravitational potential of the form: The quantity k according to Newtonian theory is k N = GM which follows from equation (75) and equation (67).This would also be the value according to General Relativity for a slowly moving test particle (see equation ( 77)).However, for a light ray we should use equation ( 80) instead and thus k E = 2k N = 2GM.The quantity l is a constant of motion (which can be thought as an angular momentum in the direction perpendicular to the plane of motion per unit mass) and is given by: The eccentricity e is given by the equation: which is dependent on another constant of motion E (which can be thought as the energy of the test particle per unit mass): A light ray coming from far away will have v = c at infinity while lim r−>∞ φ = 0 hence we may write: Thus: Consider a star is located far away from the sun at a distance r ∞ .Let the line connecting this star and the sun which is located at r = 0 coincide with the x-axis.Assume that the earth is also located along the x-axis and has an x coordinate of x earth = −r earth as depicted in figure 3. We consider a light ray orbit in which a star emits a photon which is detected on earth.At the closest point of the trajectory to the sun the photons travels at a direction which is purely in the θ direction and has no radial component.At this point the photon is a distance r p from the origin of axis and thus: This is inserted into equation (83) to yield: Hence: Now according to general relativity k = k E = 2GM: In which r s is the Schwarzschild radius defined in equation ( 60) and we assume the weak field approximation as described in equation ( 62), for the Newtonian theory we have k = k N = GM for which we obtain: In either relativistic or Newtonian theory the trajectory's shape is a rather extreme hyperbola (see Goldstein [41] p. 94).Now the starting and ending points of the trajectory are known: The photon trajectory starts at the star with the cylindrical coordinates (r, θ) = (r ∞ , 0) and ends on earth with the cylindrical coordinates (r, θ) = (r earth , π).Thus according to equation (81): It follows that we may also estimate l by: We may insert the above results into equation (81) and write We are now at a position to calculate the lensing angle given by equation (56): Now: x(θ) = r(θ) cos θ, y(θ) = r(θ) sin θ. (98) Calculating the derivative of the above quantities we have: The quantity dr dθ is deduced from equation ( 96) Inserting equation (100) into equation (99) and using some basic trigonometric identities and equation ( 93) leads to the following results: It is now easy to insert the expression of equation ( 101) into equation ( 97) and obtain a simple expression: Thus, at the star θ = 0 and thus the light ray satisfies tan α = 0, with a plausible physical solution of α = π.This means that initially the light ray is propagating parallel to the x-axis in the direction of earth.On the other hand as the light ray approaches earth, we have θ = π leading to: In which we used the fact that e 1 as follows from equation (91) and equation (92).Choosing the physical plausible positive sign and taking into account that the lensing angle is small we have: The Schwarzschild radius of the sun can be calculated from equation (60) from the sun's mass which is M sun 1.99 10 30 kg leading to r s 2950 m. r p was taken in the classical observation of Eddington [34] to be the Sun's radius, thus r p = 6.96 10 8 m, it now follows that: Eddington was able to show during a sun eclipse that the observed lensing angle is closer to α E than to α N thus adding an important empirical argument if favour of general relativity in addition to the theoretical arguments presented by Einstein.
To conclude this section we notice that the maximal h 00 in the solar system is according to equation (61): Thus justifying the linear approximation of equation ( 9) easily at least for the Solar system.

A suggested experiment
The lensing angle measured by Eddington is not the only possible lensing angle in the Solar system.To see that other angles are possible consider the case depicted in figure 4. Here, the light ray is emitted from a satellite on the same trajectory as the earth but in the other side of the sun.Consider a line connecting the earth and the satellite, and assume that this line coincides with the x-axis.The sun is located at the origin of axis and both the satellite and the earth are in equal distance from the sun.Now the starting and ending points of the trajectory are known also in this case: The photon trajectory starts at the satellite with the cylindrical coordinates (r, θ) = (r earth , 0) and ends on earth with the cylindrical coordinates (r, θ) = (r earth , π).Thus according to equation (81): The difference of the above two equations leads to: And their sum leads to: Hence l can be also evaluated as: We may insert the above results into equation ( 81) and write: The quantity dr dθ is deduced from equation ( 113) Inserting equation ( 113) and equation ( 114) into equation (99) will lead after some trigonometry to: It is now easy to insert the expression of equation ( 115) into equation ( 97) and obtain a simple expression: The ray should be launched from (r, θ) = (r earth , 0) at an angle: that is at according to equation (91) a small angle to the negative x direction of 1 e .And will arrive at earth at the angle: that is at a small angle 1 e below the negative x axis, causing the observer to see a deviation of the location of the satellite of the same angular magnitude, which is: The relativistic α E is expected rather than the Newtonian α N .Notice that in this case the angular deviation is half of the one expected for a distant star.Of course one can use a laser source in a short wavelength that is not common in the suns natural radiation and modulate the laser such that a matched filter can be easily constructed in the receiver side to minimize noise interference.

Beyond the Newtonian Approximation
So far we have considered only the case of Newtonian potential which neglects retardation phenomena.This approach seem to suffice for the solar system in which despite the sun's slow change in mass through the solar wind, retardation effect seem to be negligible.
However, in other gravitating systems such as galaxies (not to mention galaxy clusters) there mere size and the nature of their mass exchange with the environment leads to a situation in which retardation effects cannot be neglected.Indeed it was shown [15] that the peculiar shape of galactic rotation curves can be explained by the retardation phenomena.This raises the question on the effect of the same on gravitation lensing.Let us reiterate our main results.The velocity of a photon in a gravitational field is given to the first order of h 00 according to equation (52), equation (58), and equation ( 22): The duration R c for galaxies may be a few tens of thousands of years, but can be considered short in comparison to the time taken for the galactic density to change significantly.Thus, we can write a Taylor series for the density: By inserting Equations (122) into Equation ( 121) and keeping the first three terms, we will obtain: The Newtonian potential is the first term, the second term contributes only to v but no to v ⊥ , and the third term is the lower order correction to the Newtonian potential affecting v ⊥ in equation ( 121) : We recall that according to equation ( 56) it is only v ⊥ that affects the lensing angle, and thus despite the fact that second term G c ρ (1) ( x , t)d 3 x has a physical measurable effect on v it does not have any effect on the lensing phenomena to the first order in h 00 .
We underline again that we are not considering a post-Newtonian approximation in this paper (see section 3), in which matter travels at nearly relativistic speeds, but we will be considering the retardation effects and finite propagation speed of the gravitational field.We emphasize that taking v c 10 −3 << 1 (as for galaxies v 10 5 m/s see figure 1 while c 3 10 8 m/s) is not the same as taking R c << t r (with R being the typical size of a galaxy say about: R 3 10 20 m ⇒ R c 10 12 s) and t r is the typical time the mass of the same galaxy changes (to be discussed in section 9).Equation ( 58) can be rewritten in terms of a perpendicular "force" per unit mass such that: The total perpendicular force per unit mass is: Now if the lensing trajectory is far from the gravitating mass such that r ≡ | x| r ≡ | x | it follows that R r and R r ≡ x | x| .The force can now be written as: In the above M = ρd 3 x is the total mass of the gravitating body, which may be a galaxy or a cluster of galaxies, and M is the second derivative of the same.Now although M may be affected by many astrophysical processes we suggest that main contribution is the depletion of gas outsider the gravitating body (see section for a detailed discussion in section 9), hence M = −| M| < 0. It follows that according to equation (127): where the "dark matter" mass is defined as: we observe that the "dark matter" mass associated with lensing is the same as the "dark matter" mass associated with galactic rotation curves (see equation ( 106) of [15]), thus explaining the observational results of [27].

A Dynamical Model
As mass is accumulated in the galaxy or galaxy cluster, it must be depleted in the surrounding medium.This is due to the fact that the total mass is conserved; still, it is of interest to see if this intuition is compatible with a model of gas dynamics.For simplicity, we assume that the gas is a barotropic ideal fluid and its dynamics are described by the Euler and continuity equations as follows: where the pressure p(ρ) is assumed to be a given function of the density, ∂ ∂t is a partial temporal derivative, ∇ has its standard meaning in vector analysis and d dt is the material temporal derivative.We have neglected viscosity terms due to the low gas density.

General considerations
Let us now take a partial temporal derivative of equation (130) leading to: Using equation (130) again we obtain the expression: We divide the left and right hand sides of the equation by c 2 as in equation ( 126) and obtain: Since v c is rather small in galaxies and galaxy clusters it follows that 1 ∂t 2 is also small unless the density or the velocity have significant spatial derivatives.A significant acceleration ∂t resulting from a considerable force can also have a decisive effect.The depletion of available gas can indeed cause such gradients as we describe below using a detailed model.Taking the volume integral of the left and right hand sides of equation ( 134) and using Gauss theorem we arrive at the following equation: The surface integral is taken over a surface encapsulating the galaxy or galaxy cluster.This leads according to equation (129) to a "dark matter" effect of the form: Thus we obtain the order of magnitude estimation: In the above we define three gradient lengths: We can also write: in which the smallest gradient length will be the most significant one in terms of the "dark matter" phenomena.In the depletion model to be described below we assume that l ρ associated with density gradients is the shortest length scale.For galaxies we have v c 2 ≈ 10 −6 , hence the factor r l t should be around 10 6 to have a significant "dark matter" effect.A detailed model of the depletion process in galaxies leading to the desired second derivative of galactic mass is given in [15] and will not be repeated here.

Conclusions
In this paper we have deduced from general relativity a linear approximation.Under the said linear approximation we have solved Einstein field equations in term of retarded solutions.Those where used to derive the trajectory of a light ray (photon) both in the direction parallel to its original direction and perpendicular to it.Our current approach was compared to the approach based on the Schwarzschild metric [35] and was shown to be better in the sense that it is applicable to general mass distributions including ones that are changing in time.The light ray equations is the presence of a static point mass allowed us to re-derive the classical light deflection of Einstein and Eddington [24] and suggest a new experiment.This was followed by a discussion on the retardation effects on light ray trajectories, deriving an expression for "dark matter" which is the same as the one obtained in [15] for slowly moving bodies.Thus justifying results reported in the literature of the equivalence of "dark matter" for both galaxy rotation curve and gravitational lensing.This is followed by a discussion on the physical requirements needed in order to derive "dark matter" effects from retardation and a detailed model showing how those requirements are satisfied in a galactic scenario.
Lorentz symmetry invariance does not allow action at a distance potentials and forces, but retarded solutions are allowed.Retardation is significant for large distances and large second derivatives.It should be emphasized that the retardation approach does not require that velocities in, v in the gravitating body are high; in fact, galactic & galactic cluster bodies (stars, gas) move slowly with respect to the speed of light-thus the quantity v c 1. Typical velocities in galaxies are 100 km s (see Figure 1), which makes v c to be about 0.001 or smaller.However, every gravitational system, even if it consists of subluminal entities, has a retardation distance, above which retardation cannot be neglected.Natural systems, for example a star or a galaxy and even a galactic cluster, exchanges mass with its environment.The sun loses mass through solar wind and galaxies accrete matter which originate in the intergalactic medium.Thus all natural (gravitational) systems have a finite retardation distance.This leads to quantitative inquiry: what is the actual size of the retardation distance?The modification of the solar mass is quite small and thus the retardation distance of the solar system is extremely large, we can thus neglect retardation within the solar system.On the other hand, for the M33 galaxy, velocities indicate that retardation cannot be neglected.The retardation distance was calculated in [15] to be roughly R r = 4.54 kpc for M33; other galaxies of different types have shown similar results [14].We demonstrated, in Section 9, that this does not require a high velocity of gas or stars and is perfectly consistent with current observational knowledge of galactic and extragalactic dynamics.
We underline, that if extra galactic mass is abundant (or totally consumed), M 0 and the retardation force vanish.As was reported [43] for NGC1052-DF2.
We emphasize that the terms in the GR equations responsible for gravitational radiation recently discovered are also the cause for the peculiar shape of the rotation curves of galaxies and gravitational lensing.The approximation used here is not a far field approximation but a near field one.Indeed, the expansion (123), being second order, is only valid up to limited radii: This is reasonable since the extension of the rotation curve in galaxies and the distance of lensing trajectories from galaxies is the same order of magnitude as the size of the galaxy.The case in which the dimensions of the source is much smaller than the distance to the observer will result in a different valid approximation to (13), leading to the famous quadruple expression of gravitational radiation, as derived by Einstein [46] and verified (indirectly) in 1993 by Russell A. Hulse and Joseph H. Taylor.The observation of the Hulse-Taylor binary pulsar has given the first (indirect) evidence of gravitational waves [47].On 11 February 2016, the LIGO and Virgo Collaboration announced that they made the first (direct) observation of gravitational waves.The observation was made earlier, on 14 September 2015, using the LIGO detectors.The gravitational waves were caused by the merging of a binary black hole system [48].Thus, here we discuss only a near-field application of gravitational radiation contrary to previous works discussing farfield results.
Unfortunately no direct measurement of the second temporal derivative of the galactic mass is available.What is available is the remarkable fit between the retardation theoretical velocity and the observed galactic rotation curve, as can be seen in Figure 1, this constitutes indirect evidence of the total mass second derivative.Competing theories like dark matter do not supply any direct observational evidence either.Despite the work of many people and the investment of a large financial resources, there is no evidence of dark matter.Occam's razor postulates that when theories compete, the one that makes less ontological assumptions about the physical existence of an exotic form of matter, wins.Retardation theory assumes only baryonic matter and a large second temporal derivative of mass.
Problems related to dark matter, such as the core-cusp problem which refers to the difference between the dark matter density profiles of galaxies and the density profiles predicted by N-body simulations.Almost all simulations form dark matter halos, which posses "cuspy" dark matter mass distributions, with density increasing rather steeply at small radii, while the rotation curves of most dwarf galaxies suggest that they have flat central dark matter profiles ("cores").This does not occur in the retardation theory which does not require dark matter.One cannot consider flat or sharp profiles of dark matter distribution if dark matter is not there.The persistent difficulties with dark matter's dynamics strengthen the claim that dark matter is not needed and the gravitational lensing characteristics attributed to dark matter should be attributed to retardation.
To conclude, we would like to mention the significant theory of conformal gravity put forward by Mannheim [3,4].The current retardation approach leads approximately to a Newtonian potential and in addition a linear potential.Such potential types can be derived from conformal gravity.On purely phenomenological grounds, rotation curve fits for linear plus Newtonian potentials have already been published.While those fits are very good, they had to treat the coefficient of the linear potential as a variable that changed from galaxy to galaxy; this element contradicts conformal gravity in which the coefficient of the linear potential is a universal constant.This can be explained in the framework of retardation, in which a M/M depends on the dynamical particular conditions of every galaxy.Indeed, M/M has a theoretical reason, lacking in a pure phenomenological approach.The work of Mannheim [2][3][4], is related to conformal gravity, which is different from GR, and thus has to justify other results of GR (Big Bang Cosmology, etc.).We also underline that retardation per se does not contradict conformal gravity and, both effects may exist, although the principle of Occam's razor forbids us to add new universal constants if the existing ones suffice to explain observations.Retardation theory's approach is minimalistic (it satisfies the Occam's razor rule), and does not affect observations that are beyond the near-field regime and, thus, does not contradict GR theory and its observational consequences (nor with Newtonian theory, as the retardation effect is negligible for "small" distances).The perfect fit to the rotation curve and gravitational lensing is achieved with a single parameter and we do not adjust the mass to light ratio in order to improve our fit as is done by other authors.Retardation effects beyond gravity, in particular with respect to electromagnetic theory were studied in [51][52][53][54].
In this paper we study the gravitational lensing scenario showing that retardation effects lead to the same "dark matter" mass for lensing as for galactic rotation curves.
The current paper does not discuss dark matter in a cosmological context; this is left for future work.We mention, however, that on the cosmological scales it is not enough to invoke dark matter, but one must also consider dark energy as described in the ΛCDM model.
The CMB anisotropy spectrum that was observed precisely using WMAP in 2003-2012, and even to higher precision by the Planck spacecraft in 2013-2015 is in agreement with the ΛCDM model [55,56].
Moreover, small anisotropies of the homogenous universe grew gradually and condensed the homogeneous material into stars, galaxies and larger structures.Since ordinary matter is affected by radiation, which is dominant at very early times.It follows,that its density perturbations are washed out and unable to condense [57].Thus it is suggested that if there is only baryonic matter, there would not have been sufficient duration for perturbations to grow into galaxies and clusters seen.
Dark matter is assumed to provides a solution to this problem because it does not interact with electromagnetic radiation.Therefore, its perturbations can grow fast.The resulting gravitational potential attracts ordinary matter collapsing later, speeding significantly the structure formation process [57,58].
On the other hand we know [15] that the Newtonian gravity is only a part of the gravitational force and at large distances retardation forces prevail.Thus what is attributed to dark matter attraction in terms of density perturbation growth may be attributed to retardation.
The CAMB anisotropy spectrum as well as the distant Super Novae data may be explained by a more thorough perturbation analysis of the Friedman Robertson Walker metric, preliminary results in this direction are given in [59].

Figure 1 .
Figure 1.Velocity curve for M33.Observational points were obtained by Dr. Michal Wagman, a former PhD student under my supervision, using[40]; the full line describes the rotation curve, which is the sum of the dotted line, describing the retardation contribution, and the dashed line, which is Newtonian.

Figure 2 .
Figure 2. Light travelling toward the observer is bent due to the gravitational field of a massive object, thus a distant star appears to the observer at an angle α with resect to its true location.The tangent of alpha is the velocity of light perpendicular to the original direction v ⊥ divided by the velocity of light in the original direction v .