New Classes of Charged Spheroidal Models

New classes of exact solutions to the Einstein-Maxwell system is found in closed form by assuming that the hypersurface t = constant is spheroidal. This is achieved by choosing a particular form for the electric field intensity. A class of solution is found for all positive spheroidal parameter K for a specific form of electric field intensity. In general, the condition of pressure isotropy reduces to a difference equation with variable, rational coefficients that can be solved. Consequently, an explicit solution in series form is found. By placing restrictions on the parameters, it is shown that the series terminates and there exist two classes of solutions in terms of elementary functions. These solutions contain the models found previously in the limit of vanishing charge. Solutions found are directly relating the spheroidal parameter and electric field intensity. Masses obtained are consistent with the previously reported experimental and theoretical studies describing strange stars. A physical analysis indicates that these models may be used to describe a charged sphere.


Introduction
In recent years, there have been several investigations into the Einstein-Maxwell system of equations for static spherically symmetric gravitational fields with isotropic pressures in the presence of the electromagnetic field.In such study, regular interior spacetime is matched smoothly at the pressure free interface to the Reissner-Nordstrom exterior model.The models generated are useful to describe charged relativistic bodies with strong gravitational fields such as neutron stars.Gravitational collapse of a spherically symmetric distribution of matter to a point singularity may be avoided in the presence of electromagnetic field.In this situation, the gravitational attraction is counterbalanced by the repulsive Columbian force with the pressure gradient and, hence, charged fluids have a tendency to resist the gravitational collapse.This property persuades the researchers to work on charged perfect fluid distribution.Bonnor [1] has shown that charged dust solutions are expected to form a point like model of electron when its radius shrinks to zero.The presence of electromagnetic field affects the value of redshifts, luminosities, and maximum mass of a compact relativistic object (Ivanov [2], Sharma et al. [3]).Many exact solutions which satisfy the conditions for a physically acceptable charged relativistic sphere have been given by Ivanov [2], Thirukkanesh and Maharaj [4], and Gupta and Maurya [5], among others.Detailed studies of Sharma et al. [6] in cold compact objects, Sharma and Mukherjee [7] analysis of strange matter and binary pulsars and Sharma and Mukherjee [8] analysis of qark-diquark mixtures in equilibrium are of interest physically.Thomas et al. [9], Tikekar and Thomas [10], and Paul and Tikekar [11] have shown that charged relativistic matter is relevant in modeling core-envelope stellar system in which the stellar core is an isotropic fluid surrounded by a layer of anisotropic fluid.
Vaidya and Tikekar [12] proposed the geometry of the spacelike hypersurfaces generated by  = constant are of 3-spheroid to generate exact solutions since it provides a clear geometrical interpretation: the models with spheroidal geometries directly related to the physical situations.Tikekar [13] found an exact solution for a particular spheroidal geometry which could be used to model superdense neutron stars of densities in the range of 10 14 gcm −3 ; this solution has been generalized by Maharaj and Leach [14].There have been extensive studies on charged spheroidal stars by considering a particular form for the electric field in recent years [15][16][17][18][19][20][21][22][23].These charged spheroidal models contain uncharged neutron stars in the relevant limit and are consequently relevant in the description of dense astrophysical objects.Therefore, the study of charged fluid spheres in static spherically symmetric spacetimes is important in relativistic astrophysics.
The objective of this paper is to generate new classes of charged spheroidal solutions in terms of elementary function, which may be used to describe the interior of a relativistic compact sphere.In Section 2, the Einstein-Maxwell system of equations is expressed for static spherically symmetric spacetime.In Section 3, particular forms for one of the gravitational potentials with spheroidal parameter and the electric field intensity are chosen, which reduces the condition of pressure isotropy to a second order linear differential equation in the remaining gravitational potential.In Section 4, a class of solutions for a particular parameter value is first obtained.In general, the solution is obtained in series form using the method of Frobenius, and then two categories of solutions in terms of elementarily functions are derived by placing restrictions on the parameters.The physical features are illustrated graphically, and numerical values of some physical quantities are calculated for a particular example in Section 5.

