f ( R , T ) Gravity and Constant Jerk Parameter in FLRW Spacetime †

: It is well known that the universe is undergoing accelerated expansion during recent times and that it underwent a decelerated expansion in early times. The deceleration parameter, essentially the second derivative of the scale factor, can be used to describe these eras, with a negative parameter for acceleration and a positive parameter for deceleration. Apart from the standard Λ CDM model in general relativity, there are many cosmological models in various other theories of gravity. In order to describe these models, especially the deviation from general relativity, the jerk parameter was introduced, which is basically the third derivative of the scale factor. In the Λ CDM model in general relativity, the jerk parameter j is constant, and j = 1. The constant jerk parameter, j = 1, leads to two different scale factor solutions, one power law and the other exponential. The power-law solution corresponds to a model in which our universe expands with deceleration, while the exponential solution corresponds to a model in which it expands by accelerating. In this study, the cosmological consequences of such a selection of the jerk parameter on the non-minimally coupled f ( R , T ) theory of gravity (where R is the Ricci scalar, and T is the trace of the energy–momentum tensor) and the dynamic properties of these models are investigated on a ﬂat Friedmann–Lemaitre–Robertson–Walker backgfround.


Introduction
It is known that while our universe was in a period of decelerated expansion in the past, it has been in a period of accelerated expansion during recent times. The accelerated expansion has been indicated by several astrophysical observations [1][2][3][4]. Einstein's General Relativity (GR) cannot adequately explain this late-time acceleration. Therefore, various suggestions have been put forward to explain this acceleration. One of them is the concept of dark energy (DE), which connects accelerated expansion to an energy of unknown nature with negative pressure. Apart from the standard Λ cold dark matter (ΛCDM) model in GR, in this framework, various DE terms are added to the energy-momentum tensor on the right-hand side of Einstein's field equations [5][6][7][8][9][10][11][12][13]. Another attempt to explain this accelerated expansion is modified gravitation theories. Amongst these, f (R, T) theory is important because it is a modified theory through a function that depends on both the curvature (R) and the matter source T [14]. In other words, the geometry-matter coupling is considered in f (R, T) theory. In this theory, various universe models are studied by taking different forms of the f (R, T) function. However, in general, it would not be wrong to say that non-minimally coupling forms are less considered [15][16][17][18].
In the present study, we consider f (R, T) gravity on the background of the flat Friedman-Lemaitre-Robertson-Walker (FLRW) universe for a non-minimal coupling form, To solve the field equations, we assume that the jerk parameter j is constant, and j = 1. One of the two solutions to this assumption corresponds to the past period when the universe expanded by decelerating, and the other corresponds to the late period when it expanded by accelerating. In Section 2, the field equations of f (R, T) gravity are given, the solutions of the field equations and their consequences are derived in Section 3, and the conclusion is formed in Section 4.

Field Equations of f (R, T) Gravity
The f (R, T) gravity field equations are based on the gravitational action [14] where f (R, T) is an arbitrary function of the Ricci scalar R and trace T of the energymomentum tensor T ij , g is the determinant of the metric tensor g ij , and L m is the Lagrange density of the matter. The variation in the gravitational action S in Equation (1) with respect to the metric g ij gives the field equations of f (R, T) gravity: Here, ∇ i is the covariant derivative, ≡ ∇ i ∇ i is the d'Alembertian operator, and When calculating the variation, L m is taken to depend only on the metric tensor and not on its derivatives. If the matter assumed to fill the universe is considered to be a perfect fluid, then the energy-momentum tensor attached to it is where ρ and p is the energy density and pressure, respectively, and u i (u i u i = 1) is the four-velocity vector of the fluid. When the matter's Lagrangian L m = −p is taken, Θ ij becomes Substituting this into Equation (2), the field equations take the form Three different functional forms of the function f (R, T) are considered in the study in Ref. [14]. In this study, we adopt the last of the forms (which is not usually chosen): where f i (i = 1, 2, 3.) are arbitrary functions of their arguments. In this case, the field equations' Equation (6) becomes: Here, a prime ( ) denotes a derivative with respect to the argument for any f i function. Now, we customize the functional form even more, taking it as to obtain where λ is a coupling constant. If λ = 0 is put into Equation (10), it is clearly seen that the equations are reduced to the field equations of GR.
On the other hand, the usual conservation equation of GR is not valid in this theory, but it takes the following form for Equation (9) [18]: where H =ȧ/a is the Hubble parameter.

