arXiv is now an independent nonprofit! Learn more
License: CC BY 4.0
arXiv:2406.16113v1 [gr-qc] 23 Jun 2024
00footnotetext: PACS: 04.50.Kd, 98.80.jk, 98.80.cq.

:

Pseudo generalized Ricci-recurrent spacetimes with certain applications to modified gravity

Uday Chand De1 and Krishnendu De * 2 Address: 1 Department of Pure Mathematics, University of Calcutta, West Bengal, India. ORCID iD: https://orcid.org/0000-0002-8990-4609 Email address: uc$˙-$de@yahoo.com, ucde1950@gmail.com Address: 2 Department of Mathematics,
Kabi Sukanta Mahavidyalaya, The University of Burdwan.
Bhadreswar, P.O.-Angus, Hooghly,
Pin 712221, West Bengal, India.
ORCID iD: https://orcid.org/0000-0001-6520-4520
Email address: krishnendu.de@outlook.in, krishnendu.de67@gmail.com
Abstract.

In this article we introduce and characterize a pseudo generalized Ricci-recurrent spacetimes. At first, we produce an example to justify the existence of such a spacetime. Then, it is provided that a pseudo generalized Ricci-recurrent generalized Robertson-Walker spacetime represents a perfect fluid spacetime and a pseudo generalized Ricci-recurrent perfect fluid spacetime represents either a dark energy epoch of the Universe or, the velocity vector field is parallel, conservative, acceleration-free, vorticity-free, and shear-free and becomes a static spacetime. Lastly, we study the impact of this spacetime under f()f(\mathcal{R}) gravity scenario and deduce several energy conditions.

Key words and phrases: 
Pseudo generalized Ricci-recurrent spacetime; generalized Robertson-Walker spacetime; perfect fluid spacetime; static spacetime; f()f(\mathcal{R}) gravity.
*Corresponding author.

1. Introduction

Let MnM^{n} (n4n\geq 4) be a semi (or, pseudo)-Riemannian manifold and gg be its semi-Riemannian metric with signature (p,m)(p,m), such that p+m=np+m=n. MnM^{n} equipped with gg is referred to as a Lorentzian manifold [1] if gg is a Lorentzian metric with signature (n1,1)(n-1,1) or, (1,n1)(1,n-1). Time-oriented Lorentzian manifolds are Lorentzian manifolds that accept a globally time-like vector field, also physically referred to as spacetime. A generalized Robertson-Walker (shortly, GRW) spacetime is a Lorentzian manifold of dimension nn (n4)\left(n\geq 4\right) which can be shaped by a warped product I×Ψ2M-I\times_{\varPsi^{2}}\stackrel{{\scriptstyle\ast}}{{M}}, M\stackrel{{\scriptstyle\ast}}{{M}} stands for (n1)\left(n-1\right)-dimensional Riemannian manifold, II\in\mathbb{R} (set of real numbers) is an open interval, and Ψ>0\varPsi>0 is the warping function. This notion was presented by Alías et al. [2] in 19951995. In particular, GRW spacetime becomes Robertson-Walker ( shortly, RW) spacetime if we assume that M\stackrel{{\scriptstyle\ast}}{{M}} is a three dimensional Riemannian manifold of constant sectional curvature. For more details about GRW spacetimes, we refer ([2]-[7]).

A 44-dimensional Lorentzian manifold M4\mathrm{M}^{4} is described as a perfect fluid spacetime (shortly, PFS) if the Ricci tensor lk\mathcal{R}_{lk} fulfills

lk=αglk+βuluk,\mathcal{R}_{lk}=\alpha g_{lk}+\beta u_{l}u_{k}, (1.1)

α\alpha, β\beta being scalars (not simultaneously zero) and the velocity vector uku_{k} stands for a unit time-like vector {that is, ukuk=1u_{k}u^{k}=-1, uk=glkulu^{k}=g^{lk}u_{l}}. The matter field in general relativity (in short, GR) is demonstrated by TlkT_{lk}, which is called the energy momentum tensor (in short, EMT). Since heat conduction terms and stress terms corresponding to viscosity does not occur, the fluid is referred to as perfect [8]. The EMT [1] for a PFS has the following form

Tlk=(σ+p)uluk+pglk,T_{lk}=\left(\sigma+p\right)u_{l}u_{k}+pg_{lk}, (1.2)

pp, σ\sigma being the isotropic pressure and the energy density, respectively. If σ=p\sigma=p, then the PFS is named stiff matter fluid. If σ+p=0\sigma+p=0, p=0p=0 and σ=3p\sigma=3p, then the PFS is called the dark energy epoch of the Universe, dust matter fluid and radiation era [9], respectively. Without the cosmological constant, the Einstein’s field equations (shortly, EFEs) are stated as

κTlklk+2glk=0\kappa T_{lk}-\mathcal{R}_{lk}+\dfrac{\mathcal{R}}{2}\,g_{lk}=0 (1.3)

in which κ\kappa and \mathcal{R} stand for the gravitational constant and Ricci scalar, respectively. Equations (1.1), (1.2) and (1.3) infer that

α=κ(σp2)andβ=κ(σ+p).\alpha=\kappa\left(\dfrac{\sigma-p}{2}\right)\quad\mathrm{and}\quad\beta=\kappa\left(\sigma+p\right). (1.4)

In M4\mathrm{M}^{4}, the Weyl conformal curvature tensor 𝒞ijkl\mathcal{C}^{l}_{ijk} is demonstrated by

𝒞ijkl=ijkl12{gijklgikjl+δklijδjlik}+6{δklgijδjlgik}\mathcal{C}^{l}_{ijk}=\mathcal{R}^{l}_{ijk}-\dfrac{1}{2}\left\{g_{ij}\mathcal{R}^{l}_{k}-g_{ik}\mathcal{R}^{l}_{j}+\delta^{l}_{k}\mathcal{R}_{ij}-\delta^{l}_{j}\mathcal{R}_{ik}\right\}+\dfrac{\mathcal{R}}{6}\left\{\delta^{l}_{k}g_{ij}-\delta^{l}_{j}g_{ik}\right\} (1.5)

in which ijkl\mathcal{R}^{l}_{ijk} denotes the curvature tensor and kl=jkglj\mathcal{R}^{l}_{k}=\mathcal{R}_{jk}g^{lj}.

It is well-known that [10]

l𝒞ijkl=12{(kijjik)16(gijkgikj)}.\nabla_{l}\mathcal{C}^{l}_{ijk}=\dfrac{1}{2}\left\{\left(\nabla_{k}\mathcal{R}_{ij}-\nabla_{j}\mathcal{R}_{ik}\right)-\dfrac{1}{6}\left(g_{ij}\nabla_{k}\mathcal{R}-g_{ik}\nabla_{j}\mathcal{R}\right)\right\}. (1.6)

It is important to note that the Weyl geometry initially presented as a unifying theory of gravity and electromagnetic [11] that could challenge Einstein’s theory of gravity which is the main motivation behind its investigation in this article.

For the purpose of determining Einstein’s field equations, the study of spacetime symmetries is important. Symmetries are characteristic of geometry and reveal physics. There are many symmetries in matter and spacetime geometry. The metric equations are useful because they make many solutions easier. In GR, their primary utilization is that they classify solutions of Einstein’s field equations.

In 1952, Patterson [12] invented the idea of Ricci-recurrent manifolds. In 1995, De et al. [13] introduced the concept of generalized Ricci-recurrent manifolds. In a recent paper [14] Mallick et al. studied generalized Ricci-recurrent manifolds with applications to relativity. The significance of the Generalized Ricci-recurrent structure and its interaction with the modified f()f(\mathcal{R})-theory [15] and the modified Gauss-Bonnet f(,G)f(\mathcal{R},G)-theory [16] are well established. On the other hand, pseudo ZZ-symmetric spacetimes have been investigated by Mantica and Suh [17] and Ozen[18] have studied mm-projectively flat spacetimes. Also, in [19] Ψ\Psi-conformally symmetric spacetimes have been investigated. Moreover, in [20], we have studied ψ\psi-conharmonically symmetric spacetimes. As well, many authors have looked at the spacetime of general relativity in various methods; for additional information, see ( [21], [22]). Inspired by foregoing studies we introduce and characterize a novel spacetime named pseudo generalized Ricci-recurrent spacetimes.

If the covariant derivative of lk\mathcal{R}_{lk} is in the following form:

ilk=uilk+viglk+wi𝒟lk,\nabla_{i}\mathcal{R}_{lk}=u_{i}\mathcal{R}_{lk}+v_{i}g_{lk}+w_{i}\mathcal{D}_{lk}, (1.7)

where uiu_{i}, viv_{i} and wiw_{i} are non-zero covariant vectors, 𝒟\mathcal{D} is the (0,2)\left(0,2\right)-type structure tensor of the manifold obeying

𝒟lk=𝒟kl,glk𝒟lk=0and𝒟lkuk=0,\mathcal{D}_{lk}=\mathcal{D}_{kl},\quad g^{lk}\mathcal{D}_{lk}=0\quad\mathrm{and}\quad\mathcal{D}_{lk}u^{k}=0, (1.8)

then we named such an nn-dimensional manifold as pseudo generalized Ricci-recurrent manifold, denoted by P(GRn)P\left(GR_{n}\right). If wi=0w_{i}=0, then P(GRn)P\left(GR_{n}\right) reduces to a generalized Ricci-recurrent manifold (GRn)\left(GR_{n}\right) [13]. For vi=wi=0v_{i}=w_{i}=0, P(GR4)P\left(GR_{4}\right) becomes a Ricci-recurrent manifold (Rn)\left(R_{n}\right) [12]. If lk\mathcal{R}_{lk} obeys (1.7) and uiu_{i} is a unit time-like vector, then M4\mathrm{M}^{4} is called a P(GR4)P\left(GR_{4}\right) spacetime.