Field Equations
The gravitational field should be static and spherically symmetric to describe the internal structure of a charged dense compact relativistic sphere.Therefore, the generic form of the line element for describing such configuration is given by in Schwarzschild coordinates (  ) = (, , , ), where ]() and () are arbitrary function of radial coordinate .The Einstein-Maxwell system of field equations, for the metric (1), can be written in the form The energy density  and the pressure  are measured relative to the commoving fluid 4-velocity   =  −]   0 and primes denote differentiation with respect to the radial coordinate .The quantities  and  denote the electric field intensity and the proper charge density, respectively.In the system (2a)-(2d), the units used are such that the coupling constant 8/ 4 = 1 and the speed of light  = 1.This system of equations determines the behaviour of the gravitational field for a charged perfect fluid source.When  = 0 the Einstein-Maxwell system (2a)-(2d) reduces to the uncharged Einstein system.

Choosing Gravitational Potential and Electric Field Intensity
The aim is to seek solutions to the Einstein-Maxwell system (2a)-(2d) by making explicit choices for the gravitational potential  2() and electric field intensity  on physical grounds.The system (2a)-(2d) comprises four equation with six unknowns , ], , , , and  so that it is necessary to choose two of the variables to integrate the system.In this treatment,  and  are specified.A particular choice for  is made such that where  and  are real constants.The above form has been used previously by Tikekar [13] and Maharaj and Leach [14] to study the behaviour of uncharged superdense stars.Note that the choice (3) for the gravitational potential  restricts the geometry of the 3-dimensional hypersurfaces  = constant to be spheroidal for  ̸ = 0 and spherical for  = 0. Eliminating  from (2b) and (2c), for the particular form (3), one obtain the condition of pressure isotropy: The transformation reduces the condition (4) for pressure isotropy to a convenient form in terms of the new variables  and , where dots denote differentiation with respect to .
In terms of new variable , the Einstein-Maxwell system (2a)-(2d) becomes for the choice (3).Note that in (7a)-(7c), , , and  are defined in terms of .Equation ( 6) may be integrable if a particular choice of the electric field intensity  is made.For mathematical convenient one may take where  is a constant which is different from the choice of Komathiraj and Maharaj [15].On substituting ( 8) into ( 6), we obtain a second order linear differential equation in .It is expected that investigation of (9) will produce physically reasonable models of charged stars since  = 0 yields models found previously by Maharaj and Leach [14] for neutron stars which contain Tikekar [13] superdense stars as special case.

Solution
It is clear that the solution of the Einstein-Maxwell system depends on the integrability of (9).One may consider the following two cases.

Particular Case.
When  = −1, (9) becomes If we utilize the transformation equation ( 10) reduces to the equation of free oscillation for  > 0. The solutions of (12) become This becomes in terms of the variable .Thus, a class of exact solution to the Einstein-Maxwell system is generated for all positive value of .Solution ( 14) is given in simple form which is an advantage for physical analysis.

General
Case.Note that ( 9) can be transformed to a hypergeometric equation.However, it is impossible to express the solutions in terms of elementary function for all .
In general, the solution will be given in terms of special functions.Solutions in a simple form are important for a detailed physical analysis.Hence, first I attempt to obtain a general solution of (9) in a series form using the method of Frobenius and then demonstrate the possibility to extract solutions in terms of polynomials and algebraic functions by imposing restrictions on the parameters.

