Fake dark matter from retarded distortions

We push ahead the idea developed in a previous work, that some fraction of the dark matter and the dark energy can be explained as a relativistic effect. The inhomogeneity matter generates gravitational distortions, which are general relativistically retarded. These combine in a magnification effect since the past matter density, which generated the distortion we feel now, is greater than the present one. The non negligible effect on the averaged expansion of the universe contributes both to the estimations of the dark matter and to the dark energy, so that the parameters of the Cosmological Standard Model need some corrections. In this second work we apply the previously developed framework to relativistic models of the universe. It results that one parameter remain free, so that more solutions are possible, as function of inhomogeneity. One of these fully explains the dark energy, but requires more dark matter than the Cosmological Standard Model (91%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$91\%$$\end{document} of the total matter). Another solution fully explains the dark matter, but requires more dark energy than the Cosmological Standard Model (15%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$15\%$$\end{document} more). A third noteworthy solution explains a consistent part of the dark matter (63%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$63\%$$\end{document} of the total matter) and also some of the dark energy (4%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4\%$$\end{document}).


Introduction
The evidences of dark matter are of many kinds. We can roughly divide the "dark matter phenomena" in two categories: global dark matter effects, which consist in unexpected values of some cosmological parameters (the deceleration parameter [1,2], the deuterium abundance, the power spectrum of CMB anisotropies, etc.), and local dark matter effects, appearing in observations of astronomical objects (the galaxies rotation curves [3], the virial of galaxy clusters [4], gravitational lensing [5], etc.). In both cases are exclusively gravitational phenomena, anomalies of gravitaa e-mail: fre@uninsubria.it (corresponding author) tion with respect to what we expect: some of them derive directly from the study of a gravitationally interacting system, as a galaxy, or from the space-time dynamics, as the deceleration parameter; for some others the derivation is indirect, like for the deuterium abundance, which, anyway, depends on "dark matter contributions", i.e. on gravitation effects that cannot be associated to a distribution of baryonic matter. Several followed hypotheses on dark matter assume that it could also interact weakly, at least, but such "weakly dark matter effects" have not been observed up to now [6,7].
The evidences of dark energy are restricted to more specific phenomena. There are only global dark energy effects, which consist in unexpected values of some cosmological parameters, due to a negative deceleration parameter [1,2]. Again, it is a purely gravitational phenomenon, consisting in an unexpected distortion of the space-time metric. Dark energy is often described as an exotic kind of energy having a negative pressure, or as a "vacuum energy"; these assume some non gravitational interaction. However, none of such new non gravitational interactions have been observed yet, and as for the dark matter, all our knowledge about the dark energy comes just from gravitational phenomena.
In the last 20 years, several lines of research have been opened that seek to study some essentially relativistic effects in cosmology. MOND theories are an attempt in this direction, but they actually didn't give rise to a good matching with data [8][9][10]. Rather, we will consider rather the attempts to obtain unusual effects from the usual General Relativity.
A general relativistic explanation of the unexpected distortion of the space-time tensor, even without the presence of real energy-matter, would provide the justification of at least a part of dark matter and/or dark energy effects. A point is that the formalism usually used to show these effects is not truly general relativistic. Especially for local effects, it is adopted the Newtonian approximation, while for global ones, one assumes a Friedmannian model for the expan-sion, with matter and curvature assumed to be homogeneous, which could be an oversimplification. By now, we will refer to the Newtonian approximation, for galaxies, and the homogeneous Friedmannian model, for the universe expansion, as "the standard calculations".
There are three principal lines of research that try to explain (at least a fraction of) dark matter and/or dark energy as a difference between the standard calculations and a more precise general relativistic model. One of them started form a deep analysis of the coordinate systems adopted in observational general relativity [11][12][13][14][15][16][17]. Recently it had some confirmations by observations on the rotation curves of our galaxy [18,19]. Another line of research explores the backreaction effects, which means the discrepancy between the spatial curvature due to an averaged quantity of matter, and the averaged spatial curvature due to an inhomogeneous quantity of matter [20,21]; it doesn't vanish in general since the Einstein Equations are not linear. This effects needs averaging on large dominions, and arise the standard deviation of the spatial curvature.
The third line considers two main ingredients: the inhomogeneity of the matter distribution, and the delay of the space-time distortions that it generates. In the universe, the matter is inhomogeneous and anisotropic, so that the Birkhoff Theorem cannot be applied in general; the distortions due to the inhomogeneity do not cancel out. The influence of the external objects on the dynamic of galaxies was recently confirmed in [22], which allows also a quantitative comparison. The Cosmological Principle states that, at large scale, the matter distribution is near to a homogeneous one, so that we must remember that globally the matter inhomogeneities are quite small. However, the second ingredient, the retarded distortions, provides a magnification effect, for which the resulting space-time distortion is not negligible even if its source is small. This magnification is due to the fact that the distortions we feel nowadays were generated by matter sources in the past, when the matter density was greater, because of universe expansion. The actual distortion is a superposition of all retarded distortions from all past times, which predictably have a singularity at the Big Bang time. However, there is also a decrease with the distance. It needs a precise calculation to see if it prevails the magnification or the reduction, and [23,24] highlighted that the magnification dominates.
These ideas were presented in nuce in the article [23], which however had several lacks. It doesn't work in a general relativistic context. The authors adopted a linearized gravity model on a Minkowskian background, where the simulated expansion of the universe is forced by hand, and an effective FLRW metric is imposed with a "compatibility condition". In [24] we firstly developed a truly general relativistic framework apt for this line of investigation. Since the matter inhomogeneity is small according to the Cosmological Principle, we considered it as a first order perturbation on a homogeneous FLRW background metric. We derived the linearized Einstein Equations, which returns the metric distortion, at the first order. Choosing the suitable gauge, these result to be wave equations, so that the metric distortion is retarded, traveling at the speed of light. This framework was then tested on a specific universe, which was not a realistic one, but a model chosen such that the linearized Einstein Equations have constant coefficients, with the only aim to make them easier to be solved. We obtained an explicit formula for the retarded gravitational potentials. Averaging the resulting metric, we got an evaluation of the global effects of such perturbative correction, which manifested the magnification effect we hoped, and was able to explain a non negligible fraction of the dark matter and the dark energy.
In the present article, we will push forward such treatment in a substantial way, by applying the method to a general realistic background universe, with any possible combination of matter and energy components. The main consequence is that the linearized Einstein Equations will not have constant coefficients, in general. From those, we will obtain an averaged metric, which can be compared with the metric assumed by the Cosmological Standard Model. We will obtain a set of consistency conditions, which reduce the acceptability of the possible universe. For acceptable universes, the gap between our averaged metric and the standard calculation provides some quantity of "fictitious" matter and dark energy. We will apply these formulas to the case of the real universe, which depends in general on three parameters: the real quantity of radiation, of matter, and of dark energy (where the last one can eventually be zero, if all the dark energy can be explained as fictitious). Substituting the most recent cosmological data, two parameters will be fixed. The last one will leave a range of possible solutions, some of which looks of particular interest.
The paper is structured as follows. In Sect. 2 we outline the perturbative method, defining the background universe and the perturbed universe, and showing the fictitious matter and dark energy from the gap between them. In Sect. 3 we recall the linearized Einstein Equations from [24]. In Sect. 4 we average the metric perturbations, ignoring temporarily their precise development, and we get a formula for the fictitious matter and dark energy. We reduce the PDEs for the evolution of the metric perturbations to ODEs, with an averaging procedure. In Sect. 5 we solve such ODEs with a single-component approximation, finding also some "Selfconsistency Conditions" for the components that can fill the universe. In Sect. 6 we apply all the previous formulas to our universe, getting an explicit model, which will result to leave a free parameter. The model is numerically solved in Sect. 7, finding a set of solutions depending on the free parameter.