Multiplying (1.7) with glkg^{lk} and using (1.8), we reach

i=ui+4vi.\nabla_{i}\mathcal{R}=\mathcal{R}u_{i}+4v_{i}. (1.9)

The physical motivation for researching various spacetime models in general relativity and cosmology is to gain additional insight into particular phases of the universe’s evolution, which can be divided into the following three stages:
(i) The initial stage, (ii) The intermediate stage and (iii) The final stage.

The initial stage concerns viscous fluid whereas the intermediate stage concerns non-viscous fluid and both are admitting heat flux. The final stage equipped with thermal equilibrium tells about the perfect fluid stage. In our current study, we select the final stage.

The scientific world as a whole accepts the idea that our Cosmos is currently going through an accelerated phase. EFEs are not adequate to determine the late-time inflation of the cosmos without supposing the existence of certain unseen components that could account for the dark energy and dark matter origins. It is the primary source of inspiration for the extension to acquire higher order field equations of gravity.

According to GR theory, it is commonly accepted that energy conditions (ECs) are essential resources for studying black holes and wormholes in various modified gravities ([23]-[27]). The Raychaudhuri equations [28], which methodically produce the ECs, express the intriguing nature of gravity through the positivity condition lkuluk0\mathcal{R}_{lk}u^{l}u^{k}\geq 0, ulu^{l} is a null vector. The last geometric criterion is the same as the null energy condition (NEC) Tlkuluk0T_{lk}u^{l}u^{k}\geq 0. Certainly, the weak energy condition (WEC) reflects that Tlkuluk0,T_{lk}u^{l}u^{k}\geq 0, for every time-like vector ulu^{l} and preserves a positive local energy density. In addition, a spacetime fulfills the dominant energy condition (DEC) if Tlkulvk0T_{lk}u^{l}v^{k}\geq 0 holds for every two co-oriented time-like vectors uu and vv and strong energy condition (SEC) [29] if Rlkuluk0R_{lk}u^{l}u^{k}\geq 0 holds for all time-like vectors uu. An overview of the four primary ECs is given here. A number of additional, lesser-known point wise ECs exist that express different restrictions on the stress energy tensor. A trace energy condition (TEC) was also present in the past. This implies that for a perfect fluid, σ3p0\sigma-3p\geq 0. For a number of years in the 1950s and 1960s, this was thought to be a physically feasible state. Now, opinions have changed. Particularly, the TEC is violated with the finding of stiff equations of state for matter found in neutron stars [30]. We bring it up here as a specific illustration of an EC.

Interestingly, the idea of f()f(\mathcal{R})-gravity appears as a spontaneous extension of Einstein’s theory of gravity. Here, the function f()f(\mathcal{R}), in which \mathcal{R} denotes the Ricci scalar, modifies the Hilbert-Einstein action term. The aforementioned theory was invented by Buchdahl [31], and Starobinsky [32] has demonstrated its validity through research on cosmic inflation. Through the introduction of certain couplings between the geometrical quantities and the matter sector, the f()f(\mathcal{R}) theory of gravity has been further generalized. The non-minimal coupling between the matter lagrangian density and the curvature invariant has been proved in [33], which is known as the f(,Lm)f(\mathcal{R},L_{m}) theory of gravity. The corresponding Lagrangian can be modified by incorporating an analytic function of TjkTjkT_{jk}T^{jk} in this generalization procedure for the f(,Lm)f(\mathcal{R},L_{m}) theory. f(,T2)f(\mathcal{R},T^{2}) gravity or energy-momentum squared gravity is the result of selecting the corresponding Lagrangian. In 2014, Katirci and Kavuk [34] originally put forward this theory, which allows the existence of a term in the action functional that is proportional to TjkTjkT_{jk}T^{jk}. The f(,G)f(\mathcal{R},G)-gravity theory [25] was one of these modified theories. It was developed by changing the previous Ricci scalar \mathcal{R} by a function of \mathcal{R} and GG, the Gauss-Bonnet invariant. The f(,T)f(\mathcal{R},T)-gravity theory, discovered by Harko et al. [26], was another modified theory. This is an extension of f()f\left(\mathcal{R}\right)-gravity ([23], [24]) in which the trace TT of the EMT is directly linked to any arbitrary function of \mathcal{R}. Several functional forms of f()f(\mathcal{R}) were previously provided in the following works: ([35]-[39]) which shows that this modified theory has a number of cosmological applications.

The literature mentioned above makes it very clear that more attention needs to be paid to f()f(\mathcal{R}) gravity, and there are still a lot of unanswered questions. Inspired by the foregoing investigations, this article is focused to investigate P(GR4)P\left(GR_{4}\right) GRW spacetime satisfying f()f(\mathcal{R}) gravity. The f()f(\mathcal{R}) gravity model, f()=μctanh(c)f\left(\mathcal{R}\right)=\mathcal{R}-\mu\mathcal{R}_{c}\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}}) [40] satisfies local gravity constraint [41]. In this case, f>0f_{\mathcal{R}\mathcal{R}}>0 if >0\mathcal{R}>0 and f>0f_{\mathcal{R}}>0 if μ<1\mu<1. The condition f>0f_{\mathcal{R}\mathcal{R}}>0 is required for the consistency of local gravity tests [42], for the presence of the matter-dominated epoch [43] and for the stability of cosmological perturbations [44]. In this paper, we choose the model f()=μctanh(c)f\left(\mathcal{R}\right)=\mathcal{R}-\mu\mathcal{R}_{c}\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}}) in which μ\mu and c\mathcal{R}_{c} are positive constants, which is chosen to explain different ECs. The level of theoretical variety introduced by these higher curvature components is not found in more basic theories like General Relativity. This variability may open up new avenues for understanding and predicting late-time cosmic acceleration. Higher curvature components in gravity theories may result in changed dynamics at late-times, which may enable more realistic modelling of the cosmic acceleration at late-times and providing a better fit to observational data. The inclusion of higher curvature terms may be inspired by their importance in the effective action of quantum gravity theories, like string theory [45]. Examining their consequences for late-time cosmology can provide insight into the interaction between quantum physics and gravity. Investigating late-time cosmology in these structures offers a special chance to test basic physics outside of general relativity. Higher curvature terms may affect the behaviour of cosmic structures or affect cosmological perturbations, among other particular cosmological effects. Studying these implications can provide important insights into the fundamentals of gravity theory. These ideas might provide an alternate explanation for cosmic acceleration to dark energy, which could help to solve some of the long-standing cosmological problems. Higher curvature terms can have an impact on gravitational wave generation and propagation. Examining late-time cosmology in these theories may have consequences for gravitational wave signal detection and interpretation.

This article is structured as: In Section 2, an example of a P(GR4)P\left(GR_{4}\right) spacetime is illustrated. P(GR4)P\left(GR_{4}\right) GRW spacetimes and P(GR4)P\left(GR_{4}\right) PFS are investigated in next two Sections. Finally, a P(GR4)P\left(GR_{4}\right) GRW spacetime in f()f(\mathcal{R}) gravity theory is considered.

2. Example of a P(GR4)P\left(GR_{4}\right) spacetime

Choose a Lorentzian metric gg on 4\mathbb{R}^{4} described by

ds2=gijdyidyj=(dy1)2+(y1)2(dy2)2+(y2)2(dy3)2(dy4)2,ds^{2}=g_{ij}dy^{i}dy^{j}=\left(dy^{1}\right)^{2}+\left(y^{1}\right)^{2}\left(dy^{2}\right)^{2}+\left(y^{2}\right)^{2}\left(dy^{3}\right)^{2}-\left(dy^{4}\right)^{2}, (2.1)

where i,j=1,2,3,4i,j=1,2,3,4. Using (2.1), we observe that the non-vanishing components of the metric tensor are

g11=1,g22=(y1)2,g33=(y2)2,g44=1g_{11}=1,\quad g_{22}=\left(y^{1}\right)^{2},\quad g_{33}=\left(y^{2}\right)^{2},\quad g_{44}=-1 (2.2)

and the associated contravariant components are

g11=1,g22=1(y1)2,g33=1(y2)2,g44=1.g^{11}=1,\quad g^{22}=\dfrac{1}{\left(y^{1}\right)^{2}}\,,\quad g^{33}=\dfrac{1}{\left(y^{2}\right)^{2}}\,,\quad g^{44}=-1. (2.3)

Utilizing equations (2.2) and (2.3), here we determine the components (non-vanishing) of the Christoffel symbols, the curvature tensor and the Ricci tensor and they are

Γ221=y1,Γ332=y2(y1)2,Γ122=1y1,Γ233=1y2,\Gamma_{22}^{1}=-y^{1},\quad\Gamma_{33}^{2}=-\dfrac{y^{2}}{\left(y^{1}\right)^{2}}\,,\quad\Gamma_{12}^{2}=\dfrac{1}{y^{1}}\,,\quad\Gamma_{23}^{3}=\dfrac{1}{y^{2}}\,, (2.4)
1332=y2y1,12=1y1y2\mathcal{R}_{1332}=-\dfrac{y^{2}}{y^{1}}\,,\quad\mathcal{R}_{12}=-\dfrac{1}{y^{1}y^{2}} (2.5)

and the symmetric properties lead to the other components.

The covariant derivative of non-vanishing Ricci tensor is written as

12,1=1y2(y1)2and12,2=1y1(y2)2.\mathcal{R}_{12,1}=\dfrac{1}{y^{2}\left(y^{1}\right)^{2}}\quad\mathrm{and}\quad\mathcal{R}_{12,2}=\dfrac{1}{y^{1}\left(y^{2}\right)^{2}}\,. (2.6)

The one-forms we select are as follows:

ui(y)={1, when i=40, otherwise, u_{i}\left(y\right)=\begin{cases}1,&$ when $i=4\\ 0,&$ otherwise, $\end{cases} (2.7)
vi(y)={y1, when i=2y2, when i=3 0, otherwise v_{i}\left(y\right)=\begin{cases}y^{1},&$ when $i=2\\ y^{2},&$ when $i=3\\ \,0,&$ otherwise $\end{cases} (2.8)

and

wi(y)={y2, when i=1y1, when i=2 0, otherwise w_{i}\left(y\right)=\begin{cases}y^{2},&$ when $i=1\\ y^{1},&$ when $i=2\\ \,0,&$ otherwise $\end{cases} (2.9)

for all y4y\in\mathbb{R}^{4}.

We consider 𝒟ij\mathcal{D}_{ij} as follows:

𝒟ij(y)={1(y1)2(y2)2, when i=1,j=20, otherwise \mathcal{D}_{ij}\left(y\right)=\begin{cases}\dfrac{1}{\left(y^{1}\right)^{2}\left(y^{2}\right)^{2}}\,,&$ when $i=1,\,j=2\\ \quad\quad 0,&$ otherwise $\end{cases} (2.10)

for all y4y\in\mathbb{R}^{4}. It is enough to examine at the subsequent equations in order to confirm the relation (1.7):

12,1=u112+v1g12+w1𝒟12\mathcal{R}_{12,1}=u_{1}\mathcal{R}_{12}+v_{1}g_{12}+w_{1}\mathcal{D}_{12} (2.11)

and

12,2=u212+v2g12+w2𝒟12.\mathcal{R}_{12,2}=u_{2}\mathcal{R}_{12}+v_{2}g_{12}+w_{2}\mathcal{D}_{12}. (2.12)

The other cases hold trivially.

Now, right hand side of (2.11) =u112+v1g12+w1𝒟12\displaystyle=u_{1}\mathcal{R}_{12}+v_{1}g_{12}+w_{1}\mathcal{D}_{12}
=0(1y1y2)+0+y21(y1)2(y2)2\displaystyle=0\cdot\left(\dfrac{-1}{y^{1}y^{2}}\right)+0+y^{2}\cdot\dfrac{1}{\left(y^{1}\right)^{2}\left(y^{2}\right)^{2}}
=1y2(y1)2=12,1.\displaystyle=\dfrac{1}{y^{2}\left(y^{1}\right)^{2}}=\mathcal{R}_{12,1}.

By applying the same deduction, it is possible to demonstrate that equation (2.12) is likewise true.

By virtue of (2.3) and (2.7) we find

gijuiuj=1,andui(y)=gijuj(y)={1, when i=4  0, otherwise g^{ij}u_{i}u_{j}=-1,\quad\mathrm{and}\quad u^{i}\left(y\right)=g^{ij}u_{j}\left(y\right)=\begin{cases}-1,&$ when $i=4\\ \;\;0,&$ otherwise $\end{cases} (2.13)

Equations (2.10) and (2.13) together yield

𝒟ijui=0.\mathcal{D}_{ij}u^{i}=0. (2.14)

Clearly, the trace(𝒟ij)=0.\left(\mathcal{D}_{ij}\right)=0. Therefore, (4,g)\left(\mathbb{R}^{4},g\right) is a P(GR4)P\left(GR_{4}\right) spacetime.

3. P(GR4)P\left(GR_{4}\right) GRW spacetimes

Theorem A. [6] A M4\mathrm{M}^{4} is a GRW spacetime if and only if it admits a unit torse-forming time-like vector uiu_{i}:

kui=φ{gki+ukui}\nabla_{k}u_{i}=\varphi\left\{g_{ki}+u_{k}u_{i}\right\} (3.1)

and uiu_{i} is an eigenvector of ij\mathcal{R}_{ij}, that is,

ijuj=θuj\mathcal{R}_{ij}u^{j}=\theta u_{j} (3.2)

in which φ\varphi and θ\theta are non-zero scalars.

Multiplying (1.7) with uiu^{i} and using (1.8) infers

(kij)ui=ukijui+vkuj.\left(\nabla_{k}\mathcal{R}_{ij}\right)u^{i}=u_{k}\mathcal{R}_{ij}u^{i}+v_{k}u_{j}. (3.3)

Since k(ijui)=(kij)ui+ij(kui)\nabla_{k}\left(\mathcal{R}_{ij}u^{i}\right)=\left(\nabla_{k}\mathcal{R}_{ij}\right)u^{i}+\mathcal{R}_{ij}\left(\nabla_{k}u^{i}\right), (3.3) becomes

k(ijui)ij(kui)=ukijui+vkuj.\nabla_{k}\left(\mathcal{R}_{ij}u^{i}\right)-\mathcal{R}_{ij}\left(\nabla_{k}u^{i}\right)=u_{k}\mathcal{R}_{ij}u^{i}+v_{k}u_{j}. (3.4)

Equations (3.1), (3.2) and (3.4) together imply

φjk=θφgjk+θkujθukujvkuj.\varphi\mathcal{R}_{jk}=\theta\varphi g_{jk}+\theta_{k}u_{j}-\theta u_{k}u_{j}-v_{k}u_{j}. (3.5)

Interchanging jj and kk in (3.5), we obtain

φkj=θφgkj+θjukθujukvjuk.\varphi\mathcal{R}_{kj}=\theta\varphi g_{kj}+\theta_{j}u_{k}-\theta u_{j}u_{k}-v_{j}u_{k}. (3.6)

Subtracting (3.6) from (3.5), we have

θkuj=θjuk+vkujvjuk.\theta_{k}u_{j}=\theta_{j}u_{k}+v_{k}u_{j}-v_{j}u_{k}. (3.7)

Multiplying (3.7) with uju^{j}, we acquire

θk=(vjuj)uk(θjuj)uk+vk,\theta_{k}=\left(v_{j}u^{j}\right)u_{k}-\left(\theta_{j}u^{j}\right)u_{k}+v_{k}, (3.8)

that is,

θk=(f1f3)uk+vk,wheref1=vjujandf3=θjuj.\theta_{k}=\left(f_{1}-f_{3}\right)u_{k}+v_{k},\quad\mathrm{where}\quad f_{1}=v_{j}u^{j}\quad\mathrm{and}\quad f_{3}=\theta_{j}u^{j}. (3.9)

From (3.5) and (3.9), it follows that

jk=θgjk+(f1f3θφ)ujuk.\mathcal{R}_{jk}=\theta g_{jk}+(\dfrac{f_{1}-f_{3}-\theta}{\varphi})u_{j}u_{k}. (3.10)

Hence, we reach:

Theorem 1.

A P(GR4)P\left(GR_{4}\right) GRW\mathrm{GRW} spacetime represents a PFS\mathrm{PFS}.

In light of equations (1.1), (1.4) and (3.10), we acquire

κ(σp2)=θ\kappa\left(\dfrac{\sigma-p}{2}\right)=\theta (3.11)

and

κ(σ+p)=f1f3θφ.\kappa\left(\sigma+p\right)=\dfrac{f_{1}-f_{3}-\theta}{\varphi}\,. (3.12)

Equations (3.11) and (3.12) together give

pσ=f1f3θ2θφf1f3θ+2θφ.\dfrac{p}{\sigma}=\dfrac{f_{1}-f_{3}-\theta-2\theta\varphi}{f_{1}-f_{3}-\theta+2\theta\varphi}\,. (3.13)

We observe that (3.13) implies p=0p=0 for f1=f3+θ(1+2φ)f_{1}=f_{3}+\theta\left(1+2\varphi\right), σ=3p\sigma=3p for f1=f3+θ(1+4φ)f_{1}=f_{3}+\theta\left(1+4\varphi\right) and σ+p=0\sigma+p=0 for f1=f3+θf_{1}=f_{3}+\theta, respectively. Hence, we obtain:

Corollary 1.

A P(GR4)P\left(GR_{4}\right) GRW\mathrm{GRW} spacetime represents a

  1. (1)

    state equation of the form (3.13),

  2. (2)

    dust matter fluid for f1=f3+θ(1+2φ)f_{1}=f_{3}+\theta\left(1+2\varphi\right),

  3. (3)

    radiation era for f1=f3+θ(1+4φ)f_{1}=f_{3}+\theta\left(1+4\varphi\right) and

  4. (4)

    dark energy epoch of the Universe for f1=f3+θf_{1}=f_{3}+\theta.

Definition 1.

A vector field uku_{k} is Riemann compatible [46] if and only if it is Weyl compatible and

uijhuhujihuh=0.\displaystyle u_{i}\mathcal{R}^{h}_{j}u_{h}-u_{j}\mathcal{R}^{h}_{i}u_{h}=0. (3.14)

In a P(GR4)P\left(GR_{4}\right) GRW spacetime

ij=αgij+βuiuj\displaystyle\mathcal{R}_{ij}={\alpha}g_{ij}+{\beta}u_{i}u_{j} (3.15)

in which

α=θ,β=f1f3θφ.\displaystyle\alpha=\theta,~~\beta=\frac{f_{1}-f_{3}-\theta}{\varphi}. (3.16)

Multiplying equation (3.15) by gilg^{il}, we acquire

uijhuhujihuh=ui(αδjh+βuhuj)uh\displaystyle u_{i}\mathcal{R}^{h}_{j}u_{h}-u_{j}\mathcal{R}^{h}_{i}u_{h}=u_{i}(\alpha\delta^{h}_{j}+\beta u^{h}u_{j})u_{h}
uj(αδih+βuhui)uh\displaystyle-u_{j}(\alpha\delta^{h}_{i}+\beta u^{h}u_{i})u_{h}
=0.\displaystyle=0. (3.17)

Thus the equation (3.14) is verified. Hence, we write:

Theorem 2.

Every vector field of a P(GR4)P\left(GR_{4}\right) GRW spacetime is Riemann compatible.

The magnetic and electric parts of the Weyl tensor are given by

Hij=ukulC~kijl,Eij=ukulCkijl,\displaystyle H_{ij}=u^{k}u^{l}\tilde{C}_{kijl},\;\;\;\;E_{ij}=u^{k}u^{l}C_{kijl}, (3.18)

in which ukuk=1u_{k}u^{k}=-1 and C~kijl=12ϵkilmCjllm\tilde{C}_{kijl}=\frac{1}{2}{\epsilon}_{kilm}C^{lm}_{jl} is the dual [47]. These two tensors are traceless, symmetric and obey Ehkuh=0E_{hk}u^{h}=0 and Hhkuh=0H_{hk}u^{h}=0.

For n=4n=4, from the definition of Weyl tensor, we provide

Eij=hijkuhuk12(gijhkgikhj+ghkijghjik)uhuk\displaystyle E_{ij}=\mathcal{R}_{hijk}u^{h}u^{k}-\frac{1}{2}(g_{ij}\mathcal{R}_{hk}-g_{ik}\mathcal{R}_{hj}+g_{hk}\mathcal{R}_{ij}-g_{hj}\mathcal{R}_{ik})u^{h}u^{k}
+6(gijghkgikghj).\displaystyle+\frac{\mathcal{R}}{6}(g_{ij}g_{hk}-g_{ik}g_{hj}). (3.19)

By hypothesis P(GR4)P\left(GR_{4}\right) is a GRW spacetime.

Hence, using equation (3.15) in (3.19), we infer

Eij=hijkuhuk+αβ3(gij+uiuj).\displaystyle E_{ij}=\mathcal{R}_{hijk}u^{h}u^{k}+\frac{\alpha-\beta}{3}(g_{ij}+u_{i}u_{j}). (3.20)

From equation (3.1), we obtain

kjui\displaystyle\nabla_{k}\nabla_{j}u_{i} =\displaystyle= φk{gij+uiuj}\displaystyle\varphi_{k}\left\{g_{ij}+u_{i}u_{j}\right\} (3.21)
φ[φ{gik+uiuk}uj+φ{gjk+ujuk}ui]\displaystyle\varphi[\varphi\left\{g_{ik}+u_{i}u_{k}\right\}u_{j}+\varphi\left\{g_{jk}+u_{j}u_{k}\right\}u_{i}]

and hence we get

kjuijkui=φk{gij+uiuj}φj{gik+uiuk}.\nabla_{k}\nabla_{j}u_{i}-\nabla_{j}\nabla_{k}u_{i}=\varphi_{k}\left\{g_{ij}+u_{i}u_{j}\right\}-\varphi_{j}\left\{g_{ik}+u_{i}u_{k}\right\}. (3.22)

From the above, we can easily acquire

uhihjk=φk{gij+uiuj}φj{gik+uiuk}u_{h}\mathcal{R}^{h}_{i}jk=\varphi_{k}\left\{g_{ij}+u_{i}u_{j}\right\}-\varphi_{j}\left\{g_{ik}+u_{i}u_{k}\right\}

which entails

uhhijk=φk{gij+uiuj}φj{gik+uiuk}u^{h}\mathcal{R}_{h}ijk=\varphi_{k}\left\{g_{ij}+u_{i}u_{j}\right\}-\varphi_{j}\left\{g_{ik}+u_{i}u_{k}\right\} (3.23)

Multiplying the foregoing equation by uku^{k} yields

ukuhhijk=f{gij+uiuj},u^{k}u^{h}\mathcal{R}_{h}ijk=f\left\{g_{ij}+u_{i}u_{j}\right\}, (3.24)

in which f=ukφkf=u^{k}\varphi_{k}.

Using equation (3.24) in (3.20), we get

Eij={f+αβ3}(gij+uiuj).\displaystyle E_{ij}=\{f+\frac{\alpha-\beta}{3}\}(g_{ij}+u_{i}u_{j}). (3.25)

Now, from equation (3.24), we provide

ukuhhk=3f.u^{k}u^{h}\mathcal{R}_{h}k=3f. (3.26)

Also, from equation (3.15), we acquire

ukuhhk=αβ3.u^{k}u^{h}\mathcal{R}_{h}k=-\frac{\alpha-\beta}{3}. (3.27)

Hence, the last two equations jointly give

f=αβ3.f=-\frac{\alpha-\beta}{3}. (3.28)

Using equation (3.28) in (3.25), we get Eij=0E_{ij}=0.

Thus we can write

Theorem 3.

In a P(GR4)P\left(GR_{4}\right) GRW spacetime, the electric part of the spacetime vanishes.

Let P(GR4)P\left(GR_{4}\right) spacetime be a GRW spacetime. From the last two Theorems it follows that the velocity vector is Riemann compatible and the electric part of the spacetime under consideration vanishes. Hence, using Definition 1, we conclude that the velocity vector is Weyl compatible.

Also, we are aware that H=0H=0, or the magnetic component of the Weyl tensor disappears when a spacetime’s vector field is Weyl compatible [46]. Therefore, the spacetime under consideration becomes conformally flat since E=0E=0 and H=0H=0 and hence is of petrov type O. Therefore we state:

Theorem 4.

A P(GR4)P\left(GR_{4}\right) GRW spacetime is conformally flat and is of petrov type O.

4. P(GR4)P\left(GR_{4}\right) perfectfluidspacetimes\mathrm{perfect\,fluid\,spacetimes}

Considering the covariant derivative of (1.1), we infer that

kij=αkgij+βkuiuj+β{ui(kuj)+uj(kui)}.\nabla_{k}\mathcal{R}_{ij}=\alpha_{k}g_{ij}+\beta_{k}u_{i}u_{j}+\beta\left\{u_{i}\left(\nabla_{k}u_{j}\right)+u_{j}\left(\nabla_{k}u_{i}\right)\right\}. (4.1)

Using (1.1) and (4.1) in (1.7), we find

β{ui(kuj)+uj(kui)}=(αuk+vkαk)gij+(βukβk)uiuj+wk𝒟ij.\beta\left\{u_{i}\left(\nabla_{k}u_{j}\right)+u_{j}\left(\nabla_{k}u_{i}\right)\right\}=\left(\alpha u_{k}+v_{k}-\alpha_{k}\right)g_{ij}+\left(\beta u_{k}-\beta_{k}\right)u_{i}u_{j}+w_{k}\mathcal{D}_{ij}. (4.2)

Since ujuj=1u_{j}u^{j}=-1, therefore uj(kuj)+uj(kuj)=0u^{j}(\nabla_{k}u_{j})+u_{j}(\nabla_{k}u^{j})=0, which implies

uj(kuj)=0.u^{j}\left(\nabla_{k}u_{j}\right)=0. (4.3)

Multiplying (4.2) with uiu^{i} and using (4.3), we notice that

β(kuj)=pkuj,wherepk=βukβkαukvk+αk.\beta\left(\nabla_{k}u_{j}\right)=p_{k}u_{j},\quad\mathrm{where}\quad p_{k}=\beta u_{k}-\beta_{k}-\alpha u_{k}-v_{k}+\alpha_{k}. (4.4)

Again, multiplying (4.4) with uju^{j} and using (4.3), we acquire

pk=0.p_{k}=0. (4.5)

Equations (4.4) and (4.5) reflect that

β(kuj)=0,\beta\left(\nabla_{k}u_{j}\right)=0, (4.6)

which entails the following cases:

Case 1. For β=0\beta=0, equation (1.4) gives us σ+p=0\sigma+p=0. Thus, the spacetime represents a dark energy epoch of the Universe.

Case 2. For β\beta0\neq 0, (4.6) becomes

kuj=0.\nabla_{k}u_{j}=0. (4.7)

Contracting with gkjg^{kj} yields

kuk=0.\nabla_{k}u^{k}=0. (4.8)

Again multiplying equation (4.7) by uku^{k} gives

uj˙=ukkuj=0.\dot{u_{j}}=u^{k}\nabla_{k}u_{j}=0. (4.9)

The equations (4.7)-(4.9) show that the vector field uku_{k} is parallel, conservative, and acceleration-free.

The vorticity tensor μlk\mu_{lk} is described by

μlk=12(lakkal)+12(alak˙akal˙).\mu_{lk}=\frac{1}{2}(\nabla_{l}a_{k}-\nabla_{k}a_{l})+\frac{1}{2}(a_{l}\dot{a_{k}}-a_{k}\dot{a_{l}}). (4.10)

Therefore, using equations (4.7) and (4.9), we provide μlk=0\mu_{lk}=0, that is, vorticity free.

The covariant gradient kuj\nabla_{k}u_{j} in a spacetime has the following standard decomposition [8]:

kuj=iuin1(gkj+ukuj)ukuj˙+μkj+σkj.\nabla_{k}u_{j}=\frac{\nabla_{i}u^{i}}{n-1}(g_{kj}+u_{k}u_{j})-u_{k}\dot{u_{j}}+\mu_{kj}+\sigma_{kj}. (4.11)

In view of the above equations, the last equation provides σkj=0\sigma_{kj}=0, that is, shear free.

A spacetime is said to be stationary if uiu_{i} is Killing and static ([48], [49], p. 283) for irrotational vector uiu_{i}. For a smooth vector vv,

£vglk=lvk+kvl,\pounds_{v}g_{lk}=\nabla_{l}v_{k}+\nabla_{k}v_{l},

£\pounds denotes the Lie derivative. As kui=0\nabla_{k}u_{i}=0, then £ugki=0\pounds_{u}g_{ki}=0, means that uiu_{i} is Killing. Moreover, kui=0\nabla_{k}u_{i}=0 gives uiu_{i} is irrotational. Therefore, the spacetime is static.

Hence, it is possible to declare the following:

Theorem 5.

A P(GR4)P\left(GR_{4}\right) PFS\mathrm{PFS} represents either a dark energy epoch of the Universe or, the vector field uku_{k} is parallel, conservative, acceleration-free, vorticity-free, and shear-free and becomes a static spacetime.

5. P(GR4)P\left(GR_{4}\right) GRW spacetime in f()f(\mathcal{R}) gravity

Here, we will look into a few characteristics of a P(GR4)P\left(GR_{4}\right) GRW spacetime in f()f(\mathcal{R}) gravity.

Choose the Einstein-Hilbert action

S=12k2d4xgf()+d4xgLm\displaystyle S=\frac{1}{2k^{2}}\int{d^{4}x{\sqrt{-g}}f(\mathcal{R})}+\int{d^{4}x{\sqrt{-g}}L_{m}}

in which LmL_{m} is the matter Lagrangian density described as

Tij=2gδ(gLm)δgij\displaystyle T_{ij}=-\frac{2}{\sqrt{-g}}\frac{\delta(\sqrt{-g}L_{m})}{\delta g^{ij}}

In this case, κ2=8πG\kappa^{2}=8{\pi}G, GG stands for the Newton’s constant and hence the modified formula can be written as [50]

S=12k2d4xgf().\displaystyle S=\frac{1}{2k^{2}}\int{d^{4}x{\sqrt{-g}}f(\mathcal{R})}.

The field equations can be expressed in the following form by applying the variation with gijg^{ij}

f()ijf()2gij+(gijij)f()=k2Tij\displaystyle f_{\mathcal{R}}(\mathcal{R})\mathcal{R}_{ij}-\frac{f(\mathcal{R})}{2}g_{ij}+(g_{ij}{\square}-\nabla_{i}\nabla_{j})f_{\mathcal{R}}(\mathcal{R})=k^{2}T_{ij} (5.1)

where \square indicates the D’Alembert’s operator.

Taking trace of equation (5.1) provides

3f()+f()2f()=k2T.\displaystyle 3{\square}f_{\mathcal{R}}(\mathcal{R})+\mathcal{R}f_{\mathcal{R}}(\mathcal{R})-2f(\mathcal{R})=k^{2}T. (5.2)

Subtracting the expression f()2gij\frac{\mathcal{R}f_{\mathcal{R}}(\mathcal{R})}{2}g_{ij} from (5.1), we infer

f()ijf()2gij=k2Tij+k2Tij(curve)\displaystyle f_{\mathcal{R}}(\mathcal{R})\mathcal{R}_{ij}-\frac{\mathcal{R}f_{\mathcal{R}}(\mathcal{R})}{2}g_{ij}=k^{2}T_{ij}+k^{2}T^{(curve)}_{ij} (5.3)

such that

Tij(eff)=Tij+Tij(curve)\displaystyle T^{(eff)}_{ij}=T_{ij}+T^{(curve)}_{ij} (5.4)

in which

Tij(curve)=1k2[(f()f())2gij+(ijgij)f()].\displaystyle T^{(curve)}_{ij}=\frac{1}{k^{2}}\big[\frac{(f(\mathcal{R})-\mathcal{R}f_{\mathcal{R}}(\mathcal{R}))}{2}g_{ij}+(\nabla_{i}\nabla_{j}-g_{ij}{\square})f_{\mathcal{R}}(\mathcal{R})\big]. (5.5)

Using the equation (5.3), we acquire

ij2gij=k2f()Tij(eff)\displaystyle\mathcal{R}_{ij}-\frac{\mathcal{R}}{2}g_{ij}=\frac{k^{2}}{f_{\mathcal{R}}(\mathcal{R})}T_{ij}^{(eff)} (5.6)

which obeys the equations (5.4) and (5.5).

Multiplying equation (3.10) with gjkg^{jk}, we reach

f1f3θφ=4θ.\dfrac{f_{1}-f_{3}-\theta}{\varphi}=4\theta-\mathcal{R}. (5.7)

Utilizing (5.7) in (3.10), we get

jk=θgjk+(4θ)ujuk.\mathcal{R}_{jk}=\theta g_{jk}+\left(4\theta-\mathcal{R}\right)u_{j}u_{k}. (5.8)

Taking covariant derivative of the foregoing equation yields

ijk\displaystyle\nabla_{i}\mathcal{R}_{jk} =\displaystyle= θigjk+(4θii)ujuk\displaystyle\theta_{i}g_{jk}+\left(4\theta_{i}-\nabla_{i}\mathcal{R}\right)u_{j}u_{k} (5.9)
+(4θ)[ukiuj+ujiuk].\displaystyle+\left(4\theta-\mathcal{R}\right)[u_{k}\nabla_{i}u_{j}+u_{j}\nabla_{i}u_{k}].

Multiplying the above equation by uku^{k} provides

(4θ)iuj=(i3θi)ujukijk.\left(4\theta-\mathcal{R}\right)\nabla_{i}u_{j}=(\nabla_{i}\mathcal{R}-3\theta_{i})u_{j}-u^{k}\nabla_{i}\mathcal{R}_{jk}. (5.10)

Again multiplying the previous equation by uiu^{i} gives

(4θ)uj˙=(˙3f3+θf1)uj,\left(4\theta-\mathcal{R}\right)\dot{u_{j}}=(\dot{\mathcal{R}}-3f_{3}+\theta-f_{1})u_{j}, (5.11)

where uiiuj=uj˙u^{i}\nabla_{i}u_{j}=\dot{u_{j}} and uii=˙u^{i}\nabla_{i}\mathcal{R}=\dot{\mathcal{R}}.

Using ujuj˙=0u_{j}\dot{u_{j}}=0, we readily conclude that the acceleration vector

uj˙=0.\dot{u_{j}}=0. (5.12)

Now transvecting equation (1.9) by uiu^{i}, we acquire

˙+=4f1.\dot{\mathcal{R}}+\mathcal{R}=4f_{1}. (5.13)

Covariant derivative of equation (1.9) produces

ji=ui(j)+(jui)+4jvi.\nabla_{j}\nabla_{i}\mathcal{R}=u_{i}(\nabla_{j}\mathcal{R})+\mathcal{R}(\nabla_{j}u_{i})+4\nabla_{j}v_{i}. (5.14)

Interchanging ii and jj yields

ij=uj(i)+(iuj)+4ivj.\nabla_{i}\nabla_{j}\mathcal{R}=u_{j}(\nabla_{i}\mathcal{R})+\mathcal{R}(\nabla_{i}u_{j})+4\nabla_{i}v_{j}. (5.15)

Utilizing the last two equations, we get

[iujjui]+[4viuj4uivj]+4[ivjjvi]\mathcal{R}[\nabla_{i}u_{j}-\nabla_{j}u_{i}]+[4v_{i}u_{j}-4u_{i}v_{j}]+4[\nabla_{i}v_{j}-\nabla_{j}v_{i}] (5.16)

Transvecting the last equation by uiu^{i}, we obtain

vj=f1uj,v_{j}=-f_{1}u_{j}, (5.17)

since ui˙=0\dot{u_{i}}=0.

Suppose viv_{i} is parallel. Then using equation (5.17) in equation (5.16), we provide

ij=(4f1)uiuj.\nabla_{i}\nabla_{j}\mathcal{R}=(\mathcal{R}-4f_{1})u_{i}u_{j}. (5.18)

Since f()f(\mathcal{R}) is an analytic function, then we write

ijf()=f()ij+f()(i)(j).\nabla_{i}\nabla_{j}f_{\mathcal{R}}(\mathcal{R})=f_{\mathcal{R}\mathcal{R}}(\mathcal{R})\nabla_{i}\nabla_{j}\mathcal{R}+f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})(\nabla_{i}\mathcal{R})(\nabla_{j}\mathcal{R}). (5.19)

Hence, equation (5.19) yields

f()=f()+f()gij(i)(j).{\square}f_{\mathcal{R}}(\mathcal{R})=f_{\mathcal{R}\mathcal{R}}(\mathcal{R}){\square}\mathcal{R}+f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})g^{ij}(\nabla_{i}\mathcal{R})(\nabla_{j}\mathcal{R}). (5.20)

Now, using the equations (1.9), (5.18) and (5.19), we acquire

ijf()\displaystyle\nabla_{i}\nabla_{j}f_{\mathcal{R}}(\mathcal{R}) =\displaystyle= {(4f1)f()+(28f1)f()}uiuj\displaystyle\{(\mathcal{R}-4f_{1})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})+(\mathcal{R}^{2}-8\mathcal{R}f_{1})f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})\}u_{i}u_{j} (5.21)
+16f()vivj.\displaystyle+16f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})v_{i}v_{j}.

