An anisotropic stellar fluid configuration with vanishing complexity

In this work, we construct a new stellar model in the regime of anisotropic fluid pressure using the concept of vanishing complexity for spherically symmetric fluid distributions (Herrera in Phys Rev D 97:044010, 2018) and a convenient ansatz in order to close the Einstein’s field equations. The resulting model fulfills the fundamental physical acceptability stellar conditions for a specific set of compactness factor. The stability and its response against fluctuations in the matter sector is also investigated.


Introduction
In 1915, Einstein published his famous work on the General Theory of Relativity, which changed completely our view and understanding of the structure of space and time in the universe. The concepts related to this theory are realized in the Einstein's Field Equations (EFE), which describe how gravity arises from the curvature of space-time in the presence of matter or radiation. Since then, several works have been devoted to studying these equations and the construction of a wide range of models for describing phenomena from the solar system to cosmological scales. Thus, Schwarzschild [1] obtained the first EFE solution describing the space-time in the interior and exterior of a compact uniform-density sphere, and ever since, the modeling of relativistic stellar compact objects has moved from the regime of toy models to sophisticated, realistic stellar structures. For example, it is worth mentioning the seminal work of Tolman [2] in obtaining exact solutions of EFE for static fluid spheres and the pioneering work of Bowers and Liang [3] where the local anisotropy for static, and spherically symmetric distribution of matter is analysed. In this context, researchers have concentrated a e-mail: julio.andrade@espoch.edu.ec (corresponding author) b e-mail: diego.santana@espoch.edu.ec on the formation and properties of these stellar structures, which are commonly used to refer to white dwarfs, neutron stars, and black holes. Thus, despite the fact that the complete nature of interactions inside of these stellar compact objects is still unknown, several models are published in order to characterize compact stellar structures.
Researchers consider both isotropic and anisotropic pressures on the relativistic fluid that supports the compact objects in order to develop models that can describe their main characteristics. To first approximation the local pressure isotropy can be considered a valid assumption in the study of stellar compact objects, however, it is known that many physical processes can produce deviations from isotropy and/or local anisotropy fluctuations in pressure that may be caused by a large variety of physical phenomena, especially in compact objects. Theoretical studies on more realistic stellar models in [4] show that nuclear matter may be locally anisotropic at very high densities. According to these views, the radial pressure may not be equal to the tangential pressure in massive stellar objects. Some of the research in this regard is reported in [5][6][7][8][9][10][11][12]. Additionally, it has been found strong evidence suggesting that for certain ranges of density, a large number of physical phenomena can cause local anisotropy and must therefore be considered in order to describe more realistic models. For example, one possible source of anisotropy is related to the intense magnetic fields observed in compact objects such as white dwarf stars, neutrons or magnetized quark stars [13][14][15], pion and meson condensations [16], exotic phase transitions [17] and others phenomena are extensively investigated in [18][19][20][21][22].
Even so, a recent result in [23] shows that any realistic physical process of the type expected in stellar evolution always tends to produce pressure anisotropy due to the presence of dissipation, energy density inhomogeneities, and shear. As a result, any final stage of a dynamical regime in the evolution of a star should exhibit pressure anisotropy. Con-sequently, radial and tangential pressures become unequal, and the concept of local anisotropy inevitable emerges in the study of realistic compact stellar models [24][25][26][27].
It is worth noting that the problem of obtaining exact EFE solutions for static and spherically symmetric space-times supported by anisotropic fluid includes solving a robust system of coupled differential equations. The problem is that there are only three independent EFE, but five unknown quantities: two metric functions, the energy density, the radial and tangential pressures (matter sector). This issue has been addressed over time by employing widely successful strategies for obtaining precise EFE solutions [28][29][30][31][32][33][34][35][36][37][38]. Even exact solutions of the EFE for anisotropic fluid distribution on the background of a pseudo-spheroidal or paraboloidal spacetimes have been obtained in [39][40][41][42][43]. And for example, recently, the polytropic equation of state was used for developing models for compact objects with anisotropic polytropes [11,44,45], as well as modified polytropic equations of state are used to obtain anisotropic stellar models in [46,47], which result versatile in the way to describe several applications in astrophysics such as white dwarfs, neutron stars or Fermi fluids. Furthermore, it is worth noticing that in the Ref. [48] a novel approach to integrate the Lane-Emden equations for relativistic anisotropic polytropes is presented.
Particularly, in this work we obtain a new interior solution of EFE in the anisotropic regime of pressure as a result of the using as an extra condition the idea complexity for static self-gravitating spheres proposed by Herrera [49] and a convenient ansatz for the radial metric potential λ. This idea of complexity for the self-gravitating system has been widely accepted and used; it arises from the existence of a structure scalar called complexity factor Y T F that contains information about the matter content of the fluid distribution. Specifically, this factor depends on pressure anisotropy and energy inhomogeneity. It is worth mentioning the case of the called "vanishing complexity", where Y T F = 0, which can be reached when both pressure anisotropy and energy inhomogeneity vanish, which pertains to the simplest case. When pressure anisotropy and energy inhomogeneity cancel out, another more fascinating scenario emerges. This idea of complexity for self-gravitating spheres has been used as an extra condition to close the EFE in order to obtain new stellar models. Thus, we can mention several works where this definition has been used; for example, recently, anisotropic solutions of EFE describing embedding Class I compact stars have been obtained by using a vanishing complexity factor condition in the context of the Gravitational Decoupling approach (GD) (see Refs. [50][51][52][53][54][55] for details about GD approach) [56], as well as the construction of an anisotropic generalization of the Buchdahl static stellar model by implementing the method of GD via extended minimal geometric deformation and further requiring vanishing complexity in [57]. Also, it is worth mentioning that the role of GD on isotropization and complexity of selfgravitating systems under the complete geometric deformation approach and the role of complexity on self-gravitating compact stars under GD have been studied in [58,59], respectively. Even this definition has been used for the construction of new stellar models in the framework of alternative theories of gravity (see, for example, the Refs. [60][61][62][63][64][65][66][67][68][69][70][71][72][73]); for example, the vanishing complexity and gravitational decoupling approach have been used to investigate a spherically symmetric anisotropic solution in f (Q) gravity theory for the first time recently in [74].
In this work, we also use the definition of vanishing complexity for self-gravitating spheres in order to close the resulting EFE. The obtained solution is regular and meets all physical plausibility conditions, and its stability is also investigated. It is worth mentioning that the number of interior solutions in the regime of anisotropic fluid with vanishing complexity is reduced nowadays; therefore, new models with these characteristics are important in the field of theoretical astrophysics, and therefore the relevance of obtaining this kind of new solutions is well established.
The paper has been organized as follows: Sect. 2 is dedicated to reviewing briefly the EFE for anisotropic fluid distribution, in Sect. 3 we introduce and review the idea of complexity for self-gravitating fluid spheres, in Sect. 4 we obtain a new stellar model with the aim of vanishing complexity and a convenient choice of ansatz on radial metric potential, in Sect. 5 we briefly review the fundamental stellar physical acceptability conditions for a realistic stellar compact object, while in Sect. 6, contains a physical analysis of the model. Finally, in Sect. 7, we draw some conclusions about our work.