Framework and notations
We will adopt the following terminology. With "radiation" we mean any component of the universe with a pressure p = wρ with w = 1 3 . We call "matter" any component with w = 0; and "dark energy" any one with w = −1. A component with w < 0 will be called "exotic".
Moreover, we call "total matter" the quantity of matter M0 , which is required in the Cosmological Standard Model in order to explain the observed deceleration parameter; analogous for the "total dark energy" 0 . We call "dark matter" DM0 the part of total matter that is not observed, unless indirectly via gravitational phenomena. The observed part is essentially the "baryonic matter" B M0 ∼ = M0 − DM0 .

Perturbative method
Let us consider a universe filled with any choice of components, each one characterized by its constant w. ρ = w ρ w is the true matter-energy density. Only the matter component ρ M (x; τ ) can be inhomogeneous. Let us call its homogeneous part 1 Let now consider a background universe, approximating our true universe at zero order, so assuming it filled with a perfectly homogeneousρ. We assume this universe to be spatially flat, in order to keep all calculations as simple as possible. For the same reason, we assume irrotational matter. Then, for a more realistic description of the universe, we perform a first order perturbative. The perturbation of the energy-matter, which in fact consists only on matter, is Notice that this inhomogeneity is always non negative. In a symmetric way, we could also definē in such a case,ρ would be non positive. It is a matter of convenience to fix the choice with always positive or negativẽ ρ.

