Isotropic compact stars in four-dimensional Einstein–Gauss–Bonnet gravity coupled with scalar field: reconstruction of model

It has been suggested that the Einstein–Gauss–Bonnet theory coupled with a scalar field (EGBS) may allow us to obtain physically viable models of celestial phenomena such that the scalar field effect is active in standard four dimensions. We consider the spherically symmetric and static configuration of the compact star and explain the consequences of the EGBS theory in the frame of stellar modeling. In our formulation, for any given static profile of energy density ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} with spherical symmetry and the arbitrary equation of state (EoS) of matter, we can construct a model which reproduces the profile. Because the profile of the energy density determines the mass M and the radius Rs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_s$$\end{document} of the compact star, an arbitrary relation between the mass M and the radius Rs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_s$$\end{document} of the compact star can be realized by adjusting the potential and the coefficient function of the Gauss–Bonnet term in the action of EGBS theory. This could be regarded as a degeneracy between the EoS and the functions characterizing the model, which indicates that the mass–radius relation alone is insufficient to constrain the model. Here, we investigate a novel class of analytic spherically symmetric interior solutions by the polytropic EoS. We discuss our model in detail and show that it is in agreement with the necessary physical conditions required for any realistic compact star, confirming that EGBS theory is consistent with observations.


Introduction
Although Einstein's theory of general relativity (GR) is successful at present, and can forecast and elucidate increased observational data, there are strong reasons to expect that it a e-mail: nashed@bue.edu.eg (corresponding author) b e-mail: nojiri@gravity.phys.nagoya-u.ac.jp must be modified due to its shortfall in the quantization of gravity and explaining the recent observational puzzles in modern cosmology leading to the study of amended theories of gravity.
The Lovelock gravitational theories [1] are of special interest since they are Lagrangian-based theories that can give conserved covariant field equations which do not include derivatives higher than the second degree. In this regard, Lovelock's theories are the physical extensions of GR. The Gauss-Bonnet (GB) theory is considered the first physical nontrivial expansion of Einstein's GR. This theory is meaningful if its spacetime is greater than four-dimensional, in which the GB invariant can create a rich phenomenology. Through the use of Chern's theorem [2], it can be shown that in four dimensions, the GB expression is a non-dynamical term because the GB invariant becomes a total derivative. To make the GB expression a dynamical one in four dimensions, we must invoke a novel scalar field with a canonical kinetic term coupling to the GB term [3][4][5][6][7][8][9][10][11][12] as stimulated, for example, by low-energy effective actions stemming from string theory, such as the Einstein-dilaton-GB models [5,[13][14][15]. Actually, because of Lovelock's theorem, in principle, all amended gravitational theories in four dimensions will have extra degrees of freedom, which can be considered as new basic fields. The exact solutions of the gravitational system supply scientific society with a simple test of spacetime and evaluation of observable forecasts. Nevertheless, amended gravitational theories with new basic field(s) usually provide equations of motion with high intractability so that the evaluations become analytically out of the question. To address such an issue, one is forced either to apply perturbation theory, which is not well qualified in the strong gravitational field, or to defy numerical methods [16]. However, the field equations of GR coupled with matter have conformal invariance since they possess the constant Ricci scalar curvature on-shell, limiting the spacetimes and permitting analytic solutions to be easily derived. An example of such a theory that has conformal invariance and yields simple analytic solutions is the electrovacuum, whose Reissner-Nordström (Kerr-Newman) solution was the first-ever discovered static (spinning) black hole (BH) with a matter source. Another model is the gravitational theory coupled with a conformal scalar field, in which the matter action obeys the conformal invariance and has the form where R is the Ricci scalar and ξ is the scalar field. The field equations of the above action give a solution with the nohair theorem (see, e.g., Ref. [17] for a review) and the static Bocharova-Bronnikov-Melnikov-Bekenstein BH [18][19][20], which has been much debated. Gravitational theory with a conformal scalar field and its solutions have been discussed in recent years because of its compelling properties (see, e.g., Refs. [21][22][23][24][25][26][27][28][29][30][31] and references therein). As we discussed above, in four dimensions, the GB term is topological and does not yield any dynamical effect. Nevertheless, when the GB term is non-minimally coupled with any other field, such as a scalar field ξ , the output dynamics are nontrivial. Many cosmological proposals have been presented in the recent literature (  and references therein). In the frame of astrophysics, however, to the best of our knowledge, the GB theory with a non-minimal coupling of a scalar field via potential and coefficient function has not been tackled, although there are some pioneering works such as [8]. It is the aim of the present study to derive exact spherically symmetric interior solutions of this theory and discuss the physical consequences. By using our formulation, we can construct a model which reproduces any given profile of the energy density ρ for an arbitrary equation of state (EoS) of matter. The mass M and the radius R s of the compact star are determined by the profile of the energy density, and therefore we can obtain an arbitrary relation between the mass M and the radius R s of the compact star by adjusting the scalar potential and the coefficient function of the GB term in the action of the Einstein-Gauss-Bonnet gravity coupled with a scalar field (EGBS), which could be a kind of degeneracy between the EoS and the functions characterizing the model. Therefore, we find that the mass-radius relation alone is not sufficient to constrain the model.
The remainder of the paper is organized as follows: In Sect. 2, we describe the fundamentals of the EGBS, and we apply the field equation of the EGBS theory to a spherically symmetric spacetime and derive the full system of the differential equation. Here we show that we can construct a model which reproduces any given profile of the energy density ρ for an arbitrary EoS of matter. Also in Sect. 3, we give the form of a polytropic EoS as an example and the form of one of the metric potentials as an input and then derive all the unknown functions including the profile of the scalar field, the coefficient function, the potential of the scalar field, and the form of another metric potential. Section 4 discusses the physical conditions that must be satisfied for any real stellar configuration. In Sect. 5, we discuss the physical properties analytically and graphically, showing that the solutions have realistic physical properties. In Sect. 7, we discuss the issue of stability by using the adiabatic index and show that our model satisfies the adiabatic index; that is, the value of the index is greater than 4/3, which is the condition of stability. The final section is reserved for the conclusion and discussion of the present study.

