The physical acceptability conditions and the strategies to obtain anisotropic compact objects

We studied five methods to include anisotropy, or unequal stress distributions, in general relativistic matter configurations. We used nine acceptability conditions that the metric and physical variables must meet to determine if our models were astrophysically viable. Our analysis found the most effective way to introduce anisotropy while keeping a simple density profile. We also found a practical “rule of thumb” that relates the density at the boundary to the density at the centre of relativistic matter distributions. Additionally, we calculated the configuration radius and encountered that values observed by NICER for PSR J0740+6620 are consistent with several acceptable matter configurations, both isotropic and anisotropic.


Introduction
General Relativity is experiencing an extraordinary era where what was once considered a mathematical curiosity, such as black holes, and faint phenomena, like gravitational waves, have transformed into observable astrophysical entities [1,2].Significant efforts explore the properties of physically viable matter configurations that may describe general relativistic compact objects in various states: static, stationary, or undergoing collapse.Any exact solution to the Einstein Equations has certain restrictions, constraining the metric and the energy-momentum tensor to ensure that emerging space-time geometry is astrophysically reasonable.
Since the seminal result of M.S.R. Delgaty and K. Lake [3], several works have expanded the set of acceptability conditions to obtain more meaningful solutions [4,5,6,7,8].These conditions are elaborated under the assumption that there are two distinct components for the pressure, one radial and the other tangential, which yields a richer and more realistic description of the internal structure of a compact object.The consideration of local anisotropy, where the radial and tangential stresses are unequal (P ̸ = P ⊥ ), has gained recognition as a relevant concept in describing general relativistic stars.This idea can be traced back to the pioneering works of J. H. Jeans [9] and G. Lemaître [10] and has continued to be explored in both Newtonian and relativistic frameworks (see [11,12,13,14,15,16,17,18,19,20,21,22,23,24] and references therein).Notably, a recent paper [25] presents intriguing insights into the instability of isotropic pressure distribution in self-gravitating matter systems.
Throughout this work, we shall consider the first five common assumptions to model nonpascalian fluids in general relativistic matter configurations and examine, through extensive modelling, the consequences of the acceptability conditions.We identify the relevant parameters for a particular equation of state, their range and relevance.We integrate the structure equations implementing every anisotropic equation of state with the same density profile ρ(r) for all configurations.We also identify a comparable set of parameter ranges so as to compare all the physical acceptabilities of the different anisotropic modelling strategies.
Within this framework, for a particular common density distribution, we explore answers to the following two questions: • Which type of anisotropy strategy leads to more acceptable matter configurations?
• Are these acceptable models consistent with the Neutron Star Interior Composition Explorer (NICER) observations [2]?
This paper answers the above questions by organizing our subject matter into several sections.The next section describes the notation and framework of General Relativity.In Section 3, we list the acceptability conditions that our models must meet to be considered candidates for compact stellar objects.Section 4 discusses five approaches to include anisotropy in a general relativistic matter configuration.Next, in Section 5, we explore the parameter space while fulfilling several acceptability conditions and answer the above queries.Finally, Section 6 summarizes our closing remarks and conclusions.