Multiplying equation (5.21) by gijg^{ij} and making use of the equations (5.18) and (5.20), we get

f()=(4f1)f()(8f12+16)f().{\square}f_{\mathcal{R}}(\mathcal{R})=(4f_{1}-\mathcal{R})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})-(8\mathcal{R}f_{1}-\mathcal{R}^{2}+16)f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R}). (5.22)

Again, utilizing the equations (5.4), (5.5) and (5.6), we provide

f()(ij2gij)=k2Tij+f()f()2gij+(ijgij)f().\displaystyle f_{\mathcal{R}}(\mathcal{R})(\mathcal{R}_{ij}-\frac{\mathcal{R}}{2}g_{ij})=k^{2}T_{ij}+\frac{f(\mathcal{R})-\mathcal{R}f_{\mathcal{R}}(\mathcal{R})}{2}g_{ij}+(\nabla_{i}\nabla_{j}-g_{ij}{\square})f_{\mathcal{R}}(\mathcal{R}). (5.23)

Substituting the equations (5.21) and (5.22) in equation (5.23) infer

f()(ij2gij)=κ2Tij\displaystyle f_{\mathcal{R}}(\mathcal{R})\big(\mathcal{R}_{ij}-\frac{\mathcal{R}}{2}g_{ij}\big)=\kappa^{2}T_{ij}
+[f()f()2(8f12+16)f()\displaystyle+\big[\frac{f(\mathcal{R})-\mathcal{R}f_{\mathcal{R}}(\mathcal{R})}{2}-(8\mathcal{R}f_{1}-\mathcal{R}^{2}+16)f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})
(4f1)f()]gij\displaystyle-(4f_{1}-\mathcal{R})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})\big]g_{ij}
[(8f12)f()+(4f1)f()]uiuj.\displaystyle-\big[-(8\mathcal{R}f_{1}-\mathcal{R}^{2})f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})+(4f_{1}-\mathcal{R})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})\big]u_{i}u_{j}. (5.24)

