Next Article in Journal
Directed Flow in Microscopic Models in Relativistic A+A Collisions
Next Article in Special Issue
Cosmic Inflation, Quantum Information and the Pioneering Role of John S Bell in Cosmology
Previous Article in Journal
Vacuum Effective Actions and Mass-Dependent Renormalization in Curved Space
Previous Article in Special Issue
The Quantum Nature of Lorentz Invariance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Matter Growth in Imperfect Fluid Cosmology

by
Winfried Zimdahl
1,*,†,
Hermano E.S. Velten
2,† and
William C. Algoner
3,†
1
PPGCosmo & Núcleo Cosmo-UFES, Universidade Federal do Espírito Santo, Av. Fernando Ferrari, 514, Campus de Goiabeiras, CEP 29075-910 Vitória, Espírito Santo, Brazil
2
Departamento de Física, Universidade Federal de Ouro Preto (UFOP), CEP 35400-000 Ouro Preto, Minas Gerais, Brazil
3
Departamento de Física, Universidade Federal do Paraná, 81531-980 Curitiba, Paraná, Brazil
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2019, 5(3), 68; https://doi.org/10.3390/universe5030068
Submission received: 2 November 2018 / Revised: 24 February 2019 / Accepted: 25 February 2019 / Published: 4 March 2019

Abstract

:
Extensions of Einstein’s General Relativity (GR) can formally be given a GR structure in which additional geometric degrees of freedom are mapped on an effective energy-momentum tensor. The corresponding effective cosmic medium can then be modeled as an imperfect fluid within GR. The imperfect fluid structure allows us to include, on a phenomenological basis, anisotropic stresses and energy fluxes which are considered as potential signatures for deviations from the cosmological standard Λ -cold-dark-matter ( Λ CDM) model. As an example, we consider the dynamics of a scalar-tensor extension of the standard model, the e Φ Λ CDM model. We constrain the magnitudes of anisotropic pressure and energy flux with the help of redshift-space distortion (RSD) data for the matter growth function f σ 8 .

1. Introduction

The last two decades have seen a tremendous activity within the community of cosmologists and astrophysicists to find a satisfactory explanation for the results of the observations of supernovae of type Ia (SNIa), reported in [1,2,3], which suggested an accelerated expansion of the scale factor of the Robertson–Walker (RW) metric. The apparently simplest way to account for this behavior is to assume the existence of a cosmological constant Λ with a suitable value (see, e.g., [4,5]). From a purely general relativistic (GR) point of view, Λ might be seen as another gravitational constant along with Newton’s gravitational constant. On the other hand, an effective cosmological constant has been associated with the quantum vacuum (see, e.g., [6]). Taking this into account, the actually measured gravitational constant might be a combination of a purely geometric and a quantum contribution. This context gave rise to a number of discussions concerning the cosmic coincidence problem [7,8,9,10,11].
Since accelerated expansion in inflationary models of the early universe is conveniently described in terms of scalar fields, a corresponding description has also been applied to the current phase of the cosmic evolution. This is equivalent to making the cosmological “constant" a dynamic quantity. Along this line, a number of different approaches have been developed which either assume the existence of some form of exotic matter within GR or modify Einstein’s theory of gravity. In the meantime, a huge amount of data of various types have been accumulated [12,13,14,15].
Facing the multitude of models at hand, one might wish to keep an eye on unifying phenomenological aspects of different types of description. It is this intention which motivates the present paper. Our starting point is the observation that any extension of Einstein’s GR can formally be recast into an Einsteinian structure with a suitably defined effective energy-momentum tensor [16,17]. In the presence of an adequate timelike vector field, this energy-momentum tensor is then characterized by an energy density, a scalar isotropic pressure, an anisotropic pressure and an energy flux. It has the structure of the energy-momentum tensor of an imperfect fluid [18,19,20,21]. Imperfect-fluid dynamics can provide an effective description for cosmological models beyond the standard model. Anisotropic stresses and/or energy fluxes are typical ingredients of a non-standard dynamics of, e.g., scalar-tensor theories. A robust phenomenological theory may represent a unifying framework for approaches that differ in their detailed underlying microscopic dynamics. Since the anisotropic pressure determines the gravitational slip, a suitable parametrization may be used to quantify potential deviations from the standard model and, combined with observational data, to put limits on such deviations. A similar comment holds for heat–flux caused deviations.
The resulting gravitational dynamics depends on the (effective) fluid quantities directly, but it does not explicitly depend on the detailed underlying microscopic dynamics which is either unknown or beyond a straightforward and transparent analytical treatment [22]. For a scalar field, e.g., there is a dependence on this field only through the mentioned fluid type quantities, in the simplest case energy density and isotropic pressure, but there is no further direct dependence on the scalar-field details themselves. On this basis, we aim to give a phenomenological fluid-type description of the cosmological dynamics in which different models are determined by a set of phenomenological parameters such as an equation-of-state parameter and an adiabatic sound speed parameter. Additionally, and this is what we consider to be the new aspect of our study, it is necessary to introduce analogous parameters for the anisotropic pressure and the energy flux. We split the total effective energy-momentum tensor into two components, one of them being a pressureless perfect fluid, the second one an imperfect fluid. The pressureless fluid is supposed to account for some form of (dark) matter, the imperfect component is intended to provide an effective description of dark energy. While rather general, such a purely phenomenological description necessarily leaves open the microscopic origin of these quantities.
One of the tools to discriminate between different models of the cosmological dark sector is the growth rate of matter perturbations [23,24,25,26,27]. Competing models with similar behavior in the homogeneous and isotropic background will generally have different predictions for the dynamics of matter inhomogeneities. As an example, we consider the simplest possible scalar-tensor extension of the Λ CDM model. In this minimalist approach, we remain in the vicinity of the standard model at the present epoch [28]. The present paper advances a preliminary study in [29] where the matter growth rate was obtained in a simplified manner on the basis of a rough approximation which both avoided a solution of the full dynamics and neglected the heat flux. Here, we consider the full perturbation dynamics with the heat flux included.
In Section 2, we introduce the energy-momentum tensors of the cosmic medium as a whole and those of the individual components. The general conservation laws are presented in Section 3. Section 4 is devoted to the perturbation dynamics. The relevant metric and fluid quantities are introduced and a modified Poisson equation is obtained. A coupled system of two second-order equations for the entire cosmological dynamics is established from which the dynamics of matter perturbations is derived. All this is applied in Section 5 to the special case of the e Φ Λ CDM model. In Section 6, we constrain the parameters that describe deviations from the standard model by contrasting the predictions of this model with recent redshift-space distortion (RSD) data. Section 7 provides our conclusions concerning the status of the present approach.

2. Energy-Momentum Tensors

