On covariant perturbations with scalar field in modified Gauss-Bonnet gravity

We investigate cosmological perturbations of f(G) gravity in the presence of a scalar field. Using the 1 + 3 covariant formalism, we present the energy overdensity perturbation equations responsible for large scale structure formation. After applying harmonic decomposition method together with the redshift transformation technique, we obtain the fully perturbed equations in redshift space. The equations are solved to study the growth of matter overdensities contrast with redshift. For both short- and long-wavelength modes, we obtain numerical results for particular functional form f(G) models and scalar field. We find that, for this choice the energy overdensity perturbations decay with increase in redshift. However, for both short- and long- wavelength modes, the perturbations which include amplitude effects due to the f(G) models with a scalar field do differ remarkably from those in Lambda CDM. The results reduce to GR results in the limit of f(G) tends to G and in the absence of scalar field.


Introduction
Recent observational studies such as high redshift supernovae I a [1,2,3], cosmic microwave background (CMB) [4,5] and baryon acoustic oscillations (BAO) [6] have provided convincing evidence for an accelerated expansion of the Universe.The standard model of cosmology referred to as ΛCDM which includes the cosmological constant introduces inflation addresses several problems in the early universe, including the horizon problem, the flatness problem and the monopole problem, the model is also consistent with various astrophysical and cosmological observations such as the origin of cosmic microwave background radiations (CMB), the formation and distribution of large scale structure, the synthesis of the light elements in the universe and its expansion [7,8,9,10].However this ΛCDM model currently shows limitations in explaining some cosmological issues in both the early and late time Universe.These include the cosmological constant problem, the coincidence problem, the fine tunning problem, and the current problems of Hubble tensions and Sigma 8.The model still cannot explain the nature 1 arXiv:2312.15183v1[gr-qc] 23 Dec 2023 and even the source of both dark matter and dark energy problems.As a result, these shortcomings have encouraged a search for alternative approaches seeking to address some of these problems.These include the use of different scalar fields to describe the cosmic acceleration, they include the consideration of Chaplygin gas models [11,12,13,14,15,16,17,18,19].On the other hand, these latter models present severe constraints on the matter power spectrum in galaxy clusters.Another alternative is to modify gravity on large scale.Several different modified theories of gravity which are extensions of Einstein's gravity in the past have been introduced to address early and late time cosmic problems.They include Scalar-Tensor theories, f (R) models, f (T ) models and f (G) models of gravity, to name but a few, where R, T and G are the Ricci scalar, torsion tensor and Gauss-Bonnet invariant, respectively [20,21,22,23,24,25,26,27,28,29,30,31,32,33].These theories are capable of unifying the early inflationary era and the late-time era in a way similar to the ΛCDM model.It has been suggested that some of these models may have issues to address such as achieving a consistent description of neutron stars, as pointed out in [34,35], dealing with singularity problems, as pointed out in [36,37], controlling matter instability, as pointed out in [38], and satisfying solar system tests, as pointed out in [37,39,40].In some of these models such as 1  R model, it was found [36,41] that the Newtonian limit may not be recovered.For several of these 1  R models, the corrections to Newton's law do not comply with the solar system tests.
Most of these challenges appear to be in absence in f (G) as investigated in [42].As a result, there has been a marked interest in the use of modified f (G) models.For examples, It was demonstrated that scalar field coupled with Gauss-Bonnet invariant and with a potential has no ghosts and is stable [42,43,44].The work conducted by Nojiri et al. [42] investigated the modified Gauss-Bonnet gravity and demonstrated that it may lead to a unified cosmic history.Further, the authors in [45,46] demostrated that some scalar-Gauss-Bonnet gravities may be compatible with the known classical history of the universe expansion (radiation, matter dominance, transition to deceleration and acceleration).In [46], the authors proposed a unified f (R)-scalar-Gauss-Bonnet gravity for dark energy models.A reconstruction program for such models was later developed.It was further shown that Gauss-Bonnet term may play an important role in another class of gravitational models where it couples with scalar field kinetic term.
In [47], the authors focussed mainly on a model of the form R + f (G) in order to avoid the presence of primordial superluminal perturbation modes.That work finds an appropriate model of Gauss-Bonnet gravity which may describe the dark energy era,it also demonstrated that in some Gauss-Bonnet models, it is possible to provide a unified description of inflationary and dark energy era with the same model.
In the work done by [33,48], the authors presented the phenomenology of the late time universe of a given scalar-tensor theory in the context of string inspired gravity.This theory involves a scalar field which is minimally coupled with the Ricci scalar and a function of f (G) [39,41,43,49,50,51,52].Here the coupling between the scalar field and the Gauss-Bonnet invariant was ignored for simplicity reasons.The combination of both a scalar field and f (G) models has the advantage that a unified description between early and late time epochs may be achieved.One possible scenario involves a dominant scalar field during the inflationary era while in the late time era, the curvature corrections in the Gauss-Bonnet invariant can drive the accelerated expansion of the universe.The validity of these theories has been studied and can provide the viable phenomenological models in the primordial late time era [53].There is a need to study whether this combination at perturbative level can produce significant results for large scale structure formation and can be treated as the form of dark energy to source cosmic acceleration.This can best be approached through the study of 1 + 3 covariant perturbations using the combination of both scalar field and Gauss-Bonnet f (G) model.
There are two approaches to study cosmological perturbations namely metric based approach [54,55,56,57,58] and the 1 + 3 covariant approach [59,60,61,62,63].In the latter approach, the perturbations defined describe the true physical degrees of freedom and there are no unphysical modes present.This approach has been used to study the cosmological perturbations in different contexts of general relativity and different alternative theories of gravity [64,65].For instance the work done in [65] considered the covariant form of the field equations of f (R) gravity as a subclass of scalar tensor to study the linear cosmological perturbations.In [66] they analysed the evolution of density perturbations in a particular spatially flat cosmological model with a vacuum decay.In [20], the authors considered the perturbation dynamics in modified Gauss-Bonnet gravity theory at first-order and found that at perturbation level, the small scale dark matter density perturbations grow much more quickly than in the ΛCDM paradigm, which might lead to a strong scale dependent matter power spectrum.
In the previous research works [14,67,68,69,70], we applied 1 + 3 covariant perturbation in the context of f (G) gravity, and we found that the energy overdensity perturbations decay with increase in redshift.In the present work, we apply the 1 + 3 covariant approach to a mixture of the scalar field and Gauss-Bonnet fluids.We consider the fluids to be non-interacting and study the perturbation of gradient variables of scalar field and Gauss-Bonnet fluids in addition to the gradient variables of physical standard matter in both short and long wavelength modes.In this context, the consideration of two different cases is made while analysing the energy density perturbations.In the first case, we solve the whole system of energy density perturbation equations to analyse the large scale structure formation scenario, where the matter perturbation equations couple with the perturbation equations for scalar field and Gauss-Bonnet energy densities.For comparison purposes, we show that the energy density perturbations decay with increase in redshift for both GR and the considered scalar field assisted f (G) gravity model.In the second case, we consider the quasi-static approximation technique to study the implication of perturbation equations on small scale structure.In this approximation, we assume the very slow fluctuations in the scalar field and Gauss-Bonnet energy densities and their momenta, compared to matter energy density.As a result, the matter energy density decouple from the energy densities of the scalar field and Gauss-Bonnet fluids.
The rest of this paper is organised as follows: in Sec.(2), a review of cosmic dynamics in GR is presented; in Sec.(3), we present the 1 + 3 covariant form of the field equations of the scalar field assisted f (G) gravity; in Sec.(4) we define the fluid gradient variables and presents their evolution equations.We analyse the perturbation equations for matter fluctuations for both GR and f (G)-scalar field models in Sec.(5) and Sec.(6), respectively.Sec.(7) gives closing remarks.
Natural units in which c = 8πG = 1 will be used throughout this paper, and Greek indices run from 0 to 3. The symbols ∇, ∇ and the overdot .represent the usual covariant derivative, the spatial covariant derivative, and differentiation with respect to cosmic time, respectively.The (− + ++) spacetime signature will be considered.