Using equation (5.8) in equation (5.24), we acquire

Tij\displaystyle T_{ij} =\displaystyle= 1κ2[f()2θf()2(8f12+16)f()\displaystyle-\frac{1}{\kappa^{2}}\big[\frac{f(\mathcal{R})-2\theta f_{\mathcal{R}}(\mathcal{R})}{2}-(8\mathcal{R}f_{1}-\mathcal{R}^{2}+16)f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R}) (5.25)
(4f1)f()]gij\displaystyle-(4f_{1}-\mathcal{R})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})\big]g_{ij}
+1κ2[(8f12)f()\displaystyle+\frac{1}{\kappa^{2}}\big[-(8\mathcal{R}f_{1}-\mathcal{R}^{2})f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})
+(4f1)f()(4θ)f]uiuj.\displaystyle+(4f_{1}-\mathcal{R})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})-(4\theta-\mathcal{R})f_{\mathcal{R}}\mathcal{R}\big]u_{i}u_{j}.

Therefore, using equations (1.2) and (5.25), we find

p\displaystyle p =\displaystyle= 1κ2[f()2θf()2(8f12+16)f()\displaystyle-\frac{1}{\kappa^{2}}\big[\frac{f(\mathcal{R})-2\theta f_{\mathcal{R}}(\mathcal{R})}{2}-(8\mathcal{R}f_{1}-\mathcal{R}^{2}+16)f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R}) (5.26)
(4f1)f()]\displaystyle-(4f_{1}-\mathcal{R})f_{\mathcal{R}\mathcal{R}}(\mathcal{R})\big]

