The Scale-Invariant Vacuum (SIV) Theory: A Possible Origin of Dark Matter and Dark Energy

The Scale Invariant Vacuum (SIV) theory rests on the basic hypothesis that the macroscopic empty space is scale invariant. This hypothesis is applied in the context of the Integrable Weyl Geometry, where it leads to considerable simplifications in the scale covariant cosmological equations. After an initial explosion and a phase of braking, the cosmological models show a continuous acceleration of the expansion. Several observational tests of the SIV cosmology are performed: on the relation between $H_0$ and the age of the Universe, on the $m-z$ diagram for SNIa data and its extension to $z=7$ with quasars and GRBs, and on the $H(z)$ vs. $z$ relation. All comparisons show a very good agreement between SIV predictions and observations. Predictions for the future observations of the redshift drifts are also given. In the weak field approximation, the equation of motion contains, in addition to the classical Newtonian term, an acceleration term (usually very small) depending on the velocity. The two-body problem is studied, showing a slow expansion of the classical conics. The new equation has been applied to clusters of galaxies, to rotating galaxies (some proximities with Modifies Newtonian Dynamics, MOND, are noticed), to the velocity dispersion vs. the age of the stars in the Milky Way, and to the growth of the density fluctuations in the Universe. We point out the similarity of the mechanical effects of the SIV hypothesis in cosmology and in the Newtonian approximation. In both cases, it results in an additional acceleration in the direction of motions. In cosmology, these effects are currently interpreted in terms of the dark energy hypothesis, while in the Newtonian approximation they are accounted for in terms of the dark matter (DM) hypothesis. These hypotheses appear no longer necessary in the SIV context.


Introduction
This work pertains to the exploration of new ways to try to understand the dark components of the Universe.We present the basic physical assumptions of the Scale-Invariant Vacuum (SIV) theory, the cosmological properties with new results, and several tests in the weak field approximation.The results differ from those of General Relativity (GR) by a small additional acceleration in the sense of motion.The cosmological solutions are close to the Concordance Model (ΛCDM), but with some moderate time-dependent differences.
Among the numerous theoretical attempts to explain the dark energy and dark matter in the Universe, the scale-invariant vacuum (SIV) theory is likely the one based on the most "economic" hypothesis, i.e., the assumption that the macroscopic empty space is invariant to scale transformations.We may call it the most economic one because this hypothesis is first a simple minimal extension of Einstein's GR, compared to the many ones proposed to account for the dark components of the Universe, and, as further discussed below, this hypothesis is already satisfied in electrodynamics, as well in General Relativity in the absence of cosmological constant.Over the last few decades, an impressive variety of new models based on exotic particles and fields have been forwarded in order to account for the assumed dark components of our Universe.Explorations in high energy particle physics have been pursued without much advance in solving the problem; see recent reviews by Bertone & Hooper [1] and de Swart et al. [2].Simultaneously, theories of modified gravity have been developed, in particular the theory of Modified Newtonian Dynamics (MOND) developed by Milgrom [3][4][5], which well accounts for some properties of galactic dynamics; see Famaey & McGaugh [6], Kroupa [7], Kroupa [8] and McGaugh et al. [9].
As pointed out by Dirac [10], "It appears as one of the fundamental principles in Nature that the equations expressing basic laws should be invariant under the widest possible group of transformations".For a long time, there has been a particular concern about a possible scale or gauge invariance in the theory of gravitation; see, for example, Weyl [11], Dirac [10] and Canuto et al. [12].
Scale invariance means that the basic equations do not change upon a transformation of the line element of the form, ds ′ = λ(t) ds , (1) where λ(t) is the scale factor, considered here not to depend on space for reason of homogeneity and isotropy.We know that a general scale invariance of the physical laws is prevented by the presence of matter, which defines scales of mass, time, and length [13].However, empty space at large scales could have the property of scale invariance, since by definition there is nothing to define a scale in it.As a matter of fact, the scale invariance of the physical laws in empty space is a fundamental and well known property.Maxwell's equations are scale invariant in empty space.Although it is not so frequently mentioned, it is also true for General Relativity (GR) if the Einstein cosmological constant Λ E is equal to zero.However, the equations of GR in empty space are no longer scale invariant if the cosmological constant is different from zero.According to Bondi [14], this was one of the reasons for Einstein's disenchantment with Λ E .The inclusion of scale invariance, in addition to the general covariance, implies to move from the Riemann Geometry to Weyl's Geometry, which accounts for the possibility of gauge transformations.Weyl's Geometry is based on the co-tensor calculus rather than tensor calculus, since the mathematical objects such as derivatives of scalars, vectors and tensors, Christoffel's symbols, etc. must have the property of gauge invariance in addition to the usual general covariance.Moreover, we use here a particular form of Weyl's Geometry, the Integrable Weyl's Geometry, see Appendix H.In this geometry, a general scale-invariant field equation and a geodesic equation have been obtained first by Dirac [10] and Canuto et al. [12].Noticeably, in this context, it is possible to have scale invariant field and geodesic equations in the presence of a non-zero cosmological constant.
At this stage, we do not know whether, in her wisdom, Nature presents the properties of gauge invariance simultaneously with a non-zero cosmological constant.However, our role as scientists is to explore new paths, particularly in the present context where 95% of the matter-energy in the Universe is in an unknown form.The observations are the judge; they will tell whether this is true or wrong.In this context, there already are some encouraging first signs as shown below.
A limited form of scale invariance is also present in MOND, which is remarkably successful in accounting for the high observed rotation velocities of spiral galaxies.MOND implies a dilatation invariance where lengths and time are multiplied by a constant λ-factor, (r ′ , t ′ ) → (λr, λt) , (2) where λ is independent from time.Expressions (1) and ( 2) are in fact similar for comparisons concerning the present epoch; however, the results are not necessarily the same, when cumulative effects over the ages intervene.The scale-invariant vacuum (SIV) theory presents several promising features.The cosmological equations predict an acceleration of the cosmic expansion [15], and the test analyses performed are in agreement with observations.The growth of the density fluctuations appears fast enough during the matter era so that it is not necessary to call for dark matter [16].Other tests of the SIV model are performed (related to the dynamics in galaxies based on the weak field approximation) concerning the mass excess in clusters of galaxies, the growth with age of the "vertical" velocity dispersion in the Milky Way [17], and the rotation velocities of galaxies [18].
Section 2 summarizes the basic properties of the SIV theory and gives the cotensorial form of the equations.Section 3 examines the cosmological equations, their results, and properties.The relation between H 0 and the age of the Universe is studied in Section 4. Three other basic cosmological tests are analyzed in Section 5: the distances, the magnituderedshift m − z or Hubble diagram is studied with the SNIa results of the joint analysis of the SDSS-II and SNLS observations by Betoule et al. [19], and with a recent quasar and GRB data by Lusso et al. [20], the past expansion rates H(z) analyzed by de Jesus et al. [21] and finally a study of the predictions of the redshift drifts.Section 6 summarizes some basic properties and tests of the weak field approximation (Newton like), including in Section 6.3 results on the growth of density fluctuations in the early Universe.Section 7 gives the conclusions and further perspectives.