Background evolution
We adopt the most minus signature and natural units, so that c = 1. The background metric is Any quantity has associated an "unperturbed" (or "averaged") versionQ, computed using the background metricḡ μν , and a "perturbed" (or "true") version Q computed using the real metric g μν . Its "perturbation" is the differencẽ Q := Q −Q, which we consider negligible beyond the first order.
τ is the conformal time for the unperturbed metric. The dot will always denote derivation with respect to the conformal time:Q := ∂ τ Q. a(τ ) is the unperturbed expansion parameter, so thatt = a(τ )dτ is the usual (unperturbed) time. H (τ ) :=ȧ a is the Hubble parameter for the unperturbed model. τ , a and H are written without the overline, with abuse of notation, for better readability. Their perturbed versions will be a and H. τ 0 is the actual instant for the unperturbed universe, s.t. a(τ 0 ) = 1. t 0 is the present for the true model, s.t. a(t 0 ) = 1. The "0" label means evaluation in the present time, both for unperturbed and perturbed quantities. Notice that in general t 0 = t (τ 0 ), again with a little abuse of notation. a(τ ) solves the Friedmann Equations with p(τ ) = w p w (τ ) = w wρ w (τ ). From the second Friedmann Equation we know We can describe the background components as (2.6) We define a(τ ) as a maximal solution of this ODE, with maximal domain (τ I ; τ F ), eventually unbounded. The initial condition for a is provided by the request that lim τ →τ I a(τ ) = 0. The radius of the visible universe results to be R(τ ) = τ − τ I ; notice that it is always infinite if τ I = −∞.

Comparison with the cosmological standard model
Averaging itρ on the whole space, we get ρ (τ ). The average of a certain S on the whole x ∈ R 3 is defined as the average on some increasing sequence of compact domains D ⊂ R 3 filling the whole space in the limit The exact process of averaging is described in [25]. Since we need the average on the whole space, we perform the averaging prescription called J (S, ρ, V 0 , A, B s , B s ) with the thickness B s tending to infinity, and with trivial weight ρ ≡ 1. The choice of B and B s are hence irrelevant, and we take V := A := τ . Indeed, we do not need a lightlike gradient for V , as it is prescribed in [25] for averaging cosmological observables: a set of measures suffers from the delay of information, due to the speed of light, but the quantities ρ or g μν we need to average are not measures, but rather inhomogeneous fields that we want to consider as they were homogeneous, at a certain instant of the time foliation. Moreover, the relativistic delay of information is already accomplished by the study of the retarded potentials, as in Sect. 3. For these reasons, we choose a time section for V .
The a 3 factor is simplified dividing by the volume J (1, ρ, V 0 , A, B s , B s ). Thus, for a first order calculation, this cosmological average is equal to the average of S − 3 S on a usual Euclidean metric. Notice that the second term vanishes whenever S is a first order quantity, as it will be throughout this article.
We can define the "inhomogeneous matter" as a part of the total matter component The Cosmological Standard Model measures the cosmic components via the observed deceleration parameter ⎧ ⎨ ⎩ w w0 = 1 but it assumes a homogeneous ρ. Since it is not our case, a will be obtained by just an adaptation of the true space-time metric to a FLRW one This provides a distortion of the expansion law, so that in general a(t) = a(t), q 0 = −ä H 2 0 , and w0 =¯ w0 . To fit the two conditions (2.8), two more parameters are needed. Interpreting the distortion as the unexpected presence of matter and dark energy, we will see an effect of "fictitious matter and dark energy". We evaluate them as F M0 , F 0 .
Remark 1 Since they come from a global evaluation and we used a first order approximation, these fictitious components will result to be proportional to the total perturbation I M0 . Thus, these global effects do not depend on the spatial distribution of the matter inhomogeneities, but only on their total amount.
The matter and the dark energy components used in the CSM have a "true" and a fictitious part The other components are all "true". The "true" parts must be proportional to the same components of the background universe with the exception of matter, for which we have to add again the inhomogeneous part Remember I M0 can be considered as positive or negative. In the second case, the homogeneous approximationρ is a rounding up, so that H M0 > T M0 . Some of the true matter must be the baryonic matter we know to exists. If there is still some part left, it is "true dark matter" T DM0 . It is some kind of matter that actually exists, which gravitational action is not just a relativistic effect, but that is not a directly observable matter, like primordial black holes, neutron stars, and so on and so forth.

