Existence of Compact Structures in $f(R,T)$ Gravity

The present paper is devoted to investigate the possible emergence of relativistic compact stellar objects through modified $f(R,T)$ gravity. For anisotropic matter distribution, we used Krori and Barura solutions and two notable and viable $f(R,T)$ gravity formulations. By choosing particular observational data, we determine the values of constant in solutions for three relativistic compact star candidates. We have presented some physical behavior of these relativistic compact stellar objects and some aspects like energy density, radial as well as transverse pressure, their evolution, stability, measure of anisotropy and energy conditions.


Introduction
General relativity (GR) is considered as the most fruitful theory for understanding the evolution of universe and its hidden secrets, yet the evidence of dark matter (DM) and the cosmic accelerating nature of spacetime put some challenges on this [1][2][3][4][5][6][7]. The Einstein's GR explained the cosmological phenomena in a regime of weak field, while some modifications may be needed to study the strong fields in the scenario of accelerating expansion of the universe. In this direction, Qadir et al. [8] reinforced the requirement of the modified relativistic dynamics and indicated that this modification may help to settle down the problems related to DM and quantum gravity. As a result, many techniques were used like by introducing the cosmological constant as well as the modified theories from time to time.
Modified gravitational theories (MGTs) are actually the generalization of GR in which function of the Ricci scalar (R) is substituted in the Einstein-Hilbert action. These modified gravity theories are dubbed with the names, Einstein-Λ [9], f (R) [10] (R is the Ricci scalar), f (R, T ) [11] (T is the trace of energy momentum tensor), f (G) [12] (G is the Gauss-Bonnet term) and f (R, T, R ξπ T ξπ ) gravity [13]. In the recent times, Nojiri et al. [14] presented various mathematical techniques to understand burning issues of cosmos related to bouncing cosmos. They asserted that gravity mediated by f (R) and f (G) theories could be used to realize many hidden secrets of our universe. Once can observe the pity good agreement results between the cosmological models in MGTs and the observational data [15][16][17][18]. The f (R, T ) gravity is one of the MGTs, in which the f (R) is replaced with the function of R and T [19]. It is claimed that the evidence behind the dependence of T may come from the presence of imperfect fluid or it may be some kinds of quantum effects (for further reviews on DE and MGTs, see, for instance, [20][21][22][23][24][25][26][27][28][29][30][31][32]).
In f (R, T ) gravity, many cosmological applications were discussed in [33][34][35][36][37][38]. From literature, some of them are, The non-static line element for collapsing of spherical body having anisotropic fluid were discussed in [39]. The static spherical wormhole solutions were found in [40,41]. Furthermore, the perturbation techniques were used in study of spherical stars [42]. The effects on gravitational lensing due to f (R, T ) gravity were discussed in [43]. The spherical equilibrium theme of polytropic and strange stars were investigated in [44]. Houndjo [45] constructed few observationally notable cosmic models in f (R, T ) gravity for studying matter dominated era of the expanding universe. Baffou et al. [46] applied perturbation on the spacetimes of de-Sitter and power law models in order to explore some cosmic viability bounds.
Bamba et al. [47] analyzed the effects of higher degrees of freedom coming from MGT on the dynamical features of our accelerating cosmos. Bamba et al. [53] further checked the viability regimes on the parameters of f (G) gravity models and presented some mathematically consistent cosmic zones. The stability of gravitational evolving stellar bodies have been investigated in few models of f (R) gravity by [49,51]. Das et al. [52] calculated exact relativistic models of spherical interiors in MGT and discussed the physical implications of their results on compact stars.
Yousaf and his collaborators examined the role of various curvature invariant functions on the existence as well as stability of the planar [54], spherical [55] and cylindrical [56] geometries. Sahoo with his coworkers [57] studied the viability of the spatially regular cosmos along with some other cosmological aspects in f (R, T ) gravity. Moraes et al. [58] worked out the stability of some well-known compact stars by computing their corresponding hydrostatic equations in f (R, T ) = R + 2λT gravity.
The exploration on the existence of self-gravitating compact stars have always been a source of great attention among gravitational physicists. [59][60][61][62][63][64][65][66]. In this direction, many researchers reconstructed different models for the study of anisotropic relativistic compact stars. Various physical properties, like masses of compact stars, radii, stability etc. and the moment of inertia of neutron stars were studied with the comparison developed with GR and other MGTs [67]. Virbhadra et al. [68,69] discussed the relation between naked singularity and black holes formation accompanied with a well-consistent mathematical stand point. Egeland [81] performed the modeling of neutron star by examining some mass radius relations and concluded that Λ should exist to justify the vacuum density. Sharif and Yousaf [70] also found these relations and checked the existence of compact structures in the platform of MGTs. Bhatti et al. [71] calculated dark dynamical variables and checked the dynamics of compact stars with the help of these variables. Recently, Yousaf et al. [72] investigated the stability of three different compact structures in the presence of dark sources terms mediated by quadratic, cubic and exponential f (R) formulations.
Here, we use the modified f (R, T ) gravity, which can be considered as more well established theory than that of cosmological constant. To explore the formation of relativistic anisotropic compact stellar objects in f (R, T ) gravity, we use two such a viable f (R, T ) models for three candidates of strange stars, i,e., Her X-1, SAXJ1808.4-3658 and 4U1820-30. We shall label these stars with CS1, CS2 and CS3, respectively. We find the solutions of relativistic anisotropic stellar bodies in f (R, T ) gravity, that enforce some constraints on the cosmic model parameters. In very next section, we formulate f (R, T ) equation of motions, after having the background of f (R, T ) gravity, we take anisotropic matter contents for static spherical star geometry. In section 3, we demonstrate some of the physical viable models in f (R, T ) gravity. In section 4, we explore the solutions and discussed the physical properties of relativistic compact stars through graphical illustrations. At the end, we finalize our results.

f (R, T ) gravity
The action for f (R, T ) gravity can be written as where g and L matter indicate determinant of metric tensor and the matter Lagrangian, respectively. By varying the above equation with metric tensor, we get where the notations f R , ∇ γ and f T are operators for covariant and partial derivatives of their arguments, respectively, while Θ µν can be expressed through the stress-energy tensor (T γδ ) as 3) The aim of this work is to check the different properties of some of the compact stars with locally anisotropic pressure. Now, we consider that the spherical distribution of geometric system is coupled with the following relativistic fluid in which Π ≡ P r − P ⊥ and ρ is the matter energy density and V γ and X γ are the 4-vectors corresponding to fluid and radial directions, respectively. Under the non-tilted coordinate system, the four vectors satisfy V α V α = 1 and X β X β = −1 relations. Equation (2.3) can be written, after choosing L m = ρ, as Then, the field equations (2.2) boil down to where T eff ξπ is dubbed with the effective stress energy tensor for f (R, T ) gravity whose mathematical formulations is given by Now, we wish to consider the diagonally symmetric static form of spherically symmetric spacetime as where metric coefficients a and b are the r dependent functions. Our aim is to analyze the role of anisotropicity in the modeling of some stellar toy models. So, we assume that Eq.(2.7) is composed of the locally anisotropic fluid content given in Eq. (2.4). The f (R, T ) equations of motion (2.2) for the geometry (2.7) and fluid (2.4) yield , where the over prime indicates ∂/∂r operator. In order to analyse the impact of f (R, T ) gravity on the construction of stellar models, we assume the separable form of R and T in f (R, T ) model as follows This choice of separable R and T can be regarded as a possible linear corrections in the well-known f (R) theory. By choosing f i (R) from [16] along with the linear combination of g(T ), the viable f (R, T ) model can be designed. Therefore, we suppose g(T ) = ǫT in which ǫ is a very small positive number. In this context, Eqs.(2.8)-(2.10) provide where We write a and b as a combinations of radial coordinates suggested by Krori and Barua [82]. They proposed the specific forms of these functions as a(r) = Br 2 + C and b(r) = Ar 2 , where A, B and C are the three constant numbers. One can find the values of this triplet (A, B, C) by considering some appropriate boundary conditions. Then, Eqs.(2.12)-(2.14) yield (2.17)

Boundary Conditions
We consider a timelike hypersurface denoted by Ω that has differentiated interior manifold given in Eq.(2.7) and outer geometry. The exterior region is through the vacuum Schwarzschild spacetime as where M is the gravitating matter content. At the boundary surface Ω, the continuity of the metric variables g tt , g rr and the ∂g tt /∂r provide the following constraints For the smooth and continuous matching of the manifolds, the above constraints at Ω should be fulfilled. We now pick some observational values of A, B and C from the literature to check the construction as well as the stability of the compact stellar bodies. We consider, three stellar toy models Her X-1, SAXJ1808.4-3658 and 4U1820-30 labeled with CS1, CS2 and CS3 here. The masses of these star candidates are 0.88, 1.435 and 2.25 solar masses, respectively. All of these stellar structures satisfy Buchdahl Bondi bound as their 2M/R ratios are less than 8/9. We further suppose that the junction conditions at the stellar core are [84] where ρ c is the critical mass density.

Various f (R, T ) Models
The aim of this section is to study some of the notable f (R, T ) models. We want to analyze the influences of these f (R, T ) models in the construction and stability of the compact stellar candidates in the background of some observations. We shall discuss some realistic features of spherical systems like, compactness, stability, evolution of energy density with radial coordinate, the measurement of anisotropic pressure evolution and the different energy conditions. The results obtained such investigations may provide some hidden realities corresponding to the both theoretical and astrophysical regimes. Equation (2.11) provides where i = 1, 2. There has been very interesting results which reveal that the inclusion of extra dark energy/DM components mediated from alternative theories could bring some exciting results. For instance, collapse time, existence of more compact structure, stability, phenomenon of core formation and above all can be well influenced by these dark source terms unlike GR. The reexamination of GR problems in alternatives theories may be helpful to shed some light on the models viability and their usefulness on physical grounds. The cross-examinations may present both quantitative and qualitative different consequences than that of GR. It is of particular interest for many researchers to explore the predictions as well as the outcomes of modified gravitational theories like f (R, T ) theory concerning the existence of stellar structure and their stability. The f (R, T ) gravity can be treated as a mathematical tool to examine various unknown features of gravitational dynamics at large scales. Schaefer and Koyama [88] generalized their study of gravitational collapse of spherical structures in the realm of modified gravity and Birkhoff-theorem and found enhanced cluster merging rates as well as overdensed populations of relativistic structures (due to modified gravity). Capozziello et al. [89] computed modified Lané-Emden expression with metric f (R) corrections and Newtonian approximation. They get some exceptional results of density and pressure distribution in the analysis of the hydrostatic scenario of few celestial bodies. Cembranos et al. [90] studied stellar structure formation in the background of modified gravity and found comparatively higher contraction in the collapsing rate of spherical systems at its initial stages unlike GR.
Astashenok et al. [91] examined the impact of few modified gravity formulations on the existence of compact structure and inferred that there exists some models whose corrections could lead arena of having relatively more compact stars than in GR. Yousaf and Bhatti [77] studied the influences of modified dark source terms on mass-radius relationships for compact stars and concluded that some extended gravities mediated by Lagrangians αR 2 and αR 2 −βR 1+γR are likely to host supermassive relativistic systems with comparatively smaller radii than in GR. Resco et al. [92] numerically calculated apparent masses of neutron star models in modified gravity and inferred that generically some modified extra degrees of freedom likely to keep significantly massive neutron star with smaller radii than in GR. Such type of investigations could provide theoretical well-consistent way to handle and study classes of massive and super massive structures at large scales. Bamba et al. [93] claimed that modified gravitational dynamics could provide an additional platform in few regions of spacetime that could lead to stable configurations of relativistic system. In this paper, we found relatively stellar bodies with much higher densities than that found in GR [80] as seen from Fig.(1).
In the following we discuss two different f (R, T ) models.

Model 1
Firstly, we take the tanh modification of the Ricci scalar function that was proposed by [85]. In this respect, the f (R, T ) gravity model (4.1) give where α ∈ R + andR ∈ R + in which R + denotes the set of positive real numbers. On setting α = 0, the dynamics of Einsteins gravity can be discussed.

Model 2
Next, we consider another viable formulations of f (R) gravity that was first suggested by [86]. Under this context, Eq.(4.1) becomes whereR, γ and q are the free non zero and non negative parameters.
In the coming section, we will use these two viable model for further investigation of compact stars.

Some Physical Properties
This section is devoted to analyze various aspects of some compact stellar toy models. We assume three different star distributions, i.e., CS1, CS2, CS3 as well as two observationally consistent f (R, T ) models. Upon substituting this data in Eqs.(2.8)-(2.10), we get values of matter variables in the form of four parameters, α,R, q and ǫ. Then, we study some physical features to obtain the realistic configurations of compact stellar structures (shown in Table 1). The comparison between our theoretical outcomes and the observational data may provide the strong evidences for f (R, T ) models. We shall show our results with the help of plots.

Density and Pressure Evolutions
This section is devoted to analyze the matter parameters of all the three stellar bodies whose variations are depending on the radial coordinate. We check the radial variations in the anisotropic pressure, energy density and their corresponding radial derivatives. Now with the help of Eqs.(2.8)-(2.10), we plot diagrams (1) for all the three stellar models in the background of two above mentioned f (R, T ) models¿ We infer that the behavior of energy densities keep on increasing till the constraint r → 0. This indicates the distribution of ρ to be increasing with respect to the decreasing choices of r. One can say from these that our stellar toy models are of having highly compact cores in the degrees of freedom coming from Eqs. Furthermore, Figures (4), (5) and (6) are describing the variations in the r derivative of ρ and p i . One can observe from these that dρ dr < 0, dpr dr < 0 and dpt dr < 0 for the given two viable models and the three toy models. Under the constraint r = 0, one can notice that the r variations of all the matter variables vanishes, such that dρ dr = 0, dp r dr = 0.
Furthermore, the twice variations of these structural quantities have been found to be negative. These consequences are providing a seed for the high compact profiles of such star structures around their subsequent cores, thus indicating that the configurations of compact and dense stars do exists in the arena of f (R, T ) gravity.

Energy Conditions
For a realistic matter content distribution, one need to check some particular conditions, such a mathematical conditions are called energy conditions (ECs). These mathematical conditions are the coordinate-invariant. In the regime of modified f (R, T ) gravity like the effective forms of density and the anisotropic pressure, the null energy Condition (NEC) and the weak energy conditions (WEC) are devised, respectively, as while the rest of ECs, i.e., strong energy condition (SEC) and dominant energy condition (DEC) give • SEC ⇔ ρ eff + 2P eff t + 3P eff r ≥ 0 and ρ eff + P eff i ≥ 0, • DEC ⇔ ρ eff ≥ 0 and ρ eff ± P eff i ≥ 0.   It has been easily be figure out from the plots shown in Figures (7) and (8) that our under considerations of both models with CS1, CS2 and CS3 obey all ECs. This suggests that anisotropic fluid content (2.4) describes the realistic source of gravitational effects.

TOV Equation
A mathematical expression that provide some restrictions on the geometrical distributions of the relativistic spherical system, which is normally hydrostatic equilibrium and is often called as TOV equation. For the spherical line element, this equation is written as dP r dr + a ′ (ρ + P r ) 2 + 2(P r − P t ) r = 0, (5.1) The above equation contains three well-known interactions, i.e., the gravitational, (F g ), hydrostatic, (F h ) and anisotropic, (F a ) forces. Therefore, one may write above equation as The values of these forces by using of Eqs.(2.8)-(2.10) giives [4{1 + 2ǫ + B 2 r 4 ǫ + A 2 Br 6 ǫ − exp(Ar 2 )(1 + 2ǫ) + Ar 2 Figure 4: Graph describing dρ/dr with respect to r for both f (R) models. Figure 5: Graph describing dp r /dr with increasing r for both f (R) models. Figure 6: Graph describing dp t /dr with increasing r for both f (R) models.
Using above relations together with observational values of A, B and C from Table 1, we have drawn some diagrams mentioned in Figure (9). In Figure (9), the left plot is for f (R, T ) model (4.2) and right plot is for model (4.3). The role of these three types of interactions with respect to r (km) has been seen in this diagram in the modeling of compact bodies CS1, CS2 and CS3.

Stability Analysis
This section illustrates the stability of the compact relativistic structures. It would be interesting to mention that, for any arbitrary relativistic observer, only those stellar bodies are worthy to study which remain in stable position after the application of fluctuations. Thus, the stability problem is among the burning issues for the relativistic astrophysicists.

The Measurement of Anisotropy
The extent of anisotropicity as well as its magnitude within a local anisotropic fluid distributions can be defined as Using Eqs.(2.8)-(2.10), we can write ∆ as The above equation contains f i terms. We use the two different viable f i models from Eqs.(4.2)-(4.3) and get the three different equations, each incorporating a particular f (R, T ) model corrections. After using observational results of the compact stars given in Table 1, we terminated with six set of equations. We check the behavior of anisotropy in the degrees of anisotropicity in the background of relativistic compact objects, by plotting these equations. It can be seen from the plot shown in Figure (13) that the ∆ is remained greater than zero for all the three relativistic compact stellar structures which resulting that the influence of radial pressure, P r is lesser than that of transverse pressure, P t . This outward directed measure of anisotropy for these two viable models.

Summary
As the universe is accelerating, this ground reality attracts researchers towards the extension of GR and modified theories of gravity. The first modification to GR gives birth to f (R) gravity, which provides a good attention while the extension of this modification gives rise to f (R, T ) gravity and attain much attentions because of some quantum effects arising in the theory. These modifications appear in the field of low energy action for effective quantum gravity theory.
In this work, we study the anisotropic relativistic compact stellar objects with static and spherical structures. For this purpose, we consider the relativistic compact stellar objects whose interior geometry is based on anisotropy in f (R, T ) gravity. Krori and Barua [82] proposed the relativistic interiors of stellar bodies through specific metric combnations. With these techniques, we connect this interior geometry with the exterior Schwarzschild geometry and find the constant (in term of mass and radius) of interior metric over the boundary. After this, we used some observational data from which we find the numerical values of these constants. Then we put these numerical values in our calculations and plots our results.
In this paper, we have used BK solution, according to which a(r) = Br 2 + C and b(r) = Ar 2 , where A, B and C are the three constant numbers which can be evaluated depending on several physical requirements. Such solutions [82] are claimed to be singularityfree for the static spherical systems in the background of GR. Furthermore, this solution is asserted to be regular everywhere and the matter variables, like mass density, pressure etc are finite all over within the relativistic spherical system.
It is of our interest to check the outcomes and the predictions of one of the extended gravities, i.e., f (R, T ) theory regarding the existence and stability of spherical stars. Therefore, we have explored the modified version of the TOV equation in the realm of couple of f (R, T ) models. Figure 9 states that for the particular choices of f (R, T ) model parameters, there exists some eras under which system could attain equilibrium condition by keeping all the forces sum to be zero. It is worthy to mention that usual GR forces are being modified due to f (R, T ) models, thus producing some extra effects in the forces F g , F h and F a . The stability of our compact stars not only depends upon KB-solution but also on the choices of parameters involved in f (R, T ) models mentioned in Eqs.(4.2) and (4.3). The difference of the squares of sound speeds, i.e, v 2 st − v 2 sr has been found to be within [0, 1], thus describing our relativistic stars to be in stable window with certain values of parameters involved in the corresponding f (R, T ) models.
The energy density remains positive and maximum at the core of compact stars. Energy conditions holds for all these three compact stars and the radial as well as the transverse equation of state parameter are in a usual range i.g. 0 < ω i < 1, here i = r, t. This indicate that the interior structures of these relativistic compacts stellar objects are composed of normal ordinary matter. For all case, we find the anisotropy directed outward, e.g., p t > p r or ∆ > 0. Similarly, the transverse as well as the radial sound speed for all these compact stars are in stable limits that implies the stability of these spherical anisotropic compact bodies in the realm of f (R, T ) gravity.