A complete cosmological scenario from $f(R,T^\varphi)$ gravity theory

Recent elaborated by T. Harko and collaborators, the $f(R,T)$ theories of gravity contemplate an optimistic alternative to dark energy, for which $R$ and $T$ stand for the Ricci scalar and the trace of the energy-momentum tensor, respectively. Although the literature has shown that the $T$ dependence on the gravitational part of the action - which is due to the consideration of quantum effects - may induce some novel features in the scope of late-time cosmological dynamics, in the radiation-dominated universe, when $T=0$, no contributions seem to rise from such theories. Apparently, $f(R,T)$ contributions to a radiation-dominated universe may rise only from the $f(R,T^\varphi)$ approach, which is nothing but the $f(R,T)$ gravity in the case of a self-interacting scalar field whose trace of the energy-momentum tensor is $T^\varphi$. We intend, in this article, to show how $f(R,T^\varphi)$ theories of gravity can contribute to the study of the primordial stages of the universe. Our results predict a graceful exit from inflationary stage to a radiation-dominated era. They also predict a late-time cosmic acceleration after a matter-dominated phase, making the $f(R,T^\varphi)$ theories able to describe, in a self-consistent way, all the different stages of the universe dynamics.


Introduction
Plenty of efforts have been made in the theoretical framework with the purpose of explaining the accelerated regime our universe has passed through a fraction of a e-mail: moraes.phrs@gmail.com b e-mail: joaorafael@uaf.ufcg.edu.br second after the Big-Bang, named "inflationary era". It is common to describe this phenomenon via scalar fields [1]- [7], although another reputed form to do so comes from the f (R) gravity [8]- [11].
Recently, a more general theory of gravity, named f (R, T ) gravity, was proposed by T. Harko and collaborators [12]. A priori it can also contribute to inflationary era studies through its scalar field approach, named f (R, T φ ) gravity [12].
Although the late time acceleration of the universe expansion [13]- [15] has been broadly investigated in such a theory of gravity [16]- [25], the inflationary era still presents a lack of examination. The same happens for the radiation-dominated era of the universe.
For the radiation era, the lack of f (R, T ) applications is quite predictable for the following reason. The equation of state (EoS) of the universe is p = ρ/3 at this stage, with p and ρ representing the pressure and density of the universe, respectively. Such an EoS yields a null trace for the energy-momentum tensor of a perfect fluid. A null trace for the energy-momentum tensor yields the f (R) formalism for the f (R, T ) functional forms found in the literature (check, for instance, [12], [16]- [25]). Hence the study of the radiation era of the universe in f (R, T ) gravity seems to be quite tricky.
A priori, it seems reasonable to affirm that the f (R, T ) gravity does not contribute to the study of the radiationdominated universe, since the contribution coming from the trace of the energy-momentum tensor in f (R, T ) vanishes at this stage. Specifically, such a shortcoming or incompleteness has attracted attention recently [23]- [26].
In [23], the authors argued that the high redshift f (R, T ) cosmological solutions tend to recover the standard model of cosmology if the f (R, T ) functional form is linear in R precisely because when z >> 1, the radi-ation with EoS p = ρ/3 dominates the universe dynamics, making the trace of the energy-momentum tensor to vanish. Therefore, no novel contributions would come from the f (R, T ) theories of gravitation at the radiation era.
In [24], in order to make f (R, T ) gravity theory able to contribute also to a radiation-dominated universe, a compactified space-like extra dimension had to be invoked.
In [25], with the purpose of obtaining well behaved solutions at the radiation stage of the universe, the f (R, T ) gravity was generalized, by allowing the speed of light c to vary. Moreover, a scenario alternative to inflation was obtained.
By developing the f (R, T φ ) approach, which is submissive to f (R, T ) gravity, we intend, here, to make the f (R, T ) theories able to contribute also to inflationary and radiation-dominated eras, in a self-consistent way. Recall that once T = 0 at the radiation stage, it is natural to expect that f (R, T ) theories retrieve f (R) gravity or even general relativity (GR) when f (R) ∼ R, in such a way one would not expect any new informations derived from f (R, T ) gravity. To make the f (R, T ) gravity able to describe a radiation universe by itself, without necessarily recovering the f (R) formalism outcomes, will make such a theory able to contribute to all the different stages of the universe dynamics. In order to do so, we will implement the first-order formalism [27,28] to the f (R, T φ ) gravity.
We will not omit the recent dynamical features our universe has presented. In other words, a late-time cosmic acceleration, which is corroborated by Type Ia Supernovae observations [13,14] and described in ΛCDM cosmology by the presence of the cosmological constant in standard Einstein's field equations (FEs), shall also be described in our model.