Gauss-Bonnet theory coupled with scalar through f (ξ )
Now we consider the EGBS in N dimensions. This theory takes the following amended action, where ξ is the scalar field and V is the potential which is a function of ξ , f (ξ ) is an arbitrary function of the scalar field, and S M is the matter action, where we assume that matter is to couple minimally to the metric, i.e., we are working in the socalled Jordan frame. In four dimensions, i.e., when N = 4, the aforementioned action is physically nontrivial because the GB invariant term G is coupled with the real scalar field ξ through the coupling f (ξ ). Because of this coupling, the Lagrangian is not a total derivative but contributes to the field equations of the system. The variation of the action (3) w.r.t. the scalar field ξ yields the following equation The variation of the action (3) w.r.t. the metric g μν yields the following field equations Through the use of the Bianchi identities, in Eq. (5), we obtain The field equations (4) and (7) are the full system of equations describing the theory under consideration. In the fourdimensional case, i.e., N = 4, Eq. (7) yields In the present study, we assume that the scalar field ξ is a function of the radial coordinate r and therefore the function f (ξ ) depends only on r , i.e., f (r ) ≡ f (ξ (r )), because we deal with static and spherically symmetric spacetime, In the following section, we study the system of field equations (4) and (8) and try to find the analytic form of the unknown functions when T μν = 0.