We assume the cosmic medium to consist of pressureless perfect-fluid type matter (subindex (m)) and a second component, modeled as a general imperfect fluid (subindex (x)). The cosmic substratum as a whole is then an imperfect fluid as well. It is described by an energy-momentum tensor
T μ ν = T ( m ) μ ν + T ( x ) μ ν .
Greek tensorial indices run over 0 , 1 , 2 , 3 . The first part, T ( m ) μ ν , describes pressureless matter through
T ( m ) μ ν = ρ m u ( m ) μ u ( m ) ν ,
where u ( m ) μ with u ( m ) μ u ( m ) μ = 1 is the matter four-velocity which may be associated with an observer. The matter energy density is
ρ m = T ( m ) μ ν u ( m ) μ u ( m ) ν .
(For scalar quantities, we use the subindex m instead of (m)). The second part, T ( x ) μ ν , is given by the general split with respect to a timelike unit vector u ( x ) μ ,
T ( x ) μ ν = ρ x u ( x ) μ u ( x ) ν + p x h ( x ) μ ν + Π ( x ) μ ν + q ( x ) μ u ( x ) ν + q ( x ) ν u ( x ) μ ,
where u ( x ) μ u ( x ) μ = 1 and
ρ x = T ( x ) μ ν u ( x ) μ u ( x ) ν , p x = 1 3 h ( x ) μ ν T ( x ) μ ν , q ( x ) μ = h ( x ) μ ν T ( x ) ν σ u ( x ) σ ,
Π ( x ) μ ν = h ( x ) μ σ h ( x ) ν τ T ( x ) σ τ 1 2 h ( x ) μ σ h ( x ) ν τ T ( x ) σ τ + h ( x ) ν σ h ( x ) μ τ T ( x ) σ τ 1 3 h ( x ) μ ν h ( x ) σ τ T ( x ) σ τ
with
h ( x ) μ ν = g μ ν + u ( x ) μ u ( x ) ν , h ( x ) μ ν u ( x ) ν = q ( x ) μ u ( x ) μ = Π ( x ) μ ν u ( x ) μ = Π ( x ) μ μ = 0 .
(For scalar quantities, we use the subindex x instead of (x)). Here, ρ x is the energy density of component x and p x is its isotropic pressure. The heat–flux vector is denoted by q ( x ) μ and the anisotropic pressure by Π ( x ) μ ν . In general, the timelike unit vector u ( x ) μ does not coincide with u ( m ) μ . In a scalar-field representation, e.g., it may be related to the timelike gradient of a scalar field. We shall choose both four-velocities to coincide in a homogeneous and isotropic background, but they will differ, in general, on the perturbative level.
For the total energy-momentum tensor, we assume a decomposition with respect to still another timelike unit vector u μ ,
T μ ν = ρ u μ u ν + p h μ ν + Π μ ν + q μ u ν + q ν u μ ,
where
ρ = T μ ν u μ u ν , p = 1 3 h μ ν T μ ν , q μ = h μ ν T ν σ u σ ,
Π μ ν = h μ σ h ν τ T σ τ 1 2 h μ σ h ν τ T σ τ + h ν σ h μ τ T σ τ 1 3 h μ ν h σ τ T σ τ ,
with
u μ u μ = 1 , h μ ν = g μ ν + u μ u ν , h μ ν u ν = q μ u μ = Π μ ν u μ = Π μ μ = 0 .
In this decomposition, ρ is the total energy density, p is the total isotropic pressure, q μ is the total heat flux and Π μ ν is the total anisotropic pressure. The relation between u μ , u ( m ) μ and u ( x ) μ will be clarified below.

3. General Conservation Equations

We shall assume separate energy-momentum conservation of both components. The separate matter conservation equations are
u ( m ) ν T ( m ) ; κ ν κ = ρ m , α u ( m ) α + Θ m ρ m = 0
and
h ( m ) ν α T ( m ) ; κ ν κ = ρ m u ˙ ( m ) α = 0 .
Here, Θ m u ( m ) ; α α and u ˙ ( m ) α u ( m ) ; β α u ( m ) β is the matter four-acceleration. In the homogeneous and isotropic background, Equation (12) reduces to
d ρ m d t + 3 H ρ m = 0 ,
where H is the Hubble rate H a ˙ a and a is the scale factor of the RW metric. In this background, Equation (13) is identically satisfied.
The conservation equations for the x component are
u ( x ) ν T ( x ) ; κ ν κ = ρ x , α u ( x ) α + Θ x ρ x + p x + μ q ( x ) μ = 0 ,
where μ q ( x ) μ h μ α q ( x ) ; α μ and
h ( x ) ν α T ( x ) ; κ ν κ = ρ x + p x u ˙ ( x ) a + p x , ν h ( x ) α ν + κ Π ( x ) μ κ + h μ ( x ) ν q ˙ ( x ) ν + 4 3 Θ h μ κ ( x ) q ( x ) κ = 0
with Θ x u ( x ) ; α α and a four-acceleration u ˙ ( x ) α u ( x ) ; β α u ( x ) β . We have neglected here terms that will not contribute in a linear perturbation theory about a homogeneous and isotropic background with vanishing Π ( x ) μ κ and q ( x ) μ . The pressurefree matter motion is geodesic ( u ˙ ( m ) α = 0 ), while the motion of the x component is generally not. For the background dynamics, one has
d ρ x d t + 3 H 1 + w x ρ x = 0
with an equation-of-state (EoS) parameter w x that has to be specified for a given model.
The general total energy-conservation equation becomes (neglecting contributions of u ˙ μ q μ with u ˙ α u α ; β u β and σ μ ν Π μ ν with σ μ ν 1 2 u μ ; ν + u ν ; μ 2 3 h μ ν u α α which will be of higher than first order in a linear perturbation theory)
ρ ˙ + Θ ρ + p + μ q μ = 0 ,
where ρ ˙ ρ , μ u μ and Θ u ; α α . For the momentum conservation, we have (again taking into account only those terms that contribute in a linear perturbation theory)
ρ + p u ˙ μ + μ p + κ Π μ κ + h μ ν q ˙ ν + 4 3 Θ h μ κ q κ = 0 .
In the homogeneous and isotropic background, q μ = u ˙ μ = Π μ κ = σ μ κ = ω μ κ 1 2 u μ ; κ u κ ; μ = 0 . In our linear perturbation theory, we will neglect products of these first-order quantities.
In the homogeneous and isotropic background, we identify all four-velocities, i.e., u ( m ) a = u ( x ) a = u a . Then,
d ρ d t + 3 H 1 + w ρ = 0
with
w p ρ = p x ρ = w x Ω x , Ω x = ρ x ρ = 1 Ω m , Ω m = ρ m ρ .

4. Perturbation Dynamics

4.1. Metric and Fluid Quantities

We restrict ourselves to scalar perturbations, described by the line element (in the longitudinal gauge)
d s 2 = 1 + 2 ϕ d t 2 + a 2 1 2 ψ δ a b d x a d x b .
The fluid-dynamical system of perturbation equations will be established by adequately generalizing the corresponding steps in [30,31,32,33]. First-order fluid quantities will be denoted by a hat symbol. The perturbed time components (tensorial index 0) of the four-velocities are
u ^ 0 = u ^ 0 = u ^ ( m ) 0 = u ^ ( x ) 0 = 1 2 g ^ 00 = ϕ .
We define the (three-) scalar quantities v, v m and v x by (latin indices denote spatial components)
a 2 u ^ m = u ^ m v , m ,
as well as by
a 2 u ^ ( m ) a = u ^ ( m ) a v m , a and a 2 u ^ ( x ) a = u ^ ( x ) a v x , a ,
respectively. The perturbed expansion scalar becomes
Θ = u ; μ μ Θ ^ = 1 a 2 Δ v 3 ψ ˙ 3 H ϕ ,
where Δ denotes the three-dimensional Laplacian. The corresponding first-order expressions for Θ m = u ( m ) ; μ μ and Θ x = u ( x ) ; μ μ are
Θ ^ m = 1 a 2 Δ v m 3 ψ ˙ 3 H ϕ and Θ ^ x = 1 a 2 Δ v x 3 ψ ˙ 3 H ϕ ,
respectively.
With T ^ 0 0 = ρ ^ as well as T ^ ( x ) 0 0 = ρ ^ x and T ^ ( m ) 0 0 = ρ ^ m , it follows that, at first order,
ρ ^ = ρ ^ x + ρ ^ m .
With
T ^ a 0 = ρ + p u ^ a + q a , T ^ ( m ) a 0 = ρ m u ^ ( m ) a , T ^ ( x ) a 0 = ρ x + p x u ^ ( x ) a + q ( x ) a ,
the first-order relation
ρ + p u ^ a + q a = ρ m u ^ ( m ) a + ρ x + p x u ^ ( x ) a + q ( x ) a
is valid. According to our restriction to scalar perturbations, the vectors q ( x ) a and q a are represented by gradients of scalars, i.e.,
q ( x ) a = q x , a and q a = q , a ,
respectively. This leads to
v = ρ m ρ + p v m + ρ x + p x ρ + p v x + q x q ρ + p ,
i.e., in the presence of heat fluxes the relation between the different velocities is modified compared with the perfect-fluid case.
From the perturbed spatial components
T ^ a b = p ^ δ a b + Π a b , T ^ ( m ) a b = 0 , T ^ ( x ) a b = p ^ x δ a b + Π ( x ) a b ,
we have
p ^ δ a b + Π a b = p ^ x δ a b + Π ( x ) a b .
Since Π a a = Π ( x ) a a = 0 , this reduces to
p ^ = p ^ x , Π a b = Π ( x ) a b .

4.2. Modified Poisson Equation

