Class I approach as MGD generator

In this work we build a relativistic anisotropic admissible compact structures. To do so we combine the class I approach with gravitational decoupling in order to generate the deformation function f(r). As an example we have re-anisotropized two anisotropic matter distributions previously obtained by the class I procedure. To produce all the graphical study supporting this analysis, we have considered the data corresponding to the compact object 4U 1538-52, SMC X-1 and LMC X-4 for model 1 and Cen X-3 for model 2. In considering the last one, we have taken the constant parameter α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} to be {-0.3;0.1;0.3}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{-0.3;0.1;0.3\}$$\end{document}. It is found that the resulting models satisfy all the general requirement in order to represent or describe realistic compact structures such as neutron or quark stars.


Introduction
The embedding class I condition is a versatile and simple tool, which can be considered an auxiliary condition to solve Einstein's field equations. Staged by Karmarkar [1], this condition generally applies to any n-dimensional (pseudo)-Riemannian space-time. In short, this technique says that a ndimensional (pseudo)-Riemannian manifold can be embedded into a pseudo-Euclidean n + p-dimensional space, where p denotes the class of the embedded manifold. Since general relativity works on a 4-dimensional manifold a spherically symmetric and static space-time described by a Riemaniann variety with Lorentzian metric can be embedded within a 5diemensional pseudo-Euclidean space, then the submerged manifold is of class I. For the aforementioned case the derivaa e-mail: francisco.tello@ua.cl (corresponding author) b e-mail: sunil@unizwa.edu.om c e-mail: ygomez@ucn.cl tion of the class I condition made by Karmarkar reads The above relationship between the Riemann tensor components is valid only if R θφθφ = 0 [2]. As we will see in the following, the components of the Riemann tensor engaged in expression (1) lead to a beautiful and simple relationship between the metric potentials ν and λ that determine the geometry of the space-time under study. For the resolution of the field equations in the context of general relativity, this relationship is very helpful, given that in the context of compact solutions such as the modeling of neutron stars, whose material content is described by an imperfect fluid, the number of unknowns exceeds the number of equations. Specifically one has five unknowns: the geometry {ν, λ} and the matter content {ρ, p r , p t }. Then, with the imposition of a suitable ansatz for ν or λ the class I condition leads to the determination of the modeling geometry, whereby the thermodynamic variables that characterize the model are obtained by inserting the resulting geometry into the field equations. Therefore the problem is mathematically solved. Due to the great difficulty involved in solving the general relativity field equations, during the last years the conspiracy of class I has been widely used in the study on modeling compact structures with anisotropic matter distribution . As we mentioned earlier this technology solves the problem in the mathematical sense, however the challenge is even greater since such models must be physically admissible to represent a realistic situation. In a simpler context a compact configuration can be made of a perfect fluid matter distribution p = p r = p t , which greatly reduces the resolution of the field equations. Nevertheless, the work by Lake and Delgaty [38] showed that not all models built on the assumption of a perfect fluid distribution are admissible from the physical point of view. The inclusion of local anisotropies ≡ p t − p t in the stellar interior has a long history and has been thor-oughly studied . Due to extreme internal density and strong gravity within self-gravitating compact object indicates that pressure may not be equal, i.e., there exist two different kinds of interior pressures, namely the radial and tangential pressure [41,67]. Such kind of compact object may be developed to study the phase transitions and distributions containing the combination of two fluids [68,69]. Lemaître [70] first pointed out this outcome in the stellar structure and evolution of compact objects. Later on, Bowers and Liang [39] revived the interesting study of anisotropic relativistic matter distributions in general relativity. Using the generalization of the equation of hydrostatic equilibrium, they created a static spherically symmetric stellar structure and investigated the modifications in the surface redshift and gravitational mass. Through the theoretical studies, Ruderman in [71] jagged that nuclear matter tends to become anisotropic in nature at very high densities of order of 10 15 g/cm 3 . For a massive compact star model, they had considered that both pressures (radial and tangential) might be differ to each other. In this situation, many different opinions have been presented for the presence of anisotropy in the stellar compact objects such as by the existence of different kinds of phase transitions [72], the occurrence of solid core, combination of two fluids, presence of type 3A superfluid [73] or by other different physical aspects. Nowadays, several works have been done for anisotropic stellar models in different context [32,[74][75][76][77]79].
In this direction, recently a novel method was developed to generate anisotropic solutions to Einstein's field equations [80][81][82][83][84][85][86][87][88][89][90]. The so-called gravitational decoupling by minimal geometric deformation (MGD) was designed to extend solutions driven by an isotropic matter distribution to anisotropic domains. This scheme not only modifies the material content of the object, also deforms its geometry in such a way that the symmetry of the solution is preserved. These last two years the amount of works available in the literature using this method has grown considerably. Applications range from stellar interiors, black holes and modified gravity theories, to name a few [78,. What is more, the MGD inverse problem ı.e, given an anisotropic solution it is possible to know the isotropic counterpart, was developed in [117] and the extended case, that is, deformation on both metric potentials was worked in [118].
Recently, class I condition has been used within the framework of MGD [119]. In that work, a space-time has been generated using the embedding technique and then deformed by applying MGD. In this sense, the proposal of this article is not only to extend solutions to scenarios dominated by anisotropic matter. The fundamental idea is to use Karmarkar's condition as a generator of MGD, which in turn allows to determine the components of the new material sector responsible for the anisotropic behavior inside the compact structure. As examples we have considered the uncharged Adler-Finch-Skea solution [120,121] to be reanizotropized using the described methodology and a hybrid obtained by using the temporal metric potential corresponding to Kuchowicz space-time [122] which was worked in [13]. It is worth mentioning that the Adler-Finch-Skea solution was already obtained in the framework of class I scheme, as well as the hybrid formed by Kuchowicz metric potential. Furthermore, in both cases the matter distribution corresponds to an anisotropic charged one (see [28] and [29] for further details). However, in this case we have taken only the geometry of such solutions. To produce the profiles of the energy-density, radial and tangential pressures, anisotropy factor, velocities of the pressure waves, stability and balance mechanism we have used the observational data corresponding to the compact objects [123] 4U 1538-52, SMC X-1 and LMC X-4 for model 1 (Adler-Finch-Skea) and Cen X-3 for model 2 (Kuchowicz). The article is organized as follows: Sec. 2 presents in brief the class I approach and the gravitational decoupling by means of minimal geometric deformation. In Sect. 3 the anistoropic stellar interior solution is provided. Next in Sect. 4 are analyzed the main results and finally in Sect. 5 some remarks are reported.

