#### Abstract

In this paper we derive a novel cosmological model from the theory of gravitation, for which is the Ricci scalar and is the trace of the energy-momentum tensor. We consider the functional form , with being the Starobinsky Model, named , and , with and being constants. We show that a hybrid expansion law form for the scale factor is a solution for the derived Friedmann-like equations. In this way, the model is able to predict both the decelerated and the accelerated regimes of expansion of the universe, with the transition redshift between these stages being in accordance with recent observations. We also apply the energy conditions to our material content solutions. Such an application makes us able to obtain the range of acceptability for the free parameters of the model, named and .

#### 1. Introduction

The theories of gravity [1, 2] are an optimistic alternative to the shortcomings General Relativity (GR) faces as the underlying gravitational theory, such as those discussed by Padmanabhan [3]; Antoniadis et al. [4]; Demorest et al. [5]; Bull et al. [6]. They can account for the cosmic acceleration [7, 8] with no need for a cosmological constant, providing good match between theory and cosmological observations [9–11].

Particularly, in standard cosmology, derived from GR, the dark energy and dark matter should compose of the universe [12], but their nature is still dubious [13, 14].

Another crucial trouble surrounding GR is the difficulty in quantizing it. Attempts to do so have been proposed [15–17] and can, in future, provide us a robust and trustworthy model of gravity-quantum mechanics unification.

Meanwhile it is worthwhile to attempt to consider the presence of quantum effects in gravitational theories. Those effects can rise from the consideration of terms proportional to the trace of the energy-momentum tensor in the gravitational part of the action, yielding the gravity theories [18]. Those theories were also motivated by the fact that although gravity is well behaved in cosmological scales, the solar system regime seems to rule out most of the models proposed so far [19–22]. Furthermore, rotation curves of spiral galaxies were constructed in gravity, but the results did not favour the theory, as it can be checked in Chiba [23]; Dolgov & Kawasaki [24]; Olmo [25]. The structure and cosmological properties of the modified gravity starting from theory to power-counting renormalizable covariant gravity were presented by Nojiri & Odintsov [26]. The review in Nojiri et al. [27] describes the cosmological developments regarding inflation, bounce, and late-time evolution in , , and modified theories of gravity, with and being the Gauss-Bonnet and torsion scalars.

Despite its recent elaboration, gravity has already been applied to a number of areas, such as cosmology ([28–34] and astrophysics [35–42]). Particularly, solar system tests have been applied to gravity [43] and the dark matter issue was analysed by Zaregonbadi et al. [44]. The late-time behaviour of cosmic fluids consisting of collisional self-interacting dark matter and radiation was discussed by Zubair et al. [45]. Moreover, considering the metric and the affine connection as independent field variables, the Palatini formulation of the gravity can be seen in [46, 47].

By investigating the features of an or model, one realizes the strong relation they have with the functional form of the chosen functions for and , as well as with their free parameters values. In fact, a reliable method to constraint those “free” parameters to values that yield realistic models can be seen in Correa & Moraes [48] and Correa et al. [49] for these theories, respectively.

In the gravity, a reliable and reputed functional form was proposed by A.A. Starobinsky as [50, 51]which is known as Starobinsky Model (SM), with a constant. It predicts a quadratic correction of the Ricci scalar to be inserted in the gravitational part of the Einstein-Hilbert action.

SM has been deeply applied to the cosmological and astrophysical contexts in the literature. Starobinsky showed that a cosmological model obtained from (1) can satisfy cosmological observational tests [51]. On the other hand, the model seems to predict an overproduction of scalarons in the very early universe. In an astrophysical context, SM is also of great importance. In Sharif & Siddiqa [52], the authors have explored the source of a gravitational radiation in SM by considering axially symmetric dissipative dust under geodesic condition. In Resco et al. [53], it has been shown that in SM it is possible to find neutron stars with , which raises as an important alternative to some of the GR shortcomings mentioned above [4, 5]. In Astashenok et al. [54], the macroscopical features of quark stars were obtained.