Classifying possible results
We can apply this framework to a universe filled by any choice of components {w}.

Definition 1 We will call "selfconsistent" a choice for which
• all calculations return a finite result; • the linearized Einstein Equations give a unique solution; • the perturbative method is justified by small enough perturbations.
We can write the last condition as where the second condition means that the Cosmological Principle holds.

Definition 2
We will call "acceptable" a choice for which The second part states that all the baryonic matter we see is really existing, so it is included in the model.

Remark 2
Notice that the fictitious components can be negative, and such a case means that the dark matter and/or the dark energy is not explained at all, but rather its quantity is more than what is predicted by the CSM.

Definition 3
We will call "good" the choices for which both the dark matter and the dark energy are explained, at least for some fraction, i.e.
Even better choices are whose which fully explain the dark matter and/or the dark energy, i.e. T DM0 = 0, T 0 = 0.

Linearized Einstein equations
Let us summarize one of the main results of [24].
in the harmonic gauge where the "metric perturbations" A, B, C follow the PDEs The box operator denotes the flat d'Alembertian The source q(x; τ ) has nothing to do with the deceleration parameter q 0 , but it is an expression for the (irrotational) velocity field A general solution of the linearized Einstein Equations has also a wave term g μν − 2Hġ μν + 2(Ḣ − 2H 2 )g μν = 0; but we will not consider it, since we are seeking for selfconsistent choices, so we want that the linearized Einstein Equations have a unique solution. From now on, we will call as "the linearized Einstein Equations" the system (3.3).
Near τ I , the matter inhomogeneities cannot have yet generated the metric perturbations A, B, C. For this reason, as initial conditions for (3.3) we ask that A, B, C are zero at τ I .

Gauge-invariant quantities
In order to get effective parameters F M0 , F 0 with a physical meaning, we must show their gauge independence. In perturbative cosmology, the gauge represents a map which identifies points of the background manifold with those of the physical manifold, so that the metric perturbatioñ g μν := g μν −ḡ μν is well defined (see [26, §4]). Hence, the gauge is a different matter from the coordinates choice in GR; in particular, a gauge change x μ := x μ + ξ μ (x) has ξ μ of a first order in perturbation.
Since we want to use only gauge-invariant fields, we consider the Bardeen's potentials, that for a metric as (3.1) have the form Notice that the second gauge condition (3.2) guarantees that A crucial observation is that the Bardeen's potentials result to be the components of metric perturbation once it is expressed in Newtonian gauge.

Lemma 4.1 Via the transformation
where the gravitational potential coincides to both the Bardeen's potentials.
Hence, from now on we will use this gauge, without writing the primes. The gauge-invariance of its components guarantees the gauge-independence of the following calculations.
Notice that the newtonian gauge has some nice properties, for our purposes: it is already diagonal, while we want to compare its average to a FLRW metric, hence a diagonal one. With another gauge choice, e.g. the harmonic one, the offdiagonal terms would vanish just for the averaging procedure, which would be acceptable but less elegant. The framework of "fake dark matter from GR" finds dark matter effects from the comparison between a "usual" model, i.e. a Newtonian one for local DM and a Friedmannian one for global DM, and a "true" inhomogeneous model; such comparison result to be more natural if the true model has a diagonalized metric, since both the Newtonian approximation of GR and the FLRW metric are diagonal.

Averaging the metric
What we will compare with the CSM is just the average of the metric, since the metric itself is not homogeneous and never allow for an exact equivalence. Such an average depends only on time where we know Now we recall another result from [24]. Let us consider the Green functions for (3.3) where the second-order operators for gravitational waves are defined: Notice that, for such homogeneous operators, the Green functions are completely described by the dependence on x. E.g. G τ (x; τ ) represents the retarded potentials for A(x; τ ) generated by a delta source localized at the point 0 and at the instant τ .

Lemma 4.2 Let us assume the separation of variables for the matter inhomogeneitỹ
Then we can express the average of metric distortions as follows (4.10) Proof L is a linear operator, hence the field generated by a source 4π Gρ can be expressed as the superposition of all Green functions, generated from any point x ∈ R 3 of the space and at any times τ after τ I .
where we exploited the causality of wave equations in the second passage, and the isotropy of G w.r.t. r in the third passage. Now we average A on a region D. As we saw in Sect. 2.3, since it is a first order quantity its average can we written at the first order as where we exchanged the integrals in r and x.
For D → R 3 , the integral inρ 0 becomes ρ 0 and is pulled out, so that the remaining integral coincides with the definition of u A .

