Effects of Higher Order Retarded Gravity on Galaxies †

: In a recent paper, we have a shown that the ﬂattening of galactic rotation curves can be explained by retardation. However, this relies on a temporal change of galactic mass. In our previous work, we kept only second order terms of the retardation time in our analysis, while higher terms in the Taylor expansion were not considered. Here, we consider analysis to all orders and show that, indeed, a second order analysis will sufﬁce, and higher order terms can be neglected.


Introduction
The dynamics of large scale structures is inconsistent with Newtonian mechanics. This was noted in the 1930s by Fritz Zwicky [1], who pointed out that if more (unseen) mass were present, one would be able to solve the apparent contradiction. The phenomena was also observed in galaxies by Volders [2], who has shown that star trajectories near the rim of galaxies do not move according to Newtonian predictions, and later corroborated by Rubin and Ford [3][4][5] for spiral galaxies.
In a series of papers, we have shown that these discrepancies can be shown to result from retarded gravity as dictated by the theory of general relativity [6][7][8][9][10]. Indeed, in the absence of temporal density changes, retardation does not effect the gravitational force. However, density is not constant for galaxies; in fact, there are many processes that change the mass density in galaxies over time. Mass accretion from the inter galactic medium and internal processes such as super novae leading to super winds [10] modify the density. In addition to these local processes, there is a cosmological decrease in density due to the cosmic expansion. However, the later process is many orders of magnitude smaller than the former.
In previous analysis [6][7][8][9][10], the corrected gravitational force was evaluated assuming a second order approximation in the retardation time R c , neglecting higher order terms without justification. Here, we take into account all orders and show that a second order approximation is indeed sufficient.

