All spherically symmetric charged anisotropic solutions for compact stars

In the present paper we develop an algorithm for all spherically symmetric anisotropic charged fluid distributions. Considering a new source function ν(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu (r)$$\end{document} we find a set of solutions which is physically well behaved and represents compact stellar models. A detailed study specifically shows that the models actually correspond to strange stars in terms of their mass and radius. In this connection we investigate several physical properties like energy conditions, stability, mass–radius ratio, electric charge content, anisotropic nature and surface redshift through graphical plots and mathematical calculations. All the features from these studies are in excellent agreement with the already available evidence in theory as well as observations.


Introduction
Historically the possibility that self-gravitating stars could actually contain a non-vanishing net charge was first pointed out by Rosseland [1] and later on by several other researchers [2][3][4] with different view points. The general relativistic analog for charged dust stars were discussed by Majumdar [5] and Papapetrou [6]. However, in his pioneering work Bonnor [7] further discussed this issue and also several investigators considered the problem later on in detail in connection to stability and other aspects [3,[8][9][10][11][12][13][14].
Following Treves and Turolla [15], to justify the present work with a charged fluid distribution, Ray and Das [16,17] argue that even though the astrophysical systems are by and large electrically neutral, recent studies do not rule out the existence of massive astrophysical systems that are not electrically neutral. The mechanism is mainly related to acquiring a e-mail: sunil@unizwa.edu.om b e-mail: kumar001947@gmail.com c e-mail: saibal@associates.iucaa.in a net charge by accretion from the surrounding medium or even by a compact star during its collapse from the supernova stage. In this connection it is interesting to note that to study the effect of electric charge in compact stars Ray et al. [18] by assuming an ansatz have shown that in order to see any appreciable effect on the phenomenology of the compact stars, the total electric charge is to be ∼10 20 C.
It has been pointed out by Ivanov [14] that substantial analytical difficulties associated with self-gravitating, static, isotropic fluid spheres when pressure explicitly depends on matter density. However, it is also observed that simplification can be achieved with the introduction of electric charge. One is to note that charged, self-gravitating anisotropic fluid spheres have been investigated by Horvat et al. [19] in studies of gravastars and also recently Thirukkanesh and Maharaj [20] found solutions for the charged anisotropic fluid.
In connection to the stability of the stellar model Stettner [21] argued that a fluid sphere of uniform density with a net surface charge is more stable than without charge. Therefore, as pointed out by Rahaman et al. [22], a general mechanism has been adopted to overcome the singularity due to gravitationally collapsing of a static, spherically symmetric fluid sphere one is to include charge to the neutral system. It is observed that in the presence of charge several features may arise: (i) gravitational attraction is counter balanced by the electrical repulsion in addition to the pressure gradient [23], (ii) it inhibits the growth of space-time curvature which has a great role to avoid singularities [24] and (iii) the presence of the charge function serves as a safety valve, which absorbs much of the fine tuning, necessary in the uncharged case [14].
One can notice that since the breakthrough idea of a white dwarf by Chandrasekhar [25] the study of compact stars was tremendously motivated in the field of ultra-dense objects. In this line of research the other dense compact stars are neutron stars, quark stars, strange stars, boson stars, gravastars and so on. As far as composition is concerned in the compact stars the matter is found to be in a stable ground state where the quarks are confined inside the hadrons. It is argued by several workers [26][27][28][29] that if it is composed of de-confined quarks then also a stable ground state of matter, known as 'strange matter', is achievable, which leads to a 'strange star'. There are two aspects of this assumption behind the strange star: (i) theoretically to explain the exotic phenomena of gamma ray bursts and soft gamma ray repeaters [30,31], and (ii) observational confirmation of SAX J1808.4-3658 as one of the candidates for a strange star by the Rossi X-ray Timing Explorer [32].
It was Ruderman [33] who investigated the idea that nuclear matter may have anisotropic features at least in certain very high density ranges (>10 15 gm/cm 3 ), where the nuclear interaction must be treated relativistically. However, later on Bowers and Liang [34] showed specifically that anisotropy might have non-negligible effects on such parameters like maximum equilibrium mass and surface redshift. We notice that recently an anisotropic matter distribution has been considered by several authors in connection to compact stars [22,[35][36][37][38][39][40][41].
Studies have shown that at the center of the fluid sphere the anisotropy vanishes. However, for a small radial increase the anisotropy parameter increases, and after reaching a maximum in the interior of the star, it becomes a decreasing function of the radial distance [42,43]. So there are several possibilities of expressions for charge functions and pressure anisotropy. It is also indicated by Varela et al. [38] that inward-directed fluid forces caused by pressure anisotropy may allow equilibrium configurations with larger net charges and electric field intensities than those found in studies of charged isotropic fluids.
The algorithm for a perfect fluid and an anisotropic uncharged fluid has already been published by others [44][45][46]. In his work Lake [44,45] has considered an algorithm based on the choice of a single monotone function which generates all regular static spherically symmetric perfect as well as anisotropic fluid solutions of Einstein's equations. On the other hand, Herrera et al. [46] have extended the algorithm to the case of locally anisotropic fluids. Therefore, there remains a natural choice of an algorithm so as to apply to a more general case with the inclusion of charge along with an anisotropic fluid distribution.
Against the above background and with the above motivation, in the present paper, we have carried out an investigation for a relativistic stellar model with charged anisotropic fluid sphere. The schematic format of this study is as follows: we provide the Einstein-Maxwell field equations for a charged anisotropic stellar source in Sect. 2, whereas an allied algorithm has been constructed in Sect. 3. The general solutions are shown in Sect. 4, along with a special example for the index n = 1 and matching of the interior solution with the exterior Reissner-Nordström solution. In Sect. 5 we explore several interesting properties of the physical parameters which include density, pressure, stability, charge, anisotropy and redshift. Special case studies are conducted in Sect. 6 to verify (i) the mass-radius ratio and (ii) the density of the star, both of which clearly indicate that the model represents a stable configuration of a strange compact star. Section 7 is meant as a platform for providing some salient features and concluding remarks.

The field equations for charged and anisotropic matter distribution
In this work we intend to study a static and spherically symmetric matter distribution whose interior metric is given in Schwarzschild coordinates [47,48] x i = (r, θ, φ, t) as follows: The Einstein-Maxwell field equations are as usual given by where κ = 8π is the Einstein constant with G = 1 = c in relativistic geometrized unit, G and c, respectively, being the Newtonian gravitational constant and velocity of photon in vacua.
The matter within the star is assumed to be locally anisotropic fluid in nature and consequently T i j and E i j are the energy-momentum tensor of fluid distribution and electromagnetic field defined by [49] where v i is the four-velocity as v i = e ν(r )/2 δ i 4 , θ i is the unit space like vector in the direction of radial vector as θ i = e λ(r )/2 δ i 1 , ρ is the energy density, p r is the pressure in the direction of θ i (normal or radial pressure) and p t is the pressure orthogonal to θ i (transverse or tangential pressure), while T 1 1 = −p r , T 2 2 = T 3 3 = −p t , T 4 4 = ρ and r 4 . Now, anti-symmetric electromagnetic field tensor F i j can be defined by which satisfies the Maxwell equations where g is the determinant of quantities g i j in Eq. (2) defined by where A j = (φ(r ), 0, 0, 0) is the four-potential and J i is the four-current vector defined by where σ is the charged density. For static matter distribution the only non-zero component of the four-current is J 4 . Because of spherical symmetry, the four-current component is only a function of radial distance, r . The only non-vanishing components of electromagnetic field tensor are F 41 and F 14 , related by F 41 = −F 14 , which describe the radial component of the electric field. From Eqs. (7) and (9), one obtains the following expression for the component of the electric field: where √ g 44 = e ν/2 and if q(r ) represents the total charge contained within the sphere of radius r , then it can be defined by the relativistic Gauss law as From Eqs. (10) and (11), we obtain the electric charge q(r ) as For the spherically symmetric metric (1), the Einstein-Maxwell field equations may be expressed as the following system of ordinary differential equations [49]: where the prime denotes differential with respect to r . If the mass function for electrically charged fluid sphere is denoted by m(r ), then it can be defined by the metric function e λ(r ) as If R represents the radius of the fluid spheres then it can be showed that m is constant m(r = R) = M outside the fluid distribution where M is the gravitational mass. Thus the function m(r ) represents the gravitational mass of the matter contained in a sphere of radius r . The gravitational mass M of the fluid distribution is defined as where μ(R) = κ 2 R 0 ρ r 2 dr is the mass inside the sphere, ξ(R) = κ 2 R 0 σ r q e λ/2 dr is the mass equivalence of the electromagnetic energy of distribution and q(R) is the total charge inside the fluid spheres [50]. Now using Eqs. (17) and (11), we can write the mass m(r ) of the fluid spheres of radius r in terms of energy density and charge function as whereas from Eqs. (13) and (16) we obtain We suppose here that the radial pressure is not equal to the tangential pressure i.e. p r = p t , otherwise if the radial pressure is equal to the transverse pressure i.e. p r = p t , which corresponds to isotropic or perfect fluid distribution. Let the measure of anisotropy = p t − p r and is called the anisotropy factor [51]. The term 2( p t − p r )/r appears in the conservation equations T i j;i = 0 (where a semi-colon denotes the covariant derivative) which is representing a force due to anisotropic nature of the fluid. When p t > p r then direction of force to be outward and inward when p t < p r . However, if p t > p r , then the force allows construction of more compact object for the case of anisotropic fluid than isotropic fluid distribution [52].
By using Eqs. (13)- (16) and also Eqs. (18) and (19) the expression of pressure gradient in terms of mass, charge, energy density and radial pressure read where m ≡ dm dr i.e. variation of mass with radial coordinate r . The above Eq. (20) represents the charged generalization of the well-known Tolman-Oppenheimer-Volkoff (TOV) equation of hydrostatic for anisotropic stellar structure [47,48].

The algorithm for constructing all possible anisotropic charged fluid solutions
The Einstein equations (13), (14) and (15) in terms of mass function reduce to 2m Using Eqs. (21) and (22), we obtain a Riccati equation in the first derivative of ν(r ). However, after the re examination of the differential equation we come across a linear differential equation of first order in m(r ) [53].
The first order linear differential equation of m(r ) in terms of ν(r ), anisotropy = ( p t − p r ) and charge function q(r ) can be provided as follows: where Equation (24) gives the mass m(r ) as follows: where where we have used the symbol ≡ d dr . At this point we would like to construct useful algorithm to generate solutions for any known generic function ν(r ). Now from Eqs. (21) and (23), we get Note that the inequalities in (29) and (30) are to be viewed from reality or energy conditions which will impose the restrictions on ν(r ). At the center of symmetry (r = 0) the regularity of the Ricci invariants requires that energy density ρ(r ), radial pressure p r (r ) and tangential pressure p t (r ) at origin should be finite. The regularity of Weyl invariants requires that mass m(r ) and charge q(r ) at r = 0 should satisfy Now the metric function ν(0) is a finite constant, q(0) = 0 and it follows from (30) that ν (0) = 0 and ν (0) = κ 3 [ρ(0)+ 3 p r (0)] − (q (0)) 2 > 0. Since ρ ≥ 0 and continuous, and also since p r > 0 is finite, it follows that r > 2m(r ) [54,55]. With r > 2m(r ) for r > 0. It also follows from (30) for p r > 0 that ν (r ) = 0. As a result, the source function ν(r ) must be a monotone increasing function with a regular minimum at r = 0.

A class of new solutions for charged anisotropic stellar models
For a class of new anisotropic charged stellar models we consider the following suitable source function in the form of the metric potential: where n, B and C are positive integers. It is suitable in the sense that the source function given by Eq. (31) is monotonic increasing with a regular minimum at r = 0. It is to note that charged and uncharged perfect fluid of this source function with different electric intensity has already been carried out [56,57] where it was proved that the above kind of source function with increasing and non-singular behavior provides physically valid solutions.
In terms of the source function expressed in Eq. (31) we consider the electric charge distribution and anisotropic pressure distribution are in the following forms: where K , N and β are positive constants, a and b are positive real numbers and m is a positive integer. The electric field intensity and anisotropy are vanishing at the center and remains continuous, regular and bounded in the inside of the fluid sphere for certain range of values of the parameters. Also these forms of the electric intensity and anisotropy function allow us to integrate Eq. (26). Thus these choices may be physically reasonable and useful in the study of the gravitational behavior of anisotropic charged stellar models. It is observed that Durgapal and Pandey [58], Ishak et al. [59], Lake [44], Pant [60] and Maurya et al. [61] have proposed solutions via the ansatz (31) with some particular values of n. After that Maurya and Gupta [62,63] showed that the same ansatz for the metric function (1) by taking n is a negative integer, C < 0 and C > 0, 0 < n < 1 and it produces an infinite family of analytic solutions of the self-bound type (see details in the Tables 7 and 8 of Appendix). Recently Maurya and Gupta [64] have also obtained infinite family of anisotropic solutions for the same ansatz. But recently Murad [65] obtained charged stellar model for n = −2 and C < 0, however, neutral solutions of this are irregular in the behavior of d p/dρ (Durgapal and Fuloria [66], Delgaty and Lake [67], Pant [60], Maurya and Gupta [63]). Hence the solution is not suitable for application to a neutron star model because the equations of state for nuclear matter show a regular behavior of d p/dρ [66]. So in the present problem we have started with regular behavior of d p/dρ in the same ansatz by taking the value of n = 1 and 3. Recently Maurya et al. [68] argued that neutral solutions for these cases have the regular behavior of d p/dρ and it may be suitable for application to a neutron star model. By using Eqs. (31), (32) and (33), Eq. (26) gives m(r ) in the following form: where where In the absence of electric field intensity (K = 0) and pressure anisotropy (β = 0), Eqs. (21), (22) and (23) reduce to the equations obtained by Maurya et al. [68]. Corresponding solutions belongs to the solutions of Maurya et al. [62,63,68] for the values of n being all negative integers, all positive fractional values between 0 and 1 and some positive integers (n = 1, 2 and 3) and solutions for particular values of n to the well-known Tolman [47] for n = −1, Wyman [69], Kuchowicz [70], Adler [71], Adams and Cohen [72] all for n = −2, Heintzmann [73] for n = −3, Durgapal [74] for n = −4, −5 and Pant [60] for n = −6, −7 for the ansatz (31). 4.1 An example: physical parameters of charged anisotropic model for n = 1 We calculate mass of the charged anisotropic fluid sphere as where The expressions for energy density, radial pressure and tangential pressure are (by taking x = Cr 2 ) given by where is the measure of anisotropy as defined earlier and also In a similar way one can calculate the mass m(r ) of charged anisotropic model for n = 3 and other permissible cases.

Matching and boundary conditions
The metric or first fundamental form of the boundary surface should be the same whether obtained from the interior or exterior metric, guarantees that for some coordinate system the metric components g i j will be continuous across the surface. The requirements of matching condition for metric (1) that the above system of equations is to be solved subject to the boundary condition that radial pressure p r = 0 at r = R (which is the outer boundary of the fluid sphere). It is clear that m(r = R) = M is a constant and, in fact, the interior metric (1) can be joined smoothly at the surface of spheres (r = R) to an exterior Reissner-Nordström metric whose mass is the same as m(r = R) = M [75]. Thus one can get which requires the continuity of e λ(r ) , e ν(r ) and q across the boundary r = R where M and Q are called the total mass and charge inside the fluid sphere, respectively. The continuity of e λ(r ) and e ν(r ) on the boundary is e −λ(R) = e ν(R) , which gives the constant B in the following form: where X = C R 2 .
On the other hand, the arbitrary constant A will be determined from the boundary conditions by putting radial pressure p r = 0 at r = R for the case n = 1 as follows: Hence the total charge inside the star, central density and surface density can, respectively, be evaluated for the case n = 1 as follows: where

Physical acceptability conditions for anisotropic stellar models
In order to be physically meaningful, the interior solution for static fluid spheres of Einstein's gravitational-field equations must satisfy some general physical requirements. Because the Einstein field equation (2) is highly nonlinear in nature, not many realistic physical solutions are known for the description of static spherically symmetric perfect fluid spheres. Out of 127 solutions only 16 were found to be physically meaningful [67]. The following conditions have been generally recognized to be crucial for anisotropic fluid spheres [76].

Case 1
The solution should be free from physical and geometrical singularities i.e. pressure and energy density at the center  (Table 5) should be finite and metric potentials e −λ(r ) and e ν(r ) should have non-zero positive values in the range 0 ≤ r ≤ R. At origin Eq. (17) provides e −λ(0) = 1 whereas from Eq. (32) we obtain e ν(0) = B. So it is clear that metric potentials are positive and finite at the center (Fig. 1).

Case 2
The density ρ and radial pressure p r and tangential pressure p t should be positive inside the star.

Case 3
The radial pressure p r must be vanishing at the boundary of sphere r = R but the tangential pressure p t may not vanish at the boundary r = R of the fluid sphere and may follow p t > 0 at r = R. However, the radial pressure is equal to the tangential pressure at the center of the fluid sphere.

Case 4
(d p r /dr ) r =0 = 0 and (d 2 p r /dr 2 ) r =0 < 0 so that pressure gradient d p r /dr is negative for 0 ≤ r ≤ R.

Case 5
(d p t /dr ) r =0 = 0 and (d 2 p t /dr 2 ) r =0 < 0 so that pressure gradient d p t /dr is negative for 0 ≤ r ≤ R.

Case 1
Inside the fluid ball the speed of sound should be less than the speed of light i.e. 0 ≤ d p r dρ < 1, 0 ≤ d p t dρ < 1 i.e. both d p r dρ and d p t dρ are lies between 0 and 1 which can be observed from Fig. 4 as well as from Table 1.

Case 2
The velocity of sound monotonically decreasing away from the center and it is increasing with the increase of density i.e. d dr d p r dρ < 0 or d 2 p r dρ 2 > 0 and d dr d p t dρ < 0 or d 2 p t dρ 2 > 0 for 0 ≤ r ≤ R (see Fig. 4). In this context it is worth mentioning that the equation of state at ultra-high distribution has the property that the sound speed is decreasing outwards [77].

Case 3
The ratios of the pressure to density, p r /ρ and p t /ρ (as can easily be obtained from Table 1), should be monotonically decreasing with the increase of r , i.e. d     Table 1 which indicates the ratios via the data of p r , p t and ρ.

Energy conditions
A physically reasonable energy-momentum tensor has to obey the following energy conditions [78]: Now we check whether all the energy conditions are satisfied or not. For this purpose, numerical values of these energy conditions are given in Table 2 and accordingly their behavior are shown in Fig. 6. This figure indicates that in our model all the energy conditions are satisfied throughout the interior region.

Method 1
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, that the radial pressure p r must be a monotonically non-decreasing function of ρ [79].
With the energy-momentum tensor of the form (3), the relativistic first law of thermodynamics may be expressed as where p r is the radial pressure, ρ is the total energy density and ρ m is that part of the mass density which satisfies a continuity equation and is therefore conserved throughout the motion. We let the pressure change with density as From Eq. (54) we have By Eqs. (53) and (55) we have where γ is a parameter called the adiabatic index. Matter obeying these equations is stable to gravitational collapse if the pressure times the surface area increases more rapidly than R −2 . Because the density is proportional to R −3 , the force exerted by the pressure is proportional to R 2−3γ . This force increases more rapidly than the gravitational force when γ > 4/3. The latter condition is, however, necessary but not sufficient to obtain a dynamically stable model [80]. Heintzmann and Hillebrandt [81] also proposed that a neutron star with an anisotropic equation of state is stable for γ > 4/3. Also it is well known that Newton's theory of gravitation has no upper mass limit if the equation of state has an adiabatic index γ > 4/3.  The behavior of adiabatic index (γ ) is shown in Fig. 7. It is clear from the figure that the value of γ is higher than 4/3. So our model is stable.

Method 2
For this case let us write the generalized Tolman-Oppenheimer-Volkoff (TOV) equation in the following form: where M G is the gravitational mass within the radius r and is given by Substituting the value of M G (r ) in the above equation we get The above TOV equation describes the equilibrium condition for a charged anisotropic fluid subject to gravitational (F g ), hydrostatic (F h ), electric (F e ) and anisotropic stress (F a ) so that where Now, the above forces can be expressed in the explicit forms as follows: with x = Cr 2 as mentioned earlier.
We have shown the plot for TOV equation in Fig. 8. From the figure it is observed that the system is in static equilibrium under four different forces, e.g. gravitational, hydrostatic, electric and anisotropic to attain overall equilibrium. However, a strong gravitational force is counter balanced jointly by hydrostatic and anisotropic forces. The electric force, it seems, has a negligible effect on this balancing mechanism.

Method 3
In our anisotropic model, to verify stability we plot the radial (V 2 sr = d p r /dt) and transverse (V 2 st = d p t /dρ) sound speeds in Fig. 9. It is observed that these parameters satisfy the inequalities 0 ≤ V 2 sr ≤ 1 and 0 ≤ V 2 st ≤ 1 everywhere within the stellar object which obeys the anisotropic fluid models [82,83]. Again, to check whether local anisotropic matter distribution is stable or not, we use the proposal of Herrera [82], known as cracking (or overturning), which states that the potentially stable region is where the radial speed of sound is greater than the transverse speed of sound. From the left panel of Fig. 10, we can easily say that V 2 st − V 2 sr ≤ 1. Since 0 ≤ V 2 sr ≤ 1 and 0 ≤ V 2 st ≤ 1, we have |V 2 st − V 2 sr | ≤ 1 as can be seen from the right panel of Fig. 10. Hence, we can conclude that our compact star model provides a stable configuration.

Electric charge
From the present model it is observed that in units of Coulomb, the charge on the boundary is 1.5151 × 10 13 C and at the center it is as usual zero. In Table 3 we have put the data for charge q in the relativistic units km. However, to convert these values to Coulomb one has to multiply every value by a factor 1.1659 × 10 20 . The graphical plot is shown in Fig. 11 where the charge profile is such that starting from a minimum it acquires a maximum value at the boundary.
Let us now justify this feature of the charge from the available literature. It was shown by Varela et al. [38] that spheres with vanishing net charge contain fluid elements with unbounded proper charge density located at the fluid-vacuum interface and net charges can be huge (10 19 C). On the other hand, Ray et al. [18] have analyzed the effect of charge in Fig. 11 Distribution of charge with respect to fractional radial distance (r/R) for R X J 1856-37. For plotting this figure we have employed the same data set of numerical values as used in Figs. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 compact stars considering the limit of the maximum amount of charge they can hold and they have shown through numerical calculation that the global balance of the forces allows a huge charge (10 20 C) to be present in a neutron star. Thus we see that the net amount of charge has less effect to balance the mechanism of the force in our model.

Pressure anisotropy
For the present model we calculate the measure of pressure anisotropy as follows: It is in general argued that the 'anisotropy' will be directed outward for the condition p t > p r i.e.
> 0, and inward for the condition p t < p r i.e. < 0. This special feature can be observed from Fig. 12 related to our model. This kind of repulsive 'anisotropic' force allows for construction of a more massive compact stellar configuration [84].
One can also calculate variation of the radial and transverse pressures which are, respectively, given by d p r dr , as can be obtained from Eq. (20), and d p t dr = d p r dr + 2C 2 βr 8π .

Surface redshift
The effective gravitational mass in terms of the energy density can be written as where e −λ(R) is given by Eq. (46). One can therefore provide the compactness of the star as Again we define the surface redshift corresponding to the above compactness factor as follows: We plot the redshift in Fig. 13 from which it is evident that it shows a gradual decrease. This feature also can be observed from Table 3. The maximum surface redshift for the present stellar configuration of radius 6.0 km turns out to be Z = 0.3826. In this connection it is worth to mention that for the isotropic case and in the absence of the cosmological constant the surface redshift is constrained: Z ≤ 2 [85][86][87]. Again for an anisotropic star in the presence of a cosmological constant the constraint on the surface redshift is Z ≤ 5 [88], whereas Ivanov [14] put the bound Z ≤ 5.211. Based on the above discussion we therefore conclude that for an anisotropic star without cosmological constant the value for our model Z = 0.3826 is in good agreement.
6 Some case studies: comparison of present stellar model with compact stars

Allowable mass-to-radius ratio
Buchdahl [85] has proposed an absolute constraint on the maximally allowable mass-to-radius ratio (M/R) for isotropic fluid spheres of the form 2M/R ≤ 8/9 (in units c = G = 1) which states that for a given radius a static isotropic fluid sphere cannot be arbitrarily massive. Böhmer and Harko [87] 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 [89] and it was proved that By substituting the following data: mass M = 0.9693M and radius R = 6.0 km, we find that M/R = 0.238 < 4/9 and also 2M/R = 0.4760, which satisfy the Buchdahl condition of a stable configuration [85]. We also note from Fig. 14 and Table 4 that the charged stars have large mass and radius, as we should expect due to the effect of the repulsive Coulomb force with the M/R ratio increasing with charge [18]. However, unlike Ray et al.'results [18] where in the limit of the maximum charge the mass goes up to 10, which is much higher than the maximum mass allowed for a neutral compact star, our model seems very satisfactory.

Validity with strange star candidates
We have presented two tables here (Tables 5, 6) from which it can be observed that the mass and radius exactly correspond to the strange stars R X J 1856-37 and Her X-1. What we did in the tables is as follows: by considering the mass and radius of the above mentioned stars we have figured out data for the model parameters, and in the next step we evaluated the data for the different physical parameters, e.g. central density, surface density and central pressure, of those strange stars. One can observe that these data are in good agreement with the available observational data.
In this connection we would like to mention that previously Gupta and Maurya [90] showed a similar result for P S R J 1614-2230 with an isotropic fluid distribution and the charge generalization of Durgapal [66]. We also note that like the models presented by Kalam et al. [41], Hossein et al. [84] and Kalam et al. [91] our models provide significantly promising results as regards observational evidence.

Conclusion
In this work we have presented a set of new solutions for an anisotropic charged fluid distribution in the framework of General Theory of Relativity. To solve the Einstein-Maxwell field equations we construct a general algorithm for all possi-  ble anisotropic charged fluid spheres. As an additional condition which simplifies the physical system of space-time we consider a special source function in terms of the metric potential ν. We further adopt the exterior solution of Reissner-Nordström, so that our interior solution can be matched smoothly as a consequence of the junction conditions at the surface of spheres (r = R).
The solutions thus obtained exhibit regular physical behavior as can be observed from the figures and tables for different parameters. We specifically discuss (i) regularity and reality conditions (applied for metric potentials e −λ(r ) and e ν(r ) , energy density ρ, fluid pressures p r and p t , pressure gradients d p r /dr and d p t /dr , and density gradient dρ/dr ), and (ii) causality and well-behaved conditions (applied for speed of sound d p r /dρ and ratios of pressure to densities p r /ρ and p t /ρ). Beside all these general physical properties the solution set shows the desirable and essential features for energy condition, stability condition, charge distribution, pressure anisotropy and surface redshift. Among these physical parameters as a special case, regarding the electric charge distribution of our model, we note that the charge on the boundary is 1.5151 × 10 13 C and at the center it is as usual zero. Other features of charge are also available in the literature [18,38,65] in connection to a stable configuration of compact stars where it has been shown that the global balance of the forces allows a huge charge (∼10 20 C) to be present in a neutron star.
We also observe some special and interesting features for our stellar models, which are related to compact stars as follows.
1. Allowable mass-to-radius ratio: the condition of Buchdahl [85] related to the maximally allowable mass-toradius ratio for isotropic fluid spheres is of the form 2M/R ≤ 8/9. By substituting the following data: mass M = 0.9693 M and radius R = 6.0 Km, we find that 2M/R = 0.4760, which satisfies the Buchdahl condition of a stable configuration [85] as mentioned above. 2. Validity with strange stars: we have prepared several data sets from which it is observed that the mass and radius exactly correspond to the strange stars R X J 1856-37 and Her X-1. Therefore, one can note that like the mod-els of Kalam et al. [41], Hossein et al. [84] and Kalam et al. [91] our models also provide significantly promising results as regards observational evidence.
In this work we have studied the case for n = 1 only in the source function because of the fact that this value is more relevant for exploring existence and properties of strange stars. There is, however, scope for further study with other values of n also as follows: (1) 1 and p = 2). However, the specific value 3/2 is not allowed for these factors to study the solutions for the present model.
As a final comment we would like to mention that Tiwari and Ray [92] proved that any relativistic solution for spherically symmetric charged fluid sphere has an electromagnetic origin and hence it provides an electromagnetic mass model [93][94][95][96]. Therefore, it would be an interesting task to verify whether our model also represents an electromagnetic mass or not and this is to be studied elsewhere in a future project.