Ricci cosmology in light of astronomical data

Recently, a new cosmological framework, dubbed Ricci Cosmology, has been proposed. Such a framework has emerged from the study of relativistic dynamics of fluids out of equilibrium in a curved background and is characterised by the presence of deviations from the equilibrium pressure in the energy-momentum tensor which are due to linear terms in the Ricci scalar and the Ricci tensor. The coefficients in front of such terms are called the second order transport coefficients and they parametrise the fluid response to the pressure terms arising from the spacetime curvature. Under the preliminary assumption that the second order transport coefficients are constant, we find the simplest solution of Ricci cosmology in which the presence of pressure terms causes a departure from the perfect fluid redshift scaling for matter components filling the Universe. In order to test the viability of this solution, we make four different ans\"{a}tze on the transport coefficients, giving rise to four different cases of our model. On the physical ground of the second law of thermodynamics for fluids with non-equilibrium pressure, we find some theoretical bounds (priors) on the parameters of the models. Our main concern is then the check of each of the case against the standard set of cosmological data in order to obtain the observational bounds on the second order transport coefficients. We find those bounds also realising that Ricci cosmology model is compatible with $\Lambda$CDM cosmology for all the ans\"{a}tze.


I. INTRODUCTION
The discovery of the late-time accelerated expansion of the Universe was made in the late 1990s, observing fainter than previously predicted Type Ia Supernovae (SnIa) [1,2]. Since then, more and more precise observations, such as Baryon Acoustic Oscillation (BAO) [3] and Cosmic Microwave Background (CMB) anisotropies [4], have pointed in the same direction making the explanation of these observations so compelling that it has become one of the central issues of modern cosmology. The most successful model in explaining such feature of our Universe and in fitting the available data (SnIa, BAO, CMB, Big Bang Nucleosynthesis (BBN) and Large Scale Structure (LSS)) is the Standard Cosmological Model, better known as ΛCDM [5,6], in which the Universe is well described below the Planck scale by General Relativity (GR) and at scales of 100 Mpc and larger each component of the Universe can be described as a perfect fluid. The Universe is filled with baryonic matter, radiation, Cold Dark Matter (CDM) and the cosmological constant Λ, which is responsible for the late-time accelerated expansion [7,8].
Despite the success and the simplicity of the ΛCDM model, an explanation for the physical nature of CDM and Λ is still missing and related to the latter, there are two big unsolved issues: it is still unclear why Λ is so small compared to vacuum energy predicted by Quantum Field Theory (QFT) (fine-tuning problem) and why it a roberto.caroli@usz.edu.pl b mariusz.dabrowski@usz.edu.pl c vincenzo.salzano@usz.edu.pl becomes important only at very late times (coincidence problem) [9,10].
More recently, as soon as we entered in the era of precision cosmology, tensions in the data started to appear. The most important example is the Hubble tension [11]. On one hand, from the CMB anisotropy measurements from Planck [4], for the value of the Hubble parameter H 0 we have H 0 = (67.36 ± 0.54) km s −1 Mpc −1 , by assuming ΛCDM as the fiducial cosmological model. On the other hand, the local measurements of the same parameter point towards greater values, among which from the Hubble Space Telescope we have H 0 = (74.03 ± 1.42) km s −1 Mpc −1 [12], with a 4.4σ discrepancy with the early-universe value above from Planck.
Currently, this discrepancy does not seem to be due to systematic effects in either early-time or late-time measurements and may point to new physics beyond ΛCDM [13,14].
Because of these issues, in the last decades, a great variety of models, which try to explain the late-time accelerated expansion of the Universe without resorting to the cosmological constant Λ, have been proposed [15,16] which also try to address the Hubble tension problem [17]. For reviews on the experimental status of some of these theories beyond ΛCDM, see Refs. [18,19].
In most of the alternative models, the matter components filling the Universe are usually treated as perfect fluids. A slightly unconventional path to explain latetime accelerated expansion, which we consider in this paper is modifying the usual matter content by adding dissipative or viscous terms in the Energy-Momentum Tensor (EMT) describing perfect fluids. For recent reviews, see Refs. [20,21].
Most of the best explored viscous models so far are those involving bulk viscosity where a term proportional to the Hubble function is added to one or more matter components of the Universe, giving rise to a modified effective pressure for them. These models may differ from one another in the assumptions made on the bulk viscosity coefficient and/or on the matter content of the Universe. For this reason, the effective pressure arising from viscosity can be important either in the early stages of the Universe or at the present time.
In the first case, fall all those models in which the inflation is driven by viscosity [22][23][24][25][26][27][28][29][30][31][32]. A relevant feature of these models is that by admitting viscosity, one is able to avoid singularity in the Universe as well as to grow the Universe in size in oscillatory sequences.
In the second case, the models that can describe the late-time accelerated expansion of the Universe are considered [33][34][35][36][37][38][39][40][41][42][43][44][45][46]. These models have revealed to be successful in reproducing the background evolution of the Universe at late time while attempting to solve the problems of ΛCDM, unify or replace the dark components of the Universe and be compatible with the large structure formation. It is still debated whether viscous terms in the pressure may or may not contribute to relieve the Hubble tension [47][48][49].
Recently, with the purpose of describing the inflationary epoch of the Universe, a new framework, dubbed Ricci Cosmology, has been proposed [50]. This framework involves second order corrections in gradients of the metric tensor g µν to the perfect fluid EMT, which arise in out-of-equilibrium relativistic fluid dynamics theory [51].
In this paper we explore the viability of a simple solution in such a framework to describe at the background level the late-time accelerated expansion of the Universe, relieving at the same time the Hubble tension. The structure of the work is as follows. In Section II, we shortly review the relativistic fluid dynamics theory from which Ricci cosmology emerges. In Section III, we describe the features of our solution in the Ricci cosmology framework. In section IV, we present statistical analysis and the data used to put constraints on the parameters of our model. In Section V, we show and discuss our results. In Section VI, we draw conclusions about the viability of the model in light of the cosmological data taken into account.

II. NEAR-EQUILIBRIUM FLUID DYNAMICS
To understand the features of Ricci cosmology, it is useful to discuss some basics of fluid dynamics.
The first attempts to construct a consistent theory of dissipative fluids on a general curved spacetime are due to Eckart [52] and Landau and Lifshitz [53] in the first half of the twentieth century. After these first efforts which have been proven to be plagued with stability and causality problems, in the 1960s Müller [54] and in the 1970s Israel and Stewart [55,56] succeeded to solve these issues, as shown in [57].
Since then, further developments occurred which culminated in the theory of relativistic fluid dynamics without conserved charges reviewed in [51], with an equivalent formulation for a fluid with a U (1) symmetry based on an action principle, proposed by Kovtun et al. [58]. In the rest of the section, we will follow [51].

A. General construction
A perfect fluid at the equilibrium in a curved spacetime is described by the Energy-Momentum Tensor (EMT) where ρc 2 is the energy density of the fluid as seen by a comoving observer with the fluid 4-velocity u µ and h µν = g µν + u µ u ν /c 2 is the 3-spatial metric of the hypersurface orthogonal to the fluid 4-velocity. The pressure P (ρ) is the equilibrium pressure of the relativistic fluid, whose form represents the equation of state (EoS) of the fluid. In general, a perfect fluid at the equilibrium has a barotropic EoS given by where w is the constant EoS parameter depending on the nature of the fluid. The energy density ρc 2 , the fluid 4-velocity u µ and the metric g µν which fully characterise the description of the perfect fluid go under the name of hydrodynamic fields.
If the fluid is slightly out-of-equilibrium for the presence of dissipative effects or anisotropic non-stationary expansion, the perfect fluid description above is not enough to accurately characterize its dynamics, and gradients of such hydrodynamic fields are needed.
In general, the EMT can be written as where the subscripts (0), (1), (2), . . . indicate the number of gradients in each term of T µν . In order to obtain such higher order corrections to the perfect fluid EMT (1) scalars, vectors and tensors of such gradients are constructed and combined into tensors in such a way to split the contributions into a traceless part π µν , referred to as the shear stress tensor where < · · · > indicates symmetrization over the indices and subtraction of the trace, and a trace part Πh µν , with which is called the bulk stress. The form of these corrections is further constrained by the conservation equation of motion for the full EMT ∇ µ T µν = 0 (6) and by the fact that the fluid energy flux is not modified by the corrections at any order, i.e.
meaning that the energy density ρc 2 seen by an observer comoving with the fluid is not modified by the non-equilibrium terms. This choice corresponds to the so-called Landau frame. Instead, from the corrections above, it is evident that the local pressure P (ρ) departs from its equilibrium expression where in the second term there is no summation over i.
When we neglect anisotropies, that is possible on cosmological scales which we are interested in here, the last equation simplifies to which is an isotropic non-equilibrium pressure, due to internal friction in a fluid and can be responsible, in the cosmological context for the late-time accelerated expansion of the Universe. At the first-order, the construction described above gives rise to the well-known expression The term σ µν = 2∇ <µ u ν> is the well-known shear viscosity which becomes relevant in presence of anisotropies, and the second term is the bulk viscosity, which has been used in alternative cosmological models to Dark Energy and Inflation. The shear viscosity coefficient η and the bulk viscosity coefficient ζ, are collectively called first order transport coefficients.
By applying the same reasoning used for the first order correction T µν (1) , the second order correction to the perfect fluid EMT T µν (2) can be obtained, by considering all the linearly independent scalars, vectors and rank-two tensors, containing exactly two gradients of the fundamental hydrodynamic fields and of the metric g µν . The full details of the second order shear tensor π µν (2) and of the second order bulk stress Π (2) can be found in [51].
Here, we are only interested in those terms in the second order bulk pressure Π (2) appearing in the Ricci Cosmology framework proposed in [50], namely where R λρ is the Ricci tensor and R is the Ricci scalar. The coefficients ξ 5 and ξ 6 are two of the second order transport coefficients, which parametrize the response of the non-equilibrium fluid to the change of the background. These terms are the only terms at the second order which are vanishing in Minkowski spacetime and non-zero in a curved spacetime.
In the following, we will assume the second order transport coefficients and we will study their impact on the usual scaling of pressureless matter, radiation and cosmological constant with the redshift.