Series Solution.
Since  = 0 is a regular point of the differential equation ( 9), we can apply the method of Frobenius about  = 0 to obtain a series solution.Thus, we assume that is a solution of (9).To express the solution, the coefficients of the series   need to determined explicitly.Substituting (15) into (9) we obtain For the validity of ( 16), we must have It remains to obtain the coefficients   from the system (17a)-(17c).Equation (17c) is the linear difference equation governing the structure of the solution.The difference equation (17c) consists of variable, rational coefficients.First consider the even coefficients  0 ,  2 ,  4 , . . . .These coefficients generate a pattern where the symbol ∏ denotes multiplication.The odd coefficients  1 ,  3 ,  5 , . . .can be written in the form Hence, the difference equation (17c) has been solved and all nonzero coefficients are expressible in terms of the leading coefficients  0 and  1 .Thus, from (15), (18), and (19) the solution of (9) becomes where  0 and  1 are arbitrary constants.Therefore, the general solution of (9) for the choice ( 15) is given by where are linearly independent solutions of (9).
It is interesting to observe that when  = 0 the series solution (21) reduces to a simple form In this case, the electric field intensity  vanishes and there is no charge.Note that when  0 = 0 and  1 = 1, the line element (1) takes the form The above metric corresponds to the familiar isotropic uncharged de Sitter model.
When  1 = 0, the line element (1) takes the particular form The above metric corresponds to the well-known isotropic uncharged Einstein model.
Secondly, consider the polynomials of odd degree.For this case, if we set then for a fixed integer  > 0, (17c) becomes where I have set  1 = 1/.Therefore, from ( 15) and (32), we can express the polynomial in odd powers of  as Polynomial (29) and (33) comprise the first solution of ( 9) for appropriate values of .
Second Solution.Assume the second solution of (9) to be of the form where () is an arbitrary function.Substituting  in (9) we obtain which is a linear differential equation in ().
Observe that ( 9) and ( 35) are of the same type.As in Section 4.2.1, we can first find a general series solution and then two classes of polynomial solution (in even powers of  and in odd powers of ) for (35) using the above technique.Hence, I present the final form of the solution: the polynomial in even powers of  leads to the expression where the real constant  is restricted as integer such that  ≤  − 1; the polynomial in odd powers of  leads to the result for  = 2 − (2 − 1) 2 + 4, where the real constant  is restricted as integer such that  ≤  − 2 in this case.Hence, the solutions to (9) becomes: for  = 2[1 − 2 2 + (2 + 1)], where  is an integer such that  ≤  − 1; for  = 2 − (2 − 1) 2 + 4, where  is an integer such that  ≤  − 2.
The algebraic functions (38) and (39) comprise the second solution of ( 9) for appropriate values of .
Exact Solution.The solutions generated in Section 4.2.2 can be expressed in terms of two classes of elementary functions.The first category of solution for () = () is for the values The second category of solution for () = () is for the values In ( 40) and (42),  and  are arbitrary constants.The solutions (40) and ( 42) are given completely in terms of elementary functions: this has the advantage of facilitating the analysis of physical feature of the stellar interior.These solutions are applicable to a charged superdense star with spheroidal geometry.Note that this treatment has combined both charged and neutral cases for a relativistic star: by setting  = 0 one obtain the solution for neutral case directly.
From these general class of solutions (40) and (42), it is possible to regain particular solutions found in the past.For example, the solutions (40) and (42) reduce to the uncharged models of Maharaj and Leach [14] when  = 0 which contain the Tikekar [13] superdense neutron star model for  = −7.Other explicit functional forms for  are obtainable which could be useful to study dense stars.For example, if we set  = 1 and  = −4 ( = 2) then (42) becomes For this case, the line element takes a simple form in terms of the original variable .

