Self-bound embedding Class I anisotropic stars by gravitational decoupling within vanishing complexity factor formalism

We present a spherically symmetric embedding Class I solution for compact star models using the gravitational decoupling approach. We have chosen a null complexity factor condition proposed by Herrera (Phys Rev D 97:044010, 2018) in the context of a self-gravitating system and derive the anisotropic solution through a systematic approach given by Contreras and Stuchlik (Eur Phys J C 82:706, 2022). In this regard, we use the Finch–Skea model along with the mimicking of mass constraint to find fluid pressure and the matter-energy density from the Einstein Field Equations (EFE). We tested the physical viability and impact of gravitational decoupling on the anisotropic solution through the graphical representation. Moreover, the energy exchange between the fluid distributions along with the mass-radius ratio of different compact objects has been also discussed.


Introduction
Despite research over several years by numerous researchers, the exact nature of the compact stars is still somewhat unknown. With the rapid advancement in observational techniques, theoretically modeling these objects have been ever so important. Speaking of observational techniques, in recent times, fruitful insights about compact stars have been found by studying the accreting systems, where the matter is supa e-mail: 21434093@uofn.edu.om b e-mail: sunil@unizwa.edu.om (corresponding author) c e-mail: 27716095@uofn.edu.om d e-mail: 14313615@uofn.edu.om e e-mail: 14648319@uofn.edu.om f e-mail: 19212282@uofn.edu.om g e-mail: gmustafa3828@gmail.com plied to the compact star by a donor star, and the observed radiation is generated by accreted matter. The temperature increases gradually as the matter comes near the compact star. Thus the inner parts of this flow emit high-energy X-rays and UV radiation is produced from the outer parts. Alongside this, one can think of the recent detection of gravitational waves as well. For example, the properties of some nuclear parameters (see [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]) have been constrained by observations by NICER [18,19] of X-ray emissions coming out of hot spots on the surface of rotating neutron stars, observation of pulsars having mass 2M [20], LIGO and VIRGO [21] observations of gravitational waves coming out of two colliding neutron stars. For probing strong-field gravity, these neutron stars are great sources of their compactness. Several modifications to GR have been influenced by observations of binary neutron star merger GW170817 [22][23][24][25][26][27] and binary pulsars [28][29][30][31][32][33][34][35][36]. It's particularly fascinating to think about what goes on inside these objects, as they have nuclear matter in tremendously high densities and pressure. So how the fluid distribution may behave is a very important question to answer. Some of these compact objects have very strong magnetic fields, for example, pulsars and several other spinning stars. It can be assumed that these objects consist of elementary particles such as baryons, mesons, and leptons. It must be noted that the simultaneous and precise measurement of the mass and radius of a neutron star is quite difficult. The maximal mass value of a neutron star is still an open topic of research. Some observational studies estimated this limiting value as ≈ 2M . Some recent studies have found neutron stars having different masses such as PSR J1614-2230 (∼ 1.97M [37]), Vela X-1 (∼ 1.8M [38]), and 4U 1822-371 (∼ 2M [39]). However, LIGO/VIRGO collaboration reported a gravitational wave signal originating from the merging of a black hole and a compact object, whose nature is unknown. This is known as the GW190814 event which invoked debate among the scientific community about its origin. The identification of such an object can be done in several ways, such as, by looking at the mass, checking its tidal deformability or, analyzing the electromagnetic signal originating from the same binary responsible for the gravitational waves. However, for the GW190814 event tidal deformability could not be measured because of inadequate information, and the electromagnetic counterpart was also absent. The mass was found to be 2.5−2.67 M , which is a bit higher than the expected maximum mass limit of a neutron star. Although it was argued by several researchers that this object can be a black hole, in a recent paper of 2020, it was demonstrated that this secondary compact star is likely to be a neutron star [40]. Later it will be shown that in the current article, the maximum value of mass was found to be ∼ 2.64M , which backs up the theory that the unknown object of the GW190814 event can actually be a neutron star. In astrophysics, one of the main objectives is to find out the configuration and geometry of inshore and intramural substances of compact stars. Another fascinating research area is the formation of compact stars and the phenomenon of gravitational collapse. The first exact interior solution for symmetric stars was obtained by Schwarzschild in 1916 [41]. The gravitational collapse was explored by Oppenheimer and Snyder in 1939 [42]. So far, there are only a handful of interior solutions of Einstein's field equations, considering perfect fluid distribution. This shows how difficult it is, to find interior solutions to Einstein's field equations because of high non-linearity. But considering anisotropy often makes it more convenient to find interior solutions. But the question is, can the matter distribution at all be anisotropic inside a compact star? In this context, there has been a detailed analysis in a paper by Herrera and Santos [43], where it was shown that the matter distribution in compact stars is most likely anisotropic, and considering anisotropy is more realistic. In fact, even if the matter distribution was initially isotropic, eventually it may become anisotropic [44]. The idea of anisotropy in compact stars was initially proposed by Ruderman [45] in a review article on the dynamics of pulsar's structure. Anisotropy in compact stars can generate for various reasons, such as very high-density core, different condensate states such as pion and meson condensate states, rotational motion, relativistic particles, superfluid 3A, different fluids intermixed and so on. Gravitational tidal effects also seem to be resulting in anisotropy, as it causes deformation and hence results in anisotropy. So it can be said that in extremely dense matter distribution, the pressure breaks into radial and tangential components, resulting in anisotropy ( = p r − p t ). Several researchers have been studying compact stars considering different matter distributions including anisotropy [46][47][48][49][50][51][52][53][54][55][56][57][58][59][60][61].
As it was stated earlier, obtaining solutions to Einstein's field equations is a very difficult task, researchers have adopted several techniques to make the problems easier to solve mathematically. In recent times, one particular methodology has been very popular amongst researchers trying to find solutions to field equations. This is known as the gravitational decoupling method. It was first incorporated by Ovalle [62,63], in form of minimal gravitational deformation (MGD). MGD technique involves deformation along the radial component of the metric tensor. However, the MGD was further extended to incorporate transformation along both radial and temporal components of the metric tensor. The extended version of MGD is known as CGD. One great advantage of gravitational decoupling (GD) is that it allows the generalization of a simple case, and it can also be used to simplify a relatively complex problem. For example, a simple problem can be generalized by introducing an arbitrary source that can include many physical features including anisotropy, charged distribution, and so on. The first system is the Einstein system and the second system is the quasi-Einstein system. Then the field equations for both systems are then solved separately and the linear combination of the solutions yields the final solution. In order to solve the quasi-Einstein part, several more approaches such as the mimic approach, considering some well-acknowledged solution, and so on can be incorporated, which may make the problem more interesting. In a similar fashion, a complex system can be split into a combination of two sources, one being the Einstein system and another being the quasi-Einstein system and then the same methodology can be followed. The reader can go through several works by various researchers where GD has been effectively used in different contexts to obtain interior solutions for compact stars .
In this paper, another concept called the vanishing complexity condition has been used. For this, let us shed light on what is meant by "complexity" in the first place. The understanding of complexity can appear very arbitrary and subjective. This is the most challenging factor to quantify the concept such as complexity. For example, if we think in terms of information, a perfect crystal has minimum information so it can be thought of as the least complex system. On contrary, ideal gas has random motion and has maximum information. But, if we think about minimum disequilibrium, ideal gas has minimum disequilibrium whereas, disequilibrium is predominant for the perfect crystal. So one way of overcoming this issue is to express complexity as the product of information and disequilibrium. However, for the self-gravitating system, the much accepted and used definition of complexity was given by Herrera [94] where the complexity factor for static spherically symmetric fluid distributions was proposed. The notion of this concept of complexity is to take into account the complexity that may arise in a fluid's structural properties as well as its pattern of evolution. The definition of complexity factor given by Herrera [95] is dependent on active gravitational mass or, Tolman mass. This scalar quantity named the "complexity factor" depends on pres-sure anisotropy and energy inhomogeneity. So, the vanishing complexity can be achieved in two scenarios, such as, when both pressure anisotropy and energy inhomogeneity vanishes, which is the most basic scenario. Another scenario can be when pressure anisotropy and energy inhomogeneity cancel one another. So the vanishing complexity condition is very much dependent on the properties of the system itself. Several researchers have used the definition given by Herrera in various circumstances [96][97][98][99][100][101][102][103].
In the current work, anisotropic solutions of Einstein's field equations describing embedding Class I compact stars have been obtained by using a vanishing complexity factor condition in the context of the gravitational decoupling method. The Finch-Skea model along with mimic to mass function approach [104] is applied to solve the system of equations. In this regard, some works on gravitational decoupling in the context of the complexity domain can be seen in the following Refs. [105][106][107][108][109][110][111].
The paper is organized as follows: In Sect. 2, a brief review of Einstein's field equations under GD and complexity in self-gravitating systems have been discussed. Section 3 deals with complexity-free embedding Class-I solution via MGD. The energy exchange between the relativistic fluids for the sources θ i j and T i j have been analyzed in Sect. 4. Section 5 deals with the junction condition and in Sect. 6, a detailed physical analysis has been done. The concluding remarks are given in the Sect. 7.