and

σ\displaystyle\sigma =\displaystyle= 1κ2[f()(10θ2)f()2\displaystyle\frac{1}{\kappa^{2}}\big[\frac{f(\mathcal{R})-(10\theta-2\mathcal{R})f_{\mathcal{R}}(\mathcal{R})}{2} (5.27)
(16f122+16)f()].\displaystyle-(16\mathcal{R}f_{1}-2\mathcal{R}^{2}+16)f_{\mathcal{R}\mathcal{R}\mathcal{R}}(\mathcal{R})\big].

Therefore, we provide:

Theorem 6.

In a P(GR4)P\left(GR_{4}\right) GRW spacetime satisfying f()f(\mathcal{R}) gravity, pp and σ\sigma are described by (5.26) and (5.27), respectively.

If f()=f(\mathcal{R})=\mathcal{R}, then f()f(\mathcal{R}) theory reduces to Einstein theory. In this case the equations (5.26) and (5.27) takes the form

p\displaystyle p =\displaystyle= 1κ2[2θ2]\displaystyle-\frac{1}{\kappa^{2}}[\frac{\mathcal{R}-2\theta}{2}] (5.28)

and

σ\displaystyle\sigma =\displaystyle= 1κ2[310θ2].\displaystyle\frac{1}{\kappa^{2}}[\frac{3\mathcal{R}-10\theta}{2}]. (5.29)

Thus, we state:

Corollary 2.

In a P(GR4)P\left(GR_{4}\right) GRW spacetime satisfying f()f(\mathcal{R}) gravity with the condition f()=f(\mathcal{R})=\mathcal{R}, pp and σ\sigma are given by (5.28) and (5.29), respectively.

Equations (5.28) and (5.29) jointly yield

pσ=2θ310θ.\dfrac{p}{\sigma}=-\dfrac{\mathcal{R}-2\theta}{3\mathcal{R}-10\theta}. (5.30)

We observe that (5.30) implies p=0p=0 for =2θ\mathcal{R}=2\theta, σ=3p\sigma=3p for =83θ\mathcal{R}=\frac{8}{3}\theta and σ+p=0\sigma+p=0 for =4θ\mathcal{R}=4\theta, respectively. Hence, we write:

Remark 1.

A P(GR4)P\left(GR_{4}\right) GRW spacetime satisfying f()f(\mathcal{R}) gravity with the condition f()=f(\mathcal{R})=\mathcal{R}, represents a

  1. (1)

    state equation of the form (5.30),

  2. (2)

    dust matter fluid for =2θ\mathcal{R}=2\theta,

  3. (3)

    radiation era for =83θ\mathcal{R}=\frac{8}{3}\theta and

  4. (4)

    dark energy epoch of the Universe for =4θ\mathcal{R}=4\theta.

5.1. Energy Conditions

In the following subsection, we verify the ECs for the f()f\left(\mathcal{R}\right)-gravity model f()=μctanh(c)f\left(\mathcal{R}\right)=\mathcal{R}-\mu\mathcal{R}_{c}\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}}) [40], in which μ\mu and c\mathcal{R}_{c} are positive constants.

In modified gravity the ECs are demonstrated as

NEC\displaystyle\mathrm{NEC} ifandonlyifσ+p0,\displaystyle\quad\mathrm{if\;and\;only\;if}\quad\sigma+p\geq 0,
TEC\displaystyle\mathrm{TEC} ifandonlyifσ3p0,\displaystyle\quad\mathrm{if\;and\;only\;if}\quad\sigma-3p\geq 0,
SEC\displaystyle\mathrm{SEC} ifandonlyifσ+p0andσ+3p0,\displaystyle\quad\mathrm{if\;and\;only\;if}\quad\sigma+p\geq 0\quad\mathrm{and}\quad\sigma+3p\geq 0,
DEC\displaystyle\mathrm{DEC} ifandonlyifσ±p0andσ0,\displaystyle\quad\mathrm{if\;and\;only\;if}\quad\sigma\pm p\geq 0\quad\mathrm{and}\quad\sigma\geq 0,
WEC\displaystyle\mathrm{WEC} ifandonlyifσ+p0andσ0.\displaystyle\quad\mathrm{if\;and\;only\;if}\quad\sigma+p\geq 0\quad\mathrm{and}\quad\sigma\geq 0.

With the help of (5.26) and (5.27), σ\sigma and pp are given by

κσ=\displaystyle\kappa\sigma= 212μctanh(c)(5θ)[1μ(1tanh(c)2)]\displaystyle\frac{\mathcal{R}}{2}-\frac{1}{2}\mu\mathcal{R}_{c}\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})-(5\theta-\mathcal{R})[1-\mu(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})]
(16f122+16)[2μ(1tanh(c)2)2c24μtanh(c)2(1tanh(c)2)c2],\displaystyle-(16\mathcal{R}f_{1}-2\mathcal{R}^{2}+16)\big[\frac{2\mu(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})^{2}}{\mathcal{R}_{c}^{2}}-\frac{4\mu\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2}(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})}{\mathcal{R}_{c}^{2}}\big], (5.31)
κp=\displaystyle\kappa p= 2+12μctanh(c)+θ[1μ(1tanh(c)2)]\displaystyle-\frac{\mathcal{R}}{2}+\frac{1}{2}\mu\mathcal{R}_{c}\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})+\theta[1-\mu(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})]
(8f12+16)[2μ(1tanh(c)2)2c24μtanh(c)2(1tanh(c)2)c2]\displaystyle-(8\mathcal{R}f_{1}-\mathcal{R}^{2}+16)\big[\frac{2\mu(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})^{2}}{\mathcal{R}_{c}^{2}}-\frac{4\mu\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2}(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})}{\mathcal{R}_{c}^{2}}\big]
2(4f1)μtanh(c)2(1tanh(c)2)c.\displaystyle-\frac{2(4f_{1}-\mathcal{R})\mu\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2}(1-\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}})^{2})}{\mathcal{R}_{c}}. (5.32)

The ECs for the above model are now examined. The ECs for this arrangement may now be discussed using equations (5.31) and (5.32).

[Uncaptioned image] [Uncaptioned image]
Fig. 1: Development of σ\sigma with reference to \mathcal{R} and θ\theta Fig. 2: Development of p+σp+\sigma with reference to \mathcal{R} and θ\theta
[Uncaptioned image] [Uncaptioned image]
Fig. 3: Development of σp\sigma-p with reference to \mathcal{R} and θ\theta Fig. 4: Development of σ+3p\sigma+3p with reference to \mathcal{R} and θ\theta
[Uncaptioned image] Fig. 5: Development of σ3p\sigma-3p with reference to \mathcal{R} and θ\theta

Figures 11 and 22 demonstrate that, for parameters θ,[1,3]\theta,\mathcal{R}\in\left[1,3\right], the energy density and p+σp+\sigma can not be negative, and that, for larger values of θ\theta and \mathcal{R}, they are high. Because NEC is a component of WEC, NEC and WEC are fulfilled. The σρ\sigma-\rho profile for θ,[1,3]\theta,\mathcal{R}\in\left[1,3\right] is positive, as seen in Fig. 33. It is evident from Figs. 11, 22, and 33 that DEC is validated. Furthermore, we can observe that SEC is satisfied from Figs. 22 and 44, and this finding yields the late-time acceleration of the Cosmos[51]. Moreover, each result aligns with the Λ\LambdaCDM model [52]. Fig. 55 shows that TEC is satisfied.