Einstein field equations for anisotropic fluid distributions
The Einstein field equations (EFE) for vanishing cosmological constant are where R μν the Ricci tensor, R is the curvature scalar, T μν is the energy-momentum tensor and κ = 8π G c 4 . Now, in order to model a star, we must use the space-time of a static self-gravitating sphere provided by the metric where ν and λ are functions (called as metric potentials) that depend only on the radial coordinate r . In such sense in the co moving frame the physical interior of the self-gravitating object can be modelled as an relativistic anisotropic fluid given by a density energy ρ, radial pressure p r and tangential pressure p t . So in this way the energy-momentum tensor T μν is represented by whose components (ρ, p r , p t , p t ) are known as the matter sector of the interior solution and is the four velocity of the fluid and s μ another quantity define as with the properties s μ u μ = 0 and s μ s μ = −1.
So if one uses the information of Eqs. (2) and (3) in EFE (1), one arrives at where primes denote differentiation with respect to the radial coordinate r , and we assume the geometric units G = c = 1. Note that Eqs. (6)-(8) constitutes a set of three differentials equations with five unknown quantities: {ν, μ} (metric potentials) and {ρ, p r , p t }. So in order to solve this system, it is necessary to use two extra conditions, which usually can be geometrical relations (for example, the Karmakar condition [75]), equations of state (EoS) that relate the physical quantities of matter sector (for example, the polytropic equation of state, the Van der Waals equation of state, etc. [76][77][78][79]) and others. In particular, in this work we use the idea of complexity for static self-gravitating fluid spheres of Herrera [49] and a convenient choice of an ansatz for the radial metric potential λ in order to obtain a new interior solution for EFE. Furthermore, the contracted Bianchi identities ensure that the Einstein tensor is divergence-free. Then, by the Eq. (1), we can derive the covariant conservation of the energymomentum tensor as follows: which explicitly gives us the generalized Tolman-Opp enheimer-Volkoff (TOV) equation for anisotropic fluid dp r dr with Π ≡ p r − p t . If one uses the definition of mass function given by or equivalently m = 4π r 0r 2 ρ(r )dr (12) the TOV can be expressed as which accounts for the hydrostatic balance of the fluid within the interstellar object. Note that the pressure gradient is balanced by a gravitational term and a term that includes the local anisotropy distribution.