Integrable Weyl Geometry, Cotensor Expressions, and the Scale Invariant Field Equations
The SIV theory is based on the Integrable Weyl Geometry (IWG), some properties of which are given in Appendix H.In addition to a metric form ds 2 = g µν dx µ dx ν , this geometry is endowed with a scalar field Φ: where λ is the scale factor, expressing a transformation of the line element as in Equation ( 1).
The quantity κ ν is the metrical connexion1 , it is a fundamental parameter of the IWG, where the lengths can undergo gauge changes of the form, The scale factor λ could in principle depend on the 4-coordinates x ν , but, because of homogeneity and isotropy, we consider it to be a function of time only.The quantities with a prime refer to Riemann Geometry of GR, while the quantities without a prime refer to the IWG.The integrable Weyl's space is conformally equivalent to a Riemann space (defined by g ′ µν ) via the λ-mapping that represents gauge re-scaling: In the space with a prime (GR), the metrical connexion κ ′ ν = 0.Such conformal mappings have been related to space-time deformations and have been studied in connection with a large class of Extended Theories of Gravity [22,23].According to Capozziello & Stornaiolo [22], the space-time deformations act like a force that deviates the test particles from the unperturbed motions.These authors also express the corresponding geodesics and general field equations that are consistent with those by Canuto et al. [12] and Bouvier & Maeder [24] as shown below.In relation to other works in the field, we may mention that Meierovich [25] is using a longitudinal vector field (cf. the derivative in Equation ( 3) ) to represent in the framework of GR the effects of the dark components of the Universe.The acceleration is a quadratic function of this vector field there, while in our SIV context it is a linear function.
The cotensor analysis has been developed [10][11][12]26], with some further complements about the derivation of the equation of geodesics from an action principle, the study of the notion of parallel displacement, of isometries, and Killing vectors [24].Standard scalars, vectors, or tensors that transform by following the rule: upon a gauge change (4) consistent with the conformal transformation (5), respectively define co-scalars, co-vectors, or co-tensors of power n = Π(Y ν µ ).If n = 0, one has an in-scalar, in-vector or in-tensor, such objects are invariant upon a scale transformation.In general, the derivative of a scale-covariant object (co-tensor, co-vector or co-scalar) is not scale covariant; furthermore, the derivative of a scale-invariant object is not scale-invariant either.However, one can write a modified scale-covariant derivative of an object, which has the properties of scale covariance, respectively of scale invariance when n = 0.Such scale covariant derivatives are called "co-covariant derivatives" by Dirac [10].The co-covariant derivatives have the same power as the original object.For example, the co-covariant derivative (marked by an asterisk * ) of a coscalar S of power n = Π(S) is The scale factor λ is considered to be a co-scalar of power Π(λ) = −1 by imposing the requirement λ * µ = λ ,µ − Π(λ)λκ µ ≡ 0 ⇔ Equation (3), which assures the commutativity of the scaling ( 6) and the co-covariant derivative operation.Furthermore, κ µ is then a zero power Π(κ µ ) = 0 gauge vector (κ µ → κ ′ ν = κ ν + ∂ ν ln λ) assuring that the co-covariant derivatives preserve the power of the object they are applied on.For a covector of power n = Π(A), the co-covariant derivatives are: with Here, * Γ α µν is a modified Christoffel symbol, while Γ α µν is the usual Christoffel symbol.The first derivatives of a co-tensor have the following expressions: Second and higher co-covariant derivatives of scalars, vectors, and tensors can also be expressed.Pursuing the developments, a corresponding Riemann-Christoffel tensor with all symmetry properties can be obtained; see, for example, Dirac [10].The contracted Riemann-Christoffel tensor appears to be an intensor; it writes where R ′ν µ is the usual expression.The sign ";" indicates a covariant derivative with respect to the mentioned coordinate.The total curvature R in the scale-invariant context is : where R ′ is the total curvature in a standard Riemann geometry.We note that, in all expressions, the additional terms depend on κ ν .For more details on the cotensor calculus, the interested reader may read Chapter VII of "The Mathematical Theory of Relativity" by Eddington [26], as well as Dirac [10].We also point out that Canuto et al. [12] provided a short summary of the cotensor analysis.The above expressions, developed in the IWG, consistently lead to a general scale invariant field equation [12] that contains additional terms depending on κ ν , where G is the gravitational constant, taken here as a true constant.The cosmological term in Einstein's GR-equation is Λ E g ′ µν , which according to expression (5) transforms like Here, Λ = λ 2 Λ E in the SIV context; it is a coscalar of power Π(Λ) = −2.Thus, the field equation in the above form (15) consistently is a scale invariant expression.
The second member of the RHS of ( 15) is scale invariant, and so is the first one.This means that the energy-momentum tensor behaves like an in-scalar, if G is considered to be in-scalar as well, This expression has some implications on the behavior of the relevant densities.The tensor T µν being scale invariant, one may write The velocities u ′µ and u ′ µ transform as follows: and Thus, the energy-momentum tensor transforms by This implies the following scaling of p and ϱ [12], Thus, pressure and density are not scale invariant but are coscalars of power Π(ρ) = −2, like the term Λ = λ 2 Λ E .However, to avoid any ambiguity, we keep all expressions with Λ E , the Einstein cosmological constant.

The General Field Equation from an Action Principle and the Lagrangian
The generalization of the Einstein GR-field equations to (15) and the relationships of the Riemann curvature and Ricci tensor to their counterparts within the IWG could be derived via the substitution of (5) directly into the corresponding tensors or via variation of the action: Here, we use an asterisk in front of R to emphasize that ( * R) is a IWG co-scalar of non-zero power and therefore some additional objects are needed to make the whole action in-scalar object.The in-scalar properties of this action are seen easily by re-calling that λ is a co-scalar of power Π(λ) = −1, while ( * R) is of power Π( * R) = −2 since it is a contraction of the in-tensor R µν using g µν , Π(g µν ) = 2, and g = det(g µν ) is therefore a co-scalar of power Π(g) = 8 [12].The corresponding variational equations are the same as those derived by substitution of (5) directly into the corresponding tensors.However, the variational approach allows for an extension of the Lagrangian to include additional terms in λ and its co-tensor derivatives ( * µ).That is, a term c 1 λ * µ λ * µ and a term c 2 λ 4 .The c 1 term is irrelevant because λ * µ = λ ,µ − Π(λ)λκ µ = 0 due to the definition of the metrical connexion κ µ = −∂ µ ln(λ).The quartic term then is related to the cosmological constant upon the substitution: Λ ∼ c 2 λ 2 .It has been shown by Canuto et al. [12] that the equations corresponding to the variation with respect to the scale factor λ do not produce new field equations but can be derived via trace of the corresponding Einstein field equations.In our case, (15) results in: The requirement of scale-invariance of the action leads to λT µ µ = λGT µ µ = δL δλ (see Canuto et al. [12]).Thus, in the absence of matter (L = 0 ⇒ T µ µ = 0) and in the Einstein gauge (λ = 1), the cosmological constant is related to the Ricci scalar R ′ = 4Λ E ; however, in a flat IWG vacuum (R ′ = 0), one then obtains an equation about the scaling factor: 6 In general, if T µ µ is not zero, it could contribute to an effective cosmological constant.In what follows, a specific gauge is chosen which allows for utilization of the flat IWG vacuum as a foundation for the scale invariant vacuum theory.

Fixing the Gauge
As it stands, the field Equation ( 15) is undetermined due to the gauge symmetry of the equations.The same problem appears in General Relativity, where the under-determinacy of GR is resolved by the choice of coordinate conditions.Here, one needs to impose some gauging conditions to define the scale factor λ. The key hypothesis is made that the properties of the empty space, at macroscopic scales, are scale invariant as in Maeder [15].This choice is justified since the usual equation of state for the vacuum p vac = −ϱ vac is precisely the relationship permitting ϱ vac to remain constant for an adiabatic expansion or contraction [27].Under the above key hypothesis, one is left with the following condition for empty "vacuum" spacetime deduced from the field Equation (15): The geometrical terms R ′ µν and R ′ of the field equation have disappeared from (15), since the de Sitter metric for an empty space endowed with a cosmological constant is conformal to the Minkowski metric where R ′ µν = 0 = R ′ [28].In particular, if 3λ −2 /(Λ E τ 2 ) = 1, where the time parameter τ is discussed in the 4th remark after Equation (27) and in more detail in (53) and in [15].This condition is noticeably satisfied by the solution of ( 23), as shown below, if λ is a function of time.
Since λ may only be a function of time t, only the zero component of κ µ is nonvanishing.Thus, the coefficient of metrical connexion becomes We use the Minkowski metric with a positive time-component signature, having checked that it is compatible with the scale invariant field equation in empty space [15].However, the results that follow are the equivalent if one is to use mostly positive metric signature.The 0 and the 1, 2, 3 components of what remains from the field Equation (23) become respectively: Notice that (25) and (26) correspond to the two non-trivial coordinate equations from (23); in this respect, the trace of (23), given also by (22), can be obtained by multiplying the spatial coordinate Equation ( 26) by (−3) and, when added to (25), one recovers (22) for the choice of gauge (24) .From (3), one has κ 0 = − λ/λ (using c = 1 units), while (25) and (26) lead to κ0 = −κ 2 0 along with: These expressions, which relate Λ E and the scale factor λ(t), result from our key hypothesis and some important remarks can be made about them:

•
Since Λ E represents the energy of the vacuum, it also establishes a relation between the energy density of the vacuum and the scale factor λ and its possible variations.
• As noted by Durrer and Maartens [29], "Instead of a cosmological constant, one may also introduce a scalar field or some other contribution to the energy-momentum tensor..." , this is just what is happening here.The field Φ as defined in Equation ( 3) is therefore a scalar field associated with the vacuum properties at macroscopic scales in the Universe.
• In GR, Λ E and therefore the properties of the empty space associated with it are considered to not depend on the matter content of the Universe.This is evident upon considering the Einstein gauge (λ = 1) in (21) along with the fact that radiation does not contribute to the trace T µ µ and the matter contribution vanishes since ϱ m ≃ a −3 → 0 in the limit of ever increasing scale factor a. Thus, in standard GR, one is forced to accept 4Λ E = R ′ .We adopt the same assumption here.This means that the above expressions of Λ E (27) are valid whatever the matter content of the Universe.In particular, if the IWG with λ satisfying ( 27) is adopted, then the corresponding equations ( 15) would imply local no-cosmological constant Einstein GR equations for the matter fields (see (30) and ( 31)).

•
The solution of the above two equations ( 27) is of the form: λ = y follows from the first equation.Any value of the parameter b would satisfy the equations.However, the value of b has to be consistent with the cosmological equations for the expansion factor a(t) when applied to the considered case of an empty space.The results below and in Maeder [15] show that b = 0 has to be taken in this case.

•
The expression of λ = ct in current units (in calculations c = 1 is usually taken).The value λ 0 at the present time is then

Accelerating Cosmological Models
The main cosmological consequence of the scale invariant vacuum theory is that the models predict an acceleration of the expansion as a result of the scalar field resulting from the properties of the vacuum.

The Basic Equations of SIV Cosmology
The application of the FLWR metric, appropriate for an homogeneous and isotropic space, to the general field Equation ( 15) leads to some additional terms resulting from the κ-terms in the cosmological equations as obtained by Canuto et al. [12], These equations do not have known analytical solutions for the expansion factor a(t) and have not led to satisfactory cosmological models.However, if we accept "the postulate of GR that gravitation couples universally to all energy and momentum" [27], we may account for expressions (27) characterizing the properties of a homogeneous and isotropic empty space.Thus, with the two relations in (27), the equations obtained from (15) nicely simplify to [15] 8 πGϱ The combination of these two equations leads to Term k is the curvature parameter which takes values 0 or ±1; p and ϱ are the pressure and density in the IWG.Einstein cosmological constant has disappeared due to the account of Equations (27).These three equations differ from the classical ones, in each case, only by one additional term containing ȧ λ a λ , which depends on time t.If λ(t) is a constant, one gets the usual equations.Thus, at any fixed time, the effects that do not depend on time evolution are just those predicted by GR (for example, the gravitational shift in stellar spectral lines).Significant departures from GR may essentially appear in cosmological evolution over the ages.
The significance of the additional term appearing in (32), where we have λ/λ = − 1 t , is as follows.If there is an expansion with ȧ > 0, it represents an acceleration in the direction of the expansion and opposed to gravitation.On the contrary, if there is a contraction, it is an acceleration in the direction of the contraction, thus adding to gravitation.Equations ( 30) to (32) are fundamentally different from those of the ΛCDM models: a variable term replaces the cosmological constant Λ E .The new term implies an acceleration (variable with time) of the expansion or contraction.In this connection, we note that the need to have a time-dependent term has been emphasized by several authors, Sahni et al. [30], Solá et al. [31], Ding et al. [32], and Solá et al. [33].

Density, Geometry Parameters, and Conservation Laws
Let us examine some of the general properties of the above equations in the matter dominated era with p = 0. Equation (30) divided by H = ȧ a can be written, With the following definitions, considered at time t 0 , we get The term Ω λ is associated with the scale factor λ(t).Let us examine how this acceleration term behaves over the ages.Based on the model results of Equation (52), where, in the case of a flat model, we have H(t) = 2 t 2 /(t 3 − Ω m ); this implies the following behavior of Ω λ : We see that relation ( 35) is valid at the present time, while at other epochs it becomes of the form (36), the two forms being evidently compatible at t 0 = 1.Thus, as time is going on, Ω λ is increasing, while the mass density term is decreasing.This shows that the acceleration term is gaining importance over the ages.Conversely in the early stages, the role of Ω λ was smaller compared to that of Ω m -notice that the first expression in Equation (52) implies Ω λ = 0 at t in when the scale factor takes the value a(t in ) = 0.This is also confirmed by the SIV solutions for the relativistic era [34], which shows a braking of the expansion with a(t) ∼ t 1/2 in the radiative era.
Here, ϱ c is the usual value of the present critical density in Friedman models, and so does Ω m .(We make this choice in order to not introduce new confusing definitions; note that another definition of the critical density could be introduced in the SIV cosmology [15]).While in Friedman's models there is only one model density corresponding to k = 0, here, for k = 0, there is a variety of models with different Ω m and Ω λ .For all models, whatever the k-value, the density parameter Ω m remains smaller than 1.For k = −1 and k = 0, this is clear since 2 /(t 0 H 0 ) is always positive for an expanding Universe.For k = +1 (negative Ω k ), numerical models confirm that Ω k + Ω λ is always positive so that Ω m < 1.Thus, scale invariant models only exist for Ω m < 1.In other words, scale invariance is limited to densities lower than critical.That is, the presence of matter in the Universe with a density higher than the critical one is killing the scale invariance.
Let us examine the geometry parameter q 0 = − ä0 a 0 ȧ2 0 .Equation (31), in the absence of pressure and divided by H 2 0 , becomes For flat models with k = 0, we get This establishes a relation between the deceleration parameter q 0 and the matter content for the scale invariant cosmology.Using (35), we also have which provides a very simple relation between basic parameters.For Ω m = 0.30, 0.20, or 0.10, we would get q 0 = −0.20 − 0.30 or −0.40.The above basic relations are different from those of the ΛCDM, which is expected since the basic Equations ( 30)-( 32) are different.Let us recall that, in the ΛCDM model with k = 0, one has For Ω Λ = 0.70, 0.80 or 0.90 and Ω m = 0.30, 0.20 or 0.10, one gets q 0 = −0.55,-0.70, or -0.85, which shows that the second derivative of the expansion factor a(t) behaves differently with the density in the two kinds of models.
Conservation laws and invariances are intimately connected, thus the scale invariance of the empty space must influence the conservation laws.Let us use the above cosmological equations for that.We first rewrite (30) as follows and take its derivative, Equations ( 30), (32), and ( 27) lead to welcome simplifications, This can also be written and 3 We may also express it in a form similar to the usual conservation law, These last two equations give the law of conservation of mass-energy in the scale invariant cosmology.For a constant λ, we evidently recognize the usual conservation law.We now write the equation of state in the general form, p = w ϱ, with c 2 =1, where w is taken here as a constant.The equation of conservation (44) becomes 3 da a + dϱ ϱ + dλ λ + 3 w da a + 3 w dλ λ = 0, with the following simple integral, Notice that the integer powers of a and λ correspond to familiar values in the literature: such as radiation, matter, and vacuum energy w = 1/3, 0, −1.The case w = −1 is the case of a cosmological constant ϱ Λ ∝ Λ E λ 2 .Notice that such scaling of the cosmological constant is consistent with scale-invariant stress energy tensor (17) of an ideal fluid discussed earlier (19); that is, a scale-invariant vacuum energy.For radiation, w = 1/3, the corresponding trace of the stress-energy tensor is zero therefore, it has no contribution to a potential effective cosmological constant via (21).For w = 0, this is the case of ordinary matter of density ϱ m without pressure, which means that the inertial and gravitational properties (respecting the Equivalence principle) within a covolume should slowly increase over the ages with M = M ′ /λ.However, interestingly enough, this implies that the gravitational potential ϕ = G M/r of a given object is an inscalar, meaning that it stays constant with time, since r also behaves like 1/λ.A change of the inertial and gravitational mass is not a new fact; it is well known in Special Relativity, where the effective masses change as a function of their velocity.In the standard model of particle physics, the constant masses of elementary particles originate from the interaction of the Higgs field [35,36] in the vacuum with originally massless particles.In addition, in the ΛCDM, the resulting acceleration of a gravitational system does not keep the total energy of the system unchanged [37].The assumption of scale invariance of the vacuum (at large scales) makes the inertial and gravitational masses to slowly slip over the ages; however, it is by a rather limited amount in realistic cosmological models.For example, in a flat model with Ω m = 0.30, λ varies only from 1.4938, at the origin (the Big-Bang where a = 0) to 1 at present.The lower the density, the bigger the variations.

