Lorentz Symmetry Group, Retardation, Intergalactic Mass Depletion and Mechanisms Leading to Galactic Rotation Curves

The general theory of relativity (GR) is symmetric under smooth coordinate transformations, also known as diffeomorphisms. The general coordinate transformation group has a linear subgroup denoted as the Lorentz group of symmetry, which is also maintained in the weak field approximation to GR. The dominant operator in the weak field equation of GR is thus the d'Alembert (wave) operator, which has a retarded potential solution. Galaxies are huge physical systems with dimensions of many tens of thousands of light years. Thus, any change at the galactic center will be noticed at the rim only tens of thousands of years later. Those retardation effects are neglected in the present day galactic modelling used to calculate rotational velocities of matter in the rims of the galaxy and surrounding gas. The significant differences between the predictions of Newtonian instantaneous action at a distance and observed velocities are usually explained by either assuming dark matter or by modifying the laws of gravity (MOND). In this paper, we will show that, by taking general relativity seriously without neglecting retardation effects, one can explain the radial velocities of galactic matter in the M33 galaxy without postulating dark matter.


Introduction
The general theory of relativity (GR) is known to be symmetric under general coordinate modifications (diffeomorphism). This group of general coordinate transformations has a Lorentz subgroup, which is maintained even in the weak field approximation. This is manifested through the field equations containing the d'Alembert (wave) operator, which has a retarded potential solution.
From an observational point of view, it is well known that GR is verified by many observations. However, at the present time, the standard Newton-Einstein gravitational theory stands at something of a crossroads. It simultaneously has much in its favor observationally, while, at the same time, it has some very disquieting challenges. The observational successes that it has achieved in both astrophysics and cosmology have to be tempered by the fact that the theory needs to appeal to two as yet unconfirmed ingredients, dark matter and dark energy, in order to achieve these successes. The dark matter problem has not only been with us since the 1920s and 1930s (when it was initially known as the missing mass problem), but it has also become more and more severe as more and more dark matter has had to be introduced on larger and larger distance scales as new data have come online. Moreover, extensive-now 40-year-underground and accelerator searches have failed to find any of it or establish its existence. The dark matter situation has become even more dire in the last few years as the Large Hadron Collider has failed to find any super symmetric particles, not only of the community's preferred form of dark matter, but also the form of it that is required in string theory, a theory that attempts to provide a quantized version of Newton-Einstein gravity.
While things may still eventually work out in favor of the standard dark matter paradigm, the situation is disturbing enough to warrant consideration of the possibility that the standard paradigm might at least need to be modified in some way if not outright replaced. The present proposal sets out to seek such a modification. Unlike other approaches such as Milgrom's MOND, Mannheim's Conformal Gravity or Moffat's MOG, the present approach is, in a sense, the minimalist one adhering strictly to the basic scientific principle dictated by Occam's razor. It seeks to replace dark matter by effects within standard General Relativity itself.
In 1933, Fritz Zwicky noticed that the velocities of a set of Galaxies within the Comma Cluster are much higher than those predicted by the virial calculation that assumes Newtonian theory [1]. 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 entire cluster. Volders, in 1959, indicated that stars in the outer rims of the nearby spiral galaxy M33 do not move "correctly" [2]. It is the result of the virial theorem coupled with Newtonian Gravity which implies that MG/r ∼ Mv 2 , that is to say, the expected rotation curve should at some point decrease as 1/ √ r. This was well established during the seventies when Rubin and Ford [3,4] demonstrated that, for a large sample of spiral galaxies, this behavior can be considered a general feature: velocities at the our rim of the galaxies do not decrease-rather, in a general case, they attain a plateau at some velocity, which is different for each galaxy. In what follows, we will show that such effects can be deduced from GR if retardation effects are not neglected. The derivation of the retardation force described in previous publications [5][6][7] is repeated for completeness. However, a fit of the theory to the M33 rotation galaxy is given for the first time. It should be stressed that the current approach does not require that velocities v are high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal-in other words, the ratio of v c 1. Typical velocities in galaxies are 100 km/s, which makes this ratio 0.001 or smaller. However, one should consider the fact that every gravitational system, even if it is made of subluminal bodies, has a retardation distance, beyond which the retardation effect cannot be neglected. Every natural system, such as stars and galaxies and even galactic clusters, exchanges mass with its environment. For example, the sun loses mass through solar wind and galaxies accrete gas from the intergalactic medium. This means that all natural gravitational systems have a finite retardation distance. The question is thus quantitative: how large is the retardation distance? The change in the mass of the sun is quite small and thus the retardation distance of the solar system is quite large, allowing us to neglect retardation effects within the solar system. However, for the M33 galaxy, the velocity curve indicates that the retardation effects cannot be neglected beyond a certain distance smaller than 14,000 light years; similar analyses for other galaxies of different types have shown similar results [8,9]. We demonstrate, in Section 6, using a detailed model, that this does not require a high velocity of gas or stars in or out of the galaxy and is perfectly consistent with the current observational knowledge of galactic and extra galactic material content and dynamics.

