Bouncing cosmological solutions from \(f(\mathsf{R,T})\) gravity

The European Physical Journal C, May 2018

In this work we study classical bouncing solutions in the context of \(f(\mathsf{R},\mathsf{T})=\mathsf{R}+h(\mathsf{T})\) gravity in a flat FLRW background using a perfect fluid as the only matter content. Our investigation is based on introducing an effective fluid through defining effective energy density and pressure; we call this reformulation as the “effective picture”. These definitions have been already introduced to study the energy conditions in \(f(\mathsf{R},\mathsf{T})\) gravity. We examine various models to which different effective equations of state, corresponding to different \(h(\mathsf{T})\) functions, can be attributed. It is also discussed that one can link between an assumed \(f(\mathsf{R},\mathsf{T})\) model in the effective picture and the theories with generalized equation of state (EoS). We obtain cosmological scenarios exhibiting a nonsingular bounce before and after which the Universe lives within a de-Sitter phase. We then proceed to find general solutions for matter bounce and investigate their properties. We show that the properties of bouncing solution in the effective picture of \(f(\mathsf{R},\mathsf{T})\) gravity are as follows: for a specific form of the \(f(\mathsf{R,T})\) function, these solutions are without any future singularities. Moreover, stability analysis of the nonsingular solutions through matter density perturbations revealed that except two of the models, the parameters of scalar-type perturbations for the other ones have a slight transient fluctuation around the bounce point and damp to zero or a finite value at late times. Hence these bouncing solutions are stable against scalar-type perturbations. It is possible that all energy conditions be respected by the real perfect fluid, however, the null and the strong energy conditions can be violated by the effective fluid near the bounce event. These solutions always correspond to a maximum in the real matter energy density and a vanishing minimum in the effective density. The effective pressure varies between negative values and may show either a minimum or a maximum.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-5886-x.pdf

Bouncing cosmological solutions from \(f(\mathsf{R,T})\) gravity

