Spherically symmetric anisotropic strange stars

In this work, we made an extensive study about the possible presence of anisotropies in strange stars. To accomplish this task, we use three different configurations for the strange matter: the unpaired matter, a two-flavor super-conducting (2SC) strange matter, and a fully three-flavor super-conducting strange matter (CFL). For each configuration, we calculate the relevant quantities for the strange stars, such as the mass-radius relation, the dimensionless tidal parameter, the moment of inertia, and the surface curvature for different degrees of anisotropies. Whenever possible, we compare our results with constraints found in the literature, especially focusing on the existence of very massive pulsars (PSR J0952-0607), as well as very light compact objects (HESS J1731-347).


Introduction
Strange stars are self-bounded compact objects composite by deconfined quarks.The theory of strange stars is based on the so-called Bodmer-Witten conjecture [1,2].It assumes that the ordinary matter we know, composed of protons and neutrons, may only be meta-stable, while the true ground state of strongly interacting matter would therefore consist of the so-called strange matter (SM), which in turn is composed of deconfined up, down and strange quarks.For the SM hypothesis to be true, the energy per baryon of the deconfined phase (for p = 0 and T = 0) is lower than the nonstrange infinite baryonic matter, i.e., (E uds /A < 930 MeV), while at the same time, the nonstrange matter still needs to have an energy per baryon higher than the one of nonstrange infinite baryonic matter (E ud /A > 930 MeV); otherwise, protons and neutrons would decay into u and d quarks.
a e-mail: llopes@cefetmg.brb e-mail: harishchandra.das@ct.infn.it In this work, we discuss strange stars within three different possible scenarios.The normal, unpaired quark matter.Two flavors paired color-superconducting quark matter (2SC) and a full three flavors, color-conducting quark matter known as Color-Flavor-Locked (CFL) phase.The quark matter is expected to behave in a full CFL phase at very high densities.However, at intermediate densities, the behavior of the quark matter is still an open issue [3,4].
Furthermore, we study not only isotropic strange stars but also investigate the effects of local anisotropy.The local anisotropy can arise due to exotic phenomena, such as a strong magnetic field, superfluity, phase transitions, etc. (see Ref. [5] for an extensive discussion).Several models have been proposed to explore the effects of pressure anisotropy inside the star [6].However, in this study, we use the socalled Bower-Liang (BL) model [7].In this simple model, the anisotropy is gravitationally induced and non-linear in pressure.Moreover, the BL model has the advantage that the presence of local anisotropy does not break the spherical symmetry of the star, as discussed in Ref. [8].
We have already discussed anisotropic stars in our previous work [9].However, in that short study, we considered only one model for the quark matter, the unpaired one.Moreover, that study also restricted our analysis to the massradius diagram and the tidal deformability parameter.In the present study, we extend our study in all aspects such as the unpaired matter and two models of color-superconducting phases.Moreover, altogether with the mass-radius relation and the dimensionless tidal parameter, we also discuss the effect of color-superconducting and anisotropy in the inertia moment of strange stars and the surface curvature [10].This work is, therefore, complementary to other studies found in the literature.For instance, the effect of anisotropy on radial oscillation was studied in Ref. [11].The anisotropic charge in strange stars has been explored in the Tolman-Kuchowicz space-time geometry in Ref. [12].They got an interesting result that the pressure anisotropy initially dominates the Coulomb repulsive forces, and the repulsive forces dominate the anisotropic when the radius increases.On the other hand, in Ref. [13], the authors have explored the effect of the magnetic field on the anisotropic strange stars and observed that the transverse component of the magnetic field increases the mass and size of the anisotropic star and vice-versa for the radial component of the field.
We also check which of our results fulfill some constraints, mainly focusing on those that are harder to explain for a traditional neutron star.The first is the HESS J1731-347 object [14] whose mass and radius are M = 0.77 +0.20  −0.17 M ⊙ and R = 10.4 +0.86 −0.78 km respectively.The true nature of the HESS object is one of the hot topics in nuclear astrophysics.In Ref. [14], it was shown that some chiral theory can explain this object as an ordinary hadronic neutron star.The hadronic nature was also explored in Ref. [15].On the other hand, the HESS object as a strange star was studied in refs.[9,16].Finally, in Ref. [17] the authors show that the HESS object can be a hadronic neutron star with a soft EOS, or a hybrid star with an early deconfinement phase transition.The second constraint is the very small upper limit of the radius (R 1.4 < 11.9 km) of the canonical star, given by Capano et al. [18].Here again, we cannot completely rule out hadronic neutron stars, although such constraint is tight enough to rule out mostly hadronic EoSs.Finally, an extreme constraint is the speculative mass of the black widow pulsars PSR J0952-0607, M = 2.35 ± 0.17M ⊙ [19].Therefore, in this study, we provide an analysis of the anisotropic stars based on different models that could be able to fulfill all three strong constraints simultaneously.
Despite those extreme cases, we also investigate if our models are able to describe some more traditional constraints.For instance, the PSR J0740+6620 with a mass of 2.14 M ⊙ was pointed in Ref. [20].Neutron Star Interior Composition Explorer (NICER) results refined this data, and today we have the PSR J0740+6620 with a mass of 2.08 ± 0.07 M ⊙ and a radius of 12.35 ± 0.35 km [21].The canonical star, M = 1.4M ⊙ , also received great attention in the last years.Two NICER results point that the radius of the canonical stars must be in the range 11.52 < R 1.4 < 13.85 km [22], and 11.96 < R 1.4 < 14.26 km [23].Nowadays, these results were revised to 11.80 < R 1.4 < 13.10 km [21].We also investigate if our results are able to explain the mass-gap object in the GW190814 event [24], whose mass was estimated to be 2.50 − 2.67 M ⊙ .Still, there are several debates about its true nature; see Refs.[25,26].
In addition, we also take the constraints from the gravitational wave observations by LIGO/VIRGO/KAGRA.The GW170817 event [27] put the constraints on the dimensionless tidal parameter of the canonical star Λ 1.4 < 800 [27].This result was then refined in Ref. [28], to 70 < Λ 1.4 .< 580.Moreover, assuming that the mass-gap object in the GW190814 event was not a black hole implies that the dimensionless tidal parameter for the canonical star must be in the range of 458 < Λ 1.4 < 889 [24].
Another important quantity is the moment of inertia (MOI) of the compact stars.Till now, we don't have any observational data for any NS.The authors of Ref. [29] have obtained the MOI for several pulsars using the universal relations between the mass and the tidal deformability.Here, we calculate the MOI of the anisotropic SQS with different SQM models by varying the degrees of anisotropicity.One can constrain the value of MOI from the future observational data for the anisotropic compact star.There are other ways to constrain the magnitude of the MOI for different systems, such as Millisecond Pulsars (MSP), Double NS (DNS), and Low Mass X-ray Binary (LMXB), as done in Refs.[29,30].The MOI of these pulsars is expected to be measured soon; for example, the PSR J0737-3039(A) is the only known DNS up to date.Since the NS equation of state (EOS) is believed to be universal, the tidal deformability constraints from GW170817 have implications for all NSs, including PSR J0737-3039(A), the MOI has been obtained to be 1.15 +0.38  −0.24 × 10 45 g cm 2 [29].One can also estimate the MOI for other MSPs and LMXBs.
Finally, we discuss the surface curvature (SC) of the strange stars.The SC is a quantity used to measure the magnitude of the curvature made by a compact object due to its huge mass, and it is also one of the major consequences of Einstein's general theory of relativity.The SC is one of the observables that can be measured for different compact objects.Recently, the direct detection of gravity-field curvature has also been observed by using atom interferometers [31].In the future, we might have modern interferometers like LISA that could be able to answer the SC for compact objects such as neutron stars/white dwarfs.Since the SC mainly depends on macroscopic properties such as mass and radius, we can put constraints on the internal composition of the star.We explain the details of the methodology to obtain the SC in the following sections.
2 Quark matter models

Unpaired quark matter: the vector MIT bag model
We begin our study with a normal, unpaired quark matter.This implies that all the flavors and colors are degenerated.Moreover, the particle distribution in the momentum space forms a sphere, the Fermi sphere, whose radius is the socalled Fermi momentum, k f .At T = 0K, all particles are constrained inside the Fermi sphere.Consequently, the Fermi-Dirac distribution becomes the Heaviside step function [32].To model the unpaired quark matter, we use a vector-enhanced MIT bag model [33].
The vector MIT bag model is an extension of the original MIT bag model [34] that incorporates some features of the quantum hadrodynamics (QHD) [35].In its original form, the MIT bag model considers that each baryon is composed of three non-interacting quarks inside a bag.The bag, in its turn, corresponds to an infinity potential that confines the quarks.As a consequence, the quarks are free inside the bag and are forbidden to reach its exterior.All the information about the strong force relies on the bag pressure value, which mimics the vacuum pressure.
In the vector MIT bag model, the quarks are still confined inside the bag, but now they interact with each other through a vector meson exchange.This vector meson plays a role analog to the ω meson of the QHD [35].The Lagrangian density reads [36,37]: where m q is the mass of the quark q of flavor u, d or s, ψ q is the Dirac quark field, B is the constant vacuum pressure, and Θ ( ψq ψ q ) is the Heaviside step function to assure that the quarks exist only confined to the bag.Imposing meanfield approximation (MFA) and applying Euler-Lagrange to Eq. ( 1), we obtain the energy eigenvalue for the quark, as well as the expected value for the vector field.
where n q is the number density of the quark q and µ q is its chemical potential.Now, applying Fermi-Dirac statistics, the energy density is analogous to the QHD plus the bag term: γ q = 6 (3 colors × 2 spins) is the degeneracy factor.Electrons are also added as a free Fermi gas to ensure chemical equilibrium.The pressure is then easily obtained by thermodynamics relations: p = µn − ρ.
In this work we use m u = m d = 4 MeV, m s = 95, and define G V ≡ (g V /m V ) 2 , as suggested in Refs.[36,37].Although the bag constant B and the parameter G V are not fully independent, they are weakly constrained to only satisfy the stability window [38].This implies that the maximum mass of an isotropic strange star can vary from 1.61 M ⊙ for G V = 0 [36], up to 2.81 M ⊙ for G V = 0.40 fm 2 [25].
However, it is well accepted in the literature that colorsuperconducting matter is stiffer than unpaired one [4,[39][40][41][42].We then use this fact to constrain the values of G V and B by assuming that this model produces a softer EoS than the 2SC color-superconducting quark matter.Therefore, in this work, we take G V = 0.15 fm 2 and B 1/4 = 150 MeV.

2SC color-superconducting quark matter
The 2SC color-superconducting is the simplest pairing that evolves a BCS coupling, although it is also the least symmetrical one.In this model, the s−quark does not pair due to its high mass value.The lightest quarks form a Cooper pair if they present different colors and flavors (anti-symmetric coupling).As a consequence, only two of the three colors form Cooper pairs.Therefore, we have two of the three flavors and two of the three colors paring.This means that of the nine quarks, only four are paired [4,40].The presence of unpaired quarks gives rise to gap-less quasi-particles.The density of states of such quasi-particles is proportional to µ 2 and, therefore, is very large [43].Therefore, the 2SC EOS is not much stiffer than the unpaired one, as already suggested in Ref. [40].
In the 2SC color-superconducting matter, as in the unpaired one, a small quantity of electrons is needed to keep the electric charge net equal to zero.In Ref. [44], the 2SC color-superconducting was explored within the NJL model.Moreover, in Ref. [45], the authors show that the 2SC phase can be accurately described by an analytical approximation called constant-sound-speed (CSS) parametrization.The CSS model also seems to accurately reproduce the Field Correlator Method (FCM) with and without a color-superconducting phase [46,47].
In the CSS model, the pressure is a linear function of the energy density.Due to its simplicity and accuracy (see, for instance, Figs 1 and 2 in Ref. [45]), the CSS model is widely used in modern literature [45][46][47][48][49][50][51][52].The CSS EoS and total number density read We have, therefore, three free parameters: the square of the speed of sound (v 2 s = a), the energy density at p = 0 (ρ * ), which plays a role similar to the bag in the MIT bag models, and the number density at p = 0 (n * ), which in turn, plays the role of the saturation density (n 0 ) of the MIT based models.In Ref. [45], the authors freely vary the value of a in the range 0.2 < a < 0.8 and found that -depending on the NJL parametrization -the 2SC phase is well described with a < 0.33 while the CFL phase is described by a > 0.35.Therefore, following that paper, we use a = 0.302, which was explicitly used in their fig.3 to describe the 2SC quark matter.

CFL color-superconducting quark matter
If in the 2SC model, only four of the nine quarks are paired, in the CFL, all of them are paired in a fully anti-symmetric coupling in the space of colors and flavors [4].Also, there are no gap-less quark quasi-particles in the low-energy spectrum of the CFL phase due to the color Meissner effect [43].Moreover, the Cooper pairing in the CFL phase helps to enforce the equal number densities of all three quark flavors, n u = n d = n s , implying that electrons are not present in this phase.In general, the CFL phase is expected at densities above four times the nuclear saturation density, n 0 , while the 2SC phase is expected in the region 2 < n/n 0 < 4 [45].However, some studies indicate that the 2SC phase is not favored at all due to the significant free energy cost of the 2SC phase [3].In this case, the unpaired quark matter undergoes a phase transition to a full CFL phase, or even a more exotic phase, as a non-BCS pattern.The gap parameter (∆ ) that determines the pairing strength of Cooper pairs can be very high in the CFL phase, above 250 MeV [39].In this case, a very stiff EoS and consequently a very massive strange star is produced.
As we did in the 2SC phase, we also use the CSS model (Eq.4) as an analytical approximation to the CFL phase.Here, following the Ref. [45], we use a = 0.57, as it was explicitly used in their fig.3 to the CFL phase.We also use ρ * = 209 MeV/fm 3 , and n * = 0.235 fm −3 , which produce a energy per baryon E/A = 897 MeV, satisfying the Bodmer-Witten conjecture.

Hydrostatic equilibrium
In this section, we briefly review the theory of static equilibrium distribution of matter, which is spherically symmetric but whose stress tensor is, in general, locally anisotropic, as originally introduced in Ref. [7].In the isotropic case, the stress-energy tensor reads: T µν = diag(ρ, −p, −p, −p).We now introduce anisotropies without breaking the spherical symmetry by assuming the following stress-energy tensor: T µν = diag(ρ, −p r , −p t , −p t ).Spherical symmetry implies that (in canonical coordinates) the stress-energy tensor T µν is diagonal, and moreover p φ = p θ = p t [8].
We now can redefine the anisotropic energy-momentum tensor as done in Refs.[53][54][55]:  where ρ, p r , and p t are the energy density, radial pressure, and tangential pressure, respectively.k µ is the unit radial vector (k µ k µ = 1) with u µ k µ = 0.The Schwarzschild metric for this type of star having the spherically symmetric and static configuration is defined as where r, θ , and φ are the Schwarzschild coordinates.Applying Einstein's field equations, we obtain: and the contracted Bianchi identities give us: Finally, by isolating the ν ′ in Eq. 8 and replacing it in Eq. 9, we can write the equilibrium equations in the Tolman-Oppenheimer-Volkoff form [53]: where σ = p t − p r is the anisotropy parameter, 'm' is the mass enclosed within the radius r.The radial pressure is then obtained from a pre-determined EOS.On the other hand, for the case of the transverse pressure, we use the BL model in the following [7]: where the factor λ BL measures the degree of anisotropy in the fluid.There are some boundary conditions required to solve the above Eqs.(10)(11)(12) as done in Refs.[6,56].Also, different fluid conditions must be satisfied for the anisotropic stars such as (i) p r , p t , and ρ > 0, (ii) 0 < c 2 s,t < 1, (iii) p r = p t for r = 0, etc.Other conditions are mentioned in Ref. [6].We plot the mass-radius relation for different values of λ BL for unpaired, 2SC, and CFL strange stars in Figure 1.
We begin our analysis by comparing the different models of quark matter in the symmetric case.It is observed that the unpaired matter produces a maximum mass of 2.09 M ⊙ , the 2SC produces a maximum mass of 2.16 M ⊙ , and the CFL produces an impressive maximum mass of 2.76 M ⊙ .
In the case of the three main constraints, we see that the unpaired matter is not able to reproduce the radius of the HESS object [14] in the absence of anisotropies, while the 2SC and CFL easily fulfill this task.In relation to the constraint related to the radius of the canonical star 10.4 < R 1.4 < 11.9 km, presented in Ref. [18], all quark models are able to fulfill it in the isotropic case.Our last main constraint is the speculative mass of the black widow pulsars PSR J0952-0607, M = 2.35 ± 0.17M ⊙ [19].In this case, only the CFL quark matter is able to reach such high mass, although the 2SC is very close.
In relation to the other constraints, we see that no isotropic model is able to fulfill the revised NICER results for the canonical star, 11.80 < R 1.4 < 13.10 km [21] (Notice that the overlap between the results in Ref. [18] and those in Ref. [21] is only 100 m), but the CFL still fulfill its original range [22].In relation to the PSR J0740+6620 [21], we see that all quark models reach the mass range of 2.08 ± 0.07 M ⊙ , but the unpaired matter produces too-low radii, in disagreement with the range of 12.35 ± 0.35 km.Finally, the mass-gap object in the GW190814 event, with M > 2.50 M ⊙ [24], can be explained only by the CFL phase.Before we finish the analysis of the isotropic case, it is worth pointing out that most of the constraint can be fulfilled with unpaired matter by simply increasing the value of G V [25].However, in this work, we impose that unpaired matter must be softer than the super-conducting one.
We now analyze the effects of anisotropies.As expected, a positive value of λ BL increases the maximum mass, as well as increases the radius of a star with fixed mass [7,9], while negative values do the opposite.For the extreme positive value (λ BL = +2.0),we see that the maximum mass increases by around 23% in the MIT model, 25% in the 2SC, and 22% in the CFL phase.In the extreme negative value (λ BL = −2.0),we see that the maximum mass decreases by around 16% in all three quark models.
In relation to the radius of the canonical star, the extreme positive value (λ BL = +2.0)causes an increase by around 2.5% in the MIT model, 2.3% in the 2SC, and 1.0% in the CFL phase.In the extreme negative value (λ BL = −2.0),we see a decrease of around 3.2% in the unpaired matter, 3.0% in the 2SC, and 1.8% in the CFL phase.
We can conclude that the effect of the anisotropies affects more the maximum mass than the radius of the canonical star.At the same time, its effects in the maximum mass are almost independent of the quark matter EOS, but the effects on the canonical radius are strong for a softer EOS.
In relation to the constraints, we see that a positive value of λ BL will improve the results of the MIT.For instance, for λ BL > 0.5, the HESS object can be described even in the unpaired quark matter.For λ BL > 1.0, the radius range of the PSR J0740+6620 [21] is satisfied, as well as the mass of the black widow pulsars PSR J0952-0607.Additionally, for λ BL = 2.0, even the mass-gap object in the GW190814 event can be explained as an unpaired quark matter.For the 2SC, we see that λ BL > 0.5 is enough to fulfill the mass range of the black widow pulsars PSR J0952-0607, while for λ BL = 1.5, the mass of the secondary object in the GW190814 event is reached.Finally, in the case of the CFL phase, we have a different behavior.As the CFL matter presents a stiff EOS even in the isotropic case, if we increase λ BL up to values of 1.5, their curves do not cross the radius range of the PSR J0740+6620 pulsar [21].On the other hand, most of the constraints are satisfied even for extremely negative values of λ BL .

Tidal deformation
The star shape is deformed due to its presence in the external field of its companion star.The degree of deformation is measured by the parameter λ , which is defined as the tidal deformability of a star.The definition for the dimensionless tidal deformability Λ , which is a quantity observed by the LIGO/Virgo, and it has a unique relation with the tidal Love number (k 2 ), and the compactness (C) of the star is [57,58]: We use the linear perturbation in the Throne and Campolattaro metric to determine the value of k 2 , [59].In the case of anisotropic stars, we use the formalism originally introduced in Ref. [56]: The term d p t /dρ represents the change of p t with respect to energy density for a fixed value of λ BL .The internal and external solutions to the perturbed variable H at the star's surface can be matched to get the tidal Love number [57,60].The value of the tidal Love number can then be calculated using the equation in the following [57,58,61,62].
where y 2 depends on the surface value of H and its derivative where the ρ s is the energy difference between the internal and external regions.The results are presented in Figure 2.
As in the last sub-section, we begin our analyses with the isotropic case.In the case of unpaired matter, it provides the softer EOS, which predicts the lower values of Λ .For the canonical mass, we have Λ 1.4 = 383.As we consider colorsuperconducting, the EOS becomes stiffer, and the dimensionless tidal parameter increases.The values of Λ 1.4 = 622 and 787 for the 2SC and the CFL phases, respectively.In relation to the effects of the anisotropies on the dimensionless tidal parameter, we see that the higher values of λ BL predict the higher value of Λ 1.4 and vice-versa.For an extreme positive value, λ BL = +2.0,we see an increase of the Λ 1.4 by around 26% in the MIT model, 24% in the 2SC, and 14% in the CFL phase.For an extreme negative value, λ BL = −2.0, it is observed that the decrease of the Λ 1.4 by around 22% in the MIT model, 20% in the 2SC, and 11% in the CFL phase.Therefore, as in the case of the radius of the canonical mass, the effects of anisotropies in the dimensionless tidal parameter are stronger in softer EOS.
We now discuss the results in the light of two constraints: the GW170817 (70 < Λ 1.4 .< 580) [28] and the GW190814 event (458 < Λ 1.4 < 889) [24], although it is worth to emphasize that the nature of the mass-gap object in the GW190814 is not know yet, therefore such constraint is still speculative.We see that the MIT model fulfills the GW170817 constraint for all values of λ BL .Moreover, for λ BL = +2.0, the constraint of GW190814 is also satisfied.For the 2SC, we see that the GW170817 is satisfied for values up to λ BL = −1.0,and the GW190814 for all degrees of anisotropies.Finally, the CFL model fails to fulfill the GW170817 for all values of λ BL , but it can describe the GW190814 limit if λ BL < +2.0.

Moment of Inertia
For a slowly rotating NS, the equilibrium position can be obtained by solving Einstein's equation in the Hartle-Throne metric, as given in Refs.[63][64][65]: − 2ω(r)r 2 sin 2 θ dt dφ For a slowly rotating anisotropic compact star, the MOI can be calculated as in Ref. [66]: where,  [29,30].and ω = ω/Ω , where ω is the frame dragging angular frequency, ω = Ω − ω(r).J is defined as e −ν (1 − 2m/r) 1/2 .Hence Eq. ( 18), can be rewritten using Eq. ( 12) as As already pointed out in Ref. [66], a measurement of NS moment inertia is crucial because its inherent capability to constrain quite restricts the EOS of NSs at high density are insensitive to EOS, and it has a universal relation with compactness and tidal deformability.The results are displayed in Fig. 3.
At first approximation, the MOI is proportional to MR 2 .Therefore, the larger and heavier strange stars will present a significantly higher MOI than smaller and lighter ones.This is reflected in Fig. 3. Considering the isotropic case, we see that the MOI for the canonical star (I 1.4 × 10 45 g.cm −2 ) are 1.47 for the unpaired matter, 1.66, 1.78 for the 2SC and the CFL respectively.In relation to the maximally massive strange stars, the MOI (I max × 10 45 g cm −2 ) are 2.36, 2.82, and 4.92 for MIT, 2SC, and CFL, respectively.We can also see that the MOI is strongly affected by anisotropies.For positive values of λ BL the increase of the I 1.4 can reach 6.8% in the MIT, 6.6% for the 2SC, and 3.4% for the CFL.For negative values of λ BL , we have a decrease of the I 1.4 of 8.1% in the MIT, 7.8% in the 2SC, and 4.5% for the CFL.These results are expected; as pointed out earlier, the effects of the anisotropy in the canonical star are stronger in the softer EoS.On the other hand, in relation to the maximum mass of strange stars, we see for positive values of λ BL , an increase of 65 % in the MIT, 67% in the 2SC, and 62% in the CFL.For negative values of λ BL , we have a decrease of 29% in the MIT, 33% in the 2SC, and 32% in the CFL.These large variations reflect the dependence of the MOI with MR 2 .
In relation to the constraints, due to the high uncertainty in the mass, most of the results satisfy the DNS, MPS, and LMXB bounds discussed in Refs.[29,30].However, for higher values of λ BL in the 2SC case, some curves do not cross some of the DNS and MSP regions.In the CFL case, even the isotropic case does not cross some curves of the MSP bound.Nevertheless, we pay special attention to the constraint of the PSR J0737-3039(A).With a well-established mass of 1.34 M ⊙ , in Ref. [29] was able to constraint its MOI to the range 1.15 +0.38  −0.24 ×10 45 g cm 2 .For the MIT, we see that such constraint is satisfied by all values of λ BL .For the 2SC it is satisfied if λ BL < +1.5, and for the CFL, it is satisfied only if λ BL < −0.5.

Surface Curvature
In the general theory of relativity (GR), there are different quantities to measure the curvature of space-time, mainly the Ricci scalar, the Ricci tensor, the Weyl tensor, and the Riemann tensor.However, except for the Riemann tensor, all of the magnitude becomes zero outside the star.Therefore, we adopt the curvature quantity from the Refs.[10,67], that is known as the Kretschmann scalar (full contraction of the Riemann tensor) and is defined as At the surface, m → M as r → R. Except for the last term, all other terms vanish outside the star because they depend on ρ(r) and P(r), which are zero outside the star.But, there is a non-vanishing component of the Riemann tensor that does not vanish; R 1 010 = − 2M R 3 = −ξ , even in the outside of the star [10,68].Therefore, the Riemann tensor is a more relevant quantity to measure the curvature of the stars.Kretschmann scalar is the square root of the full contraction of the Riemann tensor.The vacuum value for K is as one can easily see from Eq. ( 21).Therefore, one can take K reasonable measures for the curvature within the star.The SC is defined as the ratio of curvature at the surface of the NS K (R) to the curvature of the Sun K ⊙ , SC = K (R)/K ⊙ .At the surface of the Sun, the curvature is equal to K ⊙ = 3.0 × 10 −17 km −2 , which can be considered as a small quantity.We show in Fig. 4 the quantity SC in units of ×10 14 .We first notice that, despite the model utilized, the selected model, or even the anisotropy degree, all values of the SC are in the same order of magnitude: 10 14 times the curvature of the sun.In relation to the isotropic case, we see that the SC for the canonical star is (×10 14 ) 3.93 for the unpaired matter, 3.16 for the 2SC, and 2.94 for the CFL.When anisotropies are taken into account, as a positive value of λ BL increases the radius of a fixed mass, we have a decrease of SC.For λ BL = +2.0,we have a decrease of around 8% for the MIT, 7% for the 2SC, and 4% for the CFL.On the other hand, for λ BL = −2.0,we have an increase of around 10% for the MIT and for the 2SC and 5% for the CFL.For more massive strange stars, we see that anisotropy has a stronger influence.
Although today, there is no data for the SC of neutron stars, as pointed out in Ref. [10], this quantity can potentially act as a strong constraint in the future.For instance, in Ref. [67], it was shown that the presence of dark matter in neutron star interiors can increase the curvature by one order of magnitude.0.91-1.53 The main results of this study, as well as some constraints discussed, are summarized in Tab. 1.Before we finished, we quickly analyzed some results presented in this table.As pointed out before, our main goal is to fulfill three constraints that are not easily fulfilled simultaneously in ordinary hadronic models: the small radius of the canonical star (R 1.4 < 11.9 km) presented in Ref. [18]; the existence of very light compact objects, HESS J1731-347 [14]; and the speculative mass of the black widow pulsar PSR J0952-0607 [19].
The models that satisfy simultaneously these three constraints are the MIT for λ BL ≥ +1.0, 2SC for λ BL ≥ +0.5, and the CFL for all degrees of anisotropy.If altogether with those three bounds, we impose that the model also needs to fulfill the mass and radius range of the PSR J0740+6620 [21] and the inertia moment of the J0737-3039(A) [29], therefore we must rule out the CFL for λ BL ≥ −0.5, while for the 2SC, only the values λ BL = +0.5 and +1.0 are now valid.On the other hand, if altogether with the three main bounds, we impose that the model must satisfy the constraint related to the dimensionless tidal parameter in the GW170817 event.Therefore, we must rule out all models based on super-conducting quark matter.We finish by pointing out that unpaired quark matter with λ BL = +2.0virtually fulfills every constraint presented.The only exception is the radius constraint for the canonical star.However, we must emphasize that the window comprehending the NICER (R 1.4 > 11.8 km) [21] and the Capano (R 1.4 < 11.9 km) [18] is very narrow.

Conclusion
In this work, we study the effect of the anisotropy in strange stars within three different quark models: unpaired quark matter, 2SC, and CFL.For the unpaired quark model, we use the vector MIT bag model as introduced in Ref. [36,37].For both super-conducting quark matter (2SC and CFL) we use the CSS model [45][46][47][48][49][50][51][52].The parameters for the 2SC and CFL are taken from Ref. [45], and as an additional constraint, we impose that the unpaired quark matter must be softer than the color super-conducting one, as indicated in the literature [4,39].The effects of the anisotropy are studied in the BL model [7], which has the advantage of preserving the spherical symmetry of the star [8].
We searched for models that are able to satisfy three strong constraints simultaneously: the small radius of the canonical star (R 1.4 < 11.9 km) presented in Ref. [18]; the existence of very light compact objects, HESS J1731-347 [14]; and the speculative mass of the black widow pulsar PSR J0952-0607 [19].Additional constraints, as well as additional features such as the dimensionless tidal parameter, Λ ,  .
the MOI, and SC are also presented.The main results are summarized below: -In the mass-radius diagram, we observed that the presence of the anisotropy increases the maximum mass, as well as the radius of a fixed mass star for positive values of λ BL , and vice-versa for negative ones.-The effect of anisotropies in the maximum mass is similar for all three quark matter models; however, the effect in the canonical star is stronger in the softer models.-The models that are able to fulfill the three main constraints simultaneously are the MIT for λ BL ≥ +1.0, 2SC for λ BL ≥ +0.5, and the CFL for all degrees of anisotropy.
-The anisotropy increase the tidal dimensionless parameter Λ for positive values of λ BL and reduce it for negative ones.Despite this, only unpaired quark matter can fulfill the constraints related to the GW170817 event.
-The MIT presents the smaller values of MOI, and the CFL model has the large MOI.The anisotropies change the MOI of the canonical star by around 10%.However, due to the proportionality of I with MR 2 , the anisotropy strongly affects the maximum mass of strange stars.The MOI increases more than 60% for positive λ BL values and decreases around 30% for negative values.Only for the MIT the constraints related to the J0737-3039(A), are satisfied for all values of λ BL .
-The anisotropy has a weak influence of the SC in lowmass strange stars and a moderate one in most massive ones.For all values of mass and all degrees of anisotropy, the SC is always in the same magnitude order (10 14 ).-The unpaired matter with λ BL = +2.0virtually fulfill all constraints presented, with the exception of the revised NICER constraint.-Therefore, it has been observed that pressure anisotropy has significant effects on the various strange star properties.There are still unsolved issues based on the origin of anisotropy and its effects on the various properties of compact stars, which will be discussed in future work.

Fig. 1 :
Fig. 1: Mass-radius relation for unpaired (top), 2SC (middle), and CFL (bottom) strange stars for different values of anisotropies.The hatched areas are bounds discussed in the text.

Fig. 2 :
Fig. 2: Dimensionless tidal parameter (Λ ) for unpaired (top), 2SC (middle), and CFL (bottom) strange stars for different values of anisotropies.The bounds related to the GW170817 and GW190814 events are also shown.

Fig. 3 :
Fig. 3: Moment of Inertia (I) for unpaired (top), 2SC (middle), and CFL (bottom) strange stars for different values of anisotropies.The bounds are discussed in the text and in Refs.[29, 30].

Table 1 :
Some strange star properties for unpaired (MIT), 2SC, and CFL quark matter and some constraints discussed in the text.The inertia moment is given in units of × 10 45 g.cm −2 .The scalar surface, SC = K (R)/K ⊙ , is given in unit of10 14