Our proposal in this paper is to construct a cosmological scenario from an functional form whose dependence is the same as in the SM, i.e., with a quadratic extra contribution of , as in (1). The dependence will be considered to be linear, as , with being a constant. Therefore, we will take

As far as the present authors know, the SM has not been considered for the dependence within models for cosmological purposes so far, only in the study of astrophysical compact objects [55–57] and wormholes [58, 59]. We believe this is due to the expected high nonlinearity of the resulting differential equation for the scale factor. Anyhow, the consideration of linear material corrections together with quadratic geometrical terms can imply interesting outcomes in a cosmological perspective as it does in the astrophysical level.

#### 2. The Gravity

Following the steps in Harko et al. [18], we can write the gravity total action asin which is the determinant of the metric , is the matter Lagrangian, and we are working with natural units.

By taking , with being the pressure of the universe, the variational principle applied in (3) yields the following field equations:

In (4), is the Einstein tensor, is the Ricci tensor, , is the matter-energy density of the universe, and . Still in (4), it can be straightforwardly seen that the limit recovers the GR field equations.

Also, the above choice for the matter Lagrangian is usually assumed in the literature as it can be checked in Harko et al. [18]; Moraes et al. [35, 38], among many others.

#### 3. The Cosmology

Let us assume a flat Friedmann-Robertson-Walker metric in the field equations above. Such a substitution yields the following Friedmann-like equations:where we are using the following definitions:

In the equations above, is the scale factor and dots represent time derivatives. Once again, the limit retrieves the standard formalism.

In terms of the Hubble parameter , the values of and from (5)-(6) arewith the following definitions:For (9)-(10), we considered . Moreover, and , with , are functions of and its time derivatives, expressed by the following:

We can consider, as a solution for (9)-(10), the scale factor in the hybrid expansion law form [60]:with and being constants. It can be seen that such a form for the scale factor consists of a product of power law and exponential law functions. Equation (23) mimics the power law and de-Sitter cosmologies as particular cases and can predict the transition from a decelerated to an accelerated regime of the universe expansion. It has been applied to Brans-Dicke models in Ozgur et al. [60], yielding observational constraints to and .

From (23), the Hubble and deceleration parameters aresuch that the deceleration parameter is defined in such a way that its negative values describe an accelerated expansion of the universe.

We know that the universe not always has accelerated its expansion [7, 8]. The accelerated regime of expansion is considered to be a late-time phenomenon.

In this way, one can choose the constants and in such a way that the power law dominates over exponential law in the early universe and the exponential law dominates over power law at late times. Therefore, the decelerated and accelerated regimes of the universe expansion can be, respectively, well described, as well as the transition between these regimes.

From (25) it is clear that there is a transition from deceleration to acceleration phases of the universe expansion at , with . Since the negativity of the second term leads to a negative time, which indicates an nonphysical situation, we conclude that the cosmic transition may have occurred at .

From , with being the present value of the scale factor and being the redshift, we obtain the following time-redshift relation:where denotes the Lambert function (also known as “product logarithm”).

Using (26), we can plot the deceleration parameter with respect to the redshift, which can be appreciated in Figure 1, in which the values chosen for and are in agreement with the observational constraints found in Ozgur et al. [60].

Plotting as a redshift function has the advantage of checking the reliability of the model, through the redshift value in which the decelerated-accelerated expansion of the universe transition occurs. We will call the transition redshift as and in our model it can be seen that it depends directly on the parameter . From Figure 1, the transition occurs at , corresponding to , respectively. The values of the transition redshift for our model are in accordance with the observational data, as one can check in Capozziello et al. [61, 62]; Farooq et al. [63].

Now, let us write the solutions for the material content of our model, named and . Using (24) in (9)-(10), we havewhere we are using the definitionswith and , with , being functions of only, defined as

