Consequences of energy conservation violation: late time solutions of Λ(T)CDM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda (\mathsf{T}) \mathsf{CDM}$$\end{document} subclass of f(R,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathsf{R},\mathsf{T})$$\end{document} gravity using dynamical system approach

Very recently, Josset and Perez (Phys. Rev. Lett. 118:021102, 2017) have shown that a violation of the energy-momentum tensor (EMT) could result in an accelerated expansion state via the appearance of an effective cosmological constant, in the context of unimodular gravity. Inspired by this outcome, in this paper we investigate cosmological consequences of a violation of the EMT conservation in a particular class of f(R,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathsf{R},\mathsf{T})$$\end{document} gravity when only the pressure-less fluid is present. In this respect, we focus on the late time solutions of models of the type f(R,T)=R+βΛ(-T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathsf{R},\mathsf{T})=\mathsf{R}+\beta \Lambda (-\mathsf{T})$$\end{document}. As the first task, we study the solutions when the conservation of EMT is respected, and then we proceed with those in which violation occurs. We have found, provided that the EMT conservation is violated, that there generally exist two accelerated expansion solutions of which the stability properties depend on the underlying model. More exactly, we obtain a dark energy solution for which the effective equation of state depends on the model parameters and a de Sitter solution. We present a method to parametrize the Λ(-T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda (-\mathsf{T})$$\end{document} function, which is useful in a dynamical system approach and has been employed in the model. Also, we discuss the cosmological solutions for models with Λ(-T)=8πG(-T)α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda (-\mathsf{T})=8\pi G(-\mathsf{T})^{\alpha }$$\end{document} in the presence of ultra-relativistic matter.


Introduction
Today's astrophysical measurements reveal that the Universe is experiencing an accelerated expansion phase [1][2][3][4][5][6][7][8][9][10][11]. These sets of observational data have driven the quest for convincing theoretical explanations of such a phenomenon. Among the various proposed models, the most popular one a e-mail: h.shabani@phys.usb.ac.ir b e-mail: ah.ziaie@gmail.com is the theory of general relativity (GR) modified by a cosmological constant term , which is called "the concordance" or CDM model [12]. In this model, it is assumed that the term may take over the recent eras of the dynamical evolution of the Universe after domination of what is called "Dark Matter" (DM), of which the interactions are still somewhat obscure. Observational data have discovered that at least 70% of the total energy budget of the Universe is in the form of the so-called "Dark Energy" (DE), which is regarded as a cosmic medium with unusual properties attributed to cosmological constant effects. These data show that the CDM model is in good shape [13][14][15]. In spite of its fine agreement with the observation data, there are two major concerns in this context; the first one is referred to as "the cosmological constant problem" which asks the question of the origin and the great disagreement between theoretical and expected values of the cosmological constant [16][17][18]. The other problem deals with this puzzlement of why we happen to live in a special era of evolution of cosmos where the contribution of , DM and the baryonic matter are of the same order. This is pointed out as "cosmic coincidence problem" in the literature.
These issues have motived people to seek for some other theoretical foundations or at least apply some modifications to the assumed CDM model. In this respect, cosmological scenarios with running have been proposed. The first developments, in this context, have been made by Shapiro et al. [19][20][21][22]. They have shown that there is no sturdy evidence to indicate that the cosmological constant is running or not. This fact encourages one to investigate cosmological scenarios within different theoretical backgrounds that admit running cosmological parameters. Up to now, different running cosmological constant models have been proposed, among which we can quote a time dependent cosmological constant motivated by quantum field theory [22][23][24], a running vacuum in the context of supergravity [25], (t) cosmology induced by Elko fields [26], running cosmological constant via covariant/non-covariant parametrization [27] and some more [28][29][30].
In this paper we work on a particular subclass of f (R, T) gravity, in which R and T are the Ricci scalar and the trace of EMT, respectively. Firstly, this model was introduced by Harko et al. [31] and later has been widely investigated in [32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51]. In the background of f (R, T) gravity, we have investigated a linear combination of the Ricci scalar and an arbitrary function of the trace of EMT, i.e., f (R, T) = R + (T). In this model, the Einstein gravity has been modified by a minimally coupled "trace-dependent" cosmological constant. One may find some efforts made to elaborate on the cosmological features of (T)CDM in the literature. The idea of a running cosmological constant, (T), probably dates back to the paper of Poplawski [52]. He found that (T) gravity will reduce to Palatini f (R) gravity when the pressure of the fluid is neglected. Besides, he concluded that cosmological data are consistent with (T) gravity without any knowledge about the functionality of (T) parameter. In [53], Bianchi Type-V cosmological solutions have been derived. 1 The locally rotationally symmetric (LRS) Bianchi type-I cosmological models have been considered in [55]. In most work on (T) gravity, the EMT is forced to be conserved. With this assumption, the authors of [37,40] have shown that these types of models lead to an accelerated expansion era with an undesirable present value for the EoS parameter.
Up to now, cosmological consequences of the violation of the EMT conservation have not been studied properly. It may be a good idea to consider the minimally coupled part of the f (R, T) model as a running cosmological constant and inspect its cosmological solutions. As we shall see, when the conservation of EMT is allowed to be violated, it would result in a DE era which is accompanied by an observationally allowed present value for the EoS parameter. Specifically, such a model mimics a de Sitter solution. Interestingly, a similar DE solution has been reported in [56], which arises from the violation of EMT conservation. The authors of [56] have pointed out that violation of EMT conservation can be predicted by modified quantum mechanical models or by models that utilize the causal set approach to quantum gravity. 2 In our analysis, we have applied the dynamical system approach and employed a useful method to parametrize the (T) function. In this method we have presented a way to 1 However, it seems that the authors of [53] have made a wrong assumption as regards the conservation of EMT. For more details, one can compare the paragraph above Eq. (6) in the original paper with equation (2.1) in [54]. 2 In [56], the authors have obtained an effective cosmological constant in the context of unimodular gravity. determine whether or not a given model might lead to a stable late time solution.
The present paper has been organized as follows: In Sect. 2, the field equation of f (R, T) gravity has been reviewed. Besides, the relevant dimensionless parameters which will be used in the construction of subsequent equations are introduced. A discussion of the EMT conservation will be given as well. In Sect. 3, we discuss the late time solutions of the only conserved f (R, T) model. In this case, some issues have already been illustrated, however, we present some other features. Section 4 is devoted to a description of the mentioned method. In Sect. 5, we consider models in which the (T) function shows a power-law behavior. These models will be investigated independently. Finally, in Sect. 6 we summarize our results.

Field equations of f (R, T) gravity
In this section, we present the field equations of f (R, T) modified gravity (MG) and discuss the conservation of EMT. We assume that a pressure-less fluid (DM along with baryonic matter) and ultra-relativistic matter are present. The action of f (R, T) gravity can then be written as where we have defined the Lagrangian of the total matter as In the above definitions, we have used R and T (p, u) ≡ g μν T (p, u) μν as the Ricci curvature scalar and the trace of EMT of pressure-less and ultra-relativistic fluids (for which we get these matters as the total matter content of the Universe), respectively. The letters "p" and "u" indicate the pressureless and ultra-relativistic fluids and g is the determinant of the metric, κ 2 ≡ 8π G is the gravitational coupling constant, and we set c = 1. Since the ultra-relativistic fluid has a traceless EMT i.e., Hereupon, we will drop the letter "p" from the trace of the pressure-less matter for simplicity. 3 The EMT T (p, u) μν is defined as the Euler-Lagrange expression of the Lagrangian of the total matter, i.e., The field equations for f (R, T) gravity can be obtained [31]: where (p, u) and for the sake of convenience, we have defined the following functions as derivatives with respect to the trace T and the Ricci curvature scalar R: We consider a spatially flat, homogeneous and isotropic Universe, which is described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric: where a(t) denotes the scale factor of the Universe. Therefore, the field equation (4), by assuming the metric (7), leads to as the modified Friedmann equation, and as the modified Raychaudhuri equation. In Eqs. (8) and (9), H indicates the Hubble parameter. Hereafter, we work on the Lagrangians which include minimal coupling between the trace of EMT and the Ricci scalar, i.e., where β can control the strength of the coupling. Since, for the pressure-less matter we have T (p) = −ρ (p) , in order to avoid ambiguity due to the negative sign, we work hereafter with the (−T) function. In fact, to guarantee that (T) is always a real-valued function, we consider (−T) instead. For this class of f (R, T) models we have F(R, T) = −β (−T) 4 and F(R, T) = g (R). The field equations (8) and (9) can then be rewritten as and 2Ḣ where the arguments have been dropped for convenience. In order to construct the dynamical system for the field equations (11) and (12), it is helpful to define a few dimensionless variables and parameters by where we have used Ricci scalar, R = 6(Ḣ + 2H 2 ) for metric (7), within the definition of (15). Moreover, we use the following definitions: In general, eliminating T from (20) and (21) yields n = n(s). Describing the models with n = n(s) instead of (−T), can be suitable in dynamical system analysis. In the following subsections we discuss the consequences of conservation/non-conservation of EMT, which in turn, results in a key equation that helps us to study the dynamical evolution of the models. We classify minimally coupled models as those that respect EMT conservation and those that do not.

Models which obey the conservation of EMT
In this subsection, we present the results that come from considering the EMT conservation. If we apply the Bianchi identity to the field Eq. (4) and assume that the conservation of EMT holds for pressure-less and ultra-relativistic fluids, independently, we geṫ We then find the following constraint for the pressure-less fluid 5 : After straightforward but lengthy algebra, we arrive at a specific form for (−T), as follows: where C 1 and C 2 are constants of integration. It means that solution (25) is the only subclass of f (R, T) theories of gravity with minimal coupling that respect the conservation of EMT. For this solution we obtain x 5 = x 4 /2, which reduces the space constructed from the variables of the theory.

Models which violate the conservation of EMT
These models do not generally respect conservation laws (22) and (23). Applying the Bianchi identity to the field Eq. (4) leads to the following equation between the function F(R, T), the EMT and its trace: where the argument of F(R, T) has been dropped for abbreviation. Notice that in the above equation, we have used F (p) in the corresponding terms of pressure-less fluid, since only T (p) would appear in the argument of F(R, T); we further note that the function F and its derivative are zero for the ultra-relativistic fluids. Equation (26) has a specific form such that we can consider the evolution of two fluids separately, i.e., we can write 5 See [37,40] for more details.
where p (p) = 0 has been used. From Eq. (28) we can deduce that in the minimal form of f (R, T) gravity, the conservation of EMT for ultra-relativistic fluid, i.e., Eq. (23) can be always assumed, at least, as long as mutual interactions are not taken into account. Therefore, regarding the choice (10) for f (R, T) function, we can rewrite Eq. (27) as where we have used ρ (p) = −T. Once the function (−T) is determined, the dependency of −T and thus ρ (p) on the scale factor can be calculated. More precisely, Eq. (29) can be simplified as where T 0 and a 0 denote the present values for the scale factor and EMT trace. Note that, in general, the integral on the left hand side of Eq. (30) may not be simply solved. Moreover, after the integration process, it may not be possible to clearly write the density as a function of scale factor. Let us choose the functionality of parameter as (−T) = χ 2 (−T) α , whence we obtain where C is a constant of integration. In this case, for later application, let us rewrite Eq. (29) in terms of the pressureless fluid density in the following form: or correspondingly, eliminating time gives In Sect. 3 we present an overview of the cosmological implications of the only conserved models, i.e., the models with (−T) = κ 2 √ −T. The dynamical system representation of this case has been considered in [40]. However, in this section we review the corresponding cosmological consequences of this case to complete our discussion. Moreover, we will present new details that have not been considered before. In Sect. 4 we consider cosmological behavior of models of type f (R, T) = R + β (−T) for a general (−T) function via the dynamical system approach. We will see that relaxing the EMT conservation which gives Eq. (29), leads to some interesting features; DE solutions will be achieved. In Sect. 5 we study the models with (−T) = χ 2 (−T) α when the EMT conservation has not been considered. Among non-conserved models there are two special cases that the equivalent dynamical system cannot be constructed properly. More precisely, as we will see the process of recasting the field equations into equivalent dynamical system would breaks down for cases with α = 1 and α = −1/2. In Sects. 5.1 and 5.2 we will consider these cases algebraically.

Conserved (T)CDM model in phase space
In this section we present a brief review on the cosmological solutions of the only case that respect the EMT conservation. The EMT conservation leads to the constraint Eq. (24) which gives Eq. (25) for parameter, as the only solution. To illustrate late time effects of the extra term √ −T, we put aside the ultra-relativistic fluid. Therefore, field equations (8) and (9) for the model f (R, T) = R + βκ 2 √ −T will take the following form: The energy conservation yields solution (22) for which the solution is given by 0 . We can check that differentiating Eq. (34) with respect to time together with usingρ (p) = −3Hρ (p) gives Eq. (35). This means that solutions of equation (34) would satisfy Eq. (35). Equation (34) can then be solved to give where we have set κ 2 = 1 and a(t = 0) = 0. We can also check that solution (36) reduces to the standard matter dominated era solution for β = 0 and (p) 0 = 1. From solution (36), the age of Universe can be calculated as where t (+) 0 shows the solution with positive sign and t (−) 0 the solution with negative sign between the two terms in brackets. We suppose that β = bH 0 where b is a constant, therefore we have In the MG theories we can define an effective EoS parameter as w (eff) ≡ −1 − 2Ḣ /3H 2 . For the conserved (T)CDM model, using Eqs. (34) and (35) along with the solution (36) leads to the following solution for the effective EoS: As can be seen, this solution goes to zero for early times and to −1/2 in the late times. Besides, we can obtain the fluid density for this case using ρ (p) = ρ (p) 0 a −3 , however, in this case the scale factor follows the form given in (36). In Fig. 1, we have drawn the related plots for cosmological parameters discussed above. The upper left plot presents the age of Universe for both positive and negative solutions obtained in (38). In this plot the orange line denotes the age of Universe for t (+) 0 and the green one indicates t (−) 0 . We have plotted the rest of the diagrams for b = ±0.9. For the purpose of comparison, we have employed a black line for the scale factor of the standard pressure-less fluid dominated era, i.e. a (sp) , for which a(t) ∼ t 2/3 holds.
For reconstructing the dynamical system representation of Eqs. (34) and (35) we can use the dimensionless variables (16)- (19) (remember that in this case we have x 5 = x 4 /2). In terms of these variables we obtain where we have redefined (DE) ≡ −βx 4 . Note that since x 4 is always positive [see the definition of (16)], it restricts the allowed values of β to negative values, in order that (DE) stays positive. The fixed points of this dynamical system are presented in Table 1. Also, we have drawn in Fig. 2, the phase portrait of this model in the ( (u) , (DE) ) plane.  . This family of models can be interesting as they make a minor modification to GR. In fact we can interpret these types of models as a (T)CDM theory, which implies a matter density dependent cosmological constant. Our aim here is to consider late time solutions only, hence we do not include the ultra-relativistic fluid. To reconstruct Eqs. (11) and (12) in terms of a closed dynamical system for g(R) = R, we use the definitions (16), (17), (20) and (21). The last two ones play the role of a (kind of) parametrization 6 in determining the functionality of (−T). This parametrization can be suitable at least for some well-defined models. Generally, eliminating the trace T between the definitions (20) and (21) leaves us with a relation between the n and s parameters, which results in a function n(s). In fact, each f (R, T) = R + β (−T) model can be specified by a function n(s). Models with constant s and n shall be considered in Sect. 5. Equations (11), (12) and (29) can be rewritten in terms of the dimensionless variables as follows: 2Ḣ From the definition for effective EoS parameter and Eq. (44) we get Using Eqs. (43)-(45) the following autonomous differential equations can be obtained: The above system admits three fixed points with the properties we have listed in Table 2. Table 2 shows that its elements depend on the value of parameter n, generally. Note that the DE density parameter is defined so that the relation (p) + (DE) = 1 holds and n ≡ dn/ds at the fixed point. From for which (p) and thus (DE) also depend on this parameter. Therefore, only some specific models can give rise to an accelerated expansion solution in the late times. Especially, to be consistent with the observational measurements, the value of n parameter can be much more confined. One of the eigenvalues of P (DE) 2 is zero, which shows that it is a non-hyperbolic critical point and its stability properties cannot be determined by linear approximation techniques. Hence, we focus on the solution characterizing the fixed point P The fixed points shown in Table 2 are solutions of the system dx 4 /dN = 0, dx 5 /dN = 0. From this fact we can conclude that, for any arbitrary function, namely, f (x 4 , x 5 ) we must have d f (x 4 , x 5 )/dN = 0 at the equilibrium points. Hence, for the parameter s = s(x 4 , x 5 ) we obtain Therefore, all solutions originating from the presence of the function (−T) must satisfy the conditions s = 0 or n(s) = s − 1. As can be seen, the latter condition holds for the fixed points P , the constraint equation x 5 = sx 4 = −2(n + 1)/β must hold. The matter fixed point is related to the geometrical sector of the Lagrangian, i.e., R, and therefore it is not necessary to satisfy the conditions s = 0 or n(s) = s − 1. Specifying a function (−T) with the corresponding n parameter, the condition n(s) = s − 1 leads to an algebraic equation with possible s i roots. 7 Among them, we only pick out those which exhibit DE properties. Briefly speaking, there may be some s i solutions for which the DE fixed points (or at least one) may pass necessary conditions so that a stable late time solution for the underlying model could be achieved. We note that the functions n(s) and n (s) are calculated for these s i solutions. True cosmological solutions are those that include an unstable DM fixed point which is connected to a stable DE one. Therefore, discarding the fixed point P and P (DM) are determined as follows: P (DM) is unstable provided for n < 1 n < 1−n n , and for n > 1 n < 0.
As a result, in order to have the allowed DM and DE solutions, it suffices that conditions (50) be satisfied. To complete this section, we explore the discussed method for two specific models. In the next section, we discuss the cosmological solutions for models with a power-law (−T) function.
• Models with f (R, In this case we obtain where the equation n(s) = s − 1 for (52) gives rise to the solutions s 1 = α and s 2 = −β. Therefore, a true cosmological solution can be achieved provided that the following conditions hold: • Models with f (R, n (s) = αγ The only solution for the equation n(s) = s − 1 is s * = α which leads to n (s * ) = 1 + γ /α. Hence, the conditions for a true cosmological solution read n (s * ) = α + γ α < 1, 0 < α < 1.

Late time solutions for models with power-law (−T) function
In this section, we consider a class of models of type f (R, T) = R+βκ 2 (−T) α , which violate the EMT conservation. We will investigate these types of models as dynamical systems as before. We then study, via considering the equilibrium points, the cosmological solutions provided by these models. However, for two values of α, the dynamical system approach does not work, thereby encouraging us to study these specific cases algebraically. In Sects. 5.1 and 5.2 we will study these specific cases and in Sect. 5.3 we deal with the general form of f (R, T) = R + βκ 2 (−T) α gravity.

Models of type
For α = 1, Eq. (33) reduces to the following equation: for which the solution is given by For β = 0, the above solution leads to the standard form for DM energy density and behaves as a 2 for large values of β parameter. In case in which β = 1 we have ρ (p) (a) = ρ (p) 0 . Moreover, for α = 1, modified Friedman Eqs. (8) and (9) lead to Inserting the solution (63) into Eq. (64) and solving the resultant equation, we get the scale factor as a function of time, as This solution is valid only for β < 2/3 and it leads to the standard form for the pressure-less matter dominated era for β = 0. We can also check that the solutions (63) and (66) satisfy Eq. (65). Applying the definition given for effective EoS on Eqs. (64) and (65) we see that these types of models correspond to a constant value, w (eff) = β/(2 − 3β), which for β < 2/3 gives w (eff) > −1/3. However, there is a special case; Eqs. (63), (64) and (65) yield a de Sitter solution for β = 1. As a result, power-law models with α = 1 result in a single decelerated expanding cosmological state for the Universe (at all times) without ever passing through a pressure-less matter dominated era. However, these models predict a single de Sitter state as well.

Models of type
In this case we have − 3 2 + T = 0; thus Eq. (33) reduces to the following simple form: which has the following solution for the matter density in terms of the scale factor: where we have set ρ (p) (a = 1) ≡ ρ 0 . This solution in the early times where a → 0 behaves as a −3 and in the late times where a → ∞, as (−β/2) 2/3 . The Friedman equations can then be obtained: As can be seen, the first modified Friedman equation reduces to its standard form (i.e., its form in GR), and for this reason writing down the equations as a physically consistent dynamical system breaks down. Exact solutions of Eq. (69) with energy density of DM given by (68) From the solution (71) we see that the term in square brackets reduces to the standard form for pressure-less matter, i.e., (3H 2 0 /2) 2/3 , once we set β = 0. The above solutions have been obtained so that in the present time they would be equal to unity. The solution (72) shows that in the late times, we have a de Sitter solution for −2 < β < 0 (notice that for the late time solution, i.e., (−β/2) 2/3 , Eqs. (69) and (70) give H = constant andḢ = 0, respectively). In the left panel of Fig. 3 we have plotted the scale factor for three cases; a numerical diagram has been drawn from Eqs. (69) for (68) in purple color, an asymptotic curve in the late times, i.e. the solution (72) in brown color, and the scale factor for the standard pressure-less dominated era in blue, for comparison. We used β = −1 and ρ (p) 0 = 1 for plotting purpose. In the right panel of Fig. 3, we have presented the effective EoS for the same values of parameters β and ρ (p) 0 . These diagrams show that in the model f (R, T) = R − κ 2 (−T) −1/2 , the Universe could experience an accelerated expansion in the late times, when only a pressure-less fluid is present. It is interesting to note that the unusual interaction between the only cosmological fluid and curvature of space-time, which leads to the violation of EMT conservation, has the consequence of a late de Sitter epoch in the evolution of the Universe.

Models of type
In Sects. 5.1 and 5.2 we have considered cosmological implications of two specific models that could not be analyzed by the dynamical system approach. In the present subsection we will provide a comprehensive study of these models via this approach. We will see that unlike the conserved case, there is, however, a de Sitter solution in the case of non-conserved models. For those models in which f (R, T) = R+β (−T), the field equations (8) and (9) reduce to the following equations: Substituting (−T) = κ 2 (−T) α into Eqs. (73) and (74) using (29) leads to As can be seen, Eqs. (75) and (76) reduce to Eqs. (64) and (65) for α = 1 and in the case of α = −1/2, these equations reduce to (69) and (70) in the absence of ultra-relativistic fluid. Using the definitions (16), (18) and (19), the above equations can be rewritten as Finally, the dimensionless evolutionary equations for variables (u) and x 4 are obtained: Note that in the case of the power-law dependency, we have The fixed point solutions of the system (81) and (82) are calculated in Table 3. Defining (DE) ≡ −β α + 1 2 x 4 , provided that (p) + (DE) + (u) = 1, the fixed point P (DE) can be accounted as a DE solution for which we have (DE) = 1. The condition for accelerated expansion, w (eff) < −1/3 is satisfied only for − 1 2 < α < 1. It is interesting to note that only for − 1 2 < α < 1 both eigenvalues become negative, simultaneously. Therefore, in non-conserved class of models of the type f (R, T) = R + βκ 2 (−T) α , there is a stable solution for late times with the following properties: It is noteworthy that, for α → 0 ± , we have w eff → −1 so that we can get a DE solution with observationally accepted values of the EoS parameter. Planck 2015 measurements show that the Universe is undergoing an accelerated expansion driven by DE for which the present values of the EoS parameter lie in the interval −1.051 < w DE 0 < −0.961 [58]. This fact imposes a constraint on the α parameter: −0.024 α 0.020. In addition to the existence of a DE fixed point, there are unstable solutions which indicate domination of the ultra-relativistic and pressure-less fluids, i.e., the points P (u) and P (p) , respectively. Two other solutions are Px 1 and Px 2 , which are not physically interesting, since they do not correspond to dominant cosmological epochs. 8 Nevertheless, Px 1 can be a late time solution for small values of α. Table 3 shows that the fixed point Px 1 is a stable one for small values of α. Therefore, depending on the value of α, each of these solutions may be the late time solution. To show the two possibilities, we have plotted in Fig. 4 the phase space diagrams for the values α = 0.02 and α = −0.02. The red solid circle denotes the fixed point P (DE) , the green one indicates Px 1 , the purple solid circle shows P (p) , the cyan one indicates P (u) and the orange solid circle shows Px 2 . The diagrams show that physically interesting trajectories begin from P (u) , pass along P (p) and then terminate at either P (DE) or Px 1 . Figure 4 also shows that, for α < 0, the fixed point Px 1 is for the late time solutions and for α > 0 the solution P (DE) would be chosen. These points will overlap for α = 0, that is, the GR model plus a cosmological constant, known as the celebrated CDM model. The fixed point Px 2 coincides with the fixed point P (u) for α = 4/3, otherwise it is physically meaningless. The fixed point Px 2 is always unstable, that is, the functions f (α) and g(α) never get the same negative sign or pure imaginary values. In Fig. 4 the diagrams for the evolution of the matter density parameters and the effective EoS parameter are plotted for α = 0.02, as well.

Conclusion
In this work we have investigated the cosmological consequences of violation of EMT conservation for a class of f (R, T) theories of gravity. We have considered both the ultra-relativistic fluid and DM in a spatially flat, homogeneous and isotropic background given by the FLRW metric.
We have studied models of type f (R, T) = R + β (−T) which we call minimally coupled f (R, T) models. A specific model can be considered, a (T)CDM model which allows for a density dependency for the cosmological constant.
Firstly, we have presented the field equations of f (R, T) gravity and defined some dimensionless variables. We also classified the minimal models to those that respect the conservation of EMT and those that do not. The former models have been considered elsewhere, however, to complete our study, we have briefly reviewed their cosmological solutions through the dynamical system approach. Some new results have been obtained as well. We have algebraically showed that these types of models cannot be accepted, since they have a late time solution with an undesirable EoS parameter. Their EoS parameter varies from zero to −1/2, which is not obser-vationally confirmed. Thus, considering the EMT conservation, GR theory modified by a minimal (−T) function has still the problem of incompatibility with recent observational outcomes. The latter models do not respect the conservation of EMT, for which a modified version of DM density conservation has been obtained. We have shown that, in the minimal models of f (R, T) gravity, it is always possible to consider the evolution of ultra-relativistic fluid and DM independent of each other, as long as interactions are turned off. Therefore, only a modification in the behavior of DM density can provide a different cosmological scenario, at least in the late time epochs.
To consider the cosmological consequences of the violation of EMT conservation, we presented a general method to formulate the dynamical system equations for generic minimally coupled models. We have defined two dimensionless n and s parameters constructed out of the (−T) function and its derivatives and showed that the resulted autonomous equations will depend upon these parameters. As a result, we have obtained a set of closed equations of which the solutions, and hence the stability properties of them, will be controlled by these parameters. We have discussed that, at least for well-behaving models, we can parametrize (−T) function in terms of a function n(s). We illustrated that all fixed points originating from the (−T) function must lie on the line n = s − 1. In other words, for every function n(s) all fixed point solutions must occur at the location where the n(s) curve intersects with the line n = s − 1. The fixed points, representing the ultra-relativistic and pressure-less matter domination, are not subject to the above discussion as they are solutions of GR. Briefly speaking, this method shows that there generally exist two fixed point solutions which represent accelerated expansion/DE era in the late times. We have applied the method to inspect two models specified by f (R, T) = R + a(−T) α + b(−T) −β and f (R, T) = R + a(−T) α exp b(−T) γ and discussed the validity of their solutions.
As a special case that cannot be explained by the above technique, we investigated the cosmological behavior of models with (−T) = κ 2 (−T) α . In this case, we included the ultra-relativistic fluid and discussed the late time solutions. We found that there are two different DE solutions of which the properties depend on the constant α. Nevertheless, each f (R, T) = R + β (−T) model accepts only one solution, by which, accordingly, such types of f (R, T) functions can be classified as two different models. We have argued that observationally consistent models may be constructed by small values of α. For example, Planck 2015 measurements have shown that if we believe in DE as one of the ingredients of the Universe which is presently driving the observed accelerated expansion, its EoS parameter must lie within −1.051 < w (DE) 0 < −0.961. This fact restricts us to accept −0.024 < α < 0.02. We have shown that, for the two specific models with α = 1 and α = 1/2, the dynamical system approach does not work, due to the structure of the related equations. Algebraic treatments showed that the former model, in general, indicates a single decelerating late time cosmological era. However, there is an exact single de Sitter solution, as well. For the determined values of α and β parameters, we obtained a constant value for the EoS parameter with a specific value lying within the range w α=1 > −1/3. The latter case gives a proper solution including a connected matter and DE dominated eras. This model accepts a de Sitter solution in the late times. Finally, we would like to end this article by highlighting the importance of examining the observable signals of MG theories, in order to test the physical validity of the resulted models. If experiments confirm that a modified version of GR can explain observations better than the original version, the results could shed light on some fundamental cosmological questions. Modified gravity theories have been utilized successfully to account for galaxy cluster masses [59], the velocity field of DM and galaxies [60], the cosmic shear data [61], the rotation curves of galaxies [62,63], velocity dispersions of satellite galaxies [64], and globular clusters [65]. These theories have also been used to propose an explanation for the Bullet Cluster [66] without resorting to nonbaryonic DM; see also [67] and the references therein. However, among the MG theories that have been proposed so far, Rastall's gravity touches one of the cornerstones of GR, i.e., the conservation of EMT [68] and, interestingly, this issue has been entered within the context of MG theories [69,70]. While in the present work, we have studied cosmological consequences of violation of EMT in the framework of f (R, T) gravity theory, it is of utmost importance to seek for observational evidence (such as the gold sample supernova type Ia data [71], SNLS supernova type Ia data set [72] and X-ray galaxy clusters analysis [73]) that could distinguish between the resulting models from this theory and GR. However, observational features of this theory need to be carried out with more scrutiny and dealing with this issue is beyond the scope of the present paper.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .