We have presented a class of charged superdense star models, starting with a static spherically symmetric metric in isotropic coordinates for anisotropic fluid by considering Hajj-Boutros-(1986) type metric potential and a specific choice of electrical intensity E and anisotropy factor which involve charge parameter K and anisotropy parameter . The solution is well behaved for all the values of Schwarzschild compactness parameter u lying in the range , for all values of charge parameter K lying in the range , and for all values of anisotropy parameter lying in the range . With the increase in , the values of K and u decrease. Further, we have constructed a superdense star model with all degree of suitability. The solution so obtained is utilized to construct the models for superdense star like neutron stars and strange quark stars . For and , the maximum mass of neutron star is observed as and radius . Further for strange quark stars and are obtained.

1. Introduction

Since the formulation of Einstein-Maxwell field equations, the relativists have been proposing different models of immensely gravitating astrophysical objects by considering the distinct nature of matter or radiation (energy-momentum tensor) present in them. Einstein-Maxwell field equations with anisotropic matter in isotropic coordinates have more importance over Einstein field equations for perfect fluid in curvature coordinates due to following rationale justifications.(i)The presence of some charge may avert the catastrophic gravitational collapse by counterbalancing the gravitational attraction by the electric repulsion in addition to the pressure gradient.(ii)The inclusion of charge inhibits the growth of space time curvature which has a great role in avoiding singularities (Ivanov [1]; de Felice et al. [2]).(iii)Bonnor [3] pointed out that a dust distribution of arbitrarily large mass and small radius can remain in equilibrium against the pull of gravity by a repulsive force produced by a small amount of charge.(iv)The solutions of Einstein-Maxwell equations are useful to study the cosmic matter.(v)The charge dust models and electromagnetic mass models are providing some clue about the structure of electron (Bijalwan [4]) and Lepton model (Kiess [5]).(vi)Several solutions which do not satisfy some or all the conditions for well-behaved nature can be renewed into well-behaved nature by charging them.(vii)Maharaj-Takisa [6] pointed out that the astrophysical objects have essential characteristics of rotational motion, which is caused by the presence of anisotropic parameter. Therefore, the matter in reality cannot be perfect fluid; Ruderman [7] and Sharma and Maharaj [8] also justified that the anisotropic always prevails in a certain density range ≈1015 g/cm3.(viii)Ivanov [9] pointed out that solutions in isotropic coordinates are more significant than the solutions in curvature coordinates, due to the following reasons: (a) the solutions in isotropic coordinates are simple in terms of algebraic expressions; (b) isotropic coordinates solutions can be used as seed solutions in Quasar modeling or nonstatic solutions.

Thus for realistic model it is desirable to study the insinuation of Einstein-Maxwell field equations with reference to the general relativistic prediction of gravitational collapse with anisotropic matter. For this purpose charged fluid ball models are required. The external field of such ball is to be matched with Reissner-Nordstrom solution. The solutions of Einstein-Maxwell field equations successfully explain the characteristics of massive objects like neutron star, quark star, or other superdense objects. Further, these stars are specified in terms of their masses and densities.(a)A neutron star has surface density  g/cm3 (Astashenok [10]) and mass . However, Astashenok [10] established that models with realistic equation of state of neutron star have upper limit mass and minimal radius close to 9 km.(b)A strange quark star has surface density  g/cm3 (Fatema and Murad [11]; Zdunik [12]) and possible maximum mass . However, Dong et al. [13] established that due to presence of half skyrmions in the dense baryonic matter the stable strange quark star can have upper mass limit up to .

In recent past, a considerable number of exact solutions with well-behaved nature of general relativistic field equation with anisotropic matter have been obtained; Dev and Gleiser [14], Komathiraj and Maharaj [15, 16], Thirukkanesh and Regel [17], Takisa and Maharaj [18, 19], Mak and Harko [20], Mak et al. [21], Ivanov [1], Maurya and Gupta [22, 23], Chaisi and Maharaj [24], and Feroze and Siddiqui [25] deal with curvature coordinates and some of them are charged models. By motivation of Maharaj-Takisa [6] and Ivanov [9], in this paper, we present a new class of well-behaved exact solutions of Einstein-Maxwell field equations in isotropic coordinates for anisotropic fluid assuming a particular form of one of the metric potentials and suitable choice of electric intensity and anisotropy.

2. Field Equations in Isotropic Coordinates

We consider the static and spherically symmetric metric in isotropic coordinates where and are functions of .

The Einstein-Maxwell field equations for a nonempty space-time are where is Ricci tensor, is energy-momentum tensor, is the scalar curvature, is the electromagnetic field tensor, denotes the radial pressure, is the transversal pressure, is the density distribution, is the unit space-like vector in the radial direction, and is the velocity vector, satisfying the relation Since the field is static, Thus we find that for the metric (1) under these conditions and for matter distributions with anisotropic pressure the field equation (2) reduces to the following: where prime () denotes differentiation with respect to . From (5) and (6) we obtain the following differential equation in and : Our task is to explore the solutions of (8) and to obtain the fluid parameters , , and from (5), (6), and (7).

To solve the above equation we consider a seed solution as a particular case of Hajj-Boutros [26], Murad-Pant [27] and the electric intensity of the following form: where is a positive constant defined as charge parameter. The electric intensity is so assumed that the model is physically significant and well behaved; that is, remains regular and positive throughout the sphere. In addition, vanishes at the center of the star and increases towards the boundary.

We also take where “” is the anisotropy factor whose value is zero at the center and increases towards the boundary and “” is a positive constant defined as anisotropy parameter.

3. Conditions for Well-Behaved Solution