Embedding class I
A space-time is said to be of class I ı.e, admits to be embedded into a 5-dimensional pseudo-Euclidean space, if there exist a second fundamental form symmetric tensor K σ γ = K γ σ satisfying the Gauss-Codazzi equations where = ± (according to the normal to the manifold being time-like "-" or space-like "+"), R σ γβω is the Riemann tensor and ∇ ω the affine connection associated to the metric tensor g γβ , ∇ γ g βω = 0. Regarding the spherically symmetric and static space-time given in Schwarzschild like coordinates x = {t, r, θ, φ} by the only non trivial components of the symmetric tensor K σ γ are: K tt , K rr , K θθ = sin 2 θ K φφ and K tr = K rt . By inserting these components into (2) one arrives to The Riemann tensor components compromised in the previous condition and associated to the line element (4) are R θφθφ = −e −λ r 2 sin 2 θ e λ − 1 , Next, replacing (6)- (11) in Eq. (5) we obtain the following differential equation with e λ = 1, from where being A an integration constant. Equation (12) also can be solved to express the metric potential ν in terms of λ as follows where B and C are integration constants.

Gravitational decoupling by MGD
Minimal geometric deformation (MGD from now on) approach is a novel tool useful to generate anisotropic solutions of the Einstein field equations starting from an isotropic (anisotropic) one [90]. In general, there are many ways to introduce local anisotropies. In this regard we will focus on the case where the shear term of the energy-momentum tensor is taken to be null. Hence, the anisotropic behaviour appears when p t − p r = 0. In order to produce the aforementioned anisotropies, it is necessary to introduce an extra gravitational source which, in principle, can be e.g. a scalar, vectorial or tensorial field. This extra source is coupled to the energy-momentum tensor associated to the seed solution. The effective energy-momentum tensor can be defined as follows [90,124] T μν ≡T μν + αθ μν , whereT μν corresponds to a perfect fluid given bỹ being χ μ = e −ν/2 δ μ t the time-like four velocity of the fluid satisfying χ μ χ μ = 1,ρ the isotropic energy-density and p the isotropic pressure. The new field θ μν encodes the anisotropies introduced into the system. Now, the Einstein field equations associated with the geometry (4) and the matter distribution (15) explicitly reads where the primes denote differentiation with respect to the radial coordinate r . Besides, hereinafter we shall employ geometrized relativistic units where G = c = 1. The full diffeomorphims symmetry entails the conservation of the energy-momentum tensor which reads where the function H (θ i i ) encodes the corresponding anisotropies and it is defined as: The above expression (22) is a linear combination of Eqs.
It is essential to point out that the inclusion of θ -term introduces anisotropies if θ r r = θ ϕ ϕ only. Thus the effective anisotropy is defined in the usual manner, namely: Naturally, we recover a perfect fluid when α is taken to be zero. On the other hand if the seed solution already contains anisotropies inside the matter content the energy-momentum tensor is described by an imperfect fluid distribution as follows withp r andp t being the pressure waves in the principal directions ı.e, the radial and tangential ones respectively. The four-velocity of the above fluid distribution is characterized by the time-like vector χ ν . Moreover u ν is a unit space-like vector in the radial direction (orthogonal to χ ν ). In that case (27) becomes to So, in this case the extra term α θ r r − θ ϕ ϕ introduces a stronger anisotropic behaviour into the matter distribution. This serves to improve the stability and equilibrium mechanism. In general, it is not trivial to obtain analytic solutions of the Einstein field equations in the context of interior solutions, ı.e., relativistic stars. To find a tractable exact solution, albeit recent, a popular alternatives is the gravitational decoupling via the MGD approach. The crucial point of this technique relies on the following map: in which we deform minimally the g tt and g rr components of the metric. The later maps deform the metric components by the inclusion of certain unknown functions h(r ) and f (r ). At this level, it is noticeable remark that the corresponding deformations are purely radial. The later feature remains the spherical symmetry of the solution. The so-called MGD The first case maintain the deformation in the radial component only, which means that any temporal deformation is excluded. In light of this, the anisotropic tensor θ μν is produced by the radial deformation (31). The system of differential equations can be split under the replacement (31). Thus, the field equations are naturally decoupled in two set: i) the first set satisfies Einstein field equations, and correspond to the isotropic (anisotropic) case, namely α = 0, and it is given by along with the following conservation equatioñ and ii) the second set of equations corresponds to the θ -sector which is obtained when we turn on α. Thus, the equations for the later sector are: The corresponding conservation equation associated to the θ -sector is computed to be It is important to remark that the above equation is precisely the essential point to use the MGD approach, given that it guarantees that the energy interchange is pure gravitational only. At this stage it should be noted that if the seed solution is anisotropic the pressure in the left hand side of Eqs. (33)-(34) must be replaced byp r andp t , respectively. Moreover, the conservation equation is modified tõ