III. RICCI COSMOLOGY
A. Ricci Cosmology model Let us assume that GR describes gravity in our Universe and it is well described at large scales by the FLRW metric where a(t) is the scale factor. The Universe is assumed to undergo the usual sequence of radiation, matter and Dark Energy dominated epochs in which the contributions of the other matter components to the energy budget of the Universe is small and can be safely ignored.
The effective pressure of such matter components, due to the background expansion, gets modified by the terms in Eq. (11) and it is given by where the subscript c ranges on m, r, Λ which stands for for matter, radiation and cosmological constant, respectively, and ξ 5c and ξ 6c are the second order transport coefficients assumed constant and different for every matter component. The Ricci scalar R and the projection of the Ricci tensor R αβ along the fluid 4-velocity u α for the FLRW background are given by where H ≡ȧ a is the Hubble function. The effective pressure in Eq. (13) can be seen to be analogous to the model of nonlinear bulk viscosity for the Dark Energy proposed in [34]. The two models coincide when the phenomenological parameter introduced in [34], are ζ 0 = 0, ζ 1 = −4ξ 5Λ /c 2 and ζ 2 = −4ξ 5Λ /c 2 +ξ 6Λ , while the other coefficients are zero.
The conservation equation of the energy density ρ c c 2 for such out-of-equilibrium cosmic fluid is given bẏ with the energy density related to the Hubble function via the first Friedmann equation in the epoch dominated by the matter component c. By inserting Eq. (13) in Eq. (15), using the expressions in Eq. (14), the equation for the energy density becomeṡ By replacing cosmic time derivatives with redshift derivatives and using Eq. (16) to express the Hubble function in terms of the energy density, after rearranging the terms, the last equation becomes where the reduced second order transport coefficientŝ have been defined. From the differential equation (18), we find for the energy density the following expression Hence, the squared Hubble function for a flat Universe filled with matter, radiation and cosmological constant, all having a modified redshift scaling, reads where we have defined the dimensionless energy density parameter for the generic matter component c as with the relation Ω Λ + Ω m + Ω r = 1 valid in a flat Universe. Furthermore, the deviation parameters from usual scaling for matter, radiation and Cosmological Constant, i.e. δ m , δ r and δ Λ , which appear in Eq. (21), can be easily derived from Eq. (20) and are given by and respectively.

B. Thermodynamical priors
Before entering the observational bounds on the parameters of our models, we impose some physical bounds on them by applying the Second Law of Thermodynamics which says that the entropy of a system S never decreases [59]. These bounds can then be treated as priors for further statistical considerations.
For a fluid with energy E = ρc 2 V in a volume V , with temperature T and pressure p, from the First Law of Thermodynamics the conservation equation for the energy density of such a fluid can be deriveḋ By comparing the last equation with Eqs. (13) and (15), we find a the following differential equation for the entropy S From the thermodynamic relation for enthalpy H ≡ E + pV = T S, the entropy can be expressed as By dividing the last two equations, we arrive aṫ This must be valid for each matter component filling the Universe. For a matter component with energy density ρ c c 2 , barotropic EoS parameter w c and the constant reduced second order transport coefficientsξ 5c andξ 6c , the entropy S in terms of the scale factor, after a simple calculation, reads As we have already mentioned, from the Second Law of Thermodynamics for an isolated system, for each epoch of the Universe, there must be an increase of entropy Thus, by specializing Eq.(31) to each matter components of the Universe and using the expressions for δ m , δ r and δ Λ found above, we have direct or indirect constraints on the deviation parameters.
For the cosmological constant, the entropy is given by which increases for 2ξ 5Λ −ξ 6Λ ≥ 0. For pressureless matter, the entropy reads which increases for δ m ≤ 0. Finally, for radiation, it holds which increases for δ r ≤ 0.
In the next section, we describe the statistical analysis and the data we use to put bounds on the parameters of our model.