Taking into account heat–flux effects in the perturbation dynamics modifies the usual Poisson-type equation that relates the comoving density perturbations and the gravitational potential ψ . Coupling the 0-0 and the 0-a field equations leads to
Δ ψ = 4 π G a 2 ρ ^ c 3 H q = 3 2 H 2 a 2 ρ ^ c ρ 3 H q ρ ,
where we introduced the comoving density perturbation
ρ ^ c ρ ^ 3 H ρ + p v .
Then, the (generalized) Poisson equation becomes (in the k space)
2 3 k 2 H 2 a 2 ψ = ε ,
where
ε δ 3 H 1 + p ρ v + q ρ .
The heat–flux contribution has to be included in the generalized definition of comoving energy–density perturbations. The quantity ε represents a generalized comoving energy–density perturbation which takes into account also the existence of a heat flux contribution q in addition to the velocity potential v. Equation (38) with (39) generalizes the well-known corresponding relation for perfect fluids. Either the potential ψ or the (generalized) comoving energy–density perturbation ε may be used as independent variable of the perturbation theory.

4.3. Combination of Conservation and Raychaudhuri Equations

The basic set of first-order perturbation equations can be obtained by adequately generalizing previous fluid cosmological calculations (cf. [30,31,32,33]). The essential ingredient is a combination of the total first-order conservation equations with the first-order Raychaudhuri equation. The result is
ε + 3 2 15 2 p ρ + 3 p ρ ε a 3 2 + 12 p ρ 9 2 p 2 ρ 2 9 p ρ ε a 2 1 a 2 H 2 Δ a 2 p ^ c ρ p ρ l 3 2 1 + p ρ l a + 3 2 a 2 1 p ρ l 2 a Δ Π a 2 ρ 3 1 p ρ + 2 p ρ Δ a 2 H 2 + 2 3 Δ 2 a 4 H 4 H 2 a 2 Π ρ = 0 ,
where l Θ q ρ encodes the heat–flux contribution and the prime denotes a derivative with respect to the scale factor a. So far, neither the pressure perturbations p ^ c nor the quantities l and Π are specified. In a strict sense, they are to be determined from an underlying microscopic theory along with the EoS. If such theory is not available, progress can be made through a phenomenological approach.

4.4. Relative Perturbations

With the help of the quantities
D c ρ ^ c ρ + p = δ c 1 + p / ρ , ρ ^ m c ρ ^ m + ρ ˙ m v , δ m c ρ ^ m c ρ m ,
we define the relative density perturbations (“entropy perturbations”)
S m D c δ m c ,
as the difference between total and pure matter perturbations. The perturbations of the pressure are the sum of an adiabatic part p ˙ ρ ˙ ρ ^ c and a nonadiabatic contribution p ^ n a d ,
p ^ c = p ˙ ρ ˙ ρ ^ c + p ^ n a d .
For our configuration, the nonadiabatic part explicitly becomes
p ^ n a d = p ^ x c p ˙ x ρ ˙ x ρ ^ x c + ρ m ρ x + p x ρ + p p ˙ x ρ ˙ x ρ ^ x c ρ x + p x ρ ^ m c ρ m ,
where
ρ ^ x c ρ ^ x + ρ ˙ x v and p ^ x c p ^ x + p ˙ x v .
The combination p ^ x c p ˙ x ρ ˙ x ρ ^ x c in Equation (44) accounts for the intrinsic nonadiabatic perturbations of the x component while the last term appears due to the two-component nature of the cosmic medium. As a consequence, the fluid as a whole is nonadiabatic even if each of its components is adiabatic on its own. With ρ ^ x c = ρ ^ c ρ ^ m c in the last term on the right-hand side of Equation (44), we may write
ρ ^ x c ρ x + p x ρ ^ m c ρ m = ρ + p ρ x + p x D c δ m c = ρ + p ρ x + p x S m .
Then, the nonadiabatic pressure perturbations are
p ^ n a d = p ^ x c p ˙ x ρ ˙ x ρ ^ x c + ρ m p ˙ x ρ ˙ x S m .
This implies that, through the pressure perturbations, the dynamics of the total energy–density perturbation ε (or δ c ) is coupled to the dynamics of S m (cf. Equation (40)). To obtain the dynamics of S m , we combine the conservation equations for the total medium with those for the matter component. On this basis, the following second-order equation for S m can be derived (again generalizing steps described in [30,31,32,33]),
S m + 3 2 1 p ρ S m a Δ a 2 H 2 p ^ c a 2 ρ + p + 3 a p ^ n a d ρ + p + 9 7 6 + p ρ 1 2 p ρ p ^ n a d a 2 ρ + p 2 3 1 a 2 H 2 Δ 2 Π a 4 ρ + p + 3 a 2 H 2 p ρ H Δ q a 2 ρ + p = 0 .
As in Equation (40), neither the pressure perturbations p ^ c nor the quantities l (or q) and Π are specified. The idea now is to determine the pressure perturbations as well as l and Π in terms of our basic variables ε and S m in order to establish a closed system for ε and S m from which subsequently the matter perturbations can be obtained.

4.5. Perturbations of (an-)Isotropic Pressures and Energy Flux

The pressure perturbations should be considered in the rest frame of the x component. Since the combination p ^ x p ˙ x ρ ˙ x ρ ^ x is gauge invariant, the first part on the right-hand side of Equation (47) may be written
p ^ x c p ˙ x ρ ˙ x ρ ^ x c = p ^ x p ˙ x ρ ˙ x ρ ^ x = p ^ x c x p ˙ x ρ ˙ x ρ ^ x c x ,
with
p ^ x c x p ^ x + p ˙ x v x , ρ ^ x c x ρ ^ x + ρ ˙ x v x ,
where v x is the velocity potential of component x, introduced in Formula (25).
The generalized gauge-invariant comoving energy–density perturbation of the x component is
ρ ^ x ε ρ ^ x 3 H ρ x + p x v x + q x .
This definition for ρ ^ x ε parallels the definition of the comoving energy–density perturbation ε for the medium as a whole (cf. Equation (39)).
The speed of sound is defined as the coefficient that relates pressure perturbations and energy–density perturbations within the rest frame (which in the perfect-fluid case is v x = 0 ). In the presence of a heat flux with scalar potential q x , a more appropriate definition is
p ^ x ε = c ε 2 ρ ^ x ε
with the accordingly defined pressure perturbation
p ^ x ε p ^ x 3 H p ˙ x ρ ˙ x 1 + w x ρ x v x + q x .
Then, the total (isotropic) pressure perturbation on subhorizon scales k 2 a 2 H 2 are given by
p ^ ε ρ p ^ c ρ p ˙ ρ ˙ l = c ε 2 ρ x + p x ρ + p ε + c ε 2 Ω m S m a 2 H 2 k 2 1 .
Obviously, the pressure perturbations couple the dynamics of ε to that of the relative perturbations S m . The entire dynamics is given by the coupled system of equations for ε and S m in which there appear the so far undetermined quantities l and Π . The simplest and most direct way to obtain a closed system consists of requiring l and Π to be linear combinations of ε and S m . This is achieved by a general ansatz for l,
l = 3 H q ρ = α ε + β S m ,
where α and β are phenomenological coefficients, as well as by a corresponding ansatz for Π ,
H 2 Π ρ = μ ε + ν S m ,
with phenomenological coefficients μ and ν . Notice that relations (55) and (56) are constructed to parallel the expression (54) for the isotropic pressure. The coefficients α and β quantify the rôle of the heat flux on the perturbation dynamics, the coefficients μ and ν determine the influence of the anisotropic pressure. As the sound speed parameter c x , the coefficients α , β , μ and ν should be calculable from an underlying microscopic theory. Here, these coefficients are entirely phenomenological. The perfect-fluid case is recovered for α = β = μ = ν = 0 . Phenomenological relations between anisotropic stresses and energy–density perturbations in different contexts can be found, e.g., in [34,35,36]. A different parametrization of perturbations via equations of state was put forward in [18,19,20] on the basis of a general scalar-field Lagrangian. Eliminating internal degrees of freedom, these authors introduced equations of state for the entropy perturbation and the anisotropic stress in terms of perturbations of the density, the velocity and the metric perturbations to obtain closed perturbation equations. Their basic variables are different from ours. The entropy perturbation, e.g., is one of the two basic dynamical quantities in our context and neither velocity components nor metric functions do appear explicitly in the system of equations. An imperfect fluid description of scalar-tensor theories has been recently performed also in [21].

4.6. Coupled System of Equations for ε and S m