Corollary 4.3 The separation of variables does not hold exactly for A, but let we can approximate
Then, in the same way (4.14) Another calculation of C will be shown in Sect. 5.2.
Here we use the u functions to describe the time evolution of the perturbations.
The separation of variables forρ holds when there is a single component dominating. We can express it with the density contrast E.g. when the matter dominates, it is as we know from [26][27][28].
When to dominate is dark energy, the matter structures are ripped apart with the same expansion rate of the universe When to dominate is radiation, the density contrast is well described by as [26][27][28] say again, and a R M is the value of a for which the matter starts to dominate on the radiation; thus the T function is Since we are in the Newtonian gauge, we need also the average of B. We can obtain it averaging the second gauge condition (3.2).
which proves the assertion.

Formulas for the fictitious components
The fictitious components are determined by (2.8). We can rewrite it using the auxiliary variables ract and sum, defined as in [24]. (4.17) For an evaluation of these, we need to know the perturbations of on q 0 and T w0 . The magnitude of the perturbations is determined by the comparison with the CSM metric The conditions at the present time are From the first of these, we obtain the value of Now, we can consider a as the time variable. By now, we denote with a prime the derivatives with respect to a. From so that for any given quantity Q depending on the time, we have Using the relation in (4.19), we can find firstly the perturbations of T w0 . Indeed, from definition (2.5) From the second equation in (4.19), we can compute For any w it's T w0 = H w0 , with the exception of T M0 = H M0 + I M0 . Thus As for the perturbation of q 0 , we must remember that its zeroth order part is not zero, in general, but the background has a deceleration It is distorted by the perturbation, then at first order we expect to have for some coefficient q . We can compute each of these from the third of (4.19), obtaining In particular, this means that Together with (4.20), this gives where the auxiliary quantities are 3) are a formidable mathematical task. We obtained in [24] an exact solution for the case with constant coefficients, but it seems to be impossible an analytical solution when the coefficients depend on τ . However, here we are investigating only the global effects, which depend only on the average of A, B, C, as (B.11) shows.
Let us notice again that L, L C are linear differential operators, depending on time but homogeneous in space. Hence, the average of the solution of (3.3) is equal to the solution of the average of (3.3); in other words, (3.3) are transparent to the spatial average, without backreaction terms, because all the non-linearity appears in the time dependence of the coefficients. Performing a spatial averaging procedure on the PDEs (3.3), these are replaced by simpler ODEs, depending on the time only, for A , B , C . Such ODEs admit analytical solutions.
From Lemmas 4.2 and 4.4 we have where the u functions are studied in Appendix A. From those results, we get Thevs obey equations like (A.6). Writing them for the us we have the assertion.

Single component cases
It is still impossible to solve analytically the evolution (2.6) for a and the ODEs (4.28) for a general choice of components {¯ w0 } w . Moreover, for such a general choice it's quite difficult to determine the form of the sourceρ ∝ T (τ ). However, we are able to solve exactly the equations when a single component¯ w dominates. We can approximate the general evolution as a succession of "epochs"; during each epoch, we consider just the dominant component so that each epoch has a single-component evolution. The full evolution is obtained sticking the partial functions, imposing that a(τ ) ∈ C 0 (τ I ; τ F ), since (2.6) is first order, A , C ∈ C 1 (τ I ; τ F ), since (4.28) are second order; and B ∈ C 0 (τ I ; τ F ), because u B is obtained by an integral in (4.16).

The first selfconsistency condition
Let's start solving (2.6) for a general epoch where c is an integration constant. We get immediately the coefficients of (4.28) Recalling (2.5) and that a(τ ) is increasing (at least) near τ I , we see that the epochs must be in order of decreasing w. In particular, during the first epoch it dominates w M := max{w}. Setting the initial condition By definition it is always α = 0 for definition. From the previous Theorem, we get immediately In particular, a selfconsistent universe must develop the metric perturbations as described by (4.28), with non constant coefficients.

Remark 3
In [24] we studied the constant coefficient case, filling the universe with an exotic component s.t. w = − 1 3 . This breaks the First Selfconsistency Condition, which explains the divergences we found in [24, §4.3]: it is the contribution of r v(r ; τ )| r =R(τ ) ≡ ∞. It is possible to extract finite results even when the I SC is broken, as we did with a renormalization via analytic continuation. A general renormalization method could be to always neglect the term r v(r ; τ )| r =R(τ ) ≡ ∞, using (4.28) for any w M .
As long as the I SC holds, we can fix τ I := 0 without lost of generality.