Results of Cosmological Models
We study here the SIV models for the case of ordinary matter with w = 0, for a flat Universe with k = 0 as supported by the Planck Collaboration [38].The numerical solutions for k = 0, ±1 were obtained by Maeder [15].From (30), we have the basic equation in the case k = 0 considered here, Time t is expressed in units of t 0 = 1, at which we also assume a 0 = 1.The origin, the Big-Bang, occurs when a(t) = 0 at an initial time t in (see Equation ( 52)).The constant C contains the scale factor λ. With . However, as we may see below, the exact value of λ 0 has no influence on the solution of the cosmological models.Indeed, we notice that, if we have a solution a vs. t, then (x a) vs. (x t) is also a solution.This demonstrates an overall scale invariance of (48).The above equation implies at the present time t 0 , where we will take the sign "+" since H 0 is positive.Equation ( 30) becomes (33) with Ω m = ϱ/ϱ c and gives: and with (49) For Ω m ranging from 0 to 1, C varies between 0 to ∞, according to the above relation.The reason is that, when Ω m tends towards 1, the solution a(t) tends towards a singularity: the curve a(t) vs. t is a straight vertical line tending to infinity, its derivative H 0 is infinite and the same for C according to Equations ( 49) and ( 50), which follows from the general Equation ( 30).An alternative way to see this behavior for C, from its definition in Equation (48), is to recall that the general solution of the Equations ( 27) has the form λ = λ 0 /(t − t in ) which diverges when t tends towards t in .The general trend of a(t) is well illustrated in Figure 1.The integration of the differential Equation ( 48) leads to a useful analytical relation given by de Jesus [39] for flat space (k = 0): Figure 1.Some solutions of a(t) for the models with k = 0 [15].The curves are labeled by the values of Ω m .The Einstein-de Sitter model (EdS) is indicated by a dotted line.The small circles on the curves show the transition point between braking (q > 0) and acceleration (q < 0).For Ω m = 0.80, this point is at a = 2.52.The blue curve corresponds to Ω m = 0.20.
One can see again that relation (50) between the constant C and Ω m does not depend on the exact form of λ 0 and so does the solution (51) for a(t).The initial time of the Big-Bang is given by a(t in ) = 0, and the Hubble expansion rate by H(t) = ȧ/a, For H 0 , this is evidently consistent with expression (50).All the above values are on a timescale where t 0 = 1.We have taken as a boundary condition a(t) = 0. Some models with a vector field adopt da/dt = 0 [25], which then leads to a regular initial evolution instead of a singular Big-Bang.If one would adopt such an initial condition, Equation (30) would lead to a static empty Universe for k = 0. Thus, it seems difficult to adopt a different initial condition than a(t) = 0. Achieving a regular initial evolution via da/dt = 0 may not be needed for k = 0 since the Big-Bang could be viewed as an apparent singularity within the SIV paradigm.This point has been discussed towards the end of the section "Fixing the gauge within the scale-invariant vacuum theory" in [16].Here, we recall briefly that upon the use of a logarithmic time t → τ = ln(t/t 0 ) the "now" is at τ = 0, the Big Bang is at −∞, upon cosmic time when t in = 0, and the spacetime is locally Lorentz invariant since the metric tensor is now Minkowski like.
Figure 1 shows the model results for various values of Ω m .After an initial phase of braking, there is an accelerated expansion, which goes on all the way.Only the empty model starts with a horizontal tangent; in this case, the effects of scale invariance are the largest.All other models with matter start explosively with very high values of H = ȧ/a and a positive value of q, indicating braking, before turning later to an accelerated expansion (at points marked by small circles in Figure 1).Figure 2 shows the comparison of the scale invariant and ΛCDM models both with k = 0.The curves of both kinds of models are similar with larger differences for lower density parameters.These results imply that for most observational tests the predictions of the scale invariant models are not far from those of the ΛCDM models.The ratio λ(t)/ √ 3 Λ E varies between 1/t in = 1/(t 0 Ω 1/3 m ) and 1 at the present time.For an empty space, this ratio varies enormously, between ∞ at the origin and 1 at present.As soon as matter is present, the range of the above factor falls dramatically.For Ω m = 0.30, λ(t)/ √ 3 Λ E varies only from 1.4938 to 1.0 between the Big-Bang and now.Thus, the presence of less than 1 H-atom per cubic meter is sufficient to drastically reduce the amplitude of the domain of λ-variations.For Ω m approaching 1, the effects of scale invariance disappear and the cosmological solutions tend towards those of GR.Indeed, we have seen that there is no scale invariant solution for Ω m > 1.Thus, in the line of the remark by Feynman [13] in the Introduction, we see that the presence of even very tiny amounts of matter in the Universe very rapidly tend to kill scale invariance.The point is that, for realistic values of Ω m , the effects appear to be not yet completely killed and may thus deserve the present investigation.The elliptic and hyperbolic scale invariant models have been discussed in Maeder [15].
On the whole, it is a satisfactory result that the scale invariant models, which are resting on a very simple hypothesis, consistently predict an accelerated expansion.We now turn to comparisons with observations, which are the judge of the models.

Theoretical Predictions
The relation between the Hubble constant H 0 and the estimated age τ 0 of the Universe is a fundamental test.For a given curvature parameter k, it only depends on Ω m .The test is simple in its principle, but not necessarily in practice due to the tension in the observational determinations of H 0 , as well as uncertainties in the age of the Universe.
Let us call τ the age in the usual scale of years, or seconds, while t is a time in the scale where the present age t 0 = 1, t in is the initial time in this timescale.We can define τ as where ∆ 0 is the present age of the Universe in the chosen units, typically ∆ 0 = 13.8 • 10 9 yr [40].The Hubble expansion rate in the t and τ scales are related by Thus, for k = 0, the Hubble function H(τ) becomes: with Concerning the units, let us note that a recession velocity of 70 km s −1 Mpc −1 corresponds to (1/∆ 0 ) = 2.2685 • 10 −18 s −1 , i.e., ∆ 0 = 4.408 • 10 17 s or 13.969 Gyr.This correspondence gives easy calculations of H 0 vs. age in their proper units.Figure 3

