An analysis of the Reissner-Nordstrom metric shows that a spherically symmetric distri- bution of charged dust may avoid the catastrophic gravitational collapse, a seemingly unavoidable feature of Schwarzschild’s  geometry exterior to an electrically neutral fluid sphere of mass bigger than certain critical limit. As in evidence we have Bonnor’s model   of the equilibrium ball of charged dust in contrast to the Oppenheimer- Snyder  continually contracting ball of electrically neutral dust. Though Bonnors model has been found to the unstable to small radial perturbations and also to a change in the total charge content of the system, it initiates a general interest in the study of the implications of Einstein-Maxwell field equations with reference to the general relativistic prediction of gravitational collapse. On the other hand, it is generally considered that a black hole may carry non-zero net charge, apart from its mass and angular momentum. Such an assumption may require the collapse of stellar masses of charged matter. It has been suggested by Shvartsman  that on account of interaction between a star and its surroundings, it is possible that stellar systems carrying electric charge may exist in nature. It is therefore not surprising that in recent years the problem of finding non- singular, physically meaningful solution of Einstein-Maxwell field equations for static ball of charged coherent perfect has received wide attention. The inclusion of charge seems to affect the stability of the system―the stability of Schwarzschild’s uniform density sphere increases by the introduction of net surface charge. It has been shown that the stability is more profound if the same amount of charge be distributed unifor- mally throughout within the sphere.
The search for the exact solutions is of continuous interest to researcher. Buchdahl  proposed a famous bound on the mass radius ratio of relativistic fluid spheres which is an important contribution in order to study the stability of the fluid spheres. Delgaty- Lake  studied all the then existing solutions and established that Adler  , Heintzmann  , etc. do not satisfy all the well behaved conditions and also pointed out that only nine solutions are well behaved; out of which seven in curvature coordinates (Tolman  , Finch and Skea  , Patvardhav and Vaidya  , Mehra  , Kuchowicz  , Matese and Whitman  , Durgapal’s two solutions  ) and only two solutions (Nariai  , Goldman  ) in isotropic coordinates. Ivanov      , Neeraj Pant  , Maurya and Gupta  , Pant et al.   , Pant and Sah   , Tewari, Charan and Chandra  , Sah, Chandra and Charan  studied the existing well behaved solutions of Einstein’s field equations. Some pioneer work in Relativity is given by Herrera et al.  -  , Tewari and Charan     . Nduka   , Whitman and Burch  , Tikekar  , Ivanov  , Ray et al.  , Stettner  , Krori and Barua  , Ray and Das  , Pant and Negi  , Florides  , Dionysiou  , Pant et al.  etc. gave the well behaved solution for charged fluid sphere. Pant et al.  , Pant and Tewari  , Fuloria et al.  gave charge analogue of Heintzmann, Adler, Durgapal’s relativistic exact solution respectively. Gupta and Maurya  gave charge analogue of Durgapal and Fuloria superdense star. Bijalwan and Gupta  , Gupta and Kumar  gave charge analogue of Schwarzschild’s interior solution.
In this paper, we present a new solution of Einstein-Maxwell field equations in sphe- rically symmetric coordinates which are well behaved solutions charge analogous solution of Sah and Chandra  . In our present study the paper consists of nine sections. In Section 2, Einstein’s field equations for charged fluid sphere in canonical coordinates are given. In Section 3, gravitational binding energy of a charged fluid sphere is given. Section 4 consists of boundary conditions for well behaved solutions. New class of solution of Einstein’s field equations for a charged fluid sphere in canonical coordinates is given in Section 5. Section 6 stipulates the properties of this new class of solution of Einstein-Maxwell field equations. In Section 7 the matching conditions of interior metric of the charged fluid with the exterior metric are given. For better illustration of our physically accepted solution, the relevant physical quantities are presented by tables and figures in Section 8. Finally, some concluding remarks have been made in Section 9.
2. Field Equations for a Charged Fluid Sphere in Canonical Coordinates
The Einstein-Maxwell field equations in general relativity are given by
where, , and, are Ricci mixed tensor, scalar curvature, metric tensor and the energy momentum tensor for fluid sphere respectively. is conserved quantity such that
The energy momentum tensor for a charged fluid sphere is defined as
Here is the part of the energy momentum tensor due to matter distribution of the system and, for a perfect fluid distribution, it is given by
where and are the density and isotropic pressure of the fluid element measured locally in its proper reference frame. The density gives total matter energy in proper volume V as
is element’s time-like four-velocity vector such that
is the part of energy momentum tensor due to electromagnetic character of matter within the fluid sphere and is defined by
where the electromagnetic tensor satisfies Mexwells equations
is the 4-current density for a fluid of null charge conductivity and a conserved quantity such that
and is given by
being the charge density of the fluid element and gives the total charge contained in proper volume V as
The total non gravitational energy in a proper volume V is given by
We consider a static spherically symmetric charged perfect fluid distribution. The interior space-time metric for spherically symmetric fluid distribution in canonical coordi- nate is given by
where and are functions of r only.
The electrostatic field is described by the only non-singular components,
In view of (8), (9) and (10) we obtain
where Q stands for the total charge contained within the sphere of radius r and is given as
In view of the metric (14) and energy momentum tensor (3), the field Equation (1) gives
In view of Equation (17) and Equation (18), pressure isotropy gives
The charge conservation Equation (10) is identically satisfied whereas (2) for energy- momentum gives rise to the following surviving equation
which is clearly contained in the field Equations (17) to (20).
We have three equations to determine five unknown functions, , , and. Thus, we find a two degree of arbitrariness in the general relativistic problem of the electrically charged fluid ball, and to obtain a solution we can always choose one of these functions arbitrarily and correlate other to this by certain relation. The volume field surrounding the charged sphere is described by the Reissner-Nordstrom field
where M and E are constants. We observe that whereas in Schwarzschilds field that total energy is constant, in case of Reissner-Nordstrom field the total energy
increases as r increases. At large distances the Reissner-Nordstron field
approximates to Schwarzschilds field. The junction of interior and exterior field over the boundary of the sphere is governed by the junction conditions due to Darmois (1927), Viz. the continuity of the first and second fundamental forms across the boundary, which imply the continuity of and the fluid pressure across the boundary. For the junction of electromagnetic field it is sufficient to consider the continuity of (and not its first derivatives) across the boundary. In view of (16) and (23) the continuity of gives
Thus the constant measures the total charge contained within the ball.
Also we have
The continuity of over the boundary demands
Thus, in general, the constant M can not be identified with Euclidean mass of the sphere as against the case of uncharged sphere in which case. In view of (27) we have we observe what the distribution of total energy within the charged sphere as measured by an electrically neutral test particle close to the boundary of the star is similar to that within a sphere of electrically neutral perfect fluid. Equation (27) can be rewritten as,
Thus three physical quantities contribute to viz. 1) the total matter energy within the ball distributed as if the geometry is Euclidean, 2) the total electromagnetic energy within the ball as if the geometry is Euclidean 3) the total electromagnetic energy distributed over the entire exterior space-time as if the geometry there too is euclidean.
The gravitational redshift of massive spherically symmetric ball is
which gives central and surface gravitational redshifts
3. Gravational Binding Energy of a Charged Sphere
In view of (11) the total non gravitational energy of a charged sphere is given
Also the total energy of the ball measured by an electrically neutral test particle close to the boundary is given by
Clearly the difference is the expression for the gravitational energy of the charged fluid sphere as measured by an observer close to the boundary. The negative of this quantity is the gravitational binding energy of the system. In Newtonian limit we obtain from (24), (32) and (33)
For a particle at large distances from the object total energy approximates to
. The difference between the two energy measurements is the classical
expression for the energy of vacuum electrostatic field surrounding a sphere of charge and radius. In the Reissner-Nordstrom field energy is distributed in the entire space time. As such the set of field Equations (17) to (19) can be solved under two given relations or assumptions. Physically speaking, one of them has to be the equation of state for fluid and another a law for the distribution of charge with the sphere. The non-singular solution due to Naduka   and that due to Whitman and Burch  follow the charge distribution given by = constant..
4. Boundary Conditions for Well Behaved Solution
For well behaved nature of the solution in isotropic coordinates, the following conditions should be satisfied:
1) The solution should be free from geometrical and physical singularities. Metric potentials A and B must be non-zero positive finite for free from geometrical singularities while central pressure, central density, should be positive and finite or and for free from physical singularities.
2) The solution should have maximum positive values of pressure and density at the center and monotonically decreasing towards the surface of fluid object i.e.
i) and such that the pressure gradient, is negative for
ii) and such that the density gradient, is negative for
3) At boundary pressure, must vanish.
4) The pressure, , and density should be positive.
5) Solution should have positive value of pressure-density ratio which must be less
than 1 (weak energy condition) and less than (strong energy condition) throughout
within the fluid object and monotonically decreasing as well (Pant and Negi  ).
6) The casualty condition must be satisfied for this velocity of sound should be less
than that of light throughout the model i.e.. The velocity of sound should
be monotonically decreasing towards the surface and increasing with the increase of
density i.e. or. 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.
7) For realistic matter, the adiabatic index i.e., everywhere within the
8) The red shift at the center and at the boundary should be positive, finite and monotonically decreasing in nature with the increase of r.
Under these conditions, we have to assume the one of the gravitational potential component in such a way that the field Equation (1) can be integrated and solution should be well behaved.
5. New Class of Well Behaved Solution
We present the following general analytic solution of the field Equations (17) to (20).
The isotropic pressures, matter-energy density, charge, charge density and red shift of charged fluid ball are given by
Here a, b, and d are arbitrary constants.
In order to construct a new relativistic model, we assume then we have
Here a, b, and d are arbitrary constants.
The variation in pressure, density, charge, charge density and red shift with radial distance are given as
6. Properties of the Solution
For real values of metric potentials A and B,. Figure 1 shows that the metric potentials and are positive at the center which are slightly and monotonically increasing with r for suitable choice of constants, and. For the positive central value of,
The central value of, and are given as,
Figure 1. Variation of metric potentials with.
It is clear from Equation (53) to Equation (57) that for positive central values of physical quantities, , , and are positive if, and
. In view of Equations (49) and (50), the variation in the pressure and density
with the distance from the center of fluid ball are identically zero at the center.
At the center of fluid ball the second order derivatives of pressure and density with respect to radial distance from the center of fluid ball are
The pressure is maximum at the center if i.e..
The density is maximum at the center for all constants as.
The central equation of state
must satisfies the condition which demands and.
The central value of is given by
The causality condition at the center gives
It is found that and fall monotonically from their maximum positive values at the center up to non negative values at the boundary (Figure 2) and falls mono- tonically from it’s maximum positive values at the center up to zero value at the boundary while charge increases from zero at the center to maximum positive value at the boundary (Figure 3) for different values of the parameters, , satisfying
Figure 2. Variation of energy density, charge density and red shift with.
Figure 3. Variation of pressure, pressure density ratio and charge with.
Figure 4 shows that speed of sound is less than speed of light i.e. and the ratio of speeds of sound and light falls monotonically from center to the
boundary of the fluid ball.
7. Matching Conditions of Boundary
The solution so obtained are to be matched over the pressure free boundary of fluid sphere smoothly with the Reissner-Nordstrom metric:
which requires the continuity of and across the boundary and
Figure 4. Variation of adiabatic velocity of sound with.
where;, Schwarzchild parameter.
8. Tables of Numerical Values of Physical Quantities and Their Graphs
In view of Equations (64) to (66) the values of, , and are 0.701, −0.062c,
0.0823c and −0.0753 respectively and the value of for and surface density. For better illustration of our physically accepted solu- tion, the relevant physical quantities are presented by means of Table 1, Table 2 and Figures 1-4 for these constants.
In order to construct a super dense star model, we prescribe the surface density of the star as and the values of constants, and can be evaluated for different values of for describing the well behaved solutions. Thus a compact star model can be constructed by finding mass and radius corresponding to assumed surface density. The variation in the mass and radius with schwarzschild parameters for our model of compact star is tabulated in Table 3 and the variation of surface charge, Surface Density, Schwarzschild Parameter and with showing different models is shown in Figure 5.
Table 1., , , , and for.
Table 2., , , , and for.
We have given a new solution for spherically symmetric isotropic charged fluid ball. It has been observed that the physical parameters pressure, density, adiabatic speed of sound and redshift are positive at the centre and within the limit of realistic state equation and monotonically decreasing and the causality condition is obeyed through- out the fluid ball. The charge and electric field intensity are zero at the center and monotonincally increasing towards the intervening surface. Thus, the solution is well behaved for all values of Schwarzschild parameter within the charged fluid ball. Our solution is useful to construct the models of compact star like Strange star family,
Table 3. Variation in, , , , , and with Schwarzschild parameter.
Figure 5. Variation of surface charge, surface density, Schwarzschild parameter and with showing different models.
Neutron star and many more. We have discussed a model of massive neutron star having mass and radius 14.66 km with surface density and central density. The central pressure of neutron star is
while the surface pressure of the star is zero. The electric field intensity at the center is zero and at the surface it comes. Table 3 shows that we can construct different models for neutron star having mass lies between and for different values of ranging from 0.15 to 0.259. The solution reduces to Schwarzschild interior solution for n = −1/2 for electrically neutral fluid ball.
 Schwarzschild, K. (1916) On the Gravitational Field of a Mass Point According to Einstein’s Theory. Sitzungsberichte der koniglich Preussischen Akademie der Wissenschaften Berlin (Mathematical Physics), S42, 189-196.
 Delgaty, M.S.R. and Lake, K. (1998) Physical Acceptability of Isolated, Static, Spherically Symmetric, Perfect Fluid Solutions of Einstein’s Equations. Computer Physics Communications, 115, 395-415.
 Matese, J.J. and Whitman, P.G. (1980) New Method for Extracting Static Equilibrium Configurations in General Relativity. Physical Review D, 22, 1270-1275.
 Ivanov, B.V. (2016) A Different Approach to Anisotropic Spherical Collapse with Shear and Heat Radiation. International Journal of Modern Physics D, 25, Article ID: 1650049.
 Pant, N. (2011) Some New Exact Solutions with Finite Central Parameters and Uniform Radial Motion of Sound. Astrophysics and Space Science, 331, 633-644.
 Pant, N., Fuloria, P. and Tewari, B.C. (2012) A New Well Behaved Exact Solution in General Relativity for Perfect Fluid. Astrophysics and Space Science, 340, 407-412.
 Pant, N., Fuloria, P. and Pradhan, N. (2014) An Exact Solution of Perfect Fluid in Isotropic Coordinates, Compatible with Relativistic Modeling of Star. International Journal of Theoretical Physics, 53, 993-1002.
 Tewari, B.C., Charan, K. and Chandra, P. (2015) A Well Behaved Exact Solution for Spherically Symmetric Perfect Fluid Ball. Journal of Ramanujan Society of Mathematics and Mathematical Sciences, 4, 77-86.
 Sah, A., Chandra, P. and Charan, K. (2016) Spherically Symmetric Anisotropic Perfect Fluid Ball in General Relativity. South East Asian Journal of Mathematics and Mathematical Science, 12, 81-92.
 Herrera, L., Ospino, J. and Perisco, A.D. (2008) All Static Spherically Symmetric Anisotropic Solutions of Einstein’s Equations. Physical Review D, 77, Article ID: 027502.
 Herrera, L., Ospino, J., Perisco, A.D., Fuenmayor, E. and Triconis, O. (2009) Structure and Evolution of Self-Gravitating Objects and the Orthogonal Splitting of the Riemann Tensor. Physical Review D, 79, Article ID: 064025.
 Herrera, L. Parisco, A.D., Hernandez, P. and Santos, N.O. (2004) Dynamics of Dissipative Gravitational Collapse. Physical Review D, 70, Article ID: 084004.
 Herrera, L., Parisco, A.D., Hernandez, P. and Santos, N.O. (1998) On the Role of Density Inhomogeneity and Local Anisotropy in the Fate of Spherical Collapse. Physics Letters A, 237, 113-118.
 Tewari, B.C., Charan, K. and Rani, J. (2016) Spherical Gravitational Collapse of Anisotropic Radiating Fluid Sphere. International Journal of Astronomy and Astrophysics, 6, 155-165.
 Ray, S., Espindola, A., Malheiro, M., Lemos, J. and Zanchin, V. (2003) Electrically Charged Compact Stars and Formation of Charged Black Hole. Physical Review D, 68, Article ID: 084004.
 Krori, K.D. and Barua, J. (1975) A Singularity Free Solution for a Charged Fluid Sphere in General Relativity. Journal of Physics A: Mathematical and General, 8, 508-511.
 Ray, S. and Das, B. (2004) Tolman-Bayin Type Static Charged Fluid Spheres in General Relativity. Monthly Notices of the Royal Astronomical Society, 349, 1331-1334.
 Pant, N. and Negi, P.S. (2012) Variety of Well Behaved Exact Solutions of Einstein-Maxwell Field Equations: An Application to Strange Quark Stars, Neutron Stars and Pulsars. Astrophysics and Space Science, 338, 163-169.
 Florides, P.S. (1983) The Complete Field of Charged Perfect Fluid Spheres and Other Static Spherically Symmetric Charged Distributions. Journal of Physics A: Mathematical and General, 16, 1419-1433.
 Pant, N., Mehta, R.N. and Pant, M. (2011) Well Behaved Class of Charge Analogue of Heintzmann’s Relativistic Exact Solution. Astrophysics and Space Science, 332, 473-479.
 Fuloria, P., Tewari, B.C. and Joshi, B.C. (2011) Well Behaved Class of Charge Analogue of Durgapal’s Reltivistic Exact Solution. Journal of Modern Physics, 2, 1156-1160.
 Gupta, Y.K. and Kumar, M. (2005) A Superdense Star Model as Charged Analogue of Schwarzschild’s Interior Solution. General Relativity and Gravitation, 37, 575-583.