Brief review of Einstein's field equations under gravitational decoupling and complexity in self-gravitating system
Let us write Einstein's field equations by including of extra source through a coupling constant β by assuming geometrical units G = c = 1. The line element for describing the spacetime geometry for the internal structure of the stellar model can be given by static spherically symmetric metric as where ν(r ) and λ(r ) are functions depending only on the radial coordinate r . Then we assume that the internal structure of the fluid model contains the anisotropic matter distribution governing by the energy-momentum tensor (T i j ) given by, where the index i = 0, 1, 2, 3. The tensor component u i is called the fluid's four-velocity vector while χ i denotes the radial unit space-like vector defined as u i = (e −ν/2 , 0, 0, 0) and χ i = (0, e −λ/2 , 0, 0), (4) for which χ i u j = 0 and χ i χ j = −1.
Hence, the Einstein field Eq. (1) for a spherically symmetric line element (2) for anisotropic matter distribution is given by and Bianchi identity gives the conservation equation, The mass function m(r ) can be determined as, Now our primary objective is to solve the system of Eqs. (5)-(7) containing eight unknowns, which is not an easy task due to its non-linear nature. Therefore, we apply a condition for the vanishing complexity factor given by Herrera [95] for the field Eqs. (5)- (7), On inserting pressure and density components in Y T F , we get After integration and some manipulation, we obtain the following relation [104], where A and B are integration constants. Also, we split the system of equations into two systems by using transformations in order to reduce the unknowns in the equations, where f (r ) = 0 and g(r ) = 0 are called deformation functions. After plugging the above transformations into the Eqs. (5)- (7), we get the following governing equations, and Moreover, we consider spacetime describing the solution of the system of Eqs. (16)-(18) as, whose internal mass can be obtained from, Then from Eq. (9), we can write where, Now we solve both systems individually. The first system contains five unknowns which will be solved by using the embedding Class I condition. This condition can be given by Karmarkar condition [112], subject to R 2323 = 0 [113], where the quantities given in Eq. (26) are called the Riemann components which will be calculated under the metric (22). The following differential with Y = 1. After integrating, we get the relation between X (r ) and Y (r ) where C andD are constants of integration the integration constant.