Cosmic dynamics-evolution equations in GR
On the large scale structure of the universe, the homogeneous and isotropic assumptions imply that the current universe is close to a flat geometry with radius of curvature R, same at every point in space and that the universe expansion has to be the same at every space with the scale factor a(t).The relationship betwween curvature and matter content of the universe is given by the Einstein's equation represented as where G µν is the Einstein tensor, R µν , R = g µν R µν , T µν and G N are Ricci tensor, Ricci scalar, energy momentum tensor and Newton gravitational constant respectively and g µν is the metric tensor.For a perfect fluid, the energy momentum tensor is given by where, ρ and p are the energy density and isotropic pressure respectively and u µ is the 4-velocity.
The Friedmann-Robertson-Walker (FRW) metric is given by where a(t) is the scale factor governing the expansion of the Universe and it is associated with the dynamics of the universe.For the metric of the form eq. ( 3) and considering a perfect fluid in a flat goemetry, the Friedmann, acceleration and the continuity equations are represented as where H is the Hubble parameter.The eq. ( 4) tells us that the universe containing matter has to be dynamically evolving.Knowing that the critical density of matter is given by ρ c ≡ 3H 2 8πG N and the density parameter Ω i = ρ i ρc , eq. ( 4) can be rewritten as i where Ω i is the density parameter of the matter species present in the universe.For a barotropic perfect fluid with an equation of state parameter given by w = p ρ , and from eq. ( 4) and eq.( 5), one gets [71] whereby for a universe dominated by a dust with an equation of state parameter (w = 0), or radiation (w = 1 3 ), yields a(t) ∝ (t − t 0 ) 2 3 , ρ ∝ a −3 and a(t) ∝ (t − t 0 ) 1 2 , ρ ∝ a −4 respectively, leading to a decelerated expansion of the universe as presented in eq. ( 6).The accelerated expansion (a(t) > 0) occus for w < − 1 3 .
3 The 1 + 3 covariant form of the field equations of the scalar field assisted f (G) gravity The gravitational action involving a scalar field assisted by f (G) gravity is presented as [15,20,33,71,72] where f (G) represents an arbitrary function depending on the Gauss-Bonnet invariant G given by , where R µναβ is the Riemann curvature tensor, g is the determinant g µν and κ is the gravitational constant. 12 ▽ µ ϕ ▽ µ ϕ and V (ϕ) represent the kinetic term of the scalar field and the scalar potential respectively.For the case f (G) = G, d 4 x √ −gG = 0, and there is no scalar field consideration, we recover the gravitational action In this context, eq. ( 12) produces the field equations represented as where T tot µν represents the energy momentum tensor of the total fluids.In this paper, we use the 1 + 3 covariant formalism, where the spacetime is split into temporal and spatial components with respect to the congruence.This formalism treats the slicing of four dimensional space by time and hypersurface.For congruence normal to a spacelike hypersurface, the 1 + 3 covariant formalism reduces to 3 + 1 [73].The basic structure of the 1 + 3 covariant formalism is a congruence of one dimensional curves, mostly timelike curves [67,73,74].In this case, we decompose spacetime cosmological manifold into time and space submanifold separately with a perpendicular 4-velocity field vector u a so that with the metric tensor related to the spatial component as For the metric of the form eq. ( 3) and considering a perfect fluid in a flat goemetry, the Ricci scalar and Gauss-Bonnet parameter are presented as where H = ȧ a is the Hubble parameter.The dot describes differentiation with respect to cosmic time t.The variation principle in the gravitational action with respect to the metric g µν , the scalar field ϕ, and G produces the field equations for f (G) gravity which can be presented as where ρ m and p m represent both relativistic matter (photons, neutrons) and non-relativistic matter (baryons, leptons, Cold Dark Matter ), ρ ϕ and p ϕ represents the matter contribution from the scalar field, ρ G and p G represents the contribution from the Gauss-Bonnet term.The energy density for scalar field contribution and for Gauss-Bonnet contribution and the their pressures are presented as [20,33,71] The matter follows the perfect fluid assumption for any cosmological era of interest.The pressure of the perfect fluid is presented as where m specifies relativistic or non-relativistic matter, whereas the pressure of the total fluids is given by where The continuity equations in the context of scalar field assisted f (G) gravity are given by with the equation of state parameter for matter fluid, Gauss-Bonnet and scalar field fluids are given respectively as Throughout this work, we assume w m to be a constant and w G and w ϕ dynamically changes.The Raychaudhuri equation θ = − 1 3 θ 2 − 1 2 (ρ + 3p) + ▽ a ua that governs the expansion history of the universe in this case is given by with ua as the 4-acceleration in the energy frame of the total fluids given by ua In the next section, we define the gradient variables which help us to obtain perturbation equations governing the evolution of the universe.
4 Fluid description and linear evolution equations Considering the homogeneous and isotropic expanding FRW cosmological background, spatial gradient variables namely matter energy overdensity and volume expansion of the universe can be defined as with m as the running index to specify matter.