The field equations
Let us consider the interior of a dense star described by a spherically symmetric line element written as ds 2 = e 2ν(r) dt 2 − e 2λ(r) dr 2 − r 2 dθ 2 + sin 2 (θ)dϕ 2 , ( with regularity conditions at r = r c = 0, i.e. e 2νc = constant, e −2λc = 1, and ν ′ c = λ ′ c = 0. We shall consider a distribution of matter consisting of a non-Pascalian fluid represented by an energy-momentum tensor: where ρ(r) is energy density, with P (r) and P ⊥ (r) the radial and tangential pressures, respectively.From Einstein's field equations, we obtain the physical variables in terms of the metric functions as where primes ′ denote differentiation with respect to r.Now, assuming the metric function λ(r) is expressed in terms of the Misner "mass" [35] as The Tolman-Oppenheimer-Volkoff equation (i.e.T µ r ;µ = 0, the hydrostatic equilibrium equation) for this anisotropic fluid can be written as Thus, we can identify two forces competing in compensating the pressure gradient: the "gravitational force", F g and the "anisotropic force", F a .Equation (7) together with constitute the relativistic stellar structure equations.
From equation (7), notice that the pressure gradient becomes less steep when the anisotropy is positive ∆ + = P ⊥ − P > 0, and conversely, it changes more rapidly when the anisotropy is negative ∆ − = P ⊥ − P < 0. The only possibility for negative anisotropy is that the tangential and radial pressures vanish at r = r b .
Thus, for a fixed central stiffness, σ = P c /ρ c , the compactness, C ⋆ , of the sphere increases when there is positive anisotropy ∆ + , and decreases when there is negative anisotropy ∆ − .Concerning positive anisotropy, we can adjust more massive configurations compared to isotropic ∆ 0 = 0 scenarios.If both forces balance, i.e., F g = F a , we obtain a specific matter configuration characterised by vanishing radial pressures and solely sustained by tangential stresses [36].This is because the tangential stresses support the mass shells, reducing the required radial pressure in such circumstances [37].

The physical acceptability conditions
The emerging physical variables have to comply with the various acceptability conditions [3,4,5,7,8], which are crucial when considering self-gravitating stellar models.Only acceptable selfgravitating objects are of astrophysical interest and, in this work, those models have to comply with nine requirements expressed as [8]: [38,39]: (a) That the metric potentials e λ and e ν are positive, finite and free from singularities within the matter distribution, satisfying e λc = 1 and e νc = const at the centre of the configuration.
(b) The inner metric functions match the exterior Schwarzschild solution at the boundary surface.
(c) The interior redshift should decrease with increasing of r.
C2: Positive density and pressures, finite at the centre of the configuration with P c = P ⊥c [39].

C9:
The adiabatic convective stability condition ρ ′′ ≤ 0, which is more restrictive than the outward decreasing density and pressure profiles [6].
Acceptability conditions for general relativistic spheres refer to the criteria that must be satisfied by the metric and physical variables in a relativistic matter distribution to be considered astrophysically viable and consistent within the framework of General Relativity.They are motivated by • Regularity conditions on the physical and metric variables, i.e.C1 and C2: A physically acceptable solution should exhibit regular behaviour, particularly at the centre of the sphere, avoiding singularities or divergences in physical quantities such as energy density, pressure, and metric components.
• Energy conditions and equation of state, i.e.C2, C3, C4 and C5: Relativistic matter distributions are typically required to satisfy certain energy conditions, which impose constraints on the stress-energy tensor components.These conditions ensure the energy density and pressures associated with the matter distribution are within physically reasonable bounds.
• Stability, i.e.C6, C7, C8 and C9: This involves assessing the stability of the matter distribution against perturbations or dynamic changes, ensuring that it remains in a state of equilibrium and does not collapse, cracks or other undesirable behaviours.

Anisotropy heuristic strategies
This section will introduce several assumptions and heuristic strategies to model anisotropy in relativistic matter configurations.Local anisotropy in compact objects is a hypothesis that has gained relevance over time.Nowadays, it is well understood that unequal radial and tangential stresses may increase the stability of neutron star models.However, a complete description of the complex interactions in the fluid that cause such phenomena is still unknown.The most common approaches in introducing anisotropy for modelling relativistic matter configuration are: • Anisotropy proportional to gravitational force.M. Cosenza et al. [13] inspired by the work of Bowers and Liang [12] proposed suitable models for anisotropic matter by considering the anisotropic force proportional to the gravitational one.This relationship leads to the following expression for the difference between the tangential and radial pressures: • Quasi-local anisotropy Local anisotropy can also be considered as the influence of quasilocal variables, which are quantities that are not solely dependent on the state of the fluid at a specific point in space-time [51].These variables, such as the curvature radius r or the compactness µ (= 2m/r), are employed as a quasi-local equation of state to describe anisotropy [29].Within this approach, a particular type of anisotropy is: • Anisotropy proportional to a pressure gradient.Another potential form for the anisotropic force, considering equation ( 7), is for it to be proportional to the pressure gradient.Raposo and collaborators [30] proposed an anisotropy proportional to the covariant derivative of pressure as: where f (ρ) (see appendix B for details) is a generic function of the energy density and k µ = (0, k 1 , 0, 0) is a unitary space-like vector orthogonal to the fluid four-velocity.
• Complexity factor anisotropy.This factor is a quantity defined by decomposing the Riemann tensor, which measures the level of complexity in self-gravitating systems [31,52,53].It reflexes the impact of local anisotropy and density inhomogeneity on the active gravitational mass.Consequently, systems with minimal complexity are represented by homogeneous and isotropic fluids.In the case of anisotropic fluids, satisfying the condition of a vanishing complexity factor with minimal complexity, the local anisotropy can be expressed as follows: • Karmarkar anisotropy.The Karmarkar condition [32] is a relationship among components of the Riemann tensor, given by This condition provides a geometric mechanism for incorporating anisotropy into matter configurations.To express equation ( 13) in a scalar form, we introduce a set of scalar functions known as structure scalars, obtained from the orthogonal splitting of the Riemann tensor (refer to [54] and [33] for more a detailed discussion).Hence, the scalar Karmarkar condition for spherically symmetric static configurations is with where Thus, the induced anisotropy by the Karmarkar condition, written in terms of the physical variables, is given by In the above equations ( 9), ( 10), ( 11), (12), and (17), we denoted the corresponding anisotropic parameter by C GF , C QL , C P G , C CF and C KC , respectively.The relationship between the complexity and Karmarkar anisotropies is evident when we rewrite equation (17) in terms of ∆ CF , i.e.
Where we have set C = C CF = C KC .It is worth mentioning that when ∆ CF = 0 ⇔ ∆ KC = 0 and the only matter configuration for both anisotropic strategies corresponds to the Schwarzchild homogeneous isotropic solution.
Another strategy for implementing anisotropic fluids in General Relativistic matter configurations is the gravitational decoupling approach [34].This procedure assumes that the energymomentum tensor splits into two parts as where T ν µ corresponds to the perfect fluid contribution and θ ν µ describes any other coupled form of gravitational source.Implementing the anisotropic parameter C for modelling this method is unattainable.Thus, comparing the models emerging from this strategy with those executed with all previous techniques is impossible.It deserves a more detailed consideration which will be developed elsewhere.

Anisotropy and physical acceptability
In this section, we discuss the physical acceptability of relativistic anisotropic models.We numerically integrate the structure equations ( 7) and ( 8) implementing every equation of state for anisotropy from the previous section, (i.e.∆ GF , ∆ QL , ∆ P G , ∆ CF , and ∆ KC ) and selecting a common density profile where the central density ρ c and the constant α are free parameters.This simple Tolman VII density profile [55] is not deprived of physical interest [56] and has a long tradition of modelling compact objects.It corresponds to the Gokhroo-Mehra [57] solution used in several anisotropic static spheres in General Relativity [51,58,37].Additionally, under some circumstances [59], it leads to densities and pressures that give rise to an equation of state similar to the Bethe-Börner-Sato Newtonian equation for nuclear matter [60].It also describes radiating anisotropic fluid spheres [59,61] representing the Kelvin-Helmholtz phase in the birth of a neutron star [62,63].Now, from equation ( 20), we obtain the boundary radius of the configuration as a function of the physical parameters of the problem, i.e. ρ c and κ = ρ b /ρ c as Where we have defined the following quantities where R and M are the structure boundary radius and total mass, respectively.In appendix B, we present the dimensionless expressions for the structure equations and each anisotropic EoS, including some information about the simple numerical integration techniques.

The range of the parameters
The solution of equation (7) for the density profile ( 20) is sensitive to ρ c , κ = ρ b /ρ c and the anisotropy factor, C [56].A variation of these three factors generates a parameter space, exhibiting several acceptability conditions satisfied by each model.We shall identify which anisotropy delivers more physically acceptable configurations, i.e. satisfy more acceptability conditions.Thus, we shall identify a common set of parameter variations so as to compare the physical acceptability of the different anisotropic modelling strategies.We start determining the possible variation of a common anisotropic parameter, C. Regarding the case of ∆ GF , observe that equation ( 7) leads to with h = 1 − 2C, and when h = 1 the isotropic case is recovered.Notice that condition C3 and equation ( 23) implies h > 0, therefore if ρ b ̸ = 0 we have, , and since The tangential pressure should be positive at the boundary P ⊥ b ≥ 0 within the matter distribution and from equation ( 9), it restricts the anisotropic parameter to 0 ≤ C < 1 2 for any EoS having ρ b ̸ = 0. We selected six values for the anisotropy parameter, i.e.C = 0.000, 0.050, 0.150, 0.250, 0.350, and 0.450.
In addition to the anisotropic parameter, C, there are two other significant elements: the central density, ρ c , and κ.According to typical values in compact objects/neutron stars, ρ c could go from 0.1 × 10 15 to 2.5 × 10 15 g/cm 3 .The scale variation for κ runs from 0.0 (vanishing density at the surface) to 0.9 (almost homogeneous density).Finally, we have to provide the total mass, M , of the configuration (≈ 2.08 M ⊙ , the highest reliable gravitational mass of any neutron star [64,65]) to determine the central density of each model.

The best method to introduce anisotropy
To answer the first question, we shall follow two lines of reasoning in identifying which of the above anisotropy strategies is best suited in providing more acceptable models.The next section identifies regions in the parameter space (C, κ and ρ c ), that comply with the acceptability criteria.Figure 1 displays, in a colour scale, those patches for five different values of the anisotropy factor C. For example, in the isotropic case, i.e.C = 0, we obtain 33 of these physically fully acceptable models.More acceptable matter configurations are placed below the red line in all cases shown, i.e. when As displayed in Figure 2, the second approach is to sum up the total number of models satisfying all nine-acceptability criteria.This method, discussed in section 5.2.2, complements the previous criterion because we explore the number of possible acceptable models for the whole range of variation of the anisotropic parameter.

Acceptable model distribution in a parameter space
In this section, we shall discuss the acceptable model distribution in the common parameter space defined by 0.000 ≤ C ≤ 0.450; 0.1 × 10 15 ≤ ρ c ≤ 2.5 × 10 15 g/cm 3 and 0.0 ≤ κ ≤ 0.9. Figure 1 displays, in a colour scale, this model distribution for six different values of the anisotropy factor C.
As will be clear in the following discussion, this range of variation in the parameter plane (ρ c , κ) is due to the NICER-acceptable models [64] when we considered the total mass of the configuration M ≈ 2.08M ⊙ .See Figure 3 to grasp the rationale of the parameter variation.For a low anisotropic presence (C = 0.050 displayed in Figure 1), ∆ GF and ∆ QL strategies deliver more acceptable models (yellow patches represent models satisfying the nine criteria) than in the isotropic case.Several models with anisotropy proportional to the pressure gradient, ∆ P G , become unacceptable because they do not meet the adiabatic index's stability criterium C6.On the other hand, configurations with anisotropy defined by the complexity factor have nonphysical negative pressure and positive tangential pressure gradient.Moreover, configurations with vanishing density at the boundary (ρ b = 0 ⇌ κ = 0.0) do not comply with C5, C6 or C8.The Karmarkar anisotropy scheme, ∆ KC , produces unsuitable configurations having negative tangential pressures.When κ = 0.0 the corresponding matter distributions do not comply with C4, Figure 1: Number of conditions fulfilled by the different anisotropy heuristic strategies.The strategies ∆ GF and ∆ QL lead to better modelling (more yellow squares representing models complying with nine criteria) than the isotropic case for C = 0.050, 0.150 and 0.250, This is not the case for the other anisotropic approaches: ∆ CF , ∆ KC , which have the maximum number of acceptable models in the isotropic condition.When the anisotropic parameter increases to C = 0.150, the number of acceptable models boost for ∆ GF and very little for ∆ QL and ∆ KC .For the other approaches, the number of models decreases.Increasing the anisotropy to C = 0.250 generates a maximum of acceptable models with the ∆ GF method.When C = 0.350, no ∆ GF -model satisfies the nine acceptability criteria.The ∆ QL -strategy gives a few new models, and other anisotropic methods decrease the number of acceptable configurations.Nearly all matter distributions meet the nine criteria for the considered values of κ and ρ c .Finally, when C = 0.450, only the ∆ QL -strategy provides more acceptable models than the isotropic condition.Clearly, for all the anisotropic methods, we found more acceptable models when ρ b ≤ 9  10 ρ c 1 − 2ρc 5 .C6 or C8.In general, as κ increases, the speed of sound exceeds the light speed, failing C4.
The increase in the central density oversteps the condition on the trace of the energy-momentum tensor, and the darker region in the upper right corner is due to the models' cracking (C8).Raising C to 0.150 enhances the acceptability when anisotropy is proportional to the gravitational force, and 65 out of 90 models satisfy all the conditions.On the other hand, quasi-local anisotropy does not show variation in the acceptable configurations, while the number of acceptable models with ∆ P G decreases drastically by failing with condition C6.The increased anisotropy in ∆ CF makes acceptable models, now to have positive tangential pressure gradient and negative radial pressure.One model becomes acceptable for ∆ KC by fulfilling condition C5.
Increasing the anisotropic factor C to 0.250 further causes a maximum in ∆ GF .As seen from Figure 1, almost all models, ∆ GF , meet the nine criteria for the considered values of κ and ρ c .This is also evident from Figure 2. It attains 76 fully acceptable models when C = 0.225.The remaining 12 configurations mainly do not satisfy the causality condition of radial sound speed C4.
Quasi-local anisotropy also increases, to a lesser extent, the number of acceptable models.Only one model for ∆ P G becomes unacceptable due to condition C6.Regarding ∆ CF , a few models are no longer acceptable, breaking the tangential pressure condition C3.Acceptable models with ∆ KC remain unchanged.
As it is clear from Figure 1, when C = 0.350, there is no ∆ GF -model satisfying the nine acceptability criteria; the ∆ QL -strategy allows a few new models, while ∆ P G and ∆ CF decrease in acceptable models.In contrast, ∆ KC remains unchanged.Finally, as displayed in Figure 1, i.e. for C = 0.450, only the ∆ QL -strategy provides more acceptable models than the isotropic condition.Several models under the red line, for ∆ CF , become acceptable by fulfilling condition C4.
Table 1: Values of the anisotropic parameter C and the number of acceptable models N (C) in the range of the considered κ and ρ c .C 0 represents the anisotropic parameter where no other fully acceptable model is found.C max represents the central value among the C values that yield the highest number of acceptable models for a particular anisotropy strategy.N (C max ), is greater than the number of acceptable isotropic models for ∆ GF , ∆ QL and ∆ P G strategies.The area under the curve, N T otal ∝ C 0 0 dC N (C), for each anisotropic strategy.Clearly, for the range in κ and ρ c considered, there are more anisotropic acceptable models than their isotropic counterparts.Concerning the NICER-acceptable models, we found N N ICER (C N ICER−max ) = 10, 9, 8, 8 and 8, for ∆ GF , ∆ QL , ∆ P G , ∆ CF and ∆ KC , respectively and the area under the curve for the anisotropic strategies follow the same pattern.Figure 2  This section extends our analysis by determining the total number of entirely acceptable models.To achieve this, we redefine the range for the anisotropic parameter, starting from isotropy (C = 0.000) and continuing until the last value (C = C 0 ) where no other fully acceptable model exists, i.e. satisfying all nine criteria for acceptability.
In Figure 2 and Table 1, we show the most effective methods for introducing anisotropy.The most useful approach is the quasi-local ∆ QL method [29], followed by the Karmarkar scheme [32,33], and finally the complexity factor approach [31,52,53].These strategies include varying degrees of effectiveness due to the significant range of the anisotropic parameter C: 0.000 ≤ C QL ≤ 3.345, 0.000 ≤ C KC ≤ 4.760, and 0.000 ≤ C CF ≤ 3.570, respectively.
Anisotropy leads to more acceptable models than isotropic ones.The most distinct scheme is the anisotropy proportional to gravitational force, ∆ GF .It has a narrow range of variation for the anisotropic parameter, with 0.000 ≤ C GF ≤ 0.250, and a pronounced peak with 76 models for C GF = 0.225.Around half of the models fall within 0.050 ≤ C GF ≤ 0.250.Models with C GF > 0.250 are considered unacceptable due to their positive tangential pressure gradient, which violates condition C3.Regarding the quasi-local approach, ∆ QL , the number of models satisfying Table 2: Range of κ and ρ c for the PSR J0740+6620 with 2.08 ± 0.07 M ⊙ , considering distinct anisotropy approaches.In this table, we present the range of κ for the possible values of ρ c consistent with the observational NICER data.We found various possible configurations for different ρ c , ranging from almost homogeneous density profiles, i.e. κ ≈ 0.8, to others with vanishing density at the boundary of the configuration where κ ≈ 0.0.Assuming a simple density profile, the observed radius for PSR J0740+6620, can be described by acceptable isotropic matter configurations and several anisotropic approaches.Models with the highest κ correspond to anisotropic configurations with ∆ GF strategy.

Range of κ ρ
all requirements increases with the level of anisotropy, reaching a peak of 53 for C QL = 0.910.This method has a significant range in anisotropy 0.000 ≤ C QL ≤ 3.345.The simplest but least effective method is anisotropy proportional to pressure, ∆ P G .It has a limited anisotropic range of 0.000 ≤ C P G ≤ 0.578, with a maximum of 34 models at very low anisotropy, C P G = 0.025.The following method is complexity anisotropy, ∆ CF , which has a considerable range, 0.000 ≤ C CF ≤ 3.570, but no anisotropic parameter generates more acceptable models than the isotropic case.Finally, the geometric Karmarkar strategy is ∆ KC related to the complexity anisotropy and has the broadest range, i.e., 0.000 ≤ C CF ≤ 4.760.

NICER acceptable models
The Neutron Star Interior Composition Explorer is an X-ray telescope on the International Space Station which studies the X-ray emissions from neutron stars, helping to determine their size, mass, and the properties of their dense interiors.By measuring the mass and radius of multiple neutron stars, NICER refines our understanding of the equation of state, providing valuable constraints on the properties of ultra-dense matter.NICER also detects the pulsation of neutron stars, permitting scientists to explore the dynamics of their atmospheres, unravelling the physical processes occurring in and around them.NICER has been employed to obtain the first precise (and dependable) measurements of a pulsar's size and mass and the first-ever map of hot spots on its surface (see [2] and references therein).
In Figure 3, we include a region covering observational estimates for the radius of the PSR J0740+6620, with mass 2.08 ± 0.07 M ⊙ , which is the highest reliable gravitational mass of any neutron star [64,65].In Table 1, we indicate the number of NICER-compatible models for different anisotropic strategies.
From equation (21), we calculate the configuration radius as a function of ρ c and κ.We find that the central density, 0.7 × 10 15 g/cm 3 ≤ ρ c ≤ 1.3 × 10 15 g/cm 3 , is consistent with observational data.The range for κ associated with various central densities for different anisotropy strategies are displayed in Table 2.
We also find various acceptable configurations with different ρ c , ranging from almost homogeneous density profiles, i.e. κ ≈ 0.8, to others with vanishing density at the configuration boundary where κ ≈ 0.0.Assuming a simple density (20), the observed radius for PSR J0740+6620 can be described by acceptable isotropic matter configurations and several anisotropic approaches.Models with the highest κ correspond only to anisotropic configurations within ∆ GF strategy.Figure 3: NICER acceptable models for PSR J0740+6620.There could be several descriptions (isotropic and anisotropic) for this pulsar.We consider three anisotropic methods, ∆ GF , ∆ QL , ∆ P G and the isotropic case.We present the isotropic scenario because the number of acceptable models is the same for the ∆ CF and ∆ KC .We show that the value of all possible radii ranges within 11.6 km ≤ R ≤ 13.1 km, corresponding to several NICER-acceptable models with the selected parameter space, fulfilling all the physically acceptable conditions.The NICER acceptable region overlaps with more yellow models in ∆ GF and ∆ QL plots than in the other anisotropic strategies.The red arrows indicate the possible values for the central density, 0.7 × 10 15 g/cm 3 ≤ ρ c ≤ 1.3 × 10 15 g/cm 3 and the corresponding values of κ consistent with the assumed density profile (20), which are displayed in Table 2.
Figure 3 provides a visual representation of the results detailed in Table 1.All possible radii 11.6 km ≤ R ≤ 13.1 km harmonise with various NICER-acceptable models.The red arrows denote workable values for the central density, ranging from 0.7 × 10 15 g/cm 3 to 1.3 × 10 15 g/cm 3 , coupled with their corresponding κ values that align with the assumed density profile (20).
In both the ∆ GF and ∆ QL anisotropy strategies, the NICER acceptable region encompasses more yellow models compared to other anisotropic approaches.Specifically, there are ten NICER models for ∆ GF with C GF = 0.225 and nine for ∆ QL with C QL = 0.910.However, if we consider the whole range of possible values for the anisotropy parameter C, the total number of NICER models reaches approximately N N ICER−T otal−QL ≈ 11.775 for the ∆ QL strategy (see Table 1).

Final remarks
This work introduces the most common assumptions in modelling non-pascalian fluids in general relativistic matter configurations.Local anisotropy in compact objects is a hypothesis that has gained relevance over time.So far, however, it is still not well known how unequal radial and tangential stresses may increase the stability of neutron star models.The complete description of the complex interactions in the fluid that cause such phenomena is unknown [15,23].
We explore five different heuristic methods to include anisotropy in general relativistic matter configurations.We found that the most effective approach in introducing anisotropy, with a physically meaningful density profile (20) is the quasi-local ∆ QL method [29]; followed by the Karmarkar scheme [32,33]; and last by the complexity factor approach [31,52,53].Incorporating any of the five types of anisotropy schemes considered in this study results in a significantly greater number of acceptable configurations than their isotropic counterparts within the specified range of critical parameters (C, κ and ρ c ).
Furthermore, as shown in Figure 1 and from equation ( 25), we have established a "rule of thumb" that provides a simple relationship between the density at the boundary, ρ b , and the centre, ρ c , for relativistic matter distributions.This rule can serve as a helpful tool for identifying potentially realistic and acceptable models of compact objects.By leveraging this relationship, researchers can make informed judgments about the physical viability of different matter configurations.
From equation ( 21), we calculate the configuration radius as a function of ρ c and κ.We found that the central density, 0.7 × 10 15 g/cm 3 ≤ ρ c ≤ 1.3 × 10 15 g/cm 3 , is consistent with NICERobservational data.In Table 2, we introduce the corresponding ranges for κ.All the possible radii values, 11.6 km ≤ R ≤ 13.1 km, correspond to several NICER-acceptable models within the selected parameter space, fulfilling all physically acceptable conditions.
Assuming a simple density profile (20), the observed radius for PSR J0740+6620 can be described by acceptable isotropic matter configurations and several anisotropic approaches.Models with the highest κ correspond only to anisotropic configurations with ∆ GF strategy.

Appendices A The structure equations
We determine the physical variables (ρ, m, P, P ⊥ ) and check the acceptability conditions.We compare the physical acceptability among models with the same parameters (ρ c , α, C) having different anisotropy strategies.The more acceptable models an anisotropy generates, the more it may represent observable compact objects.Now, expressing the structure equations ( 7) and ( 8) in term of dimensionless quantities we have leaving only equation ( 26) to be integrated.We have this new set of dimensionless physical variables: where R is the boundary radius of the configuration.
It is convenient to transform the parameter α (= αR 2 ) into a quantity with greater physical meaning.Evaluating the dimensionless density in (28) at the surface of the configuration, where x = 1, we can define κ = 1 − α = ρ b /ρ c (30) as the density ratio at the surface to the density at the centre.

B Dimensionless equations of state for anisotropy
The change of variables proposed in (29) to express equations in dimensionless form has the virtue of preserving the equations without additional constants.That is, we can directly put the tilde mark on the variables (and swap r for x in the case of the radial coordinate) to obtain the dimensionless version.Equations ( 26) and ( 27) are an example of what has just been stated.However, here are some simple calculations that prove it.
• Anisotropy proportional to gravitational force.Carrying out the change of variables (29) on the anisotropy proportional to gravitational force yields Rx − 2R m .
Now, rearranging the constant R to the right-hand side, we have from where we have that • Quasi-local anisotropy.Quasi-local anisotropy is a more straightforward case since compactness is a dimensionless variable.Implementing the change of variables (29) in equation (10) leads us to ∆QL Rx , and therefore ∆QL = 2C QL m P x .
• Anisotropy proportional to a pressure gradient.In this particular case, we first choose the function f (ρ) = ρ as in [30], leaving anisotropy (11) as Therefore, the anisotropy factor C P G has dimensions length cubed.Thus, substituting the change of variables into the last equation gives ∆PG

C Numerical integration
Equation ( 26) was numerically integrated with Python, implementing the RK45 method through the solve ivp function.The solution was started at the surface of the model, with initial values x b = 1 and Pb = 0, and proceeded with an adaptive step towards the centre, with final values x c = 10 −15 and P (x c ) = P c .Since x takes values between 10 −15 and 1 we can identify R as the total radius.
Additionally, the interior metric should continuously match the Schwarzschild exterior solution at the sphere's surface, r = r b = R.This implies that e 2ν b = e −2λ b = 1 − 2C ⋆ = 1 − 2M/R, where M = m b is the total mass and C ⋆ = M/R the compactness of the configuration.From now on, the subscripts b and c indicate the variable's evaluation at the boundary and the centre of the matter distribution.

Figure 2 :
Figure 2: In the left panel, we present the number of acceptable models N (C) fulfilling the nine acceptability conditions as a function of the anisotropic factor C for each anisotropic strategy.The embedded plot displays a zoom of the region near the isotropy, i.e.C ≈ 0. The black point represents the number of nine-condition-acceptable-isotropic models.There are more anisotropicacceptable models than isotropic ones.∆ GF has a limited range of anisotropy variation, 0.000 ≤ C GF ≤ 0.250, with a pronounced peak for 76 models for C GF = 0.225.The number of acceptable quasi-local models, ∆ QL , increases with anisotropy, reaching 53 for C QL = 0.910.The areas under the curves are proportional to the total number of anisotropic configurations described by the three parameters: C, ρ c and κ.The quasi-local anisotropic modelling, ∆ QL , is the most effective anisotropic strategy, followed by the Karmarkar scheme, ∆ KC , and then by the complexity factor approach, ∆ CF .The right panel illustrates the number of acceptable models representing the NICER observational estimation of the radius R for PSR J0740+6620.The area under the curve for the anisotropic strategies follows the same pattern in the NICER-acceptable models.More possible models exist for the ∆ QL approach.
illustrates these results.