Complexity-free embedding class I solution via MGD
In this section, we discuss the solution of the equation of motions (5)- (7) under the gravitational decoupling satisfying the complexity factor condition. For this purpose, we must solve both systems of Eqs. (16)- (18) and (19)- (21). In order to find the solution for the system of Eqs. (16)-(18), we need two auxiliary conditions to close the system. For this purpose, a well-behaved metric function corresponding to the Finch-Skea solution, where L denotes constant with dimension length −2 along with the mimic constraint for the mass function proposed by Contreras and Stuchlik [104], have been employed. Now substituting Y in condition (28), we get while the condition (30) gives which derives the expression for g(r ) as, Then deformed metric function e −λ(r ) is calculated as and the gravitational mass of a deformed star is After substituting potential λ(r ) into the condition (13), we obtain Hence, the other deformation function f (r ) is determined by using the transformation (14) as, After substituting the expression for ν(r ) and X (r ) in above Eq. (37), we obtain and corresponding deformed spacetime is given by, where, e ν (r ) is given by Eq. (37). Now we define the effective components of density and pressures as, whose expressions are given as, 8π where, A systematic diagram for generating the anisotropic solution from perfect fluid distribution using complexity factor and gravitational decoupling in GR is shown in Fig. 1.
Before moving next section, it is necessary to mention an important development given by Herrera et al. [114] in connection to a solution-generating technique to construct all possible types of static spherically symmetric solutions to the Einstein field equations for the anisotropic fluids distributions. They proved that this can be given in terms of two generating functions. By inspiring that methodology, we discover the above two generating functions for the present gravitationally decoupled complexity-free solution. To obtain this, let us analyze by using Eqs. (6) and (7), one finds To simplify the calculations, we present new variables as Substituting Eqs. (45) into the Eq. (44), we get where (r ) = 8 π (p eff r − p eff t ). Integrating Eq. (46) we obtain (r ) as This general formula was given by Herrera et al. [114] to find all possible solutions in terms of two generating functions (r ) and (r ). Then, the general spacetime can be written using (45) and (47) as, Therefore, in our present case the Eq. (45) turns out to be: because of the resulting algorithm, the generating functions for the present complexity-free anisotropic fluid solution under the gravitational decoupling solution are as follows [using the Eqs. (42) and (43) along with Eq. (36)]:

Energy exchange between the relativistic fluids for the sources T i j and θ i j
The idea of the energy exchange concept under extended gravitational decoupling was discussed by Ovalle [63]. Later on, this concept was introduced by Ovalle et al. [124] and Contreras and Stuchlik [125] for the first time to discuss the flow of energy exchange between the sources. It was noted that a successful decoupling of both systems under the CGD approach requires the energy exchange ( E) between both sources T i j and θ i j . Therefore, it is necessary to discuss the exchange between both sources which can be determined by the formula by, Further details for the development of this formula can be seen in the above-mentioned references. On inserting the temporal deformation function along with the seed pressure and density in the above formula, we find an expression for E as,

Junction condition
The essential boundary conditions will be discussed in this section in order to obtain the constants involved in the complexity-free anisotropic solution. In this regard, we match the spacetime for the interior solution i.e. complexityfree anisotropic solution with the Schwarzschild exterior spacetime at the boundary of star r = R. Here, M denotes the total mass of the object at r = R and then M = m(R). The necessary and sufficient boundary conditions are derived by Darmois-Israel boundary conditions [119,120] which are given by, We find the constants after utilizing of above boundary conditions as 6 Physical analysis

Regularity conditions
It is already known that any physically viable regular solution must have finite central pressures and density in order to ensure a singularity-free solution at the center. In the present complexity-free model, the central pressures and density are given as: where, i ≡ {r, t}, F 1 = A(β + 1) 2 tan −1 1 √ β . It can be observed from Eqs. (61) and (62) that both pressures and density are free from the singularity. Furthermore, for modeling any realistic object, it is required that pressure and density must be positive everywhere inside the object and the pressure-density ratio must be less than one. Using these bounds on the pressures and density, we get the following constraints on the constant L, On the other hand, the p eff (0) ρ eff (0) ≤ 1 provides After combining the inequalities (63) and (64), we get, The regular behavior of the density, pressures and anisotropy within the compact objects can be observed via density plots (Figs. 2, 3, 4 and 5). The energy density (ρ eff ) is seen to be maximum at the center and falls off gradually in the radially outward direction (Fig. 2). This means the core has more energy density than the surface. From Fig. 3 it can be seen that the radial pressure ( p eff r ) has the maximum value at the center and decreases towards the surface. It is also seen that the radial pressure approaches zero near the surface. The tangential pressure ( p eff t ) is also maximum at the center and decreases outwards and is minimum at the surface (Fig. 4). However, p eff t does not approach zero near the surface and keeps a non-zero positive value. So, the radial pressure approaches zero near the surface and tangential pressure remains non-zero positive value, which means the pressure distribution satisfies physically acceptable stellar objects. Also, as both the pressure components grow from the surface towards the center, it indicates positive hydrostatic force acting in the outward direction. The gravitational pull is balanced by the hydrostatic force and which in turn helps avoid gravitational collapse, which is one of the fundamental requirements of realistic compact objects. It must be noted, for β ≈ 0.35, the tangential pressure doesn't monotonically decrease, instead, it slightly increases with r initially, then starts decreasing. This anomaly suggests, when gravitational decoupling is not in effect, the model is not completely well-behaved. The variation in the anisotropic factor shows that the anisotropic factor ( eff ) is zero at the center and it gradually increases and is maximum at the surface (Fig. 5). This means the anisotropy is positive and thus the repulsive anisotropic force helps maintain stability inside the star. Looking into the dependence of β, it is seen that both p eff r and p eff t increase with the increase in β, but the energy density (ρ eff ) and anisotropic factor eff seem to be relatively less impacted by the influence of β as both ρ eff and eff slightly increases with the increase in β.