Complexity of self-gravitating spheres
The concept of complexity varies greatly depending on the subject of study, particularly; in this work, we shall use a definition of complexity for static and spherically symmetric selfgravitational systems that was recently proposed by Herrera [49,80,81] in the context of the general theory of relativity. This definition replaces the fundamental idea of probability distribution that appears in the definition of "imbalance" and information by the energy density of the fluid distribution (see [82]). And likewise, this definition replaces certain previous definitions of complexity for self-gravitating spheres that considered only the energy density of the fluid, ignoring the other components of the energy-momentum tensor that describes the interior of the star from the point of view of general theory of relativity. This definition is based on the structure of the fluid distribution (i.e., fluid inhomogeneity and pressure anisotropy), in such a way that the simplest system is the one with perfect fluid distribution, and that more complex systems are those that vary from this fundamental system, especially those that deviate from the regular pattern of constant energy density and pressure isotropy.
Specifically, such a definition surges from the existence of a structure scalar (denoted by Y T F ) that is connected to the orthogonal splitting of the Riemann tensor [83,84] in static and spherically symmetric space-times (for the first time, such a scalar and others were thoroughly examined in [85]). The Riemann tensor can be expressed through the Weyl tensor C ν αβμ , the Ricci tensor R μν and the curvature scalar R in the following way: and by other hand, we can only express the Weyl tensor in terms of its electric part because the magnetic part vanishes in the spherically symmetric case: It should be noted that E μν can also be written as with having the properties for E μν : It is now possible to show that the Riemann tensor can be expressed using tensors (see [84] for details) as in what is called the orthogonal splitting of the Riemann tensor. Here * denotes the dual tensor, i.e., R * μγ νδ = 1 2 η σ γ δ R σ μν and η μνλρ corresponds to the Levi-Civita tensor. T μν can be expressed in a particularly useful manner so that after some manipulations, namely and with Thus, from the tensors X μν and Y μν is possible to define four structure scalar functions in the following way: and As a result of the preceding, X T F and Y T F determine the local anisotropy of pressure: Also, it is possible to express Y T F in terms of energy density inhomogeneity and system local anisotropy by This scalar captures the concept of complexity since it measures the relationship between the inhomogeneity in the energy density and the pressure anisotropy of a static and spherically symmetric self-gravitational system. Additionally, it can be demonstrated that (33) permits one to express the Tolman mass as where the subscript Σ indicates that the quantity is evaluated on the boundary surface Σ. Equation (34) shows that this scalar includes all the alterations caused by the energy density inhomogeneity and the anisotropy of the pressure on the active gravitational mass, namely, the Tolman mass, which is a combination of its value for a zero-complexity system and two other terms related to energy density inhomogeneity and pressure anisotropy, respectively. It can be viewed as a convincing justification to define the complexity factor by means of this scalar. Thus, this scalar represents an appropriate parameter that characterizes the complexity of self-gravitating static spheres because, first and foremost, it is based on a structure scalar (which is critical because it ensures that this characteristic is founded on a quantity that is invariant for any observer) that contains all physical parameters of the matter sector of the interior of the selfgravitating sphere, and, more specifically, it is dependent on inhomogeneity in the energy density and anisotropy in the pressure. Now, it is worth noticing that if one uses the EFE (6)-(8) in (33) arrives at which represents an alternative method of calculating Y T F using knowledge of the space-time within the stellar compact source. Equation (33), in particular, has been used as a state equation to construct a limited number of new interior solutions for static self-gravitating spheres [56][57][58][59][86][87][88][89][90][91][92][93][94][95]. For example, it is significant to note that the vanishing complexity criterion, namely, when Y T F = 0, is met not only in the most straightforward instances of isotropic and homogeneous sys-tems, but also in the situations where namely, in the scenarios where the pressure anisotropy and energy inhomogeneity cancel each other.
The system of EFE can be closed using Eq. (36) as a complementing condition because it reflects a non-local equation of state (three interesting formalisms to construct solutions with such a characteristic have been developed recently in [91,96,97]).