5.2.

The field equations of f()f(\mathcal{R}) gravity are as follows:

κTlk\displaystyle\kappa T_{lk} =\displaystyle= f()lkf′′′()lkf′′()lk\displaystyle f^{\prime}(\mathcal{R})\mathcal{R}_{lk}-f^{\prime\prime\prime}(\mathcal{R})\nabla_{l}\mathcal{R}\nabla_{k}\mathcal{R}-f^{\prime\prime}(\mathcal{R})\nabla_{l}\nabla_{k}\mathcal{R} (5.33)
+glk[f′′′()mm+f′′()212f()].\displaystyle+g_{lk}[f^{\prime\prime\prime}(\mathcal{R})\nabla_{m}\mathcal{R}\nabla^{m}\mathcal{R}+f^{\prime\prime}(\mathcal{R})\nabla^{2}\mathcal{R}-\frac{1}{2}f(\mathcal{R})].

For =\mathcal{R}= constant, we acquire

lkf2fglk=κfTlk.\mathcal{R}_{lk}-\frac{f}{2f^{\prime}}g_{lk}=\frac{\kappa}{f^{\prime}}T_{lk}. (5.34)

Using equation (1.2) in equation (5.34), we infer

lk=f2fglk+κf[(σ+p)ukul+pgkl].\mathcal{R}_{lk}=\frac{f}{2f^{\prime}}g_{lk}+\frac{\kappa}{f^{\prime}}[(\sigma+p)u_{k}u_{l}+pg_{kl}]. (5.35)

Making use of equations (3.10) and (5.35), we obtain

κf(σ+p)=f1f3θφ\frac{\kappa}{f^{\prime}}(\sigma+p)=\frac{f_{1}-f_{3}-\theta}{\varphi} (5.36)

and

pκf+f2f=θ.\frac{p\kappa}{f^{\prime}}+\frac{f}{2f^{\prime}}=\theta. (5.37)

Solving the last two equations, we get

p=fθκf2κp=\frac{f^{\prime}\theta}{\kappa}-\frac{f}{2\kappa} (5.38)

and

σ=f(f1f3θθφ)κφ+f2κ.\sigma=-\frac{f^{\prime}(f_{1}-f_{3}-\theta-\theta\varphi)}{\kappa\varphi}+\frac{f}{2\kappa}. (5.39)

Therefore, we provide:

Theorem 7.

In a P(GR4)P\left(GR_{4}\right) GRW spacetime satisfying f()f(\mathcal{R}) gravity with constant Ricci scalar, pp and σ\sigma are described by (5.38) and (5.39), respectively.

5.3. Energy Conditions

In general relativity, energy conditions are vital tools to study black holes and wormholes in numerous modified gravities. To specify certain energy conditions in our current investigation of the f()f(\mathcal{R}) gravity, we must find the effective isotropic pressure peffp^{eff} and the effective energy density σeff\sigma^{eff} (see, [37]-[15]).

Equation (5.34), can be rewritten in the subsequent form

lk2glk=κfTlkeff,\mathcal{R}_{lk}-\frac{\mathcal{R}}{2}g_{lk}=\frac{\kappa}{f^{\prime}}T_{lk}^{eff}, (5.40)

in which

Tlkeff=Tlk+(ff)2κglk.T_{lk}^{eff}=T_{lk}+\frac{(f-\mathcal{R}f^{\prime})}{2\kappa}g_{lk}. (5.41)

Then equation (1.2) reduces to

Tkleff=(σeff+peff)ukul+peffgkl,T_{kl}^{eff}=(\sigma^{eff}+p^{eff})u_{k}u_{l}+p^{eff}g_{kl}, (5.42)

in which peff=p+(ff)2κp^{eff}=p+\frac{(f-\mathcal{R}f^{\prime})}{2\kappa} and σeff=σ(ff)2κ\sigma^{eff}=\sigma-\frac{(f-\mathcal{R}f^{\prime})}{2\kappa}.

In our case, using equations (5.38) and (5.39), we provide

peff=fθκ(f)2κp^{eff}=\frac{f^{\prime}\theta}{\kappa}-\frac{(\mathcal{R}f^{\prime})}{2\kappa} (5.43)

and

σeff=f(f1f3θθφ)κφ+(f)2κ.\sigma^{eff}=\frac{f^{\prime}(f_{1}-f_{3}-\theta-\theta\varphi)}{\kappa\varphi}+\frac{(\mathcal{R}f^{\prime})}{2\kappa}. (5.44)

Now we investigate the energy conditions in a P(GR4)P\left(GR_{4}\right) GRW spacetime with non-zero constant Ricci scalar obeying f()f(\mathcal{R}) gravity. With the help of equations (5.43) and (5.44), we find the NEC, TEC, DEC, WEC and SEC in this set up and they are given as follows:

Table 1
Validity of Energy Conditions in P(GR4)P\left(GR_{4}\right) GRW spacetime
Energy Condition Inequalities Conditions of validation
NEC peff+σeff0p^{eff}+\sigma^{eff}\geq 0 (f1f3θ)\mathcal{R}\leq(f_{1}-f_{3}-\theta)
TEC σeff3peff0\sigma^{eff}-3p^{eff}\geq 0 (f1f3θ4θφ)φ\mathcal{R}\leq\frac{(f_{1}-f_{3}-\theta-4\theta\varphi)}{\varphi}
WEC σeff0\sigma^{eff}\geq 0 and peff+σeff0p^{eff}+\sigma^{eff}\geq 0 (f1f3θθφ)φ\mathcal{R}\leq\frac{(f_{1}-f_{3}-\theta-\theta\varphi)}{\varphi}
and (f1f3θ)\mathcal{R}\leq(f_{1}-f_{3}-\theta)
DEC σeff0\sigma^{eff}\geq 0 and σeff±peff0\sigma^{eff}\pm p^{eff}\geq 0 f1f3θ(1+2φ)0\frac{f_{1}-f_{3}}{\theta(1+2\varphi)}\geq 0 , (f1f3θθφ)φ\mathcal{R}\leq\frac{(f_{1}-f_{3}-\theta-\theta\varphi)}{\varphi}
and (f1f3θ)\mathcal{R}\leq(f_{1}-f_{3}-\theta)
SEC σeff0\sigma^{eff}\geq 0 and σeff+3peff0\sigma^{eff}+3p^{eff}\geq 0 (f1f3θθφ)φ\mathcal{R}\leq\frac{(f_{1}-f_{3}-\theta-\theta\varphi)}{\varphi}
and (f1f3θ4θφ+φ)0(f_{1}-f_{3}-\theta-4\theta\varphi+\varphi\mathcal{R})\geq 0

6. Discussion

Because of an additional higher-order curvature terms, modified gravitational theories are thought to be the most intriguing and promising way to study the present cosmic expansion. Here, with the geometric restriction of P(GR4)P\left(GR_{4}\right) GRW\mathrm{GRW} spacetime, f()f(\mathcal{R}) gravity models are investigated and we find pp and σ\sigma are not constants. Hence, we may assert that the P(GR4)P\left(GR_{4}\right) GRW\mathrm{GRW} spacetime is consistent with the universe as it exists right now. Also, for the condition f()=f(\mathcal{R})=\mathcal{R}, the spacetime reveals dust matter era, radiation era and dark energy epoch under certain restrictions.

The initial f()f(\mathcal{R})-model, f()=+α2f(\mathcal{R})=\mathcal{R}+\alpha\mathcal{R}^{2}, (α>0\alpha>0) proposed by Starobinsky [32] was aimed at explaining cosmic inflation as a pure gravitational effect without the use of dark energy. Carroll et al. [53] presented the model f()=μ4f(\mathcal{R})=\mathcal{R}-\frac{\mu^{4}}{\mathcal{R}}, (μ>0\mu>0) to explain late-time acceleration as a scalar field. Despite their limitations, these models were still able to popularize f()f(\mathcal{R})-models in general. In [51], by choosing the model f()=α(1eα)f(\mathcal{R})=\mathcal{R}-\alpha(1-e^{-\frac{\mathcal{R}}{\alpha}}) the authors have shown that DEC, NEC and WEC have been satisfied, whereas SEC violated. Here, our findings have been assessed both analytically and graphically. Our formulation was constructed using the analytical technique, and one cosmological model, f()=μctanh(c)f\left(\mathcal{R}\right)=\mathcal{R}-\mu\mathcal{R}_{c}\tanh(\frac{\mathcal{R}}{\mathcal{R}_{c}}), were evaluated for stability. For this model, we find that NEC, WEC, DEC, SEC and TEC are satisfied. Also, this finding yields the late-time acceleration of the Cosmos[51] and each result aligns with the Λ\LambdaCDM model [52].

In [54], Capozziello et al. deduced that a GRW spacetime of dimension nn with kClijk=0\nabla_{k}C^{k}_{lij}=0 reveals a perfect fluid type EMT for any f()f(\mathcal{R}) gravity model. Hence, from Theorem 2, we state the subsequent:

A P(GR4)P\left(GR_{4}\right) GRW\mathrm{GRW} spacetime represents a perfect fluid type EMT for any model of f()f(\mathcal{R}) gravity.

A collective study of validation of energy conditions are investigated and the result is mentioned in Table 1. It is seen that the presence of exotic matter is not required for our case.

7. Declarations

7.1. Funding

NA.

7.2. Code availability

NA.

7.3. Availability of data

NA.

7.4. Conflicts of interest

The authors have no conflicts to disclose.

Acknowledgment

We would like to thank the referee and the Editor for reviewing the paper carefully and their valuable comments to improve the quality of the paper.