Causality condition
For every model to be physically acceptable, it has to follow the causality condition, i.e. neither of the radial and tangential components of the velocity of sound should exceed the speed of light. If the speed of light is scaled to unity, the causality condition becomes 0 < v 2 r < 1 and 0 < v 2 t < 1 where v r and v t are the radial and tangential components of the velocity of sound respectively. The analysis of the causality condition has been done via density plots in r − β plane (Fig. 6). Here, it is noted that v 2 t becomes negative for β ≈ 0.35 in the range 0 ≤ r ≤ 4 which is not physically acceptable. So, this shows, for β ≈ 0.35, or, when gravitational decoupling is absent, the model is not physically well-behaved and the presence of gravitational decoupling stabilizes the model. Barring the above range, it can be seen that both v 2 r and v 2 t are well within the range 0 < v 2 r < 1 and 0 < v 2 t < 1, which means, the causality condition is satisfied everywhere else. Both radial and tangential velocities increase radially, i.e. they have a minimum value at the center and increase as one travels towards the boundary. Also, it is recorded, that β increases both v 2 r and v 2 t . The maximum value of both v 2 r and v 2 t are recorded near the surface when β = 0.60.

Fig. 4
The contour diagram for equi-tangential pressure on the r − β plane for L = 0.0095 km −2

Adiabatic index
Another key aspect to check while checking the stability of any model is to check the radial component of the adiabatic index (γ ). This becomes even more significant when the fluid distribution has the presence of anisotropy. This is because the presence of anisotropy may result in a departure from the spherically symmetric distribution in order to escape gravitational collapse. The collapsing condition for perfect fluid distribution with non-relativistic distribution is γ < 4/3 [115,116]. However, when the matter distribution is anisotropic, this changes the form [117,118] to, The second and third terms of the r.h.s represent correction due to relativistic consideration and anisotropy respectively. It is obvious from the above equation that for nonrelativistic perfect fluid distribution, the previous relationship is retrieved. Moreover, as anisotropy increases, the instability decreases, so the stability condition further gets updated to γ > 4/3 considering the anisotropic factor ( = p t − p r ). The Equi-adiabatic index, as seen in Fig. 7 contour diagram in r − β plane is seen to be slowly increasing with r up to some point and then it rapidly increases near the boundary. Also, the condition of γ > 4/3 is always fulfilled everywhere within the model. Moreover, the adiabatic index (γ ) is seen to be very slightly increasing with β.
Due to the consideration of the relativistic matter, some instabilities may occur inside the compact stars [121,122]. To resolve this issue, another condition was imposed by Moustakidis [123], where a critical value of γ was introduced (γ crit ). For a model to be stable, γ ≥ γ crit condition must be satisfied. The γ crit is defined as, γ crit = 4 3 + 19M 21R . From Table  1, it was found, that the condition of γ ≥ γ crit is not satisfied for lower values of β, and from β ≥ 0.55, this condition is satisfied. So it suggests, that gravitational decoupling leads more stable model as compared to pure GR.

Mass function and compactness
The effect of gravitational decoupling on the mass m(r ) is discussed here. For a detailed clarification, we start by analyzing the consequences of gravitational decoupling on Einstein's gravity theory. As mentioned above, the mass of the object in the context of gravitational decoupling is the sum of seed mass (m G R ) and mass (m θ ) for an extra source via decoupling constant β. Since we find the deformation function using the condition (31), then the mass of the deformed spacetime becomes Long back Buchdahl [126] proposed the maximum allowable mass-radius ratio constraint M G R R ≤ 4 9 for isotropic matter distribution satisfying the decreasing density. Here, we check the effect of decoupling constant β on the total mass (M) and whether the object satisfies the Buchdahal limit under gravitational decoupling. It is obvious that the total mass of the object under gravitational decoupling will be (1 + β) times the total mass in the pure GR scenario i.e. for fixed radius R, M ≥ M G R will always hold true when β ≥ 0 and consequently the upper bound will change in presence of gravitational decoupling. The graphical analysis of mass is done by a contour diagram of equi-mass in r − β plane which shows how the mass is varying with β along the radius r . From Fig. 8 it is recorded that mass (in terms of M/M ) is zero at the center and it increases as one moves towards the surface. This means the compact object has a kind of mass distribution where it is more concentrated towards the boundary and near the core has a lower mass. Also, it is noticed that M/M also increases with the increase in β, which means, β allows the star to be more massive. The   Table 3 Mass-radius ratio of 11 observed compact stars corresponding to Table 2  Moreover, we have also predicted the radii for 11 different compact objects for different values of decoupling constant β in order to predict the influence of the gravitational decoupling on the predicted radii. These predicted radii are given in Table 2. In this regard, Table 3 shows the mass-radius ratio (M/R) to the same 11 different compact objects for different β. We can observe that the mass-radius ratio (compactness) increases when β increases but the obtained compactness values for each object lie within the range given the Buchdahal. On the other hand, the increasing value of compactness shows that the gravitational decoupling provides more compact objects as compared to the model in a pure GR scenario. and e each deformed compact object.