New stellar model
The vanishing complexity condition, Y T F = 0, yields the differential equation shown below for which is necessary the information of the radial metric potential λ in order to obtain ν. So in this work we propose an ansatz for the radial metric potential given by where A and B are constants. In principle, we propose this ansatz only because we have tested it and it is suitable for us to find a regular and adequate new solution to the EFE, and it is not associated with any particular physical motivation.
In particular, one can associate our ansatz with the particular case of the well-behaved metric function corresponding to the Finch-Skea solution [56,98] given by e −λ = 1 1+Lr 2 where L = B A . Therefore, using this ansatz in (37) we arrive at with c 1 and c 2 integration constants. Now, using the metric potentials (38)- (39) in EFE, we obtain the matter sector given by with A + Br 2 > 0 and c 1 = − Furthermore, by employing the matching conditions (see C6. in the following section), we have with A + B R 2 > 0.

Physical acceptability conditions
The fundamental physical requirements necessary to any interior solution can describe a realistic stellar compact object are (see Ref. [99] for a detailed discussion of these conditions): C1. The metrics e ν and e λ should be finite and regular in the interior of a stellar compact object; e ν should be a monotonously increasing function of radial coordinate r , and e −λ should be a monotonously decreasing function. Furthermore, e ν(0) = constant and e −λ(0) = 1. C2. The matter sector given by the density energy ρ, radial pressure p r and transversal pressure p t should be positive and regular inside the compact stellar object. They should be monotonously decreasing functions of radial coordinate r , with their maximum values at the center. Also, p t (r ) > p r (r ) inside the star, with the exception of the center, where p t (0) = p r (0). Moreover, the radial pressure should vanish at the boundary of the star. C3. The Dominant Energy Condition (DEC): ρ − p r ≥ 0 and ρ − p t ≥ 0 should be fulfilled ∀r, r ≤ R. Also, it is desirable that even the strong energy condition (SEC) ρ ≥ p r + 2 p t is satisfied. Obviously, the latter encompasses DEC. C4. The redshift function z(r ) = g −1/2 tt (r ) − 1 should decrease outward, and its value at the surface z(r = R) is less than the universal bound for solutions satisfying the DEC, namely z bound = 5.211 [100]. C5. The stability of the stellar object requires that the speed of sound should be less than the speed of light, c = 1, which leads 0 ≤ v 2 r = dp r dρ ≤ 1 and 0 ≤ v 2 t = dp t dρ ≤ 1 inside of a stellar compact object. C6. At the boundary of the stellar compact object (r = R), the metric components e ν and e −λ should match continuously with the Schwarzschild exterior solution, namely where M and r Σ = R are the total mass and radius of the star, respectively. As also is necessary that since the exterior of the star in this case is considered to be empty.

Physical analysis
In this section, we analyse the physical acceptance of the stellar model obtained in this work. In this sense, we have checked that the behavior of the whole model given by Eqs.

Metrics
It is noticeable that C1 of the previous section is fulfilled for the model in Figs. 1 and 2, as well, we can notice that e ν is a monotonously increasing function of r and e −λ is a monotonously decreasing function of the same variable r . So the space-time within the stellar compact object behaves as expected.

Matter sector
In Figs. 3, 4, and 5, the behavior of the matter sector of the solution is presented. We can note from these figures that condition C2 is fulfilled. Furthermore, from Fig. 6, we can note that Δ = −Π > 0, and the radial and tangential pressures are equal in the center of the stellar object, which is absolutely necessary in order to maintain certain stability of the stellar object.

Energy conditions and causality
On the other hand, it is also necessary to show that our model should satisfy the dominant energy condition given by C3 in order to affirm that this matter sector belongs to an acceptable physical energy-momentum source. This is effectively observed in Figs. 7 and 8. However, we have checked that the model does not fulfill the Strong Energy Condition (SEC) described in C3 for the compactness parameters 0.35 < u = M/R (see Fig. 9). Based on the foregoing, we can conclude that our model can explain realistic stars from an energy standpoint because it should obey either of the two energy conditions (DEC or SEC) [99][100][101][102]. Moreover, in Figs. 10 and 11 the profiles of sound velocities inside of the star are shown. It is noticeable that these velocities do not surpass the limit value of light velocity in vacuum that in this case is taken as c = 1. However, these results are interesting since their profiles are not monotonous functions of the radial coordinate r .

