Some new Wyman–Leibovitz–Adler type static relativistic charged anisotropic fluid spheres compatible to self-bound stellar modeling

In this work some families of relativistic anisotropic charged fluid spheres have been obtained by solving the Einstein–Maxwell field equations with a preferred form of one of the metric potentials, and suitable forms of electric charge distribution and pressure anisotropy functions. The resulting equation of state (EOS) of the matter distribution has been obtained. Physical analysis shows that the relativistic stellar structure for the matter distribution considered in this work may reasonably model an electrically charged compact star whose energy density associated with the electric fields is on the same order of magnitude as the energy density of fluid matter itself (e.g., electrically charged bare strange stars). Furthermore these models permit a simple method of systematically fixing bounds on the maximum possible mass of cold compact electrically charged self-bound stars. It has been demonstrated, numerically, that the maximum compactness and mass increase in the presence of an electric field and anisotropic pressures. Based on the analytic models developed in this present work, the values of some relevant physical quantities have been calculated by assuming the estimated masses and radii of some well-known potential strange star candidates like PSR J1614-2230, PSR J1903+327, Vela X-1, and 4U 1820-30.


Introduction
The subject of modeling relativistic compact stellar objects through the analytical solutions of Einstein's gravitational field equations has a long history and still the interest remains as one of the key issue to the present researchers. Since the work of Schwarzschild [1], Tolman [2] and Oppenheimer and Volkoff [3] the determination of maximum mass of very compact astrophysical objects has been a key issue in relativisa e-mail: mhmurad@bracu.ac.bd b e-mail: saba@daffodilvarsity.edu.bd tic astrophysics. Such findings are important in astrophysics because analytical solutions enable the distribution of matter in the interior of stellar object, under extreme conditions, to be modeled in terms of simple algebraic relations.
The central energy density of compact stellar object could be of the order of 10 15 g cm −3 , several times higher than the normal nuclear matter density, and due to the absence of reliable information as regards the behavior of matter at such ultra-high density insight into the structure can be obtained by reference to applicable analytic solutions to the equation of relativistic stellar structure [4]. The known analytic solutions of Einstein's gravitational field equations fall into two classes. The first class that describes "normal" matter neutron stars for which the density vanishes at the surface where the pressure vanishes. The Tolman VII solution with vanishing surface energy density falls into this class and hence is a useful approximation to realistic neutron star models. The class that describes stars for which density is finite, about 2-3 times the normal nuclear matter saturation density [5], at the surface where the pressure vanishes includes Tolman IV solution and the solutions discussed in [6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. This type of solutions is a useful approximation to realistic models of "self-bound" strange quark star (strange star for short) [22]. The best-known example of self-bound stars results from the Bodmer-Witten hypothesis also known as the strange quark matter hypothesis asserts that strange quark matter is the ultimate ground state of matter. Still the fundamental significance of this hypothesis remains as a serious possibility in physics and astrophysics [23][24][25][26][27].
An important distinction between quark stars and conventional neutron stars is that the quark stars are self-bound by the strong interaction, gravity just make them massive, whereas neutron stars are bound by gravity. This allows a quark star to rotate faster than would be possible for a neutron star. A quark star can also be bare. The surfaces of a bare strange star and that of normal matter neutron star have striking differences. The very properties of the quark surface, e.g., strong bounding of particles, abrupt density change from 4 × 10 14 g cm −3 to ∼0 in ∼1 fm.
In a very recent past a polytropic quark star model has been suggested [28,29] in order to establish a general framework in which theoretical quark star models could be tested by observations. The key difference between polytropic quark stars and the polytropic model studied previously for normal (i.e., non-quarkian) stars is that the quark star models with nonvanishing density at the stellar surface may not be avoidable due to the strong interaction between quarks which is relevant to the effect of color confinement. As discussed in [28] the polytropic equations of state are stiffer than the conventional realistic models (e.g., the MIT bag model) for quark matter, and pulsar-like stars calculated with a polytropic equation of state could then have high maximum masses >2M . In this framework of a polytropic model a very low massive quark star can also exist and be still gravitationally stable even if the polytropic index n > 3.
Apart from the constituents of these types of compact stars, the most fascinating distinction between a strange star and a normal neutron star is the surface electric fields associated with it. Bare strange stars possess ultra-strong electric fields on their surfaces, which, for ordinary strange matter, is around 10 18 and 10 20 V/cm for color superconducting strange matter [30][31][32]. The influence of energy densities of ultra-high electric fields on the bulk properties of compact stars was explored in [33][34][35][36][37][38]. It also has been shown that electric fields of this magnitude, generated by charge distributions located near the surfaces of strange quark stars, increase the stellar mass by up to 30 % depending on the strength of the electric field. In contrast to the strange star the surface electric field in the case of neutron star is absent [39]. These features may allow one to observationally distinguish quark stars from neutron stars.
In order to obtain a realistic charged stellar model one can start with an explicit EOS and suitable form of electric charge distribution and then integrating the equation of hydrostatic equilibrium, also known as the charged generalization of Tolman-Oppenheimer-Volkoff (TOV) equation [40], which may be obtained by requiring the conservation of mass-energy, as that determines the global structure of electrically charged stars. The integration starts at the center of the star with a prescribed central pressure and ends where the pressure decreases to zero, indicates the surface of the star. Some recent studies include [41][42][43][44]. Such input equations of state do not normally allow for closed-form solutions.
In the second approach one can have insight into such structures by solving the Einstein-Maxwell equations which represent an under-determined system of nonlinear ordinary differential equations of the second order. Due to the high nonlinearity it is difficult to obtain exact solutions to this system. For the special case of a static isotropic perfect fluid, the system of field equations can be reduced to a set of four coupled ordinary differential equations in five unknowns and arrive at exact solutions by making an ad hoc assumption for one of the metric functions or for the energy density. The EOS can then be extracted from the resulting metric. The first exact solutions of field equations, in this approach, known to have astrophysical significance, may have been discovered by Tolman [2]. Out of the different types of exact solutions obtained by Tolman, model V and VI are not considered physically viable, as they correspond to singular solutions (infinite values of the central density and pressure). Except these models, all the other solutions are known as regular solutions (finite and positive pressure and density at the origin). Models IV and VII are found physically viable in the study of compact astrophysical stellar objects. The numerous publications following Tolman's approach include [29,. As might be expected with Tolman's method, unphysical pressure-density configurations are found more frequently than physical ones.
In recent years, however, several authors followed an alternative approach to present analytical stellar models of electrically neutral/charged compact strange stars within the framework of linear equation of state based on MIT bag model together with a particular choice of metric potentials/mass function [74,[78][79][80][81][82][83][84][85]. Some works also studied the viability of nonlinear EOS based on suitable geometry for the description in the interior 3-spaces of such compact star [11,86]. This approach leads to physically viable and easily tractable models of superdense stars in equilibrium. Tikekar and Jotania [87,88], Jotania and Tikekar [89] showed that the ansatz suggested by Tikekar and Thomas [90] has these features and the general three-parameter solution based on it also leads to physically plausible relativistic models of strange stars. Several aspects of physical relevance and the maximum mass of class of compact star models, based on Vaidya-Tikekar ansatz, for the both isotropic and anisotropic pressures have been investigated in [53,[91][92][93].
Out of the 127 known analytical solutions to Einstein's equations, compiled in [94], only a few satisfy elementary tests of physical relevance and, hence, are viable in the description of relativistic compact stellar objects. For strange quark stars, the energy density does not vanish at the surface. Known applicable analytic solutions include [4,5,22]: In contrast, as far the literature is concerned known to the present authors, the charged analogs of Tolman's models (V-VI) obtained in [95][96][97][98][99][100][101] are not physically viable in the description of compact astrophysical objects as regards the infinite values of the central density and pressure. Though the Schwarzschild constant density solution is physically unrealistic, the charged analogs, obtained in [56,[102][103][104], and the charged analog of the Matese and Whitman solution, obtained in [105], may be relevant in the description of self-bound electrically charged strange quark stars. Charged analogs of the Tolman IV and VII models [106][107][108]), as the neutral ones, exhibit the physical features required for the construction of a physically realizable relativistic compact stellar structure. The charged analogs of the Vaidya-Tikekar models have been derived in [109,110]. Astrophysical consequences of the charged analogs of the Vaidya-Tikekar solutions in modeling a electrically charged compact star have been discussed in [59,[111][112][113][114].
It was shown by Bonnor [115,116] that a spherical body can remain in equilibrium under its own gravitation and electric repulsion if the matter present in the sphere carries a certain modest electric charge density. The problem of the stability of a homogeneous distribution of matter containing a net surface charge was considered by Stettner [117]. Stettner showed that a fluid sphere of uniform density with a modest surface charge is more stable than the same system without charge. The electric charge weakens gravity to the extent of turning it into a repulsive field, as happens in the vicinity of a Reissner-Nordström singularity. Thus the gravitational collapse of spherical matter distribution to a point singularity may be avoided if the matter acquires large amounts of electric charge during an accretion process onto a compact object. The gravitational attraction may then be balanced by electrostatic repulsion due to the same electric charge and by the pressure gradient [40,118]. Hence the study of the gravitational behavior of compact charged stellar object has raised the possibility of modeling such compact astrophysical objects in terms of simple algebraic relations between the matter pressure and its energy density.
Of course, no astrophysical object is entirely composed of a perfect fluid. The theoretical investigations [119][120][121][122][123][124][125][126] of more realistic stellar models show that the nuclear matter may be locally anisotropic at least in certain very high density ranges (ρ > 10 15 g cm −3 ), where the nuclear interactions in the stellar matter must be treated relativistically. According to these views, in such massive stellar objects the radial pressure may not be equal to the tangential pressure. Since the pioneering work of Bowers and Liang [127], there has been an extensive literature devoted to the study of anisotropic spherically symmetric static general relativistic configurations (see [75,76,79,82,97, and the references therein).
The principal motivation of this work is to develop some new analytical relativistic stellar models by obtaining closedform solutions of Einstein-Maxwell field equations following the approach of Durgapal [13], and of Maurya and Gupta [107,157]. Our analysis depends on several mathematical key assumptions. First, we choose a particular functional form for one of the metric potentials. The form chosen ensures that the metric function is nonsingular, continuous, and well behaved in the interior of the star. On a physical basis this is one of the desirable features for any well-behaved model. Further, we assume particular forms of electric charge distribution and pressure anisotropy. The maximum allowable mass and corresponding values of physical quantities have been determined. The solutions obtained in this work are expected to provide simplified but mathematically easy to analyze charged stellar models with non-zero super-high surface density, which could reasonably model the stellar core of an electrically charged strange quark star by satisfying applicable physical boundary conditions.
The presentation of this work is as follows. The next section, Sect. 2, is devoted to the solution of the Einstein-Maxwell field equations of an anisotropic fluid and derives the pressure and density relation. In Sect. 3 we present the elementary criteria that have to be satisfied by the obtained solution so as to present a realistic stellar model. Section 4 develops the important ratios by matching the obtained metric components with the space-time exterior to the charged object which is described by the unique Reissner-Nordström metric. A physical analysis is presented of the obtained models in Sect. 5. It is demonstrated numerically that we have a maximum compactness, redshift, and mass increase in the presence of an electric field and anisotropic pressures; this is in agreement with some other work [158]. In Sect. 6 some explicit numerical models of relativistic anisotropic stars, of possible astrophysical relevance, are also presented and we also apply our model to some well-known potential strange star candidates to calculate some physical quantities by assuming the estimated masses and predicted radii. Finally, Sect. 7 discusses and concludes the work.
(2.1.1) 1 Throughout the work we will use c = G = 1, except in the tables.
The functions ν and λ satisfy the Einstein-Maxwell field equations, where κ = 8π is Einstein's constant. The matter within the star is assumed to be locally anisotropic fluid in nature and consequently T μ ν and E μ ν are the energy-momentum tensor of fluid distribution and electromagnetic field defined by [131,159] T μ ν = P t + ρc 2 v μ v ν − P t δ μ ν + (P r − P t )χ μ χ ν , where ρ, P r , P t , v μ , denote the energy density, radial pressure, and tangential pressure of the fluid distribution respectively. v μ and F μν denote the velocity vector and antisymmetric electromagnetic field strength tensor, defined by which satisfies the Maxwell equations, where g is the determinant of quantities g μν in Eq. and A ν = (ϕ(r ), 0, 0, 0) is the four-potential and j μ is the four-current vector, defined by where ρ ch denotes the proper charge density. For a static matter distribution the only non-zero component of the four-current is j 0 . Because of spherical symmetry, the four-current component is only a function of the radial distance, r . The only non-vanishing components of the electromagnetic field tensor are F 01 and F 10 , related by F 01 = −F 10 , which describe the radial component of the electric field. From Eq. (2.1.4a) one obtains the following expression for the electric field: where q(r ) represents the total charge contained within the sphere of radius r defined by Equation (2.1.5) can be treated as the relativistic version of Gauss' law. For the metric (2.1.1), the Einstein-Maxwell field equations may be expressed as the following system of ordinary differential equations [131]: where a prime denotes the r -derivative. In analogy to the electrically uncharged case, one usually introduces a quantity m(r ) by the following expression: (2.1.12) Finally, combining (2.1.11) and (2.1.12), we get which is the charged generalization of Tolman-Oppenheimer-Volkoff (TOV) equation of a hydrostatic equilibrium for the anisotropic stellar configuration [133]. In Eq. (2.1.13) the additional term, 2(P t − P r )/r , represents the "force" which is due to the anisotropic nature of the fluid. This force is directed outward when P t > P r and inward when P t < P r . The existence of a repulsive force (in the case P t > P r ) allows for the construction of a more compact distribution when using an anisotropic fluid than when using an isotropic fluid. Instead of solving Eq. (2.1.13), for any prescribed equation of state, we are rather interested in solving Eqs. (2.1.6)-(2.1.8) with the help of the following ansatz [12,13]: where N is a positive integer and B N , C > 0 are two constants to be determined by the appropriate physical boundary conditions. Subtracting (2.1.6) from (2.1.7) one obtains the equation of the "pressure anisotropy", ) is a second order ordinary nonlinear differential equation in ν but first order linear in λ. An algorithm recently presented by Herrera et al. [149] shows that all static spherically symmetric anisotropic solutions of Einstein's field equations may be generated from Eq. (2.1.15) by two generating functions κ and ν. In our case we have one additional generating function 2q 2 /r 4 . At this moment it is convenient to introduce the following transformations: and A N is the constant of integration, which may be determined by imposing appropriate physical boundary conditions. Once the metric potential Z is obtained, the other physical variables may be expressed in terms of the generating functions and the equation of state may be extracted, parametrically, from the following equations: 2.2 Models of electric charge distribution and pressure anisotropy The "realistic" charge distribution inside the fluid sphere is not known [161], but it seems intuitively reasonable that due to electrical repulsion the charge distribution should be weighted toward the surface [42].  Table 1 (also see [162]). In this work we consider the following forms of the electric charge distribution and the pressure anisotropy: where K , δ ≥ 0, n is a nonnegative integer, and m, p, a are any real numbers. It must be emphasized that these hypothetical models of electric charge distribution and pressure anisotropy are chosen, in terms of x, in such a way that these allow us to integrate Eq. (2.1.18) rather than for any particular physical reasons. Moreover, the electric field intensity and anisotropy vanish at the center and remain continuous and bounded in the interior of the star for a wide range of values of the parameters. Thus these choices may be physically reasonable and useful in the study of the gravitational behavior of anisotropic charged stellar objects. and

Anisotropic charged stellar models
Ishak et al. [175], Lake [176], and recently Maurya and Gupta [107,157] showed that the ansatz for the metric function (2.1.14) produces an infinite family of analytic solutions of the self-bound type. Five of these were previously known (N = 1, 2, 3, 4, and 5). For example, N = 1 corresponds to the Tolman IV model, and N = 2, 3, 4, 5 correspond to Adler [8], 2 Heintzmann [7], 3 and Durgapal models, 4 respectively. The most relevant case is for N = 2, for which the speed of sound ≈ 1/ √ 3 throughout most of the star, similar to the behavior of strange quark matter [180]. The astrophysical significance and the adiabatic stability of the Wyman-Leibovitz-Adler solution (N = 2) in modeling neutron stars was first discussed in [181][182][183].
(a) Regularity conditions: (i) The solution should be free from physical and geometric singularities, i.e., e ν > 0 and e λ > 0 in the range 0 ≤ r ≤ R. (ii) The radial and tangential pressures and density are nonnegative, P r , P t , ρ ≥ 0. (iii) Radial pressure P r should be zero at the boundary r = R, i.e., P r (r = R) = 0, the energy density and tangential pressure may follow ρ(r = R) ≥ 0 and P t (r = R) ≥ 0.
(b) Stability conditions: (iv) In order to have an equilibrium configuration the matter must be stable against the collapse of local regions. This requires Le Chatelier's principle, also known as the local or microscopic stability condition: the radial pressure P r must be a monotonically non-decreasing function of ρ [129], (c) Causality condition: (v) The condition 0 ≤ √ d P r /dρ ≤ 1, 0 ≤ √ d P t /dρ ≤ 1 be the condition that the speed of sound not exceeds that of light.
(e) Monotone decrease of physical parameters: (vii) Pressure and density, should maximum at the center and monotonically decreasing toward the pressure free interface (i.e., boundary of the fluid sphere). Mathematically, dP r dr < 0, dρ dr < 0, 0 < r ≤ R.
(viii) Additionally, the radial velocity of sound may be monotonically decreasing toward the surface. In this context, it is worth mentioning that for different equations of state at ultra-high densities available in the literature [120,121,187] it is found that the speed of sound is decreasing outwards from the center of the fluid sphere, i.e., d dr dP r dρ < 0 for 0 ≤ r ≤ R. (ix) The ratio of pressure to density, P r /ρ and P t /ρ, should be monotonically decreasing with the increase of r , i.e., d dr (f) Matching condition: The interior solution should match continuously with an exterior Reissner-Nordström solution, This requires the continuity of e ν , e λ and q across the boundary r = R, and q(R) = Q, where M and Q represent the total mass and charge inside the fluid sphere, respectively. (g) Charge distribution: (xi) The electric field intensity E, such that E(0) = 0, is taken to be monotonically increasing, i.e., dE/dr > 0 for 0 < r ≤ R. (g) Pressure anisotropy: (xii) The pressure anisotropy vanishes at the center, i.e., (0) = 0 [127,188]. (h) Allowable mass-to-radius ratio: (xiii) Buchdahl [189] obtained an absolute constraint of the maximally allowable mass-to-radius ratio (M/R) for isotropic fluid spheres of the form 2M/R ≤ 8/9 (in the units c = G = 1), which states that for a given radius a static isotropic fluid sphere cannot be arbitrarily massive. Böhmer and Harko [190] proved that for a compact object with charge, Q(<M), there is a lower bound for the mass-radius ratio, The upper bound of the mass of charged sphere was generalized by Andréasson [191] and one proved that

Determination of the arbitrary constant A 2
The boundary condition P r (R) = 0 can be utilized to specify A 2 . For Case I: where X = C R 2 .

Total charge to radius ratio Q/R
Using X = C R 2 in Eq. (2.2.1) we obtain the square of the ratio Q/R, (4.4.1)

Determination of the constant B 2
The constant B 2 can be specified by the boundary condition e ν(R) = e −λ(R) , which gives (4.5.1)

Surface redshift
The surface redshift of the charged fluid sphere is given by

Relativistic adiabatic index and stability
The stability of a fluid sphere, i.e., whether it is Newtonian or relativistic, isotropic or anisotropic, is related to the adiabatic index (the ratio of two specific heats). It is well known that the collapsing condition for a Newtonian isotropic sphere is < 4/3 [195]. For an anisotropic general relativistic sphere the adiabatic index is defined by and the collapsing condition then becomes [196,197] where P r 0 , P t0 , and ρ 0 are the initial radial, tangential, and energy density in static equilibrium satisfying Eq. (2.1.13).
The first and last term inside the square brackets, the anisotropic and relativistic corrections, respectively. It is to be noted that the positive anisotropy, P t − P r > 0, increases the unstable range of [197,198].
To study the stability of anisotropic stars under the radial perturbations, Herrera [199] (also see [200]) introduced the concept of "cracking", breaking of self-gravitating spheres, which results from the appearance of total radial forces of different signs in different regions of the sphere once the equilibrium is perturbed. The occurrence of such a "cracking" may be induced by the local anisotropy of the fluid.
By this concept of cracking Abreu et al. [186] proved that the region of the anisotropic fluid sphere where −1 ≤ v 2 st − v 2 sr ≤ 0 is potentially stable, but the region where 0 < v 2 st − v 2 sr ≤ 1 is potentially unstable. The radial and tangential speeds of sound of the anisotropic sphere developed so far may be obtained from Eqs. (5.1.1)- To satisfy −1 ≤ v 2 st − v 2 sr ≤ 0 throughout the fluid distribution we require d /dρ ≤ 0. As we have dρ/dx < 0, we further require that d /dx ≥ 0, which will be satisfied as long as is an increasing function of x.

For a given radius
(a) Total mass: where the mass M is in the unit of km. 6 6 The following physical constants, with their conventional values, have been used for the numerical calculation: c = 1 = 2.997 × 10 8 m s −1 , G = 1 = 6.674 × 10 −11 N m 2 kg −2 , M = 1.486 km = 2 × 10 30 kg.
(b) Total charge: where the charge Q is in the unit of radius. 7 (c) Central density:

For a given surface density
The total mass, total charge, and the central density can be calculated by Eqs.

For a given central density
The radius of the charged fluid sphere for a prescribed central density can be calculated by the following equation: where the central energy density ρ c is given in kg m −3 and the radius in m. The total mass, total charge, and the surface density then can be calculated by Eqs.

For a given central pressure
The radius of the charged fluid sphere for a prescribed central pressure can be calculated by the following equation: where the central pressure P c is given in N m −2 . 8 The total mass, total charge, and the central and surface densities then can be calculated by Eqs. (5.3.1.1)-(5.3.1.4).

For a given electric charge
The radius of the charged fluid sphere for a prescribed total charge can be calculated by where the charge Q is given in km. Then the total mass, total charge, and the central and surface densities can be calculated by Eqs.

Physical analysis of the models
For each choice of constant parameters (K , m, n, p, δ, a), the maximum mass of the charged star depends on the corresponding set of maximum values of X = X max up to which the pressure and density and their gradients satisfy P r ≥ 0, P t ≥ 0, ρ ≥ 0, dP r /dx < 0, dρ/dx < 0, and the speeds of sound satisfy 0 ≤ √ dP r /dρ ≤ 1, 0 ≤ √ dP t /dρ ≤ 1 and may be monotonically decreasing with increasing x from center to the boundary.
Case Ia: Isotropic pressure The maximum limit of compactness parameter (2M/R) obtained by Adler [8] set by the causality condition v) for the neutral isotropic fluid sphere is about 0.7. But the neutral model fails to satisfy condition viii). To generate an isotropic charged fluid sphere we set (m, n, p, δ, a) = (10 4 , 0, 0.23, 0, 0) and K = 0.08. For this choice the range of values 0 < X ≤ 0.71 is obtained over which the fluid distribution satisfies the elementary criteria mentioned in Sect.
3. With the decrease of K , X increases. The maximum value of compactness parameter is obtained (2M/R) max = 0.8066, using (4.3.1), at K = 0.08, X max = 0.71. Corresponding to (K , X max ) the total charge-to-radius ratio, and the total charge-to-total mass ratio are found to be Q/R = 0.3641 and Q/M = 0.9029 using Eq.  Fig. 5, which also shows that the energy density is positive and monotonically decreasing inside the fluid distributions. The pressure-density profile for the sphere is plotted, parametrically, in Fig. 6. The pressure and energy density gradients are found to be strictly negative throughout the distribution. The adiabatic speed of sound is found to be less than the speed of light, in the unit c = 1, and monotonically decreasing in nature. Case Ib: Anisotropic pressure Footnote 9 continued some 14 orders of magnitude larger than the surface density of normal neutron stars. The mass-radius relation for a sequence of anisotropic charged fluid spheres generated with a particular set of input parameters m = 10 4 , n = 0, p = 0.23, δ = 0.21, a = 0.75, K = 0.228, and 0 < X ≤ 0.65 with surface density ρ s = 4.68 × 10 14 g cm −3 is in Fig. 1. The behavior of Fig. 1 reproduces that of other quark star models [38]. But a  Fig. 3 maximum value of X = X max is found of 0.319 up to which the fluid distribution satisfies the conditions dP t /dx < 0 and P r ≥ 0. Hence, the maximum value of compactness parameter is obtained (2M/R) max = 0.5978, at X max = 0.319. The charge-radius ratio, and charge-mass ratio are also found to be Q/R = 0.2652 and Q/M = 0.8874. The maximum mass is found M max = 2.1274M , Fig. 6 Pressure-density profiles, P r (ρ), for the same fluid spheres as in Fig. 3   Fig. 7 Behavior of pressure gradients for the same fluid sphere as in Fig. 4. The solid (blue) line corresponds to radial pressure gradient κd P r /Cdx and the dash-dotted (red) line corresponds to the tangential pressure gradient κd P t /Cdx with radius R = 10.49 km, central pressure is P c = 1.37 × 10 35 dyne cm −2 , the central energy density is ρ c = 1.40 × 10 15 g cm −3 , and the total charge is Q = 3.23 × 10 20 C, and z s = 0.4547. The behavior of pressure anisotropy is demonstrated in Fig. 2. For this particular choice of constant parameters is found to be maximum at   From the figure it is clear that the gradients remain strictly negative throughout the distribution. The adiabatic speeds of sound v sr and v st are shown in Fig. 9 from which it is found that the speeds are monotonically decreasing in nature. Figure 10 shows that the condition −1 < v 2 st − v 2 sr ≤ 0 is satisfied throughout the fluid configuration and Through a numerical and a graphical analysis we have demonstrated that the models obtained in Sect. 2.3 satisfy the physical requirements for a wide range of values of m, n, p, δ, a, and K , giving us a possibility for different charge variations and anisotropy within the fluid spheres. The resulting spheres can be utilized to construct physically reasonable compact self-bound charged stellar model such as a charged strange quark star.

An application of the model for some well-known strange star candidates
The analysis of very compact astrophysical objects has been a key issue in relativistic astrophysics for the last few decades. Recent observations show that the estimated mass and radius of several compact objects such as X-ray pulsar Her X-1, Xray burster 4U 1820-30, millisecond pulsar SAX J 1808.4-3658, X-ray sources 4U 1728-34, PSR 0943+10, and RX J185635-3754 are not compatible with the standard neutron  star models [202,203]. For a recent review on this the readers are referred to [24]. Based on the analytic model developed so far, to get an estimate of the range of various physical parameters of some potential strange star candidates we have calculated the values of the relevant physical quantities, such as the central/surface pressure and density, by using the refined mass and predicted radius of 12 pulsars recently reported in [204]. The values are reported in Table 4.

Concluding remarks
In this work we have studied some particular simple families of relativistic charged stellar models obtained by solving Einstein-Maxwell field equations for a static spherically symmetric locally anisotropic fluid distribution. We based our work on three ad hoc assumptions: (1) one of the metric potential, (2) the electric charge distribution, (3) and the pressure anisotropy, the analytical equation of state has been com-puted from the resulting metric. These families of analytical relativistic stellar models may be considered as anisotropic charged analogs of the Wyman-Leibovitz-Adler solution.
We have shown that the maximum limit of compactness (2M/R) = 0.7 [8] for the neutral isotropic fluid sphere, by the causality condition, may be changed significantly by the insertion of charge (K = 0) without anisotropy (δ = 0). We have also observed that both the isotropic and the anisotropic neutral models fail but charged models exhibit the monotonically decreasing adiabatic sound speeds.
A wide range of values of constant parameters are allowed to specify the maximum mass of charged fluid spheres. At this point we want to make some remark on the particular choice of stellar surface density ρ s . Various authors usually have chosen ρ s = 2 × 10 14 g cm −3 to calculate the mass and radius of the charged fluid spheres, which have given rise to the stellar configuration as massive as 4-6M with much lower central density. Such a massive configuration may not serve as a realistic model for a self-bound star. This choice is, therefore, not a physical one. Modeling a compact (quark) star requires the use of a higher surface density. Certainly, the value of the surface density affects the calculated value of the stellar mass-to see this, we observe that by the method employed in the present work one can obtain an arbitrarily large maximum mass just by inserting a vanishing small surface density (e.g., 0.1 − 1 g cm −3 to model a thin crust). In our model calculation, the density at the stellar radius was chosen within the range 4-10 × 10 14 g cm −3 [205] and drops abruptly to zero, as with all stellar models matching an interior metric to the external Reissner-Nordström form. This sharp drop in density is a reasonable model approximation, since the thickness of the "quark surface" is of order 1 fm, a negligibly small dimension compared to the stellar radius.
In the construction of fluid spheres we assumed a positive anisotropy, P t > P r , which increases the region of instability, and have shown that the upper bound on the maximum mass decreases in the presence of a positive anisotropy in order to satisfy the condition −1 < v 2 st − v 2 sr ≤ 0 and get the potentially stable sphere. Moreover, for the models the radial speed of sound is obtained ≈ 1/ √ 3 at the center, and it remains almost the same throughout most of the fluid sphere. This behavior is somewhat like the MIT bag model. An analytical stellar model with such physical features could play a significant role in the description of the internal structure of superdense compact astrophysical objects like electrically charged bare strange quark stars. Nevertheless, it would also be interesting to study the behaviors of strange stars within the framework of a linear equation state, like the MIT bag model, in Wyman-Leibovitz-Adler space-time, which the authors hope to do in the near future.