Four-dimensional spherically symmetric interior solution in EGBS
For the metric in (9), the (t, t)-component of the field equation Eq. (8) has the following form where the (r, r )-component is given by and the (θ, θ )and (φ, φ)-components are The field equation of the scalar field (4) takes the following form Here, ρ is the energy density and p is the pressure of matter, which we assume to be a perfect fluid and which satisfies an EoS, p = p (ρ). The energy density ρ and the pressure p satisfy the following conservation law The conservation law is also derived from Eqs. (10), (11), (12), and (13). Here we have assumed that ρ and p depend only on the radial coordinate r . Other components of the conservation law are trivially satisfied. If the EoS ρ = ρ( p) is given, Eq. (14) can be integrated as Because Eq. (14) and therefore (15) can be obtained from Eqs. (10), (11), (12), and (13), as long as we use (15), we forget one equation in Eqs. (10), (11), (12), and (13). In the following, we do not use Eq. (13). Inside the compact star, we can use Eq. (15), but outside the star, we cannot use Eq. (15). Instead of using Eq. (15), we may assume the profile of a = a(r ) so that a(r ) and a (r ) are continuous at the surface of the compact star. By Eq. (10) + Eq. (11), we obtain On the other hand, Eq. (10) − Eq. (11) gives Furthermore, Eq. (10) − Eq. (12) gives which can be regarded with the differential equation for f and therefore for f if a = a(r ), a 1 = a 1 (r ), ρ = ρ(r ), and p = p(r ) are given and the solution is given by Here, c 1 and c 2 are constants of the integration. Let us assume the r -dependencies of ρ and a 1 , ρ = ρ(r ) and a 1 = a 1 (r ). Then by using the EoS p = p(ρ), we find the r -dependence of p, p = p(r ) = p (ρ (r )). Furthermore, by using (15), we find the r -dependence of a, a = a(r ). However, Eq. (15) is not valid outside the compact star because ρ and p, of course, vanish there. Then outside the compact star we may properly assume the profile of a(r ) so that a(r ) and a (r ) are continuous at the surface, that is, the boundary of the compact star, and coincide with a(r ) and a (r ) obtained from (15). Therefore, by using (19), we find the r -dependence of f , f = f (r ), and by using Eqs. (16) and (17), we find the r dependencies of V and ξ , V = V (r ) and ξ = ξ(r ). By solving ξ = ξ(r ) with respect to r , r = r (ξ ), we find f and V as , which realize the model which has a solution given by ρ = ρ(r ) and a 1 = a 1 (r ).
We should note, however, that the expression of ξ in (17) gives a constraint, 2 a 1 (ρ + p) so that the ghost can be avoided. If Eq. (20) is not satisfied, the scalar field ξ becomes purely imaginary. We may define a new real scalar field ζ by ξ = iζ i 2 = −1 , but because the coefficient in front of the kinetic term of ζ becomes negative, ζ is a ghost, that is, a non-canonical scalar field. The existence of the ghost generates the negative norm states in the quantum theory, and therefore the theory becomes inconsistent. When we consider compact stars like neutron stars, we often consider the following EoS 1. Energy-polytrope with constants K and n. It is known that for the neutron stars, n can take the value 0.5 ≤ n ≤ 1.

Mass-polytrope
where ρ m is rest mass energy density and K m , N are constants.
Now let us study the case of the energy-polytrope (21) in detail, in which we can rewrite the EoS as follows For the energy-polytrope, Eq. (15) takes the following form Here, c is a constant of the integration. Similarly, in the case of the mass-polytrope (22), we obtain Here,c is again a constant of the integration. Under one of the above EoS, we may assume the following profile of ρ = ρ(r ) and a 1 = a 1 (r ), just for an example, Here, r 0 is a constant, ρ c is a constant expressing the energy density at the center of the compact star, R s is also a constant corresponding to the radius of the surface of the compact star, and M is a constant corresponding to the mass of the compact star, When r → ∞, a 1 behaves as a 1 (r ) ∼ 1 − 2M r , and therefore M can be regarded as the mass of the compact star. Equation (27) gives the M-r relation, that is, the relation between the mass and the radius of the compact star when r = R s . We also note that we need to choose r 0 large enough that a 1 is positive. In order for a 1 in (26) to be positive, we require We should also note that when r → 0, a 1 behaves as a 1 (r ) ∼ 1 − 2Mr 2 r 3 0 . 1 Therefore, a 1 (r ) vanishes at the center r = 0, a 1 (r = 0) = 0, and thus there is no conical singularity.
As an example, we use the energy-polytope as the EoS by choosing n = 1 just for simplicity. Then Eq. (24) gives which gives a = 8e c Kρ c r Outside the star, we assume that a(r ) = a 1 (r ) in (26), and therefore Because a(r ) and a (r ) should be continuous at the surface r = R s , we obtain 1 It is well known that the junction conditions for the matching of two spacetime manifolds have further restrictions in the EGBS gravity [88]. With regard to a static configuration, this is not a real problem, since the interior will match to vacuum, and so the pressure will still vanish at the surface r = R s , as we will show below.
By deleting e c in the two equations in (32), we obtain where we have used Eq. (27) when r = R s . Because r 0 should be positive, we find Then, by using (19), we find the r -dependence of f , f = f (r ), and by using Eqs. (16) and (17), the r dependencies of V and ξ , V = V (r ) and ξ = ξ(r ), are determined. If we can solve ξ = ξ(r ) with respect to r , r = r (ξ ), we find f and V as functions of ξ , f (ξ ) = f (r (ξ )), V (ξ ) = V (r (ξ )).
Inside the compact star, by using Eqs. (26) and (29), we find that the GB term G behaves as Equation (36) shows that the GB term does not vanish, and it depends on the mass of the star. Now we calculate the form of f (r ) using the data given in Eqs. (29), (32), and (26). The explicit form of f (r ) is displayed in Appendix A. The form of ξ(r ), by using the data given in Eqs. (26), (29), and (32), is also displayed in Appendix A. Finally, we calculate the explicit form of V (r ) using the data given in Eqs. (29), (32), and (26), and list the results in Appendix A.
To complete our study, we solve Eq. (A4) asymptotically and obtain The above equation is valid provided that the constant C 12 < 0. Now, using Eq. (37) in (A2), we obtain f (ξ ) as Also using Eq. (37) in (A6), we obtain V (ξ ) as A final remark that we should stress is that using Eqs. (26), (29), (A2), and the constraints (35) with R s = 1, one can easily show that the inequality (20) holds. We have four differential equations for seven unknown functions, as shown in Eqs. (10), (11), (12), and (13), that is, ρ, p, V , ξ , f , a, and a 1 . As a result, we need to require three additional conditions to close such a system. One of these extra conditions is the continuity equation given by Eq. (14). The second condition is the polytropic EoS given by Eq. (21). The third is the profile of the energy density of matter given by Eq. (26). When these additional conditions are combined with Eqs. (10), (11), (12), and (13), the system is in a closed form, allowing all seven unknown functions to be explicitly fixed.

Ingredient requirements for a real physical stellar configuration
For a physically reliable isotropic stellar model, the solution has to satisfy the conditions inside the stellar configurations as follows: • The metric potentials a(r ) and a 1 (r ), and the energymomentum components ρ and p should be well defined at the center of the star and should have a regular behavior and have no singularity in the interior of the star. • The density ρ must be positive in the stellar interior, i.e., ρ ≥ 0. Moreover, its value at the center of the star must be finite, positive, and decreasing to the boundary of the star, i.e., dρ dr ≤ 0. • The pressure p should have a positive value inside the fluid configuration, i.e., p ≥ 0. In addition, the derivative of the pressure should yield a negative value inside the star, i.e., d p dr < 0. At the surface of the star, r = R s , the pressure p should vanish. • The causality condition which should be satisfied to obtain a realistic model, i.e., the speed of sound should be less than 1 (provided that the speed of light is c = 1) in the interior of the star, i.e., 1 ≥ d p dρ ≥ 0. • To obtain a stable model, the adiabatic index must be greater than 4 3 .
It is time to analyze the above conditions to see whether we have a real isotropic star.

Physical behavior of our model
To test whether our model given by Eqs. (22) and (24) agrees with a real stellar construction, we discuss the following issues: 5.1 Non-singular model 1. The metric potentials of this model satisfy a(r → 0) = e c (1 + Kρ c ) 4 and a 1 (r → 0) = 1 , (40) which yields that the metric potentials have finite values at the center of the star configuration. Additionally, the derivatives of these metric potentials vanish at the center of the star, i.e., a (r → 0) = a 1 (r → 0) = 0. If the derivatives do not vanish even if they are finite, there appear conical singularities at the center. The above constraints yield that the metric is regular at the center and that the metric has good behavior in the interior of the star. 2. Density (26) and pressure (21), at the center, have the form The above Eq. (41) clearly shows that the density and pressure at the center of the star always have positive values if ρ c > 0 and K > 0; otherwise they become negative.
3. The gradients of density and pressure of our model are given respectively as Here, ρ = dρ dr and p = d p dr . Equation (42) shows that the derivatives of density and pressure are negative. Furthermore, because they vanish at the center of the star, the conical singularities do not appear. 4. The velocity of sound using relativistic units, i.e., (c = G = 1), are derived as [89] v r Now we are ready to plot all the above conditions to examine their behaviors using the numerical constraints listed in Eq. (35). In Fig. 1a and b, we present the behavior of metric potentials. As Fig. 1 shows, the metric potentials assume the values a 1 (r → 0) = 1 and a(r → 0) = 0.5 for r = 0, which ensure that both of the metric potentials have finite and positive values at the center of the star. Now we plot the energy density and pressure, listed by Eqs. (21) and (26) in Fig. 2. Figure 2 shows that the energy density and pressure are positive, which is in agreement for a realistic stellar configuration. Additionally, as Fig. 2a and b indicate, the density and pressure have high values at the center and decrease toward the boundary, which is relevant for a realistic star. Figure 3 shows that the derivatives of density and pressure have negative values, which ensure the decrease in density and pressure throughout the stellar configuration.
In Fig. 4, we plot the speed of sound and the mass-radius relation. As Fig. 4a shows, the speed of sound is less than 1, which confirms the non-violations of causality condition in the interior of the stellar configuration. Moreover, Fig. 4c shows that the compactness of our model is constrained by 0 < C < 0.55, where C = M r in the stellar configuration. As Fig. 4a shows, the causality condition is satisfied, which is one of the advantages in this study due to the procedure we follow, although in the frame of GR, it is shown that this condition is not satisfied [90]. We may infer that the procedure used in this study is responsible for the correction in the behavior of the causality condition. Moreover, also as in [90], it is shown that the maximum mass lies in the range 0.2 M . In our model, however, due to the procedure we follow in this study, the maximum mass is about 0.25 M , as shown in Fig. 4b, which could be used to be compared with the recent data.  (26) and (29) vs. the radial coordinate r using the constraints (35) Fig. 2 Plot of the energy density and pressure of (21) and (26) vs. the radial coordinate r using the constraints (35) Figure 5 shows the behavior of the energy conditions. In particular, Fig. 5a-c indicate the positive values of the NEC, WEC, and SEC energy conditions, which ensure that all the conditions are verified through the stellar configuration as it should be for a physical stellar model.
In Fig. 6, we plot the EoS. As Fig. 6a shows, the EoS is not linear. It was shown in [91] that the EoS of neutral compact stars is almost a linear one, in contrast to the EoS presented in this study, which shows a nonlinear form due to the form of the pressure given by Eq. (21).