Physical Analysis
It is easy to observe from the solutions (40) and (42) that the spheroidal parameter increases with  so that the electric field intensity affect the spheroidal parameter.For vanishing of pressure at the boundary  =  in the solutions ( 14), (40), and (42), we require () = 0, that gives the condition The generated interior metric matches smoothly with the exterior Reissnar-Nordstrom metric: across the boundary  = , where  and  are the total mass and charge of the sphere.This generate the relationship between , , ,  and the arbitrary constants involved in the solutions  as follows: The conditions (46) and (48a)-(48b) place the restrictions on the parameters involved in these equations.However, there are sufficient free parameters to satisfy the necessary condition that arise for the model under investigation.Now, we shall show that the obtained models are physically reasonable by plotting for a particular solution (45) as illustrated in Figures 1-7 by assuming the parameter values  = −10,  = −1, and  = 16 km in the interval 0 ≤  ≤ 1, where  is in km.The software package Mathematica was utilized to generate the plots for  2] ,  2 , , ,  2 , /, and  versus .Gravitational potentials are plotted in Figures 1 and 2 which are nonsingular at the center and increasing from center to the boundary.The behaviour of the energy density is plotted in Figure 3, which is positive and monotonically decreasing towards the boundary of the stellar object.The behaviour of matter pressure  is plotted in Figure 4, which is regular, monotonically decreasing throughout the stellar interior.The behaviour of the electric field intensity is described in Figure 5, which is continuous and well behaved.The derivative / is plotted in Figure 6, which is 0 ≤ / ≤ 1 throughout the interior of the stellar object.Therefore, the speed of the sound is less than the speed of the light and causality is maintained.The pressure  verses the density  is plotted in Figure 7, and this looks like a straight line so that the equation of state approximate to a linear relationship between  and .Thus, we have demonstrated that the particular solution satisfies the requirements for a physically reasonable stellar interior in the context of general relativity.Numerical values for the physical quantities for the metric (45) for various radius are given in Table 1, where the constant of integration  = 0 is assumed for simplicity.The radius of the star  is chosen to be compared with experimental values reported for realistic stars.In all cases, it is noted that at the surface the pressure () = 0 (up to 8 decimal places), the density () > 0 and the charge () > 0; at the centre the pressure (0) > 0 and the density (0) > 0. The mass of the sphere of radius 7.07 km obtained in Table 1 is comparable with the experimentally determined value 1.44  ⊙ of the strange star SAX J1808.4-3658(SS1)having the same radius [24].Moreover, similar value of mass was theoretically obtained for SAX J1808.4-3658(SS1) by Dey et al. [25], Sharma and Maharaj [26] and Tikekar and Jotania [27].A similar mass value for the SAX J1808.4-3658(SS2) was reported by Tikekar and Jotania [27].Experimental observations suggest a star model for Her.X-1, which is estimated to have a mass between 1.1-1.8 ⊙ and the radius between 6-7.7 km [28].It is noted that the mass obtained for radius 7.7 km and  = 9.445343 km in Table 1 corroborates with this experimentally estimated result for Her.X-1.The mass obtained for 4U 1820-30 is  similar to the value reported in the past [27].However, for this particular example of parameter values considered in Table 1, the casuality condition is not satisfied as for the results obtained by Kumar and Gupta [23] and the argument given [23] for this scenario is notable.

Conclusion
New classes of solutions are generated to the Einstein-Maxwell system by assuming that the hypersurface  = constant is spheroidal.For a particular choice of electric field intensity  ( = −1), a class of solution is for all positive spheroidal parameter values .In general, a series solution is generated to the Einstein-Maxwell system and is demonstrated that two classes of solutions ( 40) and (42) can be extracted in terms of elementary functions.These class of solutions contain de Sitter model, Einstein universe, Tikekar superdense star [13] and Maharaj and Leach neutron star model [14] in the limit of vanishing electric field intensity.For particular parameter values, it is shown geometrically that the model satisfy the necessary physical requirements in the description of a compact object with isotropic matter distribution.Also, the mass values obtained are comparable with experimentally estimated values of realistic stars such as SAX J1808.4-3658(SS1),SAX J1808.4-3658(SS1),Her.X-1, and 4U 1820-30.We believe that the general class of exact solutions found in this paper may assist in more detailed studies of relativistic compact objects.
Energy density in terms of nuclear density   .