Decoupling
As we say in Lemma 4.2, for general coefficients of (3.3) we have just an approximated solution of C . This is due to the coupling between C and A. Another advantage of the single component evolution is to allow the decoupling the PDEs of A and C Let α = − 1 2 . 2 Then it is convenient to define the auxiliary field D := A + (2α + 1)C, (5.6) which must satisfy the PDE All the results in Sect. 5 hold true for D, so that From these we get an exact formula for C Remark 4 Notice that in the dark energy epoch α = −1 and the ODE for u D is free of source. However, this doesn't imply that u D is zero, thus in general C = A .

Solving the ODEs
To solve (4.28) for a general w, we need the form of T (τ ). We will assume δ M ∝ a(τ ) n , (5.10) with n(w) a regular function, of which we know n(0) = 1 and n(−1) = 0. This assumption does not certainly hold for the radiation epoch (w = 1 3 ), when Let us start by solving for u A . In general, it has a term u I A generated by the source −a(τ ) 2 T (τ ) = −a(τ ) n−1 , and a term u H A without sources. They result to be The exponent of u H A is 14) It has an imaginary part if and only if Because of the arbitrariness of the integration constants c 1 and c 2 , we can write in general where ξ := 3α 2 + 3α − 1 4 . The solution for D is simpler.
. (5.16) Using this, we get for C The evolution of B is determined by Lemma 4.4.

Particular components
In the following sections, we will need the single-component solutions for some particular components. The dark energy has w = −1 < w − . The perturbations evolve as The matter has w = 0 ∈ (w − ; w + ). The perturbations evolve as For the peculiar evolution during the radiation epoch, we don't use T = a n . The perturbations evolve as (5.24)

