Anisotropic strange stars in the Einstein-Maxwell spacetime

We present here a detailed analysis on the effects of charge on the anisotropic strange star candidates by considering a spherically symmetric interior spacetime metric. To obtain exact solution of the Einstein-Maxwell field equations we have considered the anisotropic strange quark matter (SQM) distribution governed by the simplified MIT bag equation of state (EOS), $p=\frac{1}{3}\left( {\rho}-4\,B \right)$, where $B$ is the bag constant and the distribution of the electrical charge is given as $q(r)=Q\left({r}/{R}\right)^3=\alpha {r^3}$, where $\alpha$ is a constant. To this end, to calculate different constants we have described the exterior spacetime by the Reissner-Nordstr{\"o}m metric. By using the values of the observed mass for the different strange star candidates we have maximized anisotropic stress at the surface to predict the exact values of the radius for the different values of $\alpha$ and a specific value of the bag constant. Further, we perform different tests to study the physical validity and the stability of the proposed stellar model. We found accumulation of the electric charge distribution is maximum at the surface having electric charge of the order ${{10}^{20}}~C$ and electric field of the order ${10}^{21-22}~V/cm$. To study the different physical parameters and the effects of charge on the anisotropic stellar system we have presented our analysis graphically and in the tabular format by considering $LMC~X-4$ as the representative of the strange star candidates.

stant. To this end, to calculate different constants we have described the exterior spacetime by the Reissner-Nordström metric. By using the values of the observed mass for the different strange star candidates we have maximized anisotropic stress at the surface to predict the exact values of the radius for the different values of α and a specific value of the bag constant. Further, we perform different tests to study the physical validity and the stability of the proposed stellar model. We found accumulation of the electric charge distribution is maximum at the surface having electric charge of the order 10 20 C and electric field of the order 10 21−22 V /cm. To study the different physical parameters and the effects of charge on the anisotropic stellar system we have presented our analysis graphically and in the tabular format by considering LM C X − 4 as the representative of the strange star candidates.
Keywords General Relativity; anisotropic fluid; strange stars.

Introduction
The theoretical possibility of the existence of hypothetical strange quark stars were first speculated in Refs. [1,2,3,4]. According to the strange quark matter hypothesis [5,6,7] the strange quark matter (SQM), made of equal number of up, down and strange quarks can be considered as the absolute ground state for the confined state of hadrons [2,3,8,9]. Although strange stars form a distinct hypothetical branch of compact stars but these heavier members have masses and radii quite similar to the neutron stars. However, strange stars are not part of the the continuum of equilibrium configurations like white dwarf and neutron stars [10,11,12]. In this context it is worth mentioning that strange matter equation of state (EOS) appears as the suitable EOS to explain observed compactness of the compact astrophysical objects like 4U 1820−30, SAX J 1808.4−3658, 4U 1728 − 34, Her X − 1, RX J185635 − 3754 and P SR 0943 + 10 [3,4,13,14,15,16], whereas neutron star EOS failed to explain those estimated compactness.
To maintain global charge neutrality strange stars which made of approximately equal number of up, down and strange quarks, should include smaller number of electrons. Alcock et al. [3,8] and Usov et al. [17,18] in their study showed that high electric fields in the order of 10 18−19 V /cm is expected to present on the surface of strange stars and presence of electrons play a significant role to the formation of the electric dipole layer at the surface. Such strong electric fields have values on the order of the energy density of SQM and it should be included in the stress-energy tensor which describes strange stars. The presence of the charge affects the relativistic stellar system in the following ways: (i) it causes the space-time curvature, (ii) it produces Coulomb interaction by introducing an extra term in the relativistic hydrodynamic equilibrium equation, and (iii) the energy density associated with the electric field has significant role in producing the gravitational mass of the relativistic stellar system. In this line several literature [19,20,21,38,23,41,25,26,27,28] can be referred to understand the effects of the electric charge on the relativistic compact stellar system.
In his pioneering work Ruderman [29] first introduced the idea of pressure anisotropy and showed that the high density of the nuclear matters which interact relativistically are the key reason of the formation of the anisotropy. Here, by anisotropy we are addressing the difference between the radial component, p r (r) and the angular component, p θ (r) = p φ (r) ≡ p t (r) of the pressure. Clearly, p θ (r) = p φ (r) is the consequence of the assumed spherical symmetry of the stellar system. An extensive study by Bowers and Liang [30] showed that in the presence of complex strong interactions anisotropy in the spherically symmetric stellar system may be arising due to the presence of superconductivity and superfluidity of the ultradense matter. Later, Herrera and Santos in their detailed review [31] discussed the possible reasons behind the formation and existence of the local anisotropic stress in a self gravitating system and also studied their effect on a static spherically symmetric stellar system. Dev and Gleiser [32,33,34] in their series of work studied the significant effect of anisotropy on the redshift and maximum mass. They also showed that the presence of anisotropic stress enhances stability of the relativistic stellar system compared to the isotropic cases and predicted that for the lower adiabatic index values too anisotropic systems are stable. In this line several authors as in Refs. [35,36,37,17,38,39,40,41,42,43,44,45,46,47] have studied the effect of anisotropy on the spherically symmetric compact stellar system and examined it's effect on the different physical properties of the stellar system. The present work is the charged generalization of the earlier work done by Deb et al. [45], where they presented an unique anisotropic model for the strange stars and showed the typical mass-radius relation for the strange stars by solving the Einstein field equations. To this end, they assumed simplified MIT bag EOS and showed that maximum anisotropy at the surface of the ultra dense strange stars is their inherent property. Using the motivation of the earlier work [45], in the present article we have studied charged and anisotropic spherically symmetric stellar systems for the strange stars by considering a specific form of the electric charge distribution, q(r). We also presented exact solutions for the Maxwell-Einstein field equations. It is interesting to note that though there are several literature which separately studied the effect of anisotropy or charge on the strange stars. But we found there is no other literature, which has studied the combined effects of anisotropy and charge on the compact stellar system by providing typical mass-radius relation for the strange stars in the framework of the Maxwell-Einstein gravity. However, in the present study by considering the combined effects of anisotropy and charge on the stellar system we attempt to present the exact solutions for the Maxwell-Einstein field equations by providing the typical mass-radius relation for the strange stars. Further, we have also examined the physical validity of the obtained solutions.
The outline of our study is as follows: In Sect. 2 we have presented the basis of using the MIT bag EOS and the chosen form of the electric charge distribution. The basic equations to describe the anisotropic charged stellar system are presented in Sect. 3. In Sect. 4 we have derived the solutions for the Maxwell-Einstein field equations and presented expressions for the different physical parameters. Further, in Sect. 5 to show physical acceptibility of the stellar system on the basis of the obtained solutions we have performed different tests like Energy conditions 5.1, mass-radius relation 5.2, compactification factor and redshift 5.3, and the stability of the system 5.4. Finally, in Sect. 6 we have concluded our study by discussing in detail the effects of the electric charge distribution on the anisotropic stellar system.

The MIT Bag equation of state and the electric charge distribution
In the present article we consider MIT bag model EOS [48] to describe the SQM distribution. in MIT bag model to maintain all the corrections due to energy and pressure functions of SQM an ad hoc bag function has been introduced. For the simplicity we assume that the up (u), down (d) and strange (s) quarks are are massless and non-interacting in nature. Hence, the quark pressure, p r is defined as where p f is the pressure due to individual quark flavors viz. u, d and s. B is the vacuum energy density and usually known as 'Bag constant'. The relation between p f and energy density due to each quark flavors reads p f = 1 3 ρ f . Hence the energy density, ρ due to deconfined SQM distribution inside the bag is defined as Hence, substituting relation between p f and ρ f into Eq.
To study the effects of charge on the relativistic stellar system Felice et al. [58,59] in their literature considered an specific form of electric charge distribution q(r) given as q(r) = Q(r/R) n . Following Felice et al. [58,59] in the present study we choose this specific simplest form of q(r), for the parametric values of n = 3 as follows where Q and R are the total charge and the total radius of the stellar system, respectively and α is a constant which can be defied as α = Q/R 3 .

Basic stellar structure equations
To describe interior spacetime of the ultra dense spherically symmetrical stellar system in Schwarzschild-like coordinates [60,61] we use metric as follows where the metric potential ν and λ are the functions of the radial coordinate r only. Now, to obtain hydrostatic stellar structure of the charged sphere we have to solve the Einstein-Maxwell field equations provided as where we assumed G = 1 = c in the relativistic geometrized unit. Here T i j and E i j represents stress-energy tensor for the locally anisotropic fluid distribution and the present electromagnetic field, respectively and they are defined as [62] where v i and u i are the four-velocity and radial fourvector, respectively, ρ is the energy density of SQM, p r represents pressure in the direction of u i , known as radial pressure and p t represent component of pressure normal to u i , known as tangential pressure. Here, F ij is the anti-symmetric electromagnetic field tensor and can be defined as where, A j = (φ(r), 0, 0, 0) is the four-potential. F ij satisfies the covariant Maxwell equations, where J i is the electromagnetic four-current vector defined as where σ = e ν/2 J 0 (r) represents the charged density and g is the determinant of the metric g ij defined by For a static spherically symmetric stellar system J 0 is the only non vanishing component of the the electromagnetic four-current J i which is a function radial coordinate, r only. F 01 and F 10 are the only non zero components of the electromagnetic field tensor and they are related by F 01 = −F 10 . F 01 and F 10 are the radial component of the electric field. Using Eqs. (11) and (12) the expression for the electric field is given as If q(r) represents the total charge of a spherical system of radius r then following the relativistic Gauss's law the electric charge q(r) can be defined as Using Eqs. (7), (8), (11), (12), (14) and (15) the stress-energy tensor for the anisotropic charged matter distribution can be written as where the electric charge and the electric field are related by q 2 (r)/8πr 4 = E 2 (r)/8π. Substituting Eq. (16) into Eq. (6) we have the explicit form of the Einstein field equation for the anisotropic charged spherically symmetric stellar system as follows [62] e −λ λ ′ r − 1 4 In the analogy of the uncharged case, let we define the mass function of the spherically symmetric charged stellar system as follows where ρ ef f = ρ + E 2 8π . To describe the exterior spacetime of our system we consider the exterior Reissner-Nordström metric given as Now using Eqs. (4), (20) and (21) we find from Eqs. (17) as follows Following Mak and Harko [63] to obtain singularity free monotonically decreasing SQM density function ρ, we define where ρ c and ρ 0 are the central and surface density, respectively.
To obtain hydrostatic equilibrium equation for the anisotropic charged stellar system we perform covariant divergence of the electromagnetic stress-energy tensor, i.e., ∇ a T a b = 0, which leads to the equation of energy conservation as follows For q = 0 in Eq. (24) where λ 1 , λ 2 , ν 1 , ν 2 , ν 3 , ν 4 , ν 5 and ν 6 are constants and their expressions are shown in Appendix 6. We featured variation of the physical parameters, viz. e λ , e ν , ρ, p r and p t with respect to the radial coordinate r/R in Figs. 1 and 2.
The anisotropy for our system is given as The variation of the anisotropic stress (∆) with respect to the radial coordinate r/R is shown in Fig. 3. We find from Fig. 3 anisotropy is zero at the center and maximum at the surface as predicted by Deb et al. [45].  Table 1 and following Deb et al. [45] we shall maximize the anisotropic stress, ∆(r) at the surface r = R to predict the exact value of the radius, R for the different strange stars. To this end, we consider that the value of the bag constant is B = 83 M eV /f m 3 [57] and the chosen values of α are 0, 0.0005, 0.0010 and 0.0015 km −2 . Clearly, ∆ ′ (R) = 0 will yield several values of R and we will choose only that value of R for which the Buchdahl conditions [64] will be satisfied.
In Fig. 4 we have featured the variation of the electric charge distribution q(r) and electrical energy density E 2 (r)/8π with respect to the radial coordinate r/R in the upper and lower panel, respectively. Fig. 4 clearly suggests that both the distribution of the electric charge and electrical energy density is minimum, i.e., zero at the center and maximum at the surface. In this section to discuss physical validity of the achieved solution we will study some salient physical features of the stellar system as follows:

Energy conditions
To satisfy energy conditions, viz., Null Energy Condition (NEC),Weak Energy Condition (WEC), Strong Energy Condition (SEC) and Dominant Energy Condition (DEC) the anisotropic charged stellar system have to be consistent with all the inequalities simultaneously as follows

Mass-radius relation
Andréasson [65] predicted the upper bound of the massradius ratio for the charged spherically symmetric stellar system, which was generalization of the Buchdahl limit [64] that provides upper limit for the allowed massradius ratio in the uncharged case. Hence, in the present system the upper bound [65] is given as The mass function for our system is provided as follows We have presented variation of the total mass, M (normalized in solar mass, M ⊙ ) with respect to the total radius, R due to different parametric values of α in Fig. 6, where we chose that the bag constant is B = 83 M eV /f m 3 [57]. We find the maximum mass of the system increases as the value of α increases, which is clearly shown in the lower panel of Fig. 6.  Fig. 7 show the variation of R with respect to ρ c . We find, R Mmax , the radius corresponding to M max

Compactification factor and redshift
The compactification factor for our system is defined as Hence, the surface redshift, Z s corresponding to the compactification factor u is given as The variation of the redshift function, Z(r) with respect to the radial coordinate r/R is presented in Fig. 8. Clearly, in a spherically symmetric anisotropic charged stellar system as the value of α increases the values of the surface redshift gradually decreases.

The stability of the system
To examine stability of our system we will study i) Generalized TOV equation and ii) Herrera cracking concept as follows