Eur. Phys. J. C Bouncing cosmological solutions from f (R, T) gravity Hamid Shabani 1 Amir Hadi Ziaie 0 0 Department of Physics, Kahnooj Branch, Islamic Azad University , Kerman , Iran 1 Physics Department, Faculty of Sciences, University of Sistan and Baluchestan , Zahedan , Iran In this work we study classical bouncing solutions in the context of f (R, T) = R + h(T) gravity in a flat FLRW background using a perfect fluid as the only matter content. Our investigation is based on introducing an effective fluid through defining effective energy density and pressure; we call this reformulation as the “effective picture”. These definitions have been already introduced to study the energy conditions in f (R, T) gravity. We examine various models to which different effective equations of state, corresponding to different h(T) functions, can be attributed. It is also discussed that one can link between an assumed f (R, T) model in the effective picture and the theories with generalized equation of state (EoS). We obtain cosmological scenarios exhibiting a nonsingular bounce before and after which the Universe lives within a de-Sitter phase. We then proceed to find general solutions for matter bounce and investigate their properties. We show that the properties of bouncing solution in the effective picture of f (R, T) gravity are as follows: for a specific form of the f (R, T) function, these solutions are without any future singularities. Moreover, stability analysis of the nonsingular solutions through matter density perturbations revealed that except two of the models, the parameters of scalar-type perturbations for the other ones have a slight transient fluctuation around the bounce point and damp to zero or a finite value at late times. Hence these bouncing solutions are stable against scalar-type perturbations. It is possible that all energy conditions be respected by the real perfect fluid, however, the null and the strong energy conditions can be violated by the effective fluid near the bounce event. These solutions always correspond to a maximum in the real matter energy density and a vanishing minimum in the effective density. The effective pressure varies between negative values and may show either a minimum or a maximum. 1 Introduction Today, the standard cosmological model (SCM) or the bigbang cosmology has become the most acceptable model which encompasses our knowledge of the Universe as a whole. For this reason it is called also the “concordance model” [1]. This model which allows one to track the cosmological evolution of the Universe very well, has matured over the last century, consolidating its theoretical foundations with increasingly accurate observations. We can numerate a number of the successes of the SCM at the classic level. For example, it accounts for the expansion of the Universe (Hubble law), the black body nature of cosmic microwave background (CMB) within the framework of the SCM can be understood and the predictions of light-element abundances which were produced during the nucleosynthesis. It also provides a framework to study the cosmic structure formation [2]. However, though the SCM works very well in fitting many observations, it includes a number of deficiencies and weaknesses. For instance some problems which are rooted in cosmological relics such as magnetic monopoles [3,4], gravitons [5–9], moduli [10–13] and baryon asymmetry [14,15]. Despite the self-consistency and remarkable success of the SCM in describing the evolution of the Universe back to only one hundredth of a second, a number of unanswered questions remain regarding the initial state of the Universe, such as flatness and horizon problems [16–19]. Moreover, there are some unresolved problems related to the origin and nature of dark matter (DM) [20–22]. Notwithstanding the excellent agreement with the observational data there still exists a number of challenging open problems associated with the late time evolution of the Universe, namely the nature of dark energy (DE) and cosmological constant problem [23–25]. Though the inflation mechanism has been introduced to treat some of the mentioned issues such as, the horizon, flatness and magnetic monopole problems at early Universe [26–28], the SCM suffers from a more fundamental issue, i.e., the initial cosmological singularity that the existence of which has been predicted by the pioneering works of Hawking, Penrose and Geroch in 1960s, known as the singularity theorems [29–37] and their later extensions by Tipler in 1978 [38,39] and by Borde, Vilenkin and Guth in the 1990s [40–45] (see also [46] for a comprehensive study). According to these theorems, a cosmological singularity is unavoidable if spacetime dynamics is described by General Relativity (GR) and if matter content of the Universe obeys certain energy conditions. A singular state is an extreme situation with infinite values of physical quantities, like temperature, energy density, and the spacetime curvature from which the Universe has started its evolution at a finite past. The existence of such an uncontrollable initial state is irritating, since “a singularity can be naturally considered as a source of lawlessness” [47]. A potential solution to the issue of cosmological singularity can be provided by “non-singular bouncing cosmologies” [48,49]. Beside a huge interest in the solutions that do not display singular behavior, there can be more motivations to seek for non-singular cosmological models. The first reason for removing the initial singularity is rooted in the initial value problem since a consistent gravitational theory requires a well-posed Cauchy problem [50,51]. However, owing to the fact that the gravitational field diverges at a spacetime singularity, we could not have a well formulated Cauchy problem as we cannot set the initial values at a singular spatial hypersurface given by t = const. Another related issue is that the existence of a singularity is inconsistent with the entropy bound S/E = (2π R)/c h¯, where S, E , R, h¯ and c being entropy, proper energy, the largest linear dimension, Planck’s constant and the velocity of light, respectively [47]. During the past decades, models which describe bouncing behavior have been designed and studied as an approach to resolve the problem of initial singularity. These models suggest that the Universe existed even before the big-bang and underwent an accelerated contraction phase towards reaching a non-vanishing minimum radius. The transition from a preceding cosmic contraction regime to the current accelerating expansion phase (as already predicted in SCM) is the so called “Big Bounce”. From this perspective, the idea that the expansion phase is preceded by a contraction phase paves a new way towards modeling the early Universe and thus, may provide a suitable setting to obviate some of the problems of the SCM without the need to an inflationary scenario. Although an acceptable model can be considered as the one being capable of explaining the issues that have been treated by inflationary mechanism, e.g., most inflationary scenarios can give the scale-invariant spectrum of the cosmological perturbations [52], problems of the SCM may find solutions in the contracting regime before the bounce occurs. The horizon problem, for example, is immediately resolved if the far separated regions of the present Universe were in causal connection during the previous contraction phase. Similarly, the homogeneity, flatness, and isotropy of the Universe may also be addressed by having a smoothing mechanism in the contraction phase, see e.g., [53–55] for more details. Moreover, though the fine-tuning is required to keep a stable contracting regime, the nonsingular bounce succeeds in sustaining a nearly scale-invariant power spectrum [56,57]. For several years, great effort has been devoted to the study of bouncing cosmologies within different frameworks. The resulted cosmological models could be obtained at a classical level or by quantum modifications. Most of the efforts in quantum gravity are devoted to reveal the nature of the initial singularity of the Universe and to better understand the origin of matter, non-gravitational fields, and the very nature of the spacetime. Non-singular bouncing solutions generically appear in loop quantum cosmology (LQC) [58– 62], where the variables and quantization techniques of loop quantum gravity are employed to investigate the effects of quantum gravity in cosmological spacetimes [63–65]. The recent large amount of works done within the loop quantum cosmology (LQC) show that when the curvature of spacetime reaches the Planck scale, the big-bang singularity is replaced by a quantum Big-Bounce with finite density and spacetime curvature [66]. Another approach, based on the de Broglie–Bohm quantum theory, utilizes the wave function of the Universe in order to determine a quantum trajectory of the Universe through a bounce [67,68]. In the framework of LQC the semi-classical Friedmann equations receive corrections as [65,69,70] ( 1 ) ( 2 ) HL2QC = 8π G 3 ρ ρ 1 − ρmax , H˙LQC = −4π G(ρ + p) 1 − 2 ρ ρmax . where, ρmax ≈ 0.41ρPl and ρPl = c5/hG2 being the Planck ¯ density [65,71–73]. We note that, the relative magnitude of ρ and ρmax enables one to distinguish classical and quantum regimes. By a short qualitative inspection of the above equations the general feature of the bouncing behavior will be revealed: initially, the Universe were in contracting phase at which the matter density and curvature are very low compared to the Planck scale. As the Universe contracts more, the maximum density is reached so that the quantum evolution follows the classical trajectory at low densities and curvatures but undergoes a quantum bounce at matter density ρ = ρmax, where we have HLQC = 0 and also H˙LQC = 4π G(ρ + p). The quantum regime then joins on to the classical trajectory that was contracting to the future. Therefore, the quantum gravity effects create a non-singular transition from contraction to expansion and thus the big-bang singularity is replaced by a quantum bounce. Furthermore, we see that for all matter fields which satisfy the weak energy condition (WEC) we have H˙LQC > 0. These two results are accounted for the general conditions for the existence of a bouncing solution. Moreover, nonsingular bouncing scenarios have been also reported in nonlocal gravity model [74] where effective Friedmann equation with a slight difference to Eq. ( 1 ) has been proposed exhibiting a bouncing-accelerating behavior. See also [75–80] for probing the issue of singularity avoidance in nonlocal gravitational theories. Another type of theories are called non-singular “matter bounce” scenarios which is a cosmological model with an initial state of matter-dominated contraction and a nonsingular bounce [81]. Such a model provides an alternative to inflationary cosmology for generating the observed spectrum of cosmological fluctuations [53–55,82–85]. In these theories some matter fields are introduced in such a way that the WEC is violated in order to make H˙ > 0 at the bounce. From Eq. ( 2 ), it is obvious that putting aside the correction term leads to negative values for the time derivative of the Hubble parameter for all fluids which respect WEC. Therefore, in order to obtain a bouncing cosmology it is necessary to either go beyond the GR framework, or else to introduce new forms of matter which violate the key energy conditions, i.e., the null energy condition (NEC) and the WEC. For a successful bounce, it can be shown that within the context of SCM the NEC and thus the WEC, are violated for a period of time around the bouncing point. In the context of matter bounce scenarios, many studies have been performed using quintom matter [86–88], Lee–Wick matter [89], ghost condensate field [90], Galileon fields [91,92] and phantom field [93–97]. Cosmological bouncing models have also been constructed via various approaches to modified gravity such as f (R) gravity [98–102], teleparallel f (T) gravity [103,104], brane world models [105], Einstein–Cartan theory [106–113], Horava–Lifshitz gravity [114], nonlocal gravity [115,116] and others [117]. There are also other cosmological models such as Ekpyrotic model [118,119] and string cosmology [120–124] which are alternatives to both inflation and matter bounce scenarios. In the present work we study the existence of the bouncing solutions in the context of f (R, T) gravity. These type of theories have been firstly introduced in [125] and later, their different aspects have been carefully studied and analyzed in [126–144]. In the current work we use an effective approach that we have introduced previously in [145]. In this method one defines an “ effective fluid” endowed with the effective energy density ρ(eff) and effective pressure p(eff) allowing thus, to reformulate the f (R, T) gravity field equations. In a class of the minimally coupled form, i.e., f (R, T) = R + h(T), one usually presumes h(T) functions and solves the resulted field equations. Instead, using the effective fluid description we obtain the h(T) function which corresponds to an EoS defined as p(eff) = Y(ρ(eff)). We therefore observe that the effective fluid picture may at least be imagined as a mathematical translation of gravitational interactions between the actual matter fields and the curvature to an overall behavior attributed to a mysterious fluid with p(eff)(T) and ρ(eff)(T). In the current article we discuss three different classes of models specified by three different effective EoSs or equivalently three different h(T) functions. We shall see that in these cases we obtain ρ(eff)(T) = β1T + β2Tγ + ρ˜(eff), p(eff)(T) = λ1T + λ2Tγ + p˜(eff), where, βi , λi , γ , ρ˜(eff) and p˜(eff) are some constants. Eliminating the trace between the effective density and pressure leaves us with an EoS as follows λ1 (ρ − ρ˜) − β1 ( p − p˜) λ2 (ρ − ρ˜) − β2 ( p − p˜) γ where the subscript “eff ” has been dropped. In view of relation ( 5 ), we may conclude that in f (R, T) = R + h(T) gravity, interactions of a perfect fluid with the spacetime curvature can be mapped effectively onto the behavior of an exotic fluid obeying equation of state ( 5 ). It is quiet interesting that a reduced form of ( 5 ) has been introduced in [146] and further studied in [147]. Such a complicated EoS has been introduced to study the cosmological implication of a model with a mixture of two different fluids, i.e., effective quintessence and effective phantom. Additionally, one can find other related works investigating some “exotic” fluids which follow various subfamily of ( 5 ). These theories may be called “modified equation of state” (MEoS) models. This branch of research presumes a mysterious fluid(s) specified by an unusual EoS with the hope of dealing with some unanswered questions in the cosmological realm. For example some relevant works in the literature can be addressed as follows; in [148] the author has employed an EoS of the form p = −ρ + γρλ in order to obtain power-law and exponential inflationary solutions. The case with λ = 1/2 has been analyzed in [147,150,151] to focus on the future expansion of the Universe. Emergent Universe models have been studied in [149] by taking into account an exotic component with p = Aρ − Bρ1/2 and in [152] with A = −1. Different cosmological aspects of DE with more simple form of EoS, i.e, pDE = α(ρDE − ρ0) have been investigated in [153] and the study of cosmological bouncing solutions can be found in [52]. Therefore, recasting the f (R, T) field equations into the “effective picture” may provide a bridge to the cosmological models supported by MEoS. Via this connection the problem of an exotic fluid turns into the problem of a usual fluid with exotic gravitational interactions. However, contrary to the former, in the latter case we start with a predetermined Lagrangian, i.e., f (R, T) gravitational Lagrangian. The importance of the effective picture becomes more clear ally chooses a mathematically suitable f (R, T) Lagrangian, then obtains the corresponding field equations, and finally tries to solve them. In [145], we introduced a novel point of view to dealing with the field equations of f (R, T) gravity. This approach is based on reconsidering the field equations in terms of a conserved effective fluid. One of the benefits of using this method is to obtain a form of f (R, T) function for a physically justified condition on the effective fluid. Therefore, instead of mathematical arbitrariness in selecting different f (R, T) functions, we have physically meaningful Lagrangian forms. Hence, in this paper we make use of those f (R, T) functions that we obtained in our previous work [145]. The action integral in f (R, T) gravity theories is given by [125] 1 2κ2 S = √−gd4x f R, T + L(m) , where R, T ≡ gμν Tμν , L(m) denote the Ricci scalar, the trace of energy momentum tensor (EMT) and the Lagrangian of matter, respectively. The determinant of metric is shown by g, κ2 ≡ 8π G is the gravitational coupling constant and we have set c = 1. We assume that the matter Lagrangian L(m) depends only on the metric components therefore the following form for the EMT can be defined Tμν ≡ − √ 2 −g δ √−gL(m) δgμν . By varying action ( 6 ) with respect to the metric components gαβ we obtain the following field equation [125] when one considers the energy conditions in f (R, T) gravity. As discussed in [154], in f (R, T) gravity the energy conditions would be obtained for effective pressure and effective energy density. Therefore, it is reasonable to define a fluid as a source with effective pressure and energy density. As we shall see, the bouncing solutions in f (R, T) gravity (using only one perfect fluid in a flat FLRW background) in the framework of our effective fluid approach, exhibit nonsingular properties such that in a finite value of the bounce time tb, non of the cosmological quantities would diverge. More exactly, as t → tb we observe that the scale factor decreases to a minimum non-vanishing value, i.e., a → ab, H t→tb→ 0, fρorte→, tnb→on ρobf,tρh(eefff)uttu→r etb →sin0gualnadritpie(esffw)to→ultbd→appp(eeaffr).bA.Tlshoe,reinall cases we have W t→tb→ −∞, where W being the effective EoS parameter and subscript “b” stands for the value of quantities at the time at which the bounce occurs. We then observe that if we want to describe nonsingular bouncing solutions in f (R, T) = R + h(T) gravity using a minimally coupled scalar field, a phantom field should be employed. These solutions show a violation of the NEC in addition to the strong energy condition SEC. Such a behavior is predicted in GR for a perfect fluid in FLRW metric with k = −1, 0 [155]. The current research is planned as follows. In Sect. 2 we briefly present the effective fluid picture. Sect. 3 is devoted to the bouncing solutions with asymptotic de Sitter behavior before and after the bounce. We first analyze models with constant effective pressure in Sect. 3.1, they are called models of type A. We then proceed to investigate the corresponding bouncing solutions, the energy conditions, the scalar field representation and finally the stability of these type of solutions. In Sect. 3.2 models which correspond to two different EoSs assuming p(eff) = Y(ρ(eff)) will be discussed. These models are named as B, C, D and E models. An example of the matter bounce solution is considered in Sect. 4 which is labeled as model E. The connection of A-E models with MEoS theories will be presented through the effective picture. Section 5, is devoted to study of scalar-type cosmological perturbations. In Sect. 6 we give a brief review of singular models in the context of f (R, T) gravity and obtain a class of solutions exhibiting singular behavior. Finally, in Sect. 7 we summarize our conclusions. 2 Reformulation of f (R, T) field equations in terms of a conserved effective fluid In the present section we review the field equations of f (R, T) gravity theories and rewrite them in terms of a conserved effective fluid. This reformulation would allow us to better understand the properties of bouncing solutions as well as classifying them. In f (R, T) gravity, one usu( 6 ) ( 7 ) ( 8 ) ( 9 ) gμν − μ ν F (R, T) 1 F (R, T)Rμν − 2 f (R, T)gμν + = κ2 − F (R, T) Tμν − F (R, T) μν , where we have defined μν ≡ gαβ δTαβ δgμν = −2Tαβ + gαβ L(m), and F (R, T) ≡ ∂ f (R, T) ∂T and F (R, T) ≡ ∂ f (R, T) ∂R where a(t ) is the scale factor of the universe and d 2 is the standard line element on a unit two sphere. Applying metric ( 11 ) to field equation ( 8 ) together with using the definition for μν for a perfect fluid leads to where “ ” denotes derivative with respect to the argument. Also, from Eq. ( 14 ) we obtain f (R, T) − F (R, T)R + 3F˙ (R, T)H 1 3H 2 F (R, T) + 2 = κ2 + F (R, T) ρ + F (R, T) p, as the modified Friedmann equation, and 2F (R, T)H˙ + F¨ (R, T) − F˙ (R, T)H = − κ2 + F (R, T) (ρ + p), as the modified Raychaudhuri equation, where H indicates the Hubble parameter. Note that, we have used L(m) = p for a perfect fluid in the expression ( 9 ). Applying the Bianchi identity to the field equation ( 8 ) gives following covariant equation (κ2 + F )∇μTμν + 21 F ∇μT + Tμν ∇μF − ∇ν ( pF ) = 0, where we have dropped the argument of F (R, T) for abbreviation. Equation ( 14 ) tells us that the conservation of EMT is not generally respected in f (R, T) gravity theories. There is Only a narrow class of solutions for which the conservation of energy is still preserved [137]. In this work we consider a specific class of models for which the Lagrangian is written as a minimal coupling between the Ricci curvature scalar and a function of trace of EMT, i.e., f (R, T) = R + ακ2h(T). Also, we study bouncing solutions when the Universe contains a single perfect fluid with a barotropic EoS given as p = wρ. In the next section, we seek for non-singular cosmological solutions in f (R, T) gravity with a baratropic perfect fluid. As we shall see, a useful approach for choosing the functionality of h(T) is to rewriting Eqs. ( 12 ) and ( 13 ) in terms of an effective fluid that respects the conservation of EMT. For the Lagrangian, Eqs. ( 15 ), ( 12 ) and ( 13 ) are simplified to 3H 2 = κ2 1 + (1 + w)αh T αh 3w − 1 − 2 , and 2H˙ = −κ 2 w + 1 (1 + αh )T, 3w − 1 ( 12 ) ( 13 ) ( 14 ) ( 15 ) ( 16 ) ( 17 ) α 1 + 2 (3 − w)h + α(1 + w)T h +3H (1 + w) 1 + αh T = 0. T ˙ As a matter of fact, one can directly solve Eqs. ( 16 )–( 18 ) to obtain the scale factor or the Hubble parameter, once the h(T) function is determined. Alternatively, one can also define the pressure and energy density profiles of an effective fluid along with imposing the energy conditions in order to obtain the functionality of h(T). We then proceed in this way and rewrite Eq. ( 16 ) as ( 18 ) ( 19 ) ( 20 ) ( 21 ) ( 22 ) 3H 2 = κ2ρ(eff)(T), where, ρ(eff)(T) ≡ 1 + (1 + w)αh T αh 3w − 1 − 2 With this definition the acceleration of expansion of the Universe can be obtained as follows a¨ κ2 a = − 6 ρ(eff)(T) + 3 p(eff)(T) , where, we have defined the effective pressure as w αh p(eff)(T) ≡ 3w − 1 T + 2 . Therefore, the original field equations of f (R, T) gravity can be recast into the usual Friedman form with an effective fluid. This fluid is characterized by an effective energy density and effective pressure which in turn are determined in terms of the EMT trace. Therefore, once a property/relation for energy density and pressure components is established, a first order differential equation for the function h(T) could be reached. Solving the differential equation leads to an h(T) function that conveys the specified property/relation. Hence, in this way we can obtain a minimal f (R, T) model based on the conditions on energy density and pressure profiles of the effective fluid, instead of choosing the functionality of f (R, T) based on ad hoc mathematical terms. It is straightforward to verify that Eq. ( 18 ) is turned to the usual conservation equation in terms of ρ(eff) and p(eff), i.e., ρ˙(eff) + 3H (ρ(eff) + p(eff)) = 0, ( 23 ) where the arguments are dropped for the sake of simplicity. Consequently, an effective EoS parameter can be defined for this effective fluid as 3 Asymptotic de-Sitter bouncing solutions in f (R, T) gravity In this section we study different bouncing cosmological solutions of f (R, T) gravity. We extract those solutions which correspond to some properties of the effective fluid. Such an approach may help us to understand how these solutions can emerge in f (R, T) gravity. Furthermore, there can be obtained more bouncing solutions, however, to show that f (R, T) gravity theories are capable of describing a nonsingular pre-Big Bang era, we are restrict ourselves to study only few examples. We set κ2 = 1 in the rest of the work. 3.1 Type A models: solutions which correspond to a constant effective pressure, p(eff)(T) = P Let us begin with bouncing solutions which are obtained by assuming an effective fluid with constant pressure. We show that these type of models lead to a de-Sitter era at late times [145]. From definition ( 22 ) for a constant effective pressure we obtain 2 hA(T) = α w P + 1 − 3w T . Substituting the above function into the modified conservation equation ( 18 ) together with using T = (3w − 1)ρA for a perfect fluid we get ρA = ρ0aA−3. where ρ0 is an integration constant. Substituting solutions ( 25 ) and ( 26 ) back into the Friedman equation ( 19 ) leaves us with the following differential equation for the scale factor 3 aA(t ) 2 aA(t ) ρ0 w2 − 1 − (3w − 1)aA(t )3 + P = 0. For P = 0, a particular solution to the above equation is found as aA(t ) ∝ t 2/3. Moreover, in case in which the middle term is absent, the solution represents a de-Sitter phase for P < 0. From the Friedman equation ( 27 ), we observe that the effective pressure P plays the role of a cosmological constant. Other forms for the Friedman equation similar to Eq. ( 27 ) have been found in the literature. For example in [153] authors have worked on a perfect fluid, dubbed as DE with a linear EoS, in a flat FLRW background. They used pDE = α(ρDE − ρ0) ≡ pα + p , where α and ρ0 × B + (w2 − 1)Pρ0 sinh √ 2(3w − 1) + cosh √ −3Pt − 1 2/3 A = P2 2B − w2 − 1 Pρ0 3w − 1 3 2 B = a0 P − a0 P 3 3 3 a0 P + P t − 3 −3Pt 1/3 ρ0(1 − w2) 3w − 1 being some constants. The DE pressure is then decomposed to a dynamical part which corresponds to pα as well as a constant part which corresponds to p . These assumptions lead to an equation that will reduce to ( 27 ), when only p is considered. In the upcoming subsections when we study the model C, we find that in the context of f (R, T) gravity one can obtain the same results as those given in [153]. Equation ( 27 ) admits three different general solutions which we will consider only one of them. A general solution of Eq. ( 27 ) is obtained as 1 aA(t ) = A cosh P t − 3 − sinh Note that once we set the integration constant in Eq. ( 27 ) such that aA(0) = 0, we get B = 0 and furthermore, if we consider the case w = 0, we obtain A = Pρ01/3. Thus the solution ( 28 ) reduces to the familiar form a(t ) = (−ρ0/P)1/3 sinh2/3(√−3P/2t ). However, the integration constant in solution ( 28 ) is fixed such that we have a(t = 0) = a0. The most important feature of solution ( 28 ) is the appearance of cosine hyperbolic function which allows us to have a non-singular behavior for P < 0. This solution describes the de-Sitter expansion of the Universe which is initially dominated by dust [156]. In Fig. 1 the thick black curve shows the scale factor solution ( 28 ). This solution behaves exponentially in the far past and far future from the bounce and reaches a nonzero minimum value at the bounce time, given by a(min)A = A−1 3 2(1 − 3w)B (1 − w2)Pρ0 − 1 2(3w −1)B2 × 2(3w −1)B − w2 − 1 Pρ0 − w2 −1 Pρ0 3w − 1 The Hubble parameter is obtained as HA(t) P ⎡ 2(1 − 3w)B+ w2 −1 Pρ0 sinh √−3Pt +cosh √−3Pt + 1 ⎤ = − 3 ⎣ 2(3w − 1)B + w2 −1 Pρ0 sinh √−3Pt +cosh √−3Pt − 1 ⎦ . 2/3 Fig. 1 Different cosmological quantities related to the model A which are characterized by solution ( 28 ). The black solid line indicates the scale factor. The black middle sized dashed line shows the Hubble parameter, the solid blue curve indicates the effective energy density, the black dot-dashed curve shows the pressure, The long-dashed blue curve shows the effective EoS and the black dotted curve shows the effective energy density. We see that WA < −1 always, signaling that the scalar representation of model A should be defined by a phantom field. We have set w = 0.6, ρ0 = 1, P = −2.8 and a0 = 1.12. The horizontal axis represents the cosmic time Figure 1 shows that the Hubble parameter tends to constant values before and after the bounce and vanishes at the bounce. On the other hand, the Hubble parameter and its time derivatives diverge at an imaginary time ts = i π/√−3P − tb, where tb is the (real) time at which the bounce occurs. Therefore, we observe that the Hubble parameter behaves in a well-defined way (without any future singularity) so that it respects the bounce conditions. We can also obtain the effective energy density and the effective EoS parameter. To this aim, we substitute ( 25 ) into Eqs. ( 20 ) and ( 24 ), respectively. Thus we obtain w2 − 1 ρ0 ρ(eff)A(a) = (3w − 1)aA3 − P, 3 (1 − 3w)PaA WA = (3w − 1)PaA3 − w2 − 1 ρ0 . As can be seen from Fig. 1, the effective density reduces from a constant value and tends to zero near the bounce. From Eq. ( 19 ) we see that the vanishing of Hubble parameter at the bounce demands that the effective density becomes zero. Also for the same reason, the effective EoS diverges at the bounce. Such behaviors are common for all bouncing models that we shall present in the framework of minimal f (R, T) gravity. the matter energy density itself increases from small values to a maximum value near the bounce. Based on the exchange of energy between gravitational field and matter constituents (that the mechanism of which is explained in [128]) one may explain the process of bouncing behavior; the interaction of the real fluid with curvature leads to some transformations of energy from gravitational field to matter before the bounce where the spacetime curvature is dominant ( 31 ) ( 32 ) in comparison to matter energy density. Such a transmutation, that the start of which is triggered at the time far past the bounce, gives rise to an increase in the energy density as the bounce event is approached. At the bounce time the energy density of matter grows to a maximum value after which the process of transmutation is reversed until the density falls back to zero (the post-bounce regime). Note that the effective energy density remains constant in the de-Sitter era. However, some physics is needed in order to explain the process of matter production from curvature component which disturbs the stability of de-Sitter era to enter the bounce event. Informations from Fig. 1 can help us to discuss the energy conditions. In GR the well-known energy conditions are the NEC, WEC, SEC and the dominant energy condition (DEC). In a modified gravity theory with defined effective energy density and pressure, these conditions can be written as [154] WEC ⇔ ρ(eff) ≥ 0, and ρ(eff) + p(eff) ≥ 0, NEC ⇔ ρ(eff) + p(eff) ≥ 0, SEC ⇔ ρ(eff) + 3 p(eff) ≥ 0 and ρ(eff) + p(eff) ≥ 0, DEC ⇔ ρ(eff) ≥ 0, and ρ(eff) ± p(eff) ≥ 0. For the model A, we obtain the following results w2 − 1 ρ0 ρ(eff) + p(eff) = (3w − 1)aA3 , w2 − 1 ρ0 ρ(eff) + 3 p(eff) = (3w − 1)aA3 + 2P. ( 33 ) ( 34 ) ( 35 ) ( 36 ) ( 37 ) ( 38 ) It is obvious that the fulfillment of NEC and thus WEC requires (w2 − 1)/(3w − 1) ≥ 0 which gives −1 ≤ w ≤ 13 . From ( 31 ) we see that the effective energy density tends to −P at late times and vanishes at the time of bounce. Therefore, we obtain 3P ≤ ρ(eff) + 3 p(eff) ≤ 2P, (remember that P < 0). Thus, since only negative values are valid for the effective pressure P, the SEC is always violated. However, the validity of NEC depends upon the value of w. We plot the diagrams for w = 0.6 in Fig. 1. This figure shows that the NEC and SEC are violated in this case. Our studies show that the bouncing behavior is achieved from solution ( 28 ) for w > 1/3. Note that, as we have mentioned before, solution ( 28 ) is only one of the three possible solutions of Eq. ( 27 ). Investigating other solutions may validate the cases with w < 1/3 from energy conditions point of view. One may be tempted to reinterpret the source of matter as that of a scalar field. Such a representation is also used in similar works [52,148,152,157]. In the case of constant effective pressure, the particular solution ( 28 ) (which is valid for 1/3 < w < 1) corresponds to WA < −1 which can be realized from relation ( 32 ), see also the long-dashed blue curve in Fig. 1. Therefore, if we want to translate mutual where, the subscript Ph denotes the phantom field and “;” indicates covariant differentiation. We then get φ˙(2eff)A = −(ρ(eff) + p(eff)) = (13w−−w21)aρA30 , 1 V (φ(eff)A) = 2 w2 − 1 ρ0 ρ(eff) − p(eff) = 2(3w − 1)aA3 − P. ( 42 ) A straightforward calculation reveals that φ(eff)A = −4 × arctan ⎢⎢ ⎣ P A (1 − 3w)(w − 1)A 3C ⎡ B(3w−1) tanh √−(w + 1)Pρ0C −136p t − 1 + w2 −1 P, ρ0 ⎤ ⎥⎥ , ⎦ where we have defined C = (w − 1) 2(1 − 3w)B + w2 − 1 Pρ0 . We can also obtain φ(eff)A in terms of the scale factor by solving solution ( 28 ) for time t , which gives 1 t = √−3P log ⎢⎣⎢⎢⎢⎢ ⎡ ⎛ ⎞ ⎤ ⎢ 2aA3A3(1−3w)2 ⎝ 1− aCA3(Aw3+(11−)P3wρ)02 +1⎠ +(w + 1)Pρ0C ⎥ w2 − 1 2 P2ρ02 ⎥⎥⎦⎥⎥⎥ . We have plotted φ(eff)A and V(eff)A in Fig. 2 for the same parameters of Fig. 1. Another important issue that needs to be treated is to investigate the stability properties of solution of Eq. ( 27 ). Substituting solution ( 31 ) in ( 19 ) and also ( 26 ) in ( 17 ) and taking H and ρ as dynamical variables, we arrive at the following dynamical system β H˙ = 2 ρ , ρ˙ = −3Hρ , 3H 2 = −βρ − P, interactions of perfect fluid and curvature as the behavior of an effective scalar field, we should employ a phantom scalar field. Therefore, we can define 1 LPh = − 2 φ(eff);μφ(eff);μ − V (φ(eff)), 1 2 ρ(eff) = − 2 φ˙(eff) + V (φ(eff)), 1 2 p(eff) = − 2 φ˙(eff) − V (φ(eff)). ( 39 ) ( 40 ) ( 41 ) ( 43 ) ( 44 ) ( 45 ) ( 46 ) ( 47 ) ( 48 ) where β = (1 − w2)/2(3w − 1) and we have dropped subscripted A . Note that the validity of solution ( 28 ) requires that β > 0. The system ( 46 )–( 48 ) has two critical points far from the bounce; P(±) = (±√−P/3, 0) which corresponds to the eigenvalues λ(±) = (∓√−P/3, 0). We see that the stability properties of solutions are independent of w. Due to the appearance of the zero eigenvalues, one may not decide about the stability properties of these fixed points, however, by inspecting Eqs. ( 46 ) and ( 47 ) it is possible to figure out the nature of the fixed points. For fixed point P(−) Eq. ( 47 ) becomes ρ˙(−) = √−3Pρ(−). Therefore, in the vicinity of P(−) within the phase space, by indicating the values of ρ and H on vertical and horizontal axises, respectively, we have H˙ > 0 & ρ˙ > 0 for all points with ρ > 0 and H˙ < 0 & ρ˙ < 0 for points with ρ < 0. These show that when ρ → ρ + δρ for t → t + δt , the solution at P(−) will not stay stationary and hence it is a repulsive fixed point. On the other hand, for P(+) we have ρ˙(+) = − √−3Pρ(+). ( 49 ) ( 50 ) Therefore, in this case we have H˙ > 0 & ρ˙ < 0 for all points with ρ > 0 and H˙ < 0 & ρ˙ > 0 for points with ρ < 0 in the vicinity of P(+). Hence it is a stable fixed point. The bounce corresponds to the point P(b) = (0, −P/β) for which we have H˙ = −P/2 > 0. Therefore, at this point the tangent vectors on the phase space trajectories are directed toward the right side. We plot a typical trajectory in the phase space in Fig. 3. 3.2 Solutions which correspond to a general effective EoS, p(eff) = Y(ρ(eff)) These class of models can be constructed by imposing a particular condition on the effective profiles. This approach can Fig. 3 A typical phase trajectory in the phase space plane (H, ρ) corresponding to model A. This plot shows that the evolution of the Universe has started from an unstable state and terminates in a stable phase. The bounce event is an unstable (of saddle type) fixed red point be viewed as a sort of classification of f (R, T) gravity models based on the properties of the effective quantities. Generally, one can obtain a class of h(T) functions for a determined property which is specified by an effective EoS. In the following sections we consider two subclasses based on conditions on the effective densities. We find that each class of h(T) solutions that exhibit bouncing behavior correspond to an effective EoS which is already introduced or obtained for an exotic fluid in the literature [52,152,153,157,158]. 3.2.1 Type B models: solutions which follow the relation dρ(eff)/dT = [n/(1 + w)T](ρ(eff) + p(eff)) Applying this condition together with using the definitions for the effective energy density and pressure, we arrive at a differential equation for h(T) which can be solved as follows [137] 2 B(w + 1) 2n+3w−1 2(n − 1) hB(T) = 2n + 3w − 1 T 2(w+1) − α(2n + w − 3) T + where, B and B are integration constants and n is an arbitrary constant. Substituting the relation dρeff/dT = [n/(1 + w)T](ρeff + peff) into Eq. ( 18 ) gives ρB = ρ0a− 3(wn+1) . Next we proceed to find a non-singular bouncing solution by solving the modified Friedmann equation ( 16 ) for solutions ( 51 ) and ( 52 ). We first try to obtain solutions of the form aB(t ) = R(cosh[(t − t0)/R] − S) where R and S are constants. This type of solution has been discussed in [52] under assumption of an MEoS and has the following form of the Friedmann equation 3 aB(t ) 2 aB(t ) − 3(S2 − 1) aB(t )2 6S 3 − RaB(t ) − R Applying ( 51 ) and ( 52 ) in Eq. ( 16 ) gives the Friedmann equation for arbitrary constants w and n. We can check that, there are only two cases which correspond to Eq. ( 53 ) and thus to the scale factor aB as the solution; when w = −1/5, n = 12/5 and w = −1/5, n = 6/5. However, the latter leads to similar physics to the former. The physical quantities constructed out of the bouncing solution for w = −1/5, n = 12/5 are given as follows aB(t ) = R cosh − S , t R 6 − α B , R = 3 S = − 10 Rρ0, B, ( 51 ) ( 52 ) HB (t ) = , 3 aB + R(S − 1) aB + R(S + 1) 2 2 aBR R R − S(4aB + SR) − 3aB2 2 2 aBR The behavior of the above quantities are depicted in Fig. 4. In order that the ansatz aB(t ) satisfies Eq. ( 53 ) we must have 1 B = 32α2 27 B . Solution ( 54 ) shows that the Universe shrinks from an infinite size to a minimum radius which is equal to R(1 − S). The size of Universe at the time of bounce is controlled by the constant B as well as the coupling constant α. As can be seen, in Eq. ( 53 ) and the subsequent solutions, these two constants appear as a multiplied form. This means that if either of these constants become zero, the bouncing solution would disappear. Also, expression ( 15 ) together with solution ( 54 ) ( 55 ) ( 56 ) ( 57 ) ( 59 ) ( 51 ) indicate that the bouncing solution disappears in model B when α = 0. The expressions for effective energy density and pressure, i.e. ( 56 ) and ( 57 ), can be rewritten as a sum of three densities and pressures given by From these expressions, we could suppose that our effective fluid consists of a combination of three perfect fluids with EoSs wi = ρi / pi . In this view, there is a DE component which corresponds to a cosmological constant with w1 = −1, a quintessence with w2 = −2/3 and a fluid which drives an expanding Universe with zero acceleration, with w2 = −1/3. Such a description has been presented in [52]. In the framework of f (R, T) gravity this decomposition may be translated as follows; the effects of unusual interactions of matter (here a perfect fluid with w = −1/5) with curvature, could produce the same behavior as the case of GR with three different fluids. Eliminating the scale factor from solutions ( 56 ) and ( 57 ) leads to an effective EoS of the form p(eff)B = − × 3R2 ρ(eff)B 2S 3 + 3R2 S2 − 1 1 S2 − 1 ρ(eff)B + 9 2 + R In the cosmological applications, a general type of EoS p = βρ + γ f (ρ) is ascribed to some exotic or dark fluid which can determine the evolution of the Universe. Different choices for function f (ρ) are studied in the literature. Such an equation has been already discussed in [158] where the authors considered cosmological consequences of an EoS of the form p = −ρ − ρq + 1. Substituting expressions ( 56 ) and ( 57 ) in the NEC and SEC conditions ( 34 ) and ( 35 ) leads to the following conditions However, expressions ( 63 ) and ( 64 ) are never satisfied; far from the bounce where a → ∞, the second term of ( 63 ) which always has negative sign dominates and in the limit a → ab = R(1 − S) the expression for NEC becomes −2/R2(1 − S) which is also negative because S < 0. The same line of reasons can be used to prove the violation of SEC. Thus, in model B the NEC and SEC are never satisfied. Note that, the NEC is violated only near the bounce, because expression ( 63 ) tends to zero as the scale factor gets large values. By considering the above discussion for energy conditions, we find that again, a phantom scalar field can be used for modeling the behavior of bouncing solution. In the case of model B we have 2 φ˙(eff)B = − V(eff)B = 2 aBS + S2 − 1 R , 2 aBR Our studies show that the behavior of the above solutions are similar to those of model A which in Fig. 2 a typical example is demonstrated. By choosing HB and ρB as dynamical variables one can rewrite the Friedmann equation as follows ( 65 ) ( 66 ) ( 67 ) ( 68 ) ( 69 ) ( 70 ) H˙ = 1 −ρ02S2 ρ2 + 130 ρ , ρ˙ = −Hρ , H 2 = 1 −ρ02S2 ρ2 − 53 ρ + 1 2 , R where we have dropped the subscript B. This system admits two critical points PB(±) = (±√1/R, 0) far from the bounce. In these situations we have ρ → 0, hence, the dynamics of Eq. ( 67 ) is determined by the second term. Also, near the point of the bounce as specified by Hb = 0, ρb = 3ρ02/10R(R + 1), we always have H˙ > 0. Therefore, the stability properties of the system is similar to the model A. The Friedmann equation ( 16 ) for general function ( 51 ) takes the form 3H 2 2αn(w + 1) B = (3w − 1)(2n + 3w − 1) n(w − 1)T + (3w − 1)(2n + w − 3) − T wn−+21 + 23 α B 2 . Seeking for a general solution demands that one substitutes the solution ( 52 ) into ( 70 ) (using the fact that T = (3w − 1)ρ) and solves for the resulting differential equation to find the scale factor. However, the resulting equation cannot be solved analytically for arbitrary values of w and n. Nevertheless, for particular values of these parameters a non-singular solution can be obtained as given in expressions ( 54 )–( 58 ). But, we are still able to find more general solutions. As the third type of solutions named as C, we work on a bouncing solution for which the governing differential equation is given by C = aC(t ) HC(t ) = − ( 72 ) 2n 3n+3 , ( 73 ) . ( 74 ) = a0 (Q + 1) cosh 3α C (n + 1) 2 2n t ± cosh−1(Q) 2 where Q is an arbitrary constant and for the Hubble parameter we have α C tanh 6 By substituting ( 52 ) for w = 1 within the definitions ( 20 ), ( 22 ) and ( 24 ) along with using relation ( 72 ) we get the effective quantities as follows1 Fig. 4 Cosmological parameters of the bouncing solution for the model B. We have set α = 20, B = −0.3 and ρ0 = 1 3 aC(t ) aC(t ) 2 where we have set w = 1. Choosing a relation between and ρ0, a0, n and C as the following C (n + 1) C √2ρ0a0−3/n −n−1 n(Q + 1) , Equation ( 71 ) leads to the following solution for the scale factor ( 75 ) ( 76 ) ( 77 ) ( 78 ) 1 We note that as we are concerned with the solutions that respect the conservation equation ( 18 ), expression ( 52 ) is used for these class of solutions. 2 In the original paper they used α and ρ0 instead of γ and θ , respectively. We changed this character in the present work to prevent ambiguity. Fig. 5 The cosmological quantities related to the bouncing solution of the model C for, n = 3, a0 = 1.3 ρ0 = 2, C = 3, α = 1 and Q = 0.5 ρ(eff)C(a) = α p(eff)C(a) = −α WC(a) = 1 C ⎢⎢ 2 − ⎣ ⎡ C ⎢ ⎡ ⎢ ⎣ − 3(nn+1) ⎤ a a0 Q + 1 − 3(nn+1) a a0 n(Q + 1) 2 + n(Q + 1) n 2 − (Q + 1) a a0 a a0 ⎥⎥ , ⎦ ⎤ We have plotted the cosmological parameters ( 73 )–( 77 ) in Fig. 5. Unlike the solution B, here, we see that matter creation in the time of bounce leads to a decrease in the effective pressure. Eliminating the scale factor between expressions ( 75 ) and ( 76 ) gives the effective EoS which can be viewed as the characteristic equation for the model C. We therefore get 1 p(eff)C = n ρ(eff)C − α 2 C (1 + n) . Some of the cosmological properties of model C have been investigated in [153]. The authors have considered a model of DE for which an EoS of the form pDE = γ (ρDE − θ )2 is assumed for a perfect fluid. We therefore observe that if we apply ρ(eff) → 3ρDE within Eq. ( 19 ), we obtain the same Friedman equation as the one given in [153]. Also, by redefining the parameters as n → 1/γ and α(1 + n) C/2 → θ in ( 73 ) we will obtain the corresponding solution for the scale factor. These considerations show that the problem of dark fluid with an unusual EoS (which may not clearly correspond to a definite Lagrangian) can be explained in the framework of f (R, T) gravity. 3.2.2 Type D models: solutions which are consistent with the relation dρ(eff)/dT = m Applying this condition on the definition of effective energy density, i.e., definition ( 20 ), leads to a second order differential equation for h(T) function; the solution then reads hD(T) = 3w − 1 + D, 2 D(w + 1) 3w−1 T 2(w+1) + 2[m(1 − 3w) + 1] T α(w − 3) ( 79 ) where, m is an arbitrary constant and D, D are integration constants. Substituting ( 79 ) into the conservation equation ( 18 ), we obtain a first order differential equation for the matter energy density in terms of the scale factor. However, since the mentioned equation cannot be solved for an exact general solution for arbitrary values of w and m, we proceed with particular cases. Note that, further investigations may give other exact solutions or even numerical simulations can be utilized to study other solutions. At the present, we work on a particular case w = 1. The conservation equation ( 18 ) then yields ρD(a) = ζ − α DaD3 2 8aD6m2 , where we have again used T = (3w − 1)ρ. In the limit a → ∞ one obtains D2α2/8m2 from solution ( 80 ). As can be seen, in the model D the matter energy density evolves from a non-vanishing initial value far from the bounce event. This property cannot be seen in the previous models. To obtain the solution for the scale factor, we substitute ( 79 ) in the Friedmann equation ( 19 ) for w = 1, which gives Dζ ω − ω cosh 3αω t − ϒ sinh 4m 3αω t 4m ζ = α D ± 2m#2ρ0, ω = α D2 − 2m D, ζ − 2αa03 D , From solution ( 81 ), we can obtain the time at which the bounce occurs as well as the radius of the Universe at the moment of bounce, as aD(t) = where ϒ = = 6 ζ a0 + αω Dζ − ωa03. t D(b) = log − ωϒ + ωϒ , ⎡ $ 1 %% ζ − α DaD3 2 p(eff)D = 4 ⎢⎣ 2α D& m2a6 D ( 85 ) , ( 86 ) ( 87 ) ( 88 ) ( 89 ) ( 80 ) WD = In view of what we discussed in the paragraph right after Eq. ( 62 ), model D corresponds to a generalized EoS with β = 1. This may be interesting since the most bouncing models have been obtained for β = −1. Note that, the behavior of obtained cosmological quantities for two values of ζ as given in ( 82 ), are similar to those of model C. By choosing the model parameters so that the term including si nh in ( 81 ) disappears, we arrive at a different model. In this case, the behavior of cosmological quantities is the same as the case for which ϒ = 0. However, the evolution of matter energy density and the effective pressure are different. We typically plot these quantities for both situations in Fig. 6. The thick curves belong to the solution ( 81 ) and the thin ones show the solution with ϒ = 0. The energy condition considerations show that models of type D lead to the violation of NEC near the bounce event. Also the phase space and the scalar field representation of this model are the same as model A. Fig. 6 Mutual behaviors of matter energy density and effective pressure related to the bounce solution in model D. The thin curves correspond to solution ( 81 ) when ϒ = 0. We have set α = −2.7, D = −0.1, m = −0.8, ρ0 = 2, a0 = 1.2, D = 1.2 for thick curves and α = −3, m = 0.2, ρ0 = 0.1, a0 = 1.1, D = −0.2 for dashed curves 4 Matter bounce solutions in f (R, T) gravity In this section, we deal with the well-known matter bounce scenario which can be established through the models that obey ( 51 ) and ( 52 ). We specify these class of solutions as type E models. A branch of the matter bounce scenarios have been discussed with the characteristic scale factor aE(t ) = Qt 2 + Z M , where Q, Z and M are positive constants. Note that Q = uρmax, with u = 2/3, 3/4, 4/3, Z = 1 and M = 1/3 [100,101,103,104] and M = 1/4 [159] has been used in the literature. The scale factor ( 91 ) gives the following expression for the Hubble parameter ( 91 ) ( 92 ) 2MQt HE(t ) = Qt 2 + Z . Note that in model E the Hubble parameter and its time derivatives never diverge since all of them are proportional to negative powers of Qt 2 + Z. By eliminating the time parameter between ( 91 ) and ( 92 ), we get the Hubble parameter in terms of the scale factor and thus the Friedmann equation can be obtained as 1 2 3HE2(a) = 12QM2 a− M − ZaE− M . E From another side, substituting ( 51 ) into ( 16 ) together with using ( 52 ), the Friedmann equation in terms of the scale factor reads ( 94 ) ( 95 ) ( 96 ) ( 97 ) ( 98 ) for ( 99 ) ( 101 ) ( 102 ) 3H 2(a) = As can be seen, in the context of f (R, T) gravity, there can be found a matter bounce solution for every values of w (except for w = 1/3, 1). The value w = −1 can be accessed for large values of M (see relations ( 95 )). From ( 52 ) and ( 95 ) we see that for models of type E we have ρE = ρ0aE−1/M. The effective quantities are then obtained as ρ(eff)E = p(eff)E = − 12M2Q Z 2 a− M E 4MQ a− M2 Z E 1 aEM − 1 = 12M2Q2t 2 Z Qt 2 + 1 2 , ( 100 ) 1 (3M − 1)aEM − 3M + 2 = − 1 WE = 3 4MQ (3M − 1)Qt 2 + 1 Z Qt 2 + 1 2 1 M(1 − a M ) 1 = −1 + 3M 1 1 − Qt 2 . 1 + M − 3 ( 93 ) Model E corresponds to an effective EoS as p(±eff)E = 2 2Q 3M − 1 ρ(eff)E ± Z M2 − 3ZQ ρ(eff)E 2MQ − Z ter energy density for model E. The model parameters have been set as, M = 0.6, α E = −0.011 and Z = 1 Estimating effective pressure ( 101 ) in the limiting times t → 0 and t → ±∞, indicates that the effective EoS follows the expression p(−eff)E far from the bounce and obeys p(+eff)E near the bounce. In t → 0 pressure ( 101 ) gives p(eff)E = −4MQ/Z which is consistent with p(−eff)E and in t → ±∞ we have p(eff)E = 0 which can be explained only with p(+eff)E. Figure 7 shows the behavior of different quantities. As is seen in the left panel, the scale factor decreases till reaching a minimum non-zero value at t = tb where the Hubble parameter vanishes. The Universe experiences four phases during its evolution from pre-bounce to post bounce. Before the bounce occurs, the Universe has been within an accelerated contracting regime till the first inflection point (t1inf < tb) is reached at which the accelerations vanishes. At this point the Hubble parameter maximizes in negative direction and correspondingly the effective energy density gains a peak value. The Universe then enters a decelerating contracting regime so that its velocity decreases in negative direction. The collapse of Universe halts at the bounce time after which the Universe goes into an accelerating expanding phase where both a¨E > 0 and HE > 0. Once the second inflection point (t2inf > tb) is reached, the Universe enters a decelerating expanding regime so that the speed of expansion decreases at later times. Let us now check the behavior of NEC(eff) and SEC(eff) which for model E take the following form as NEC(eff)E = p(eff)E + ρ(eff)E 4MQ Qt 2 − 1 = 2 , Z Qt 2 + 1 SEC(eff)E = p(eff)E + 3ρ(eff)E = 12MQ (1 − 2M)Qt 2 − 1 Z Qt 2 + 1 2 . It is obvious that in the expression ( 104 ), the sign of Qt 2 − 1 determines the validity of NEC. We therefore find that near the bounce, NEC is violated within the range −1/√Q < t < 1/√Q, and out of this range it is preserved. The SEC is violated within a larger time interval, i.e., in −1/√Q(1 − 2M) < t < 1/√Q(1 − 2M) for 0 < M < 1/2 and is always violated for M > 1/2. In Fig. 8 we have plotted the expression of NEC and SEC for two different values of M, i.e., for M = 0.3 and M = 0.7. We see again that, in the background of f (R, T) gravity a bouncing behavior corresponds to violating the NEC. We then may conclude that scalar representation type E models can be constructed by a phantom field within the time interval where NEC is violated (notice that from solution ( 102 ) we see that for the same intervals in which the NEC ( 104 ) ( 105 ) ( 112 ) ( 113 ) ( 114 ) ( 115 ) ( 116 ) ( 117 ) ( 106 ) ( 107 ) ( 108 ) ( 109 ) ( 110 ) passing through an unstable bounce phase and finally reaches a stable state. 5 Stability of the bouncing models In this section we verify the wholesomeness of the bouncing models which has been introduced in the previous sections through considering the possibility of occurring serious instabilities. We therefore examine the evolution of scalartype perturbations in the discussed models within the metric formalism. Since f (R, T) gravity introduces unusual coupling of matter to curvature part of its action, the evolution of matter density perturbations (specially, as the effect of bounce event on the evolution of matter perturbations is not obvious) can be problematic. In order to study such type of perturbations we consider the matter density perturbations in f (R, T) = R + ακ2h(T) models for a fat FLRW metric in the longitudinal gauge ds2 = −(1 + 2 )dt 2 + a(t )2( 1 − 2 )δi j d xi d x j , (111) where the metric scalar perturbations and are functions of four coordinates (t, x , y, z), generally. In the current work we shall obtain necessary equations for models including a barotropic perfect fluid with equation of state p = wρ and a general h(T) function. In this respect, the authors of [126] have already considered the matter perturbations in a narrow class of f (R, T) models3 for a pressure-less perfect fluid. The perturbations of EMT in the longitudinal gauge are given by [160] δTtt = −δρm , 1 δTit = a (1 + w)ρv,i , δTti = −a(1 + w)ρv,i , δTij = wδi j δρ , where, v is a covariant velocity perturbation [161]. Using the background Eqs. ( 16 ) and ( 17 ) we obtain the following equations for the scalar perturbations in Fourier space The analysis of these equations is similar to those of the previous models. Firstly, note that only for 1/6 < M < 1/2, equation of ( 110 ) leads to the standard Friedmann equation when the correction term is absent. The system ( 108 ) and ( 109 ) have two fixed points with coordinates P(E±) = (±0, 0) which correspond to the limit ρ → 0. The bounce event occurs at ρE = ρ0 at which we have H˙E = −3(2M − 1)/6M − 1 > 0. The stability of the system can be analyzed similar to the way used in the previous sections. A simple study shows the stability properties are analogous to the other models; the evolution of the Universe begins from an unstable state, then k2 2 a2 H Fig. 9 The effective phantom scalar field and the corresponding potential for the model E. As before, the horizontal axis shows the time values. In the model E, the NEC is violated only near the bounce (within the period −1/√Q < t < 1/√Q) in contrast to the other discussed models. A phantom scalar field representation can be used to equivalently describe the model which can be valid when the NEC is violated is violated, we have WE < −1). In this cases, the kinetic energy of scalar field is negative. Thus, for a phantom scalar field in model E we obtain 2H˙E = 3 2M − 1 2ρ 6M − 1 ρ 1 − ρ0 H ρ˙E = − M ρ , , 3HE2 = −9M 2M − 1 ρ 6M − 1 ρ 1 − ρ0 . 8 M φ(eff)E = − 3 Z arcsin #Qt −√2 arctan 2Qt 1 − Qt 2 , V(eff)E = 2MQ (6M − 1)Qt 2 + 1 Z Qt 2 + 1 2 . In Fig. 9 we have presented typical diagrams for the effective phantom field and its corresponding potential. To see the stability properties of model E, we rewrite the field equations as follows + 6H (H 1 + ˙ ) = δ tt + 2 F δT, as the ADM energy constraint (Gtt component of the field equation, if we rewrite ( 8 ) as Gνμ = νμ), 1 + ˙ = − 3 δ ti d xi , 3 Paper [126] has considered models in wh√ich the conservation of EMT is respected. These modes accept h(T) = T [137]. as the ADM momentum constraint (Gti component). Moreover, we have and the perturbed Ricci scalar for the FRLW metric can be obtained as = 0, as the ADM propagation equation (G ij − 1/3δ ij Gll component), as the perturbed version of Raychaudhuri equation (Gii − Gtt component), δR = − (δ + 2F δT) , as the trace equation (Gμμ = μμ), δ˙ + η H δ + ξ k2 a2 v − 3 ˙ = 0, as the time component of perturbed EMT conservation and finally v˙ + 3H σ v − + λδ = 0, as the spatial component of perturbed EMT conservation.4 Note that Eqs. ( 116 )–( 122 ) are the most general equations describing scalar perturbations in minimal f (R, T) gravity for condition F = 1 when a barotropic perfect fluid is included. These equations are not independent, so that it is possible to obtain one equation from another one; for example ( 120 ) follows from multiplying ( 116 ) by 2 then adding it to ( 119 ) and using ( 127 ). In the above equations and relations we have used the following definitions for the source terms, which appear in the right hand sides of field equation ( 8 ), its trace and in Eq. ( 14 ) when is written as ∇β Tβα = α, respectively ν μ = κ2 + F Tνμ − wρF gνμ, = κ2 + F 1 α = κ2 + F T − 4wρF , 1 w∇α (ρF ) − 2 F ∇ αT − ∇β F Tβα . δ¯R = −2 3 ¨ + 12H ˙ + 3H ˙ + 6H˙ λ = w +1 1 21 (1 − w) κ2 F+ F − w , Nη = −(1 + w)F F + (1 + w)(3w − 1)ρF 2 − κ2(w + 3) 2 κ2 + 21 (3 − w)F + (1 + w)(3w − 1)ρF F − (1 + w)(3w − 1)(κ2 + F )ρF , , 2 where Dη = and ( 118 ) ( 119 ) ( 120 ) ( 121 ) ( 122 ) ( 123 ) ( 124 ) ( 125 ) ( 126 ) ( 127 ) ( 128 ) −1 , ( 129 ) ( 130 ) ( 131 ) = 0. ( 133 ) ( 134 ) The gauge-invariant density contrast in the longitudinal gauge is defined as From Eqs. ( 120 ) and ( 127 ) we obtain a dynamical equation for the perturbed potential , as δ¯ρ δ = ρ + 3H v 4 For more details and also the utilized terminology, see [162]. 1 Nσ = 2 (1 − w)F − w(κ2 + F ) + (3w2 + 2w − 2)ρF Dσ = κ2 + 21 (3 − w)F + (1 + w)(3w − 1)ρF , ( 132 ) and prime denotes derivative with respect to the trace. Note that wherever is needed we have used ρ˙ which can be obtained from the EMT conservation equation ( 14 ). Now, using equations ( 121 ) and ( 122 ) we arrive at the evolutionary equation for the matter perturbation, as follows, δ¨ + D1δ˙ + D2δ + ξ k2 −3 ¨ − 3(3σ + 2)H ˙ + a2 k2 2 3 ¨ + 15H ˙ + (6H˙ + 12H 2 + a2 ) + θ δ = 0, Fig. 10 Evolution of δ and for models A, B and C. The behavior of these quantities are similar for the first three model. The black curves are plotted for initial values δ(0) = 0, δ (0) = 1, (0) = 0 and (0) = 1 and the gray ones are drawn for δ(0) = 1, δ (0) = 0, (0) = 1 and (0) = 0 where we have used solution ( 118 ) and defined the following coefficients ξ D1 = (2 + 3σ )H + η − ξ˙ , D2 = ξ˙ k2 η˙ − ξ η + 3σ η H − λξ a2 + 2η H , ( 135 ) ( 136 ) θ = κ2 + (3 − 5w)F − (1 + w)(3w − 1)F ρ ρ . ( 137 ) Hence, we have two differential equations ( 133 ) and ( 134 ) along with relation ( 116 ) to be solved for δ and . Obviously, the coefficients of these equations are some complicated functions of model properties, H, a, ρ , T, F , F , F , and thus it may not be possible to obtain exact solutions. However, one can resort to numerical methods or even in the case of stiff equations, obtain approximated solutions. We have plotted the evolution of δ and in Fig. 10 for models A, B and C. Numerical simulations show that the behavior of perturbations are typically similar in A, B and C models for the same initial values. We have sketched two set of plots in Fig. 10. These models show a zero value for the fluctuations when the initial values δ(0) = 0, δ (0) = 1, (0) = 0 and (0) = 1 are assumed. In this case, the fluctuations tend to zero (left panel) and constant values (right panel) in the regime of large times, see the black curves in Fig. 10. As another case, the fluctuations increase from zero to a maximum finite value in the period of bounce if the initial data are set as δ(0) = 1, δ (0) = 0, (0) = 1 and (0) = 0, see the gray curves. As can be seen from the evolution of scalar perturbations across the bounce, though small temporary fluctuations, no instability occurs during the bounce, nor does it happen in the limit of large times for these models. Unfortunately, the system of differential equations ( 133 ) and ( 134 ) become stiff for models D and E so as it is not possible to plot reasonable diagrams for δ and . In this case, we proceed to obtain approximate solutions. For models E and D (in case in which ϒ = 0), Eqs. ( 133 ) and ( 134 ) at times near the bounce will take the following forms = 0, k2 2k2 a02 + 12H 0 = 0, for which the solutions up to second order can be found as δ = δi + δi t + − 3H 0 + a02 = i + i t + ( 138 ) ( 139 ) ( 140 ) ( 141 ) where the superscript “0” denotes the values of quantities in the limit t → 0 and the subscript i shows the initial values required for integrations. As we see, the solutions are stable in the period of bounce. Note that other coefficients except those which are shown in Eqs. ( 138 ) and ( 139 ) vanish in the limit t → 0. In the limit of large times, we have numerically plotted the evolution of the matter contrast δ and the potential in Fig. 11. As can be seen, from ( 140 ) and ( 141 ), it is obvious that there happens no instability at the period of bounce in models D and E, however, far away from the bounce point, both δ and increase dramatically in model D, see Fig. 11. For model E, in the limit of large times, t → ∞, we get δ¨ − 3ξ ∞ ¨ = 0, ¨ = 0, ( 142 ) ( 143 ) Fig. 11 The evolution of δ and for model D. The Left panels are plotted for initial values δ(0) = 0, δ (0) = 1, (0) = 0 and (0) = 0. The right panels are drawn for δ(0) = 1, δ (0) = 0, (0) = 1 and (0) = 0. The black curves are the solutions of equations approximated for large times and the brown ones are solutions in the limit of small times approximated around the bounce time with the solutions δ = δi + δi t , = i + i t . ( 144 ) ( 145 ) where Q1 and Q2 are some constants. In [137] the authors have analytically shown that for the case of pressure-less perfect fluid, i.e., for w = 0, one finds aSIIN(t ) = , Therefore, depending on the initial values, the perturbations in model E can grow before and after the bounce. Therefore we conclude that, though the scalar-type perturbations in models D and E behave regularly at the bounce point, these solutions are unstable asymptotically. 6 Singular solutions In this section we seek for possible solutions that exhibit singular behavior, specially the big-bang singularity. Some singular solutions have been already studied in the literature which we suffice to give a short discussion for them. Presuming the EMT conservation (i.e., supposing T = T0a−3(1+w)), Eq. ( 18 ) can be solved to give the following solution which shows a singular behavior. This singular solution is the only one which respects the EMT conservation. Relaxing such a constraining condition, one can obtain other kind of solutions. Note that function ( 146 ) is the only form that respects the EMT conservation, that is, to obtain another solution one should assume some suitable form for h(T). For example, in [137] for f (R, T) = R + ακ 2T the authors have obtained 6(α−1) ρ (a) = ρ0a 2−3α , aSIIIN(t ) = ( 147 ) ( 148 ) ( 149 ) h(T) = 2Q1 w + 1 T 23 − w+11 + Q2, 3w + 1 ( 146 ) for a pressure-less matter. In this case singular solution can be found for some valid range of values of the model constant α. Also, the following solution has been obtained for models of form f (R, T) = R + ακ2T−1/2 ρ(a) = 2−2/3 aSIIIIN(t ) = As a new class of (big-bang) singular models, we examine the models which admit the following solution aSIVIN(t ) = a0 t t0 , where being a real positive number and t0 being the time at which the scale factor gets the present value a0. Any solution has to satisfy two of the three equations ( 16 )–( 18 ), for a presumed scale factor function. In this case, two unknown functions T(a) and h(T) should be obtained via solving the resulted equations. To proceed further we make use of the following anzats h(T) = C1Tμ + C2T, where constants Ci ’s and μ should be fixed according to the following considerations. Substituting ( 152 ) and ( 153 ) in ( 18 ) and solving for T(a) we get T(a) = ⎨⎪⎪⎪⎧ 1 ⎫⎪ 1−μ αC1μ(w − 3)T0 ⎪⎪⎬ a − 2μ6((μw−+11))(−w3+w1+) 1 +T0(1−w) ⎪⎪⎪⎭ ⎪⎪⎪⎩ αC1μ(w−3)T0μ +(w − 1)T0 a0 , where an integration constant has been set so that T(a = a0) = T0. To ensure that three functions ( 152 ), ( 153 ) and ( 154 ) provide an analytic solution, they must satisfy at least one of the equations ( 16 ) and ( 17 ). Substituting these functions in ( 16 ) shows that we have a solution only for the case of stiff fluid, i.e., for w = 1. Therefore, a singular solution can be obtained provided that 1 2 (3 − 2)T03 −2 t02ακ2 C1 = − 1 C2 = − α , 1 μ = 2 − 3 , w = 1. 1 3 , T0 1 2−3 ( 156 ) . ( 157 ) Therefore, besides the non-singular solutions obtained in the present paper, one can still find a set of singular solutions. In this brief section in addition to addressing some previous results we obtained a new singular model, though a coherent study can be performed in order to deal with possible conditions for which big-bang singularity would occur; however working on this issue is beyond the scope of the present paper and comprehensive studies on this subject will be reported elsewhere. It is worthwhile to mention that some studies have been already made to consider other forms of singular solutions, e.g., [163]. Beside the above results, Bianchi type I cosmological model with magnetized strange quark matter in the framework of f (R, T) gravity have been investigated and it is found that the model begins with big-bang and ends with big rip [164]. Using Lie point symmetry analysis method, the authors of [165] have shown that for a Bianchi type I spacetime, both singular (big-bang) and nonsingular solutions could exist subject to the type of specified symmetry. Recently, the authors of [166], have considered some cosmological features of f (T ) gravity (where here T denotes torsion scalar) using the dynamical system approach both generally and for some specific forms of f (T ) functions. The core of their studies is taking the advantage of this fact that the torsion scalar can be used interchangeably with the Hubble parameter (i.e., T = −6H 2). Thus, the field equations reduce to a single equation (in the case of pressure-less matter) in the form of H˙ = F (H ), since the matter density can also be rewritten as a function of the Hubble parameter. Briefly, they have shown that in f (T ) gravity a single equation (which can be interpreted as a simple one dimensional dynamical system) can govern the dynamics of field equations. Benefiting this useful result they investigated phase space portraits of various cosmological evolutions such as, singular and non-singular solutions. Likewise, one may be motivated to utilize such an approach in order to investigate the cosmological solutions of f (R, T) gravity (especially, in the case of present work, i.e., the function given in ( 15 )) through phase portrait diagrams. However, looking at Eqs. ( 17 ) and ( 18 ) one finds that it is impossible to obtain an equation like H˙ = F (H ) so that it reflects full information of the field equations. In this case for an assumed function h(T), we have a two dimensional dynamical system without any further reduction. Thus, the procedure proposed in [166] would be generally failed in f (R, T) gravity. B C D E Models a(t) *The subscripts A,..,E and “eff” are dropped for abbreviation 7 Concluding remarks In the present work we studied classical bouncing behavior of the Universe in the framework of f (R, T) = R + h(T) gravity theories. We assumed a single perfect fluid in a spatially flat, homogeneous and isotropic FLRW background. Having obtained the resulted field equations, we employed the concept of effective fluid (which is firstly introduced in [145]) via defining an effective energy density and pressure and also reformulating the field equations in terms of these fluid components. In this picture, one could recast the field equations of f (R, T) gravity for a real perfect fluid into GR field equations for an effective fluid. It is also shown that in a modified gravity model the energy conditions are usually obtained by using the effective EMT, not the one for real fluids. In f (R, T) gravity, the definitions for effective energy density and pressure have already been used to obtain the energy conditions [154]. The effective fluid has an EoS of the form, p(eff) = Y(ρ(eff)), which corresponds to an h(T) function. In this method one firstly specifies an effective EoS or a condition on the effective components and then obtains the corresponding h(T) function and other cosmological quantities. It is also possible to make a link between f (R, T) gravity in effective picture and models which use some exotic or dark component with unusual EoS. These models which have been widely discussed in the literature (to deal with some cosmological issues) are also called theories with “generalized EoS”. The mathematical representation of effective components provides a setting within which unusual interactions of a real perfect fluid with gravitational field can be translated as the presence of an exotic fluid which admits the EoS of the form p(eff) = Y(ρ(eff)). In this paper we have shown that it is possible to recover generalized EoS models which have been previously studied in the literature ( see e.g., [52,152,153,157,158]), in the framework of f (R, T) P = cons. − 3 + 2S√3R2(S2−1)ρ+9 ρ 2 3R2(S2−1) + R2(S2−1) n1 ρ − α 2C (1 + n) ρ − √αm ρ + α2 + α 32M − 1 ρ ± Z 2MQ 2Q M2 − 3ZQ ρ − Z h(T) 2 α P + 1−w3w T 2T2 − 47α T + 2 T n+21 T n+1 − α + 2 √T − (2m−α1)T + 2 (w+1) T wn−+21 + 23 2n+32w(n−−11) − α(2n+w−3) T + ρ ρ0a−3 ρ0a−1 ρ0a−6/n ζ−α a3 2 8a6m2 ρ0a−3(w+1)/n w w > 1/3 −1/5 1 1 gravity. Therefore, the problem of exotic fluid in the context of generalized EoS models which are mostly without a determined Lagrangian may be discussed in a Lagrangian based theory of gravity like f (R, T) gravity. In the current research, we discussed four different bouncing models in f (R, T) gravity. We labeled them as the models A, B, C, D and E and briefly mentioned their main properties in Table 1. Each model can be specified either by an h(T) function or by an effective EoS. Models A-D mimic an asymptotic de Sitter expansion in the far past and future of the bounce. The model A corresponds to a constant effective pressure, p(eff) = P; for the model B we have p(eff)B = −ρ(eff)B/3 + #bBρ(eff)B + dB + eB, the model C is specified by p(eff)C = jCρ(eff)C + eC, the model D corresponds to p(eff)D = ρ(eff)D + #bDρ(eff)D + dD + eD and finally the model E obeys the EoS, p(eff)E = aEρ(eff)E + #bEρ(eff)E + dE + eE, where the constants b, d, e and j are written in terms of model parameters. In all models the matter density grows to a maximum value at the bounce which corresponds to a minimum for the scale factor. The effective density varies from zero at the bounce to a positive value in the far past and future of the bounce. The effective pressure varies between negative values; in model A, it is a constant, in the model B it increases at the bounce, in the models C and E it decreases and the model D admits both behaviors. The effective EoS has the property −∞ < W < −1 when the bounce point is approached. The Hubble parameter satisfies H (t ) = 0 and d H/dt > 0 at the event of bounce and also all its time derivatives have regular behavior for all models. Therefore, these bouncing solutions do not exhibit future singularities which are classified in the literature of cosmological solutions. We can consider the inherent exoticism hidden behind f (R, T) gravity in another way. As already we mentioned, this issue can be described as an unusual interaction between gravitational field and normal matter or introducing an effective fluid. From the point of view of the energy conditions, in all discussed models the SEC and NEC are violated (note that for a normal fluid NEC is not violated [47]) near the bounce and the effective density gets minimized to zero. Such a result has been previously predicted in GR [155]. As discussed in [47], the exoticness can be understood as a minimization in the effective pressure. In the other words, a minimum in the effective energy density corresponds to a minimum in the scale factor. Such a behavior is permitted provided that W < −1. Note that, for a normal matter, a minimum (maximum) compression leads to a minimum (maximum) energy density. Thus, in f (R, T) gravity an abnormal or effective fluid which leads to an uncommon balance in the density and pressure can be responsible for the bouncing behavior. An interesting feature of the bouncing solution in f (R, T) gravity is that one can construct solutions in which the SEC is respected by the real perfect fluid. Such solutions cannot be found in GR [155]. Also note that the real perfect fluid with w > −1 never violates the NEC. Therefore, we have solutions without the future singularities and all energy conditions can be respected by a real perfect fluid. By this discussion, one may use the definition of an (effective) phantom scalar field if one asks for the matter source to be reinterpreted as that of a scalar matter field. We obtained the equivalent scalar field φ(eff)(t ) and its corresponding potential V(eff)(t ) in each case. Moreover, we have studied the dynamical system representation of these models. We found that the evolution of the Universe can be displayed by trajectories which initially start from an unstable state, passing through an unstable fixed point (the bounce event) and finally are absorbed by a stable point. The initial and final states are de-Sitter era in models A, B, C and D and the decelerated expanding Universe in model E. Another important issue discussed in this work is related to the study of stability of bouncing solutions through scalar-type cosmological matter perturbations in the bouncing universe. Our numerical analysis of density perturbations for models A, B and C revealed that, though a slight jump (depending on the initial conditions) at the bounce point, the amplitude of matter density perturbation (δ) and perturbed potential ( ) behave regularly throughout the bounce phase. Therefore, since the time interval during which the fluctuations that occur within density contrast and perturbed field is short, the instabilities do not have enough time to grow to a significant magnitude. However, this case does not happen for the two remaining models. As the final remarks we should emphasize that our models were obtained by indicating different conditions on the effective density and pressure which led to different h(T) functions. This means that the models A, B, C, D and E are not the only possible models for the bouncing behavior. It is obvious that one can still choose other h(T) functions or consider other assumptions on the effective density and pressure to obtain new bouncing solutions (with even new features). Our aim was to show the existence of varieties of bouncing solutions in f (R, T) gravity and study their properties. Especially, our study was confined to the Lagrangians of type f (R, T) = R + h(T) though other forms of Lagrangians can be investigated. The other issue is that our study was performed in the effective picture. In case such an approach is not taken seriously, one can think of it as only an alternative mathematical method. One can still investigate a nonsingular cosmological scenario without employing the equations which are written in terms of the effective quantities. In this case it is enough to assume a Lagrangian and solve the field equations to inspect for a bouncing solution. However, cosmological solutions for the f (R, T) gravity model presented here are not singularity free and as we observed under certain conditions, a class singular solutions could be obtained. 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 SCOAP3. 1. J.P. Ostriker , P.J. Steinhardt , Cosmic concordance . arXiv:astro-ph/9505066 2. P. Coles , F. Lucchin , Cosmolgy: The Origin and Evolution of Cosmic Structure (Wiley, Chichester, 2002 ) 3. R.A. Liddle , D.H. Lyth , The Primordial Density Perturbation: Cosmology, Inflation and the Origin of Structure (Cambrige University Press, Cambrige, 2009 ) 4. E.W. Kolb , M.S. Turner , The Early Universe. Frontiers in Physics (Avalon Publishing, Redwood City , 1994 ) 5. M.Y. Khlopov , A.D. Linde , Is it easy to save the gravitino? Phys. Lett. B 138 , 265 ( 1984 ) 6. J. Ellis , J.E. Kim , D.V. Nanopoulos , Cosmological gravitino regeneration and decay . Phys. Lett. B 145 , 181 ( 1984 ) 7. M. Kawasaki , T. Moroi, Gravitino production in the inflationary universe and the effects on big bang nucleosynthesis . Prog. Theor. Phys . 93 , 879 ( 1995 ) 8. MYu. Khlopov, YuL Levitan, E.V. Sedelnikov , I.M. Sobol , Nonequilibrium Cosmological Nucleosynthesis of Light Elements: Calculations by the Monte Carlo Method . Phys. At. Nuclei 57 , 1393 ( 1994 ) 9. M. Kawasaki , K. Kohri , T. Moroi, Hadronic decay of latedecaying particles and Big-Bang nucleosynthesis . Phys. Lett. B 625 , 7 ( 2005 ) 10. G.D. Coughlan , W. Fischler , E.W. Kolb , S. Raby , G.G. Ross , Cosmological problems for the Polonyi potential . Phys. Lett. B 131 , 59 ( 1983 ) 11. J.R. Ellis , D.W. Nanopoulos , M. Quiros , On the axion, dilaton, Polonyi, gravitino and shadow matter problems in supergravity and superstring models . Phys. Lett. B 174 , 176 ( 1986 ) 12. T. Banks , D.B. Kaplan , A.E. Nelson , Cosmological implications of dynamical supersymmetry breaking . Phys. Rev. D 49 , 779 ( 1994 ) 13. B. de Carlos, J.A. Casas , F. Quevedo , E. Roulet, Modelindependent properties and cosmological implications of the dilaton and moduli sectors of 4-D strings . Phys. Lett. B 318 , 447 ( 1993 ) 14. A.D. Dolgov , M.V. Sazhin , I.A.B. Zeldovich , Basics of modern cosmology (Editions Frontieres , France, 1990 ) 15. Z. Xing , S. Zhou , Neutrinos in Particle Physics, Astronomy and Cosmology (Springer, Berlin, 2011 ) 16. T. Padmanabhan, Cosmology and Astrophysics Through Problems (Cambridge University Press, Cambridge, 1996 ) 17. S. Dodelson, Modern Cosmology (Academic Press, Cambridge, 2003 ) 18. G.F.R. Ellis , R. Maartens , M.A.H. MacCallum , Relativistic Cosmology (Cambridge University Press, Cambridge, 2012 ) 19. A. Liddle , An Introduction to Modern Cosmology (Wiley, New York, 2015 ) 20. D.W. Sciama, Modern Cosmology and the Dark Matter Problem (Cambridge University Press, Cambridge, 1993 ) 21. C. Bambi , A.D. Dolgov , Introduction to Particle Cosmology: The Standard Model of Cosmology and its Open Problems (Springer, Berlin, 2015 ) 22. K. Freese , Status of dark matter in the universe . Int. J. Mod. Phys. D 26 , 1730012 ( 2017 ) 23. S. Weinberg, The cosmological constant problem . Rev. Mod. Phys. 61 , 1 ( 1989 ) 24. L. Amendola , S. Tsujikawa , Dark Energy: Theory and Observations (Cambridge University Press, Cambridge, 2010 ) 25. J. Gleyzes , Dark Energy and the Formation of the Large Scale Structure of the Universe (Springer , Berlin, 2016 ) 26. A.H. Guth , Inflationary universe: a possible solution to the horizon and flatness problems . Phys. Rev. D 23 , 347 ( 1981 ) 27. K. Sato , First-order phase transition of a vacuum and the expansion of the universe . Mon. Not. R. Astron. Soc . 195 , 467 ( 1981 ) 28. A.D. Linde , A new inflationary universe scenario: a possible solution of the horizon, flatness, homogeneity, isotropy and primordial monopole problems . Phys. Lett. B 108 , 389 ( 1982 ) 29. S.W. Hawking , R. Penrose , The singularities of gravitational collapse and cosmology . Proc. R. Soc. Lond. A 314 , 529 ( 1970 ) 30. R. Penrose , Gravitational collapse and space-time singularities . Phys. Rev. Lett . 14 , 57 ( 1965 ) 31. S.W. Hawking , Occurrence of singularities in open universes . Phys. Rev. Lett . 15 , 689 ( 1965 ) 32. S.W. Hawking , The occurrence of singularities in cosmology . Proc. R. Soc. Lond. A 294 , 511 ( 1966 ) 33. S.W. Hawking , The occurrence of singularities in cosmology . II. Proc. R. Soc. Lond. A 295 , 490 ( 1966 ) 34. R.P. Geroch , Singularities in closed universes . Phys. Rev. Lett . 17 , 445 ( 1966 ) 35. S.W. Hawking , The occurrence of singularities in cosmology. III Causality and singularities . Proc. R. Soc. Lond. A 300 , 187 ( 1967 ) 36. S.W. Hawking , G.F.R. Ellis , The Large Scale Structure of SpaceTime (Cambridge University Press, Cambridge, 1975 ) 37. J.M.M. Senovilla , Singularity theorems and their consequences . Gen. Relativ. Gravit . 29 , 701 ( 1997 ) 38. F.J. Tipler , General relativity and conjugate ordinary differential equations . J. Differ. Equ . 30 , 65 ( 1978 ) 39. F.J. Tipler , Energy conditions and spacetime singularities . Phys. Rev. D 17 , 2521 ( 1978 ) 40. A. Borde , Geodesic focusing, energy conditions and singularities . Class. Quantum Gravity 4 , 343 ( 1987 ) 41. A. Vilenkin , Did the universe have a beginning? Phys. Rev. D 46 , 2355 ( 1992 ) 42. A. Borde , A. Vilenkin , Eternal inflation and the initial singularity . Phys. Rev. Lett . 72 , 3305 ( 1994 ) 43. A. Borde , Open and closed universes, initial singularities and inflation . Phys. Rev. D 50 , 3692 ( 1994 ) 44. A. Borde , A. Vilenkin , Singularities in inflationary cosmology: a review . Int. J. Mod. Phys. D 5 , 813 ( 1996 ) 45. A. Borde , A.H. Guth , A. Vilenkin , Inflationary spacetimes are incomplete in past directions . Phys. Rev. Lett . 90 , 151301 ( 2003 ) 46. G. Calcagni, Classical and Quantum Cosmology, Graduate Texts in Physics (Springer, Berlin, 2017 ) 47. M. Novello , S.E. Perez Bergliaffa , Bouncing cosmologies . Phys. Rep . 463 , 127 ( 2008 ) 48. V.F. Mukhanov , R.H. Brandenberger , A nonsingular universe . Phys. Rev. Lett. 68 , 1969 ( 1992 ) 49. R.H. Brandenberger , V.F. Mukhanov , A. Sornborger , Cosmological theory without singularities . Phys. Rev. D 48 , 1629 ( 1993 ) 50. Y. Choquet-Bruhat , General Relativity and the Einstein Equations (OUP , Oxford, 2008 ) 51. M. Bojowald , Essay: initial conditions for a universe . Gen. Relativ. Gravit . 35 , 1877 ( 2003 ) 52. F. Contreras , N. Cruz , G. Palma, Bouncing solutions from generalized EoS . arXiv: 1701 .03438 [gr-qc] 53. R.H. Brandenberger , Introduction to early universe cosmology . arXiv:1103.2271 [astro-ph.CO] 54. D. Battefeld , P. Peter , A critical review of classical bouncing cosmologies . Phys. Rep . 571 , 1 ( 2015 ) 55. Y.-K.E. Cheung , X. Song , S. Li , Y. Li , Y. Zhu , The CST bounce universe model-a parametric study . arXiv:1601 .03807 [gr-qc] 56. J.E. Lidsey , D. Wands , E.J. Copeland , Superstring cosmology . Phys. Rep . 337 , 343 ( 2000 ) 57. B.K. Xue , D. Garfinkle , F. Pretorius , P.J. Steinhardt , Nonperturbative analysis of the evolution of cosmological perturbations through a nonsingular bounce . Phys. Rev. D 88 , 083509 ( 2013 ) 58. G. Date, G.M. Hossain , Genericness of a big bounce in isotropic loop quantum cosmology . Phys. Rev. Lett . 94 , 011302 ( 2005 ) 59. J. Mielczarek , T. Stachowiak , M. Szydlowski , Exact solutions for a big bounce in loop quantum cosmology . Phys. Rev. D 77 , 123506 ( 2008 ) 60. P. Singh , Are loop quantum cosmos never singular? Class . Quantum Gravity 26 , 125005 ( 2009 ) 61. A. Ashtekar , Singularity resolution in loop quantum cosmology: a brief overview . J. Phys. Conf. Ser . 189 , 012003 ( 2009 ) 62. Y.-F. Cai , E. Wilson-Ewing, Non-singular bounce scenarios in loop quantum cosmology and the effective field description . J. Cosmol. Astropart. Phys . 03 , 026 ( 2014 ) 63. M. Bojowald , Loop quantum cosmology . Living Rev. Relativ . 8 , 11 ( 2005 ) 64. K. Banerjee , G. Calcagni, M. Martin-Benito , Introduction to loop quantum cosmology . SIGMA 8 , 016 ( 2012 ) 65. A. Ashtekar , P. Singh , Loop quantum cosmology: a status report . Class. Quantum Gravity 28 , 213001 ( 2011 ) 66. A. Ashtekar , T. Pawlowski , P. Singh , Quantum nature of the big bang: improved dynamics . Phys. Rev. D 74 , 084003 ( 2006 ) 67. P. Peter, E.J.C. Pinho , N. Pinto-Neto , A non inflationary model with scale invariant cosmological perturbations . Phys. Rev. D 75 , 023516 ( 2007 ) 68. N. Pinto-Neto , J.C. Fabris , Quantum cosmology from the de Broglie-Bohm perspective . Class. Quantum Gravity 30 , 143001 ( 2013 ) 69. P. Singh , A. Toporensky , Big crunch avoidance in k = 1 semiclassical loop quantum cosmology . Phys. Rev. D 69 , 104008 ( 2004 ) 70. J.E. Lidsey , D.J. Mulryne , N.J. Nunes , R. Tavakol , Oscillatory universes in loop quantum cosmology and initial conditions for inflation . Phys. Rev. D 70 , 063521 ( 2004 ) 71. P. Singh , K. Vandersloot , G.V. Vereshchagin , Non-singular bouncing universes in loop quantum cosmology . Phys. Rev. D 74 , 043510 ( 2006 ) 72. M. Bojowald , Quantum gravity in the very early universe . Nucl. Phys. A 862-863 , 98 ( 2011 ) 73. A. Ashtekar , Loop quantum gravity and the Planck regime of cosmology . Fundam. Theor. Phys . 177 , 323 ( 2014 ) 74. G. Calcagni, L. Modesto , P. Nicolini , Super-accelerating bouncing cosmology in asymptotically-free non-local gravity . Eur. Phys. J. C 74 , 2999 ( 2014 ) 75. T. Biswas, T. Koivisto , A. Mazumdar , Towards a resolution of the cosmological singularity in non-local higher derivative theories of gravity . JCAP 1011 , 008 ( 2010 ) 76. J.F. Donoghue , B.K. El-Menoufi , Non-local quantum effects in cosmology 1: quantum memory, non-local FLRW equations and singularity avoidance . Phys. Rev. D 89 , 104062 ( 2014 ) 77. Y. -D. Li , L. Modesto , L. Rachwal, Exact solutions and spacetime singularities in nonlocal gravity . JHEP 2015 , 1 ( 2015 ) 78. L. Amendola , N. Burzilla , H. Nersisyan , Quantum gravity inspired nonlocal gravity model . Phys. Rev. D 96 , 084031 ( 2017 ) 79. L. Modesto , L. Rachwal, Nonlocal quantum gravity: a review . Int. J. Mod. Phys. D 26 , 1730020 ( 2017 ) 80. L. Modesto, Super-renormalizable quantum gravity . Phys. Rev. D 86 , 044005 ( 2012 ) 81. Y.-F. Cai , D.A. Easson , R. Brandenberger , Towards a non-singular bouncing cosmology . J. Cosmol. Astropart. Phys . 1208 , 020 ( 2012 ) 82. R. Brandenberger , The matter bounce alternative to inflationary cosmology . arXiv:1206.4196 [astro-ph.co] 83. R. Brandenberger , Alternatives to the inflationary paradigm of structure formation . Int. J. Mod. Phys. Conf. Ser . 01 , 67 ( 2011 ) 84. R. Brandenberger , Cosmology of the very early universe . AIP Conf. Proc. 1268 , 3 ( 2010 ) 85. E. Wilson-Ewing, The matter bounce scenario in loop quantum cosmology . arXiv:1211 .6269 [gr-qc] 86. Y.F. Cai , T. Qiu , Y.S. Piao , M. Li , X. Zhang, Bouncing universe with quintom matter . JHEP 0710 , 071 ( 2007 ) 87. Y.F. Cai , T. Qiu , R. Brandenberger , Y.S. Piao , X. Zhang, On perturbations of quintom bounce . J. Cosmol. Astropart. Phys . 0803 , 013 ( 2008 ) 88. Y.F. Cai , X. Zhang, Evolution of metric perturbations in quintom bounce model . J. Cosmol. Astropart. Phys . 0906 , 003 ( 2009 ) 89. Y.-F. Cai , T. Qiu , R. Brandenberger , X. Zhang, A nonsingular cosmology with a scale invariant spectrum of cosmological perturbations from Lee-Wick theory . Phys. Rev. D 80 , 023511 ( 2009 ) 90. C. Lin , R. Brandenberger , L.P. Levasseur , A matter bounce by means of ghost condensation . J. Cosmol. Astropart. Phys . 1104 , 019 ( 2011 ) 91. T. Qiu , J. Evslin , Y.F. Cai , M. Li , X. Zhang , Bouncing Galileon cosmologies. J. Cosmol. Astropart. Phys . 1110 , 036 ( 2011 ) 92. D.A. Easson , I. Sawicki , A. Vikman , G-bounce. J. Cosmol. Astropart. Phys . 1111 , 021 ( 2011 ) 93. M.G. Brown , K. Freese, W.H. Kinney , The phantom bounce: a new oscillating cosmology . J. Cosmol. Astropart. Phys . 0803 , 002 ( 2008 ) 94. V. Dzhunushaliev , V. Folomeev , K. Myrzakulov , R. Myrzakulov , Phantom fields: bounce solutions in the early Universe and Sbranes . Int. J. Mod. Phys. D 17 , 2351 ( 2008 ) 95. K. Nozari , S.D. Sadatian , Bouncing universe with a nonminimally coupled scalar field on a moving domain wall . Phys. Lett. B 676 , 1 ( 2009 ) 96. E.N. Saridakis , S.V. Sushkov , Quintessence and phantom cosmology with non-minimal derivative coupling . Phys. Rev. D 81 , 083510 ( 2010 ) 97. A. Banijamalia , B. Fazlpour , Phantom behavior bounce with tachyon and non-minimal derivative coupling . J. Cosmol. Astropart. Phys . 01 , 039 ( 2012 ) 98. C. Barragán , G.J. Olmo , H. Sanchis-Alepuz , Bouncing cosmologies in Palatini F (R) gravity . Phys. Rev. D 80 , 024016 ( 2009 ) 99. V.K. Oikonomou , Loop quantum cosmology matter bounce reconstruction from F (R) gravity using an auxiliary field . Gen. Relativ. Gravit . 47 , 126 ( 2015 ) 100. S.D. Odintsov , V.K. Oikonomou , Matter bounce loop quantum cosmology from F (R) gravity . Phys. Rev. D 90 , 124083 ( 2014 ) 101. S.D. Odintsov , V.K. Oikonomou , Bouncing cosmology with future singularity from modified gravity . Phys. Rev. D 92 , 024016 ( 2015 ) 102. V.K. Oikonomou , Superbounce and loop quantum cosmology Ekpyrosis from modified gravity . Astrophys. Space Sci . 359 , 30 ( 2015 ) 103. Y.-F. Cai , S.-H. Chen , J.B. Dent , S. Dutta , E.N. Saridakis , Matter bounce cosmology with the f (T ) gravity . Class. Quantum Gravity 28 , 215011 ( 2011 ) 104. K. Bamba , G.G.L. Nashed , W. El Hanafy, ShK Ibraheem, Bounce inflation in f (T ) cosmology: a unified inflaton-quintessence field . Phys. Rev. D 94 , 083513 ( 2016 ) 105. A. Kehagias , E. Kiritsis, Mirage cosmology . J. High Energy Phys . 11 , 022 ( 1999 ) 106. M. Gasperini , Repulsive gravity in the very early universe . Gen. Relativ. Gravit . 30 , 12 ( 1998 ) 107. T. Stachowiak, M. Szydlowski , Exact solutions in bouncing cosmology . Phys. Lett. B 646 , 209 ( 2007 ) 108. S.D. Brechet , M.P. Hobson , A.N. Lasenby , Classical big-bounce cosmology: dynamical analysis of a homogeneous and irrotational Weyssenhoff fluid . Class. Quantum Gravity 25 , 245016 ( 2008 ) 109. N.J. Poplawski , Cosmology with torsion: an alternative to cosmic inflation . Phys. Lett. B 694 , 181 ( 2010 ) 110. N.J. Poplawski , Big bounce from spin and torsion . Gen. Relativ. Gravit . 44 , 1007 ( 2012 ) 111. N.J. Poplawski , Nonsingular, big -bounce cosmology from spinortorsion coupling . Phys. Rev. D 85 , 107502 ( 2012 ) 112. J. Magueijo , T.G. Zlosnik , T.W.B. Kibble , Cosmology with a spin . Phys. Rev. D 87 , 063504 ( 2013 ) 113. H. Hadi , Y. Heydarzade , M. Hashemi , F. Darabi , Emergent cosmos in Einstein-Cartan theory . Eur. Phys. J. C 78 , 38 ( 2018 ) 114. R. Brandenberger , Matter bounce in Horava-Lifshitz cosmology . Phys. Rev. D 80 , 043516 ( 2009 ) 115. T. Biswas , A.S. Koshelev , A. Mazumdar , S.Y. Vernov , Stable bounce and inflation in non-local higher derivative cosmology . J. Cosmol. Astropart. Phys . 08 , 024 ( 2012 ) 116. B. Dragovich , On nonlocal modified gravity and cosmology , in Lie Theory and Its Applications in Physics. Springer Proceedings in Mathematics and Statistics , ed. by V. Dobrev (Springer, Tokyo, 2014 ) 117. S.D. Odintsov , V.K. Oikonomou , E.N. Saridakis , Superbounce and loop quantum ekpyrotic cosmologies from modified gravity: F (R) , F (G) and F (T ) . Ann. Phys . 363 , 141 ( 2015 ) 118. J. Khoury , B.A. Ovrut , P.J. Steinhardt , N. Turok , The ekpyrotic universe: colliding branes and the origin of the hot big-bang . Phys. Rev. D 64 , 123522 ( 2001 ) 119. J.-L. Lehners , Ekpyrotic and cyclic cosmology . Phys. Rep . 465 , 223 ( 2008 ) 120. M. Gasperini , G. Veneziano, Pre-big bang in string cosmology . Astropart. Phys. 1 , 317 ( 1993 ) 121. T. Biswas , A. Mazumdar , W. Siegel, Bouncing universes in stringinspired gravity . J. Cosmol. Astropart. Phys . 0603 , 009 ( 2006 ) 122. C. Kounnas , H. Partouche , N. Toumbas , Thermal duality and nonsingular cosmology in d-dimensional superstrings . Nucl. Phys. B 855 , 280 ( 2012 ) 123. I. Florakis, C. Kounnas , H. Partouche , N. Toumbas , Non-singular string cosmology in a 2d hybrid model . Nucl. Phys. B 844 , 89 ( 2011 ) 124. R.H. Brandenberger , C. Kounnas , H. Partouche , S.P. Patil , N. Toumbas , Cosmological perturbations across an S-brane . J. Cosmol. Astropart. Phys . 03 , 015 ( 2014 ) 125. T. Harko , F.S.N. Lobo , S. Nojiri , S.D. Odintsov , f (R, T ) gravity . Phys. Rev. D 84 , 024020 ( 2011 ) 126. F.G. Alvarenga , A. de la Cruz-Dombriz, M.J.S. Houndjo , M.E. Rodrigues , D. Sáez-Gómez, Dynamics of scalar perturbations in f (R, T ) gravity . Phys. Rev. D 87 , 103526 ( 2013 ) 127. H. Shabani , M. Farhoudi , f (R, T ) cosmological models in phasespace . Phys. Rev. D 88 , 044048 ( 2013 ) 128. T. Harko, Thermodynamic interpretation of the generalized gravity models with geometry-matter coupling . Phys. Rev. D 90 , 044048 ( 2013 ) 129. H. Shabani , M. Farhoudi , Cosmological and solar system consequences of f (R, T ) gravity models . Phys. Rev. D 90 , 044031 ( 2014 ) 130. C.P. Singh , V. Singh , Reconstruction of modified f (R, T ) gravity with perfect fluid cosmological models . Gen. Relativ. Gravit . 46 , 1696 ( 2014 ) 131. M. Sharif , M. Zubair , Cosmological reconstruction and stability in f (R, T ) gravity . Gen. Relativ. Gravit. 46 , 1723 ( 2014 ) 132. M.E.S. Alves , P.H.R.S. Moraes , J.C.N. de Araujo , M. Malheiro, Gravitational waves in f (R, T ) and f (R, T φ ) theories of gravity . Phys. Rev. D 94 , 024032 ( 2016 ) 133. G. Sun, Y.-C. Huang, The cosmology in f (R, T ) gravity without dark energy . Int. J. Mod. Phys. D 25 , 1650038 ( 2016 ) 134. R. Zaregonbadi , M. Farhoudi , N. Riazi , Dark matter from f (R, T ) gravity . Phys. Rev. D 94 , 084052 ( 2016 ) 135. R. Zaregonbadi , M. Farhoudi , Cosmic acceleration from mattercurvature coupling . Gen. Relativ. Gravit . 48 , 142 ( 2016 ) 136. H. Shabani , A.H. Ziaie , Stability of the Einstein static universe in f (R, T ) gravity . Eur. Phys. J. C 77 , 31 ( 2017 ) 137. H. Shabani , A.H. Ziaie , Consequences of energy conservation violation: Late time solutions of ( T )C D M subclass of f (R, T ) gravity using dynamical system approach . Eur. Phys. J. C 77 , 282 ( 2017 ) 138. H. Shabani , A.H. Ziaie , Late-time cosmological evolution of a general class of f (R, T ) gravity with minimal curvature-matter coupling . Eur. Phys. J. C 77 , 507 ( 2017 ) 139. H. Shabani, Cosmological consequences and statefinder diagnosis of non-interacting generalized Chaplygin gas in f( R, T) gravity . Int. J. Mod. Phys. D 26 , 1750120 ( 2017 ) 140. P.H.R.S. Moraes , W. de Paula, P.A.C. Correa , Charged wormholes in f (R, T) extended theory of gravity . arXiv:1710 .07680 [gr-qc] 141. P.H.R.S. Moraes , P.K. Sahoo , Modeling wormholes in f( R, T) gravity . Phys. Rev. D 96 , 044038 ( 2017 ) 142. D. Deb , F. Rahaman , S. Ray , B.K. Guhaa , Strange stars in f (R, T ) gravity . J. Cosmol. Astropart. Phys . 03 , 044 ( 2018 ) 143. P.H.R.S. Moraes , P.A.C. Correa , Evading the non-continuity equation in the f (R, T ) formalism . Eur. Phys. J. C 78 , 192 ( 2018 ) 144. P.K. Sahoo , P.H.R.S. Moraes , P. Sahoo , Wormholes in R2-gravity within the f (R, T ) formalism . Eur. Phys. J. C 78 , 46 ( 2018 ) 145. H. Shabani , A.H. Ziaie , Interpretation of f (R, T ) gravity in terms of a conserved effective fluid . Int. J. M. Phys. A 33 , 1850050 ( 2018 ) 146. H. Štefancˇic´, Dark energy transition between quintessence and phantom regimes: an equation of state analysis . Phys. Rev. D 71 , 124036 ( 2005 ) 147. S. Nojiri , S.D. Odintsov , Inhomogeneous equation of state of the universe: phantom era, future singularity, and crossing the phantom barrier . Phys. Rev. D 72 , 023003 ( 2005 ) 148. J.D. Barrow , Graduated inflationary universes . Phys. Lett. B 235 , 40 ( 1990 ) 149. S. Mukherjee , B.C.D. Paul , A. Beesham, Emergent universe with exotic matter . Class. Quantum Gravity 23 , 46927 ( 2006 ) 150. H. Štefancˇic´, Expansion around the vacuum equation of state: sudden future singularities and asymptotic behavior . Phys. Rev. D 71 , 084024 ( 2005 ) 151. S. Nojiri , S.D. Odintsov , S. Tsuijikawa , Properties of singularities in the (phantom) dark energy universe . Phys. Rev. D 71 , 063004 ( 2005 ) 152. F. Contreras , N. Cruz , E. Gonzàles, Generalized equations of state and regular universes . J. Phys. Conf. Ser . 720 , 012014 ( 2016 ) 153. E. Babichev , V. Dokuchaev , Yu. Eroshenko, Dark energy cosmology with generalized linear equation of state . Class. Quantum Gravity 22 , 143 ( 2005 ) 154. M. Sharif , M. Zubair , Energy conditions constraints and stability of power law solutions in f (R, T ) gravity . J. Phys. Soc. Jpn . 82 , 014002 ( 2013 ) 155. C. Molina-París , M. Visser , Minimal conditions for the creation of a Friedman-Robertson-Walker universe from a bounce . Phys. Lett. B 455 , 90 ( 1999 ) 156. J.A. Frieman , M.S. Turner , D. Huterer , Dark energy and the accelerating universe . Annu. Rev. Astron. Astrophys . 46 , 385 ( 2008 ) 157. P.-H. Chavanis , A cosmological model based on a quadratic equation of state unifying vacuum energy, radiation, and dark energy . J. Gravity 2013 , 682451 ( 2013 ) 158. K. Bamba , S. Capozziello , S. Nojiri , S.D. Odintsov , Dark energy cosmology: the equivalent description via different theoretical models and cosmography tests . Astrophys. Space Sci . 342 , 155 ( 2012 ) 159. Y.F. Cai , T. Qiu , R. Brandenberger , Y.S. Piao , X. Zhang , A CDM bounce scenario . J. Cosmol. Astropart. Phys . 03 , 006 ( 2015 ) 160. S. Tsujikawa , K. Uddin , R. Tavakol , Density perturbations in f (R, T ) gravity in metric and Palatini formalisms . Phys. Rev. D 77 , 043007 ( 2008 ) 161. K.A. Malik , D. Wands , Adiabatic and entropy perturbations with interacting fluids and fields . J. Cosmol. Astropart. Phys . 02 , 007 ( 2005 ) 162. J.-C. Hwang , H. Noh , Gauge-ready formulation of cosmological kinetic theory in generalized gravity theory . Phys. Rev. D 65 , 023512 ( 2001 ) 163. M.J.S. Houndjo , C.E.M. Batista , J.P. Campos , O.F. Piattel , Finitetime singularities in f (R, T ) gravity and the effect of conformal anomaly . Can. J. Phys . 91 , 547 ( 2013 ) 164. P.K. Sahoo , P. Sahoo , B.K. Bishi , S. Aygn , Magnetized strange quark model with big rip singularity in f (R, T ) gravity . Mod. Phys. Lett. A 32 , 1750105 ( 2018 ) 165. A.K. Yadav , A.T. Ali , Invariant Bianchi type I models in f (R, T ) gravity . Int. J. Geom. Methods Mod. Phys . 15 , 1850026 ( 2018 ) 166. A. Awad , W.E. Hanafy , G.G.L. Nashed , E.N. Saridakis , Phase portraits of general f (T ) cosmology . J. Cosmol. Astropart. Phys . 02 , 052 ( 2018 )


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-5886-x.pdf

Hamid Shabani, Amir Hadi Ziaie. Bouncing cosmological solutions from \(f(\mathsf{R,T})\) gravity, The European Physical Journal C, 2018, 397, DOI: 10.1140/epjc/s10052-018-5886-x