Realistic stellar anisotropic model satisfying Karmarker condition in f(R, T) gravity

The present study explores the f(R,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathcal {R},\mathcal {T})$$\end{document} modified gravity on the basis of observational data for three different compact stars with matter profile as anisotropic fluid without electric charge. In this respect, we adopt the well-known Karmarker condition and assume a specific and interesting model for grr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {g}_{rr}$$\end{document} metric potential component which is compatible with this condition. This choice further leads to a viable form of metric component gtt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {g}_{tt}$$\end{document} by utilizing the Karmarkar condition. We also present the interior geometry in the reference of Schwarzschild interior and Kohler–Chao cosmological like solutions for f(R,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathcal {R},\mathcal {T})$$\end{document} theory. Moreover, we calculate the spacetime constants by using the masses and radii from the observational data of three different compact stars namely 4U 1538-52, LMC X-4 and PSR J1614-2230. In order to explore the viability and stability of the obtained solutions, some physical parameters and properties are presented graphically for all three different compact object models. It is noticed that the parameters c and λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda $$\end{document} have some important and considerable role for these solutions. It is concluded that our obtained solutions are physically acceptable, bearing a well-behave nature in f(R,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\mathcal {R},\mathcal {T})$$\end{document} modified gravity.


Introduction
In last few decades, various extended theories of gravity have been proposed as the most promising candidates of mysterious dark energy for exploring the accelerated expansion aspects of our cosmos. There are several independent observations in the context of astrophysics that provide evidences about the accelerated expanding nature of space. These astroa e-mail: gmustafa3828@gmail.com b e-mails: mzubairkk@gmail.com; drmzubair@cuilahore.edu.pk c e-mail: swaheed@pmu.edu.sa d e-mail: xiatc@shu.edu.cn physical experiments include the outcomes acquired from the supernova type Ia (SNIa) [1][2][3][4][5][6][7][8], large scale construction surveys [9], X-ray brightness from galaxy erect [10], cosmic microwave background radiation (CMBR) [11][12][13], weak lensing and the baryon acoustic oscillation (BAO) surveys [14]. This phenomenon of accelerating cosmic expansion is regarded as an outstanding critical riddles of contemporary physics. It is argued that the expansion of cosmos is accelerated due to the presence of some uncertain dominant source of energy labeled as dark energy (DE). In order to incorporate this unusual motive of dynamical cosmos, some modifications in the Lagrangian density of Einstein's general relativity have been presented in the literature. One of the way to deal with the problem is to modify the matter profile of the density by adding some DE terms as scalar field or cosmological constant, k-essence, canonical kinetic scalar term, quintessence, and different versions of chaplying gas etc. [15][16][17][18][19][20]. While in the alternative approach, researchers extended the gravitational part by adding some extra degrees of freedom there, which provided the group of extended theories of gravity. In this respect, some leading examples include f (T, T G ) theory, f (G) gravity, f (R) framework, braneworld scenarios, Kalb-Ramond background, Gauss-Bonnet gravity and Brans-Dicke gravity theories etc. [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36]. Based on these modifications, many other remarkable generalized modifications are also available, which help to explore various cosmic aspects successfully. These all modifications have passed different astrophysical and necessary solar system constraints and are regarded as viable candidates.
After the formulation of general theory of relativity (GR), the f (R) theory is regarded as one of the most interesting and viable extension of GR. Later on, its different modified versions have been presented by researchers that are also proved as successful in various respects. Particularly, its recent extension namely f (R, T ) gravity, proposed by Harko et al. [37] almost 6 years ago, has attracted many researchers. This modification is based on a governing function f (R, T ) depending on the Ricci scalar R and the trace T of anisotropic like energy-momentum tensor. They also derived the corresponding field equation by utilizing the metric potentials formalism and discussed the significance of this modification. Further, the authors have presented different models for f (R, T ) function in separable form, i.e., f (R, T ) = f 1 (R) + f 2 (T ). It is considered as a very interesting modification because the resulting field equations have not much difficult form or unhandled order. There is a sufficient literature available where numerous cosmological aspects of this theory have been explored like expansion of universe due to accelerated matter, Birkhoff's theorem, scalar field reconstructions, stability using cosmological perturbation, large scale structure, thermodynamical laws and its relevant features, constrains regarding solar system, stellar equilibrium configurations of compact stars, neutron stars, gravitational collapse phenomenon [38][39][40][41][42][43][44][45][46][47][48].
The study of compact stars formation and phenomenon of gravitational collapse is regarded as one of the most fascinating subjects in modern cosmology and astronomy. In 1916, Carl Schwarzschild presented the exact solution of interior of symmetric star spherically by using the uniform density based matter profile [49]. In 1939, Oppenheimer and Snyder [50] contributed to explore the gravitational collapse with homogeneity based dust sphere. In literature [51][52][53][54][55], hundreds of static analytic models representing the relativistic stars have been constructed by introducing bulk viscous effects, anisotropic pressures, charge, multilayered fluids and equations of state etc. The discussion of ultra densities for matter profile of order 10 15 g cm −3 provides a new idea that radial pressure profile and tangential pressure profile are unequal. In this context, Ruderman [56] was the pioneer who predicted that in the high density, relativistic interaction of nuclear matter yields anisotropy as an inherent feature. Lobo [57] used barotropic EoS for exploring some compact object models which is defined as p = ρ; − 1 < < −1/3, for the dark stars existence. He also presented a modified version of Mazur-Mottola gravastar model by utilizing the junction conditions to match Schwarzschild vacuum solution with static line element. Egeland [58] explained Neutron stars by taking mass-radius relationship into account. In another study [59], Mak and Harko proposed a spherically symmetric model exhibiting the characteristics of strange stars. Further, Rahaman et al. [60] discussed the possibility of compact stars formation by using Krori-Barua model in the presence of Chaplygin gas matter distribution. In modified gravitational frameworks like f(R) gravity and scalartensor theories, much work has already been done for modeling the massive as well as neutron stars [61][62][63][64][65]. Hossein et al. [66] have developed some compact stars models with anisotropic matter and variable cosmological constant where they assumed a linear equation of state. In different collaborations [67,68], Herrera discussed some interesting anisotropic solutions for static as well as non static sources. In the context of GR and modified gravitational frameworks of f(T) and f(G) theories, much work has already been done by the researchers [69][70][71][72][73][74][75][76]. Zubair and Abbas [77] have explored some interior compact star models in f(R) extension of GR by including Krori and Barua solutions and discussed different physical features as well as the stability of the obtained models.
Naidu and Govender [78] investigated two new stellar models with same radii and mass but distinct pressure distributions. In another study, Zubair et al. [51][52][53][54][55] investigated the possible existence of compact stars in f(R,T) framework by taking the analytic models of Krori and Barua line element into account and presented a detailed analysis of the models by discussing some physical features. In 2013, Herrera and Barreto [79] acquired a mathematical procedural technique to produce a relativistic-polytropes with radial and tangential pressure profiles by using curvature coordinates. Recently, a considerable approach for deriving the analytic solutions of Einstein field equations, representing the compact objects, has been proposed by the researchers namely the Karmarkar condition [80][81][82][83][84][85][86][87]. This condition was firstly proposed by Karmarkar [88] and is regarded as a compulsory condition for a spherically symmetric spacetime to be of embedding class-I. It is basically a mathematical tool which helps us in obtaining the exact solutions of field equations. In literature [89][90][91][92][93], this condition has been used by numerous researchers for discussing the compact stars models. In the present paper, we shall adopt the Karmarkar condition to develop the analytic solutions representing compact objects in f (R, T ) gravity.
Being motivated from the above literature, in the present manuscript, we will discuss the formation of compact stars in f (R, T ) modified gravity by considering three different models of compact stars. For this purpose, we will consider static spherically geometry filled with anisotropic matter contents and also take the well-known Karmarkar condition. In the up-coming section, we will define the mathematical structure of f (R, T ) gravity and formulate its field equations. In the same section, we will present the Karmarkar condition briefly and calculate isotropic like Class-I solutions with two different cases of f (R, T ) gravity. In Sect. 3, we will discuss the set up for a new family of embedding class-I solutions along with the physical boundary conditions. In Sect. 4, we will present the analysis of obtained solutions using three different models for compact stars like 4U 1538-52, LMC X-4 and PSR J1614-2230 by investigating different physical properties analytically and graphically. In the same section, we will discuss the stability of obtained model using dif-ferent viable measures. In last, we will conclude the whole discussion by focusing on the main achievements.

f (R, T ) gravity field equations and anisotropic fluid distribution
In this segment, we will describe the basic formulation of f (R, T ) modified theory of gravity and develop its field equations for a static spherically symmetric spacetime geometry. Here we will also briefly explain the idea of Karmarker condition and some assumptions taken for this work. The modified Einstein-Hilbert action defining the well-known f (R, T ) gravitational framework is given by [37] where L m denotes the Lagrangian density of ordinary matter and f (R, T ) represents a generic function depending on two variables namely the scalar curvature R and the energy-momentum tensor trace T μν , which is defined as T = T μν T μν . The variation of the above modified action with respect to metric tensor leads to the following set of field equations: 8π In the above set of equations, the notation ∇ refers to the covariant derivative, while the symbol denotes the d'Alembert operator, i.e., = ∂ μ ( √ −gg μν ∂ ν )/ √ −g. Moreover, the function f (R, T ) with subscripts R and T represent the respective derivatives, i.e., f R (R, T ) = ∂ f (R, T )/∂R and f T (R, T ) = ∂ f (R, T )/∂T . Also, the term μν is given by The energy-momentum tensor representing the anisotropic matter profile within a star is defined by where ζ μ and ξ μ represent the 4-velocity vectors defined by the relations ζ μ = e −a δ μ 0 and ξ μ = e −b δ μ 1 and also satisfy the condition ζ μ ζ μ = −ξ μ ξ μ = 1. Here, the terms ρ, p t and p r all are radial coordinate dependent functions and represent the energy density, radial and tangential components of pressure, respectively. It is also interesting to mention here that the function p t is orthogonal in the direction of ξ ν , while p r is orthogonal in the direction of ξ μ . By taking the varia-tion of indices along with energy-momentum tensor (3), the set of field equations (2) take the following form The line element representing the geometry of a static object exhibiting spherically symmetry is given by where d 2 = (dθ 2 +sin 2 θ dφ 2 ) and defines the g θθ and g φφ components. For the simplicity in calculations, we choose an interesting simple model of f (R, T ) given by [38][39][40][41][42][43][44][45][46][47][48] where λ is an arbitrary constant. Using metric (5) in the Eq.(4) and then by re-arranging the resulting field equations for the functions ρ, p r and p t , we obtain Here the last equation for 8π( p t − p r ) gives a measurement of anisotropy which is denoted by , i.e., = 8π( p t − p r ). For having isotropic pressures condition, this equation should be vanished, = 0. Now we will explain the well-known Karmarker condition briefly which is the most important tool considered in the present study. The foundation of Karmarkar condition depends upon the Riemannian space of class-I. Eisenhart [80] presented a necessary and sufficient condition which is based on a symmetric tensor of second order β μη and the Riemann curvature tensor and is defined in terms of the following equations: where ; denotes the covariant derivative and = ±1. These values refer to a time-like or space-like manifold depending on the sign taken as − or +, respectively. Now, by utilizing the non-zero Riemann curvature components for the spacetime geometry (5) along with the nonzero components of second order of symmetric tensor β μη in Eq.(11), we get a relationship of the form This relation represents the well-known Karmarkar condition where R 2323 = 0. Here the expression R 2323 = 0 represents the Pandey-sharma condition [94]. Pandey and Sharma presented an argument about the Karmarkar condition that this condition is not a sufficient condition for a symmetric space time for being a class-I model. Further, by plugging all the respective values of the Riemann curvature components in Eq.(13), we get a differential equation, given by where e ν(r ) = 1. The integration of Eq.(14) leads to a relationship between g rr and g tt components of the line element given by where K is introduced as a constant of integration. Now, we calculate the derivative of ν(r ) with respect to radial coordinate r , which can be written as where B 2 = 1/K . Using these value of e ν(r ) and ν (r ) in Eq. (10), the anisotropy function for f (R, T ) takes the following form where λ = −1. In the upcoming sections, we will talk about the solutions representing stellar models for isotropic as well as anisotropic cases respectively.

Isotropic like class-I solutions for f (R, T ) theory
In this segment, we shall discuss the case of isotropic pressures and its consequent outcomes. The isotropic case, where p r = p t , yields the anisotropy measure function as zero, i.e., = 0. Assuming μ (r ) = 0, Eq. (17) provides two interesting results given by or or the product of both expressions equal to zero. Here Eq. (18) corresponds to Schwarzschild's constant like density model, which is already presented by Bhar et al. [95], while the second result given by Eq. (19) leads to the Kohler-Chao cosmological like solution.

Schwarzschild interior solution for f (R, T ) theory
In order to find the Schwarzschild interior solution [49], we integrate Eq. (18), which yields the following solution: where c is a constant of integration. Now, by utilizing Eq. (20) in the Eq. (15), we get a differential equation given by By solving this differential equation for μ(r ), we obtain a final solution for the g tt component of the line element given by where A is another constant of integration. The above solution (22) is a well-known Schwarzschild interior solution model that describes an incompressible static sphere with uniform energy density. Now, we present the expressions for energy density and both pressure components, which are cal-culated in Eqs. (7) and (8). In the present case, these expressions take the following form: This Schwarzschild solution is very interesting as well as simple solution and have already been studied by several researchers in literature.

Kohler-Chao cosmological like solution for f (R, T ) theory
Here we discuss the Kohler-Chao cosmological model like solution [96]. For this purpose, we integrate Eq. (19) which leads to the following solution for the reciprocal of g rr component of the line element given as follows where A denotes a constant of integration. Now, by using the value from Eq. (26) of e μ in the Eq. (15), we acquire the following solution for reciprocal of g tt component: The above solution (27) corresponds to the Kohler-Chao cosmological model like solution. Also, the corresponding expressions for energy density and pressure components take the following forms: Here the parameter B is assumed to be a positive quantity.

A new setup for a family of embedding class-I models
In this section, we shall study new stellar models which are based on anisotropic fluid and are compatible with the Karmarker condition. It is a considerable point that the pressure anisotropy, i.e., has an important role in the discussion of gravitational collapse process. Recently, Naidu and Govender [78] studied the collapse dynamics which is related to its radial pressure and density function of radial coordinate under the stellar fluid. Further, they have assumed a linear form of equation of state within a static configuration which can be presented by the expression p r = ρ − ε, with and ε as constants. They argued that the succeeding collapse is responsive to the correlation of energy density function and radial pressure. Further, they discussed the impact of parameter on the temperature description of collapsing configuration.
For having a family of solutions compatible with the Karmarker condition and involve anisotropic matter contents, we assume a specific model for the component g rr which is acceptable for all the gravitational aspects and is given by where a, c and b are parameters which can be calculated from some physical boundary conditions. This models has been extensively used and explored in many directions related to astrophysics and cosmology. In 1999, Raychaudhuri and Dadhich discussed its significance without utilizing the Big Bang singularity to find the oscillating cosmological model. Another important characteristic of this model is that it permits the divination of blue shifts beyond violating the GR basic aspects. By utilizing the Eq. (31) in Eq. (15), we get a differential equation of the following form: Solving Eq. (32) for the expression of μ(r ), we get a relationship for the metric component g tt as follows Further, we use the values of g tt and g rr components in the field Eqs. (7)-(10), we get the expressions of the following quantities: +a Bλ sin 2 br 2 + c +B cos br 2 + c a 2 (2λ + 1) a Ab(2λ + 1) sin br 2 + c +a Bλ sin 2 br 2 + c +bB −ar 2 cos 2 br 2 + c + ar 2 + 2λ + 4 +2abB(2λ + 1)r 2 sin br 2 + c cos 2 br 2 + c , = √ ar 2 a sin 3 br 2 + c − 2b cos br 2 + c −4 √ a Ab sin br 2 + c + a B sin 2 br 2 + c + 4bB 16π(λ + 1) ar 2 sin 2 br 2 + c + 1 2 √ a B cos br 2 + c − 2 Ab .

Physical boundary conditions
In this portion, we shall describe some physically acceptable boundary conditions by matching interior solution to the exterior solution. Here we consider the exterior solution as Schwarzschild space-time which is defined as where Here we have enforced a condition on the radial coordinate given by r > r ς , where r ς is termed as Schwarzschild radius. It is interesting to mention here that without enforcing this restriction, we may get some black hole solutions. Now, at the boundary r = R, the junction conditions for the continuity of metric components g tt , g rr and ∂g tt ∂r are given by where − and + correspond to the interior and exterior solutions, respectively. By plugging all the corresponding values, we obtain the relations Further, by utilizing these boundaries conditions (40) and (42), we get the following relations for unknowns:  where the parameters b and c are treated as free parameters. Moreover, M and R represent the mass and radius of the object, respectively. Here we consider three different models for the compact stars namely 4U 1538-52, LMC X-4 and PSR J1614-2230 for the discussion of the new obtained stellar model of the present study. Since in the present study, we have used a linear model of generic function given by f (R, T ) = R + λT which demonstrates the coupling of the Ricci scalar and the trace of the energy-momentum tensor, so it is very important to explore the impact of parameter λ on the obtained solutions. Here we shall investigate the impact of this parameter by taking few small and large values of parameter λ into account. For graphical analysis and further calculations, we have considered its three small values from the range 0 < λ < 1, i.e., λ = 0.35, λ = 0.45, and λ = 0.55 and three large values such that λ > 1, i.e., λ = 2.00, λ = 4.00, and λ = 6.00. In the first place, we have considered its small values and the calculated values of different unknowns against these λ values are provided in Table 1. Further in Table 2, we have provided the values of theses unknowns against three assumed slightly large values of λ. From both tables, it can be observed that the values of parameters A and B increase with the increasing values of λ.

Analysis of the physical properties of the f (R, T ) stellar model
In this section , we explore our results in more detail by focusing on some physical aspects and necessary properties of the obtained f (R, T ) stellar configuration using three different stellar objects observation data values. For this purpose, we shall present the discussions about some physical measures analytically and graphically by taking different values of parameter λ. Here we shall consider three different models for stellar objects like 4U 1538-52, LMC X-4 and PSR J1614-2230.

Evolution of metric functions, energy density and pressure components
Here we first discuss the behavior of metric function in the stellar objects, i.e., g tt = e μ and g rr = e ν . It can be noticed from the Eqs. (31) and (33) that e ν(r =0) = 1 and e μ(r =0) = 0 which indicates that this model is physically realistic and acceptable. Moreover, from the left penal of Figs. 1 and 4, it is observed that the metric potential components are regular and monotonically increasing in the inner side of the stellar object. Here the green graphs show the behavior of metric components for 4U 1538-52 model, the blue graphs represent the metric components for LMC X-4 model while the red graphs describe the metric components for PSR J1614-2230 model. Now we calculate the energy density and pressure components at the central values r = 0, i.e., ρ(r = 0), p r (r = 0) and p r (r = 0) which are given as follows It can be easily observed from the Eqs. (46) and (47) that the ratio p rc /ρ c is consistent with the Zeldovich's condition. In the present scenario, the Zeldovich's condition for f (R, T ) theory can be written as We have calculated the values of different physical parameters like metric potentials, energy density and pressures at the center for both small and large choices of λ. The obtained results are summarized in form of Tables 3 and 4. From  Table 3, it can be observed that the metric components satisfy the important condition at the center, i.e., e ν(r =0) = 1 and e μ(r =0) = 0, for all three small λ choices. We have also calculated that the values of energy density and pressure components at the centre, i.e., ρ(r = 0), p r (r = 0) and p r (r = 0) which indicate that the values of energy density, pressure components at the center, energy density at the boundary, and the ratio p r c /ρ 0 = p t c /ρ 0 (Zeldovich's condition) gradually decreasing as the λ increases.
This behavior can also be observed from the Figs. 1, 2, 3, 4, 5 and 6. Figure 1 shows that the energy density function exhibits positive and gradually decreasing behavior for small increasing values of λ for 4U 1538-52, LMC X-4, and PSR J1614-2230 stellar models. From Fig. 2, it can be seen that both pressure components are also gradually decreasing against increasing small values of λ. The existence of nonzero anisotropy is considered as a significant aspect in the relativistic stellar configurations without an electric charge. In the present case, the non-zero anisotropy played a special role for exploring the compact stars solutions. Now we describe the graphical development of anisotropy function, i.e., (r ) and the derivatives of ρ, p r and p t with respect to radial coordinate r (gradients) in the current scenario. For 4U 1538-52 model, the anisotropy function exhibits positive behavior for increasing λ through small values. For LMC X-4 model, the anisotropy function is observed positive for λ = 0.35, and λ = 0.45 while it exhibits negative behavior for λ = 0.55. Further for PSR J1614-2230 model, it shows positive behavior for λ = 0.35 only. It is seen that in contrast to small λ values, the anisotropy function exhibits negative behavior for all large choices of this parameter. Clearly, in such cases where the anisotropy function remains positive with increasing behavior, it supports to the condition p t > p r . Further, it Table 3 Calculated values of different physical properties at center and boundary for different three values of λ p t p t < 0 (near to boundary) p t < 0 (near to boundary) p t < 0 (near to boundary) (r ) (r ) < 0 (throughout) (r ) < 0 (throughout) (r ) < 0 (throughout)

Gradients
Gradients< 0 G r a d i e n t s < 0 G r a d i e n t s < 0   Table 5. From Table 5, it can be concluded that  Table 2 under the condition of non-zero anisotropy, our obtained solutions are stable only when λ = 0.35. It is evident from the graphs that all derivatives exhibit negative and decreasing behavior and hence satisfy the necessary conditions given by dρ/dr < 0, dp r /dr < 0 and dp t /dr < 0. It is also observed that these derivatives are vanished at the central radius, i.e., dρ(r = 0)/dr < 0, dp r (r = 0)/dr < 0 and dp t (r = 0)/dr < 0. Finally, we can conclude that all these gradients satisfy the required conditions in f (R, T ) theory. Likewise, a similar behavior can also be noticed from the Figs. 4, 5 and 6 where the large values of λ has been taken into account. In this case, it can be concluded that these physical parameters exhibit a quite similar behavior but their values are decreasing more rapidly as compared to the small values of λ. Mathematically, these results are also summarized in Table 4.

Equilibrium condition
In this part, we shall investigate the stability of obtained solutions by utilizing the equilibrium condition for three different models of stellar system in f (R, T ) gravity. For this purpose, we take the well-known Tolman-Oppenheimer-Volkov (TOV) equation [38][39][40][41][42][43][44][45][46][47][48] given as follows The above Eq. (49) defines the equilibrium condition for a configuration by taking three forces into account namely hydrostatic, gravitational and anisotropic forces. The following are the mathematical expressions for these forces and consequently, Eq.(50) can be rewritten in the form given by The graphical illustration of these forces are provided in the left penal of Figs. 7 and 8 for increasing small and large values of λ. Here the green, blue and red curves provide the forces for 4U 1538-52, LMC X-4 and PSR J1614-2230 model, respectively. From these graphs, it can be noticed that the gravitational forces are displayed below the x-axis, the anisotropic forces are observed almost on the x-axis, while the curves representing hydrostatic forces are above the xaxis for all three stellar models with small increasing values of λ. So it can be concluded that these forces almost balance each other's effect and hence leave the configuration stable with necessary mass and radii as described in Tables 1 and 2 in f (R, T ) theory of gravity

Energy conditions
Energy condition bounds have many cosmological applications in the context of GR and modified gravity theories. In general, these conditions are defined in terms of four constraints namely DEC, S EC, W EC and N EC and are given as Here we shall discuss the behavior of energy constraints namely N EC, W EC, S EC and DEC versus radial coordinate which is provided in the right penal of the Figs. 7 and 8 where we have considered the small as well as large values of λ. In these Figures, the green, blue and red curves correspond to the energy conditions for 4U 1538-52, LMC X-4 and PSR J1614-2230 models, respectively. It is seen from these graphs that these constraints are satisfied for the present model as these constraints show positive but gradually decreasing behavior against the increasing values of parameter λ. This is considered as one of the necessary Fig. 7 Illustrates the graphical behavior of forces namely F a , F h and F g and energy conditions versus r for 4U 1538-52, LMC X-4 and PSR J1614-2230 along with the values of parameters from Table 1.
Here the solid, dashed and small dashed curves refer to the choices λ = 0.35, λ = 0.45 and λ = 0.55, respectively  Tables 1 and 2. Here the impact of λ on the obtained solutions is if we increase its value then these conditions decrease more rapidly but remain satisfied. It is interesting to mention here that the behavior of other constraints have already been presented in the Figs. 1, 2, 3, 4, 5 and 6.

Mass function, compactness parameter and gravitational red-shift function
Here we will discuss the relation of mass and radius of a compact star, which is defined by the following equation where m represents the mass function depending on the radial coordinate. Also, the compactness parameter of a compact star is defined by the following expression From Eq. (43) the gravitational red-shift of a stellar system is given by Here we shall explore the impact of parameter λ on the obtained solutions by discussing the behavior of corresponding mass function, compactness parameter and gravitational red-shift function versus radial coordinate graphically. It can be observed from the Figs. 9, 10, 11 and 12 that mass function and compactness parameter exhibit a regular increasing behavior with the increasing values of small or large λ choices while the gravitational red-shift function indicates the decreasing behavior in both cases of parameter λ. Here it can also be seen that the λ has an impact on the obtained solutions, i.e., if we consider large λ values, then function are showing gradually increasing or decreasing behavior as compared to its small values. These behavior show that all these physical parameters satisfy the required condition for compact stars in f (R, T ) theory.  Table 2 with λ = 2.00 (solid), λ = 4.00 (dashed) and λ = 6.00 (small dashed) In this part, we formulate the equations of state parameter for radial and tangential pressures denoted by w r and w t , respectively. These parameters are defined as follows Here we discuss the graphical development of two different ratios, i.e., w r and w t versus radial coordinate. It is observed from the right penal of Figs. 11 and 12, both ratios p r ρ and p t ρ exhibit monotonically decreasing behavior less than 1 against different small and large values of parameter λ which is regraded as one of the necessary condition for  Table 2. Here λ = 2.00 (solid), λ = 4.00 (dashed) and λ = 6.00 (small dashed) Table 6 The development of Herrera cracking criterion, i.e., −1 ≤ v 2 t − v 2 r ≤ 0 for different three values of λ, under three different well-known compact stars 4U 1538-52, LMC X-4 and PSR J1614-2230  Table 2 where λ = 2.00 (solid), λ = 4.00 (dashed) and λ = 6.00 (small dashed)

Conclusion
In the present paper, we have described a new set up for a new family of embedding class-I models in f (R, T ) theory by using the well-known Karmarkar condition along with the Pandey-Sharma condition for three different compact stars namely 4U 1538-52, LMC X-4 and PSR J1614-2230. Here, for exploring physical validity of the obtained solutions, we  Table 2. Here λ = 2.00 (solid), λ = 4.00 (dashed) and λ = 6.00 (small dashed) have discussed all the compact stars related physical parameters analytically and graphically. For the present study, we have employed the matter contents as anisotropic fluid without any net electric charge for a spherically symmetric spacetime. The physical behavior of all relevant parameters have been presented in the Figs. 1, 2, 3, 4, 5, 6, 7 and 8 for the considered stellar models. All achieved results can be summarized as follows Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: There is no observational data related to this article. The necessary calculations and graphic discussion can be made available on request.] 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 .