Generalized TOV equation
The generalized form of the TOV equation in the present anisotropic charged system reads where M g denotes the effective gravitational mass and given as follows Eq. (39) features that the system is completely stable under the equilibrium of the different forces, i.e., F g + F h + F e + F a = 0, where F g , F h , F e and F a represents gravitational, hydrodynamic, electric and anisotropic force, respectively. We have presented variation of the different forces with respect to the radial coordinate r/R due to different values of α in Fig. 9. The figure features that the attractive gravitational force F g , which acts toward the inward direction along the system is counterbalanced by the combined effects of the forces F h , F e and F a .

Herrera cracking concept
To examine stability of the system in terms sound speeds the systems have to be consistent with the i) causality condition and ii) Herrera cracking concept. To be consistent with the causality condition, the square of the radial (v 2 sr ) and tangential (v 2 st ) sound speeds should satisfy the inequalities 0 ≤ v 2 sr ≤ 1 and 0 ≤ v 2 st ≤ 1 simultaneously. According to the concept of Herrera's cracking [66,67] for a potentially stable region v 2 sr should be greater than v 2 st and the difference of the square of the sound speeds should maintain it's sign same through In Fig. 10 we have shown the variation of the square of the sound speeds (upper panel) and |v 2 st − v 2 sr | (lower panel) with respect to the radial coordinate r/R due to different parametric values of α. The figure clearly features that both our system satisfies both the causality condition and Herrera cracking concept. Hence, our system is completely stable.

Discussion and conclusion
In this literature we have presented a detailed study on the effect of the electrical charge on the spherically symmetric anisotropic stellar system, which is made of SQM and governed by the MIT bag EOS. Assuming a simplified form of the electrical charge distribution given as q (r) = Q (r/R) 3 ≡ α r 3 , we have obtained exact solutions for the Maxwell-Einstein field equations. Further, using the exterior Reissner-Nordström metric we have presented expressions for the different physical parameters in Eqs. (25)- (30). We have presented the obtained solutions and studied their physical validity in terms of the star LM C X − 4 of mass 1.29 M ⊙ , by considering it as the representative of strange star candidates. Throughout the study we have considered bag The profile of the metric potentials e ν , e λ are shown in Fig. 1, which shows that at the center both the metric potentials are finite. It confirms that our system is free from any sort of singularities, i.e., physical or geometrical singularities. The variations of ρ, p r and p t are shown in the upper, middle and lower panel in Fig. 2, respectively. We find that density and pressure functions are maximum at the surface and decrease monotonically through out the system to reach the minimum value at the surface and confirms regularity of the achieved solutions. We have predicted different values of the central density, ρ c and central pressure p c for the different strange star candidates in Table 1. We find that the densities and radial pressures of the different strange stars are in the order of 10 14 gm/cm 3 and 10 34 dyne/cm 2 , respectively. Due to the strange star candidates as mentioned in Table 1 we find density is much higher than the normal nuclear density ρ normal = 2.3 × 10 14 gm/cm 3 , which confirms that the stars are made of SQM. The variation of the anisotropic stress for the different values of α is shown in Fig. 3 and it confirms the prediction by Deb et al. [45] that for an anisotropic strange star the anisotropic stress should be maximum at the surface. The profile of the electrical charge q(r) and electrical energy density E 2 (r)/8π is featured in the upper and lower panel in Fig. 4, respectively. We find that the total charge Q and the associated electric field E are in the order of 10 20 C and 10 21 − 10 22 V /cm, respectively. Our study clearly reveals that the electric charge has a significant effect on the different physical parameters and the stability of the anisotropic spherically symmetric system. Both Fig. 2 and Table 2 shows that as the charge increases the density and pressures of the stellar system decreases gradually. Interestingly, Fig. 3 features that the effect of anisotropy on the stellar system is maximum when the system is neutral. However, the anisotropic stress of the system decreases consequently with the increasing effect of the charge. We perform different physical tests, viz., energy conditions, mass-radius relation, generalized TOV equation and Herrera cracking concept, etc. In Fig. 5 we have shown that our system is consistent with all the energy conditions. We have featured variation of M (normalized in M ⊙ ) with respect to R for the different values of α in Fig. 6. The solid circles in Fig. 6 denotes the maximum mass points due to different values of α. We found as the charge increases both M max and R Mmax increases gradually. For α = 0.0015 the values of M max and R Mmax increases 4.1% and 0.86%, respectively, than the uncharged case. In the upper and lower panel in Fig. 7 we have presented variation of M and R with respect to ρ c , respectively. For α = 0.0015 km −2 the maximum mass point is achieved for ρ c = 7.024 ρ nuclear , which is 6.07% lower than the value of ρ c as in uncharged case. We have also presented variation of the total charge Q with respect to ρ c due to different values of α in Fig. 11. The figure reveals that as the value of α increases the total charge (Q max ) corresponding to M max is achieved for the lower value of ρ c . The variation of the redshift function with respect to r/R is shown in Fig. 8. To examine stability of the system we have studied Generalized TOV equation which predicts that for our system sum of the forces F g , F a , F e and F h is zero and variation of the forces due to different values of α is shown in Fig. 9. Further, Fig. 10 features that our system is consistent with the Herrera cracking concept by satisfying all the inequalities simultaneously given as 0 ≤ v 2 sr ≤ 1, 0 ≤ v 2 st ≤ 1 and |v 2 st − v 2 sr | ≤ 1. In Table 1 we have predicted a detailed data sheet of the different physical parameters for the different strange star candidates due to α = 0.0010 km −2 and B = 83 M eV /f m 3 . Further, with the motivation to discuss the effects of the increasing value of the electric charge, we have predicted numerical values of the different physical parameters for the strange star candidate LM C X − 4 in Table 2. The high redshift value (0.2824 − 0.2853) supports that the proposed model is suitable to study strange star candidates. Both Tables. 1 and 2 feature that due to different values of α the predicted values of mass to radius ratio for the different strange star candidates are well with in the upper limit of the mass-radius ratio provided by Andréasson [65].
In a summery, in this article we have presented an anisotropic charged spherically symmetric stellar model which is suitable to study ultra-dense strange stars.