Surface redshift
Furthermore, the surface redshift z s of the compact objects under the gravitational decoupling can be given from the following formula, The Eq. (68) clearly shows that the decoupling parameter β affects the surface red-shift value due to the change in compactness u = M R . This was well supported in the result found in Table 1. Here for fixed radius, it was found that z s changes for different values of mass. So it essentially varies with β and as M/R ratio gets changed due to change in β, z s also varies. The surface redshift is studied via a contour diagram in r − β plane and it is observed that (Fig. 9) z s increases as one moves from the center to the surface. Moreover, it doesn't vary much with β. However, it seems to be slightly increased with the increase in β. The highest value of z s is attained near the surface and at β = 0.60.

Energy exchange
Here, we present the energy exchange between the relativistic fluids for sources θ i j and T i j based on the discussion proposed in the works [124,125]. They have done the analysis of the energy distribution between fluid distributions according to the positive and negative values of E. The negative value of E shows the source θ i j takes the energy from the environment (seed fluid) while E > 0 indicates that the new source gives the energy to the seed fluid or environment. In order to find the flow of energy exchange, we plotted the equicontour diagram for energy exchange ( E) on r − β plane against the free constants ratio A1/B1 governing the seed solution. The contour diagrams for E are plotted in Figs. 10 and 11 for A1/B1 = 1, 2 − 0.03, and − 1. The left panel (A1/B1 = 1) of Fig. 10 shows that E is negative on r − β plane within the object i.e. when r ∈ (0, 10). This shows that the source θ i j takes the energy from the environment while there is no energy exchange between the sources at the center and boundary of the star. It is highlighted that there is a huge amount of energy exchange near the core when β < 0.40. If we see the right panel (A1/B1 = 2), the same situation is happening as the left panel but the amount of energy exchange is less as compared to the left panel. This implies that when A1/B1 increases then the amount of energy exchange from the environment to the new source decreases. Furthermore, when 0.35 ≤ β ≤ 0.37 and 0.5 ≤ r ≤ 1.5, then the value of E = −1.0 at A1/B1 = 0.1 and the value of E is more negative when A1/B1 → 0. However, the value of E is ruled out for some range of r near the core, and this range shifts towards the boundary when A1/B1 < 0, which can be observed from the left and right panels of Fig. 11. As we can see from the left panel of Fig. 11 when 0.2 < r ≤ 1.8, then E is positive and E is negative for r > 2.5 for all 0.35 ≤ β ≤ 0.60 while the value of E is ruled out for 1.8 < r ≤ 2.5. The right panel of Fig. 11 shows the same trend for different ranges of r . If 0.0 < r ≤ 6.2, then E is positive while E is negative for r > 7 for all 0.35 ≤ β ≤ 0.60. However, the value of E is ruled out 6.2 < r ≤ 7. By comparing of left and right panels of Fig. 11, we observe that if the value of A1/B1 moves left from 0, then E has both positive and negative values for different ranges of r which implies that the new source gives