The f (R, T φ ) Gravity
Early and late time accelerations of the universe expansion can be explained in the framework of scalar field models [1]- [7], [29]- [35]. The consideration of f (R, T ) gravity in the case of a self-interacting scalar field φ yields the f (R, T φ ) gravity [12].
Here we will work with the case f (R, T φ ) = −R/4+ f (T φ ), for which R stands for the Ricci scalar, while T φ is the trace of the energy-momentum tensor of the scalar field, so that is the action of such a field, with g being the determinant of the metric with signature (+, −, −, −) and throughout this article we will assume units such that 4 π G = c = 1. Moreover, we are going to deal with a standard Lagrangian density for a real scalar field, whose form is for which it has been considered a self-interacting potential V (φ). Furthermore, since the energy-momentum tensor of the scalar field is given by its trace is The f (R, T φ ) = −R/4 + f (T φ ) gravity model is defined in the following form: for which we have inserted the gravitational part of the action. By minimizing the above action, we obtain where G µν is the Einstein tensor and f ′ is the derivative of f (T φ ) with respect to T φ . For the action (5), a flat Friedmann-Robertson-Walker universe has the following Friedmann equations: in which φ = φ(t), H =ȧ/a is the Hubble parameter, a is the scale factor and dots denote derivation with respect to time.
We are also able to derive the equation of motion for such a system, which is given by where it was used the notation V φ ≡ dV (φ)/dφ.

First-Order Formalism
Let us implement the first-order formalism to the model above, by following the recipe presented in [27,28]. Firstly it is straightforward to observe from (7) and (8) that the Hubble parameter obeys the differential equatioṅ Then, the initial ingredient that we may consider in order to establish the first-order formalism is based on the definitioṅ with W φ ≡ dW (φ)/dφ and W is an arbitrary function of the field φ. Such a procedure is commonly used in cosmological scenarios coupled with scalar fields, as we may see in [29] and references therein. This assumption implies that (10) is rewritten aṡ If we establish that where c is a real constant, then, Eq.(10) requires the following constraint to the function h: Moreover, the Friedmann equations impose that the potential V (φ) has to obey the relation

Examples
We can apply the first-order formalism by considering a specific form for the function f (T φ ). As in [28], let us deal with f (T φ ) = α (T φ ) n , with α as a constant and n as an integer. Note that such a functional form for f (T φ ) straightforwardly yields the GR case when α = 0. Moreover, it is the analogous of the form f (T ) = αT n in the f (R, T ) gravity with no scalar field case, firstly suggested in [12] for deriving an accelerated cosmological scenario and then applied to a number of other f (R, T ) well-behaved cosmological models, such as [18,19], [24,25], among many others. Such an assumption yields In this study, we are going to consider that the function W is given by with a 1 being a real constant. Such a form for W (φ) is a sine-Gordon type of model, and it has been broadly studied in the literature, specially in investigations related with classical field theory, as one can see in [36] and references therein. Very recently, this model was applied in the context of braneworld scenarios for the f (R, T ) gravity [28].
Here, for a matter of simplicity, we focus on the case n = 1, which leads to Then, by substituting f ′ in Eq. (14), we find which means that Moreover, it is straightforward to obtain that the po- From Eq.(18) one can see that the first-order differential equation for this scenario has the forṁ whose analytical solution is with b 1 being a real constant. The previous solution allows us to rewrite Eq. (21) as which is presented in details in the upper panel of Fig.1 below, where we can see that H exhibits a kink-like profile as a direct consequence of the model defined in Eq. (18).
We can use the last result to determine the expansion parameter (or scale factor) a(t), whose explicit form is given by Moreover, the analytical EoS parameter is and its features can be observed in the central panel of Fig.1. We can verify from the analytical cosmological parameters that small deviations on the values of the constants produce cosmological scenarios which are similar to those presented in Fig.1. Therefore, the time evolution of our parameters are not so strongly dependent on the initial conditions of the model. Another interesting model is whose first-order differential equation is such thaṫ which is satisfied by where b 1 is a real constant. This last first-order differential equation represents the so-called φ 4 type of model, which is also broadly investigated in classical field theory, as one can see in [37] and references therein.
In the present case, we are dealing with f (T φ ) = α T φ again. Therefore, we can combine the previous ingredients to determine the following form for the function h(φ): The last equation results in the following relation for the Hubble parameter which is presented in details in the upper panel of Fig.2 below. As in the previous example, it is possible to choose the arbitrary constants in order to obtain a kinklike profile for H. We can also show that the potential V for this analytical model is Withal, the expansion parameter for such a case is given by the following while the EoS parameter reads The latter is plotted in the central panel of Fig.2. Again, if we take small deviations of these parameters we determine similar cosmological scenarios.