General Relativity
The general theory of relativity is based on two fundamental equations, the first being Einstein equations [10][11][12][13]: G µν stands for the Einstein tensor, T µν indicates the stress-energy tensor, 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 α (s) are the coordinates of the particle in spacetime, s is a typical parameter along the trajectory that for massive particles is chosen to be the length of the trajectory, u µ = dx µ ds is the µ-th component of the 4-velocity of a massive particle moving along the geodesic trajectory s (increment of x per unit length s) 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, p 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: and the affine connection: The affine connection serves to 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
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 [5]. 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 µν = 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 ( [10] page 75, exercise 37, see also [11][12][13]) 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 rotation curve frame of reference for galaxies, matter is approximately at rest. Equation (12) can always be integrated to take the form [14] (For reasons why the symmetry between space and time is broken, see [15,16] The factor before the integral is small: 4G c 4 3.3 × 10 −44 ; hence, in the above calculation one can take T µν , which is zero order in h αβ . Let us now calculate the affine connection in the linear approximation: The affine connection has only first order terms in h αβ ; hence, to the first order Γ α µν u µ u ν appearing in the geodesic, u µ u ν is of the zeroth order. In the zeroth order: For non relativistic velocities: Hence, we will not be considering a post-Newtonian approximation in this paper, 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 underline that taking v c << 1 is not the same as taking R c << 1 (with R being the typical size of a galaxy) since: and, in galaxies, R v is a very big number ( R v 10 15 seconds); thus, v c can be neglected but not R c , in which R c 10 12 seconds. By inserting Equations (14) and (16) in the geodesic equation, we arrive at the approximate form: Let us now look at Equation (3). We assume ρc 2 p and, taking into account Equation (16), we arrive at T 00 = ρc 2 , while other tensor components are significantly smaller. Thus,h 00 is significantly larger than other components ofh µν . One should notice that it is not possible to deduce from the magnitudes of quantities that such a difference exists between their derivatives. In fact, by the gauge condition in Equation (12):h α0, Thus, the zeroth derivative ofh 00 (which contains a 1 c ) is the same order as the spatial derivative ofh 0a and the zeroth derivative ofh 0a (which appears in Equation (18)) is the same order as the spatial derivative ofh ab . However, it we can compare spatial derivatives ofh 00 andh ab and conclude that the former is larger than the later. Using Equation (11) and taking into account the above consideration, we may write Equation (18) as: Thus, φ is the gravitational potential of the motion which can be calculated using Equation (13): and F is the force per unit mass. If the mass density ρ is static, we are in the realm of the Newtonian instantaneous action at a distance. We point out that it is unlikely that ρ is static, as a galaxy will obtain mass from the intergalactic medium.

Beyond the Newtonian Approximation
The retardation time R c 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 (22) into Equation (21) and keeping the first three terms, 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: The total force per unit mass: While the Newtonian force F N is always attractive, we notice that the retardation force F r can be attractive or repulsive. The Newtonian force decreases as 1 R 2 , but the retardation force is independent of distance as long as the Taylor approximation given in Equation (22) is valid. Below a certain distance, the Newtonian force will be dominant, but for higher distances the retardation force becomes larger. Newtonian force should be neglected for distances significantly larger than the retardation distance, defined as: ∆t is the typical duration in which the density ρ changes; this will be defined more accurately later. On the other hand, for R R r , the retardation effect can be neglected and only Newtonian forces should be considered; this is probably the situation in the solar system. It should be stressed that the existence of R r does not require that velocities v are high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal. In other words, the ratio of v c 1. Typical velocities in galaxies are 100 km/s, which makes this ratio 0.001 or smaller. However, one should consider the fact that every gravitational system, even if it is made of subluminal bodies, has a retardation distance, beyond which the retardation effect cannot be neglected. Every natural system, such as stars and galaxies and even galactic clusters, exchanges mass with its environment. For example, the sun loses mass through solar wind and galaxies accrete gas from the intergalactic medium. This means that all natural gravitational systems have a finite retardation distance. The question is thus quantitative: how large is the retardation distance? This question will be further addressed in Section 5. For large distances, r = | x| → ∞, such thatR x | x| ≡r; thus, we obtain: As the galaxy attracts intergalactic gas, its mass becomes larger and thereforeṀ > 0; however, as the intergalactic gas is depleted, the rate at which the mass is accumulated must decrease and thereforeM < 0. Thus, in the galactic case: and the retardation force is attractive.

Rotation Curves
To calculate the non-asymptotic rotation curve, we note that the square of the azimuthal (orbital) velocity, v θ divided byr, is proportional to the derivative of the gravitational potential byr, as shown in Equation (68). We use the cylindrical coordinatesr, θ, z, in which z = 0 is the galactic plane. The gravitational potential will be evaluated using Equation (23), which is composed of both the Newtonian φ N and retardation φ r parts, such that: To do this, we first need a density profile describing the mass distribution in and near the galactic plane.

General Considerations
We assume a density distribution of the form: Although the density profiles ρ a and ρ b need not have similar forms, we will assume, for simplicity, that they do. Hence, ρ a = ρ b and, by defining g(t) = 1 + g (t), we obtain a density distribution of the form: hence, which leads to the result: here, we define: Inserting Equations (31) and (34) into Equation (24) leads to the following form of the retardation potential: Now, we introduce the dimensionless quantities: in terms of a typical density ρ c and a typical scale R s . From now on, we will consider all time-dependent quantities as being given at the time the galaxy is being observed. This has nothing to do with the present situation in a certain galactic system that may be millions or billions year in the future with respect to the data available to us, which of course is also the result of retardation. The dependence of the temporal variable t will thus be omitted. Using Equations (37) and (32) will take the form: We now define the dimensionless constant: In terms of which: Hence, φ r can be written as: χ is a dimensionless function that satisfies: Similarly, one can write the Newtonian potential as: ψ is a dimensionless function that also satisfies:

M33 Density Profile
In order to calculate the retardation and Newtonian gravitational potentials and hence the rotation curve through Equation (68), we must know the density distribution in the galactic plane. This is usually done as follows: for the radial distribution, we will use the luminosity distribution and assume a proportionality constant known as the mass to light ratio. Such radial density distribution is given in Figure 1, which is based on the work of Corbelli [17] (see also Rega and Vogel [18]) and its logarithm is fitted to a polynomial of order five, as described in Equation (45): For the z direction distribution orthogonal to the galactic plane, we will assume a Gaussian profile with a typical width of σ = 0.2 kpc. The mass distribution can be written as: in which we assume cylindrical symmetry. The three-dimensional galactic mass distribution is depicted in Figure 2.

M33 Rotation Curve
To obtain the rotation curve we need to evaluate the functions ψ(r) and χ(r), which describe the non-trivial contribution of the Newtonian and retardation potentials, respectively, see Equations (41) and (43). We will assume, for simplicity, that the galaxy is cylindrically symmetric and that rotation curves are evaluated at the galactic plane. Thus, according to Equation (43), we need to numerically evaluate the integral: in the above: and:r The result of the numerical evaluation of ψ(r) is given in Figure 3; we notice that, due to the cylindrical symmetry, it is enough to evaluate the function for the azimuthal angle θ = 0 as the result is similar for each azimuthal angle. As can be easily seen, the ψ function of M33 converges to one for large distances, as expected from Equation (44). If we ignore the retardation potential contribution Equation (29), then Equation (43) will lead to the rotation curve depicted in Figure 4, which decreases for large distances. In the above, we assumed a galactic (baryonic) mass of 10 10 solar masses (2 × 10 40 kg) [17]. As the true galactic rotation curve does not decrease, we need to evaluate the retardation potential Equation (41) as well. For this, we evaluate the χ function as follows: The result of the evaluation is given in Figure 5, where it is obvious that the χ function approaches unity for large distances from the galactic center. The retardation potential and thus the retardation contribution to the velocity cannot be calculated without knowledge of the retardation distance R r . However, this can be obtained easily by fitting the observational galactic rotation curve, as demonstrated in Figure 6, which yields a best fit for R r = 4.54 kpc and a retardation time of t r = 14, 818.7 years. Looking at Figure 6, this seems quite reasonable. We deduce from the above and from Equation (35) that the galactic mass second derivative is as follows: |M| = M t 2 r 9.12 × 10 16 kg/s 2 We underline that this is based only on the current estimation of the second derivative of mass; we make no claim about the past or future values ofM, nor is there any claim in this section on the value ofṀ at any time or the value of M in the past or the future. It is obvious that such questions involve an understanding of the mass exchange between the galaxy and the intergalactic medium, as described in Section 6.

A Dynamical Model
The mass accumulation model described in the previous section is based on a fitting of the second derivative of the galactic mass to the galactic rotation curve. It is intuitively obvious that, as mass is accumulated in the galaxy, it must be depleted in the intergalactic 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. For simplicity, we assume axial symmetry; hence, all variables are independent of θ. Moreover, the mass influx coming from above and below the galaxy is much more significant compared to the influx coming from the galactic edge. This is due to the large difference in the galaxy surfaces perpendicular to the z axis compared to the area of its edge. The area of the surface of the galaxy which is perpendicular to the z axis is: in which S z is the total surface area of the galaxy perpendicular to the z axis, S z+ is the upper area of the surface of the galaxy perpendicular to the z axis, S z− is the lower area of the surface of galaxy perpendicular to the z axis and r m is the galactic radius (see Figure 7). The area of the surface of the galactic edge with thickness ∆z is: Thus, the ratio of the surface area is: Typical values of ∆z are about 0.4 kilo parsec and r m is about 17 kilo parsec (for M33), giving an area ratio of about 1%. In such circumstances, the edge mass influx is less important and we can assume a velocity field of the form: z andθ are unit vectors in the z and θ directions, respectively. The influx is described schematically in Figure 8. In this case, the continuity Equation (52) will take the form: By defining the quantity: and using the above definition, Equation (58) takes the form: Assuming, for simplicity, that v z is stationary and defining the auxiliary variable t z : we arrive at the equations: ∂γ ∂t This equation can be solved easily, as follows: for the function f (x), which is fixed by the density initial conditions and the velocity profile.
Let us now turn our attention to the Euler Equation (53); for stationary flows, it takes the form: According to Equation (57): Now, by writing Equation (64) in terms of its components, we arrive at the following equations: It is usually assumed that the radial pressure gradients are negligible with respect to the gravitational forces and thus we arrive at the equation: v 2 θ r ∂φ ∂r , As for the z-component part of the Euler Equation (64), it can be easily written in terms of the specific enthalpy w(ρ) = dP ρ in the form: ∂ ∂z We recall that ρ depends on v z through Equations (59) and (63): As both the specific enthalpy and the gravitational potential are dependent on the density, Equation (69) turns into a rather complicated nonlinear integral equation for v z . However, many galaxies are flattened structures; hence, it can thus be assumed that the pressure z gradients are significant as one approaches the galactic plane. We will thus assume, for the sake of simplicity, that the pressure gradients balance the gravitational pull of the galaxy and thus v z is just a function of r, in which case the convective derivative of v z vanishes. The above assumption holds below and above the galactic plane, but not at the galactic plane itself. This suggests the following simple model for the velocity v z (see Figure 8): in which |v z | is a known function ofr. The velocity field is discontinuous at the galactic plane due to our simplification assumptions, but, of course, need not be so in reality. We also assume, for simplicity, that the velocity field |v z | is constant forr < r m and vanishes forr > r m . According to Equation (63), the time-dependent density profile is fixed by the initial density conditions. In this section, we will deal with the density profile outside the galactic plane and will leave the discussion of the density profile in and near the galactic plane to the previous section. We consider an initial density profile as follows: in which the rectangular function re(z) keeps the exponential function from diverging. The density profile is depicted in Figure 9. We assume that ρ 2 is negative and thus the density becomes dilute at distances far from the galactic plane. As v z is constant both above and below the galactic plane, t z = z v z up to a constant. Therefore, it is easy to deduce from Equation (63) the functional form of f (β) (β = −t z ): Hence, according to Equation (70), the time-dependent density function for matter outside the galactic plane is obtained: The density of matter outside the galactic plane will vanish for t > t m = z i |v z | ; hence, we will discuss only the duration of t < t m . Let us look at the mass contained in the cylinder defined by the galaxy (see Figure 10) and let us assume that the total mass in that cylinder is M T . Now, the mass outside the galactic disk will be: Hence, the mass in the galactic disk is: Moreover, the galactic mass derivatives are: in which: Then, calculating the second derivative of M o (t) and using Equation (77) leads to the result: We then denote: Thus: It should be stressed that the current approach does not require that the velocity |v z | is high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal-in other words, the ratio of |v z | c 1. The typical velocities in galaxies are 100 km/s, which makes this ratio 0.001 or smaller. However, we stress again the fact that every gravitational system, even if it is made of subluminal bodies, has a finite retardation distance, beyond which the retardation effect cannot be neglected. As demonstrated above, a galaxy exchanges mass with its environment. This means that a galaxy has a finite retardation distance. The question is thus quantitative: how large is the retardation distance? For the M33 galaxy, the velocity curve indicates that retardation effects cannot be neglected beyond a certain distance, which is calculated in Section 5 to be roughly 14,000 light years; similar analyses for other galaxies of different types have shown similar results [8]. We demonstrate here, using a detailed model, that this does not require a high velocity of gas or stars, in or out of the galaxy, and is perfectly consistent with the current observational knowledge of galactic and extragalactic material content and dynamics. Equation (81) means that we must have λ 2 < 0 according to Equation (28) in order to assure an attractive force. Next, we calculateṀ(t); using Equations (77) and (78), we obtain: Hence:Ṁ Thus, λ 1 is: Inserting Equations (85) and (81) into Equation (83) leads to: Finally, by combining Equations (76), (78), (81) and (85) and noticing that: we arrive at:

The Mass Formula
A naive observer may speculate thatM, given in Equation (51), is constant over time; however, this is contradicted both by the dynamical model as presented in Equation (88) and by the following simple argument.
Let us assume thatM is constant. This leads to: Let us assume that there is no significant mass accumulation at the current time T, that is,Ṁ(T) = 0. In this case: Assuming that the age of the galaxy is about the age of the universe: T 13 × 10 9 years = 4.1 × 10 17 s, the initial mass accumulation rate is aboutṀ(0) 3.74 × 10 34 kg s −1 . By integrating Equation (89), we arrive at an equation for the mass of the galaxy at any time: Assuming that, initially, there was no significant amount of mass M(0) = 0 and taking into account Equation (90) (as for Equation (90), it should be once again recalled that the second derivative of M by t is negative) we arrive at: By plugging in the mass accumulation decrease rate, we arrive at M(T) 7.66 × 10 51 kg, which is clearly 11 orders of magnitude greater than the known mass of the galaxy. On the other hand, if we assume that M(T) = M is the current mass of the galaxy, this will lead to |M| 2.37 × 0 5 kg s −2 , which is 11 orders of magnitude less than what is needed to explain the current galactic rotation curves. This indicates that |M| has increased considerably over time due to the depletion of the surrounding gas and that retardation forces are less significant in the early stages of galactic formation. Indeed, in young galaxies, the effect of retardation seems insignificant [19], while, for older galaxies, it seems like somebody is pressing hard on the brakes of mass accumulation. The changes, by many orders of magnitude, in |M| suggest exponential growth in agreement with Equation (88).
If we assumeṀ(t) > 0, it follows from Equation (86) that there is a maximal duration for the above model to be valid: This is the time it takes the intergalactic gas to deplete (see Section 6) and thus mass accumulation stops. We shall now partition the mass into a linear and exponential growing part as follows: By dividing Equation (94) by |M(t)| and using the definition of Equation (35), we have: We define the slow retardation time as: the square of this retardation time is non-linearly dependent on time. Hence: It is obvious that t rs > τ and it is also obvious that for any t < t max1 , then t rs e − t 2τ > τ because the total mass is not reduced, only mass accumulation slows down (in other words, we assume that the galaxy does not eject matter, thus reducing its mass). Equation (98) shows that there is no simple relation between the retardation time of the galaxy and the typical time τ associated with the second derivative of mass.
We conclude this section by deriving constraints from the assumption that the apparent luminosity and rotation curve have not changed in a detectable way since observations of the M33 galaxy have begun.
If the mass to light ratio is assumed to be constant, this means that the mass of the galaxy has not changed significantly during the said period. Furthermore, the retardation time, which has a significant effect on the galactic rotation curve, is a square root of the ratio of the galactic mass by the second derivative of the same mass. This means that if the retardation time (and retardation distance) have not changed considerably during the duration of observations, then the second derivative must also be approximately constant during that period. The M33 Triangulum Galaxy was probably discovered by the Italian astronomer Giovanni Battista Hodierna before 1654. In his work, De systemate orbis cometici; deque admirandis coeli caracteribus ("About the systematics of the cometary orbit, and about the admirable objects of the sky"), he listed it as a cloud-like nebulosity or obscuration and gave the cryptic description, "near the Triangle hinc inde". This is in reference to the constellation of Triangulum as a pair of triangles. The magnitude of the object matches M33, so it is most likely a reference to the Triangulum galaxy [20]. This observation took place 366 years ago; however, there was no measurement of the apparent or real luminosity of M33 and of course no measurement of Doppler shifts were available at that time. It was only Hubble that established the distance to this object and thus declared it an independent galaxy far away from the Milky Way [21], thus enabling the calculation of the real luminosity and mass of M33. Doppler shifts of spectral lines were obtained even later, thus enabling the deduction of rotation curves. In what follows, we shall assume an observation duration of: AsM(t) is approximately constant during the said duration, it follows that: Taking into account Equation (82) gives a lower bound to τ: This is not a strong bound as τ can be anything from 100,000 years to ten times the age of the universe or even larger. Similarly, the fact that the mass of the galaxy has not changed considerably during the observation duration indicates that: The upper bound on the mass derivative is thus aboutṀ max 1.6 × 10 30 kg/s, in which we have taken the baryonic mass of the M33 galaxy to be about 2 × 10 40 kg [17]. Using Equation (86), we arrive with an upper bound toṀ(0):Ṁ Notice, however, that without a handle on τ and |M(0)|, this does not yield any definite information.

Conclusions
The need to satisfy the Lorentz symmetry group prevents the weak field approximation of GR from allowing action at distance potentials and thus only retarded solutions are allowed. Retardation is manifested more strongly when large distances and large second derivatives are involved. It should be stressed that the current approach does not require that velocities, v are high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal-in other words, the ratio of v velocities in galaxies are 100 km s (see Figure 6), which makes this ratio 0.001 or smaller. However, one should consider the fact that every gravitational system, even if it is made of subluminal bodies, has a retardation distance, beyond which the retardation effect cannot be neglected. Every natural system, such as a star or a galaxy and even a galactic cluster, exchanges mass with its environment. For example, the sun loses mass through solar wind and galaxies accrete gas from the intergalactic medium. This means that all natural gravitational systems have a finite retardation distance. The question is thus quantitative: how large is the retardation distance? The change in the mass of the sun is quite small and thus the retardation distance of the solar system is huge, allowing us to neglect retardation effects within the solar system. However, for the M33 galaxy, the velocity curve indicates that the retardation effects cannot be neglected beyond a certain distance, which was calculated in Section 5 to be roughly R r = 4.54 kpc; similar analyses for other galaxies of different types have shown similar results [8,9]. We demonstrated, using a detailed model, in Section 6, that this does not require a high velocity of gas or stars in or out of the galaxy and is perfectly consistent with the current observational knowledge of galactic and extragalactic material content and dynamics.
We point out that, if the mass outside the galaxy is still abundant (or totally consumed),M 0 and retardation force should vanish. This was reported [19] for the galaxy NGC1052-DF2.
We note that the same terms in the gravitation equation that are responsible for the gravitational radiation recently discovered are also responsible for the rotation curves of galaxies. The expansion given in Equation (23), being a Taylor series expansion up to the second order, is only valid for limited radii: This means that current expansion is related to the near field case; this is acceptable since the extension of the rotation curve in galaxies is the same order of magnitude as the size of the galaxy itself. An opposite case in which the size of the object is much smaller than the distance to the observer will result in a different approximation to Equation (13), leading to the famous quadruple equation of gravitational radiation, as predicted by Einstein [22] and verified indirectly in 1993 by Russell A. Hulse and Joseph H. Taylor, Jr., for which they received the Nobel Prize in Physics. The discovery and observation of the Hulse-Taylor binary pulsar offered the first indirect evidence of the existence of gravitational waves [23]. On 11 February 2016, the LIGO and Virgo Scientific Collaboration announced that they had made the first direct observation of gravitational waves. The observation was made five months earlier, on 14 September 2015, using Advanced LIGO detectors. The gravitational waves originated from the merging of a binary black hole system [24]. Thus, the current paper involves a near-field application of gravitational radiation while previous works discuss far-field results.
We regret that direct measurement of the second temporal derivative of the galactic mass is not available. What is available is the remarkable fit between the retardation model and the galactic rotation curve, as can be seen in Figure 6, which constitutes indirect evidence of the galactic mass second derivative. The reader is reminded that competing theories like dark matter do not supply any observational evidence either. Despite the work of thousands of people and the investment of billions of dollars, there is still no evidence of dark matter. Occam's razor dictates that when two theories compete, the one that makes less assumptions has the upper hand. In the case of retardation theory, only baryonic matter and a large second temporal derivative of mass are assumed.
Problems that inflict dark matter theory, such as the cuspy halo problem (also known as the core-cusp problem) refers to a discrepancy between the inferred dark matter density profiles of low-mass galaxies and the density profiles predicted by cosmological N-body simulations. Nearly all simulations form dark matter halos, which have "cuspy" dark matter distributions, with density increasing steeply at small radii, while the rotation curves of most observed dwarf galaxies suggest that they have flat central dark matter density profiles ("cores"). This problem does not occur in the retardation model which denies the existence of dark matter. One cannot discuss flat or sharp profiles of dark matter if dark matter does not exist. The inherent problems with dark matter's dynamics further strengthen the claim of this work that dark matter does not exist and the rotation curve characteristics attributed to dark matter should be attributed to retardation.
How can one reach an erroneous conclusion regarding the existence of "dark matter"? By ignoring retardation effects and assuming that radial velocities are a result of some mysterious substance. We obtain for a spherically symmetric distribution [25]: where v c is the speed of a particle of unchanging radius r and M d (r) is the dark matter within the radius r.
Comparing Equations (105) and (28), we deduce that the "dark matter" mass can be calculated as follows: Then, since: it follows that: This is consistent with the observational data of [26] who concluded that the "dark matter" density decreases as r −1.3 for M33.
An additional method to explaining galactic rotation curves is to postulate that either the laws of dynamics or the laws of gravitation (GR) should be changed. This is the case in an approach championed by Milgrom (modifying the laws of gravity (MOND)-modified Newtonian dynamics) [27]. In this approach, the classical law of gravity is modified: In the above µ, is the interpolation function that should be one for a 0 a. Let us assume: if a 0 a, µ a a 0 . A particle revolving in a unchanging radius will have a centrifugal acceleration of a = v 2 r and thus: For the constant v at a far away distance, this expression is similar to the retardation force and thus: Milgrom found a 0 = 1.2 × 10 −10 ms −2 to be most fitting to the data. The velocity at 15.33 kpc from the center of the galaxy is 135,640 ms −1 . We thus obtain |M| 4.94 × 10 16 kgs −2 and a retardation time of: t r = M |M| 6.35 10 11 s (113) This amounts to a typical accumulation acceleration time scale of t r 20,129 years and a retardation distance of: R r = ct r 20, 129 light years, which seems reasonable according to our fitting estimates. Hence, despite the fact that modified Newtonian dynamics theory is phenomenological and is in contradiction to general relativity it can still serve as tool for estimating retardation theory quantities. We have thus shown that "dark matter" and "MOND" effects can be explained in the framework of standard GR as effects due to retardation, without assuming any exotic matter or modifications of the theory of gravity.
To conclude this section, we would like to mention the remarkable theory of conformal gravity put forward by Mannheim [29,30]. The current retardation approach leads to a Newtonian potential plus a linear potential. Such potential types can also be derived from the completely different theoretical considerations of conformal gravity. Indeed, on purely phenomenological grounds, such fits have essentially already been published in the literature. While those fits looked very good, they had to treat the coefficient of the linear potential as a variable that changed from galaxy to galaxy; this element is not in line with conformal gravity in which the coefficient of the linear potential is a new universal constant of nature. This can be explained much more easily in the framework of retardation theory, in which a variableM/M seems reasonable depending on the dynamical conditions of various galaxies. Indeed, M/M has a theoretical basis, while a pure phenomenological approach does not. As far as we understand the work of Mannheim [28][29][30], it 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 theory does not contradict conformal gravity and, in reality, both effects may exist, although Occam's razor forbids us to add new universal constants if the existing ones suffice to explain observations. Retardation theory's approach is minimalistic (in the sense that it satisfies the Occam's razor rule), does not affect observations that are beyond the near-field regime and, therefore, does not clash with GR theory and its observations (nor with Newtonian theory, as the retardation effect is negligible for small distances). We underline that a perfect fit to the rotation curve is achieved with a single parameter and we do not adjust the mass to light ratio in order to improve our fit as other authors do.
Retardation effects in electromagnetic theory were discussed in [31][32][33]. Finally the paper does not discuss dark matter in a cosmological context or in a gravitational lensing scenario; this is left for future works. Another important application is the study of galaxy clusters. This issue led Zwicky [1] to suggest dark matter in the first place.
Author Contributions: This paper has a single author, which has done all the work presented.
Funding: This research received no external funding.
University. Special thanks are due to Jiri Bicak for inviting me to present the theory at Charles University in Prague. I would like to thank Philip Mannheim and James Obrien for our discussions during the recent IARD meetings and for supplying some relevant data. I have benefited from discussions with a very long list of distinguished scientists and ask their forgiveness for not mentioning them all.

Conflicts of Interest:
The author declares no conflict of interest.