Concluding remarks
In this paper, a new interior solution for compact stars has been obtained by solving Einstein's field equations, using the gravitational decoupling method alongside the vanishing complexity condition, which provides a bridge equation between metric functions. To solve this equation, we use the Finch-Skea model along with the mimic-to-mass functions approach. After obtaining the solutions, a thorough physical analysis has been done including thermodynamic properties, causality, and stability analysis in terms of adiabatic index, mass, mass-radius ratio, surface redshift, and energy exchange. To begin with, the thermodynamic properties were studied. The energy density ρ eff has a high value at the core and then decreases along the radially outward direction (Fig. 2). It essentially means that the core has much more energy density than the outer portion of the compact star. The radial and tangential pressure was studied and it was found that both radial and tangential pressure have their maximum value at the center and then gradually decrease in the outward direction (Figs. 3 and 4). However, at the surface, the radial pressure approaches zero, while the tangential pressure remains non-zero and positive. The study of the anisotropic factor ( eff ) revealed that it is zero at the center and then keeps increasing as one travels toward the boundary (Fig. 5). So it can be said that the core has high energy density, low anisotropy, and the highest magnitude of pressure, and as one moves away from the core, the energy density and magnitude of pressure components decrease, and anisotropy increases. Now, let's look into how β impacts the thermodynamic parameters, which would suggest the impact of grav-itational decoupling. It is noted that increase in β increases both p eff r and p eff t significantly, and ρ eff , eff mildly increase with the increase in β. Also, for β ≈ 0.35, some anomaly in the nature of tangential pressure is found, as it doesn't monotonically decrease along the radius, instead, it slightly increases up to some point and then starts to decrease again. Another anomaly for β ≈ 0.35 is found while studying the causality condition (Fig. 6). The tangential component of velocity (v 2 t ) becomes negative for β ≈ 0.35 in the range 0 ≤ r ≤ 4, which is not physically relevant. This shows, when the gravitational decoupling is not present, the model is not well-behaved. And as gravitational decoupling takes place, the model stabilizes and becomes well-behaved. Apart from β ≈ 0.35 and in the 0 ≤ r ≤ 4, both radial and tangential velocity components satisfy the causality condition 0 < v 2 r < 1 and 0 < v 2 t < 1. If we look into the nature of the curves of v 2 r and v 2 t , it can be noticed that both of them have minimum values at the center and increase in the outward direction. Also, the increase in β increases both radial and tangential velocities.
The stability in terms of the adiabatic index (γ ) was studied and it was found that the stability condition γ > 4/3 is maintained throughout the model (Fig. 7). This means, the model is stable in terms of the adiabatic index, which is a very essential thing to check, especially when the fluid distribution is anisotropic. Also, γ is seen to be slightly increasing with r up to some point, and then increases drastically towards the surface. β is seen to be slightly increasing the value of γ . Table 1 showed that the stability condition based on the critical value of the adiabatic index is not satisfied for lower values of β and only gets satisfied for β ≥ 0.55 onwards. This essentially shows that the model is more stable at higher values of β, i.e. when the effect of gravitational decoupling is predominantly more.
The study of mass via a contour diagram of equi-mass (Fig. 8) revealed that mass is minimum at the core and gradually increases towards the surface. So this essentially suggests that the value of mass is more on the surface and minimum in the center, where the energy density is maximum. Moreover, it was noticed that M/M increases with the increase in β. The maximum value of mass was found to be 2.6M and it is achieved near the surface for β = 0.60. Table 1 revealed that the mass-radius ratio (M/R) for 11 different compact stars increases with β. This shows, the compactness increases with β, so it's evident that more compact objects are supported by gravitational decoupling as compared to pure GR. Even though the compactness increased with β, it was always found to be staying within the given range by Buchdahal.
The surface redshift (z s ) was studied and it was found that z s slightly increases with β. This was also suggested by Eq. (68) that z s gets affected by β for the changes in compactness u = M R . The variation of z s shows, it increases with r and its highest value is obtained near the surface at β = 0.60 ( Fig. 9).
The energy exchange ( E) was studied in Figs. 10 and 11 for A1/B1 = 1, 2 − 0.03, and − 1. From the left panel (A1/B1 = 1) of Fig. 10, it was found out that energy is being taken from the environment by the source θ i j for r ∈ (0, 10), while at the center and boundary, the energy exchange between the sources is non-existent. Notably, it was found that when β < 0.40, there is a great amount of energy exchange near the core. For A1/B1 = 2 plot, in the Fig. 10 right panel, it was seen that the scenario is quite similar to the one before, only here the energy exchange is lesser than the previous case. So, it suggests that the increase in A1/B1 decreases the energy exchange from the environment to the new source. If we look at Fig. 11, we observed different scenarios under the different range of r when A 1 /B 1 is negative. In this case, the new source gives the energy to the environment as well as takes energy from the environment for different ranges of r while some value of E is ruled out for some range of r .
Finally, we would like to mention here that the complexityfree conditions along with the gravitational decoupling approach open a new tool for generating the series of solutions which is not limited to GR only but also applicable to modified gravity theory.
Acknowledgements All authors acknowledge that this work is carried out under TRC Project (Grant No. BFP/URG/CBS/22/133). The authors are also thankful for continuous support and encouragement from the administration of University of Nizwa.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The current study is developed for theoretical stellar models and no novel data is generated. The unique parametric space used in the article to produce the plots is stated in the text.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.