For well-behaved nature of the solution in isotropic coordinates, the following conditions should be satisfied (Mak and Harko [20] and Maurya and Gupta [22]).(i)The solution should be free from physical and geometrical singularities, that is, finite and positive values of central pressure, central density, and nonzero positive values of and .(ii)The radial pressure must be vanishing, but the tangential pressure may not vanish at the boundary of the sphere. However, the radial pressure is equal to the tangential pressure at the centre of the fluid sphere.(iii)The density and pressures , should be positive inside the star.(iv) and so that pressure gradient is negative for .(v) and so that pressure gradient is negative for .(vi) and so that density gradient is negative for .Conditions (iv)–(vi) imply that pressure and density should be maximum at the centre and monotonically decreasing towards the surface.(vii)Inside the static configuration the casualty condition should be obeyed; that is, the speed of sound should be less than the speed of light; that is, and . In addition to the above the velocity of sound should be decreasing towards the surface; that is, or and or for ; that is, the velocity of sound is increasing with the increase of density. In this context it is worth mentioning that the equation of state at ultrahigh distribution has the property that the sound speed is decreasing outwards (Canuto and Lodenquai [28]).(viii)A physically reasonable energy-momentum tensor has to obey the conditions and .(ix)The central red shift and the surface red shift should be positive and finite; that is, and and both should be bounded.(x)Electric intensity , such that and is taken to be monotonically increasing.(xi)The anisotropy factor should be zero at the center and increasing towards the surface.

4. A New Class of Solutions

Equation (8) is solved by assuming the seed solution as a particular case of Hajj-Boutros [26], Murad-Pant [27] and the electric intensity and the anisotropy factor in such a manner that the solution can be obtained and physically viable. Thus we have On substituting the above in (8), we get the following Riccati-differential equation in : which yields the following solution: where , , , and are arbitrary constants and where is real for .

The expressions for density, radial pressure, and transversal pressure are given by where

5. Properties of the New Solution

The central values of pressure and density are given by The central values of pressure and density will be nonzero positive definite if the following conditions are satisfied: The pressure to density ratios are given by Subjecting to the condition that the ratio of pressure to density should be positive and less than 1 at the centre, that is, , the following inequality has to be satisfied. All the values of which satisfy (21) will also lead to the condition .

Differentiating (16) with respect to , we get Thus it is found that maxima of occur at the centre; that is, Thus the expression of the right-hand side of (26) is negative for all values of satisfying condition (21), showing thereby that the pressure () is maximum at the centre and monotonically decreasing.

Differentiating (17) with respect to , we get Thus it is found that maxima of occur at the centre; that is, Thus the expression of the right-hand side of (28) is negative for all values of satisfying condition (21), showing thereby that the transversal pressure is maximum at the centre and monotonically decreasing.

Now differentiating equation (15) with respect to we get Thus the maxima of occur at the centre; that is, Thus, the expression of the right-hand side of (30) is negative showing thereby that the density is maximum at the centre and monotonically decreasing.

The square of adiabatic sound speed at the centre, , is given by The causality condition is obeyed at the centre for all values of constants satisfying condition (21). Due to cumbersome expressions the trend of pressure-density ratios and adiabatic sound speeds are studied analytically after applying the boundary conditions.

6. Boundary Conditions in Isotropic Coordinates

For exploring the boundary conditions, we use the principle that the metric coefficients and their first derivatives in interior solution () as well as in exterior solution () are continuous up to and on the boundary . The continuity of metric coefficients of and on the boundary is the known first fundamental form. The continuity of derivatives of metric coefficients of and on the boundary is the known second fundamental form.

The exterior field of a spherically symmetric static charged fluid distribution is described by Reissner-Nordstrom metric given by where is the mass of the ball as determined by the external observer and is the radial coordinate of the exterior region.

Since Reissner-Nordstrom metric (33) is considered as the exterior solution, we will arrive at the following conclusions by matching first and second fundamental forms: Equations (34) to (37) are four conditions, known as boundary conditions in isotropic coordinates. Moreover (35) and (37) are equivalent to zero pressure of the interior solution on the boundary.

Applying the boundary conditions from (34) to (37), we get the values of the arbitrary constants in terms of Reissner-Nordstrom parameter “,” Schwarzschild parameters , and radius of the star ; where we define a new parameter called Reissner-Nordstrom parameter “” given by whose value lies within for .

Surface density is given by

Central red shift is given by

The surface red shift is given by

7. Discussions and Conclusions

From Figures 1 and 2 it has been observed that the physical quantities , , , , , , and are positive at the centre and within the limit of realistic state equation and monotonically decreasing while the quantities are increasing for all values of , , and lying in the ranges , , and , respectively. For the pressure is negative. With the increase in the value of from 0.04 to 0.111 the Schwarzschild parameter “” increases; hence the mass increases, but the value of has to be decreased to 0. With we recover the isotropic model.

By increasing above 0.111 the causality condition is obeyed throughout within the ball, but the trend of adiabatic sound speed (transversal) is erratic. Thus, the solution is well behaved for all values of satisfying the inequality for up to 0.111 and for up to 0.016. From Figure 3 it is clearly shown that the mass of the superdense star has a linear dependence on its radius.

In Tables 1 and 2 we present a model of superdense neutron star and quark star based on the particular solution discussed above. By assuming surface density,  g/cm3 corresponding to , for which , the resulting well-behaved solution has a maximum mass and radius  km (for quark star) and by assuming the surface density  g/cm3 the obtained maximum mass is and radius (for neutron star).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors acknowledge their gratitude to Major General Ashok Ambre, SM, Deputy Commandant, NDA, for his motivation and encouragement. They also extend their gratitude to Professor O. P. Shukla, Principal, NDA, for his encouragement. They are grateful to the anonymous referees for their rigorous review, constructive comments, and useful suggestions.