With relations (55) and (56) together with Formula (54), the set of Equations (40) and (48) becomes a closed system for ε and S m . The result, after transforming to the k-space ( Δ k 2 ), is
ε + 3 2 15 2 p ρ + 3 p ρ + A 1 ε a 3 2 + 12 p ρ 9 2 p 2 ρ 2 9 p ρ k 2 a 2 H 2 c ε 2 Ω x 1 + w x 1 + w B 1 ε a 2 + A 2 S m a + k 2 a 2 H 2 c ε 2 Ω m + B 2 S m a 2 = 0 ,
where
w = p ρ = p x ρ = w x Ω x ,
A 1 2 k 2 a 2 H 2 μ 3 2 1 + p ρ α ,
A 2 2 k 2 a 2 H 2 ν 3 2 1 + p ρ β ,
B 1 3 1 p ρ + 2 p ρ k 2 a 2 H 2 2 3 k 4 a 4 H 4 μ 9 4 1 + p ρ 1 p ρ α ,
B 2 3 1 p ρ + 2 p ρ k 2 a 2 H 2 2 3 k 4 a 4 H 4 ν 9 4 1 + p ρ 1 p ρ β
and
S m + 3 2 1 p ρ + 3 c ε 2 1 + w Ω m S m a + M 1 ( k , a ) S m a 2 + E 1 ( k , a ) ε a + E 2 ( k , a ) ε a 2 = 0 ,
with
M 1 ( k , a ) 1 1 + w k 2 a 2 H 2 c ε 2 Ω m + 3 c ε 2 3 p ρ + 1 2 + 3 p ρ 9 2 p ρ Ω m 2 3 k 4 a 4 H 4 ν ,
E 1 ( k , a ) = 3 1 + w c ε 2 p x ρ x 1 + w x 1 + w Ω x
and
E 2 ( k , a ) 1 1 + w k 2 a 2 H 2 c ε 2 + 3 c ε 2 p x ρ x 1 2 + 3 p ρ 9 2 p ρ 3 Ω m 1 + w p x ρ x a p x ρ x 1 + w x 1 + w Ω x 2 3 k 4 a 4 H 4 μ .
The system (57) and (63) describes the entire linear cosmological perturbation dynamics. It represents the central set of equations for our analysis. As we shall demonstrate in the following, its solution determines also the gravitational slip and the growth rate of matter fluctuations.

4.7. Anisotropic Pressure and Gravitational Potentials

The spatial part of Einstein’s equation can be written as
G ^ b a 1 3 δ b a G ^ m m = 8 π G T ^ b a 1 3 δ a b T ^ m m ,
where the left-hand side is determined by the difference ψ ϕ ,
G ^ b a 1 3 δ b a G ^ m m = 1 a 2 a b 1 3 δ a b Δ ψ ϕ ,
while the right-hand side of Equation (67) is related to the anisotropic pressure,
T ^ b a 1 3 δ a b T ^ m m = Π b a = 1 a 2 a b 1 3 δ a b Δ Π .
It follows that
ψ ϕ = 8 π G Π ,
the difference ψ ϕ is directly proportional to the anisotropic pressure. Notice that, in the present gauge, the quantities ψ and ϕ coincide with the Bardeen potentials Ψ and Φ , respectively. The relation between the comoving total density perturbation ε and the gravitational potential ψ is given by Equation (38). Combining Equation (38) with Equations (56) and (70), we obtain the relation
ϕ = 1 + 2 μ k 2 a 2 H 2 ψ 3 ν S m
between ψ and ϕ . In the limit μ = ν = 0 , equivalent to a vanishing anisotropic pressure, we recover ψ = ϕ . In the large-scale limit k 2 a 2 H 2 1 and assuming the S m contribution to be small, one has ψ ϕ as well, but, on smaller scales, both potentials may be different. Using Equation (38) to write relation (71) in the form
ϕ = 1 + 2 k 2 a 2 H 2 μ + ν S m ε ψ ,
it is obvious that knowledge of the relation between ψ and ϕ requires the solution of the entire coupled system of ε and S m .

4.8. Matter Perturbations

The final aim is to find the matter perturbations from the solution of the system for ε and S m . From the definition of S m (cf. Equation (42)) together with δ c = ε + l , it follows that the matter perturbations are determined by
δ m c = ε 1 + w + l 1 + w S m .
These are the perturbations with respect to the rest frame of the cosmic fluid as a whole, characterized by the quantity v. It is desirable, however, to calculate the matter-density perturbations in the matter rest frame, characterized by the matter velocity potential v m . Denoting this perturbation by δ m c m , the relation between both quantities is
δ m c m = δ m Θ v m = δ m c + Θ v v m .
Restricting ourselves to sub-horizon scales again, the matter density perturbations in the matter rest frame are
δ m c m = ε 1 + w x Ω x S m a 2 H 2 k 2 1 .
These matter perturbations are directly given by the solution of the coupled system (57) and (63) for ε and S m , respectively.
The entire setup so far is completely general and does not use any specific model for the x component. It is also valid for any homogeneous and isotropic background. Since any generalized or modified (compared with GR) gravitational theory can formally be rewritten as an effective Einsteinian theory, our formalism is expected to be applicable for a broad range of models. In the following, we apply this general scheme to a previously established scalar-tensor extension of the Λ CDM model, called e Φ Λ CDM model, which provides us with a non-standard background dynamics. This completes a preliminary study in [29], where the matter growth rate was obtained in a simplified manner on the basis of a rough approximation without solving the full dynamics given by Equations (57) and (63) and without including the heat flux.

5. e Φ Λ CDM Cosmology

5.1. Jordan–Brans–Dicke Theory

Our example originates from Jordan–Brans–Dicke (JBD) scalar-tensor theory [37,38,39], which was inspired by earlier ideas of Mach. Generally, the gravitational interaction in scalar-tensor theories is mediated by a scalar field in addition to the GR-type interaction through the metric tensor. Our aim is to find an equivalent GR description of such extended gravitational theory by mapping the additional (compared with GR) geometric degrees of freedom onto an effective fluid component. The energy-momentum tensor of this effective fluid will, in general, be of the structure of an imperfect fluid, i.e., it gives rise to effective anisotropic stresses and energy fluxes which are absent in a perfect-fluid based GR description. With an effective imperfect fluid description at hand, the perturbation analysis of the previous sections applies straightforwardly. In particular, it will be possible to calculate the growth rate of matter perturbations. What is needed, however, is a solution for the homogeneous and isotropic background which determines the coefficients in the set of first-order perturbation equations. The background is no longer that of the standard Λ CDM model, but it has to be the background of the extended gravitational theory. The first part of this section deals with the issue of how to obtain a suitable solution for a homogeneous and isotropic dynamics that deviates from the Λ CDM model.
Starting from a JBD type action (see, e.g., [40,41,42,43])
S ( g μ ν , Φ ) = 1 2 κ 2 d 4 x g Φ R ω ( Φ ) Φ Φ 2 U ( Φ ) + S m g μ ν ,
where S m is the matter part, the Jordan-frame gravitational field equations for scalar-tensor theories are
Φ R μ ν 1 2 g μ ν R = κ 2 T ( m ) μ ν + ω ( Φ ) Φ μ Φ ν Φ 1 2 g μ ν Φ 2 + μ ν Φ g μ ν Φ 1 2 g μ ν U ,
with κ 2 8 π G and
Φ = 1 2 ω ( Φ ) + 3 κ 2 T d ω ( Φ ) d Φ Φ 2 + Φ d U d Φ 2 U ,
where T is the trace of the matter energy-momentum tensor
T ( m ) μ ν = 2 g δ S m δ g μ ν .
The choice of the symbol T ( m ) μ ν indicates that we shall identify this quantity with the matter energy-momentum tensor of Section 2. In order to relate the entire formalism to the quantities of Section 2, we introduce a total effective energy-momentum tensor (cf. Equation (1)) T μ ν = T ( m ) μ ν + T ( x ) μ ν , where T ( x ) μ ν is an effective part that describes geometric “matter” and has the structure
T ( x ) μ ν 1 Φ 1 T ( m ) μ ν + 1 κ 2 Φ ω ( Φ ) Φ μ Φ ν Φ 1 2 g μ ν Φ 2 + μ ν Φ g μ ν Φ 1 2 g μ ν U .
Again, we have used the symbol T ( x ) μ ν , introduced in Section 2, to indicate identification with expression (80). The same is true for the sum T μ ν = T ( m ) μ ν + T ( x ) μ ν . With this definition, the field Equation (77) acquires the Einsteinian form,
R μ ν 1 2 g μ ν R = κ 2 T μ ν .
Having written JBD theory as an effective GR theory allows us to use techniques developed within the latter to make statements concerning the former as well. Our focus will be here on the background dynamics. In a spatially flat homogeneous and isotropic model, one may associate a Hubble rate to the scalar-tensor dynamics [43]:
H 2 = κ 2 3 ρ m Φ + 1 3 Φ 1 2 ω ( Φ ) Φ Φ t 2 3 H Φ t + 1 2 U ,
where H = 1 a d a d t is the Hubble rate of the Jordan frame and a is the Jordan-frame scale factor. Furthermore, assuming the matter component to be pressureless, one has
d H d t = κ 2 2 ρ m Φ 1 2 Φ ω ( Φ ) Φ Φ t 2 H Φ t + d 2 Φ d t 2 ,
and
d 2 Φ d t 2 + 3 H d Φ d t = 1 2 ω + 3 κ 2 ρ m d ω d Φ d Φ d t 2 Φ d U d Φ + 2 U ,
as well as the matter conservation (14).
The structure of Equation (82) suggests the definition
ρ x ( Φ ) = ρ m 3 1 Φ 1 + 1 3 κ 2 Φ 1 2 ω ( Φ ) Φ Φ t 2 3 H Φ t + 1 2 U
such that Equation (82) takes the form of an effective Friedmann equation 3 H 2 = κ 2 ρ m + ρ x ( Φ ) . Obviously, ρ x ( Φ ) is determined by the dynamics of the scalar field, this is indicated by the additional subscript ( Φ ) .