IV. STATISTICAL ANALYSIS & DATA
In order to assess the viability of the model described in the previous section and determine the relative importance of the two contributions to the out-of-equilibrium, we test against data four special cases of the model described above, with four different ansätze on the transport coefficients.
In the Table I, we report the assumptions made on the constant second order transport coefficients with physical priors on the deviation parameters derived from the physical requirement of the increase of entropy, we use in our statistical analysis, together with the usual priors on the cosmological parameters in common with ΛCDM (0 < Ω b < Ω m < 1, 0 < h < 1).
In the first two ansätze, we have that the Ricci scalar and the time-time component of the Ricci tensor in the effective pressure have the same effect in the change of the scaling of the matter components filling the Universe. The difference between the two cases is in the way radiation is treated: in ansatz 1, radiation is assumed to be unaffected by the modifying pressure terms, while in ansatz 2, radiation fluid deviates from conformality due to the pressure terms.
In the last two ansätze, instead, each fluid differs from the other for the transport coefficients which characterise its response to both Ricci scalar and time-time component of the Ricci tensor. Analogously to the previous ansätze, we consider two different behaviours for radiation.
In the fits of the special cases of our model, we combine the following data sets we have on the background evolution at large scales of the Universe: Type Ia Supernovae (SNeIa) from the Pantheon sample [60], the Mayflower sample of Gamma Ray Bursts (GRBs) [61], the Early-Type Galaxies (ETG) used as Cosmic Chronometers (CC) [62][63][64], the data on the Hubble parameter H 0 from the H0LiCOW collaboration [65,66], Baryon Acoustic Oscillations (BAO) from different surveys [67][68][69][70][71], and the last Planck release for the Cosmic Microwave Background (CMB) [72].
In order to evaluate the performance of our model in relieving the Hubble tension, we consider two different combinations of our data sets: the late time data set including only late time measurements of cosmological observables (SNeIa, CC, H0LiCOW, GRBs and BAO from WiggleZ), and the full data set which combines these observations with early time observations (CMB and BAO data from SDSS).
To find the best fit parameters of our model for each ansatz, we minimise the total χ 2 given by the sum of all the χ 2 associated with the data from the probes listed above by using our own implementation of a Monte Carlo Markov Chain (MCMC). Then, to see under which assumptions on the transport coefficients our model can fit the data better than ΛCDM, we have fitted ΛCDM to the same data sets and we have compared each of the ansätze on our model to it by means of the Bayes factor. Given a generic model M i , with π(θ i |M i ), the prior probability of its set of parameters θ i , and its likelihood function L i (D|θ i , M i ) ∝ e −χ 2 /2 , the Bayesian evidence E i is defined as the probability of the data D given the model M i with a set of parameters θ i Then, to compare the fit of this model to that of another model M j , tested against the same set of data, in general depending on a different set of parameters θ j , we compute the Bayes factor defined as where E j is the Bayesian evidence of the model M j . The last model M j , in our case, is ΛCDM while the model M i is one of the four special cases of our model. Then, we can evaluate the performance of our model with respect to ΛCDM in fitting the data for each ansatz considered, by using the Jeffreys' Scale [73]: if ln B i j < 1, the evidence in favour of model M i is not significant; if 1 < ln B i j < 2.5, it is substantial; if 2.5 < ln B i j < 5, it is strong; if ln B i j > 5, it is decisive, while negative values of ln B i j can be instead interpreted as evidence against model M i and thus, in favour of model M j .

A. Type Ia Supernovae
The Pantheon catalogue contains 1048 Type Ia Supernovae (SnIa), used as standard candles, in the redshift interval 0.01 < z < 2.26 [60].