References

  • [1] B. O’Neil, Semi-Riemannian Geometry with Applications to the Relativity, Academic Press, New York-London, ISBN: 0125267401, 1983.
  • [2] L. Alias, A. Romero, M. Sanchez, Uniqueness of complete spacelike hypersurfaces of constant mean curvature in generalized Robertson-Walker spacetimes, Gen. Relativity Gravitation 27 (1995) 71-84.
  • [3] B.Y. Chen, Pseudo-Riemannian Geometry, δ\delta-invariants and Applications, World Scientific, 2011.
  • [4] B.Y. Chen, A simple characterization of generalized Robertson–Walker spacetimes, Gen. Relativity Gravitation 46 (2014) 1833.
  • [5] C.A. Mantica, L.G. Molinari, On the Weyl and Ricci tensors of Generalized Robertson-Walker space-times, J. Math. Phys. 57 (2016) 102502.
  • [6] C.A. Mantica, L.G. Molinari, Generalized Robertson Walker spacetimes-A survey, Int. J. Geom. Methods Mod. Phys. 14 (2017) 1730001.
  • [7] M. Sánchez, On the geometry of generalized Robertson-Walker spacetimes: geodesics, Gen. Relativity Gravitation 30 (1998) 915-932.
  • [8] S.W. Hawking, G.F.R. Ellis, The Large Scale Structure of Space-Time, Cambridge University Press, London, 1973.
  • [9] P.H. Chavanis, Cosmology with a stiff matter era, Phys. Rev. D 92 (2015) 103004.
  • [10] L.P. Eisenhart, Riemannian Geometry, Princeton University Press, ISBN: 9780691023533, 1949.
  • [11] H.J. Schmidt, Fourth order gravity: equations, history, and applications to cosmology, Int. J. Geom. Methods Mod. Phys. 04 (2007) 209.
  • [12] E.M. Patterson, Some theorems on Ricci-recurrent spaces, J. London Math. Soc. 27 (1952) 287–295.
  • [13] U.C. De, N. Guha, D. Kamilya, On generalized Ricci-recurrent manifolds, Tensor (N.S.) 56 (1995) 312–317.
  • [14] S. Mallick, A. De, U. C. De, On generalized Ricci recurrent manifolds with applications to relativity, Proc. Nat. Acad. Sci. India Sect. A 83 (2013) 143–152.
  • [15] A. De, T.-H. Loo, R. Solanki, P. K. Sahoo, A conformally flat generalized Ricci recurrent spacetime in f(R)f(R)-gravity, Phys. Scr. 96 (2021) 085001.
  • [16] A. De, T.-H. Loo, R. Solanki, P.K. Sahoo, How a conformally flat (GR)4 impacts Gauss-Bonnet gravity? Fortschritte der Physik 69 (2021) 2100088.
  • [17] C.A. Mantica, Y.J. Suh, Pseudo-Z symmetric space-times, J. Math. Phys. 55 (2014) 042502.
  • [18] F.O. Zengin, m-Projectively flat spacetimes, Math. Reports 14 (2012) 363-370.
  • [19] F. Mofarreh, K. De, U.C. De, Characterizations of a spacetime admitting Ψ\Psi-conformal curvature tensor, Filomat 37 (2023) 10265–10274.
  • [20] K. De, U.C. De, Some geometric and physical properties of pseudo ψ\psi-conharmonically symmetric manifolds, Quaestiones Mathematicae 46 (2022) 939-958.
  • [21] S. Güler, S.A. Demirbağ, A study of generalized quasi-Einstein spacetimes with applications in general relativity, Internat. J. Theoret. Phys. 55 (2016) 548–562.
  • [22] P. Zhao, U.C. De, B. Unal, K. De, Sufficient conditions for a pseudosymmetric spacetime to be a perfect fluid spacetime, Int. J. Geom. Methods Mod. Phys. 18 (2021) 2150217.
  • [23] O. Bertolami, C.G. Böhmer, T. Harko, F.S.N. Lobo, Extra force in f()f\left(\mathcal{R}\right) modified theories of gravity, Phys. Rev. D 75 (2007) 104016.
  • [24] S. Capozziello, S. Nojiri, S.D. Odintsov, The role of energy conditions in f()f\left(\mathcal{R}\right) cosmology, Phys. Letter B 781 (2018) 99-106.
  • [25] E. Elizalde, R. Myrzakulov, V.V. Obukhov, D. Sáez-Gómez, Λ\LambdaCDM epoch reconstruction from F(,G)F\left(\mathcal{R},G\right) and modified Gauss–Bonnet gravities, Classical Quantum Gravity 27 (2010) 095007.
  • [26] T. Harko, F.S.N. Lobo, S. Nojiri S.D. Odintsov, f(R,T)f\left(R,T\right)-gravity, Phys. Rev. D 84 (2011) 024020.
  • [27] M.D. Laurentis, M. Paolella, S. Capozziello, Cosmological inflation in f(,G)f\left(\mathcal{R},G\right) gravity, Phys. Rev. D, 91 (2015) 083531.
  • [28] A.K. Raychaudhuri, S. Banerji, A. Banerjee, General relativity, astrophysics, and cosmology, Springer-Verlag New York, Inc., 1992.
  • [29] K.L. Duggal, R. Sharma, Symmetries of Spacetimes and Riemannian Manifolds, Springer New York, NY, 1999.
  • [30] C. Barcelo, M. Visser, Twilight for the energy conditions?, Int. J. Mod. Phys. D 11 (2002) 1553.
  • [31] H.A. Buchdahl, Non-linear Lagrangians and cosmological theory, Mon. Not. Roy. Astron. Soc. 150 (1970).
  • [32] A. A. Starobinsky, A new type of isotropic cosmological models without singularity, Phys. Letter B 91 (1980) 99-102.
  • [33] T. Harko, F. S. N. Lobo, f(R,Lm)f(R,L_{m})-gravity, Eur. Phys. J. C 70 (2010) 373–379.
  • [34] N. Katirci, M. Kavuk, f(R,TμνTμν)f(R,T_{\mu\nu}T^{\mu\nu}) gravity and Cardassian-like expansion as one of its consequences, Eur. Phys. J. Plus 129 (2014) 163.
  • [35] S. Capozziello, V. F. Cardone, V. Salzano, Cosmography of f(R)f(R) gravity, Phys. Rev. D 78 (2008) 063504.
  • [36] S. Capozziello, R. D’Agostino, O. Luongo, Extended Gravity Cosmography, Int. J. Mod. Phys. D (2019) doi:10.1142/S0218271819300167
  • [37] A. De, T.-H. Loo, S. Arora and P.K. Sahoo, Energy condition for a (WCLOSE(WRSOPEN)4)_{4} spacetime in f(R)f(R)-gravity, Eur. Phys. J. Plus, 136 (2021) 218.
  • [38] K. De and U.C. De, Investigations on solitons in f()f(\mathcal{R})-gravity, Eur. Phys. J. Plus 137 (2022) 180.
  • [39] U.C. De, K. De, F.O. Zengin, S.A. Demirbag, Characterizations of a spacetime of quasi-constant sectional curvature and ()\mathcal{F}(\mathcal{R})-gravity, Fortschritte der Physik 71 (2023) 2200201.
  • [40] S. Tsujikawa, Observational signatures of f(R)f(R) dark energy models that satisfy cosmological and local gravity constraints, Phys. Rev. D, 77 (2008) 023507.
  • [41] L. Amendola, S. Tsujikawa, Phantom crossing, equation-of-state singularities, and local gravity constraints in f(R)f(R) models, Phys. Letter B 660 (2008) 125-132.
  • [42] A.D. Dolgov, M. Kawasaki, Can modified gravity explain accelerated cosmic expansion?, Physics Letters B, 573 (2003) 1-4.
  • [43] L. Amendola, R. Gannouji, D. Polarski, S. Tsujikawa, Conditions for the cosmological viability of f(R)f(R) dark energy models, Phys. Rev. D 75 (2007) 083504.
  • [44] Y.S. Song, W. Hu, I. Sawicki, Large scale structure of f(R)f(R) gravity, Phys. Rev. D 75 (2007) 044004.
  • [45] D.G. Boulware, S. Deser, String-Generated Gravity Models, Phys. Rev. Lett. 55 (1985) 2656.
  • [46] C.A. Mantica, L.G. Molinari, Weyl compatible tensors, Int. J. Geom. Methods Mod. Phys. 11 (2014) 1450070.
  • [47] E. Bertschinger, A.J.S. Hamilton, Lagrangian evolution of the Weyl tensor, Astroph.J. 435 (1994) 1-7.
  • [48] M. Sánchez, On the geometry of static spacetimes, Nonlinear Analysis: Theory, Methods & Applications 63 (2005) 455-463.
  • [49] H. Stephani, D. Kramer, M. Mac-Callum, C. Hoenselaers, E. Herlt, Exact Solutions of Einstein’s Field Equations, Cambridge University Press, Cambridge, 2009.
  • [50] T.P. Sotiriou, V. Faraoni, f(R)f(R) theories of gravity, Modern Physics, 82 (2010) 451-497.
  • [51] T.-H. Loo, A. De, S. Mandal, P.K. Sahoo, How a projectively flat geometry regulates F(R)-gravity theory? Phys. Scr. 96 (2021) 125034.
  • [52] N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A.J. Banday, R.B. Barreiro, N. Bartolo, S. Basak, R. Battye, Planck 2018 results-VI. Cosmological parameters, Astron. Astrophys. 641 (2020) A6.
  • [53] S.M. Carroll, V. Duvvuri, M. Trodden, M. S. Turner, Is cosmic speed-up due to new gravitational physics?, Phys. Rev. D 70 (2004) 043528.
  • [54] S. Capozziello, C.A. Mantica, L.G. Molinari, Cosmological perfect fluid in f(R)f(R) gravity, Int. J. Geom. Methods Mod. Phys. 16 (2019) 1950008.