Anisotropic stellar interiors
In this section we provide two examples of anisotropic stellar interiors. Specifically, we re-anisotropize the Adler-Finch-Skea model previously obtained in [28] and second one was worked in Maurya et al. [29] (although in these works the solutions include electric charge we will take only the geometrical description).
To do so the embedding class I technology is combining with the MGD machinery to generate the deformation function f (r ). Putting together Eqs. (12) and (31) one arrives to In principle this first order non-linear differential equation (41) in f (r ) looks to complicated, however it is possible to integrate this equation for specific choice of ν(r ) and μ(r ). It should be remembered that the seed solution corresponds to a solution of the Einstein field equations with anisotropic matter distribution which was obtained using the Karmarkar condition (13). A (almost complete) list of such solutions can be found in [6]. Although some solutions have a fairly complex geometry, the differential equation (41) can be solved in most cases, except for some exceptions where the seed solution is described by very complex metric functions, such as hypergeometric functions. Moreover in the works [125][126][127][128], it was proved that in the modelling of realistic compact objects, the gravitational potential ν(0) = finite and constant, ν (0) = 0 and ν (0) > 0. On the other hand, the radial pressure and the energy density must be positive and continuous within the compact objects which yields r > 2 m(r ) [129,130]. Then from p r ≥ 0 with r > 2 m(r ), it follows that ν (0) = 0. This implies that generic function ν(r ) attains its regular minimum at centre and increasing monotonically function of r . Also we ensure that the another obtained gravitational potential e λ(r ) should be the form e λ(r ) = 1 + O(r 2 ) near at r = 0. By keeping all above mathematical and physical points in our mind we have chosen two different space-time to modeling compact structures, which are given below.