B. Gamma Ray Bursts
We consider 79 Gamma Ray Bursts (GRBs) from the Mayflower sample in the redshift range 1.44 < z < 8.1 [61]. Analogously to the case of the Type Ia Supernovae, the observable for the GRBs is the distance modulus µ GRB , defined as in Eq. (40), with the corresponding offset parameter µ 0 marginalized over as for the Type Ia Supernovae.
The total χ 2 GRB for this sample is given by where ∆µ GRB = µ GRB − µ obs GRB .

C. Cosmic Chronometers
Passively-evolving Early-Type Galaxies (ETGs), i.e. galaxies with low star formation rate and old stellar populations, can be used as standard clocks or Cosmic Chronometers (CC). Differently from other probes, they give us direct information about the Hubble function over redshift ranges, from the formula where dz dt is inferred from observations. The sample used in our data analysis covers the range 0 < z < 1.97 [62][63][64]. The χ 2 H reads as where σ H (z i ) are the statistical errors on the measured values of the Hubble function H obs (z i ).

D. H0LiCOW
The H0LiCOW collaboration [65] uses 6 strong gravitationally lensed quasars with multiple images to put constraints on the value of the Hubble parameter H 0 [66]. In order to constrain the parameters of our model, we use the so-called time-delay distance, which is defined by where z L is the lens redshift, and D S , D L and D LS represent the angular diameter distances from the source to the observer, from the lens to the observer, and between source and lens, respectively, with the angular diameter distance given by .
The χ 2 for the H0LiCOW data D∆t,i (48) where σ D∆t,i are the statistical errors on the measured time-delay distances D obs ∆t,i .

E. Baryon Acoustic Oscillations
In the following, we report the five BAO data sets used in our analysis.
From the WiggleZ Dark Energy Survey [75] the physical observables taken into account are the acoustic parameter with ω m = Ω m h 2 , and the Alcock-Paczynski distortion parameter where D A is the angular diameter distance defined in Eq.(47) and D V is the volume distance given by From the SDSS-III BOSS DR12 [69], the quantities are considered, where D M denotes the comoving distance and r s is the sound horizon which is evaluated at the dragging redshift z d in Eq.(52).
In the same expression, r f id s (z d ) is the sound horizon evaluated at the dragging redshift z d , considering a fiducial cosmological model.
The dragging redshift z d can be estimated numerically as [76] z d = 1291 respectively. The sound speed for coupled photons and baryons, appearing in Eq. (54), in ΛCDM is given by whereR b is baryon-to-photon density ratio parameter, given byR with the CMB temperature T CMB = 2.726 K, while for our model, it must be modified as due to the modified scaling of matter and radiation parametrized by δ m and δ r , respectively. The other measurements for BAO used in the data analysis are the following: • from the combination of void-galaxy crosscorrelation with BAO and galaxy RSD in the CMASS galaxy catalog of the BOSS DR12 [71] D A (z = 0.57) r s (z d ) = 9.383 ± 0.077, H(z = 0.57)r s (z d ) = (14.05 ± 0.14)10 3 km/s (63) • from eBOSS DR14, a spherically-averaged BAO distance [68] D V (z = 1.52) = 3843 ± 147 r s (z d ) • from eBOSS DR14, by combining the quasar Lyman-α autocorrelation function with the quasar Lyman-α cross-correlation measurement [70,77] c H(z = 2.34)r s (z d ) = 9.00 ± 0.22 (66) For each BAO probe, the χ 2 is given by where ∆X BAO = X BAO − X BAO obs is the difference between the predicted and observed values for the observables of each probe.

F. Cosmic Microwave Background
The last data set we use is given by the Cosmic Microwave Background (CMB).
The CMB data we include in our data set are the shift parameters [78] from the last Planck data release [72]: the physical baryon density parameter ω b , the angular scale of the sound horizon at recombination and the scaled distance to recombination with the sound horizon r s and D M evaluated at the recombination redshift z * , which is given by the fitting formula [79] with the factors g 1 and g 2 given by respectively. Therefore, the χ 2 for the CMB data is where ∆X CMB = X CMB − X CMB obs is the difference between the predicted and observed values for the quantities in the vector X = {ω b , l a , R}.

