Impact of Collisional Matter on the Late-Time Dynamics of f ( R , T ) Gravity

We study the cosmic evolution of non-minimally coupled f (R, T) gravity in the presence of matter fluids consisting of collisional self-interacting dark matter and radiation. We study the cosmic evolution in the presence of collisional matter, and we compare the results with those corresponding to non-collisional matter and the Λ-cold-dark-matter (ΛCDM) model. Particularly, for a flat Friedmann–Lemaître–Robertson–Walker Universe, we study two non-minimally coupled f (R, T) gravity models and we focus our study on the late-time dynamical evolution of the model. Our study is focused on the late-time behavior of the effective equation of the state parameter ωe f f and of the deceleration parameter q as functions of the redshift for a Universe containing collisional and non-collisional dark matter fluids, and we compare both models with the ΛCDM model. As we demonstrate, the resulting picture is well accommodated to the latest observational data on the basis of physical parameters.


Introduction
Recent observations of the accelerated cosmic expansion [1,2] coming from Type Ia supernovae, and the confirmation of the B-mode power spectrum during the early-time inflationary era [3], have utterly changed the perspective of theoretical cosmologists with regard to the evolution of the Universe.The main aim and challenge is to find a unified description of the early and late-time eras, within the same theoretical framework.Modified gravity in general has provided several models that can consistently harbor the early-and late-time evolution of our Universe [4] (see also the reviews [5,6]).With regard to the late-time era, observations support a flat universe formed from pressureless matter with a non-zero cosmological constant Lambda Λ with equation of state (EoS) parameter ω = −1 [7].Additionally, the latest Planck data [8] favor the Λ-cold-dark-matter (ΛCDM) model for our present epoch, and it seems that the so-called dark energy dominates the evolution of the Universe at the present epoch at a percentage of 69.8% and the luminous matter contributes to the total energy density of our Universe only 4.9% while the remaining 26.8% is controlled by dark matter.The confirmation of the late-time acceleration is further supported by the temperature anisotropies in baryon acoustic oscillations (BAO) [9] and cosmic microwave background (CMB) [10].As we already mentioned, a consistent way to model the dark energy and the early-time acceleration era is to use the modified gravity theoretical framework [5,6,[11][12][13][14][15][16].However, alternative approaches can also describe these acceleration epochs [17][18][19].Modified gravity models are generally modifications of the general relativity (GR) picture.In the context of GR, one needs a negative pressure fluid to describe the dark energy era.However, without the use of the modified gravity theoretical framework, one frequently needs phantom fluids to describe an accelerated expansion.There are various modified gravity models in the literature, such as f (R) gravity [20], f (R, T) gravity [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37], and f (R, T, R µν T µν ) gravity [38][39][40][41], where R, T, R µν , and T µν stand for the Ricci scalar, the trace of the "energy momentum tensor"(EMT), and the Ricci tensor, respectively.
In this paper we use the f (R, T) gravity framework to study the late-time evolution of several cosmological models, in the presence of a collisional matter fluid.We focus our study on the behavior of the EoS parameter ω e f f and the deceleration parameter q, as functions of the redshift z.From the latest observational data, we know that the Universe made the deceleration to acceleration transition at a low redshift value z t = 0.46 ± 0.13 [42].This transition from a decelerated era to an accelerated one is usually achieved in modified gravity models (e.g., [4,[43][44][45][46][47]), and the aim in this paper is to investigate if a viable late-time phenomenology can be achieved for f (R, T) models of gravity in the presence of collisional matter.Similar studies in the context of f (R) gravity have been performed in [48,49], so in this paper we extend these in the context of f (R, T) gravity.Also, in Ref. [50][51][52], the dark energy dominance and the deceleration to acceleration transition was studied in the presence of cold dark matter for f (T) models.In addition, in Ref. [53] the late-time behavior of f (R, T) gravity models was compared to the ΛCDM model, with good agreement (see also Refs.[28, [54][55][56][57][58][59][60][61] for additional studies in the field).In this work, the collisional matter fluid mainly refers to the presence of a logarithmic term in the EoS parameter of the dark matter fluid.The f (R, T) gravity part will drive the late-time acceleration, so we mainly aim to see whether the presence of collisional matter alters the late-time behavior of the f (R, T) model, and to which extent the model can be comparable to the ΛCDM.
This paper is organized as follows: In Section 2, we briefly review the f (R, T) theoretical framework by presenting the general formalism of the model.In Section 3, we incorporate a collisional matter component in the EoS of dark matter and we express the field equations as functions of the redshift, and we study the behavior of the EoS parameter and of the deceleration parameter as functions of the redshift parameter z, for two characteristic models of f (R, T) gravity.Finally, the conclusions follow at the end of the paper.

f (R, T) Gravity: Basic Formalism
In this section, we discuss the basic formulation of f (R, T) gravity and its corresponding field equations for a Friedmann-Lema ître-Robertson-Walker (FLRW) spacetime in the presence of ordinary matter.The line element of the FLRW metric is where a(t) is the scale factor and the corresponding Ricci scalar is where H is the Hubble parameter and "dot" denotes differentiation with respect to cosmic time t.
The action of f (R, T) gravity is [21] S = 1 where f (R, T) is an arbitrary function of Ricci scalar R and of the trace of the EMT T = g µν T µν .In addition, κ 2 = 1/M 2 p , where M p is the Planck mass scale.Here L m represents the Lagrangian density of ordinary matter.By varying the above action with respect to the metric tensor, we have the following set of equations: (3) By contracting the above equation, we acquire a relation between Ricci scalar R and the trace T of the EMT: These two equations contain the covariant derivative and the d'Alembert operator denoted by ∇ and 2, respectively.Furthermore, f R (R, T) and f T (R, T) correspond to the functional derivatives of f (R, T) with respect to R and T, respectively.Additionally, the term Θ µν is defined as follows: The EMT for an isotropic perfect fluid is given by, where u µ is the four-velocity of the fluid.Here, we choose L m = −P m , which leads to following expression for Θ µν : With the help of Equation ( 5), the field Equation ( 3) can be written in the following form: where is the effective gravitational constant, and is the effective EMT.
Applying the covariant divergence to the field Equation (3), one can find [62]: It is important to see that any modified theory which involves non-minimal coupling between geometry and matter does not obey the ideal continuity equation.f (R, T) also involves this type of non-minimal coupling, so it also deviates from standard behavior of continuity equations.Here, non-minimal coupling between matter and geometry induces extra force acting on massive particles, whose equation of motion is given by [62]: where The matter energy density ρ m satisfies the following continuity equation: To obtain the standard continuity equation, we need to have an additional constraint by setting the right side of the above equation equal to zero.For the EoS P m = ωρ m , Equation ( 9) becomes: with constraint In [61], authors reconstructed some well-known models of f (R, T) theory of gravity corresponding to cosmological evolution including ΛCDM, phantom or non-phantom eras and possible phase transition from accelerating to decelerating by keeping the energy momentum tensor covariantly conserved.They utilized the additional constraint (11) in the comprehensive reconstruction scheme.We opted to follow two generic models of the form .
The significant aspects of this choice: these models represent the non-minimal coupling of curvature R and matter components T (as compared to previous studies [53]), and are reconstructed for the conserved energy momentum tensor.In the next section we will consider effective fluid to develop a generic equation involving the role of collisional fluid.

Collisional Matter Model within f (R, T) Theory and Late-Time Dynamics
The analytical results of high-energy particle detectors such as ATIC, PAMELA, and WMAP indicate that the production of the electron-positron in the Universe is larger compared to that observed by cosmic ray collisions and by supernova SNIa explosions [63].This issue helps us to find the total destruction rate of weakly interacting matter particles, and this process is collisional [64].The effects of collisional matter in the Universe's evolution were considered in Ref. [65] in the context of Einstein-Hilbert gravity, and the same problem was considered in Refs.[48,49] in the context of f (R) gravity.As was shown in Refs.[48,49], a transition from deceleration to acceleration can be achieved, and in this paper we shall investigate whether this transition can occur in the context of f (R, T) gravity, and whether the resulting picture can be comparable with the ΛCDM model.Mostly, self-interacting collisional matter models are quantified in terms of a perfect fluid in which total mass-energy density denoted by m depends on two terms, as follows: where Π stands for where ρ m 0 and Π 0 are present day values.In this case, the pressure of the collisional matter fluid is where the EoS parameter of the collisional dark matter fluid is represented by ω and 0 < ω < 1.Using Equations ( 12) and ( 13), we obtain the following relation of the total energy density of the Universe: The energy momentum conservation is quantified by the following equation: and the EMT becomes as follows: where ds is the four-velocity and satisfies the equation u µ u µ = 1.Using the FLRW line element given in (1), the energy-momentum conservation equation has the following form: Combining Equations ( 14) and ( 15), we obtain the following result: where a 0 is the current value of the scale factor.Here we can see that collisional matter can be described by Equations ( 15) and (19).In addition, Π 0 has the value and its numerical value is Π 0 = 2.14169.There is an important case that needs to be discussed, related to the fact that our Universe is filled with both matter and radiation.To this end, we assume that apart from collisional matter, the Universe is also filled with radiation, so the total energy density of the matter fluids is: where ρ r0 represents the present energy density for the radiation fluid.The total pressure in this case is given by where p m is the pressure for collisional matter and p r is the pressure from radiation.By using Equations ( 15) and ( 19), we can rewrite Equation ( 21) in the following form: We can also rewrite Equation ( 23) as follows: where χ is defined as χ = ρ r0 ρ m0 , its numerical value is χ = 3.1 × 10 −4 .In addition, g(a) is defined as follows: In the above cases, if we set ω = 0, Π 0 = 0, we recover the standard cold dark matter case, and in addition if we take χ = 0, the radiation fluid makes no contribution to the energy density of our Universe.Now, let us proceed in the calculation of the total EoS parameter of the f (R, T) model and of the deceleration parameter q(z) as functions of the redshift parameter defined in terms of the scale factor as follows 1 + z = 1/a.The field equations are rewritten as follows: where H = ȧ a is the Hubble parameter.We rewrite the right hand side of Equations ( 26) and ( 27) in terms of effective energy density ρ e f f and of the effective pressure P e f f as follows: Here and ρ e f f , P e f f have the following form: The energy density of dark energy is equal to ρ DE = ρ e f f − m , and the corresponding pressure is P DE = P e f f − P m .The energy momentum conservation is Using Equations ( 28) and (29), the above equation takes the form Using the redshift parameter z = 1 a − 1, the above equation can be written as follows: In the rest of this section, we will study the late-time evolution of the Universe filled with collisional matter and radiation, for two characteristic models of f (R, T) gravity.To this end we will numerically solve the differential Equation (34) for H(z), and we shall calculate the corresponding total equation of state parameter and deceleration parameter as functions of the redshift z.With regards to the deceleration parameter as a function of the redshift, we shall use the following formula: Accordingly, we shall compare the collisional matter f (R, T) gravity with the non-collisional matter late-time evolution.We shall consider two different values for the EoS parameter of the collisional matter, namely ω = 0.5 and ω = 0.8.Moreover, the observational value for present Hubble parameter is taken H 0 = 68.3, and finally the fractional energy density of dark matter at present day is Ω m0 = 0.24.

f (R, T) Model I
Let us first consider a pure non-minimally coupled f (R, T) model gravity of the form f (R, T) = α 1 R γ 1 T γ 2 + α 2 T. In Ref. [61], the authors reconstructed this model and investigated the instability of the model against density matter perturbations and the Dolgov-Kawasaki instability criterion.
In order for this model to be viable, the following condition must hold true: so we need to have γ 1 > 1 and α 1 > 0. For the numerical analysis, we choose the free parameters values as follows: α 1 = 20, α 2 = 15, H 0 = 68.3, and Ω m = 0.3183.We employ the numerical approach to find H(z) by solving Equation (34).We investigate the behavior of the total EoS parameter, ω e f f , and of ω DE for non-collisional matter and collisional matter.We present two different cases depending on the exponents of R and T in the above model.For the first model, we set γ 1 = 10 and γ 2 = −0.7,so that f (R, T) takes the form In the left plot of Figure 1 we present the evolution of the deceleration parameter for four different cases in terms of the redshift z.The blue curve represents the standard ΛCDM model, the green curve is for non-collisional matter, the black curve represents collisional matter, and the red curve represents collisional matter in the presence of radiation.We can clearly see from Figure 1 that the behavior of the green curve and red curve correspond to the ΛCDM, but the black curve is lower than the blue curve and changes its behavior.We can also observe the deceleration to acceleration transition point in the left plot of Figure 1 for the aforementioned cases.In the case of non-collisional matter, collisional matter, the transition from deceleration to acceleration is almost equal to the one corresponding to the ΛCDM model, and it is equal to z t = 1.9.For the collisional matter case, the transition point is at z t = 7.6 and this is larger than observations.In the right plot of Figure 1, we present the evolution of the total EoS parameter for non-collisional matter, collisional matter, and for collisional matter plus radiation.In all cases, the total EoS parameter ω e f f approaches the phantom divide line value ω e f f = −1 but it does not cross it, for all the models, at least up to the present value of the redshift (which is (z = 0).Also in Figure 2, we present the behavior of the Hubble rate as a function of the redshift (left plot) and of the dark energy density as a function of the redshift (right plot).As can be seen in both figures, the collisional matter phenomenology is qualitatively viable.

f (R, T) Model II
Let us now consider the following f (R, T) model: This model is of the form f (R, T) = f (R) + f (T), and the second derivative of this model with respect to the Ricci scalar is The validity of this model crucially depends on second-derivative f RR which should be greater than zero, which means µ > 1, ν > 1, β 1 > 0, and β 2 > 0. We chose the free variables values as follows: (β 1 , β 2 , β 3 , β 4 , µ, ν) = (5,10,15,20,25,30), and Ω m = 0.3183.We numerically study this f (R, T) model by solving the differential Equation (34).In Figure 3 we present the results of our study, and specifically in the left plot we present the deceleration parameter as a function of the redshift and in the right plot we present the total EoS parameter as a function of the redshift.The conventions for the colored curves are the same as in the previous model-that is, the blue curve represents the results for the ΛCDM model, the green curve represents the non-collisional matter, and the black curve stands for the collisional matter.Finally, the red curve represents the collisional matter when radiation is also considered.As can be seen, this model is in better agreement with the ΛCDM model, and also the transition point from deceleration to acceleration for all cases containing collisional matter are closer to the ΛCDM model.Moreover, in Figure 4 we plot the Hubble rate as a function of the redshift (left plot) and the EoS parameter for the dark energy sector.As can be seen, the results are in good qualitative agreement with the late-time phenomenological picture.In addition, the value of the Hubble rate today was calculated to be 68.4,which is a very good approximation of the actual present day value H 0 = 68.3.

Another Perspective of the Late-Time Cosmological Evolution of Dark Energy with Collisional Matter Fluid
In this section we specialize the dark energy study by using appropriately chosen variables.In this way we also discuss the issue of dark energy oscillations, which haunts several modified gravity models (see for example [5] for a review of this issue).The field Equation (26) for the flat FLRW metric takes the following form: Equation ( 39) can also be written as follows: We introduce the following variables [49] in order to make the study of the late-time cosmological evolution of f (R, T) theory more focused: Here ρ DE denotes the energy density of the dark energy, m2 is the mass scale given by [48] m2 = k 2 ρ m 0 3 .Using Equations ( 40)-( 42), we obtain the following set of differential equations: The curvature scalar can also be written in the form Now making use of Equations ( 43) and ( 44), we obtain Rearranging the terms, we finally obtain the following differential equation: By making use of redshift z = 1 a − 1 in conjunction with the relations we can rewrite the differential Equation (47) as follows: ( The focus in this section is on numerically solving the differential Equation ( 50) and thus finding the exact behavior of the function y H ≡ ρ DE ρ m 0 .This will make apparent the existence or not of dark energy oscillations in any of the models.We also compare the behavior of the collisional matter f (R, T) gravity in the presence and absence of radiation, with the non-collisional matter case.Consider first the model I which we studied previously, and after solving the differential Equation (50), we present the results of our study in Figure 5. Particularly, in the left plot of Figure 5, we plot the behavior of the function y H (z) as a function of z, and in the left plot we plot the dark energy EoS parameter ω DE (z) = P DE ρ DE , which in terms of y H (z) is equal to Also in all the plots we assume the following values for the free parameters of the models: γ 1 = 10, γ 2 = −0.7,α 1 = 20, α 2 = 15, H 0 = 68.3, and Ω m = 0.3183.As in the previous cases, the red curves in the plots of Figure 5 correspond to collisional matter with radiation, the green curve corresponds to collisional matter without radiation, and the black corresponds to non-collisional matter.As can be seen, the presence of collisional matter on the cosmological model slightly changes the quantitative picture of the model, but the qualitative behavior is not altered.Thus, by using the formalism we presented in this section, we demonstrated that the presence of collisional matter provides us with qualitatively viable phenomenology.

Figure 5 .
Figure 5.The evolution of y H (z) as a function of the redshift (left), and the evolution of the dark energy EoS parameter as a function of the redshift (right).The free parameter values were chosen to be, α 1 = 20, α 2 = 15, γ 1 = 10, γ 2 = −0.7,and Ω m = 0.3183.