Gauss-Bonnet fluids
Analogous to the 1+3 covariant perturbations, the definition of key gradient variables resulting from the spatial gradient variables connected to the Gauss-Bonnet fluids in f (G) gravity can be done as where G and G characterize the perturbations due to Gauss-Bonnet invariant G and its momentum Ġ and describe the inhomogeneity in the Gauss-Bonnet fluids.

Scalar field fluids
The gradient variables responsible to the perturbations due to scalar field ϕ and its momentum φ are presented as These defined gradient variables eq. ( 32) through to eq. ( 37) can be used to develop the system of cosmological perturbations of the scalar field assisted f (G) gravity in the 1 + 3 covariant formalism.Using scalar decomposition method, the scalar gradient variables can be extracted from eq. ( 32) through to eq. ( 37) and presented as The time derivative of each gradient variable helps to know how each evolves and help us to get linear perturbation equations responsible for large scale structure formation in the context of the scalar field assisted f (G) gravity.

Derivation of linear perturbation equations
The time derivative of eq. ( 38) and making use of eq. ( 30), eq. ( 31) and eq.( 79) where the parameters such as b 1 , τ are presented in the appendix.We consider the scalar potential V (ϕ) of the exponential form [71,75] where V 0 is the constant, ϕ is the scalar field given by [71] 1 , and Mp , where M p is the Planck mass, m giving rise to the power-law expansion of the scale factor a(t), so that This potential can represent the late-time acceleration phase of the Universe.For analysis purpose, we apply first the harmonic decomposition method [14, 62, 67, 76, 77, 78] where k is the wave number and Q k (⃗ x) is the eigen functions of the covariant derivative.The wave number k specifies the order of harmonic oscillator related with cosmological scale factor as k = 2πa λ , where λ is the physical wavelength of perturbation.Then we appy the redshift transformation method using [16,33,62,68,69,70,79] so that eq. ( 39)-( 47) with little algebra can be represented in redshift space as and The Eq. ( 57) through to Eq. ( 59) remain key equations for analysing the growth of energy density fluctuations for the purpose of explaining the large scale structures formation.For pedagogical purpose and for the sake of simplicity, we consider three different f (G) models, namely model1, model2 and model32 , respectively considered in [33], [41,42,46] and [37,40], where α 1 , α 2 , β, f 0 , c 1 c 2 , d 1 and d 2 are constants and β > 0, n > 0 and G ̸ = 0, and the scalar field (eq.49) with its derivatives (l 1 , l 2 , l 6 ), the Hubble parameter and the Gauss-Bonnet parameter (eq.17) with its derivatives (l 3 , l 4 , l 5 ) presented in redshift space as with and l 6 = −H to analyse the growth of matter fluctuations in the context of the scalar field assisted f (G) gravity and GR limits.For different realistic f (G) models, the work presented in [36,37,80,81] are recommended.

Matter density fluctuations in GR
For the case f (G) = G, ( when α 1 ∼ 0 and α 2 ∼ G 1 2 , f 0 ∼ 1 and β ∼ 1, c 1 = d 2 = 1, n = 1 and d 1 = c 2 = 0, from eq. ( 60)) and without considering the scalar field contribution in Eq. ( 57) through to Eq. ( 59), we obtain which coincides with GR limits.In this case, the matter (relativistic and non-relativistic) is involved in the growth of the energy density fluctuations.The evolution of Eq. ( 57) through to Eq. ( 59) coincides with the GR limits.We analyse the growth of energy density fluctuations for both dust and radiation fluids in GR limits (eq.64-66).

Dust-dominated epoch
We assume the universe is only dominated by dust fluid, with equation of state parameter w = 0, therefore eq. ( 64)-( 66) lead to where the subscript d stands for dust.Eq. (67) admits the numerical results which are presented in fig.(1), where we depict that the matter energy density fluctuations decay with increase in redshift.

Radiation-dominated epoch
By assuming a universe dominated by radiation fluid, with equation of state parameter w = 1 3 , eq. (64-66) lead to where r stands for radiation.Eq. ( 68) admits the numerical solution presented in fig.
(2).During numerical results computation, the relation 4 is useful, and the normalised energy density fluctuations can be defined as where ∆(z in ) is the initial value of ∆(z) at z in = 1100.Since the variation of CMB temperature detected observationally is in the order of 10 −5 at z in = 1100, we use the initial conditions for ∆(z) = 10 −5 and ∆ ′ (z) = 0. From the plots, one can depict that, the matter energy density fluctuations for both dust-and radiation-dominated epochs are decaying with increase in redshift for different values of k.  67) for dust-dominated universe.
Figure 2: Plot of energy density perturbations versus redshift of eq. ( 68) for radiationdominated universe.
6 Matter density fluctuations of the scalar field assisted f (G) gravity In this section, we analyze the growth of energy density fluctuations of the perturbation equations in both short-and long-wavelength modes by considering the case where the universe is dominated by the mixture of dust-scalar field-Gauss-Bonnet fluids and the case where the mixture of radiation-scalar field-Gauss-Bonnet fluids is dominating the universe.

Short-wavelength mode
The growth of fractional energy density perturbations is analyzed within the Hubble horizon, where k 2 a 2 H 2 ≫ 1 in the whole system to present numerical results of energy density fluctuations.The work done in [82] for GR, [76] for f (R), [62] for f (T ) and in [14,67,68] for f (G) considered this limit.

Quasi-static approximation
In the quasi-static approximation, the time fluctuations in the perturbations of the scalarfield energy density and its momentum, and that of the Gauss-Bonnet energy density and its momentum are assumed to be constant with time.In this case, one sets Ġ = Ġ = Φ = Ψ = 0 under this approximation, the linear perturbations equations eq. ( 57) through to eq. ( 59) reduce to where the parameters such as E 1 , l 1 etc are presented in the appendix.In the following subsections, we consider the cases, where cosmic medium is dominated by a mixture of non interacting fluids, namely dust-scalar field-Gauss-Bonnet fluid and that of radiation-scalar field-Gauss-Bonnet fluid to analyse the energy density perturbations within the Hubble horizon.We set the equation of state parameter for dust to be equal to w = 0 in the eq.( 70), to get The matter energy density perturbations in eq. ( 71) do not couple with the energy density perturbations resulting from the scalar field and the Gauss-Bonnet fluids.Eq. ( 71) admits the numerical solutions presented in figs.

Perturbations in radiation-scalar field-Gauss-Bonnet fluids dominated Universe
For the case of radiation-dominated universe, we set the equation of state parameter to w = 1 3 in eq. ( 70) to get The matter energy density perturbations in eq. ( 72) do not couple with the energy density perturbations resulting from the scalar field and the Gauss-Bonnet fluids.Eq. ( 72) admits the numerical solutions presented in figs.( 6), ( 7), (8) for f (G) model 1, model 2, model 3, respectively.The energy density perturbations decay with increase in redshift.The implications of different models on the energy overdensity contrast are presented in table (2) for radiationscalar field-Gauss-Bonnet fluid mixture in the short wavelength mode.
Table 2: Illustration of the features of energy overdensity contrast for radiation-scalar field-Gauss-Bonnet fluids system in short wavelength mode.Models Ranges of parameters for f (G) models behavior of δ(z) decay with increase in z observationally supported model3 1 < n < 2 decay with increase in z observationally supported