Linear GR
Only in cases of extreme compact objects (black holes and neutron stars) and the very early universe do we consider the solution of the full non-linear Einstein equations [6]. In most cases of astronomical interest (including the galactic case), a linear approximation to those equations around the flat Lorentz metric η µν is used, such that One thendefines the quantitȳ h µν = h µν for non diagonal terms. For diagonal terms, It was shown ( [11] page 75, exercise 37, see also [12][13][14]) that for a proper gauge, the Einstein equations are Equation (4) can be solved such that [15] h µν ( In [16][17][18][19], we explain why the symmetry between space and time is broken, which justifies the use of different notations for space and time. 4G c 4 3.3 × 10 −44 is a tiny number; hence, in the above calculation, one can take T µν to the zeroth order in h αβ . We now evaluate the affine connection in the linear approximation: Notice that the affine connection has first order terms in h αβ ; hence, to the first order Γ α µν u µ u ν appearing in the geodesic equation, u µ u ν must be taken to the zeroth order, in which For velocities much smaller than the speed of light in vacuum, Hence, the current paper does not discuss the post-Newtonian approximation, in which matter travels at speeds close to the speed of light, but we do consider the retardation effects which are due to the finite propagation speed of the gravitational field. We emphasize the assumption that v c << 1 is not the same as stating R c << 1 (with R being the typical size of a galaxy), since Now, since in galaxies, R v is huge ( R v 10 15 s), it follows that v c can be dismissed but not R c , for which R c 10 12 s. Inserting Equations (6) and (8) in the geodesic equation, we arrive at the approximate equation: Taking a standard matter T µν , assuming ρc 2 p, and taking into account Equation (8), we arrive at T 00 = ρc 2 , while the remaining tensor components are much smaller. Therefore, h 00 is larger than the other components ofh µν . Notice that it is not possible to deduce from the magnitudes of quantities that a similar difference exists between the derivatives of those quantities. Gauge conditions in Equation (4) Thus, the zeroth derivative ofh 00 (which contains a 1 c ) is of similar order to the spatial derivative ofh 0a . In addition, the zeroth derivative ofh 0a (see Equation (10)) is of similar order to the spatial derivative ofh ab . However, we can compare the spatial derivatives of h 00 andh ab and conclude that the former is larger. Taking into account Equation (3) and the above consideration, we write Equation (10) as dv a dt c 2 4h 00, Thus, the gravitational potential φ can be estimated using Equation (5): and F is the force per unit mass. In the case that the mass density ρ does not depend on time, we may use the Newtonian instantaneous action at a distance. Note that it is improbable that ρ is static for a galaxy, as it accretes intergalactic medium gas.

Retardation Effects Beyond the Newtonian Approximation
The duration R c may be tens of thousands of years but may be short with respect to the duration in which the galactic density changes considerably. Thus, we write a Taylor series for the density: By inserting Equation (14) into Equation (13), we will obtain The Newtonian potential is the first term, the second term has null contribution, and the third term is the lower order correction to the Newtonian theory: We will show later that n > 2 terms can be neglected; thus, the total force per unit mass can be approximated as F N , first introduced by Newton, is attractive; however, the retardation force F r can be either attractive or repulsive. Newtonian force decreases as 1 R 2 ; however, the retardation force does not depend on distance as long as the Taylor approximation given in Equation (14) holds. Below a certain distance, the Newtonian force dominates, but for larger distances, the retardation force has the upper hand. Newtonian force can be neglected for distances significantly larger compared to the retardation distance: ∆t is a duration associated with the second order derivative of the density ρ. For R R r , retardation can be neglected, and only Newtonian forces need to be considered; this is the situation in the solar system. As the galaxy accretes intergalactic gas, the galactic mass becomes larger, and thusṀ > 0; however, the intergalactic gas is depleted, and thus the rate at which the mass is accreted decreases resulting inM < 0. Hence, we have an attractive retardation force.

Higher Order Terms
Comparing Equation (31) to Equation (82) of [9], it follows that Hence: And also Thus, according to Equation (15), we have the following correction to the retardation potential: The deviation from the second order approximation is more pronounced for large r, for which R r, which is the case we consider here; thus: Now, using the well-known identity: We may write Equation (23) as a closed expression instead of an infinite sum: For r cτ, it is quite clear that the term in the parenthesis of Equation (25) vanishes, since Hence, φ (n>2) can be neglected if indeed r cτ for the relevant measurements of the M33 rotation curve, that is, up to about r < 20 kpc. Now, τ is dependent, according to Equation (81) of [9], on the density gradient of the inter galactic medium (IGM) and the typical velocity in this medium. Although these values are not known precisely, we may assume that v z ∼ 100 km /s, and the typical gradient is the same as the gradient of the optical disk luminosity, that is, 1 k ∼ 0.1 kpc. Thus, τ ∼ 10 6 years, and τc ∼ 300 kpc, making the second order approximation used so far reasonable.

Conclusions
The phenomena of retardation is ubiquitous in physics and follows directly from the Lorentz symmetry group. Hence, any system that is invariant under the Lorentz transformation will exhibit retardation phenomena. These include physical systems related to classical electromagnetism [20][21][22][23] and general relativity [6][7][8][9], but also to other Lorentz invariant theories such as conformal gravity [24][25][26].
Dark matter being a major candidate for explaining galactic rotation curves has only a slim chance of being found, given that accelerator experiments, such as using the Large Hadron Collider, have been unable to find any super symmetric particles, not only of the community's favorite form of dark matter, but also of the form that is mandated in string theory, a theory that also suggests a quantized version of Einstein gravity.
We have shown that at least on the galactic scale, dark matter is not needed [6][7][8][9][10], as the needed dynamics can be explained by a retarded gravitational potential when a near field approximation is used. We remark that the analysis of far field leading to gravitational waves [27] was corroborated in recent years by observations [28,29].
A justification for the second order Taylor series approximation which we used in previous works is given here for the first time, showing that indeed higher order terms can be safely neglected.
Finally, we mention other approaches to the galactic rotation curves problem, which suggest overcoming the problem by changing the laws of gravity. Such approaches are Milgrom's MOND [30] and Mannheim's conformal gravity [24][25][26]. Unfortunately, these approaches seem to contradict both general relativity (which is supported by a large body of observational evidence) and recent observations. Indeed van Dokkum et al. [31] have shown that there are galaxies with Newtonian rotation curves, excluding the possibility of a universal modification to the laws of gravity which prevail in every galaxy. Thus, either van Dokkum's galaxy is devoid of "dark matter", or it has a small retardation depletion effect, as the gas around it has not yet depleted (or fully depleted).
Given the negative results from accelerator experiments regarding dark matter, retardation theory seems to be the only valid option.