Redshift
In Fig. 12, the redshift z is plotted. In this figure, we can observe that the redshift is a monotonously decreasing function of the radial coordinate r , having its maximum value at the center of the stellar compact object. Also, it is noticeable We have checked that for the obtained solution fulfills the fundamental physical conditions for those compactness parameters mentioned before. However, since the anisotropy  Δ and sound velocities have a non-monotonous decreasing behavior inside the stellar compact object, it is possible that the model is unstable. As a result, in this work, we go a step further and investigate the stability of the new solution in the next four subsections.

Stability against convection criterion
The stability of a spherical stellar compact object to convection implies the buoyancy principle: any fluid element pushed downward must float back to its original position [103][104][105]. This implies that in any hydrostatic matter configuration, pressure and energy density must decrease outwards. Any fluid that supports a self-gravitating sphere should satisfy this principle in order to remain stable during convection, then in such sense it was demonstrated in [106] that if the energy density of any fluid distribution is such that then effectively it fluid distribution fulfills such principle. Thus, in order to study this condition, in Fig. 13, we show the profile of ρ for our model. From this figure, it is observable that the model fulfills this stability condition in the inner regions of the star, while the outer regions are unstable. To some extent, it could be argued that the outer layers can tolerate this instability while the inner layers can not; to some extent, this could be tolerated. But nevertheless, if the behavior were the opposite, that is, if Eq. (48) is true for the outer layers and not for the inner ones, this would imply that each element of mass tends towards the center of the star, potentially causing a collapse.

Stability against gravitational collapse criterion
To investigate the resistance of the model to collapse, it is necessary to examine the behavior of the adiabatic index Γ in the radial direction given by which should satisfy with This is because a spherically symmetric system is only affected in the radial direction against eventual gravitational collapse in the presence of local anisotropies. The aforementioned relationship given by Eq. (50) accounts for relativistic changes to the adiabatic index, which can cause instabilities within the star. The second term on the right side of Eq. (51) represents relativistic corrections to the Newtonian perfect fluid and anisotropy contribution. This stability condition (50) thus applies to any relativistic compact object supported by an anisotropic fluid. Furthermore, this condition asserts the existence of a critical value for the adiabatic index Γ crit . This critical value is determined by the amplitude of the Lagrangian displacement from equilibrium and u = M/R is the compactness factor (see Refs. [107][108][109] for a more in-depth discussion of this point).
In this sense, we have plotted the behavior of Γ in function of radial coordinate r in Fig. 14. We have checked that the model presents instability against collapse for the compactness factors shown in the caption of Fig. 14.

Harrison-Zeldovich-Novikov stability criterion
Also known as the "static stability criterion", it establishes that any relativistic stellar fluid will be stable and static under radial perturbations if the mass of the system is increasing as a function of its central energy density ρ 0 , namely, d M dρ 0 > 0. The total mass as a function of its central density for this model is from which we obtain that which is of course positive since R > 0. This means that the new model satisfy the static stability criterion. Therefore the anisotropic compact objects described by our new model are stable. However, they are unstable when their internal gases undergo an adiabatic process related to gravitational collapse because the stability against gravitational collapse criterion is not satisfied, which is the first point associated with extreme changes in gravitational force inside stars.
The above result leads us to analyse the behavior of the matter sector of these models against small perturbations in the next subsection.