Long-wavelength modes
In the long wavelength modes, we assume k 2 a 2 H 2 ≪ 1.All cosmological fluctuations begin and remain inside the Hubble horizon.In this case, we solve the whole system of perturbation equations without considering the quasi-static approximation and present the numerical results in both dust-scalar field-Gauss-Bonnet and radiation-scalar field-Gauss-Bonnet fluid mixtures.72) using model 1 for different values of m of the radiation-scalar field-Gauss-Bonnet fluids system in a short wavelength mode.

Perturbations in dust-scalar field-Gauss-Bonnet fluids dominated Universe
Considering the dust-dominated universe, where w = 0, eq. ( 57) through to eq. ( 59) reduce to The eqs. ( 73)-( 75) are the ones responsible for large scale structure formation for the universe made of dust-scalar field-Gauss-Bonnet energy density perturbations in long wavelength modes.From eqs. ( 73)-( 75), the matter energy density perturbations couple with the energy density perturbations for the scalar field and Gauss-Bonnet fluids.The numerical solutions of eqs.( 73)-(75) were computed using eqs.(60-63) together with l 1 − l 6 .We chose the initial conditions ∆(z) = 10 −5 and ∆ ′ (z) = 0, together with G(z) = 10 −5 , G ′ (z) = 0, Φ(z) = 10 −5 and Φ ′ (z) = 0, and the arbitrary constants as present in different f (G) models.The numerical solutions are presented in figs.( 9), ( 10), (11) for three different models, respectively.From the plot, the energy density perturbations decay with redshift.The implications of different models on the energy overdensity contrast are presented in table (3) for dust-scalar field-Gauss-Bonnet fluid mixture in the long wavelength mode.Table 3: Illustration of the features of energy overdensity contrast for dust-scalar field-Gauss-Bonnet fluids system in long wavelength mode.Models Ranges of parameters behavior of δ(z) Observation model1 α 1 = 100, α 2 = 0.001 decay with increase in z observational supported model2 0 < β < 1 2 decay with increase in z observational supported model3 1 < n < 2 grow with increase in z differ from the ΛCDM predictions