Other selfconsistency conditions
Recalling our definition of a "selfconsistent"universe, the First Selfconsistency Condition ensures that there exist finite solutions for A , B , C . We must require also that these solutions are unique and that they describe small enough perturbations. The initial conditions for ( For the I SC one has α(w M ) > 0, so that the first limit implies the second one.

Theorem 5.2 (Second Selfconsistency Condition) A selfconsistent choice of components must be such that n(w M
Remark 5 For a monotonically increasing n(w), the II SC is equivalent to w M > w 0 (5.26) for some limit value w 0 . We can estimate it with a linear interpolation n(w) ∼ = 1 + w, that gives w 0 ∼ = − 1 4 . With more generality, remembering from [24] that 3 n(− 1 3 ) ∈ (φ; 1) and that n(0) = 1, we have The II SC is not necessary if w M = 1 3 , for which always and this is enough to forbid all the values w M / ∈ (w − ; w + ). For the allowed values w M ∈ (w − ; 1], the integration constants for u C are fixed to zero as well, since (5.17) has the same functional form of (5.15). From (5.19) and (5.25) we see that also u H B is fixed to zero, so that the metric perturbations are unique.
The Three Selfconsistency Conditions we proved allow only a "selfconsistency interval" for the component dominating near the Big Bang: (5.30)

Remark 6
Our universe contains certainly radiation and matter as homogeneous components, and probably dark energy. The biggest w is that of radiation, and w M = 1 3 is included in the selfconsistency interval. This is not obvious. Some universes, as the "constant coefficient universe" studied in [24], break the Three Selfconsistency Conditions. The selfconsistency of our universe provides some empirical reinforcement to our model. When I and III SC hold, the requirement of selfconsistency is reduced to asking that the perturbations are small enough to neglect orders higher than the first. This constitutes a last Condition.

Lemma 5.5 (Fourth Selfconsistency Condition) A selfconsistent choice of components must have an inhomogeneous matter such that I M0
T M0 , and such that ∀t ∈ [0; t 0 ] : Proof The first requirement on I M0 is the same we asked in Sect. 2.4. The other requirement is evident from (4.4), where 2 is the perturbation of the metric, and must be smaller than 1.
This is no more a Condition on w, but on I M0 , so that the selfconsistency interval remains the same. Indeed, A , B , C are proportional to I M0 , and so is : the IV SC defines a maximum value M I M0 for the inhomogeneity.

Remark 7
Notice that the IV SC does not imply the II SC, since in the limit case w M = w 0 , doesn't tend to zero, but, nevertheless, it could be small.

The 1-manifold of possible universes
Until now, our computations concerned a general choice of components {¯ w0 } w , for which we found the Selfconsistency Conditions. Now we will apply this general method to our universe.
It contains just three components: radiation R0 , matter Notice that these are not independent, since¯ R0 +¯ M0 + 0 := 1. We have only two independent constraints from 2sum +ract 3 However, we have three unknown parameters: the inhomogeneity I M0 and other two among¯ R0 ,¯ M0 and¯ 0 . This means that the components of our universe are not completely determined by (6.2), but we will find more possible solutions, when a parameter changes. We choose¯ M0 ∈ [0; 1] as parameter, with is determined by the last independent constraint of (6.2). We will have to check which of these values of¯ M0 gives selfconsistent (i.e., if for them hold the IV SC and that I M0 T M0 ), acceptable and eventually good solutions.

Epochs of evolution
Applying (2.5), So we can get the values of a for which the matter starts to be more than the radiation, and the same for other couples The evolution of the universe until now is for 0 ≤ a ≤ a 0 ∼ = 1. During this time, there may have been three or two epochs, depending on the values¯ R0 ,¯ M0 and¯ 0 .  Proof Since the radiation exists, we know¯ R0 > 0, so that a R M > 0 and a R > 0 for any values of¯ M0 ,¯ 0 . Thus, we have always a radiation epoch, which is the first one after the Big Bang. The presence of other epochs depends on our parameter: the quantity of homogeneous matter¯ M0 . We will not consider the case with only the radiation epoch, because it would mean that a R M , a R ≥ a 0 ∼ = 1, which happens for high values of¯ R0 ; but we know from the measures [29] that the radiation is far more less than the matter. Moreover, if the homogeneous matter would be so little, it would mean that I M0 ∼ = T M0 , that is not selfconsistent.
Let us consider the two cases with a matter epoch. From (5.1) we get the background evolution where the continuity determines On the other hand, in the case such that there is no matter epoch, the evolution is where the continuity determines

Searching for good solutions
For any value of the free parameter¯ M0 , we can get a numerical solution of¯ R0 . Following Sect. 6.2 we have, for any chosen value, the evolution of a(τ ), and from the formulas of Appendix B the quantities ract and sum, and thus I M0 , F M0 and F 0 . Imposing (6.2), there could be one or more solutions for R0 , or no one, depending on¯ M0 . For any solution, we have to check if it is acceptable. The selfconsistency checking will require to compute the evolution of , since we have to find that its maximum is less than 1 2 . Getting a set of selfconsistent and acceptable solutions, we will seek if some of them are also good.
Applying this planwork with a numerical algorithm, we find that for a generic¯ M0 there are up to two acceptable values of¯ R0 . E.g. we can find The set of solutions with¯ R0 ∼ 10 −4 have a radiation density quite near to the value of the CSM. We can call them the "principal" solutions, and "secondary" solutions the others. Indeed, following these solutions with continuity, for M0 = 0.315 ∼ = M0 we find trivially The secondary solutions are not selfconsistent, since all of them have I M0 > 99% · T M0 , so that they break the Cosmological Principle. Moreover, the secondary solutions have quite big perturbations 2 max > 0.5: they are smaller than 1 anyway, but not small enough. On the other hand, the principal solutions are selfconsistent. I M0 T M0 becomes greater as¯ M0 runs away from M0 , but it is always less than 45%. It is the same for 2 max , which is always smaller than 0. 28 1. However, most of the principal solutions are not good. We find just a little interval around¯ M0 ∼ = 0.2 for which are explained some fraction of both dark matter and dark energy (Fig. 1). For the values¯ M0 = 0.2,¯ R0 ∼ = 10.03 × 10 −5 ,¯ 0 ∼ = 0.8 it is For¯ M0 > 0.2 we start soon to have F M0 < 0, so that the solutions are no more good. For¯ M0 < 0.2 it starts vice versa to be F 0 < 0, and the solutions are no more good as well.

Searching for solutions without dark energy or dark matter
We can seek if there is a selfconsistent and acceptable solution which fully explains the dark energy as fictitious. From the last paragraph, we know it would require an high¯ M0 , for which F M0 < 0 and the dark matter is more than in the CSM.
The condition of nonexistence of dark energy is¯ 0 := 0, to that it is automatically fixed¯ R0 = 1 −¯ M0 (Fig. 2). The In such a case wefind On the opposite, we can seek if there is a selfconsistent and acceptable solution that fully explains the dark matter as fictitious. From the last paragraph, we know it would require a small¯ M0 , for which F 0 < 0 and the dark energy is more than in the CSM. The condition of nonexistence of dark matter is¯ T M0 := M B0 (Fig. 3). The corresponding value of¯ R0 is fixed by (6.2), which we solve numericallȳ We summarize with the following graphics the numerical results about the principal solutions (Figs. 4,5,6,7,8). In this paper we developed the framework of [24], managing to apply it to a model of our universe, complete with all the components. The large number of variables leaves a free parameter, depending on which we found a one-dimensional set of possible solutions. Within this interval, more dark matter is explained less as less dark energy is, and vice versa. At an end of the range, dark matter is fully explained as a relativistic effect, but the same effect caused an underestimation of dark energy in the Cosmological Standard Model. At the other end the numbers are analogous, with dark matter and dark energy exchanged. For a particular value, both dark matter and dark energy found a partial explanation. An additional measure would determine which is the right solution, but anyway a correction of the parameters of the CSM is required.
Better measures of the density parameters will improve our estimations of F M0 and F 0 , but they can't fix the right parameter¯ M0 . The difference between¯ R0 and the measured R0 , e.g., is not matter of measure precision, but of the factor H 0 H 0ã 0 2 , which concerns the background universe and is not measurable. Rather, a measure of the actual gravitational force ∇ (x; t 0 ) could put the restraint we need. Another possible measure could be the estimation of the mat-ter inhomogeneity at large scale I M0 , i.e. the deviation from the exact Cosmological Principle.
We approximated our calculations in many points. To overcome them would be an improvement of the framework. Solving numerically the evolution a(τ ) it would not be necessary any sticking, which presumably would return more regular graphics than in Sect. 8.3; but recall that this would require the form of T (τ ) for the multi-component case. Even if we found 2 max always far smaller than 1, it could not be considered fully negligible, so that an higher order calculation could provide some relevant corrections. Moreover, we assumed a spatially flat background metric and an irrotational matter, which is not the most general framework.
In the present article we considered the global dark matter effects only. Our cosmological model requires also the calculation of the local effects, to be empirically verified. This needs to overcome the averaging of g μν , and the distribution of fictitious dark matter would depend on the spatial distribution of inhomogeneitiesρ(x; t 0 ). A study of such distribution could start from the fractal properties of the matter structures at large scales [31,32]. The fluctuations of the resultant potential (x; t 0 ) should be compared to the dark matter halo of the galaxies.
The study of local dark matter effects would provide corrections to the standard Newtonian approximations for the dynamic of galaxies and clusters. For such calculations, we cannot assume an irrotational matter as we did here. The rotation of galaxies could provide a rotational term for the non diagonal components of the metric B, which contributes to fictitious gravitational effects [18,19]. Finally, the total amount of the local fictitious effects could be compared to the global fictitious effect F M0 we found here, and the equivalence between them could be the additional restraint we need to fix uniquely the parameter¯ M0 .

Data Availability Statement
The manuscript has associated data in a data repository. [Authors' comment: The manuscript has associated data, see https://arxiv.org/abs/2006.00319.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendices A ODEs for the averaged metric components
Here we study the quantities expressed with the Green functions for non constant parameters We seek for simple ODEs which regulate such u functions.

A.1 Reduction of the dimensions
Since the Green functions are symmetric under spatial rotation, we can reduce the spatial dimensions to one (A.1) and the same for G C τ (r ; τ ). This allow us to express in another way the terms as . Now let us define the auxiliary field and similar for v AC and v C . Then, we can prove Lemma A. 1 The metric perturbations evolve as where the v fields solve the 2D PDEs τ (r ; τ ) satisfies a wave equation, so it holds a causality principle ∀|r | > τ − τ : τ (r ; τ ) ≡ 0.
so that the boundary terms ofv A andv A must vanish. Now, we can check by substitution Notice that this PDE is symmetric under r → −r , so that v A (r ; τ ) = v A (−r ; τ ). This allows us to write The proof is analogous for v AC and v C .
We will not try to write explicitly these τ (r ; τ ) or v(r ; τ ). Instead, we will exploit their properties to get simple characterizations for their averages. For an evaluation ofv A (ω; τ ) when ω goes to infinity, we can manipulate the corresponding ODÊ A solution iŝ v A (ω; τ ) ∼ ω→±∞ 0; which proves the assertion. The proof is analogous for the others vs.
For the other term, we don't need the Fourier transform.
Putting the last lemmas all together, we obtain the us as solutions of the ODEs of Theorem 2.5.5.

B Explicit evolution and fictitious components
From the results of Sect. 6, now we can get the formulas for the evaluation of fictitious dark matter and dark energy. They are different for the three cases described in Lemma 6.1. Let us start with the case where we have all the threeepochs.