Comparisons with Observational Data
The Saga of the determinations H 0 and of the age of the Universe could fill several books.Let us be short and first consider H 0 , always expressed below in km s −1 Mpc −1 .
There is a bi-polarization around two values: a high value from the cosmic ladder formed mainly by Cepheids and Ia supernovae and a low value from the CMB data obtained by the Planck collaboration.On the high side, Frieman et al. [40] gave a value H 0 = 72 ± 5 (in the same units), Freedman & Madore [43] obtained H 0 = 73 ± 2 (random) ±4 (systematic), a value later supported by Riess et al. [44] from HST data for SNIa which lead to H 0 = 73 ± 1.75.Recently, with a new anchor based on LMC Cepheids for the SNIa ladder, Riess et al. [41] obtained a slightly larger value H 0 = 74.03± 1.42; this is the value we retain in Figure 3 for H 0 based on the distance ladder.We note that some higher values have been given recently.Tully et al. [45] obtained values of H 0 = 75(±2) or 76.2 from various distance estimators.From the mass distribution and time delay of a strongly lensed quasar, Shajib et al. [46] have given 74.2 (+2.7−3.0).From distant lensed galaxies, Chen et al. [47] found a value of 76.8 (±2.6).Let us turn to the lower values, based on CMB analysis and BAO oscillations.The Planck Collaboration [38] gives H 0 = 67.3± 1.2 within the six-parameter ΛCDM cosmology.This result is supported by further studies by the same team [42] obtaining H 0 = 67.4± 0.3.The combination of BAO and SN Ia data into an inverse distance ladder leads to a value H 0 = 67.3± 1.1 [48], rather similar results were obtained by L'Huillier & Shafieloo [49] with H 0 = 68.49± 1.53.
Let us now turn to the estimates of the age of the Universe, independent of the cosmological models.The age estimates from globular clusters range between 12 and 15 Gyr, with a general reference age of 13.8 Gy [40], which we consider here.However, some interesting claims for a possible higher age have been made.VandenBerg et al. [50] studied several old halo subgiants with accurate parallaxes and high quality spectra and obtained a value of 14.3 Gyr for HD 140283; adding 0.4 Gyr for the time interval between the Big-Bang and the first star formation would lead to an age of 14.7 Gyr.Such an age of about one Gyr greater than the standard reference value is also supported by Roukema et al. [51], on the basis on Galactic Bulge microlensed low-mass stars.
Figure 3 compares the various data and SIV models.Clearly, the reference age together with distance ladder H 0 support a value of Ω m ≃ 0.20.A higher age of, say, 14.7 Gyr, would push Ω m down to about 0.14.The above age and H 0 does not depend on any cosmological model; this is not the case for CMB and BAO data.Even if they do not explicitly rely on the ΛCDM models, these determinations assume that the oscillations responsible for the CMB peaks and BAO oscillations obey a gravitation law, which is different from what is assumed by the SIV theory, as given by the field and geodesic Equations ( 15) and (72).Thus, in the present context, a determination based on this data, despite their high accuracy, would unfortunately not be a theoretically consistent process.
On the whole, the H 0 vs. age test supports a density parameter Ω m ≃ 0.20, maybe down to 0.14.It is not surprising to find a value lower than the current ΛCDM value of about 0.30, since, in the context of the scale invariant theory, the assumption of dark matter is not a necessary one.However, we note that the final word concerning the value of Ω m could only be given by models of the cosmological nucleosynthesis in the appropriate context.