V. RESULTS & DISCUSSION
The results of the fits of our model for the four ansätze are reported in the Tables II, III, IV and V. At 1σ level, we note that the values of the cosmological parameters are indistinguishable from those corresponding to the standard ΛCDM and so we have not found a relief of the Hubble tension with any of our ansätze.
From the fits with late time data sets, we have found milder upper or lower bounds on the deviation parameters of our model, with respect to the bounds obtained by using the full data set comprising late and early-time cosmological data which show more constraining power. The only exception is represented by δ Λ for the ansätze 1 and 2 which is compatible with zero in the fits with only late time data sets.
Regarding the second order transport coefficients, the fits of the model with the first two ansätze imply that ξ 50 is compatible with zero, while ξ 60 < 0 from the physical requirement of the increase of entropy for both models. The bound is tighter when the full data set is considered.
For the last two ansätze, only upper or lower bounds can be put on the transport coefficients with significantly milder bounds coming from the fit with only late time data set taken into account. Furthermore, from the negative values of the Bayesian factors reported in the tables, it can be concluded that none of the considered cases of our model has a better fit than ΛCDM to our set of data.
An analogous Hubble function to that considered in this paper has been previously studied in Refs. [80,81]. In those papers, the Hubble function arised in the framework of Quantum Field Cosmology proposed by Weinberg [82]. In his framework, the gravitational constant varies with redshift and consequently, the cosmological constant acquires a dependence on the redshift.
The differences between our results and the precedent findings stem from different aspects. In the proposed "varying Λ"CDM model (ΛCDM), the scaling of matter and radiation is modified by the same deviation parameter δ G that comes from the redshift dependence of the Newton constant G, and the deviation parameter for the cosmological constant δ Λ is related to δ G by for δ Λ , δ G ≪ 1, with the consistency relation δ G δ Λ > 0. In our model, instead, we have δ m = δ r , and the Eq.(74) does not hold so that the consistency relation is not compatible with our physical requirements for the deviation parameters, which are not constrained to be much smaller than unity. Our results are also not fully comparable when δ Λ and δ G are taken as two independent parameters in the so-called extended "varying Λ"CDM model (eΛCDM), with the last parameter still describing the deviation from usual scaling for both matter and radiation. Because of these differences, we have obtained different constraints on our model parameters with respect to those obtained in [80,81].
Nevertheless, analogously to what we have found for our model, the two modelsΛCDM and eΛCDM, tested against the combined set of CMB distance prior data from Planck 2018, BAO and SNIa Pantheon compilation, partially overlapping our more extended full data set, result to be compatible with ΛCDM. Similarly as in our model, these models can not fit the data better than ΛCDM and cannot relieve the Hubble tension, unless the local measurement by the SH0ES team [12] is added to the data set.

VI. CONCLUSION
In this paper we have explored for the first time the physical consequences of the recently proposed framework of Ricci Cosmology on the late time Universe. We have derived a simple solution of Ricci Cosmology under   the preliminary assumption that the second order transport coefficients are constant. Such an assumption of the non-zero transport coefficients induced some small deviations from perfect fluid scaling of the Universe matter components (dust, radiation, cosmological term). Further, the basic constraints (priors) for the parameters made of transport coefficients and describing deviations from the standard cosmology have been obtained from the physical requirement of the increase of entropy according to the second law of thermodynamics.
Then, we have fitted the full and the late time cosmological data sets reported in section IV (supernovea, GRBs, cosmic chronometers, H0LiCOW, BAO, CMB), we have found the bounds on the parameters of Ricci cosmology realising that for our simple model it is compatible with standard ΛCDM cosmology which statistically still fits better to the data. In order to answer the question if Ricci cosmology may give some stronger effects on the evolution of the universe it is perhaps advisable to release the assumption of the constancy of the transport coefficients. This investigation is left for future work.    The assumptions on the reduced second order transport coefficients arê ξ 5Λ =ξ 5m =ξ 5r ≡ξ 50 ,ξ 6Λ =ξ 5m =ξ 6r ≡ξ 60 , (A16) with the deviation parameters given by (A21) The inequality (A20) can thus be rewritten in terms of δ m and δ r as which for δ m ≤ 0, gives us a lower bound for δ r δ r > −1 + δ m .
Furthermore, from the first inequality in Eq.(A19), and the lower bound in Eq.(A23), we obtain Thus, the physical bounds on the deviation parameter δ r is given by For the reduced transport coefficients, from the first and third inequalities in Eq.(A19), we arrive at which combined with the first inequality in Eq.(A19) and Eq.(A20), finally gives us − |ξ 60 | 2 ≤ξ 50 ≤ |ξ 60 | 2 . (A27)
(A33) The priors from the increase of entropy are given by δ r ≤ 0, δ m ≤ 0, andξ 5Λ ≥ 0, which imply the following bounds for the deviation parameterδ Λ and reduced transport coefficientsξ 5r and ξ 5m