Stability of the model
Now we are ready to test the stability issue on our model using the adiabatic index. The stable equilibrium of a spherically symmetric spacetime can be investigated through the adiabatic index, which is an ingredient tool to test the stability criterion. The adiabatic perturbation, i.e., the adiabatic index , is defined as [92][93][94] = ρ + p(r ) A Newtonian isotropic sphere has a stable equilibrium if the adiabatic index > 4 3 [95]. If = 4 3 , the isotropic sphere is in neutral equilibrium.
Using Eq. (44), we obtain In Fig. 7, we depict the adiabatic index . As is clear from Fig. 7, the value of is greater than 4 3 throughout the stellar interior, and therefore the stability condition is satisfied.

Discussion and conclusions
In the present research, we considered the spherically symmetric and static configuration of the compact star by using the EGBS. In our formulation, for any given spherically symmetric and static profile of the energy density ρ and for an arbitrary EoS of matter, we can construct the model which reproduces the profile. Because the profile of the energy density determines the mass M and the radius R s of the compact star, an arbitrary relation between the mass M and the radius R s of the compact star can be realized by adjusting the potential V (ξ ) and the coefficient function f (ξ ) of the GB term in (3). This could be regarded as a degeneracy between the EoS and the functions V (ξ ) and f (ξ ) characterizing the model, which indicates that the mass-radius relation alone is insufficient to constrain the model.
As a concrete example, by using the polytrope EoS (21) and assuming the profile of the energy density ρ(r ) in (26), we have constructed a model and discussed the properties. The derived analytic solution is investigated analytically and graphically using different tests to assess the physical relevance of the derived solution.
In this regard, we discovered that the energy density and pressure decrease as the radial coordinates approach the surface of the star Fig. 1. This clearly indicates that the center of the star is highly compact and the model under consideration is valid for the region outside the center of the star. Additionally, we have explained analytically and graphically in Fig. 5 that all the energy conditions are verified throughout the interior of the stellar configuration. According to Herrera [89], any stable solution must yield a square of sound speed, v 2 , to lie in the interval v 2 ∈ [0, 1]. In this model, we have shown that the speed of sound lies in the required interval, confirming that the solution obtained in our model is stable. Also, the calculation of the adiabatic index of our model is in excellent agreement with the stability condition as shown in Fig. 7 (right panel). We have depicted the mass-radius relation as shown in Fig. 4 (middle panel). As this figure shows, the mass M takes a positive value through the interior of the star. Additionally, it is easy to prove that as r → 0, we obtain M → 0, which ensures that M is regular at the core of the star. We also showed that the procedure used in this study can significantly enhance the mass, corroborating recent observations of some massive two-solar mass neutron stars. Moreover, as Buchdahl [96] has shown, for static spherically symmetric isotropic matter content, the ratio between the mass and the radius should be M R < 4 9 . In this study, the ratio M R s = 1 4 (see Fig. 4, middle panel) shows that the Buchdahl condition is satisfied. The compactification C = M R s has been depicted in Fig. 4 (right panel), which shows that the compactness should be 0 < C < 0.55. In Fig. 6 (right panel), we have shown that the profile of the surface redshift is less than 2 as required for an isotropic model without a cosmological constant. It has been shown that the upper limit of surface redshift is 2, which is in agreement with our stellar configuration [96][97][98].
In the present study, we have assumed that a physical energy density is given by Eq. (26), as it has a finite value at the center of the star ρ c and it is finite at the surface of the star, which is consistent with realistic compact stars. Also, the metric potentials of this construction are physical because they are singularity-free, as r → 0, and have finite values at the surface of the star. Additionally, the mass of the star in the model under consideration has a finite value at the center as well as at the surface of the star. Moreover, the constructed model yields a consistent form of the GB term, and the scalar field ξ , the potential V (ξ ), and the coefficient function f (ξ ) have finite value, as r → 0. Also, we have shown that the model under consideration is stable and its adiabatic index is greater than 4 3 , which is consistent with observations. Remarkably, NICER (Neutron star Interior Composition Explorer) observations of PSR J0030+0451 and PSR J0740+6020 offer indications against the more squeezable models. The latter has significantly more mass than the former, although they are nearly the same size. So it is reasonable to suppose some processes to rationalize the nonsqueezability of a neutron star as its mass increases. On the other hand, the presence of high-mass pulsars∼ 2M such as PSR J0740+6020 is known to prefer violation of the upper sound speed conformal limit v 2 ≤ 1/3, posing another challenge for theoretical models even in low-density cases, as demonstrated by Bedaque and Steiner [99] (see also [100,101]). In their study of the pulsar PSR J0740+6020, Legred et al. [102] concluded that the conformal sound speed is strongly violated at the neutron star core, where v 2 = 0.75 with density 3.60ρ nuc. . It is important to mention that such an issue does not appear in our constructed model, as shown in Fig. 4a.
To conclude, to the best of our knowledge, that this is the first study to derive an analytic isotropic spherically symmetric interior solution in the frame of EGBS theory. From the above analysis, we ensure that the derived solution in this study met all the physical requirements of any isotropic stellar configuration in the frame of this theory. An isotropic model in the frame of Rastall's theory is derived using the technique of conformal killing vectors [103]. In this model, the authors showed that the maximum value of the compactness in their model was 0.028742 and the redshift was 0.09444. If we compare our results with those presented in [103], we see that the compactness and redshift of our model are greater than the ones presented in [103]. This means that the nonlinear form of the EoS has a greater effect on the structure of the model than the conformal killing vector. An isotropic model is also constructed in the framework of F(R, T ), where R is the Ricci scalar and T is the trace of the energy-momentum tensor. It was shown that the model constructed in [104] suffers from a violation of the dominant energy condition (DEC), whereas it is satisfied in the model under consideration. Moreover, it was shown that in the model constructed in [104], its energy density configuration is nonuniform, which corresponds to a quasi-constant density configuration, but our model did not possess a such defect.

Data Availability Statement
This manuscript has data included as electronic supplementary material. The online version of this article contains supplementary material, which is available to authorized users.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.