Cosmological Interpretations
In Section 3, we were able to depict the time evolution of H and ω derived from the first-order formalism applied to the f (R, T φ ) theory. Remind that our aim in this work is to make f (R, T ) gravity able to induce a complete cosmological scenario, which includes the primordial eras of the universe, since the theory seems to present some shortcomings at these stages. Therefore, it is worth checking if the cosmological parameters derived previously indeed resemble, for small values of time, what is expected for the primordial universe dynamics.
In the previous section we have also plotted the evolution of the trace of the energy-momentum tensor of the scalar field, which shall be interpreted below.
The inflationary model states that early in the history of the universe, its expansion has accelerated. From standard Friedmann equations, this happens if the universe is dominated by a component with EoS ω < −1/3 [38,39]. According to standard cosmology, during inflation, the energy density of the universe was dominated by a constant, say, Λ ι , in such a way the Friedmann equations read H 2 ι ∼ Λ ι . The Hubble parameter was, then, constant during inflationary era. Indeed, one may write, for inflation, a ι ∼ e Hιt [38,39].
From the upper panels of Figs.1-2, the Hubble parameter values predicted from our model indeed remain approximately constant during a small period of time after the Big-Bang, in agreement with inflationary scenario. In fact,Ḣ ∼ −0.003 and ∼ −0.0003 at t = 0.01 for Figs.1 and 2, respectively.
After inflation, H must evolve as ∝ t −1 H , with t H being the Hubble time [38,39], in accordance with Figs.1-2.
Still in the upper panels of Figs.1-2, one should note that for high values of time, H behaves once again as a constant. From standard cosmology, it is well known that for high values of time, H is proportional to ρ Λ (the density of the cosmological constant), which is, indeed, a constant. Therefore, the recent cosmic acceleration our universe is undergoing can be described by a scale factor which evolves as e H0t , with the constant behaviour of H being predicted in Figs.1-2 for high values of time.
In the central panels of Figs.1-2, we have depicted the evolution of the EoS parameter ω with time. By analysing it, one can see that for small values of time, ω ∼ −1. This is in agreement with some constraints that have been put to inflationary EoS recently [40]- [42].
After representing inflation, ω smoothly evolves to ∼ 1/3. This is the maximum value the EoS should assume during the universe evolution. Moreover, 1/3 is precisely the value of the EoS which describes the radiation-dominated era of the universe [38,39], i.e., the model predicts a smooth transition from inflationary stage to radiation-dominated era.
As time passes by, the central panels of Figs.1-2 show that ω → −1. According to recent observations of anisotropies in the temperature of cosmic microwave background radiation [15], this is the present value of the universe EoS, and is responsible for the recent cosmic acceleration.
The fact that nowadays the universe is passing through its second phase of accelerated expansion justifies such an evolution for ω. Since an accelerated universe is described by a negative EoS, it is expected that during the whole universe evolution, the value of ω, coming from ∼ −1, which is related to the inflationary stage, first increases, tending to a decelerating phase, then decreases, returning to an accelerated phase.
Furthermore, the lower panels of Figs.1-2 show the time evolution of the trace of the energy-momentum tensor of the scalar field for the models.
Firstly, one can note that the values assumed by such a quantity are always ≥ 0, as they, indeed, should be.
The high values of density and temperature which are known to characterize the Big-Bang, together with the inflationary EoS mentioned above, make us to expect to have the maximum values of the trace of the energy-momentum tensor of the scalar field during inflation, as corroborated by Figs.1-2.
At the radiation era, T φ must vanish. We zoomed in on Figs.1-2 lower panels in order to explicit such a predicted feature. Note that after a period assuming its maximum values, the trace of the energy-momentum tensor T φ → 0 at the same time scale ω → 1/3 in Figs.1-2 central panels. Such a property reinforces the good behaviour of our solutions, specifically at the primordial stages of the universe evolution.
At matter-dominated era, T φ must be = 0. The same happens for the second period of acceleration. Indeed, the plots of T φ show that after the radiation era, T φ increases its values. Still, such values are small, due to the low density of the universe at these stages, specially when compared to the ρ values at inflation.
Remarkably, T φ remains constant for high values of time. Since T φ ∝ ρ, this characterizes the late-time universe dynamics to be dominated by a constant. Recall that in standard cosmology, for high values of time, the cosmological constant Λ dominates the universe dynamics. As argued above, ρ Λ remains, indeed, constant as time passes by. Such an important feature of standard cosmology is predicted from the f (R, T φ ) approach here presented.