5.2. JBD Inspired Effective Background Model

Here, we recall the basic elements of the previously established homogeneous and isotropic scalar-tensor extension of the Λ CDM model [28]. Its main characteristic is an explicit analytical expression for the Hubble rate in which deviations from the standard model are described by a single constant parameter which also governs the effective scalar-field dynamics. The background relations are found through a specific solution of the effective fluid dynamics in the Einstein frame which subsequently is converted into the Jordan frame via a conformal transformation. To make our presentation self-contained, we summarize the main steps of the corresponding derivation. The starting point for the background analysis is the Einstein frame. The reason for this apparent detour is that it is the Einstein frame in which it is possible to obtain a solution of the dynamics.

5.2.1. General Einstein-Frame Dynamics

Generally, the Einstein-frame dynamics with a metric g ˜ μ ν follows from the (so far considered) Jordan-frame dynamics with a metric g μ ν through the conformal transformation
g μ ν = 1 Φ g ˜ μ ν = e 2 b ( φ ) g ˜ μ ν
together with a redefinition of the potential term,
V ( φ ) = U ( Φ ) 2 κ 2 Φ 2 .
Moreover, one has
1 4 Φ 2 d Φ d φ 2 = d b d φ 2 = κ 2 4 ω ( Φ ) + 6 .
For an RW metric and assuming u ( m ) μ = u μ , the general relations simplify considerably. With the Einstein-frame scale factor a ˜ and the Einstein-frame time coordinate t ˜ the Einstein-frame Hubble rate H ˜ = 1 a ˜ d a ˜ d t ˜ becomes
H ˜ 2 = κ 2 3 ρ ˜ m + 1 2 d φ d t ˜ 2 + V ˜ ,
where ρ ˜ m is the matter-energy density in the Einstein frame. Quantities with a tilde have their meaning in the Einstein frame, those without a tilde refer to the Jordan frame. Different from the Jordan-frame dynamics, the matter part is not separately conserved but obeys the equation
d ρ ˜ m d t ˜ + 3 H ˜ ρ ˜ m = d φ d t ˜ d b d φ ρ ˜ m
through which it couples to the following scalar-field dynamics:
d 2 φ d t ˜ 2 + 3 H ˜ d φ d t ˜ + V ˜ , φ = d b d φ ρ ˜ m .
Here, V ˜ , φ = V ˜ φ . The time coordinate t, the scale factor a and the matter energy density ρ m of the Jordan frame are related to their Einstein-frame counterparts by
d t = e b d t ˜ , a = e b a ˜ and ρ m = e 4 b ρ ˜ m ,
respectively.

5.2.2. Interacting Fluid Approach in Einstein-Frame Dynamics

We associate an effective energy density ρ ˜ φ and an effective pressure p ˜ φ to the scalar field by
ρ ˜ φ = 1 2 d φ d t ˜ 2 + V ˜ and p ˜ φ = 1 2 d φ d t ˜ 2 V ˜ ,
respectively. Equation (91) then takes the form
d ρ ˜ φ d t ˜ + 3 H ˜ 1 + w ˜ ρ ˜ φ = Q d φ d t ˜ d b d φ ρ ˜ m ,
where w ˜ = p ˜ φ ρ ˜ φ is the Einstein-frame EoS parameter for the scalar field. Equation (90) can be written as
ρ ˜ m = ρ ˜ m 0 a ˜ 3 f ( a ˜ ) d ρ ˜ m d t ˜ + 3 H ˜ ρ ˜ m = ρ ˜ m f d f d t ˜ , f = e b ( φ ) ,
where the function f encodes the effects of an interaction between matter and (Einstein-frame) scalar field. Assuming a power-law behavior of the interaction function f a ˜ ,
f a ˜ = a ˜ 3 m ,
the explicit solution of Equation (94) for a constant EoS parameter then is
ρ ˜ φ = ρ ˜ φ 0 a ˜ 3 1 + w ˜ ρ ˜ m 0 a ˜ 3 1 + w ˜ a ˜ 0 a ˜ d a ˜ d f d a ˜ a ˜ 3 w ˜ .
To obtain the explicit solutions (95) with (96) and (97) was the main motivation for considering the Einstein frame. This implies the expression
Q = 3 m H ˜ ρ ˜ m
for the interaction term Q, defined in (94). In the simple case of a linear dependence
b = b ( φ ) = K φ , K = κ 2 4 ω + 6 ,
it follows that [28]
Φ = e 2 K φ = e 2 b = a ˜ 6 m = a 6 m 1 + 3 m .
with Formula (100), we found an explicit expression for the scalar-field variable without having solved the basic scalar-field equation. We have solved the system of energy-balance Equations (90) and (94) under the assumptions (96) and (99). This solution implies an explicit scale-factor dependence of Φ (or φ ) which does not necessarily have to be a solution of the original scalar-field Equation (78). Instead, it obeys an alternative effective second-order equation with an alternative effective potential that does not coincide with U [28]. The point here is that the dynamics on the level of the fluid energy densities do not require the exact solution of the scalar-field Equation (78). On the other hand, the specific features of our fluid dynamics imply the existence of the effective scalar field Φ of the form of (100). Under such condition, the dynamics for w ˜ = 1 are explicitly solved, which results in the expression [28]
H ˜ 2 H ˜ 0 2 = e K φ Ω ˜ m 0 a ˜ 3 1 m + 1 Ω ˜ m 0 1 m , Ω ˜ m 0 = 8 π G ρ ˜ m 0 3 H ˜ 0 2
for the Einstein-frame Hubble rate. Having used the Einstein frame to explicitly solve the background dynamics, we now turn back to the Jordan-frame dynamics.

5.2.3. Effective Hubble Rate