Distances
Distances intervene in many cosmological tests.The distance of an object with coordinate r 1 (ϑ 1 = 0, φ 1 = 0) depends on the evolution of the expansion factor a(t).There are several definitions of distances: the proper motion distance d M , the angular diameter distance d A , and the luminosity distance d L , The relations between these three distances are model independent and d M is: where c is the speed of light in vacuum.One can express H as a function of z by using Equation ( 51), (t With the expressions H 0 = 2/(1 − Ω m ), a 0 /a = 1 + z, and with a 0 = 1, one gets: which can also be written as Then, the proper motion distance is: which can be integrated numerically.
Figure 4 shows the angular distance d A as a function of redshift for scale invariant models of various Ω m .Up to redshift z = 2, the SIV models give almost the same curve whatever Ω m .At large z, the differences remain relatively limited, compared to the ΛCDM models.This shows the need of accurate data for determining of values Ω m and for discriminating between the SIV and ΛCDM models [15].[15], with the scales t 0 = 1, a 0 = 1.

The Magnitude-Redshift Diagram of Distant Standard Candles
The Hubble or magnitude-redshift diagram has been a major tool for the discovery of the accelerated expansion [52,53].The flux f received from a distant standard candle, like a supernova of type Ia, with a luminosity L goes like f = L/(4 π d 2 L ) and the distance modulus m − M of a source of coordinate r 1 is m − M = const.+ 5 log(a 0 r 1 ) + 5 log(1 + z) , (64) where the proper motion distance a 0 r 1 is obtained from (63).We choose the constant so that m − M = 33.22 at redshift z = 0.01.Such calibration is based on the recent (m − M) vs. z diagram from the joint analysis of type Ia supernova observations obtained in the SDSS-II and SNLS collaboration given in Figure 8 in Betoule et al. [19].At this low redshift, the various models, whether ΛCDM or SIV, have a similar value of the above constant.Relations (m − M) vs. z are calculated for a few scale invariant models for different density parameters Ω m .Table 1 shows some results.We see that significantly different results between models of various Ω m only appear for redshift higher than z = 1. Figure 5 compares SNIa results of the above-mentioned joint analysis of the SDSS-II and SNLS data [19] with the scale invariant models.This analysis is based on a sample of 740 spectroscopically confirmed SNIa with high-quality light curves.The sample benefits from a considerable improvement in the calibration of the luminosity scale.The best fit ΛCDM model represented by a black line is obtained for Ω m = 0.295 ± 0.034.The flat SIV models with Ω m = 0.10 and 0.20 are indicated by red lines.We see, consistently with the plot of distances, that, for all considered redshifts, the differences between models of various Ω m are very small.A scale invariant model with Ω m in the range 0.10-0.20 encompasses the best fit ΛCDM model, thus suggesting a Ω m value of about 0.20, in good agreement with what we derived from the comparison of H 0 and the age of the Universe.Thus, the scale invariant models are consistent with the Hubble diagram for the available SNIa data, which show the acceleration of the expansion.The major difference, however, is that the SIV models do not imply some unknown form of dark energy.[19].This small difference is, however, of little importance in the comparisons below.The quasar sample contains 1598 objects with high quality UV and X-ray measurements for redshifts z in the range of 0.04 to 5.1.The distances are obtained from luminosity estimates based on a nonlinear relation between the UV and X-ray flux for the quasars [55].In Figure 6, the quasars have been binned in narrow ranges of redshifts; nevertheless, their scatter remains important.
For the GRBs, the sample considered by Lusso et al. consists of a collection of data of GRBs established by Demianski et al. [56].For the distance determinations, the relation between the peak photon energy and the isotropic equivalent radiated energy is applied according to the analysis by Demianski et al. [56], who have shown that this relation appears to not depend on redshifts.Quasars and especially GRBs offer a most remarkable extension of the Hubble diagram to extreme distances, far beyond the range presently covered by SNIa data.

Scale invariant
W m = 0.10 0.30 From Lusso et al. 2019 Figure 6 shows the comparison of SNIa, quasar, and GRB data with various curves as indicated in the plot.The two red lines show the SIV models of Table 1 for Ω m = 0.10 and 0.30.This last model lies very close to the ΛCDM model with Ω m = 0.30, which means that, for such density parameters, the Hubble diagram will have difficulties to discriminate between the ΛCDM and SIV models.This also implies that the SIV models lie rather close to the observations.Lusso et al. [20] claim that there is a significant difference between the ΛCDM models and the data, thus suggesting new physics.These authors use the cosmographic method to express the luminosity distance as a function of redshift.In this approach, which is model independent, the expansion factor a(t)/a(t 0 ) is developed in a Taylor series of terms, (65) Then, the luminosity distance is also developed in a Taylor series as a function of z up to z 5 .The coefficients of this development are themselves functions of the parameters H 0 , q 0 , j 0 , s 0 , l 0 , the last three being called the jerk, snap, and lerk.The coefficients of the development are all defined at z = 0. We see from Figure 6 that the third, fourth, and fifth order cosmographic developments, supposed in fact to best fit the data, appear to not well represent the data at high redshifts, say above z = 4.The authors themselves disregard two observational values between z = 7 and 9, with the argument that the Taylor series are problematic for such high redshifts; however, we do think that the problems occur already at lower values between z = 4 and 7.At these high redshifts, the ΛCDM model with Ω m = 0.30 better reproduce the observations than the Taylor series supposed to do it.The same is true for the SIV with Ω m in the range of 0.30 to 0.10.On the whole, we conclude the scale invariant models, with Ω m between 0.30 and 0.10, are consistent with high redshift GRB data, although a detailed statistical analysis could perhaps give more precise indications on the value of Ω m suggested by these observations.

The Past Expansion Rates H(z) versus Redshift
The expansion rate varies with time in a way depending on the models, thus the study of H(z) may provide new constraints on cosmological models.The function H(z) in the SIV theory, for k = 0, is given by the analytical relation (62) seen above.The first method to study H(z) vs. z was the cosmic chronometer.It is independent of the cosmological models; however, it depends on assumptions on star formation and evolution.This method has been developed and used by a number of authors [39,[57][58][59][60][61][62][63].The method is based on the simple relation: obtained from a 0 /a = 1 + z and the definition of H = ȧ/a.The ratio dz/dt is estimated from the luminosity and colour evolution of a sample of passive galaxies (with ideally no active star formation) of different redshifts which provides some age estimates based on evolution models of stellar populations.The assumption on the total absence of star formation (SF) remains a point in discussion, as well as the initial mass function.It is clear that, if there is some undetected SF in a galaxy, its age differences with respect to other galaxies will be biased and the same for H(z).
Another method is based on the Baryons Acoustic Oscillations (BAO).In brief, the oscillations result from the propagation of sound waves in the early Universe, where photons and baryons are strongly coupled.These waves have a characteristic scale of about 150 Mpc, the traces of which remain in the distribution of radiation and matter, observable in the CMB and in the clustering of matter respectively.The characteristic scale of BAO provides a standard ruler (at the time of recombination).After decoupling, the BAO participate in the expansion of the Universe, but their size can be known and thus the observations of their extension along and transverse to the line of sight provides the distance.The BAO were detected mainly by the two point correlation of quasars and to a smaller extent by the application of such methods to luminous red galaxies (LRG).Their distances, redshifts, and H(z) become accessible.The first detections were performed in 2005 by the SDSS and 2dFGRS collaborations [64,65]; they were followed by a number of further studies.The BAO at the highest redshifts of around z = 2.3 have been studied by Blomqvist et al. [66], Delubac et al. [67], Font-Ribera et al. [68].Among the recent statistical analyses of BAO data from the cosmic chronometer method, quasar and LRG distributions, we may quote de Jesus et al. [21] who found Ω m = 0.256 ± 0.014 in the context of the flat ΛCDM and Ryan et al. [69], [70] who find that the data mildly favor closed spatial hypersurfaces.A variety of different models have been compared to BAO oscillations data.For example, from combined data, Solá et al. [71] find some support for the Running Cosmic Vacuum model, which assumes that the vacuum energy density (thus Λ) and G can be a function of time.
Figure 7 shows the observed expansion rates H(z) as a function of z, compared to the best fit flat ΛCDM as given by de Jesus [39].Flat SIV models with different Ω m are also represented by red lines.For our internal consistency, we have chosen the same H 0 = 74 km s −1 Mpc −1 as in Figure 3.However, the low value of H 0 could also be chosen and would give a fine adjustment for a slightly lower value of Ω m .There we enter in the well known debate about the tension in the H 0 -determinations.Globally, a scale invariant model with Ω m ≃ 0.20 (or slightly lower) appears to be in agreement with the data from observations.The three accurate observations of BAO oscillations near z = 2.3 well fit both the ΛCDM and SIV models.We may remark that the values of H(z) are not completely independent of model assumptions.For example, the value of the standard ruler is fixed by the theory.In the SIV theory, the basic equations are different as we have seen it for the growth of density fluctuations [16] and the same situation may occur in the study of the oscillations in the initial plasma and their further evolution.The difference are likely small; nevertheless, this may introduce a possible small bias in this test.
On the whole, the results of the different tests in Figures 3, 5-7 all point towards a similar value of Ω m ≃ 0.20.The most important point is that these fundamental cosmological tests do not lead to a rejection of the SIV models, which have the major advantage to require neither an unknown source of dark energy, nor of dark matter as shown by the dynamical tests mentioned below.

The Redshift Drift in the SIV Theory
The redshift drift, i.e., the change of redshift z of a given distant object over the years, may become a major new cosmological test in the coming decades.The development of a new generation of spectrographs like ESPRESSO by Pepe et al. [72] may make such a determination possible.
The time shift of the redshift z of a distant galaxy in the observer restframe is according to McVittie [73] and Liske et al. [74], where H(z) is the expansion factor ȧ/a at redshift z = (a 0 /a) − 1.We note that Equation ( 67) is also verified in the scale invariant context, since the scale factor λ(t) applies to both time and space coordinates.Thus, the relation ∆t ∆t 0 = a a 0 used by McVittie in the demonstration is also valid here.We have seen in Section 5.1 that the expansion factor H(z) as a function of redshift can be written, as in Equation ( 62) Thus, the expression of the shift of redshift given by (67) becomes As such, this equation gives the shift of z in the timescale where the present time is t 0 = 1.To obtain it for example in years, one has to express the present Hubble constants in the appropriate units We obtain the final expression for the shift There, h 70 is the ratio of the chosen value of H 0 with respect to 70 km s −1 Mpc −1 .Figure 8 shows the curves of the shifts as a function of the redshift z predicted in the SIV cosmology for different values of Ω m .They are compared to some standard models of different values of Ω Λ and Ω m by Liske et al. [74].We notice the relative proximity of the standard and scale invariant curves in the case of Ω m = 0.30, which could make the separation of models difficult.However, the expected value of Ω m in the SIV cosmology is likely significantly smaller than in the ΛCDM models, being of the order of Ω m = 0.20 or even less; this makes the differences of the z-drifts between the two kinds of cosmological models possibly observable by very accurate observations in the future.

Properties and Tests in the Newton Like Approximation
In this section, we briefly mention the consequences of the SIV hypothesis in the weak field approximation and summarize some positive observational tests.Indeed, this shows the similarity of the effects in cosmology and in the Newtonian approximation.In both cases, the SIV hypothesis results in an additional acceleration in the direction of the motions.The reason why will be explicated after the basic Equation (74).In cosmology, these effects are currently interpreted in terms of dark energy, while, in the Newtonian approximation, they are accounted for in terms of the dark matter (DM) hypothesis.

The Basic Dynamical Equation and the Two-Body Problem
The SIV theory leads to important dynamical effects.The starting point is the geodesic equation, obtained by Dirac [10] and also demonstrated by an action principle by Bouvier & Maeder [24]; see also [17,75], with the velocity u µ = dx µ /ds-in the weak field non-relativistic approximation, with a metric given by g i i = −1 , for i = 1, 2, 3 and g 00 = 1 + (2Φ/c 2 ).For the potential, we take that of a central mass point M at a distance r.The Christoffel symbols become: With ds ≈ cdt and the velocities u i ≈ v i c = dx i cdt and u 0 ≈ 1, one obtains From Equation ( 19), the density is a coscalar of power Π(ρ) = −2, thus M r 3 = M ′ r ′3 λ 2 and with r ′ = λ r, we have M ′ = λM.The number of particles forming M is evidently not changing, but the inertial properties are changing, a situation also common in Special Relativity.However, the gravitational potential is a scale invariant quantity and is thus left unmodified, with κ(t) = 1/t.In case of expansion, there is an additional outwards acceleration opposed to the Newtonian gravity; in case of infall, the additional term adds its effect to the Newtonian attraction.This equation has the same meaning, but in a different context, as Equation (32) where we also see that the total acceleration is the sum of two terms: the Newtonian one and an additional one in the direction of the velocity.To explain it simply, let us say that the additional acceleration is the result of the negative pressure exerted by the scale invariant empty space.A positive pressure would inhibit the motions, a negative one is favoring them.The scale invariance of the equation of motion (geodesics) implies a time dependence of this negative pressure.The additional term is very small, the inverse of the age of 13.8 Gyr is 2.295 × 10 −18 s −1 or 70.85 km s −1 Mpc −1 , close to the current Hubble constant H 0 .
It is interesting to notice that the ratio x of the additional acceleration to Newtonian gravity behaves to first order like [17,18], where ξ = H(t)/κ(t) is of order 1 and ϱ c is the critical density in (34).For relatively dense systems with ϱ > 10ϱ c , the equality in the above expression is closely satisfied.The two-body problem has been studied by Maeder & Bouvier [75] and Maeder [17].The equation of motion (74) writes in the two polar coordinates r − r θ2 = − G M r 2 + κ(t) ṙ , and r θ + 2 ṙ θ = κ(t) r θ .
The second of these equations implies the following law of angular momentum conservation This equation with the first one in (76) leads to a modified form of Binet equation, the solution of which are conics, with a constant eccentricity and a secular variation of the semi-major axis (see above references).From this last relation, we see that the velocity of circular motions behaves like, Since the lengths increase with time, this expression implies that the circular velocity in the case of the two-body motion is a scale invariant, a fact which is consistent with the inscalar property of the potential found in Section 3.2.Indeed, the circular velocity being an inscalar will write Thus, the circular radius will behave like The orbital radius of circular motions, and the same for the semi-major axis A, not to be confused with the scale factor a, of the conics, experiences a secular increase so that the motions are in reality slow outwards spiraling motions.Due to the conservation law (77) of the angular momentum, the rotation period T is also changing the same way [75], If we consider a test mass around a central object of mass M, the additional term in (74) produces an acceleration proportional and along to the rotation velocity with no component along the Newtonian gravity.This acceleration leads to the above expansion of the orbit.At the same time, the above acceleration together with some slowing down due to the angular momentum properties leads to a conservation of the orbital velocity.Indeed, we verify in Equations ( 81) that the semi-major axis and the rotation period are varying in the the same way, consistently with a constant orbital velocity and also with the expression of the gravitational potential (Section 3.2).

Tests on Stellar Dynamics and Galaxies
We briefly summarize here the main points about the tests already performed, as applications of the modified Newton law (74).They concern clusters of galaxies, galaxy rotation, and the "vertical" velocity dispersion of stars in the Milky Way.

Clusters of Galaxies
Clusters of galaxies show mass to light ratios M/L which are currently much larger than the value expected for a typical stellar population (∼ 10 M ⊙ /L ⊙ ).As an example, for 600 clusters in the range of 10 14 to 10 15 M ⊙ , M/L ratios in the range of 300-500 M ⊙ /L ⊙ have been found [76].The problem within the SIV context has been discussed by Maeder [17], where further references and details are given.Writing (74) for a given galaxy in a cluster and multiplying by υ i leads to Integrating and taking the mean leads to with q ′ a numerical factor.The ratio of the last term to υ 2 is of the order of 2(R/t) | υ |.If the clusters are not too distant, t is of the order of the age of the Universe.We see that with Equation ( 83) a smaller mass is associated with a given velocity dispersion v 2 than in the standard theory where the last term is absent.The effects reach an order of a magnitude and more depending on the mean cluster density.
6.2.2.Galaxy Rotation and the RAR Spiral galaxies rotate faster in their outer layers than would be expected from Newton's theory for the mass present in them [77].This well-known problem, at the origin of the dark matter paradigm, is beautifully illustrated by the radial acceleration relation (RAR), which compares the dynamical gravity acceleration g obs = υ 2 /R with the matter gravity g bar = GM/R 2 in different points R of late, early and dwarf spheroidal galaxies, see in particular McGaugh et al. [9] and Lelli et al. [78].There, υ is the circular velocity at a distance R from the center, M is the mass inside radius R. In the SIV theory, a relation between g obs and g bar in the case of circular motion has recently been obtained by Maeder and Gueorguiev [18] in the case of galactic circular motions, This equation predicts two limiting cases: the Newtonian limit and a limit when g bar tends towards zero: (1) If k 2 ≪ g bar then g obs → g bar . ( The parameter k 2 (potentially a constant ) appears as a background limit of the dynamical gravity for a vanishing Newtonian gravity, independently of the radius considered.The background acceleration, k 2 is logically related to some cosmological properties of the Universe, in particular to the average matter density [18], Expression (84) may also be simplified for gravities g bar sufficiently larger than k 2 .One may write g obs → g bar + g bar k 2 .This may be compared to MOND, which also has two limits [3,4], the Newtonian one, and g → √ g N a 0 , for g N ≪ a 0 .There g N corresponds to g bar and g to g obs , a 0 = 1.2 • 10 −10 m s −2 , which is the so-called deep MOND limit.Figure 9 shows the comparison of MOND and SIV results.Over the range from the large elliptical galaxies to the low mass spirals down to the above-mentioned limit, the agreement between the two curves, as well as with the observations is excellent.The curves are different for very low gravities, with larger values of g obs for the SIV results, which show agreement with the dwarf spheroidals for the indicated value of k 2 .The MOND predictions and the observations of the dwarf spheroidals do not fit, a point also mentioned by Dutton et al. [79].Thus, the dwarf spheroidals, which are the objects where the amount of dark matter with respect to baryons is the highest, appear to have a critical role in constraining theories.The above results have been disputed by Banik & Kroupa [81] on the basis of the study by Hees et al. [82] on the motion of the Cassini spacecraft around Jupiter.The remark by Banik & Kroupa [81] is poorly founded and lacks a valid scientific basis.Indeed, the data reduction by Hees et al. does not apply to the SIV theory.The orbit of Saturn and the external field effects (EFE) have been calculated with the MOND equations and not with the appropriate SIV equation of motion (74).Specifically, for the determination of the orbital parameters of Saturn, the SIV theory should be used in the data reduction; otherwise, the test is inconsistent.In addition, Equation (84) for galaxies cannot be applied as such to the Solar System with the same constant k 2 .This term represents a background acceleration, which is evidently not the same for galaxies and the Solar system.
We have mentioned in Section 2.1 the introduction of a vector field by Meierovich [25] , which can represent the properties of dark matter well.This approach has similarities with ours and also results in a satisfactory representation of the rotation curves of spiral galaxies [83], and thus finding some justification for the MOND theory.According to Meierovich [84], the rotation curves also show some oscillations with a dependence of the plateau velocity on the mass of the central black hole.Such dependence is yet to be investigated within the SIV model and/or to be confirmed observationally.

The Vertical Dispersion of Stellar Velocities in the Galaxy
The velocity dispersion of stars in the Milky Way, in particular for the W-component (perpendicular to the Galactic plane), considerably increases with the age of the stars considered; see, for example, Seabroke & Gilmore [85].Continuous processes, such as spiral waves, collisions with giant molecular clouds, etc... are active in the disk plane and may effectively influence the stellar velocity distributions.However, the stars spend most of their lifetime out of the galactic plane and, in order to have some significant effect, the stars should receive some "heating" also away from the Galactic plane, as shown by these authors.This problem has a long history [86] and no consensus has been found on the origin of the effect [87].
A study of the oscillations of stars around the Galactic plane in the SIV context shows that an equation like (74) amplifies the oscillations.The velocity W(t in ) of a star born at time t in , when crossing the plane, is [17,88], where t 0 is the present time.The trend for the velocity dispersions follows that of the velocities.We take a value of 10 km s −1 for the present velocity dispersion σ W (cf. Figure 10).
As an example for a group of stars with a mean age of 10 Gyr, for an age of the universe of 13.8 Gyr, the velocity dispersion is estimated to be about 10 km s −1 × 13.8 3.8 = 36.3km s −1 .Figure 10 compares the corresponding model predictions obtained in this way (continuous red curve) with the data from Seabroke & Gilmore [85].We see that the theoretical curve well corresponds to the trend shown by the observations.[89].The continuous red curve shows the predictions of the scale invariant dynamics according to relation (87) for an age of the universe of 13.8 Gyr.The broken red curve accounts for the fact that the Galaxy formed about 400 Myr after the Big Bang.A vertical dispersion of 10 km s −1 is assumed at the present time.

The Growth of the Density Fluctuations
The growth of density fluctuations in the early Universe is a fundamental cosmological problem and one of the main reasons, along with the flat rotation curves of spiral galaxies, for the introduction of dark matter DM [29,90].At recombination (z ∼ 10 3 ), the observed fluctuations of the CMB are of the order of δT/T ≈ 10 −5 .At this time, the fluctuations of the baryonic density δρ/ρ are still of the same order, since, up to recombination, matter and radiation are closely coupled.From this time onward, the growth of the density fluctuations δ ≡ δρ/ρ is determined by gravitation and expansion [91].In the standard model, δ goes like t 2/3 ∝ a(t) ∝ (1 + z) −1 .Thus, starting with an amplitude of about 10 −5 , a growth by a factor of 10 3 would lead to fluctuations of about 10 −2 at the present epoch.This is by several orders of magnitude smaller than the large nonlinear structures presently observed in the Universe [92].
According to the standard model, during the radiative era, the DM perturbations (not subject to radiation pressure) would be growing, while the baryonic perturbations do not amplify.After recombination, the baryons fall into the potential wells of the previously assembled supposed dark matter, thus giving an extraordinary boost to the growth of the baryonic density fluctuations.This scenario of the cold dark matter (CDM) is a basis of the standard ΛCDM model.
Classically, the equations of Euler and Poisson, as well as the continuity equation, form the basis for the study of the evolution of the density fluctuations.The scale invariant form of the Euler equation becomes [16] The density ϱ is the baryon density.It contains an additional term proportional to the velocity, as in Equation (74).This term plays an important role in the enhancement of density fluctuations in the SIV theory.The corresponding continuity equation is: The second member can be viewed as a diffusion term.The potential being scale invariant, the Poisson equation ⃗ ∇ 2 Φ = △Φ = 4πGϱ stays the same.These three equations applied to a density perturbation δ = δϱ/ρ in the linear regime lead to the basic equation for the growth of density perturbations in the SIV theory [16], Here, H is the Hubble expansion rate at the time considered, n is related to the slope of the density profile of the perturbation, n = d ln δ/d ln x, where x is the distance to the perturbation .The NFW profile, which in the outer layers of a cluster is very close to the ϱ ∼ 1/r 2 profile [93] would have a value n = 2.The terms with κ = − λ/λ are absent in the standard models (there λ = 1) where one has the well-known equation, Figure 11 shows some results.The main result is that the growth of the density fluctuations δ is fast enough in the SIV cosmology to account for the formation of galaxies long enough before the present.There is no need for the dark matter hypothesis to favor structure formation.Another point we notice is that the growth of δ is faster for lower Ω m values.This result, surprising at first glance, is consistent with the fact that the scale invariant effects are larger at lower densities as shown by Equation (75).The results of Figure 11 are obtained for a value n = 2. Higher values such as n = 3 or 5 produce a faster growth of the density fluctuations [16].For example, for Ω m = 0.10 and n = 2, δ = 1 is reached at z = 10, at z ≈ 25 for n = 3, and at z ≈ 80 for n = 5, which is certainly too high.A recent observation of a gravitationally lensed galaxy shows signatures of star formation (OII lines) implying an active star formation at about z = 15 [94].Thus, also depending on the exact value of Ω m , a slighltly steeper concentration than that of the typical n = 2 model may lead to an early star formation in the Universe.This is also a positive feature of the SIV models.

Conclusions and Perspectives
The SIV hypothesis that the macroscopic empty space is scale invariant leads to an additional acceleration of motions, which in cosmology produces effects currently attributed to the dark energy, while, in the Newtonian-like mechanics, it gives effects attributed to the dark matter.At this stage, there is already a number of favorable tests of the SIV theory both in cosmology and in the study of astrophysical motions.These results do not constitute a proof of the theory.However, they are sufficiently promising to justify further investigations of the SIV theory and its properties.

Figure 2 .
Figure 2. Comparisons of the a(t) functions of the ΛCDM and scale invariant models for given values of Ω m [15].
illustrates with black lines the values of H 0 in km s −1 Mpc −1 vs. the density parameter Ω m as predicted by Equation (57) for three different values of the age of the Universe: ∆ 0 = 12.6, 13.8, and 15 Gyr.
. W m and age in SIV cosmology

Figure 3 .
Figure 3.The relation between the Hubble constant H 0 and Ω m predicted by the SIV theory (Equation (57)) for three different values of the age of the Universe.The result from the distance ladder with a new calibration of the LMC Cepheids of 74.03 (±1.42) km s −1 Mpc −1 by Riess et al. [41] is indicated by a thick red line, with its uncertainty.The thick blue line indicates the value of 67.4 (±0.5) km s −1 Mpc −1 obtained by the Planck Collaboration [42].

Figure 4 .
Figure 4. Comparison of the angular diameter distances for ΛCDM and SIV models of different density parameters Ω m[15], with the scales t 0 = 1, a 0 = 1.

Figure 5 .
Figure 5.The magnitude-redshift diagram by Betoule et al. [19] based on the joint analysis of SNIa observations in the SDSS-II and SNLS collaborations.The upper right corner of the plot is enlarged in an insert.The scale invariant models with Ω m = 0.10 and 0.30 are indicated by red lines.The black line is the best fit ΛCDM model for Ω m = 0.295.A new study of the Hubble diagram extended up to redshifts z ≃ 7 has been recently published by Lusso et al. [20].It is based on a sample of SNIa, quasars, and gamma-ray bursts (GRBs).The SNIa sample consists of 1048 supernovae (The Pantheon sample) of type Ia up to z = 2.26 established by Scolnic et al. [54] on the basis the PS1 Medium Deep Survey, together with HST samples and distances estimates from SDSS and SNL surveys.The reference of the magnitude scale used by Scolnic et al. appears from his plots to correspond to 33.17 at redshift z = 0.01, instead of 33.22 in the calibration by Betoule et al.[19].This small difference is, however, of little importance in the comparisons below.The quasar sample contains 1598 objects with high quality UV and X-ray measurements for redshifts z in the range of 0.04 to 5.1.The distances are obtained from luminosity estimates based on a nonlinear relation between the UV and X-ray flux for the quasars[55].In Figure6, the quasars have been binned in narrow ranges of redshifts; nevertheless, their scatter remains important.For the GRBs, the sample considered by Lusso et al. consists of a collection of data of GRBs established by Demianski et al.[56].For the distance determinations, the relation between the peak photon energy and the isotropic equivalent radiated energy is applied according to the analysis by Demianski et al.[56], who have shown that this relation appears to not depend on redshifts.Quasars and especially GRBs offer a most remarkable extension of the Hubble diagram to extreme distances, far beyond the range presently covered by SNIa data.

Figure 6 .
Figure 6.The Hubble diagram for SNIa, quasars (binned), and GRBs from the samples collected by Lusso et al. [20].The various models considered by Lusso et al. are indicated.The two red lines show the flat scale invariant models with k = 0 and Ω m = 0.10 and 0.30.

Figure 7 .
Figure 7.The observed expansion rates as a function of redshifts, compared to the flat ΛCDM Ω m = 0.256 as given by de Jesus et al.[21].The scale invariant models with various values of Ω m are also represented by red lines.For consistency within the SIV, and also because it is model independent, we take the value of H 0 based on the distance ladder: H 0 = 74 km s −1 Mpc −1 as in Figure3, while de Jesus et al.[21] have H 0 = 69.5 (±2.5) km s −1 Mpc −1 .

Figure 8 .
Figure 8.The drifts of redshifts dz/dt as a function of redshift in the scale invariant theory (red curves).The values of Ω m (usual definition) are indicated.The black broken lines give the results for some standard models of different couples (Ω Λ , Ω m ) by Liske et al. [74].

Figure 9 . 2 ]
Figure 9. Observed values of g obs and g bar for the 240 galaxies studied by Lelli et al. [78], forming the radial acceleration relation (RAR).The big green hexagons represent the binned data of the dwarf spheroidal galaxies.The blue curve shows the relation predicted by (84) with a value of k 2 = 10 −10.85 m s −2 .The red curve gives the MOND relation [80] g = g N [1 − e −(g N /a 0 ) 1 2 ] −1 with

Figure 10 .
Figure 10.The vertical dispersion σ W as a function of the age of the stellar populations.The blue triangles with the error bars result from the analysis by Seabroke & Gilmore [85] of the observations by Nordstrom et al.[89].The continuous red curve shows the predictions of the scale invariant dynamics according to relation(87) for an age of the universe of 13.8 Gyr.The broken red curve accounts for the fact that the Galaxy formed about 400 Myr after the Big Bang.A vertical dispersion of 10 km s −1 is assumed at the present time.

Figure 11 .
Figure 11.The growth of density fluctuations for different values of parameter Ω m for n = 2.The values δ = 1 are reached in the interval of redshifts (z + 1) = 4.0 to 29.3.From[16].

Table 1 .
Values of (m − M) as a function of redshift z for different density parameters Ω m