Model 1
To find the decoupler function f (r ) we impose the following class I seed space-time From Eqs. (42)-(43) it can be clearly observed that this spacetime is fulfilling the physical and mathematical requirements mentioned above (see Fig. 1). Moreover, the matter content inside the stellar interior described by the above geometry is given bỹ where the constants B and C have units of [length] −2 and A is dimensionless. So, the decoupler function obtained from (41) using Eqs. (43)- (42) is being F an integration constant with units [length] 2 . As it is observed at r = 0 the deformation function f (r ) is zero. furthermore, it is dimensionless. So, the deformed space-time reads Now, by inserting Eq. (47) into the set of Eqs. (36)- (38) provide the following expressions for the θ -sector components So, by virtue of Eqs. (45), (46), (50) and (51) the anisotropic factor (29) is given by As can be seen at the center r = 0 of the compact star (52) is vanishing ı.e, (0) = 0 as it is required.

Model 2
In order to find the decoupler function f (r ) for the model 2, we impose the following another class I seed space-time Again it is observed that from From Eqs. (53)-(54) the spacetime is satisfying the physical and mathematical requirements mentioned above (see Fig. 5). The matter content inside the stellar model 2 can be described using the above geometry as follows, ρ = C e A+Br 2 (3 + 2Br 2 + C r 2 e A+Br 2 8π 1 + C r 2 e A+Br 2 2 , being F an integration constant with units [length] 2 . As it is observed at r = 0 the deformation function f (r ) is zero, and has no dimension. Then the deformed space-time for the Model 2 can be read as, Now, by plugging Eq. (55) into the set of Eqs. (36)- (38) which yield the following expressions for the θ -sector components Then, from Eqs. (56), (61) and (62) the anisotropic factor (29) is given by It can be clearly noted that the anisotropy, given in Eq. (63), vanishes at the center of the compact star i.e. (0) = 0 at r = 0 which is required for physical acceptability.

Physical analysis
The feasibility of any model describing the interior of a compact object representing realistic structures such as neutron stars must satisfy some requirements to be physically and mathematically admissible [40,56] • The metric potentials, namely e ν and e λ must be free from singularities, finite and monotone increasing functions with increasing radius. Besides, e ν(0) > 0 and e λ(0) = 1. • The main thermodynamic variables, namely {ρ, p r , p t } should be strictly positive functions at every point inside the configuration. • The behavior of the thermodynamic observables responds to a monotonically decreasing one ı.e, their maximum values are attained at the center of the object reaching their minimum at the surface. • Both the radial and tangential pressure coincide at the center. Moreover, at the boundary the radial pressure must vanish and the tangential one not necessarily is. Fig. 1 The trend of metric potential against the radial coordinate r/R for C = 0.2 km −2 , α = 0.008km and different values mentioned in Table 1 for compact objects 4U 1538-52, SMC X-1 and LMC X-4 • Inside the compact static structure the velocity of the pressure waves in the principal direction of the sphere must be less than the speed of light c = 1 in order to preserve causality condition: 0 ≤ v 2 r = dp r dρ < 1 and 0 ≤ v 2 t = dp t dρ < 1. • The energy-momentum tensor has to satisfy simultaneously the following conditions: ρ − p r − 2 p t ≥ 0 and ρ + p r + 2 p t ≥ 0.
In addition to meet the above requirements, at the boundary : r = R the model describing the stellar interior M − should be joined in a smoothly way with the corresponding exterior space-time M + . As we are dealing with an uncharged anisotropic fluid sphere, in principle the outer manifold is described by the vacuum space-time ı.e, external Schwarzschild solution. Nevertheless, the inclusion of the θsector into the matter field could in principle modify both the geometry and the matter content of the outer manifold. In this situation the compact structure will be not immersed in vacuum space-time anymore. In this case the compact object can remain embedded into a vacuum space-time whether the contributions coming from the θ -sector are assumed to be confined within the stellar interior only [90]. Thereby the outer manifold M + is given by [131] being M Sch the Schwarzschild mass which coincides with the total mass M of the object at the boundary . Then, to match the inner manifold M − with the external one M + we apply the Israel-Darmois junction conditions procedure [132,133]. This procedure dictates the continuity of the first and second fundamental form across the surface . The first fundamental form says that the intrinsic geometry described by the metric tensor g μν induced by M − and M + on the interface meets g − tt | r =R = g + tt | r =R and g − rr | r =R = g + rr | r =R (65) and the second fundamental form related with the extrinsic geometry described by the extrinsic curvature tensor k i j (Latin indices run over on spatial coordinates) says and The expression (66) is obtained from the continuity of k rr while (67) comes from the continuity of k θθ and k φφ .

Junction conditions model 1
Next from Eqs. (48), (64) and (65) one obtains and from expressions (45), (50) and (66) After some algebra Eqs. (68)- (70) can be combined to provide Equations (71)- (73) are the necessary conditions to determine the constants parameters that characterize the model. In Table 1 (25), (44), (45), (85) and (50) one obtains From (74) one gets or in order to assure a positive defined density throughout the compact configuration. Now, from (75) and Zeldovich's condition p r (0)/ρ(0) ≤ 1 one has Based on the previous discussion the inequality (78) is true iff or On the other hand, if the θ -sector is turned off ı.e, α = 0 which implies F = 0 the seed solution is recovered. In particular, for this model from Eq. (24) it is observed that C must be a strictly positive quantity in order to have a physical relevant density ρ describing a compact object. Of course, if α and F are zero the values of the constants A and B will change in magnitude but they still have the same sign. So, in conclusion despite the solution was modified the sign of the constant C does not change. This means that in general the deformed solution accepts negative values for the α parameter restricted to the condition (79) or (80), being the resulting solution admissible from the physical and mathematical point of view.

Junction conditions model 2
Similarly to the model 1 we use the boundary conditions as before, obtaining the following equations for the first fundamental form  Table 1 for compact objects 4U 1538-52, SMC X-1 and LMC X-4 αC r 2 e A+Br 2 1 + C r 2 e A+Br 2 α + C e A F + FC 2 r 2 e 2 A+Br 2 and for the second fundamental form one gets In this case we have taken as free parameters the constant B and α. Besides, the mass M and radius R were fixed by using the numerical data corresponding to the compact star  Table 1 for compact objects 4U 1538-52 (top left), SMC X-1 (top right) and LMC X-4 (bottom) Cen X-3. In Table 3 are depicted the resulting values for the remaining constant parameters that characterize the solution.
It should be noted that the constant A depends only on the mass M and radius R. This is because Eq. (81) corresponds to the continuity of the temporal metric potential across the boundary which remains the same under MGD due to the deformation enters via the radial metric potential into the space-time.
From Figs. 1, 2 (upper panels), 5 and 6 (top panels) it is evident that the inner geometry is completely regular throughout the stellar interior and the main salient thermodynamic variables satisfy the aforementioned requirements to describe a realistic compact structure from the astrophysical point of view. Furthermore, the radial pressure p r vanishes at the surface of the structure and the tangential p t one dominates at all points. As it is observed in Figs. 2 and 6 (top panels) both  Table 3 and α = −0.3, 0.1 and 0.3 for compact object Cen X-3  Table 3 and different α = −0.3, 0.1 and 0.3 for compact object Cen X-3 p r and p t drift apart towards the surface inducing an and anisotropic behaviour in the stellar interior. What is more the anisotropy factor increases monotonically towards the boundary remaining finite and continuous in the interior as can be seen from Figs. 2 and 6 (lower left panel). This fact reveals the presence of an attractive force in nature within the compact object. This force helps to counteract the gravitational gradient which sustains the stability and balance of the system against radial disturbances. In Figs. 2 and 6 (lower right panels) it is appreciated that the matter distribution respects causality condition and it is described by a well behaved energy-momentum tensor (see left panels in Figs.  3, 7). The former is an important subject regarding the study of compact structures. Indeed a bounded and finite sound speed of the pressure waves in the principal directions of the fluid sphere says that any sign travelling inside the structure cannot exceed the speed of light. On the other hand the trend of the constraints imposed on energy-momentum tensor are satisfied everywhere. This corroborates that the matter distribution threading the stellar interior is positive defined. So, in conclusion the resulting model could serve to describe realistic compact objects such as neutron stars.

Stability and hydrostatic balance
To complement the above requirements it is also important analyze the stability and hydrostatic equilibrium of the system. To study the former we have considered to analyze the behaviour of the relativistic adiabatic index in the radial direction. This is so because under the presence of local anisotropies a spherically symmetric system is affected only in the radial direction against an eventual gravitational collapse.
In the arena of classical isotropic matter distribution (Newtonian fluid spheres) the collapsing condition corresponds to < 4/3 [40,48]. In the framework of relativistic fluid spheres the situation involves some extra terms which can be seen as corrections to the previous condition [49,50], where ρ 0 , p r 0 and p t0 are the initial density, radial and tangential pressure when the fluid is in static equilibrium. The second term in the right hand side represents the relativistic corrections to the Newtonian perfect fluid and the third term is the contribution due to anisotropy. It is clear from (84) that if we have a non-relativistic perfect fluid matter distribution the bracket vanishes and we recast the collapsing Newtonian limit < 4/3. In this regard, Heintzmann and Hillebrandt [40] showed that in the presence of a positive an increasing anisotropy factor = p t − p r > 0, the stability condition for a relativistic compact object is given by > 4/3, that is so because positive anisotropy factor may slow down the growth of instability. Nevertheless, relativistic correction to the adiabatic index could introduce some instabilities inside the star [134,135]. To overcome this issue in [136] was proposed a more strict condition on the adiabatic index . This condition claim the existence of a critical value for the adiabatic index crit . To have a stable structure, this critical value depends on the amplitude of the Lagrangian displacement from equilibrium and the compactness factor u ≡ M/R. The amplitude of the Lagrangian displacement is characterized by the parameter ξ , so taking particular a form of this parameter the critical relativistic adiabatic index is given by where the stability condition becomes ≥ crit . To compute the adiabatic relativistic index one has the following expression = ρ + p r p r dp r dρ . The other important study is related with the hydrostatic equilibrium under different forces, namely the hydrostatic F h , the gravitational F g and the anisotropic F a forces. To accomplish this analysis we have considered the following modified Tolman-Oppeneheimer-Volkoff (TOV) equation, It is clear that in the case α = 0 the familiar TOV [137,138] in the context of relativistic isotropic models is recovered. From Figs. 4 and 8 it is observed that the system is in equilibrium under the mentioned forces. As said earlier the anisotorpic force is repulsive in nature. Then the gravitational gradient is counteracts by the action of the hydrostatic and anisotropic forces. This prevents the system to collapse below its Schwarzschild radius onto a point singularity. Moreover, after some point inside the stellar interior the anisotropic force dominates the hydrostatic one showing the preponderance that local anisotropies have in the balance on the configuration.

Concluding remarks
In this paper, we have combined the class I approach and gravitational decoupling methodology by means of minimal geometric deformation to generate anisotropic interior solutions. The main component of this technology is that we have implemented the class I condition as generator of minimal geometric deformation functions f (r ). To obtain the decoupler function f (r ) we have tested our approach by setting two different class I space-time solutions, namely Adler-Finch-Skea solution (model 1) and Kuchowicz solution (model 2). As it is well known in the framework of gravitational decoupling the θ -sector is completely determine once the ν(r ) metric potential and the decoupler function f (r ) are specified. Although the generating equation (41) of decoupler functions looks quite complex, it can be solved once the seed geometry has been proposed. This approach offers a new possibility to find deformation functions f (r ) in the light of gravitational decoupling in addition with previous works [90,[92][93][94][95][96] where the mimic constraint scheme was used to generate f (r ) or imposing a suitable f (r ) as was done in [103,119].
It is worth mentioning that we have tested the viability of the resulting models to describe compact objects supporting by an anisotropic matter distribution, analyzing all the necessary criteria that any admissible compact structure must satisfy in order to represent a realistic neutron or quark star. In this concern we have checked the behaviour of the geometric structure as well as of the main thermodynamic variables inside of the compact configuration, causality condition, the trend of the energy-momentum tensor, stability by means of relativistic adibatic index and hydrostatic balance by using the modified Tolman-Oppenheimer-Volkoff equation. All this analysis is supported by Figs. 1, 2, 3, 4, 5, 6, 7 and 8 for both, model 1 and 2 where it is clear that the obtained solutions satisfy all the requirements to be an acceptable models capable to describe realistic compact structures. Moreover, as Tables 1, 2 and 3 exhibit we have used real data to obtain the constant and physical parameters that characterized the solution. Specifically, we have imposed the mass and radius for some known compact configurations, namely 4U 1538-52, SMC X-1 and LMC X-4 for model 1 and Cen X-3 for model 2.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There are no external data associated with this manuscript.] 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 .