The transformation to the Jordan-frame Hubble rate via
H = e b 1 + 3 m H ˜
leads to the explicit effective Hubble rate [28]
H 2 H 0 2 = A Ω m 0 a 3 Φ + 1 A Ω m 0 Φ ,
where
Ω m 0 = Ω ˜ m 0 1 + 3 m 2 , A 1 + 3 m 2 1 m
and Φ is explicitly given in terms of the scale factor by Relation (100). Formula (103) represents an explicit analytic solution for the Hubble rate of our e Φ Λ CDM model. The scalar Φ modifies the cosmological dynamics compared with the GR based Λ CDM model. For Φ = 1 , equivalent to m = 0 , we recover the Λ CDM model. For any Φ 1 , equivalent to m 0 , the expression (103) represents a testable, alternative model with small ( | m | 1 ) deviations from the Λ CDM model. It is Formula (103) which allows us to apply the general perturbation dynamics of the previous sections to a specific class of non-standard models. To motivate the result, (103) has been the main purpose of this section. The dependence of the Hubble rate (103) on the scalar Φ changes the relative contributions of matter and the dark-energy (DE) equivalent compared with the Λ CDM model. Deviations from the Λ CDM model are entirely encoded in the constant parameter m with | m | 1 . The fractional matter contribution is
Ω m = ρ m ρ = Ω m 0 a 3 A Φ 1 Ω m 0 a 3 + 1 A Ω m 0 Φ ,
the geometric “matter” part contributes with Ω x = 1 Ω m . For different values of m, the fractional abundances Ω m and Ω x are visualized in Figure 1. Postulating a conservation equation ρ ˙ x + 3 H 1 + w x ρ x = 0 , where ρ x ρ ρ m now replaces Expression (85), this corresponds to an effective, time-varying EoS parameter w x of the geometric DE,
w x ( a ) = 1 + 2 m 1 + 3 m 1 A Ω m 0 Φ + Ω m 0 a 3 1 + m 1 + 3 m A Φ 1 1 1 A Ω m 0 Φ + Ω m 0 a 3 A Φ 1 1 .
For m = 0 , it reduces to the Λ CDM value w x = 1 . At a high redshift, one has
w x 1 + 1 + m 1 + 3 m A Φ 1 1 A Φ 1 1 ( a 1 ) .
This value may be close to zero, i.e., the geometric DE may mimic dust in this limit, but the effective energy density ρ x will be negative for m > 0 . It crosses ρ x = 0 in the redshift range 10 z 4 for the values of m chosen in Figure 1. This behavior reflects that fact that the x-component is very different from a conventional fluid. The total EoS is well behaved throughout as is demonstrated in Figure 2. At the present time, the effective EoS parameter is
w x = 1 + 2 m 1 + 3 m + 3 m Ω m 0 1 Ω m 0 ( a = 1 ) .
For small | m | , this remains in the vicinity of w x = 1 . In the far future, w x approaches
w x 1 + 2 m 1 + 3 m ( a 1 ) .
From a statistical analysis using Supernovae data, data from the differential age of old galaxies that have evolved passively and baryon acoustic oscillations, we found a best-fit value [28] of m = 0.004 0.011 ( 1 σ ) 0.017 ( 2 σ ) + 0.011 ( 1 σ ) + 0.017 ( 2 σ ) . This is compatible with the Λ CDM model, but it leaves room for small deviations. Even a very small non-vanishing value of | m | will modify the standard scenario of structure formation. The quantitative analysis to follow will rely on the use of the effective Hubble rate (103) for the background coefficients of the perturbation equations.

6. Growth of Matter Perturbations

Now, we combine the general imperfect-fluid perturbation dynamics, established in Section 4, with the e Φ Λ CDM background model of the previous section. The variable of interest, the fractional matter perturbation δ m c m in Formula (75), is obtained via the solution of the coupled system (57) and (63) for ε and S m , respectively. Since on sub-horizon scales gauge issues are less important, we shall omit from now on the superscript c m and denote this quantity simply by δ m . It is convenient to use also the linear growth rate f defined by f = d ln δ m / d ln a . In most cases, observational data are provided for the combination f σ 8 where σ 8 is the root-mean-square mass fluctuation in spheres with radius 8 h 1 Mpc [44]. In the linear regime, one has [23]
σ 8 ( z ) = δ m ( z ) δ ( z = 0 ) σ 8 ( z = 0 )
and
f σ 8 ( z ) = ( 1 + z ) σ 8 ( z = 0 ) δ m ( z = 0 ) d d z δ m ( z ) .
For the matter distribution variance today, we assume σ 8 ( z = 0 ) = 0.8 , which is compatible with standard fiducial cosmology as determined by the Planck satellite [12]. This value is biased with respect to the variance in the distribution of galaxies but the combination f σ 8 is independent from the bias factor [44]. The data of RSD measurements of f σ 8 used in our analysis are listed in Table 1. In the following, we study the influence of some of the model parameters on the growth of matter perturbations. Differences to the Λ CDM model may already occur if there are isotropic pressure perturbations, described by a non-vanishing sound speed parameter. Our analysis shows that values of the order of c ε 2 > 10 4 are necessary to produce a noticeable impact (Figure 3).
Figure 4 shows the fractional matter energy density perturbations δ m as a function of the scale factor if only the parameter m is modified (different values of m for α = β = μ = ν = 0 ) with respect to the Λ CDM model ( m = 0 ). The dependence of the growth rate f σ 8 on the redshift for these values together with the data points is shown in Figure 5. Positive values of m enhance the matter growth, negative values diminish it. Figure 6 visualizes the scale-factor dependence of the fractional matter density perturbations for different values of the anisotropic stress parameter μ for m = α = β = ν = 0 . The corresponding growth rate f σ 8 for different values of μ is shown in Figure 7. Deviations from the Λ CDM model require values of μ of the order of | 10 10 | . Finally, we demonstrate the impact of a non-vanishing heat flow, represented by the parameter α , on the growth rate (see Figure 8). The heat flux basically causes a shift of the Λ CDM curve (upward for α > 0 , downward for α < 0 ).

7. Conclusions

We have established a general phenomenological scheme for implementing (effective) non-equilibrium effects in a fluid description of the cosmological dark sector. This comprises both “true" dissipative effects within Einstein’s GR, assuming that the cosmic substratum behaves less simply than taken for granted in the usually applied perfect-fluid approach, and those features which originate from a re-interpretation of geometric terms in modified gravitational theories in terms of effective fluid quantities. Using a combination of the fluid conservation equations with the Raychaudhuri equation for the expansion scalar, the complete first-order, scalar perturbation dynamics have been reduced to a manifestly gauge-invariant coupled system of two second-order differential equations for the total and the relative energy–density perturbations. We clarified how a heat flux (effective or “true”) modifies the Poisson-type equation for the gravitational potential. A characteristic feature of our approach consists in the introduction of phenomenological parameters with the purpose to make the fluid dynamical equations a closed system. The relevant relations by which these coefficients are introduced are inspired by the manner that the speed of sound is conventionally introduced on a phenomenological basis. In a standard perfect-fluid description of the Universe, such relation between the perturbations of (isotropic) pressure and energy density is required to close the system of perturbation equations. Here, we are generalizing this procedure by adding relations of a similar type which take into account anisotropic pressure and heat flux. As in the case of a phenomenologically introduced sound speed, a derivation from an underlying fundamental theory is left open. Even for the sound speed parameter, an analytic microscopic justification does exist only in special cases. This obvious shortcoming of the phenomenological theory is the price to pay for obtaining a robust, transparent and in large part analytical description of the inhomogeneous dynamics.
Our analysis is preliminary since it so far gives only a very rough account of the relevance of different (effective) dissipative phenomena on the growth of matter inhomogeneities during the cosmic history. At this point, also in view of the large error bars of the data, only order-of-magnitude estimates are possible. Our approach allows for deviations from the standard model as long as these deviations are small. Additional information is also needed to decide whether deviations from the standard model which are likely to be tiny, can be attributed to deviations from Einstein’s GR or to a “real” non-equilibrium nature of the cosmic substratum within GR.
In a general scalar-tensor theory, all the effective fluid quantities energy density, isotropic pressure, anisotropic pressure and heat flux are given in terms of the perturbations of the scalar field [16,21]. However, this dependence, while exact, is rather involved and, to obtain observationally relevant quantities, it needs numerical implementation at a much earlier stage compared with the scheme presented here. Of course, a final justification of this scheme will require a sound microscopic foundation, a problem we hope to deal with for specific cases in future work.

Author Contributions

The authors contributed as a team to this work. Observations, numerical and statistical analysis: H.E.S.V. and W.C.A., original draft: W.Z.

Funding

This research was funded by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), FUNDAÇÃO ESTADUAL DE AMPARO À PESQUISA DO ESTADO DO ESPÍRITO SANTO (FAPES) and COORDENAÇÃO DE APERFEIÇOAMENTO DE PESSOAL DE NÍVEL SUPERIOR (CAPES).

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the writing of the manuscript, and in the decision to publish the results.

