A simple protocol for anisotropic generalization of Finch–Skea model by gravitational decoupling satisfying vanishing complexity factor condition

A physically reasonable anisotropic stellar model is constructed with the help of the gravitational decoupling via complete geometric deformation (CGD) technique under the condition of vanishing complexity factor [Contreras and Stuchlik in Eur Phys J C 82:706 2022; Herrera, in Phys Rev D 97:044010, 2018]. The source splits into a perfect fluid and an anisotropic distribution. The Finch Skea metric proves a useful seed solution to solve the Einstein sector while the condition of vanishing complexity is invoked to solve the remaining anisotropic system of equations. A comprehensive battery of tests for physical significance is imposed on the model. Through a careful choice of parameter space, it is demonstrated that the model is regular, stable, and contains a surface of vanishing pressure establishing its boundary. Matching with the exterior metric is also achieved. Finally, the energy flows between the two sectors of the source fluid are studied graphically.


Introduction
The most successful theory of gravity is Einstein's general relativity (GR), which after decades of careful scrutiny, can a e-mail: sunil@unizwa.edu.om (corresponding author) b e-mail: abdelghani.errehymy@gmail.com c e-mail: mahmoodkhalid@unizwa.edu.om (corresponding author) d e-mail: hansrajs@ukzn.ac.za e e-mail: nfhariby@uqu.edu.sa f e-mail: amabdelaty@ub.edu.sa describe a broad range of phenomena at solar system scales and beyond. This theory extends Newton's notion of gravitational force and describes how the mass-energy distribution causes the geometry of space-time to deform. Indeed, a plethora of new discoveries and extremely precise experiments, including SNeIa, LSS, CMB and BAO [1][2][3][4], reveal that the expansion of the universe is accelerating. Furthermore, the unknowable nature of the majority of the universe's substance was further proven by various recent cosmic observations at astronomical scales [5,6].
A deeper understanding of the dynamics of the universe may also be gained by studying large-scale structures such as stars, galaxies, and their clusters. The intricate nature of these stellar structures has a significant effect on associated physical variables including energy density, pressure, and heat flow. In this context, researchers have concentrated on the formation and properties of these stellar structures, which are commonly applied to refer to white dwarfs, neutron stars, and black holes. It was long ago realized by Schwarzschild, who obtained the first exact solution to Einstein's field equation for the interior of a compact stellar structure. The most likely compact stellar structures have been found in pulsars, which are rapidly rotating stars with intense magnetic fields made completely of quark matter. Eventually, Rosat surveys [7], performed in 2006, identified compact stellar structures based on their X-ray emission, revealing that the gravitational energy generated is emitted in X-rays. Hewish et al. [8,9] revealed the hypothetical emergence of pulsars, which generate a beam of electromagnetic radiation that is truly continuous but pulsed. Therefore, the discovery of pulsars and X-ray bursts motivates researchers to construct theoretical models of compact stellar structures such as neutron stars and quark stars. Several models are occasionally put out to characterize compact stellar structures, despite the fact that the makeup of the particles and the nature of interactions are still unknown.
Currently, considerable efforts are being expended on studying relativistic compact stellar structures. We investigate the static spherically symmetric Einstein's equations solutions for object modeling under various physical circumstances. These solutions can be described as dust, ideal fluid, and anisotropic fluid. In this context, the majority of the approaches frequently employed by many researchers to study these novel analytic solutions of the field equations impose symmetry requirements such as spherical symmetry, an equation of state relating the pressure and energy density of the stellar fluid, the behavior of the pressure anisotropy or isotropy, vanishing of the Weyl stresses, spacetime dimensionality, and so on, in the attempt to explore exact solutions characterizing static compact stellar structures [10][11][12][13][14][15][16][17]. In this regard, the concept presented by Riemann known as Riemannian geometry is also a useful tool for analyzing the basic geometrical properties of compact stellar systems. These hypotheses make the challenge of obtaining appropriate solutions to the field equations easier and mathematically solvable. Moreover, the relaxation of the isotropy condition boosts the probability of discovering accurate solutions, with the caveat of incorporating imperfect fluids. It is generally recognized that deviations of the isotropy and fluctuations of the local anisotropy in pressures can be induced by a wide range of physical phenomena that we would anticipate to occur in compact stellar objects (see Refs. [18], for an in-depth discussion on this point). In addition, the presence of physical factors such as dissipative fluxes, and/or energy density inhomogeneities, and/or the emergence of shear in the fluid flow will always tend to generate pressure anisotropy, even if the system is originally claimed to be isotropic. These physical factors were pointed out by Herrera in 2020 [19]. However, unequal principle stresses, also termed as anisotropic fluids, might be expected when the selfgravitational fluids' densities are typically higher than the density of nuclear matter. The concept of anisotropy arises in self-gravitational compact stars due to the occurrence of exotic phase transitions [20], electromagnetic fields, rotations, superfluids or type-A fluids [21], pion and meson condensations [22], core formation and other phenomena have also been investigated extensively by Herrera et al. [18,[23][24][25]. This suggests that the self-gravitational systems have two different kinds of pressure components, namely the radial component ( p r ) and the tangential component ( p t ). Consequently, radial and tangential pressures become unequal ( p t = p r ), and the concept of local anisotropy emerges in the study of self-gravitational fluids. In this context, Herrera and collaborators [13,[26][27][28] explored static anisotropic stars by looking into the effects of Newtonian and general relativistic regimes.
A plethora of exact solutions are known in the literature for nonstatic, radiating stars, including acceleration-free collapse, Weyl-free collapse, vanishing of shear, collapse from/to an initial/final static configuration, and anisotropic collapse models. According to the Buchdahl constraint of ideal fluid distributions, the compactness factor u = 2M/R should always be ≤ 8/9 in order to avoid gravitational collapse and the formation of a singularity or black hole. Andréasson [29,30] generalizes this upper constraint in u by including charge, anisotropy, and even cosmological constant. He considered a straightforward inequality relating pressure and density, p r + 2 p t < ρ, as the starting point for deriving the new generalized upper constraint in u. More impressively, a generic method described in [31][32][33][34][35][36][37][38][39][40][41] can be attractively applied to find an algorithm for all static and anisotropic solutions with and without charge to Einstein's equation for the spherically symmetric line element.
A novel method for generating anisotropic solutions to Einstein's field equations was recently developed [42][43][44][45][46][47]. This novel method known as gravitational decoupling via minimal geometric deformation (MGD) which was intended to enable solutions generated by an isotropic matter distribution to reach anisotropic domains. This approach deforms the object's geometry while simultaneously modifying its material content, conserving the symmetry of the solution. The number of studies employing this methodology that are now available in the literature has significantly increased during the recent past. To highlight a few instances, applications include studying stellar interiors, black holes, and modified gravity theories 78,79]. Furthermore, in [56] the MGD inverse problem has been developed, i.e. knowing the isotropic counterpart of an anisotropic solution.
The MGD approach aids in the study of the key properties of compact stellar configurations, however, it has significant limitations. For example, geometric deformation can only be achieved if the interplay involving matter sources is purely gravitational. The metric component's minimal transformation, which only affects the radial coordinate while maintaining the temporal metric potential as an unchanging entity, is likely it has some drawbacks, in particular, it can not explain a stable black hole having a well-defined horizon. In this regard, by utilizing the deformation on both radial and temporal metric functions, Casadio et al. [46] proposed an extended version of the MGD technique and generated a novel solution for spherically symmetric spacetime to overcome this problem. However, the conservation law does not apply in the presence of matter, hence this extension is limited to studying vacuum solutions. Therefore, this extended approach does not allow for analysis of the interior structure or intrinsic characteristics of self-gravitating objects. By manipulating both metric functions (g rr , g tt ), Ovalle [51] put forth the innovative concept of extended geometric deformation or complete geometric deformation (CGD), which remains valid throughout all of the spacetime regardless of the selection of matter distribution. In order to extract the exterior charged BTZ model from the appropriate vacuum solution, Contreras and Bargueno [77] effectively decoupled the field equations in (1 + 2)-dimensional gravity using the CGD approach. In this connection, some notable works using CGD approach can be seen in Refs. [80][81][82][83].
For several years, the notion of complexity of a system has been the topic of extensive study in all fields. Different factors are involved in analyzing the complexity of any system. The main idea is to measure the entropy and information of the configuration contained inside a system. The concept of the complexity of the self-gravitating system is widely analyzed in the studying of heavy configurations. Whereas in physics, isolated gases (which reveal disorder and the greatest quantity of information) are somehow sophisticated systems with vanishing complexity when considering a perfect crystal (which exhibits periodic behavior and is symmetrically dispersed). The concept of disequilibrium was developed by López-Ruiz et al. [84] to analyze the complexity of the system. In general, it is a measurement of "distance" from the system's achievable form's equally likely dispersion. They argued that the concept of complexity dissipated in the scenario of an ideal gas and perfect crystal by considering complexity as a mixture of both concepts of "disequilibrium and information". Next, Herrera [85] established a novel concept of complexity based on fluid constituents including energy density, pressure, and others after observing deficiencies in existing notions of complexity when studying the self-gravitating system. In summary, it is connected to all-inclusive aspects of the fluid's content. With the help of the complexity factor, which is one of the structural scalars obtained from the orthogonal division of the intrinsic curvature, the complexity is in this case constructed. This concept for dissipative fluid content was extended by Herrera et al. [86]. They did more than only analyze the system's complexity; they also established the prerequisites for the progression design with the least amount of complexity. They explored that there are several solutions and that the fluid is shearing and geodesic in terms of dissipation.
On the other hand, Herrera et al. [87] were using the axially symmetric geometry to further investigate the effect of complexity on various geometries and identified three different categories of complexity. They showed how complexity and symmetry are interrelated. In this peculiar case, they also obtained some analytical responses. The evolution of spherically symmetric non-static geometry is also investigated by Herrera et al. [88] using the concept of complex-ity, both in terms of dissipation and non-dissipation. They constructed some models and determined their appropriate implications for understanding evolution by adopting the quasi-homologous condition, which is a relationship between areal radius velocity and areal radius. By using this approach, Contreras and Fuenmayor [89] examined the stability of selfgravitating configurations in terms of gravitational cracking and performed a thorough study to look into the effects of the source's compactness and changes to the decoupling parameters on the radial force. Taking into account the Misner-Sharp mass and the Tolman mass as well as the involvement of various structural scalars gained via the orthogonal division of the intrinsic curvature, Herrera et al. [90] extended the concept of the complexity factor on the geometry with hyperbolic symmetry. They came to the conclusion that Tolman's mass in this scenario reflects negative nature. Furthermore, many researchers have effectively employed this complexity approach not only in GR [91][92][93][94][95], but also in the analysis of findings for various geometries in the context of modified theory of gravity [96][97][98].
In this paper, we develop a simple protocol for an anisotropic generalization of the Finch-Skea model by gravitational decoupling satisfying vanishing complexity factor condition. We are using the complete geometric deformation (CGD) scheme, which transforms the gravitational potentials g rr and g tt for exploring the physical influence of the generic source i j on seed source T i j and decomposes the system of non-linear field equations into two arrays. One set of these arrays corresponds to the seed source, and the other set provides extra source terms. We would also highlight that the energy transfer between fluid distributions corresponding to the new source ( i j ) and original source (T i j ) were analyzed based on the positive and negative values of the energy exchange E. To be more specific: (i) if E is positive, then the new source is supplying energy to the environment, (ii) if E is negative, then the original/seed matter distribution is supplying energy.
The article is organized according to the following arrangement. In Sect. 2, we derive the basic field equations for a static sphere containing two sources which are decoupled through the CGD technique. by using the well-behaved Finch-Skea metric potential for the seed spacetime geometry, which ensures a well-defined horizon-free spacetime, we are able to obtain a gravitationally decoupled solution in Sect. 3. In Sect. 4, where we match the decoupled interior solution determined by the anisotropic matter distribution to the exterior Schwarzschild solution at an appropriate boundary, the external spacetime and matching conditions are addressed. The behavior of density, pressures, and anisotropy inside the self-gravitating stellar system, stability of anisotropic solution using adiabatic index, Harrison-Zeldovich-Novikov stability analysis, and energy exchange, which are each covered in Sects. 5.1, 5.2, 5.3, and 5.5, were all considered in Sect. 5 to establish how each solution behaves physically and how it is viable and stable. In Sect. 6, the concluding remarks are given. The Appendix now contains some pertinent lengthy expressions of physical quantities.

Einstein's field equations in the framework of gravitationally decoupled system
In this section, we consider a static fluid with spherical symmetry. This fluid is hypothesized to be anisotropic and constrained by the surface designated by the symbol . The analogous line element is represented in Schwarzschild-like ansatz as where ν(r ) and λ(r ) are functions depending only on the radial coordinate r . Furthermore, the above-mentioned metric complies with the gravitationally decoupled Einstein field equations provided by where in which a new source induced by gravitational decoupling is denoted by the symbol i j . We consider that the space's physical content consists of an anisotropic distribution of matter with energy density , radial pressure P r , and tangential pressure P ⊥ , in order to characterize and explain the underlying structure of the self-gravitating system associated with the source T i j . Furthermore, the matter sector may be expressed by the following energy-momentum tensor, where the index i = 0, 1, 2, 3. Moreover, u i is called the fluid's four-velocity vector, and χ i denotes the radial unit space-like vector that satisfies such that ξ i u j = 0 and χ i χ j = −1. Following that, using the spherically symmetric line element (1), the energymomentum tensor components become as Thus, the gravitationally decoupled Einstein field equations (2) are stated as where indicates the derivative with respect to the r coordinate. Next, using the conservation law ∇ i T i j = 0 or by employing Einstein's equations (7)-(9), we can readily derive the generalized Tolman-Opphenheimer-Volkoff (TOV) (hydrostatic equilibrium) equation [99,100] for anisotropic matter distribution which reads dp eff Alternatively, using the formula for ν as Here, the mass function m(r ) for the spherically symmetric distribution can be defined by or, equivalently [101] Furthermore, it is necessary to define the nature of original energy-momentum tensorT i j whether this describes either perfect fluid matter distribution or anisotropic matter distribution. Here, we assumeT i j denotes a perfect fluid matter distribution with spacetime geometry μ and ξ and matter variables ρ (energy density) and p (pressure), is defined aŝ with α i (four-velocity vector) and ζ i are given by where α i ζ i = 0 and ζ i ζ i = −1.
Then effective quantities are written as and the effective anisotropy depends on, It is clearly noticed from the above equation that anisotropy can only be introduced in the system when β = 0. Now, we employ the extended gravitational decoupling approach [51] to modify the gravitational potentials e λ and e ν for exploring the physical influence of the generic source i j on seed source T i j . The metric functions e ν and e λ are effectively deformed via a linear transformation given by where f (r ) and h(r ) stand for the geometric deformation functions subject to the radial and temporal coordinates, respectively. The coupling constant β is a real number. Furthermore, the aforementioned transformation includes both the radial and temporal components of the line element and is the extended case of minimal geometric deformation (MGD), also known as an extended geometric deformation or a complete geometric deformation (CGD). The CGD requires that both deformation functions be non-zero, i.e., f (r ) = 0 and h(r ) = 0. Therefore, using these transformations, we get the two sets of equations. The first set of equations for the seed source is whose solution can be described by the following line element, The second set of equations for the extra gravitational source can be given by, The corresponding hydrostatic equilibrium equations take the following forms dp r dr and − d 1 The mass function m s is determined by Which implies Moreover, it should be mentioned that the energy exchange between the two sourcesT i j and i j is necessary, to succeed in decoupling between them. The equivalent energy exchange between these sources, denoted by E, may be provided by [51] Next, the gravitationally decoupled mass function m(r ), as defined by Herrera [85], may be expressed in terms of the homogeneous energy density and the change caused by density inhomogeneity as Then from Eq. (19), we can get the following relation, where The formulation of the mass function for a spherically symmetric static spacetime characterizing the energy content of a fluid stellar structure was also introduced by Tolman [103] as It may be also expressed 1 as Now we move on to the process of finding the solution for both systems, which requires some points that need to be clarified about what we have so far: (i). We have two systems that depend on the nine unknowns {ξ, μ, f, h, ρ, p, 0 0 , 1 1 & 2 2 } with five independent equations. Therefore, we need four auxiliary conditions to solve the system completely. , only two auxiliary conditions must be specified in order to close the system of differential equations. These two conditions could be: The mimic constraints for the pressure ( 1 1 = p r ), density ( 0 0 = ρ), an equation of state (EoS) operating between the -sector, or any other wellknown condition. The first three conditions are widely applied to build the interior solutions under the gravitational decoupling system. Furthermore, it is mentioned that the solution of differential equations under the constraints 1 1 = p r is more difficult here because this mimic constraint incorporates the high non-linear term in h and involves both functions f and h together. Therefore, this condition is not suitable to build the solution. Due to this, we need a particular EoS that gives a linear differential equation and involves only the radial deformation function. The best condition is the mimic constraint of the density i.e ρ = 0 0 approach. (iii). Finally, we need another condition to close the second system completely. For this purpose, we use the recent condition known as the vanishing complexity factor condition proposed by Herrera [85]. In accordance with Herrera's concept [85], we define the complexity factor (Y T F ) formula in the context of Einstein's field equations (7)-(9) system as Then, using Einstein's set of Eqs. (7)-(9), we can express, Y T F , as Then vanishing of the complexity factor condition, i.e., In terms of total derivative, the aforementioned equation may thus be expressed as which leads to log ν where log A 1 denotes the integration constant. In this case, the simplified aforementioned equation becomes ν e ν/2 = A r e λ/2 .
According to the aforementioned equation, ν and λ have the following relationship e ν = A 1 re λ/2 dr + B 1 2 (44) with B as an integration constant. It is worth mentioning that the Eq.(43) enables us to derive the deformation function h in terms of the radial deformation function f and the seed metric functions {ξ, μ} using equations (19) and (44) as [102] In order to accomplish this, we must identify the radial deformation function by a particular EoS that only includes f and its derivatives i.e. the mimic constraint of the density (ρ = 0 0 ) approach, as discussed earlier. It is emphasized that the choice of the seed solution influences the analytical solutions of Eq. (45). In light of the reasoning above, we will employ the following process to obtain the anisotropic interior solutions in the presence of vanishing complexity under gravitational decoupling via extended MGD as: (i). Select the most appropriate seed metric functions: ξ and μ. (ii). Start executing the mimic constraint for the density based on a specific EoS ρ = 0 0 . (iii). By substituting μ and f into Eq. (45), one may get the solution for the variable h. once we get f and h, then we can find e −λ = μ + β f and ν = ξ + βh directly, which is our final aim.
Now we are going to present the whole approach mentioned above to find the new anisotropic solution to Einstein's field equations (7)-(9) having zero complexity by taking the wellknown Finch-Skea perfect fluid model and energy exchange between relativistic matter distributions in the next section:

Anisotropic generalization of Finch-Skea perfect fluid model
When there are no handy evidences for the source and nature of particle interactions, one must develop a model of the stable content of relativistic compact stellar configurations by formulating the analytical solutions of Einstein's field equations characterizing the static internal core of relativistic stellar configurations. However, due to the highly non-linear second-order field equations, generating the exact solutions to the gravitational field equations is not at all an easy task.
To overcome the issue in this case, a variety of appropriate techniques are frequently applied. In this regard, we specify an explicit metric function corresponding to Finch and Skea [104] spacetime which was employed to model the interiors of relativistic stellar configurations. In this connection, when Duorah and Ray [105] first developed this type of ansatz in 1987, they did not guarantee that it would be fine-tuned to satisfy Einstein's field equations governing the content of stellar relativistic astrophysical models. Such a form of Finch-Skea solution has aroused great curiosity in the construction of relativistic compact stellar stars since the findings have been well-proven and validate all of the basic and adequate requirements of physical viability [106]. The astrophysical and cosmological backgrounds in GR and as well as higher-order modified gravity theories have been discussed in several outstanding explorations relating to this ansatz [107][108][109][110][111][112]. Motivated by the preceding discussion, we consider the spacetime geometry corresponding to Finch-Skea perfect fluid solution as We can easily get the formula for thermodynamic variables like density and pressure by employing these metric functions μ(r ) and ξ(r ) as ρ(r ) = C Cr 2 + 3 Let us now proceed to the subsequent stage, where we first apply the appropriate EoS, which only contains the deformation function f (r ) as an unknown function, in order to identify the deformation function f (r ). The mimic constraint for density is, mathematically, the most straightforward EoS involving the deformation function f (r ), This gives rise to the differential equation shown below The foregoing equation yields the solution for the deformation function f (r ) as where D is an arbitrary constant of integration. This arbitrary integration constant D was selected to be zero in order to achieve a non-singular solution and fulfill the constraint f (0) = 0. The new metric potential, denoted by the symbol e λ(r ) , is then known as the deformed metric potential and may be expressed as Now, plugging the Eq. (53) into Eq. (44) and integrating, we arrive at the following form of potential e ν , e ν(r ) Now the generalized spacetime geometry for the Finch-Skea model having vanishing complexity factor can be given as where, ν(r ) is given by Eq. (54). Now from Eq. (45), we can directly find the second deformation function h(r ) as +2 Ay(r ) Cr 2 + 1 sin The flow chart of this simple protocol for generalizing the perfect fluid solution to the anisotropic domain under the vanishing complexity factor is mentioned in the Fig. 1.

Matching conditions
Besides the above, on the boundary surface r = R, the interior metric should be smoothly connected to the exterior metric, hence we need the continuity of the first and second fundamental forms over the boundary surface. The exterior spacetime is described by the exterior Schwarzschild solution, which is (60) In order to achieve the above requirements, the first and second fundamentals forms (Israel-Darmois junction conditions [113,114]) are mathematically expressed as, The matching of two spacetime metrics and vanishing of effective pressure ( p eff ) at the boundary surface r = R, the aforementioned requirements (61)-(63) are not independent. These are however sufficient to determine all necessary integration constants involved in the solution. Using the above conditions, we find the expressions for the following constants,

Physical behavior of anisotropic solution
It is generally recognized that a theoretically well-behaved compact star model must meet some fundamental physical and mathematical criteria. The salient characteristics of the anisotropic stellar solution discovered by incorporating the well-known Finch-Skea perfect fluid model with vanishing complexity factor are very helpful in describing the composition of the relativistic compact star. This will be employed to analyze the necessary criteria in the subsections as follows: 5.1 The behavior of physical quantities such as density, pressures, and anisotropy inside the self-gravitating stellar system In this subsection, based on the graphical plots shown here, we give an in-depth physical analysis of our findings with an emphasis on viability related to the anisotropic generalization of Finch-Skea contributions satisfying the vanishing complexity factor via gravitational decoupling. Figure 2 illustrates the energy density's behavior. In this graph, we can see how the density profile changes as the deformation parameter β rises. The density, as we can observe, is a monotonically decreasing function with respect to the radial coordinate, r . It is evident that the density of the compact stellar configuration grows gradually when the deformation parameter, β, is increased from 0.05 to 0.18. The radial pressure at every inside point of the stellar configuration is shown in Fig. 3. As one goes from the core towards the stellar surface, the radial pressure falls off smoothly until it vanishes at the boundary, as we would expect since there is no energy flux to the surrounding spacetime. The effect caused by the deformation parameter β demonstrates that the radial pressure increases as β grow. The same finding holds true for the tangential pressure with a little change when β assumes tiny values. Additionally, we note that, throughout the stellar configuration, the tangential pressure consistently outweighs its radial counterpart with very small deviations. In Fig. 4, we present the behavior of the anisotropy parameter with respect to the radial coordinate, r . We note that at every inside point of the stellar structure, the anisotropy, eff , is positive. A repulsive force comes from anisotropy when the tangential pressure dominates the radial stress. By preventing the stellar configuration from being pulled inward by the gravitational force, this repulsive force aids in stabilizing it. Additionally, we emphasize that the degree of anisotropy can be governed by the deformation parameter β, where an increase in β is accompanied by an increase in anisotropy, eff . Generally, we see that when one gets closer to the compact stellar configuration's surface layers, the anisotropy is the greatest.

Stability of anisotropic solution using Adiabatic Index
We are mainly interested in discussing the stability of the compact stellar object model using the adiabatic stability criterion defined by Chandrasekhar for isotropic pressure gradients (see [115,116]). The formula for this adiabatic stability criterion is = 1 + ρ p dp dρ S , where dp dρ is the sound speed and the subscript S indicates a constant specific entropy. It establishes that > 4/3 for stellar configurations with isotropic pressures, p. It has been shown by Herrera et al. [117,118] that this criterion modifies when pressure anisotropy is involved, and takes the following form, where differentiation with respect to radial coordinate, r is shown by the prime. The Newtonian limit, < 4 3 , for unstable areas, emerges from the vanishing of the second component in (67), which originates from relativistic contributions. The adiabatic index can be modified by radial heat flux dissipation or the existence of density inhomogeneities. In this context, stability versus radial perturbations is maintained whenever > crit , where the critical value for the adiabatic index, crit is defined as crit = 4 3 + 19 21 u [119], and u = M/R indicates the stellar model's compactness. The behavior of the stability criterion for our model is shown in  Table 1 for all chosen values of β. It is noteworthy to note that an increase in the deformation parameter β tends to stabilize the stellar configuration as can be ascertained from Fig. 6. This demonstrates that our model meets the Chandrasekhar stability criterion and is stable under radial adiabatic infinitesimal perturbations, since the adiabatic index ( ) is growing and greater than 4/3 at the stellar surface for all chosen values of β.

Harrison-Zeldovich-Novikov stability analysis
We also need to investigate the Harrison-Zeldovich-Novikov (HZN) stability criterion in order to confirm the stability of the stellar model. To that end, we analyze the stability of the stellar model corresponding to the relationship between total mass M and central energy density ρ eff 0 at the point when d M/dρ eff 0 = 0. According to the HZN stability criterion's definition, [120,121], which reveals to be satisfied along with every relativistic compact stellar object. Figure 7 demonstrates that, for our model, d M showing that we have a stable stellar configuration when the deformation parameter is moving from 0.05 to 0.18. In contrast, we show the stellar configuration in Fig. 8 as a function of stellar mass M and the effective central density ρ eff 0 , which increases with increasing values of β in accordance with the gradient profile shown in Fig. 7. Furthermore, the stellar configuration on the segments d M/dρ eff 0 > 0 is always stable

Compactness and surface redshift
According to the formula u(r ) = m(r )/r , the compactness factor [122,123] of a stellar configuration is the ratio of the active gravitational mass to the constraining radius. The upper bound for the maximum possible value permitted for the compactness factor is defined by Buchdahl [124], which is u(r )| max ≤ 8/9 for a gravitationally confined spherically symmetric fluid. In this connection, Fig. 9 describes the u(r ) profile; it is intriguing to notice that u(r ) rises with the equilibrium radius, r, for each branch of the deformation parameter, β ∈ [0.05, 0.18], and u(r ) stays significantly under the Buchdahl bound. The surface redshift, Z s (R), at r = R is then derived using the radial component of the line element as, The redshift spectrum is shown in Table 1, it is clearly seen that the surface redshift increases significantly with any increase in β and also Buchdahl bound [124][125][126], Z s < 2 holds for the surface, r = R for each branch of the deformation parameter, β ∈ [0.05, 0.18]. Furthermore, we plotted the equi-mass contour diagram (10) on r −β plane to observe the effect of β on mass with r . As we can see from Fig.10 that if we fix radius R between 0 to 7.5 and increase β, then no effect in mass is observed. But if R ≥ 7, then we found that mass (M/M ) is increasing with β and the maximum increment is found near the boundary. On the other hand, if Finally, we conclude that the decoupling constant (β) introduces an extra packing of mass to the stellar models.

Energy exchange
In this section, we discuss the energy exchange between the sources T i j and i j . Now using Eq. (31), we find the expression for Energy exchange, with E 1 (r ) = 1 − A 2 Cr 2 + 1 cos where constant A 2 denotes A/B. The interesting physical feature of an anisotropic solution is the exchange of energy between the relativistic fluids corresponding to the new source ( i j ) and original source (T i j ). The energy transfer between fluid distributions can be analyzed based on the positive and negative values of E. More specifically: (i) If E is positive, then the new source is giving energy to the environment. (ii) If E is negative, then the original/seed matter distribution is giving the energy.
Since the constant A 2 is independent of the gravitationally decoupled solution but it appears in expression E. Therefore, it is necessary to see the effect of this constant on the Energy exchange along with the decoupling constant β. For this purpose, we plot Figs. 11 and 12 to observe the transitions of the energy exchange between the sources on the β − r plan for different values of constant A 2 . The left panel of Fig. 11 is plotted for equi − E contour diagram on the β − r plane for A 2 = 0.5. It can be observed that the higher value of E lies between 7 ≤ r ≤ 10 and 0.02 ≤ β ≤ 0.03, which shows that the new source is giving a high amount of energy to the environment in that region. But when we move near the core i.e. 1.5 ≤ r ≤ 2.5, the E achieve the highest negative value for the same range of β. This implies that the perfect fluid matter distribution is giving a high amount of energy. On the other hand, if we look at the right panel of the Fig. 11, we observe that the higher value of E is shifting to 3 ≤ r ≤ 6 for β = 0.21 when we increase A 2 from 0.5 to 1. Furthermore, the positive range of E is shifting towards the boundary when A 2 is increasing. Now we move to Fig. 12, which has been plotted for the negative values of A 2 = −0.5 (left panel) and A 2 = −0.75 (right panel). From the left panel, we can observe that E is positive between 0.0 < r ≤ 8.5 for all values of β ∈ [0.02, 0.2] which implies that the new source is giving the energy with the above range. But when r ≥ 8.5, E stats negative, and values of E are ruled out after r ≈ 9.8 at β = 0.02. But when β increases, the negative value of E starts shifting near to boundary but the highest negative value exists when β ≤ 0.065. This implies that the perfect fluid matter distribution is giving a high amount of energy near the boundary for all values of β ≤ 0.065. Furthermore, the highest positive value of E lies between 3 < r ≤ 6 at β = 0.02. On the other, the right panel of this figure shows that the pattern of energy exchange is the same as the left panel but the negative value of E is ruled out for the wide range of r near the surface.

Concluding remarks
We now summarise the principal findings of our investigation. A physically viable static anisotropic spherically symmetric stellar model was constructed. The method of gravitational decoupling of the metric components enabled the splitting of the source into a standard Einstein system and an additional system of equations. It was deemed prudent to make use of the well-studied Finch-Skea metric to satisfy the Einstein sector which was assumed to be a perfect fluid. A further condition was necessary to fully determine the new source with pressure anisotropy. To accomplish this it was assumed that the fluid distribution displayed a vanishing complexity factor in the sense of Herrera. The even- tual model was carefully examined to check if it satisfied stringent regularity and stability conditions. In particular, the adiabatic stability criteria of Chandrasekar as well as the Harrison-Zeldovich-Novikov condition were all found to be satisfied. Graphical plots using specified parametric values verified that the mass-radius relationships comported with known physical behavior. In addition, a study of the energy exchange between the standard fluid and the anisotropic was exhibited in the form of plots. Moreover, all integration constants that appeared along the way were settled by matching the interior spacetime with the exterior Schwarzschild metric. This experiment has demonstrated the value of the decoupling approach in devising astrophysical models that harmonize with observed phenomena.