Perturbations in radiation-scalar field-Gauss-Bonnet fluids dominated Universe
In the radiation-dominated universe, the equation of state parameter w = 1 3 is used so that eq. ( 57) through to eq. ( 59) reduce to ( The eqs. ( 76)-( 78) are the energy density perturbations for a radiation-scalar field-Gauss-Bonnet fluid mixture.From eqs. ( 76)-( 78), the matter energy density perturbations couple with the energy density perturbations for the scalar field and Gauss-Bonnet fluids.The numerical results of eqs.( 76)-( 78) were computed using different initial conditions and constants and are presented in figs.(12), ( 13), ( 14) for three different f (G) models respectively.From the plots, the energy density perturbations decay with redshift and show some oscillation features.The implications of different models on the energy overdensity contrast are presented in table (4) for radiation-scalar field-Gauss-Bonnet fluid mixture in the long wavelength mode.

Discussion and Conclusion
In the current paper, we present a detailed analysis of cosmological perturbations for the scalar field assisted f (G) gravity.Using the 1 + 3 covariant gauge-invariant formalism, we derived the linear covariant evolution equations for a flat FRW spacetime background.In doing so, we have defined gradient variables for matter fluids, scalar field and Gauss-Bonnet fluids and derived their vector evolution equations.Since our main interest extends to large scale structure formation, we extracted only scalar parts of the linear evolution equations using scalar decomposition method, understood to play a key role in structure formation.
After getting the scalar perturbation equations, we applied harmonic decomposition method together with the redshift transformation technique to get a set of simplified equations in redshift space responsible for large scale structure formation.For further analysis, we consider two different cases basing on the wavenumber k.For short wavelength modes k 2 a 2 H 2 ≫ 1, we apply the quasi-static approximation to analyze the small scale fluctuations for a system made of both the mixture of dust-scalar field-Gauss-Bonnet fluids and the mixture of radiation-scalar field-Gauss-Bonnet fluids.For long wavelength modes k 2 a 2 H 2 ≪ 1, we solve the whole system of energy density perturbation equations without considering the quasi-static approximation.In both cases, for pedagogical purpose, we consider a logarthmic scalar field ϕ = ln(t), an exponential potential V (ϕ) = V 0 e −αϕ and a polynomial models and by using the initial conditions to solve the perturbation equations.For both cases, we get the numerical results which show that the energy density perturbations decay with redshift for different values of m.Some of the specific highlights of the present paper include: This work presents the energy overdensity perturbation equations Eq. ( 57)-Eq.( 59) in the context of the scalar field assisted f (G) gravity responsible for large scale structure formation.The obtained matter energy density perturbation equation couples with the energy density fluctuations resulting from the mixture of scalar field and Gauss-Bonnet fluids in the long wavelength modes while it decouples in the short wavelength modes and GR limits.The perturbation equations were solved numerically and the results are presented in figs.to study cosmological perturbations.For instance, one can take into consideration different scalar fields in the matter lagrangian in f (R), f (T ) actions.Another application is that, one can use other types of realistic modified gravity models such as one presented in [42] to investigate their impact on the large structure formation.The corresponding analysis may lead to some significant qualitative outcomes in comparison with the ΛCDM.It will be investigated elsewhere.

A Useful Linearised Differential Identities
For all scalars f , vectors V a and tensors that vanish in the background, S ab = S ⟨ab⟩ , the following linearised identities hold: ∇a

B Useful relations
We define whereas for the quasi-static approximation case, we used The relations defined below were used in different parts of the work

4. 1
Definition of vector gradient variables for the evolution of the universe 4.1.1Matter fluids

Figure 3 :
Figure 3: Plot of energy density perturbations versus redshift of Eq. (71) using model 1 for different values of m of the dust-scalar field-Gauss-Bonnet fluids system in a short wavelength mode..

Figure 4 :
Figure 4: Plot of energy density perturbations versus redshift of Eq. (71) using model 2 for different values of β of the dust-scalar field-Gauss-Bonnet fluids system in a short wavelength mode..

Figure 5 :
Figure 5: Plot of energy density perturbations versus redshift of Eq. (71) using model 3 for different values of n of the dust-scalar field-Gauss-Bonnet fluids system in a short wavelength mode.

Figure 6 :
Figure 6: Plot of energy density perturbations versus redshift of Eq. (72) using model 1 for different values of m of the radiation-scalar field-Gauss-Bonnet fluids system in a short wavelength mode.

Figure 7 :
Figure 7: Plot of energy density perturbations versus redshift of Eq. (72) using model 2 for different values of β of the radiation-scalar field-Gauss-Bonnet fluids system in a short wavelength mode.

Figure 8 :
Figure 8: Plot of energy density perturbations versus redshift of Eq. (72) using model 3 for different values of n of the radiation-scalar field-Gauss-Bonnet fluids system in a short wavelength mode.

Figure 9 :
Figure 9: Plot of energy density perturbations versus redshift of Eq. (73-75) using model 1 for different values of m of the dust-scalar field-Gauss-Bonnet fluids system in a long wavelength mode.

Figure 10 :
Figure 10: Plot of energy density perturbations versus redshift of Eq. (73-75) using model 2 for different values of β of the dust-scalar field-Gauss-Bonnet fluids system in a long wavelength mode.

Figure 11 :
Figure 11: Plot of energy density perturbations versus redshift of Eq. (73-75) using model 3 for different values of n of the dust-scalar field-Gauss-Bonnet fluids system in a long wavelength mode.

Figure 12 :
Figure 12: Plot of energy density perturbations versus redshift of Eq. (76-78) using model 1 for different values of m of the radiation-scalar field-Gauss-Bonnet fluids system in a long wavelength mode.

Figure 13 :
Figure 13: Plot of energy density perturbations versus redshift of Eq. (76-78) using model 2 for different values of β of the radiation-scalar field-Gauss-Bonnet fluids system in a long wavelength mode.

Figure 14 :
Figure 14: Plot of energy density perturbations versus redshift of Eq. (76-78) using model 3 for different values of n of the radiation-scalar field-Gauss-Bonnet fluids system in a long wavelength mode.

Table 4 :
Illustration of the features of energy overdensity contrast for radiation-scalar field-Gauss-Bonnet fluids system in long wavelength mode.