Modified Field Equations in the Flat FLRW Background
The homogeneous and isotropic flat FLRW metric is where a(t) is the time-dependent scale factor. For the flat FLRW model, the modified field equations defined in Equation (10) lead to the following two independent equations The trace of Equation (4) gives T = ρ − 3p; hence, we haveṪ =ρ − 3ṗ andT =ρ − 3p. Substituting these expressions into Equation (14), we obtain (15) As can be seen, this last equation includes not only ρ and p, but also their first and second derivatives with respect to time. To remove this mathematical difficulty, by adopting the barotropic EoS p = ωρ, we can rewrite Equation (11) aṡ where ω is the EoS parameter. One more differentiation of the last equation with respect to time gives Here, R = −6(Ḣ + 2H 2 ),Ṙ = −6(Ḧ + 4HḢ) andR = −6( ... H + 4Ḣ 2 + 4HḦ), and the Hubble parameter H is defined as Now, the field Equations (13) and (15) become Using these two equations, one obtains The last equation is known as the generalized Raychaudhuri equation. One can solve ρ from the generalized Raychaudhuri equation in terms of its first and second time derivatives as Now, we need to know the expression of the Hubble parameter to obtain the explicit form of ρ. For this task, we limit ourselves to the assumption of a constant jerk parameter (j = 1). The jerk parameter is basically the third derivative of the scale factor. Hence, our assumption is Note that in the ΛCDM model of GR, the jerk parameter j is constant, and j = 1. The integration of Equation (23) leads to two different scale factor solutions; one is power law, and the other one is exponential [19], as follows: where c, α, and β are the constants of integration. The power-law solution (24) is important for explaining the early universe, and the exponential solution (25) is important for the late universe. Now, regarding Equation (18) for the Hubble parameter and the definition of the deceleration parameter q, we obtain: In the remainder of our work, we consider these two solutions separately. The power-law solution Equation (24) yields The sign of the deceleration parameter is related to the shape of the expansion of the universe. A positive q indicates a decelerated expansion, and a negative q indicates an accelerated expansion. In this model, q equals 1/2, indicating that the power-law solution depicts a decelerated expanding universe model.
Using Equations (16) and (17), for the time derivatives of ρ, one can obtain the explicit expression of ρ. Nevertheless, the expression is very lengthy and complicated. Therefore, we consider Equation (16) with Equation (27) Upon integration, we obtain: where ρ 0 is a constant of integration. Equation (25), the exponential solution, yields: Since q is negative, the exponential solution leads to the accelerated expansion model. For this case, following the same procedure as in the previous model, Equation (16) becomeṡ and its integration yields (32) Figure 1 shows the temporal changes of the energy density for the decelerating (left panel) model for the choice of the integrating constants c = 0, c 1 = 0, and ρ 0 = 1; for the accelerating model (right panel) for the choice of the integrating constants β = 1 and ρ 0 = 1; and for different values of the EoS parameter ω. We can see that ρ is a positive value for all t. These are necessary conditions for the physicality of the model.

Conclusions
We have studied dark energy with j = 1 in f (R, The simpler form f 2 (R) f 3 (T) = λT is usually studied in the literature. It is possible to obtain a viable model with a transition from deceleration to acceleration. We have not considered graphs of the geometrical parameters nor the pressure and equation of the state, or observational constraints. We will report on these aspects elsewhere.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.