Discussion and perspectives
In this article we established a first-order formalism to determine two analytical models related to the f (R, T φ ) theory. Such a formalism was constructed from the definition observed in Eq. (11).
The mentioned expression together with Eq. (7) and (13) gave us a constraint for the potential V , which can be viewed in (15). Moreover, the arbitrary function h(φ) was introduced and its form is directly related to f ′ . The first analytical model was generated by considering W = a 1 sin φ and the specific form for f (T φ ) showed in (16). Such definitions leaded us to determine h(φ), besides the potential V , the Hubble parameter H, as well as other cosmological quantities like the expansion and EoS parameters. In the second case, we adopted W = a 1 (φ 3 /3 − φ) , and one more time we determined the function h and the cosmological parameters for the model.
We can directly see that the field φ(t) satisfies the equation of motion (9) for both models. Furthermore, this first-order formalism may help us to obtain other one-field analytical scenarios and also may lead us to hybrid analytical models, where it is considered the coupling between an f (R, T φ,χ ) function with a two scalar field Lagrangian. Such a scenario shall be derived and presented in near future.
Another further work may rise from an alternative form for f (R, T φ ), such as f (R, T φ ) = g(R)h(T φ ), with g(R) and h(T φ ) being functions of R and T φ , respectively. Such a functional form suggests a high coupling between matter and geometry, and therefore may imply some novelties in gravitational and cosmological perspectives. We can conjecture that by coupling such kind of theories with a scalar field lagrangian may lead us to a new set of Friedmann equations and to different equations of motion. Once we have these ingredients we may try to apply the first-order formalism and then compute the constraints that the superpotential will obey. In fact, the analogous of that model for the f (R, T ) gravity with no scalar field case was investigated in [16].
We have obtained from the first-order formalism applied to the f (R, T φ ) gravity, two models able to describe all the different dynamical stages of the universe, from inflation, to radiation, matter and dark-energy dominated eras.
For the first time, an f (R, T ) cosmological model was able to describe all the dynamical stages of the universe continuously, including the primordial ones. There was no need of invoking different EoS or functional forms for f (T φ ) in different stages in order to be able to account for all the transient eras of the universe dynamics, as made in [26], [43] and [44], for instance. In fact, there was no need of invoking any EoS at all. Such a quantity was obtained simply by dividing the solutions p by ρ.
It is remarkable the fact that the inflationary scenarios of our models have smoothly decayed to a radiationdominated era. This phenomenon is sometimes called "graceful exit" and it is not trivial to be obtained (check [45,46] for the graceful exit achievement through decaying vacuum models). In fact, an ideal inflationary scenario must naturally decay to a radiation-dominated era [47,48], as in the present case.
Also remarkable is the fact that the models are wellbehaved at radiation stage. The match between the time scales for ω = 1/3 e T φ = 0 in Figs.1-2 reinforces that.
Not just the present models were able to describe a smooth transition between inflationary and radiationdominated eras, but they were also able to predict the matter dominated phase, with ω = 0, and the recent cosmic acceleration the universe is undergoing [13]- [15], from an EoS parameter whose late-time values agree with WMAP cosmic microwave background observations [15].
Therefore, departing from a cosmological scenario derived from GR, we have obtained an EoS which varies with time, allowing, through its evolution, the universe dynamics to be dominated by different components in a natural form. Moreover, in order to obtain a recent cosmic acceleration, we did not have to insert a cosmological constant in the model FEs, in this way, evading the cosmological constant problem [49,50,51].
The f (R, T φ ) cosmological scenario derived in this work is somehow similar to a decaying vacuum Λ(t) model, such as those presented in [45,46], for instance. Firstly, it should be stressed that recently in [26], the f (R, T ) theories have been interpreted as a generalization of the Brans-Dicke model [52]- [56]. As argued by the authors in [26], the extra terms in the f (R, T ) FEs may play the role of additional terms of effective gravitational and cosmological "constants" G ef f and Λ ef f . Those extra terms in Λ ef f are the responsible for a vacuum decay in Λ(t)-like models.
Second, in both decaying vacuum and f (R, T ) models, the covariant derivative of the energy-momentum tensor ∇ µ T µν does not vanish. In decaying vacuum models, this requires necessarily some energy exchange between matter and vacuum, through vacuum decay into matter or vice versa.
The vacuum-decay models presented in [45,46], for instance, are capable of describing the complete history of the universe in an effective unified framework, in the same way the present f (R, T φ ) model does. In those models, the vacuum always decays in the dominant dynamical component of each stage. However, in order to be able to describe the various phases of the universe dynamics, the functional form taken for Λ(t) is proposed on phenomenological grounds, by extrapolating backwards in time the present available cosmological data. On the other hand, the f (R, T φ ) scenario is free of ansatz or phenomenology.
Note also that the Λ(t) model presented in [45] was mimic through a scalar field model (check Section 4 of such a reference). This corroborates the argument that it might have some fundamental analogy between scalar field and decaying vacuum models.
In this article, we have also formally brought to the scientific community attention the shortcoming surrounding the f (R, T ) theories of gravity in the regime T = 0. It is, indeed, straightforward to realize that in such a regime of f (R, T ) gravity one automatically recovers the f (R) theories, or GR when f (R) is linear in R [23]. We have shown, as an unprecedented outcome, that it is possible to describe radiation in f (R, T ) gravity not necessarily by the recovering of GR, but through extra dynamical terms coming from a scalar field in the f (R, T φ ) approach.
If one focus specifically on the radiation era predicted by the f (R, T φ ) model here presented, one is able to obtain a scale factor which evolves as the one of standard cosmology purely from the f (T φ ) contribution, which is a complementary novel feature of the f (R, T ) gravity that we will show below.
By taking T φ = 0 in (4) yieldsφ 2 = 4V (φ). Since in this case f (T φ ) = α(T φ ) n = 0, one obtains, from (10), From Eqs. (11), (15) and (36), one can writeḢ = −2H 2 , which integrated for the scale factor a yields a ∼ t 1/2 . From standard cosmology [38,39], this is exactly the time proportionality of the scale factor of a radiationdominated universe. Therefore the formalism presented in this article surpasses the f (R, T = 0) issue since it makes f (R, T ) gravity able to describe a radiationdominated universe purely from its dependence on T φ . Besides some questions exposed in the Introduction, very recently, the drawback in f (R, T ) theory was raised in [26], in which the authors have searched for different forms for f (T ) in different stages of the universe. In the radiation-dominated stage, they found that f (T ) should be a non-null constant. As the authors have considered the f (R, T ) theories as a generalization of Brans-Dicke gravity, the f (T ) ∼ cte case leads to G ef f and Λ ef f to be, indeed, constants, as it may be seen in [26]. Therefore no contributions to a T = 0 universe would come from the f (R, T ) formalism on this approach, departing from the results of the present work.
One could argue that the f (R) formalism retrieval in f (R, T ) theories when T = 0 is not a shortcoming, being just an intrinsic property of such theories. However, below we bring to the reader's attention some fundamental areas of knowledge which will have to be neglected if the f (R, T ) theories, indeed, simply recovers the f (R) formalism when T = 0.
Many contributions to Physics, Astrophysics and Cosmology will come from the imminent detection of gravitational waves (GWs) [57,58]. Their detection will clarify some persistent and important issues, such as the absolute ground state of matter [59], and also contemplates a powerful tool in estimating parameters of compact binary systems [60], constraining the equation of state of neutron stars [61] and brane cosmological parameters values in braneworld models [62], and distinguishing GR from alternative theories of gravity [63].
The study of the propagation of GWs is made once they have been generated [64]. The energy-momentum tensor T µν is then set to zero in order to obtain the wave equation for vacuum [65]- [67]. In the same way the EoS ω = 1/3 yields a null trace for the energymomentum tensor, this will happen for vacuum, with T µν = 0. Therefore it is expected that the study of GWs, as the derivation of their primordial spectrum, for instance, in the f (R, T ) formalism yields the same predictions as those obtained via f (R) gravity [68].
A powerful tool to study the properties of GWs was developed in [69]. It consists in characterizing the polarization states of GWs in a given theory 1 . In order to do so, one evaluates the Newman-Penrose (NP) quantities [70,71]. In [63], from the NP quantities calculation, the authors have found extra polarization states for the f (R) gravity. Those extra states may be corroborated by experiment once the advanced GW detectors, such as LIGO [57] and VIRGO [58], start they next run. Since the NP formalism is applied to the vacuum FEs of a given theory, such an application to the f (R, T ) theory is expected to yield the same extra polarization states as those of the f (R) formalism [63].
The Kerr metric describes the geometry of an empty space (T µν = 0) around a rotating uncharged black hole (BH) with axial symmetry. Kerr BHs studies have been constantly seen in the literature (check, for instance, [72]- [75]). Since BHs are characterized by strong gravitational fields in quantum length scales, it is expected that their physics bring up some new insights about quantum theories of gravity [76,77]. Once again, one would not obtain new contributions to the study of the space-time around a Kerr BH from the f (R, T ) formalism, in such a way the f (R) theories outcomes should be expected [78]- [81].
Hence, one can see that the lack of contributions from f (R, T ) gravity will appear not only in a radiationdominated universe, but also in studies made in the vacuum regime, since T µν = 0 → T = 0. In other words, in the absence of the f (R, T φ ) approach presented here, one would expect the f (R) formalism features retrieval in the study, not only of a radiation-dominated universe, but also in those of GWs and space-times surrounding BHs.
For instance, in the search for new polarization states of GWs mentioned above, one works with the vacuum FEs of a theory. It is expected that if one considers the vacuum on the rhs of standard f (R, T ) FEs (check Eq.(11) of [12]), one gets, after the NP formalism application, the same polarization states of the f (R) theory [63]. On the other hand, by taking T µν = 0 (and con-sequently T = 0) in the present model FEs (6) 2 , one expects the third term inside the brackets to survive in the vacuum regime. This term should naturally yield some departure from both GR and f (R) gravity in the study of GWs polarization states. Moreover, it could be responsible for novelties in other astrophysical areas submissive to the GW subject, such as the characterization of the GW spectrum [64], the evolution of cosmological perturbations [82] and the GW background from cosmological compact binaries [83].
Furthermore, another vacuum studies in f (R, T ) gravity will be allowed to be investigated, such as the Tolman-Oppenheimer-Volkoff-de Sitter equation [84] and the vacuum polarization by topological defects [85], among many others. Those points above justify the search for evading the f (R, T = 0) issue, which was presented in this article.
The contributions from the present work are able to motivate new branches in the study of the radiation universe and vacuum in f (R, T ) gravity. The present formalism makes possible to make the most of the f (R, T ) theories. Although the f (R) formalism retrieval in the T = 0 regime would not necessarily be a critical feature of the f (R, T ) theories, it would make us overlook plenty of applications of those theories, which may generate some remarkable outcomes. Some of those, as possible extra polarization states of GWs, may be inferred by experiment soon.