Gravitational cracking criterion
In this section we study the behavior of the fluid distribution of our model just after its departure from equilibrium, when total non-vanishing radial forces of different signs appear within the system. This analysis is based on the idea of cracking produced in a spherical fluid distribution, which was first developed by Herrera in 1992 [110] and then fine tuned in later works [111][112][113]. The analysis focuses on studying when the radial force in the interior of stellar compact object is directed inward in the inner part of the sphere and reverses its sign beyond some value of the radial coordinate, we say there is a cracking, even if only after the fluid deviates from equilibrium. In the opposite case, when the force is directed outward in the inner part and changes sign in the outer part, we shall say that there is an overturning.
It is important to mention that cracking is closely related to the problem of compact object structure formation, but only at time scales smaller than or equal to the hydrostatic time scale [114][115][116]. In such a sense, the cracking analysis is just take a photography of the state of total radial force within the stellar compact object when it departures from equilibrium. Likewise, it should be clarified that the appearance of fractures directly and drastically affects the future structure and evolution of the compact object.
So to carry out this analysis, it is necessary to note that the TOV (10) or (13) can be rewritten as where R is the total radial force within the stellar compact object per unit of volume on each fluid element. When the stellar system is in hydrostatic equilibrium, this total force R is null, however, when the system is facing some disturbance, this force can deviate from that value. So it is interesting to study those scenarios where these deviations occur.
In such a sense, we shall assume that the perturbation is done for the whole matter sector of the stellar model except that the radial pressure remains unperturbed, namely, we have to change or perturb the parameters of the model in such way that R = 0. The idea is if the model has the parameters {α, β} the "perturbed" TOV can be written as (up to first order) Now, if the systems experiments a cracking or overturning is necessary that at some pointR changes of sign in some interior point of the stellar object, which is traduced in the existence of some real root forR = 0. Thus it can be seen how the existence of some real value Ξ such as δβ = −Ξδα, and moreover that To apply this analysis to our model, we must first define the dimensionless quantities listed below in terms of which we can write for which we also use the condition (43). Now, we proceed to perturb the matter sector through the variation of the parameters {α, β} where the tilde indicates that the quantity has been perturbed. The total radial force, on a fluid element, in terms of {α,β} results in where the first term is zero given that it corresponds to the unperturbed values ofR. Of course, after perturbation, the TOV is different from zero, so the system is no longer in hydrostatic equilibrium. Note that if cracking (or overturning) occurs,R must have a zero in the interval x ∈ (0, 1). In Figs. 15, 16, and 17, we representR as a function of x for different values of Ξ , α and β, respectively. For the realization of these figures we have taken into account the condition A + B R 2 > 0, which can be written as α + β > 0. In Fig. 15, we show the profile ofR with the values of β = 0.7 and α = 1 and different values of Ξ . Note that for the value of Ξ = −0.7 the radial force isR = 0 for every point inside the star, but for the lower values of Ξ the system experiments overturning and for the upper ones the fracture corresponds to cracking. Interestingly, the perturbed total radial force experiments a kind of transition between scenarios with overturning to scenarios with cracking when the value of Ξ increases, and thus the radius where cracking/overturning occurs coincide for all the values of Ξ considered.
On the hand, in Fig. 16, the profile ofR is shown for the values of β = 0.7, Ξ = −0.7 and different values of α. It is noticeable that under this configuration of parameters values the system experiments a transition between configurations with cracking to configurations with overturning when we increase the value of α. Note that the fracture (cracking or overturning) moves to inner shells with the reduction of α. Also, it is worth noting that there is a configuration (red line) where the system experiments the hydrostatic equilibrium, namely, in such caseR = 0 inside of the stellar compact object.
In Fig. 17, instead, we plot the behavior ofR for the parameter values of α = 1.0, Ξ = −0.7 and for different values of β. We can observe that the system for this configuration can experiment a kind of transition between scenarios with overturning and scenarios with cracking when the β value is increased. Also, it is important to note that there is a configuration (red line) whereR = 0. We can see that if β value increases, the fracture (overturning or cracking) moves to the inner shells of the stellar compact object.
In light of what has been mentioned before, we can conclude that the system can experience internal fracture (cracking or overturning) when it is subjected to small perturbations in its matter sector. If the value α of is increased, these fractures move to the outer regions, whereas β increasing causes the opposite behavior. Also, as shown in Figs. 15, 16 and 17, the stellar compact object can experiment transitions between scenarios with cracking and scenarios with overturning by varying the values of Ξ , α and β, respectively. It is worth also notice that there are configurations with hydrostatic equilibrium even after the perturbation action (see red lines in Figs. 15, 16 and 17). In addition, if one observes again the Fig. 6, namely, the behavior of the local anisotropy of stellar compact object can note that there is a region in x = r/R ∈ (0.4; 0.5) where the anisotropy changes smoothly in its growth with respect to the radial variable r , which in certain measure approximately coincides with maximum values of the perturbed radial forceR in Figs. 15, 16 and 17, demonstrating a link between the presence of local anisotropy and the formation of fracture (cracking or overturning) inside the stellar compact object. Furthermore, we can see from these figures that fractures inside stars tend to occur in the outer regions of the stellar compact object, while anisotropy decreases in these regions, which can be interpreted as showing how this descent of anisotropy in these outer regions can be related to the appearance of these fractures, which is consistent with the work [111], where it was shown that cracking results only if, in the process of perturbation leading to departure from equilibrium, the local anisotropy is perturbed, suggesting that fluctuations of local anisotropy may be crucial in the occurrence of cracking (see the works [117][118][119] where it is shown that anisotropy play an important role in the appearance of cracking).
Given the preceding results, it becomes vital to apply a basic condition founded by Abreu in [101] to determine zones of stability against gravitational cracking given by In this way we plot the profile of v 2 t − v 2 r in the Fig. 18. We can observe that the condition (66) is satisfied by the inner shells of the star, namely, the inner shells are stable regions against gravitational cracking, and the outer regions are unstable regions, which corresponds to Figs. 15, 16, and 17 because the fractures appear precisely in these outer regions. Also, is worth mentioning that the transition between the stable and unstable regions in Fig. 18 is around at x = r/R ∈ (0.4; 0.5), which corresponds to those regions where the local anisotropy undergoes a change in its growth with respect to the radial variable r and with the maximum values of the perturbed radial forceR. Namely, the outer regions are unstable when we add small perturbations in the matter sector of our model. Thus, this result reinforces the idea of the existence of direct relation between the presence of local anisotropy and the apparition of cracking inside the stellar compact object.
In order to avoid leaving only qualitative comments on the important role that local anisotropy plays in the appearance of fractures in the interior of the compact stellar object, we perform the same analysis but with the caveat that now we will consider the fluctuations over the entire material sector of the solution except for radial pressure and local anisotropy; that is, we consider the previous analysis but with the consideration that the radial pressure and the local anisotropy remain undisturbed. Thus, from these considerations, we have obtained the following Figs. 19, 20, and 21. Note that we have labelled the new radial forceR in order to differentiate it fromR. The first that we note from these figures is that the totalR has the opposite sign with respect toR for the same configuration of parameters (see Figs. 15,16,and 17), but also and most importantly thatR has not any fracture (cracking or overturning). Furthermore, we have checked that if the anisotropy function Δ is not considered inR a similar result is obtained also. Also note that the configurations with null total radial forceR remain unchanged (red lines) inR . So in this sense, we show that effectively the presence of local anisotropy or its fluctuations causes the presence of cracking inside the stellar object for our model, which is in accordance with previous works. Now, since we have shown the model is unstable against gravitational collapse and stable about the Harrison-Zeldovich-Novikov stability criterion, it would seem interesting to observe what happens to this total radial force if we perturb it solely on its radial pressure. Therefore, we use the same analysis of cracking, but this time we only perturb the radial pressure, while the other physical quantities remain unchanged. In this case, we obtain Figs. 22, 23, and 24. In this case we labelled the resultant total radial force likeR for distinguish from the previous ones. We notice that the total radial force shoots towards very high values as it approaches the center of the star, which can indicate that if we perturb the total radial force just in the radial pressure a huge force in the direction of the stellar center is induced. At this point, we can associate this with the fact that this model is unstable against gravitational collapse since these small disturbances raise big distortions in the radial direction as causes of gravitational collapse. But, on the other side, the Harrison-Zeldovich-Novikov stability criterion is also satisfied; therefore, we can not affirm with all reliability that these great forces actually cause the gravitational collapse. Even one possibility can be that the huge forces generated by radial disturbances are big but finite, so in this sense, the stellar compact object can support these forces without breaking its static stability and does not collapse under radial disturbances. But since cracking only shows us the trend just when the hydrostatic equilibrium is broken and not beyond its evolution for sufficiently long times to know the final state of the star, we can not affirm whether these disturbances generate the collapse or not.
Also, it is interesting to note that the system also experiences fractures in regions further from the center of the star center and that the configurations with hydrostatic equi-librium (red lines) remain stable despite the perturbations. Despite the great instability of this model, it is possible to find configurations where it is totally stable even under the action of disturbances in its matter sector, which is reinforced from the fact that the static stability criterion is fulfilled.

Conclusions
A new stellar model in the regime of anisotropic fluid has been developed. For such construction, we introduce a convenient ansatz for the radial metric component of the spacetime of the system, and the temporal metric component was obtained using the well-established definition of vanishing complexity for self-gravitating spheres. After getting the static and spherically distributed space-time metric, we used it to obtain the matter sector. The physical behavior was analyzed, finding that the model satisfies the physical requirements to represent realistic stellar compact objects, in this sense, the whole physical conditions of Sect. 5 are fulfilled for a certain set of compactness factors of 0.4023 ≤ u = M R ≤ 0.4104.
On the other hand, the model is stable against convection in its inner regions but unstable in the outer regions of the stellar compact object. In this sense, it could be conjectured about the existence of a central core inside the star that does not allow the outermost layers of the star to collapse inside due to convective movement, which must be common due to the extreme conditions in temperature and densities inside a star. However, the model presents instability against gravitational collapse for the set of compactness parameters mentioned before. In such a sense, we can mentioned that the possible tendency to collapse in this case is not associated with the convective motion inside the stellar compact object. At the same time, we checked that the model satisfies the Harrison-Zeldovich-Novikov, which in a certain sense is opposed to instability with respect to gravitational collapse, so that it could be thought at first that the stars represented by this model have great stability with respect to their hydrostatic balance in the face of radial perturbations, but nevertheless, if they are very intense, there could be a possible collapse. In any case, we cannot affirm this with absolute certainty.
In this regard, the response of the model against small perturbations in its matter sector was analyzed in Sect. 6.8, and the gravitational cracking criterion proposed by L. Herrera was used. Namely, whenever the total non-vanishing radial force appearing after the perturbation goes out of hydrostatic equilibrium, having the possibilities that it is directed inward in the inner part of the sphere and change its direction at some r less than the radius of the star, and the opposite case named as overturning. We analyze specifically, the state of the system just after it departures from equilibrium. Firstly, we have evidence that the system experiments a type of tran-sition between scenarios where the cracking is presented and scenarios where the overturning is presented with the variation of the parameters involved in the model, specifically the fracture point (cracking or overturning) appears in the inner regions of the star if the values of Ξ and β are increased, while the fracture appears in the outer regions with the increase of values of α, namely, the inner shells of the stellar compact object are stable against cracking and the outer ones are unstable.
It is also worth noting that the system can experiment with scenarios in which the total radial perturbed forceR is equal to zero, i.e. scenarios in which the system is in equilibrium despite perturbations in the matter sector. As a result, we can argue that the presence of internal fractures is not required in all cases for this model, as there are scenarios or configurations with absolute hydrostatic equilibrium despite the action of the model's small perturbations in the matter sector. Moreover, we have shown explicitly the direct relation in the appearance of any fracture inside the stellar compact object with the presence of the local anisotropy and its fluctuations, namely, without these two characteristic our model is free of any internal fracture (cracking or overturning even when it is perturbed).
Thus, if the system only experiences fluctuations in its radial pressure, it causes a triggered increase in its total radial force directed to the center of the star, which we believe are scenarios associated with gravitational collapse in principle, but which could also be scenarios of high central density but static stability because the model responds to the Harrison-Zeldovich-Novikov stability criterion. In this sense, we are not sure that the appearance of these large forces can cause the collapse; in fact, there could be the possibility that these huge forces are compensated by the presence of local anisotropy if we also take into account the presence of disturbances in the entire material sector, and that the collapse can occur in extreme cases associated with a large adiabatic compression of the baryonic gas inside the star by huge gravitational forces. That is, we can not make a conclusive statement about whether these forces actually cause a gravitational collapse.
Also, it is worth mentioning the presence of fractures in the external parts of the compact object at the same time, which is to say that in addition to the apparition of huge forces directed to the center of the stellar compact object, there would be a dynamic in the external layers. These results mentioned before show us theoretically in a certain measure that there are a wide number of plausible physical processes giving rise to deviations from local isotropy, which take on relevance since the occurrence of these types of fracture has direct implications on the structure and evolution of the compact object at time scales that are smaller than or equal to the hydrostatic time scale and are associated with a wide number of physical processes that are present in a star, such as intense magnetic fields observed in compact objects, neu-trino viscosity, exotic phase transitions, and the superposition of perfect fluids, among others. In any case, the appearance of these fractures can be speculated about and associated with the origin of several important stellar processes such as star quakes and gravitational collapse, the emergence of inner core fractures will undoubtedly change the conditions for outer mantle ejection in supernova events, as well as the occurrence of glitches and bursts of X-rays and gamma rays in neutron stars.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical work. No experimental data was used.] 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.