The evolution of the energy density, pressure, and corresponding equation of state (EoS) parameter with are shown in Figures 2–4, in which the time units are Gyr.

#### 4. Energy Conditions

Energy conditions (ECs), in the context of a wide class of covariant theories including GR, are relations one demands for the energy-momentum tensor of matter to satisfy in order to try to capture the idea that “energy should be positive”. By imposing so, one can obtain constraints to the free parameters of the concerned model.

The standard point-wise ECs are [64–66] Null energy condition (NEC): . Weak energy condition (WEC): , . Strong energy condition (SEC): . Dominant energy condition (DEC): .

Generally the above ECs are formulated from the Raychaudhuri equation, which describes the behaviour of space-like, time-like, or light-like curves in gravity.

In the present model the energy-momentum tensor has the form of a perfect fluid. So, we will use the above relations for analysing the ECs in theory. The behaviour of the ECs with , and is given in Figures 5–8.

#### 5. Discussion and Conclusions

In this article we have proposed an cosmological model. For the functional form of the function , we investigated a quadratic correction to , as in the well-known SM, together with a linear term on . The substitution of such an function in the gravitational action generates extra terms in the field equations (4) and, consequently, in the Friedmann-like equations (5) and (6). Those extra terms bring some new information regarding the dynamics of the universe, as we discuss below.

Our solution for the scale factor is a hybrid expansion law, well described in Ozgur et al. [60]. From (23), we have obtained the cosmological parameters, namely, Hubble factor and deceleration parameter. Specifically, for the deceleration parameter, we could separate it in two phases: one describing the decelerated and the other describing the accelerated regime of the universe expansion. This can be well-checked in Figure 1, in which the transition redshift between these two stages agrees with observational data. Moreover, one can also note that, remarkably, the present () values for , named for and for , are also in agreement with observational data [12].

We have also obtained solutions for the material content of the universe, named and (see (27)-(41)). In Figures 2–4 we plot the evolution of , , and , the EoS parameter. The values chosen for and respect the energy conditions outcomes presented in Section 4. In Figure 3 we see that the pressure of the universe starts in positive values and then assumes negative values. In standard model of cosmology, a negative pressure fluid is exactly the mechanism responsible for accelerating the universe expansion. In the present model, such a behaviour for the pressure was naturally obtained.

It is also worth stressing that the EoS parameter shows a transition from a decelerated to an accelerated regime of the expansion of the universe. This can well be seen in Figure 4, by recalling that, from standard cosmology, the latter regime may happen if [67]. Moreover, it can be seen from such a figure that as time passes by, , in accordance with recent observational data on the cosmic microwave background temperature fluctuations [12].

Furthermore, in Figures 5–8 we plotted the ECs from the material solutions of our cosmological model. Those figures were plotted in terms of , for fixed . They show the validation of WEC, NEC, and DEC, with a wide range of acceptable values for .

On the other hand, Figure 7 shows violation of SEC. However, as it has been deeply discussed by Barcelo & Visser [68] an early and late-time accelerating universe must violate SEC.

As a further work, we can look for the gravity universe at very early times, particularly investigating the production of scalarons in this model. Moraes & Santos [30] have shown that the trace of the energy-momentum tensor contribution in the theory is higher for the early universe when compared to the late-time contribution. According to Starobinsky [51], some mechanism should work in the early universe to prohibit the scalaron overproduction within SM. The high contribution of the terms proportional to in the early universe may well be this mechanism.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

P. H. R. S. Moraes thanks São Paulo Research Foundation (FAPESP), grant 2015/08476-0, for financial support. P. K. Sahoo acknowledges DST, New Delhi, India, for providing facilities through DST-FIST Lab, Department of Mathematics, where a part of this work was done. R. A. C. Correa is partially supported by FAPESP (Foundation for Support to Research of the State of São Paulo) under grants numbers 2016/03276-5 and 2017/26646-5.