References

  1. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef]
  2. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of Ω and Λ from 42 High-Redshift Supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  3. Bahcall, N.A.; Ostriker, J.P.; Perlmutter, J.P.; Steinhardt, P.J. The cosmic triangle: Revealing the state of the universe. Science 1999, 284, 1481–1488. [Google Scholar] [CrossRef]
  4. Amendola, L.; Tsujikawa, S. Dark Energy: Theory and Observations; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  5. Ellis, G.F.R.; Maartens, R.; Maccallum, M.A.H. Relativistic Cosmology; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  6. Martin, J. Everything you always wanted to know about the cosmological constant problem (but were afraid to ask). C. R. Phys. 2012, 13, 566–665. [Google Scholar] [CrossRef] [Green Version]
  7. Zlatev, I.; Wang, L.M.; Steinhardt, P.J. Quintessence, cosmic coincidence, and the cosmological constant. Phys. Rev. Lett. 1999, 82, 896. [Google Scholar] [CrossRef]
  8. Steinhardt, P.J.; Wang, L.M.; Zlatev, I. Cosmological tracking solutions. Phys. Rev. D 1999, 59, 123504. [Google Scholar] [CrossRef] [Green Version]
  9. Malquarti, M.; Copeland, E.J.; Liddle, A.R. k-essence and the coincidence problem. Phys. Rev. D 2003, 68, 023512. [Google Scholar] [CrossRef]
  10. Barreira, A.; Avelino, P.P. Anthropic versus cosmological solutions to the coincidence problem. Phys. Rev. D 2011, 83, 103001. [Google Scholar] [CrossRef]
  11. Velten, H.E.S.; vom Marttens, R.F.; Zimdahl, W. Aspects of the cosmological “coincidence problem”. Eur. Phys. J. C 2014, 74, 3160. [Google Scholar] [CrossRef]
  12. Akrami, Y.; et al. [Planck Collaboration] Planck 2018 results. I. Overview and the cosmological legacy of Planck. arXiv, 2018; arXiv:1807.06205. [Google Scholar]
  13. Abbott, T.M.C.; et al. [DES Collaboration] Dark Energy Survey Year 1 Results: Constraints on Extended Cosmological Models from Galaxy Clustering and Weak Lensing. arXiv, 2018; arXiv:1810.02499. [Google Scholar]
  14. Bull, P.; Camera, S.; Kelley, K.; Padmanabhan, H.; Pritchard, J.; Raccanelli, A.; Riemer-Sørensen, S.; Shao, L.; Andrianomena, S.; Athanassoula, E.; et al. Fundamental Physics with the Square Kilometre Array. arXiv, 2018; arXiv:1810.02680. [Google Scholar]
  15. Abbott, B.P.; et al. [LIGO Scientific Collaboration and Virgo Collaboration] GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Madsen, M.S. Scalar fields in curved spacetimes. Class. Quantum Grav. 2017, 5, 627. [Google Scholar] [CrossRef]
  17. Pimentel, L.O. Energy-momentum tensor in the general scalar-tensor theory. Class. Quantum Grav. 1989, 6, L263. [Google Scholar] [CrossRef]
  18. Battye, R.A.; Pearson, J.A. Effective action approach to cosmological perturbations in dark energy and modified gravity. J. Cosmol. Astropart. Phys. 2012, 2012, 019. [Google Scholar] [CrossRef]
  19. Battye, R.A.; Pearson, J.A. Parametrizing dark sector perturbations via equations of state. Phys. Rev. D 2013, 88, 061301(R). [Google Scholar] [CrossRef]
  20. Battye, R.A.; Bolliet, B.; Pearson, J.A. f(R) gravity as a dark energy fluid. Phys. Rev. D 2016, 93, 044026. [Google Scholar] [CrossRef]
  21. Faraoni, V.; Coté, J. Imperfect fluid description of modified gravities. Phys. Rev. D 2018, 98, 084019. [Google Scholar] [CrossRef]
  22. Sawicki, I.; Saltas, I.D.; Amendola, L.; Kunz, M. Consistent perturbations in an imperfect fluid. J. Cosmol. Astropart. Phys. 2013, 2013, 004. [Google Scholar] [CrossRef]
  23. Nesseris, S.; Perivolaropoulos, S. Testing LCDM with the Growth Function δ(a): Current Constraints. Phys. Rev. D 2008, 77, 023504. [Google Scholar] [CrossRef]
  24. Basilakos, S.; Pouri, A. The growth index of matter perturbations and modified gravity. Mon. Not. R. Astron. Soc. 2008, 423, 3761–3767. [Google Scholar] [CrossRef]
  25. Huterer, D.; Kirkby, D.; Bean, R.; Connolly, A.; Dawson, K.; Dodelson, S.; Evrard, A.; Jain, B.; Jarvis, M.; Linder, E.; et al. Growth of Cosmic Structure: Probing Dark Energy Beyond Expansion. Astropart. Phys. 2015, 63, 23–41. [Google Scholar] [CrossRef]
  26. Nesseris, S.; Sapone, D. Accuracy of the growth index in the presence of dark energy perturbations. Phys. Rev. D 2015, 92, 023013. [Google Scholar] [CrossRef]
  27. Alam, S.; Ata, M.; Bailey, S.; Beutler, F.; Bizyaev, D.; Blazek, J.A.; Bolton, A.S.; Brownstein, J.R.; Burden, A.; Chuang, C.-H.; et al. The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmological analysis of the DR12 galaxy sample. Mon. Not. R. Astron. Soc. 2017, 470, 2617–2652. [Google Scholar] [CrossRef]
  28. Algoner, W.C.; Velten, H.E.S.; Zimdahl, W. Scalar-tensor extension of the ΛCDM model. J. Cosmol. Astropart. Phys. 2016, 2016, 034. [Google Scholar] [CrossRef]
  29. Zimdahl, W.; Velten, H.E.S.; Algoner, W.C. Matter growth in extended ΛCDM cosmology. arXiv, 2017; arXiv:1706.06143. [Google Scholar] [CrossRef]
  30. Hipólito-Ricaldi, W.S.; Velten, H.E.S.; Zimdahl, W. Non-adiabatic dark fluid cosmology. J. Cosmol. Astropart. Phys. 2009, 2009, 016. [Google Scholar] [CrossRef]
  31. Hipólito-Ricaldi, W.S.; Velten, H.E.S.; Zimdahl, W. Viscous dark fluid universe. Phys. Rev. D 2010, 82, 063507. [Google Scholar] [CrossRef]
  32. del Campo, S.; Fabris, J.C.; Herrera, J.C.; Zimdahl, W. Cosmology with Ricci dark energy. Phys. Rev. D 2013, 87, 123002. [Google Scholar] [CrossRef]
  33. Romero Fuño, A.; Hipólito-Ricaldi, W.S.; Zimdahl, W. Matter perturbations in scaling cosmology. Mon. Not. R. Astron. Soc. 2016, 457, 2958–2967. [Google Scholar] [CrossRef] [Green Version]
  34. Kunz, M.; Sapone, M. Dark Energy versus Modified Gravity. Phys. Rev. Lett. 2007, 98, 121301. [Google Scholar] [CrossRef] [PubMed]
  35. Cardona, W.; Hollenstein, L.; Kunz, M. The traces of anisotropic dark energy in light of Planck. J. Cosmol. Astropart. Phys. 2014, 2014, 032. [Google Scholar] [CrossRef]
  36. Blas, D.; Floerchinger, S.; Garny, M.; Tetradis, N.; Wiedemann, U.A. Large scale structure from viscous dark matter. J. Cosmol. Astropart. Phys. 2015, 2015, 049. [Google Scholar] [CrossRef]
  37. Jordan, P. Zum gegenwärtigen Stand der Diracschen kosmologischen Hypothesen. Z. Physik 1959, 157, 112–121. [Google Scholar] [CrossRef]
  38. Brans, C.; Dicke, R.H. Mach’s principle and a relativistic theory of gravitation. Phys. Rev. 1961, 124, 925. [Google Scholar] [CrossRef]
  39. Dicke, R.H. Mach’s principle and invariance under transformation of units. Phys. Rev. 1962, 125, 2163. [Google Scholar] [CrossRef]
  40. Agarwal, N.; Bean, R. The Dynamical viability of scalar-tensor gravity theories. Class. Quantum Grav. 2008, 25, 165001. [Google Scholar] [CrossRef]
  41. Batista, C.E.M.; Zimdahl, W. Power-law solutions and accelerated expansion in scalar-tensor theories. Phys. Rev. D 2010, 82, 023527. [Google Scholar] [CrossRef]
  42. Clifton, T.; Ferreira, P.; Padilla, A.; Skordis, C. Modified Gravity and Cosmology. Phys. Rep. 2012, 513, 1–189. [Google Scholar] [CrossRef]
  43. Chiba, T.; Yamaguchi, M. Conformal-Frame (In)dependence of Cosmological Observations in Scalar-Tensor Theory. J. Cosmol. Astropart. Phys. 2013, 2013, 040. [Google Scholar] [CrossRef]
  44. Song, Y.-S.; Percival, W.J. Reconstructing the history of structure formation using redshift distortions. J. Cosmol. Astropart. Phys. 2009, 2009, 009. [Google Scholar] [CrossRef]
  45. Nesseris, S.; Pantazis, G.; Perivolaropoulos, L. Tension and constraints on modified gravity parametrizations of Geff(z) from growth rate and Planck data. Phys. Rev. D 2017, 96, 023542. [Google Scholar] [CrossRef] [Green Version]
  46. Gil-Marín, H.; Guy, J.; Zarrouk, P.; Burtin, E.; Chuang, C.-H.; Percival, W.J.; Ross, A.J.; Ruggeri, R.; Tojerio, R.; Zhao, G.-B.; et al. The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: Structure growth rate measurement from the anisotropic quasar power spectrum in the redshift range 0.8<z<2.2. Mon. Not. R. Astron. Soc. 2018, 477, 1604–1638. [Google Scholar]
  47. Hou, J.; Sánchez, A.G.; Scoccimarro, R.; Salazar-Albornoz, S.; Burtin, E.; Gil-Marín, H.; Percival, W.J.; Ruggeri, R.; Zarrouk, P.; Zhao, G.-B.; et al. The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: Anisotropic clustering analysis in configuration-space. Mon. Not. R. Astron. Soc. 2018, arXiv:1801.02656480, 2521–2534. [Google Scholar] [CrossRef]
  48. Zhao, G.-B.; Wang, Y.; Saito, S.; Gil-Marín, H.; Percival, W.J.; Wang, D.; Chuang, C.-H.; Ruggeri, R.; Mueller, E.-M.; Zhu, F.; et al. The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: A tomographic measurement of cosmic structure growth and expansion rate based on optimal redshift weights. Mon. Not. R. Astron. Soc. 2018, 482, 3497–3513. [Google Scholar] [CrossRef]
Figure 1. Matter fraction Ω m and geometric energy fraction Ω x for negative (left) and positive (right) values of m.
Figure 1. Matter fraction Ω m and geometric energy fraction Ω x for negative (left) and positive (right) values of m.
Universe 05 00068 g001
Figure 2. Total EoS parameter w = p ρ for various negative (left) and positive (right) values of m.
Figure 2. Total EoS parameter w = p ρ for various negative (left) and positive (right) values of m.
Universe 05 00068 g002
Figure 3. Dependence of f σ 8 ( z ) on z for a non-vanishing sound-speed parameter ( m = α = β = μ = ν = 0 ).
Figure 3. Dependence of f σ 8 ( z ) on z for a non-vanishing sound-speed parameter ( m = α = β = μ = ν = 0 ).
Universe 05 00068 g003
Figure 4. Matter growth for various values of the parameter m with α = β = μ = ν = c ε 2 = 0 .
Figure 4. Matter growth for various values of the parameter m with α = β = μ = ν = c ε 2 = 0 .
Universe 05 00068 g004
Figure 5. Dependence of f σ 8 ( z ) on z for different values of the parameter m with α = β = μ = ν = c ε 2 = 0 .
Figure 5. Dependence of f σ 8 ( z ) on z for different values of the parameter m with α = β = μ = ν = c ε 2 = 0 .
Universe 05 00068 g005
Figure 6. Matter growth in the presence of anisotropic stresses ( m = α = β = ν = c ε 2 = 0 ).
Figure 6. Matter growth in the presence of anisotropic stresses ( m = α = β = ν = c ε 2 = 0 ).
Universe 05 00068 g006
Figure 7. Dependence of f σ 8 ( z ) on z in the presence of anisotropic stresses ( m = α = β = ν = c ε 2 = 0 ).
Figure 7. Dependence of f σ 8 ( z ) on z in the presence of anisotropic stresses ( m = α = β = ν = c ε 2 = 0 ).
Universe 05 00068 g007
Figure 8. Dependence of f σ 8 ( z ) on z in the presence of heat fluxes ( m = β = μ = ν = c ε 2 = 0 ).
Figure 8. Dependence of f σ 8 ( z ) on z in the presence of heat fluxes ( m = β = μ = ν = c ε 2 = 0 ).
Universe 05 00068 g008
Table 1. Data points for our analysis. The first 18 entries represent the “Gold” set of [45]. The entries 19 and 20 are taken from [46,47], respectively. The last three values have been reported in [48].
Table 1. Data points for our analysis. The first 18 entries represent the “Gold” set of [45]. The entries 19 and 20 are taken from [46,47], respectively. The last three values have been reported in [48].
IndexData setz f σ 8 ( z ) YearNotes
16dFGS + SnIa0.020.428 ± 0.04652016( Ω 0 m , h, σ 8 ) = (0.3, 0.683, 0.8)
2SnIa + IRAS0.020.398 ± 0.0652011( Ω 0 m , Ω K ) = (0.3, 0)
32MASS0.020.314 ± 0.0482010( Ω 0 m , Ω K ) = (0.266, 0)
4SDSS-veloc0.100.370 ± 0.1302015( Ω 0 m , Ω K ) = (0.3, 0)
5SDSS-MGS0.150.490 ± 0.1452014( Ω 0 m , h, σ 8 ) = (0.31, 0.67, 0.83)
62dFGRS0.170.510 ± 0.0602009( Ω 0 m , Ω K ) = (0.3, 0)
7GAMA0.180.360 ± 0.0902013( Ω 0 m , Ω K ) = (0.27, 0)
8GAMA0.380.440 ± 0.0602013
9SDSS-LRG-2000.250.3512 ± 0.05832011( Ω 0 m , Ω K ) = (0.25, 0)
10SDSS-LRG-2000.370.4602 ± 0.03782011
11BOSS-LOWZ0.320.384 ± 0.0952013( Ω 0 m , Ω K ) = (0.274, 0)
12SDSS-CMASS0.590.488 ± 0.0602013( Ω 0 m , h, σ 8 ) = (0.307115, 0.6777, 0.8288)
13WiggleZ0.440.413 ± 0.0802012( Ω 0 m , h) = (0.27, 0.71)
14WiggleZ0.600.390 ± 0.0632012
15WiggleZ0.730.437 ± 0.0722012
16Vipers PDR-20.600.550 ± 0.1202016( Ω 0 m , Ω b ) = (0.3, 0.045)
17Vipers PDR-20.860.400 ± 0.1102016
18FastSound1.400.482 ± 0.1162015( Ω 0 m , Ω K ) = (0.270, 0)
19SDSS-IV1.520.420 ± 0.0762018( Ω 0 m , Ω b h 2 , σ 8 ) = (0.26479, 0.02258, 0.8)
20SDSS-IV1.520.396 ± 0.0792018( Ω 0 m , Ω b h 2 , σ 8 ) = (0.31, 0.022, 0.8225)
21SDSS-IV1.230.385 ± 0.0992018( Ω 0 m , σ 8 ) = (0.31, 0.8)
22SDSS-IV1.5260.342 ± 0.0702018
23SDSS-IV1.9440.364 ± 0.1062018

Share and Cite

MDPI and ACS Style

Zimdahl, W.; Velten, H.E.S.; Algoner, W.C. Matter Growth in Imperfect Fluid Cosmology. Universe 2019, 5, 68. https://doi.org/10.3390/universe5030068

AMA Style

Zimdahl W, Velten HES, Algoner WC. Matter Growth in Imperfect Fluid Cosmology. Universe. 2019; 5(3):68. https://doi.org/10.3390/universe5030068

Chicago/Turabian Style

Zimdahl, Winfried, Hermano E.S. Velten, and William C. Algoner. 2019. "Matter Growth in Imperfect Fluid Cosmology" Universe 5, no. 3: 68. https://doi.org/10.3390/universe5030068

APA Style

Zimdahl, W., Velten, H. E. S., & Algoner, W. C. (2019). Matter Growth in Imperfect Fluid Cosmology. Universe, 5(3), 68. https://doi.org/10.3390/universe5030068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop