An EGD model in the background of embedding class I space–time

This work is devoted to the study of relativistic anisotropic compact objects. To obtain this class of solutions of the Einstein field equations, we have developed a general scheme to generate the metric of the space–time describing the interior of the compact structure. This approach is based on the class I space–time and the extended gravitational decoupling by means of an extended geometric deformation (EGD). The class I condition provides a differential equation relating both metric potential ν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu $$\end{document} 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}, whilst the EGD translates the metric potentials to ν=ξ+βh(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \nu =\xi +\beta \,h(r)$$\end{document} and λ=-ln[μ+βf(r)]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \lambda =-\ln [\mu +\beta \,f(r)]$$\end{document}, where h(r) and f(r) are the deformation functions and β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} a dimensionless constant. In this case the pair {ξ,μ}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\xi ,\mu \}$$\end{document} represents the seed solution satisfying the class I condition without any deformation. Once the deformed metric potentials are inserted into the class I, the main task is to obtain h(r) or f(r). So, in this case a particular ansatz for h(r) is considered in conjunction with β=0.5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =0.5$$\end{document} to get f(r). In order to check feasibility of our model, we have performed a thoroughly physical, mathematical and graphical analysis.


Introduction
Last few years, the embedding theorem attracted the attention among the researchers working in the area of mathematical and physical sciences. Through this embedding theory, it can be found a link between the classical general relativity and higher dimensional spaces that can explain the inner symmetry clusters of the elements. The researchers [1][2][3] have linked different manifolds by embedding the 4dimensional Einstein field equations into a 5-dimensional flat space-time. This methodology has been used in string theory, supergravity, induced matter theory and membrane theory, for many viable aspects. In this connection, in the following a e-mail: sunil@unizwa.edu.om (corresponding author) b e-mail: francisco.tello@ua.cl (corresponding author) c e-mail: mahmoodkhalid@unizwa.edu.om Refs. [4][5][6], it was shown that an m-dimensional Riemannian space V m can be locally and isometrically embedded into a pseudo-Euclidean space (V n ), where n = m(m + 1)/2. In this case, n − m = m(m − 1)/2 denotes the "embedding class" of the immersed Riemannian space. This extra dimension is required for minimal embedding, although for symmetric cases, any m dimensional spherically symmetric space-time (V m ) needs at least a (m + 1) dimensional flat space-time (V m+1 ) to be immersed. It is well-known that the interior Schwarzschild space-time solution is of class I while the exterior Schwarzschild solution represents a space-time of embedding class II. On the other hand, the Friedmann-Lemaître-Robertson-Walker (FLRW) and the de-Sitter cosmological models are only of class I. In the present paper, we study the static spherically symmetric class I space-time case. It worth mentioning that, on the background of embedding class I technology for spherically symmetric spacetime, there are only two perfect fluid solutions, namely (i) Schwarzschild solution [7] and (ii) Kohler and Chao [8] solution. The class I condition, provides an extra differential equation (also known as Karmarkar condition [9]) in the static spherically symmetric case that connects both metric potentials corresponding to r −r and t −t metric components. This condition provides an extensive tool for solving Einstein's field equations and investigate new relativistic astrophysical compact stellar structures. Recently, several solutions were obtained in the context of charged and anisotropic matter distributions with well defined compact structures for class I space-time [10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26] (and references contained therein).
In the presence of anisotropy matter distribution in the stellar structures, it was shown by Dev and Gleiser [27,28] that the mass, compactness as well as red-shift are affected by this anisotropic behavior. The Ruderman [29] and Canuto et al. [30][31][32] theoretical studies, pointed out that a stellar structure may create the anisotropy inside it if the matter density is greater than the nuclear density i.e. p t = p r when ρ m > ρ n . In this direction, an inventive contribution has been done by Bowers and Liang [33] the properties of anisotropic relativistic stellar configurations, for the static spherically symmetric space-time case. Moreover, as Mak and Harko have resisted [34][35][36], anisotropies can be unified in several backgrounds such as: by the presence of type 3A superfluid [37], pion condensation [38] or numerous kinds of system transitions [39]. Although the above ingredients, have been taken into account to explain the presence of local anisotropies inside compact astrophysical objects, a recent study by Herrera showed [40] that the anisotropic state of the matter distribution in the stellar interior seems to be a natural state regarding the compact configurations in static hydrostatic equilibrium. Moreover, even if the system is initially described by an isotropic matter distribution, it naturally falls into an anisotropic phase. As it well known that the system attains a static equilibrium stage after a dynamical process, where the matter distribution threading the stellar interior is susceptible to suffer dissipative fluxes, energy density inhomogeneities and/or the appearance of shear in the fluid flow, ending with an anisotropic fluid distribution. Furthermore, departure from the isotropic pressure condition leads to get a more stable system. Last several years, a large number of researchers have been investigated effectively the role and impact of the anisotropy within compact stellar objects . Interestingly, recently a simple, powerful and versatile tool denominated gravitational decoupling by minimal geometric deformation (MGD) [78][79][80][81][82][83][84][85][86][87][88] was developed to introduce anisotropic behavior into the matter distribution, describing a wide range of solutions within the framework of Einstein gravity theory and beyond it . In this connection, some pioneering works on MGD have been done by da Rocha and his collaborators in different aspects [113][114][115][116][117][118]. Furthermore, the inverse problem i.e., given a minimally deformed and anisotropic space-time: what is its non-deformed and isotropic counterpart? has been worked in [119].
The gravitational decoupling by MGD is based on the introduction of a new material source coupled to the energymomentum tensor via a dimensionless constant parameter. Specifically it reads where θ i j is the mentioned new piece, which in principle could be a scalar, vector or tensor field, and β is the dimensionless parameter.  1 In the simplest case, when the seed solution is representing an isotropic fluid, this scheme translates it to an anisotropic domain, and as was discussed earlier an anisotropic material content recreate a more realistic scenario from the astrophysical point of view. So, in that case at general level the new set of equations has seven unknowns, namely the isotropic densityρ, the isotropic pressurep, the geometry {e ν , e λ } and the pieces {θ t t , θ r r , θ ϕ ϕ }. On the other hand, when the kernel solution already contains anisotropies the number of unknowns rises from seven to eight since there are a radialp r and tangentialp t pressures. In such a case the re-anisotropization of the fluid introduces a strong anisotropic behavior which helps to support the stability and equilibrium of the system. Nevertheless, solve this intricate problem is not an easy task. To tackle it the gravitational decoupling by means of MGD is employed. This is a purely geometry argument where the radial component of the metric tensor e λ is deformed, being f (r ) the decoupler function. As can be seen (2) introduces some modifications on the mass function m(r ) since both are connected. Once the above map is inserted into the global set of equations, it is separated into two system of equations: (i) one is completely determined by the seed space-time and (ii) the second one contains the θ -sector and the decoupler function. An interesting point to be noted here, is that each separate system satisfies its own conservation equation (null divergence ofT i j and θ i j ). This means that both sources only interact gravitationally. The main point here now, is how to close the θ -system. In this concern, several proposals have been worked in the literature in order to solve the extra equations, yielding to an admissible solution. For example the more common is the mimic constraint approach [88,89,108], the imposition of a suitable deformation function f (r ) [96,106] and more recently by using the so-called Cosenza-Herrera-Esculpi-Witten ansatz [110]. However, a more interesting situation arises when the above methodology includes deformation on both metric potentials. Developed in [120] the extended case or extended geometric deformation (EGD from now on) contemplates the same modification introduced on the energy-momentum tensor but this time the deformation (2) is extended to Despite a new degree of freedom h(r ) appears into the θsector, still is possible to solve the problem at least mathematically. So, taking into account the precedent antecedents and considering that the study of compact structures such as neutron or quark stars, compromises and active research field, we have employed the so-called class I methodology in conjunction with the EGD to obtain anisotropic compact configurations (for further details see [121]). In the seminal work [120] to close the θ -system was assumed that the extra piece θ i j is equal to a static electric field. In this case we want to obtain the concrete form of the source θ i j , hence we have imposed an adequate temporal deformation h(r ) to find out the complete θ -sector and the function f (r ). Moreover, we have employed Karmarkar condition as an auxiliary equation to find the decoupler function f (r ). From the mathematical point of view the class I constraint yields to a non-linear ordinary differential equation in f (r ) or a second order one in h(r ). To test the feasibility of this approach, we have taken the seed space-time solution to be the hybrid Adler-Finch-Skea solution. It is worth mentioning that the seed solution must satisfies the field equations under class I condition in the limit β → 0. The complete transformed space-time, has been tested by performing an exhaustive analysis of the main features that any compact structure describing an aniso-tropic matter distribution should fulfill. To recreate a more realistic scenario, we have checked all these properties by using real observational data corresponding to some compact stars. It should be noted that deformation on both metric functions, modifies the complete relativistic hydrostatic balance equation. Since the gravitational gradient is proportional to the first derivative of the g tt potential, its magnitude is changed.
In this respect, we have checked the hydrostatic balance of the resulting model, concluding that the system is in complete balance under the hydrostatic, gravitational and anisotropic gradients. This fact is very important, since the fulfillment of the modified relativistic hydrostatic equilibrium equation, tells us that the conservation law of the energy-momentum tensor is also preserved when new sources are present. Some recent works on EGD model can bee seen in Refs. [122,123]. So, the article is organized as follows. In Sect. 2 gravitational decoupling by EGD and class I methodology are revisited. In Sects. 3 and 4, are discussed the junction condition process in this context and the class I condition is presented as generator of the deformation functions, respectively. Section 5 presents the most general requirements that any compact anisotropic solution, should satisfy in order to describe an admissible solution. In Sect. 6 the model is presented and in Sect. 7 some implications on the main macro physics parameter, such as the total mass M compactness factor u and surface gravitational red-shift z s are discussed. Finally, in Sect. 8 some conclusions and remarks for the reported study are provided. The generating function for this solution has been provided at the last in the Appendix A. An extended gravitational decoupling by means of extended geometric deformation (EGD) technique, is an innovative approach to create or generalize solutions of Einstein's field equations, for anisotropic or electric charge matter distributions [88]. As it is well-known, the matter distribution contains aniso-tropy if p r = p t . There are several methods to introduce local anisotropies in self-gravitating systems.
For this purpose, we introduce an extra gravitational source which is coupled to the energy-momentum tensor associated with the seed solution. Then the modified energy-momentum tensor can be given by [88,125], whereT i j denotes the energy-momentum tensor corresponding to the seed solution for the compact star model. Then, we consider a static spherically symmetric space-time in Schwarzschild like coordinates x = {t, r, θ, φ} as It is well-known that the class I isotropic fluid solution for the above spherically symmetric metric, is compatible with only two perfect fluid solutions namely, the Schwarzschild [7] interior solution (or its special cases of the de-Sitter Universe or the Einstein universe) and the Kohler-Chao [8] solution. Nevertheless, as we argued before the main aim is to describe compact structures subject to (5) with an anisotropic content. Therefore, in order to obtain it, we shall consider that the matter distribution threading the interior embedding class I seed space-time, is described by the most general form of an imperfect fluid distribution, expressed bŷ being χ i = e −ν/2 δ i t the time-like four velocity of the fluid, satisfying χ i χ i = 1, whereρ is the energy-density,p r andp t are the radial and tangential pressures, respectively. However, the new field θ i j also introduced anisotropies into the self-gravitating system. Then the anisotropic behavior within the compact configuration becomes stronger. Now, the Einstein field equations corresponding to modified energymomentum tensor T i j are given as follows, Then the explicit form of Einstein's field (7) equations associated with the line element (5) and the matter distribution (4) can be read as, where primes denote differentiation with respect to the radial coordinate r . From now on, we shall employ relativistic geometrized units where Newton's gravitational constant G and the speed of light c are set 2 to: G = c = 1 and modified density and pressures can be written as, Besides, the modified energy-momentum tensor T i j must satisfy the conservation law leading tô where the function L(θ i i ) is given by The above Eq. (15) is a linear combination of the expressions (8)- (10). It is necessary to point out that the presence of θterm leads an extra anisotropies if θ r r = θ ϕ ϕ only. Thus, the effective anisotropy is defined as In the present situation, the extra term β(θ r r − θ ϕ ϕ ) produces a stronger anisotropic behavior into the self-gravitating system that helps to increase the stability and equilibrium mechanisms. 2 This convention will be assumed throughout the article in order to reduce the notation. However, for presenting of the numerical values of the main salient physical quantities such as density, radial and tangential pressures, these constants will be restored in the correct units form in order to express the above mentioned physical quantities. In the Table 2, the used numerical values of constant G and c as are follows: G = 6.673 × 10 −8 cm 3 /gs 2 and c = 2.997 × 10 10 cm/s. The introduction of the new piece changes the usual mass function definition by where m 0 (r ) denotes the mass function corresponding to pure general relativity (GR hereinafter) case, while m β (r ) is the extra contribution in the original GR due to gravitational decoupling. Normally, it is not an easy task to find out interior analytical solutions of the Einstein field equations for relativistic compact objects. Therefore, in order to find exact solutions, a famous choice is the extended gravitational decoupling via the EGD approach. The essential facts of this procedure depend on the resulting deformed metric potentials given by Eq. (3), which shows the deformation of the metric components g rr and g tt consisting in two unknown functions h(r ) and f (r ) namely, the radial and temporal deformation functions. However, the remaining features will be the same due to the spherically symmetry of the solution. Under this situation, the components of the θ i j sector are formed by the functions h(r ) and f (r ) and their derivatives, in distinction with what happens in the MGD case (2), where the θ i j tensor is constructed only with the radial deformation i.e., f (r ) and its derivatives.
Splitting the system of differential equations (8)-(10) and taking into account Eq. (3) the set (8)-(10) is divided into two different systems: (i) the first system satisfies Einstein field equations for an anisotropic system, in absence of β parameter i.e. β = 0, being the general solution, with a mass function given by and subject to the following conservation equation On the other hand, the second system of equations for the θ -sector are determined by turning on β. So, one has with the expressions Ψ (r ), In this case the mass function m β (r ) is expressed by It is important to remark that the m β (r ) strongly depends on the sign of β and f (r ). The consequence is that the extra piece m β (r ) will increase the total seed mass or decrease. The respective conservation law for the above field equations (25)- (27) associated with the θ -sector is [120], The above equation is a linear combination of Eqs. (25) and (27) through the coupling constant β.

Embedding class one space-time
It is well known that a space-time describes an embedding class I space-time if the 4-dimensional space-time is embedded into 5-dimensional pseudo-Euclidean space. This can happen if there exists a second fundamental form symmetric tensor K i j = K ji , that satisfies the Gauss-Codazzi equations, here takes the value 1 or −1 corresponding to signature required for the manifold, R hi jk is the Riemann tensor and ∇ k the affine connection allied to the metric tensor g i j such that ∇ i g jk = 0.
The only non-zero components of the symmetric tensor K hi are: K rr , K tt , K θθ = sin 2 θ K φφ and K tr = K rt . By plugging these components into (31), we get where R r θr θ = 0 [64] and the non-zero components of Riemann tensor R hi jk associated with the condition (33) in the context of static spherically symmetric line element (5) are given as, Now by inserting the above components (34)-(37) into class I condition (33), we derived the following ordinary differential equation The above equation can be solved for both unknowns λ(r ) and ν(r ) and can also be expressed in terms of each other as, whereÃ,B andC denote arbitrary integration constants.

Junction conditions for anisotropic stellar models
The study of the junction conditions for anisotropic stellar models at the surface r = R, is an essential measure that can be described by joining the interior space-time M − (r < R) and the exterior one M + (r > R) at the boundary of the stellar compact object. In our current scenario, the interior stellar space-time can be described by the following line element in the context of the extended geometric deformation approach as where is internal mass of the anisotropic stellar structure. The interior space-time (40) must be joined in a smoothly way with the empty exterior space-time i.e. ani-sotropic pressures p + i and density ρ + must be zero. In general, the outer spacetime is described by the vacuum space-time i.e., exterior Schwarzschild solution. Nevertheless, the presence of the θsector into the matter distribution could in principle change both the geometry and the matter content for the outer spacetime. In this scenario, the stellar compact object will be not immersed in vacuum space-time anymore. In this situation, the compact stellar model can remain embedded into a vacuum space-time whether the contributions coming from the θ -sector are considered to be confined within the stellar interior only [88]. Therefore, the exterior Schwarzschild spacetime can be given by being M Sch the Schwarzschild mass which coincides with the total mass M of the object at the boundary Σ. For smooth joining of inner geometry with outer geometry at the boundary of the stellar interior at r = R, we apply the well known Israel-Darmois junction conditions [65,66]. These conditions are known as the continuity of the first and second fundamental forms across the boundary Σ. The continuity of the first fundamental form at the boundary of the stellar model Σ provides yielding to where the total mass M of the compact object can be determined as Next, the continuity of the second fundamental says being r j a unit vector. From Eqs. (4) and (12) one obtains thuŝ wherep r (R) denotes the radial pressure at the boundary (r = R) corresponding to interior matter distribution when β = 0. The Eq. (49) describes the general expression for the second fundamental form associated with the Einstein field equations given by Eq. (7). Now, using Eq. (26) we arrive tô here ν ≡ ∂ r ν − . In order to find out (θ r r ) + (R) in (50), we employ Eqs. (26), (44) and (45), obtaininĝ where, f * (R) and h * (R) denote the geometric deformation functions for the exterior Schwarzschild solution (42). It is essential to note that, if the exterior space-time is described by a vacuum space-time, f * (R) and h * (R) must be null. Therefore, Eq. (51) becomeŝ The condition (52) determines the size of the compact stellar object. Besides, it is worth mentioning that conditions (44), (45) and (52) describe the necessary and sufficient conditions to determine all the constant parameters.

The embedding class one generator for anisotropic stellar object
In this section we develop a generating equation in f (r ) and h(r ) for the embedding class I space-time via gravitational decoupling by means of EGD methodology. For this purpose, we substitute the extended deformed gravitational potentials ν(r ) and λ(r ) from Eq. (3) into embedding class I condition (38), obtaining a differential equation, that involves 4-unknowns, namely f (r ), g(r ), ν(r ) and λ(r ), The above Eq. (53) generalizes, all known spherically symmetric embedding class I solutions in the context of extended gravitational decoupling via a extended geometric deformation (EGD) approach. It is observed that Eq. (53) can be written either as a first order non-linear ordinary differential equation (ODE) for f (r ) or as a second-order non-linear ODE for h(r ). Besides, it is possible to generalize some previously known embedding class I solutions (an almost complete list of such solutions can be found in [21]) through this scheme. However, there are exceptions with some solutions which have a fairly complex geometry. So the differential equation (53) cannot be solved analytically.

Physical conditions of the anisotropic stellar object
The viability of any stellar compact object, describing a realistic stellar structure such as neutron stars, must satisfy some physical and mathematical requirements [67] -The gravitational functions, namely e ν(r ) and e λ(r ) must be positive, finite, free from singularities, and monotone increasing functions with increasing r . Also, e ν(r ) and e λ(r ) must satisfy e ν(0) > 0 and e λ(0) = 1. -In considering the main thermodynamic variables namely {ρ, p r , p t }, it is desirable that they are positive everywhere within the stellar interior for the stability reasons.
Since we are within a classical 3 context, then the condition of positivity must be strictly fulfilled for ρ, not for p r and p t . However, the presence of negative pressure inside the system can be also allowed [69]. Besides, the presence of this negative pressure leads to the appearance of unstable states. Clearly, these states would impair the hydrostatic balance of the system. -The behavior of the thermodynamic observables, should be monotonically decreasing one i.e., they keep their maximum value at the center of the stellar model and their minimum at the boundary. -The radial and tangential pressures, must be matched at the center. On the other hand, on the surface, the radial pressure must vanish (this condition is not necessary to be fulfilled for the tangential pressure). -Within the compact stellar structure, the subliminal sound velocities of the fluid in the principal directions i.e., along the radial and tangential directions, must be less than the speed of light c = 1, in order to satisfy causality condition. Mathematically it can be denoted as: 0 ≤ v 2 r = dp r dρ < 1 and 0 ≤ v 2 t = dp t dρ < 1. -The energy-momentum tensor has to fulfill simultaneously the following conditions: DEC: ρ − p r ≥ 0 and ρ − p t ≥ 0, and SEC: ρ + p r + 2 p t ≥ 0.

Extended gravitational decoupling solution via embedding class I generator f(r)
Determining the solution of an ordinary first-order differential equation is usually easier than finding a solution of a second-order differential equation. Therefore, we will solve the Eq. (53) for f (r ) by specifying a well-motivated ansatz for h(r ). The Eq. (53) can be written in the following form The above Eq. (54) is called a generating equation in f (r ) (known as an embedding class I generator of the system induced by EGD approach) describing all possible solutions for the present self-gravitating system. In principle, it is not a trivial task to integrate the above first order non-linear differential equation (54) in f (r ) for all kind of ξ(r ), μ(r ) in conjunction with h(r ), otherwise, a numerical approach must be adopted. Nevertheless, it might be possible to solve exactly this equation for a particular choice of ξ(r ), μ(r ) and h(r ). On the other hand, we must choose the ξ and μ based on the facts proposed by [68,70,71] i.e., the gravitational potential ξ(r ) must satisfy ξ(0) = finite and constant, ξ (0) = 0 and ξ (0) > 0 for any realistic model. Moreover, the radial pressure and the energy density should be non-negative and continuous inside the stellar model, then it follows that r > 2 m(r ) [72,73]. So, p r ≥ 0 with r > 2 m(r ) leads to ξ (0) = 0. From all constraints of ξ(r ) in the said discussion, we observe that generic function ξ(r ) must be increasing monotonically function of r , and should attain its minimum at the centre of the stellar model. On the other hand, we should ensure that the gravitational potential μ for the radial component must be of the form μ = 1 + O(r 2 ) near at r = 0. Therefore, by keeping all the above points in our mind we have chosen the following simple embedding class I space-time as seed solution (β = 0) as, where μ(r ) = (1 + Ar 2 ) −1 , ξ(r ) = 2 ln B + Cr 2 , and the corresponding thermodynamic observables like pressures and density are given bŷ and the corresponding mass function m 0 (r ) is given by, where the constants A and C have units of [length] −2 and B is dimensionless. Now we need to solve the embedding class I generating Eq. (54) for determining the components of the θsector. For this purpose, we specify a particular form of h(r ), specifically h(r ) = D r 2 . This choice is physically motivated and simple, since it is zero at the centre and increasing function with r . By plugging h(r ), ν(r ) and λ(r ) into Eq. (54) we get, where, The above Eq. (60) is a first order non-linear differential equation in f (r ), which can be integrated in closed form only if β = 1/2, otherwise we will get a complex solution. Therefore, after integrating Eq.
where F is an integration constant with units of [length] 2 . We would like to mention that the signature of the deformation function f (r ) and coupling constant β plays an important role to decide the behavior of the metric function e λ(r ) given by Eq. (3). Since the metric function e λ(r ) must be increasing function within the star, in order to preserve the increasing behavior of the mass function. Therefore, f (r ) should satisfy some constraints 1. f (r ) > 0 (increasing function) and β > 0: In this situation, the growth of μ(r ) must be faster than the decoupler function f (r ), in order to keep e λ(r ) and the mass function m(r ) positive and increasing for all r ∈ [0, R]. 2. f (r ) < 0 (decreasing function) and β > 0: In this scenario, the metric function e λ(r ) will increase throughout the compact object. 3. f (r ) > 0 (increasing function) and β < 0: In order to preserve the increasing behaviour of ν(r ), the growth of ξ(r ) must be faster than the growth of h(r ), however the function e λ(r ) and mass function m(r ) will increase automatically. 4. f (r ) < 0 (decreasing) and β < 0: In this case, the growth of ξ(r ) and μ(r ) must be faster than the growth of h(r ) and f (r ), respectively.
Based on the above points, first we plot the deformation function f (r ) in Fig. 1 to see its behavior. From this Fig. 1, we observe that the deformation function f (r ) is negative and decreasing. Since β = 1/2 is positive and temporal deformation function h(r ) is already positive and increasing.
Therefore, according to the above discussion it is suitable to describe a realistic model. Now, the embedding class I solution in the framework of extended gravitational decoupling by EGD approach can be given by the following line element, From the above Eq. (62), it is obvious that gravitational potential ν(r ) is positive, increasing and free from singularity, although we need to see the variation of potential λ(r ) of the Eq. (62). For this purpose we plot the Fig. 2 for e λ(r ) and observe that it is 1 at the centre and increasing throughout Fig. 2 The trend of metric function e λ(r ) versus radial coordinate r/R. We have described same set of values for different curves as motioned in Fig. 1 within the stellar object. Now using Eqs. (18), (29) together with Eq. (59), the mass function m(r ) for EGD can be given as, As can be seen, we have generalised the seed mass m 0 (r ) through the gravitational decoupling because the constant β is positive and the deformation function f (r ) is negative throughout the stellar interior, and then this extra piece m β (r ) will increase the total mass M of the object. Next, we determine the expressions for components of θ -sector namely θ t t , θ r r and θ where, Then, the effective physical quantities like pressures, ( p r ) and ( p t ), density (ρ), and pressure anisotropy (Δ) for the energy-momentum tensor T i j can be given as, Now, to determine the constants, we use the junction conditions (44) and (45) with coupling constant β = 1/2, yielding to Then from Eqs. (68) and (69) together with second fundamental form (52), we have determined the values of the constants B, C and F (which are necessary to describe the complete structure of the stellar model) as, The variation of pressures ( p r and p t ), density (ρ) and anisotropy factor (Δ) are shown in Figs. 3 and 4.

Causality and hydrostatic equilibrium
In order to fulfill the physical requirements for realistic models, it is necessary to examine the causality and hydrostatic equilibrium of the present self-gravitating system. First, we discuss the causality condition of the model which says that the velocity of sound must be less than the velocity of light everywhere within the object. Mathematically we can write as 0 < v 2 r = dp r dρ < 1 and 0 < v 2 t = dp t dρ < 1, where v 2 r and v 2 t represent the square of radial and tangential velocity, respectively. As we can see from Fig. 5, both v 2 r and v 2 t are less than the speed of light (c = 1) and decreasing throughout the structure. Moreover, the radial velocity dominates the tangential one at each point inside the object. Therefore, it is necessary to analyze whether the given model is stable or not. To examine this situation, we use Abreu's criteria [74] which was initially created by Herrera, namely Herrera cracking concept [43]. The Abreu's criterion for the stability analysis of compact objects states that the model is stable if t < 0 provides the unstable regions. From the Fig. 5 it is clear that our anisotropic model is stable.
On the other hand, the analysis related to the hydrostatic equilibrium under different forces can be analyzed by following modified Tolman-Oppeneheimer-Volkoff (TOV) equa- Fig. 3 The behavior of the radial pressure ( p r ) tangential pressure ( p t ), and energy density (ρ) versus radial coordinate r/R. The different curves describe the following compact objects namely the red curve (solid) for 4U 1538-52, black curve (long-dashed) for SAX J1808.4-3658, green color (small-dashed) for SMC X-1, vermilion color (dotdashed) for LMC X-4 Fig. 4 The trend of the anisotropic factor (Δ) versus radial coordinate r/R. The description and numerical values for the free parameters are same as used Fig. 3 tion, The equilibrium equation can be spread into three different forces namely: the hydrostatic F h , the gravitational F g and the anisotropic F a forces in order to achieve this hydrostatic balance. Moreover, the explicit form of these forces can be written as It is important to note that the well-known TOV equation [128,129] in the context of standard relativistic anistropic self-gravitating system can be recovered, if β = 0. From Fig. 6, it is clear that the anisotropic system is in equilibrium under the said forces. As mentioned before, the anisotropic force is repulsive in nature, then the gravitational force is counterbalanced by the joint action of the hydrostatic and anisotropic forces. This avoids the system to collapse below its Schwarzschild radius onto a point singularity. Moreover, it is noted that the anisotropic force dominates the hydrostatic one after some points inside the stellar model which is showing that local anisotropies have a majority role in the balance of the stellar configuration.

The compactness factor and surface red-shift
In this section, we discuss the influences induced by extended gravitational decoupling by EGD approach on the compact-  ness factor u = m(r ) r within the stellar interior. In order to provide a detailed explanation, we proceed through the mass function in the EGD case as defined by Eq. (41). In the present situation, the constant β is positive while deformation function f (r ) is negative and decreasing outward. Therefore, the mass function in EGD scenario will dominate the mass function in GR case, which can be clearly observed from Fig. 7. Now the surface red-shift function of the compact star model can be given by the following formula It is well-known that the impact of anisotropies on the surface red-shift z s has been widely discussed by several authors. Bowers and Liang [33] proposed a hypothetical model having a constant density ρ = ρ 0 (incompressible fluid) with a particular choice of the anisotropy factor Δ and they determined that the surface red-shift can reach its maximum value up to z s = 4.77 when Δ = 0. In addition to above, the investigation, Ivanov [77] pointed out that the maximum value of the surface red-shift can be z s = 3.842 if the aniso-tropic star model (without cosmological constant) satisfies the strong energy condition (SEC) with mass-radius ratio 0.957, while for models obeying the dominant energy condition (DEC) having a mass-radius ratio 0.974 reach up to z s = 5.211 [36]. In this case, the surface red-shift function is depicted in Fig. 8. As it is observed, z s is increasing with r and attains its maximum value at the boundary. It is interesting to note that the growth of the surface red-shift function, is faster corresponding to the compact star having less value of the compactness factor. On the other hand, from Table 1 one can see that the surface red-shift value is increasing with an increasing compactness factor. This implies that the surface red-shift cannot be arbitrarily large due to the boundedness of the compactness factor, which means that any massive realistic compact object must satisfy the Buchdhal limit i.e. M R < 4 9 [75]. So, based on the above discussion the surface red-shift for extended gravitational decoupling model will be more than the pure GR model due to the compactness factor in the GR+EGD scenario is greater than the mass-radius ratio in the pure GR case (Fig. 7).

Discussion and conclusion
In the present article, we have developed a generating scheme to obtain solutions of Einstein field equations for modelling compact stars. To achieve it, we have employed the embedding class I generator, induced by an extended gravitational decoupling via Extended geometric deformation (EGD) approach. In general this systematic approach can be summarized as follows: 1. First, we define the modified energy-momentum tensor T i j , which is a combination of two sources namelyT i j and θ i j . HereT i j is the energy-momentum tensor for the standard anisotropic matter distribution while θ i j denotes an extra energy-momentum tensor that is coupled toT i j via a dimensionless coupling parameter β. 2. Next, we write Einstein field equations for a spherically symmetric line element (with metric functions e λ(r ) for the r − r component and e ν(r ) for the t − t component) subject to T i j . 3. Then, we split the resulting system of equations into two subsystems, by using the following transformations ν(r ) → ξ(r )+β h(r ) and e −λ(r ) → μ(r )+β f (r ), where f (r ) and h(r ) denote the radial and temporal deformation functions, respectively. In this way, the first subsystem is called the Einstein system (which corresponds to set β = 0) corresponding to matter distributionT i j and sec- ond subsystem corresponding to θ -sector is known as a quasi-Einstein system. 4. Now, in order to solve these systems, the embedding class I condition (38) is adopted. After plugging the transformations given by Eq. (3) into (38), we arrive to a non-linear differential equation (53) in terms of the unknown functions ν(r ), λ(r ), f (r ) and h(r ). The Eq. (53) is a generating equation that can generalize any known embedding class I seed solution corresponding to compact star models, by taking a compatible ansatz for anyone generator f (r ) or h(r ). Nonetheless, in general, to obtain the solution of a second-order non-linear differential equation is more difficult than the first order. So, in this case, we solved a first-order non-linear in f (r ), choosing a known embedding class I seed solution and a particular ansatz for h(r ). It is worth mentioning that the equation can be solved in closed form only for β = 1/2. For this specific model, we have obtained a negative and decreasing deformation function f (r ) (see Fig. 1). Nonetheless, its behavior does not damage the monotonic increasing behavior of the e λ(r ) metric potential. In effect, as can be appreciated from Fig. 2, the radial metric potential behaves as desired. Then, we can conclude that the resulting function f (r ) is physically viable. So, the output model constitutes a generalized embedding class I space-time solution, driven by an anisotropic matter distribution, represented by the energy-momentum tensor T i j .
In considering the thermodynamic observables, namely effective radial pressure ( p r ), effective tangential pressure ( p t ) and effective energy density (ρ), they are shown by Figs. 3 and 4, respectively. It is clear from both figures that the radial pressure, tangential pressure and energy density behaviors are positively, well-defined and decreasing inside the stellar configuration, as well as both pressures are equal at the centre. Moreover, all the above physical quantities are attaining their maximum values at the center of the structure and reaching their minimum values at the boundary r = R of the stellar configuration. The variation of the anisotropy factor Δ is presented in Fig. 4. As it is observed, the trend of these curves show that the anisotropy is directed outward at all points within the object, i.e. Δ(r ) > 0 for 0 < r ≤ R.
The said features of the model indicate the physical validity of the solutions.
Another important analysis, concerns the so-called causality condition. For this purpose, we have plotted both v 2 r and v 2 t in Fig. 5. From this figure, it is clear that v 2 r and v 2 t are less than the speed og light (c = 1) and decreasing throughout the stellar model. Moreover, as can be seen, v 2 r dominates v 2 t at each point inside the model. Therefore, according to Herrera's cracking concept as well as Abreu's criteria [74], the obtained solution is stable against radial perturbations introduced by local anisotropies. Besides, the stability of the model is quite involved with the hydrostatic balance of the stellar system. In order to verify it, we have used the modified Tolman-Oppenheimer-Volkoff (TOV) equation (75). In this case the TOV equation is composed by three forces namely F g , F h and F a (given by Eqs. (74)-(76)) As Fig. 6 illustrates, the F g is balanced by combining F h and F a , what is more the sum of all these components is zero i.e. F g + F h + F a = 0, which shows that the equilibrium of the system is achieved.
From the observational point of view, two very important quantities are the mass-radius ratio and the surface red-shift. As it is well-known, these two quantities are quite involved and help to determine valuable information of stellar bodies. In this regard, we have explored the incidence of the extended gravitational decoupling by EGD on these important quantifiers. From Fig. 7, we notice that the compactness for the GR+EGD model is more than in the pure GR scenario, what is more each star model is satisfying the Buchdhal limit. Consequently, the surface red-shift will be greater in the context of GR+EGD than in the case of pure GR. The Fig. 8 shows that the growth of the surface red-shift function, is faster for those stars having less value of the compactness factor.
To confirm and contrast the previous discussion, we have obtained the numerical data of the physical parameters, reported in Table 1. To obtain this information, the following procedure has been used: Since the system has only three free constant parameters namely R, M and D, therefore in order to determine the values of all physical quantities, the constant D and mass M (corresponding to real compact objects) are fixed. After that, the radius R has been fitted in such way that the model must satisfy all the physical and mathematical requirements of the solution. It should be noted that, Table 2 The  the obtained radii matches closely with the values reported in M-R curve as proposed by Demorest et al. [133]. By using this information and with the help of Eq. (77), the surface red-shift and the compactness factor u were computed. For the present model, we have obtained the following results: z s = 0.17121 for 4U 1538-52, z s = 0.17590 for SAX J1808.4-3658, z s = 0.19523 for SMC X-1, and z s = 0.22757 for LMC X-4 (see Table 1). It is remarkable to note that the above values are within the range as proposed in Refs. [76,77]. Also, we would like to mention that the surface red-shift for the GR+EGD model will be more than the pure GR model because the compactness factor in the GR+EGD case is always larger than the compactness in pure GR scenario. On the other hand, the values of central pressure ( p c ), central density (ρ c ) and surface density (ρ s ) in CGS unit for each model is presented in the Table 2. From the Table 2, we see that the density of the model is larger than the nuclear density (ρ N = 2.3 × 10 14 g/cm 3 ).
Finally, we would like to mention that the present generating scheme via EGD approach is very useful to generalize embedding class I solutions for self gravitating compact objects, within the framework of GR. Moreover, it would be interesting to include new ingredients such as electric charge, or translate the methodology into the context of modify gravity theories such as f (R), f (R, T ) etc.