Stable and self-consistent compact star models in teleparallel gravity

In the framework of Teleparallel Gravity, we derive a charged non-vacuum solution for a physically symmetric tetrad field with two unknown functions of radial coordinate. The field equations result in a closed-form adopting particular metric potentials and a suitable anisotropy function combined with the charge. Under these circumstances, it is possible to obtain a set of configurations compatible with observed pulsars. Specifically, boundary conditions for the interior spacetime are applied to the exterior Reissner-Nordstr\"om metric to constrain the radial pressure that has to vanish through the boundary. Starting from these considerations, we are able to fix the model parameters. The pulsar $\textit {PSR J 1614--2230}$, with estimated mass $M= 1.97 \pm 0.04\, M_{\circledcirc},$ and radius $R= 9.69 \pm 0.2$ km is used to test numerically the model. The stability is studied, through the causality conditions and adiabatic index, adopting the Tolman-Oppenheimer-Volkov equation. The mass-radius $(M,R)$ relation is derived. Furthermore, the compatibility of the model with other observed pulsars is also studied. We reasonably conclude that the model can represent realistic compact objects.


I. INTRODUCTION
Soon after the formulation of General Relativity (GR), several theories were constructed in view of fixing as many issues as possible related to the gravitational field. Among these theories, there is the one formulated by H. Weyl which tried to unify gravitation and electromagnetism in 1918 [1]. Einstein himself, in 1928 [2], adopted the same philosophy by Weyl adopting the Weitzenböck geometry. In this formulation, one has to introduce tetrad fields to describe dynamics, unlike GR whose dynamical variable is the metric. The tetrad field has 16 components which made Einstein think that the extra 6 components, with respect to the metric, could describe the components of the electromagnetic field. Nevertheless, it was shown that these extra 6 components are linked to the local Lorentz invariance of the theory [3,4].
Despite the failure of the Weyl and Einstein attempts, the new approaches supplied the notion of gauge theory and thus the search for a gauge gravity started [5,6]. In 1979, K. Hayashi and T. Shirafuji [7] proposed a gravitational theory, called "New General Relativity" that is a gauge theory for the translation group. This theory involves three free parameters to be determined by the experiment.
In this situation, a good approach could be testing different formulations of gravitational interaction by searching for signatures discriminating among concurrent theories. Systems, in strong field regime, are natural candidates to this aim. Signatures can come from compact objects like neutron stars or black holes, or from polarizations of gravitational waves [10,[23][24][25].
From a physical point of view, compact objects are, in general, stars exhausting their nuclear fuel. They can give rise to stable compact objects or to black holes. For example, neutron stars are compact objects which are boosted by their neutron degeneracy pressure against the attraction of gravity. Another type of compact stars are the white dwarfs, boosted by the electron degeneracy pressure against the gravity.
From a theoretical viewpoint, it is well know that the first exact vacuum solution of GR was the Schwarzschild one [26]. Thereafter many solutions investigating compact stars have been proposed. In particular, searching for interior solutions, describing realistic compact objects, became a fertile domain to probe GR and other theories. Using the equation of state (EoS), one can study the stability structure of compact stellar objects and then guess on their internal composition. Specifically, the EoS is useful to investigate the physical behavior of the stellar structure through the Tolman-Oppenheimer-Volkov (TOV) equation which is the general relativistic equation of the stellar interior coming from the internal Schwarzschild solution.
In order to study stellar configurations, one can assumes, at the beginning, that the distribution of matter is isotropic which means that the radial and tangential pressures are equal. However, in realistic cases, such assumption does not hold and one finds that the two components of radial and tangential pressures are not equal: the differences between them create the anisotropy. Such stellar configurations have unequal radial and tangential pressures. Lemaître, in 1933, was the first which proposed anisotropic models [27]. Moreover, it possible to show that, in order to reach stable configurations, at the maximal star surface, the radial pressure has to decrease and, at the center, it has to vanish [28,29].
There are several factors to be taken into account to study anisotropic stars; among these, the high density regime where the nuclear interactions have to be relativistically treated [30,31]. Moreover the existence of a solid core or a 3A type superfluid may cause the star to be anisotropic [32]. There is another source that makes the star anisotropic: it is a strong magnetic field [33,34]. The slow rotation can be considered another source of anisotropy [35]. Letelier showed that combinations between perfect and null fluids can give rise to anisotropic fluids [36]. Other reasons can be taken into account to generate anisotropies like pion condensation [37], strong electromagnetic field [38] and phase transition [39].
Dev and Gleiser [40,41] and Gleiser and Dev [42] investigated the operators that induce the pressure to be anisotropic. It was shown that the effect of shear, electromagnetic field, etc. on self-bound systems can be neglected, if the system is anisotropic [43]. Systems that consist of scalar fields like boson stars possess anisotropy [44]. Gravastars and wormholes can be considered as anisotropic models [45][46][47][48]. An application of anisotropic model to stable configurations of neutron stars has been discussed in [49]. They showed that anisotropy might have non-negligible effects on the equilibrium mass and on the surface red-shift.
A nice study describing the origin and the effects of anisotropy can be found in [50,51]. Super dense and anisotropic neutron stars have been considered and a conclusion is that there is no limiting mass of such stars [52].
Supermassive neutron stars in alternative gravity are considered, for example, in [53][54][55][56]. The presence of torsion field in neutron star models is studied in [57]. The issue of stability of anisotropic stars has been analyzed in [58] and a conclusion is that the stability of such systems is similar to that of isotropic stars.
There are many anisotropic models that deal with the anisotropic pressure in the energy-momentum tensor. Several exact spherically symmetric solutions of interior stellar have been developed .
Beside this, adding charge effects in compact objects is an issue widely considered in literature, especially for neutron and quark stars [96][97][98][99][100][101]. An important result is the fact that the distribution of net charge can improve the maximum mass of compact stars. For example for white dwarfs, there are several results related to the charge effects affecting the structure (see, e.g. [102]). For this reason, we consider here charged compact star models showing that the net charge effects can increase the maximum mass.
The aim of the present paper is to derive a novel charged anisotropic solution in the framework of TEGR and compare it with realistic stellar configurations using physical assumptions on the form of metric potential and the combination of charge and anisotropy. In particular, we consider the pulsar PSR J 1614 -2230 which estimated mass M = 1.97 ± 0.04M ⊚ and radius is R = 9.69 ± 0.2 km [103]. This peculiar system escapes the standard GR explanation of neutron stars because it is too massive to be stable unless one assumes exotic EoS or alternative gravities. We discuss the physical parameters of such a pulsar considering our solution derived in the framework of TEGR. From our point of view, this can constitute a possible test for the theory.
The set up of the paper is the following. In Sec. II, we sketch the TEGR theory in presence of the electromagnetic field. Sec. III is devoted to the discussion of charged compact stars in TEGR. The requirements for a physically consistent stellar model in TEGR are discussed in Sec. IV. The physical properties of the model are considered in Sec.V. The model is matched with realistic compact stars, in particular with PSR J 1614-2230, in Sec. VI. In Sec.VII, the stability of the model is taken into account with respect to the TOV equation and the adiabatic index. In particular, we report the (M, R) relation and how it is affected by the electromagnetic field. Discussion and conclusions are reported in Sec. VIII.

II. TELEPARALLEL EQUIVALENT OF GENERAL RELATIVITY AND THE ELECTROMAGNETIC FIELD
In TEGR theory, at each point of the spacetime manifold, M , can be defined a tangent Minkowski spacetime η ij = (+1, −1, −1, −1) 1 . As is well-known that the dynamical fields of TEGR are the four linearly independent vierbeins (tetrads) from where we can define the metric and its inverse as with η ab being the flat Minkowski metric of the tangent space and l b ν , L b ν are the covariant and contra-variant tetrad fields. Using Eq. (1) we can define the following orthogonal conditions From the above definitions, we can define the torsion tensor as 2 where ω a bµ is the spin connection which we can set equal zero due to TEGR. Therefore the torsion tensor (3) reduces to TEGR with electromagnetic field is recovered from the action where κ 2 = 8π, |l| = det(l a µ ) = √ −g is the determinant of the tetrad field, F is a gauge-invariant scalar defined as [104,105]. The torsion scalar T is defined as where S a µν is the superpotential defined as in terms of the contortion tensor The metric components are functions of the tetrads. The variation with respect to the tetrad l a µ yields the field equations (see [106] for details) The stress-energy tensor T a µ consists of two terms where 1 Latin indices represent spacetime coordinates. Greek ones describe tangent space indices. 2 Square brackets means anti-symmetrization, i.e., is the energy-momentum tensor of an anisotropic fluid, and is the energy-momentum tensor of the electromagnetic field. Here u µ is the time-like vector defined as u µ = [1, 0, 0, 0] and ξ µ is the unit space-like vector in the radial direction defined as ξ µ = [0, 1, 0, 0] such that u µ u µ = −1 and ξ µ ξ µ = 1.
In this paper, ρ is the energy-density, p r and p t are the radial and the tangential pressures respectively. Furthermore, the electromagnetic tensor F µν satisfies the Maxwell equations where J µ = σu µ is the current density and σ is the charge density. In the next section, we are going to apply the field equations (9) and (13) to a spherically symmetric tetrad space and try to solve the resulting system of differential equations.

III. CHARGED COMPACT STARS
Let us begin with the following spherically symmetric metric using the spherical coordinates (t, r, θ, φ) where G(r) and H(r) are two unknown functions depending on the radial coordinate r. The above metric (14) can be reproduced from the following covariant tetrad field [107] The tetrad (15) is the output of the product of a diagonal tetrad and a local Lorentz transformation which can be written as It is worth noticing that the diagonal tetrad can be applied to the field equations of f (T ) TEGR models [107]. Here we shall take into account the effect of local Lorentz transformation given by Eq. (16) and see if the results are different from those presented in [108].
Using Eq. (15) into Eq. (6), the torsion scalar takes the form Despite the fact that the diagonal tetrad in Eq. (16) and tetrad (15) reproduce the same metric, their torsion scalars are different. The reason for this difference is the Local Lorentz transformation (LLT) in (16). The torsion scalar of Eq. (17) goes to zero when both the unknown functions G(r) and H(r) tend to one, which is the necessary condition to achieve the asymptotic flatness. Nevertheless, the torsion scalar presented in [108] does not vanish as it should when we apply the same condition. So the existence of LLT makes our torsion scalar more physically motivated. Using Eq. (17) into field Eqs. (9) and (13), we get where the prime indicates derivatives w.r.t the radial coordinate r and ∆ = 8π(p t − p r ) is the pressure difference. We used the energy-momentum tensor T µ ν to get the form where E 2 (r) is the electric field defined as with q(r) being the electric charge inside a sphere of radius r and σ is the charge density (for details see, e.g., [109,110]).
Here ∆(r) is the anisotropic parameter of the stellar system. Field Eqs. (18) coincide with those given in [108], and can be recast also in the field equations presented in [111] and in [112]. The above differential system are five independent equations in seven unknown functions, G, H, ρ, p r , p t , E(r) and ∆(r). Therefore, we need two extra conditions to solve the above system. One of these extra conditions is assuming the metric potential g rr having the form where k is a constant having the dimension of the length. It will be determined from the matching conditions. Eq. (21) shows that, for r = 0 and H(r) = 1, H(r) is finite at the center. Also the derivative of H is finite at the origin. The second condition comes from the use of Eq. (21) in ∆(r) which takes the form where we have put E 2 (r) = k 2 1 r 2 . Eq. (22) is a second order differential equation in the unknown G function. In order to solve Eq. (22), we assume the anisotropic ∆ to have the form Eq. (23) has not been used before for charged stars in GR and TEGR. The present aim is to study the effect of anisotropy, given in (23), in presence of electric charge, on the neutron star structure. Inserting Eq. (23) into Eq. (22), we get The solution of differential equation (24) is Clearly, G goes to a constant value for r → 0 and r → ∞. From Eqs. (21) and (25) of the system of differential Eqs. (18), we get the remaining unknown functions If the constant k 1 = 0, our star is neutral. In this case, Eqs. (26) will be identical with those presented in [108]. Eqs (26), at the boundary of the star, take the form where R is the radius at the boundary of the star. From the second equation of (27), we get the condition for the vanishing of radial pressure, that is For R = 8.51M ⊙ , which is the boundary radius of the star HerX − 1 , c 2 takes the form If we substitute k = 33.43508140, k1 = 0.008680552960 in Eq. (30), we get c 2 = 0.005. In other words, if we substitute R = 8.51M ⊙ , k = 33.43508140 and k1 = 0.008680552960 in the second equation of (27), we get a vanishing radial pressure for the pulsar Her X-1. It is important to mention that the anisotropic force is defined as 2∆ r and it is attractive for p r − p t > 0 and repulsive for p r − p t < 0. The mass contained within a radius r of the sphere is defined as Using Eq. (26) in Eq. (31), we get The compactness parameter of a spherically symmetric source with radius r takes the form [108] u(r) = 2M (r) r = r 2 (20k 6 + 20k 2 r 4 − 5r 6 − 30k 4 r 2 + 2k 8 k 1 2 r 2 ) 5k 8 .
In the next section, we are going to discuss the physical requirements to derive viable stellar structures and to see if model (26) satisfy them or not.

IV. REQUIREMENTS FOR A PHYSICALLY CONSISTENT STELLAR MODEL
A physically viable stellar model has to satisfy the following conditions throughout the stellar configurations: • The gravitational potentials G(r) and H(r), and the matter quantities ρ, p r , p t have to be well defined at the center and regular as well as singularity free throughout the interior of the star. • The energy density ρ has to be positive throughout the stellar interior i.e., ρ ≥ 0. Its value at the center of the star should be positive, finite and monotonically decreasing towards the boundary inside the stellar interior, that is dρ dr ≤ 0.
• The radial pressure p r and the tangential pressure p t must be positive inside the fluid configuration i.e., p r ≥ 0, p t ≥ 0. The gradient of the pressure must be negative inside the stellar body, i.e., dpr dr < 0, dpt dr < 0. At the stellar boundary r = R, the radial pressure p r has to vanish but the tangential pressure p t may not be zero at the boundary. At the center both pressures are equal. This means the anisotropy has to vanish at the center, that is ∆(r = 0) = 0.
• For an anisotropic fluid sphere, the fulfillment of the energy conditions refers to the following inequalities in every point inside the fluid sphere: (i) Weak energy condition (WEC): p r + ρ > 0, ρ > 0.
• Causality condition has to be satisfied to get a realistic model i.e. the speed of sound must be smaller than 1 (assuming the speed of light c = 1) in the interior of the star, i.e. 1 ≥ dpr dr ≥ 0, 1 ≥ dpt dr ≥ 0. • The interior metric functions have to match smoothly the exterior Schwarzschild metric at the boundary.
• For a stable model, the adiabatic index should be greater than 4 3 . • Herrera method [113] to study the stability of anisotropic stars suggests that a viable model should also satisfy where v r and v t are the radial and transverse speed respectively. Now we are going to analyze the above physical requirements in details to see if our model satisfy them.

V. PHYSICAL PROPERTIES OF THE MODEL
Let us test the model (26) and see if it is consistent with realistic stellar structures. To this aim, we discuss the following issues: A. Non-singular model i-The metric functions of this model satisfy, which means that the gravitational potentials are finite at the center of the stellar configuration. Moreover, the derivatives of these potentials are finite at the center, i.e., G ′ r=0 = H ′ r=0 = 0. The above conditions means that the metric is regular at the center and has a well behavior throughout the interior of the stellar.
ii-Density, radial and tangential pressures of (26), at the center, have the form Eqs. (35) show that the density is always positive and the anisotropy is vanishing at the center. The radial and tangential pressures have a positive value as soon as c 2 > 0 otherwise they become negative. Moreover, the Zeldovich condition [114] states that the radial pressure must be less than or equal to the density at the center i.e., pr(0) ρ(0) ≤ 1. Using the Zeldovich condition in Eq. (35), we get iii-The derivative of energy density, radial and tangential pressures of the model are respectively: where ρ ′ = dρ dr , p ′ r = dpr dr and p ′ t = dpt dr . Eqs. (37) show that the gradients of density, radial and tangential pressures are negative as we will show when we plot them. iv-The radial and transverse velocity of sound (c = 1) are obtained as v r 2 = dp r dρ

B. Matching conditions
We assume that the exterior spacetime for a not-rotating star is empty and described by the exterior Reissner-Nordström solution that is a solution of vacuum teleparallel gravity as shown in [115][116][117]. It has the form where M is the total mass, Q is the charge and Q 2 (r−2M) < 1. We have to match the interior spacetime metric (14) with the exterior Reissner-Nordström spacetime metric (39) at the boundary of the star r = R. The continuity of the metric functions across the boundary r = R gives the conditions in addition to the fact that radial pressure approaches to zero at a finite value of the radial parameter r which coincides with the radius of the star R. Therefore, the radius of the star can be obtained by using the physical condition p r (r = R) = 0. From the above conditions, we get the constraints on the constants c 1 , c 2 and k. The constant k, from these conditions, is where The constants c 1 and c 2 are lengthy and useless to be reported here. We shall write their numerical values when confronting with observational data.

VI. MATCHING THE MODEL WITH REALISTIC COMPACT STARS
Let us consider now the previous physical conditions of the model derived to test it by using masses and radii of observed pulsars. In order to support our model, we will study the pulsar PSR J 1614 -2230 whose estimated mass and radius are M = 1.97 ± 0.04M ⊚ and R = 9.69 ± 0.2 km, respectively [118]. We can use the maximal values M = 2.01M ⊚ and R = 9.89km as input parameters. The boundary conditions are adopted to determine the constants k = 28, c 1 = 174908.5108 and c 2 = 0.0003. Adopting these constants, we can plot the physical quantities. The regular behavior of these one can be assumed as a first requirement to fit a realistic star model.  Density, radial and tangential pressures, given in Eq. (26), are plotted in Fig. 2. Fig. 2 shows that density, radial and tangential pressures are positive as required for realistic stellar configuration. Moreover, as Fig. 2 (a) shows, the density is high at the center and decreases far from it. Fig. 2 (b) shows that the radial pressure goes to the zero at the boundary while the tangential pressure remains non-zero at the boundary. Also this feature is relevant for a realistic model. Fig. 3 (a) shows that the anisotropy is vanishing at the center and increases at the surface of the star. Specifically,   Fig. 3 (b) shows that the anisotropic force ∆ r is positive. This means that it is repulsive due to the fact that p t −p r > 0. Fig. 4 shows that the gradients of density, radial and tangential pressures are negative which confirms the decreasing of density, radial and transverse pressures through the stellar configuration.
In Fig. 5, the radial and tangential speeds of sound are reported. They are positive and both of them are less than one. This result confirms the non-violations of causality condition in the interior of the star. Fig. 6 represent the behavior of the energy conditions. In particular, Figs. 6 (a), (b), (c) and (d) show the positive values of the WEC, NEC, SEC and DEC energy conditions. Therefore, all the energy conditions are satisfied throughout the stellar configuration as required for a physically meaningful stellar model.
The matching of the interior and exterior metrics at the boundary are shown in Fig. 7. Fig. 7 (a) represents the smooth matching between G(R) and 1 − 2M R + Q 2 R 2 . Fig. 7 (b) shows the smooth matching between H(R) and    (1 − 2M R + Q 2 R 2 ) −1 . In Fig. 8, we have plotted radial and tangential EoS. As Fig 8 (a), and (b) show, the EoS is not linear. In Ref. [119], authors derive the EoS for neutral compact stars and show that it is almost a linear one. Here, both the radial and tangential EoS show non-linear form which is due to the contribution of the electric field.
The mass function given by Eq. (31) is plotted in Fig. 9 which shows that it is a monotonically increasing function of the radial coordinate and M (r = 0) = 0. Moreover, Fig. 9 shows the behavior of the compactness parameter of star which is increasing.
The radial variation of the surface red-shift is plotted in Fig. 10. Böhmer and Harko [72] constrained the surface red-shift to be Z ≤ 5. The surface redshift of this model is calculated according to P SRJ1614 − 2230 and found to  be 0.186927610.
Using Eqs. (26) and (32), it is possible to derive the Mass-Radius (M, R) relation of the model for a given value of the surface density (ρ(r = R) = 1.5 ≈ 10 15 g cm 3 ). In Fig. 11, it is reported considering also the compactness-radius relation. As Fig. 11 (c) shows, our model has a maximum mass 3M ⊚ which is well beyond the recently reported values of 2.50-2.67M ⊚ recently reported by the LIGO collaboration [120]. This means that anomalous compact objects can be addressed in the framework of TEGR.

VII. STABILITY OF THE MODEL
In this section we are going to discuss the stability issue using two different techniques; the Tolman-Oppenheimer-Volkoff equations and the adiabatic index. In this subsection we are going to discuss the stability of the model. To this goal, we assume hydrostatic equilibrium through the Tolman-Oppenheimer-Volkoff (TOV) equation. Using the TOV equation [121,122] as that presented in [123], we get the following form with M g (r) being the gravitational mass at radius r, as defined by the Tolman-Whittaker mass formula which gives Inserting Eq. (43) into (42), we get where , F a = 2(pt−pr) r , F h = − dpr dr and F e = σE √ H are the gravitational, the anisotropic, the hydrostatic and the electromagnetic forces respectively. The behavior of the TOV equation for the model (26) is shown in Fig. 12. The four different forces are plotted in Fig. 12. It shows that hydrostatics, anisotropic and electromagnetic forces are positive and dominated by the gravitational force which is negative to keep the system in static equilibrium.

B. Stability in the static state
It was shown by Harrison, Zeldovich and Novikov [114,124] that, for stable compact stars, the gradient of central density, with regard to the mass increasing, must be positive, i.e., ∂M ∂ρ r 0 > 0. If this condition is verified, we have stable configurations. Specifically, stable or unstable regions are separated when we have a constant mass i.e. ∂M ∂ρr 0 = 0. Now we are going to apply this condition to our solution (26). For this purpose, we calculate the central density for solution (26) and get Using Eq. (45) in Eq. (32), we get With the help of Eq. (46), we have From Eq. (47), it is crear that solution (26) has a stable configuration since ∂M ∂ρr 0 > 0. The behavior of (46) and (47) are show in Fig. 13. Fig. 13 shows that mass increases as the energy density increases and the mass gradient decreases as energy density increases.
In Fig. 14 we have reported Γ, Γ r and Γ t respectively. As it is clear from these plots, it can be seen that the values of Γ r and Γ t are greater than Γ throughout the stellar interior and hence the stability condition is fulfilled.   Besides PSR J 1614-2230, a similar analysis can be developed for other pulsars. In Table I and II, we report the results for other observed systems.

VIII. DISCUSSION AND CONCLUSIONS
An important remark is in order at this point. It is well known that TEGR theory is equivalent to GR up to a total derivative term [128][129][130]. In this paper, we considered a combination of anisotropy and charge in TEGR equations. This situation gives rise to the effects of enhancing mass and modifying the (M, R) relation of GR.
For a class of metric potentials and anisotropy functions, we derived an exact solution capable of figuring out realistic compact star configurations. The regularity conditions of the solution at the origin as well as at the surface of the star show a well behavior throughout the stellar structure. This is different from the results reported in [108] for GR. In that study, the authors show that pressure increases outward which is a non-physical situation. The difference between our results and those reported in [108] is due to the anisotropy given in Eq. (23), the presence of charge, and the non-vanishing of the radial pressure. In [108], they assumed a vanishing radial pressure. In our case, we showed that density, radial and tangential pressures behave regularly according to the observational data of the pulsar PSR J 1614-2230, the first reported very massive neutron star, whose existence ruled out many EoS [103]. In order to explain such a system, exotic matter such as hyperons and kaon condensates, alternative theories of gravity and other hypotheses have been invoked (see e.g. [53][54][55]). The approach seems to work for other systems, as reported in Tables I and II.
Furthermore, we show that the anisotropy has a positive value which can be interpreted as a repulsive force. This fact is because the tangential pressure is greater than the radial pressure, i.e., p t > p r [131]. The issue of stability is studied and we showed that the derived model is stable against the different forces (gravitational, hydrostatic, anisotropic and electromagnetic) acting on it. We also calculated the sound of speed and showed that it is consistent with realistic compact stars in contrast to the analogue charged models formulated in GR, where an imaginary sound speed is derived [108]. Finally we calculated the adiabatic index of our model and showed that also it represents a realistic physical star. It is worth noticing that the adiabatic index presented in [108] has a negative value which is not consistent with realistic stellar models. This indicates, in a clear way, that our assumption of the metric potential (21) and the anisotropy form (23) are physical assumptions that makes the resulting stellar model consistent with real stellar objects. We tested the model over a wide range of reported observed values of masses and radii of pulsars (Tables I and II). The conclusion is that the fit is good also in these cases. Finally, we drew the mass-radius relation and showed the effect of electric field on it.
It is shown that the electric charge plays a central role in improving the results compared with the neutral case. Among these improvements, we have: i-In the neutral case, one gets a maximum mass as 2.788M ⊙ while, in the charged case, the maximum mass becomes 3M ⊙ . ii-In the neutral case, one get an increasing radial pressure as reported in [108] while in the charged case, one gets a decreasing one which describes a consistent compact star. iii-In the neutral star, one gets an imaginary sound speed [108] while, in the charged case, we get a real physical sound speed as shown in Fig. 5.
The approach can be summarized as follows: we used a non-diagonal form of the tetrad field that gives a null value of the torsion as soon as the metric potentials approach to 1. This is a necessary condition for any physical tetrad field as reported in various studies [132][133][134][135][136][137][138].
We can conclude that the comparison of our exact solution with the physical parameters of pulsars gives indications that the model can realistically represent observed systems. Furthermore, the approach can be extended to a large class of metrics and anisotropies, if the above physical requirements are satisfied. However, a detailed confrontation with observational data is needed. This will be the argument of a forthcoming paper.