During the final stages of stellar evolution when the nuclear fuel inside the star is almost exhausted, the star deviates from its hydrostatic equilibrium stage and starts to collapse. During contraction, gravitational force dominates all other existing forces i.e. pressure gradient forces and radiation forces. The collapse will continue until some quantum-mechanical phenomena such as electron degeneracy pressure, neutron degeneracy pressure generates inside the star thereby providing stability to it against gravitational collapse. Compact objects such as white dwarf stars, neutron stars, strange stars and black holes represent the final stages of a star’s evolution. The stars at this stage become so dense that relativistic effects come into picture and have significant role in describing the crucial features of compact stars. The possibility of existence of Quark star and strange star at such high densities can not be ignored. Neutron star is composed of neutrons, while strange star is composed of u, d, and s quarks. The formation of strange matter can be categorized in two ways: the quark hadron phase transition in the early universe and conversion of neutron stars into strange stars  . When a massive star undergoes a supernova explosion, neutron star and strange quark star are expected to be formed at the inner core of the star. A phase transition between hadronic and strange quark matter may occur at such a higher densities  . Recent developments in observation cosmology have explored interesting facts regarding evolution of these compact objects and revealed many of their features from the study of their emission spectra. But still some of the parameters like mass, internal composition, radii, etc. are not clearly known which cannot be inferred from direct observational data. However many theoretical investigations have successfully explained the various characteristics of compact stellar objects. The study of relativistic models of compact stars like neutron stars and strange stars have been a field of active research in recent years. To understand the nature and complex composition of such immense gravity objects, we require an exact solution to the Einstein’s field equations. The first exact solution of the Einstein’s field equations was obtained by K Schwarzschild in 1916 describing the exterior of a spherically symmetric matter distribution  . Later on the interior solution for a uniform density sphere was also obtained by Schwarzshild  . Of course no astrophysical configuration has a perfect fluid distribution. Observations of stars and the understanding of Particle Physics within compact stars have emphasized the search for more realistic solutions of the field equations. It is interesting to note that by relaxing the condition of a perfect fluid and allowing for pressure anisotropy and charge within the interior of stellar configurations, more realistic structure of compact stellar objects can be realized. It has been suggested that in the presence of charge the gravitational collapse of a spherically symmetric astrophysical object to a point singularity may be averted. The presence of charge in a charged fluid has tendency to oppose the gravitational collapse. This property persuaded the research workers to work on the charged perfect fluid configurations.
Ivanov  and Bonner  investigated a model for charged perfect fluid sphere in which formation of singularity is avoided. Many static fluid spheres which do not satisfy some or all the physical requirements become relevant after the inclusion of charge in them   . The central density of most of compact objects could be several times higher than the normal nuclear matter density. The theoretical investigations      on compact stellar models strongly suggest that the matter distribution in such massive stellar objects may be locally anisotropic, the radial pressure may not be equal to the tangential pressure, at least in certain very high density ranges (>1015 gm∙cm−3). During the post main-sequence of a star, the surface magnetic field reaches up to 1012 - 1013G when neutron star forms. Due to very intense magnetic field of a newly born NS, the matter inside the star may generate pressure anisotropy  . The observational evidence for the neutron star is obtained from the pulsars and its glitches. The sudden change in the time period of pulsars, called glitches is due to neutron superfluidity. Neutron superfluidity also affects the neutrinos emission rate through modified Urca process  . At a density of the order of 1015 gm∙cm−3 nuclear matter may be anisotropic and it needs to be treated relativistically  . Poin condensation at nuclear densities may also cause anisotropy  . Sokolov  suggested that phase transition from normal states of pions to superconducting states can generate pressure anisotropy. Thus it is desirable to study the implications of Einstein-Maxwell field equations with the incorporation of anisotropy for the proper understanding of compact stars. Dev and Gleiser  have shown that pressure anisotropy affects the physical properties and structure of stellar matter and also provides more stability to stellar objects against radial adiabatic perturbations as compared to isotropic matter. Bowers and Liang  have studied anisotropic spheres with uniform energy density and suggested that anisotropy could also play an important role in describing the high redshift objects like quasars. Anisotropy in pressure could significantly affect the physical parameters like maximum compactness, mass and radius of star. A large number of anisotropic models are available in the literature  -  . Local anisotropy in self-gravitating systems were studied by Herrera & Santos  They conjectured that an anisotropic model can be stable. Charged anisotropic models of stellar objects have been discussed by many Astrophysicists i.e. Thirukkanesh and Ragel  , Maurya and Gupta  , Pandya et al  , Bhar et al.  and Murad  . Compact stellar models for neutral/charged static spherically symmetric configuration with anisotropic pressure have been obtained in numerous works. Some of them include; Herrera et al.   , Chaisi and Maharaj  , Maurya and Gupta   , Maharaj et al.  . Pant et al.  have also studied a class of relativistic anisotropic charged stellar models in isotropic coordinates. There have been several recent investigations of static fluid configuration incorporating charge and pressure anisotropy: Pradhan & Pant  , Newton Singh et al.  , Pant et al.  and Maharaj et al.  . Being motivated by the aforementioned recent works we develop some new analytical relativistic anisotropic charged stellar models by using Durgapal’s metric potential  . Our analysis depends on several mathematical key assumptions. The form of metric potential ensures that the metric function is nonsingular, continuous, and well behaved in the interior of the star. This is one of the desirable features for the model on physical grounds. Further, we assume some particular forms of electric charge distribution function and pressure anisotropy so that we may construct a physically reasonable models of stellar objects.
Our whole investigation is divided into 7 sections. In first we give a brief introduction on which our present study is based on. The next section contains the complete procedure to solve Einstein Maxwell’s field equations for anisotropic fluid charge distribution. In a third section we give the elementary criteria to be satisfied by the interior solution in order to present a realistic model. In the last section we conclude all results obtained after complete analysis of newly obtained solution.
2. Einstein Maxwell Field Equations of Anisotropic Charged Fluid Distribution
Let us consider a spherical symmetric metric in curvature coordinates as
For a charged and anisotropic fluid distribution the Einstein Maxwell field equation is given as
Here is the energy momentum tensor and is given by
Here is Ricci tensor, is the energy momentum tensor, R is the scalar curvature, is electromagnetic field tensor, are the radial and transverse pressures respectively, is the energy density, is the velocity vector and is the unit space vector in radial direction.
In view of metric (2.1), the field Equation (2.2) gives
Here is the proper charge density and q is the charge enclosed by a sphere of radius r.
Here B is the constant which has positive value. In view of Equation (2.8), Equation (2.4) and Equation (2.5) reduces to
Here is anisotropy factor.
3. Physical Acceptability Conditions
For well behaved nature of the solution in curvature coordinates, the following conditions should be satisfied:
1) The solution should be free from physical and geometric singularities, i.e. it should yield finite and positive values of the central pressure, central density and nonzero positive value of and .
2) The solution should have positive and monotonically decreasing expressions for pressure and density (p and ) with the increase of r. The ratio of pressure and density should be positive and less than 1(weak energy condition) and less than 1/3 (strong energy condition) throughout within the star.
3) The casualty condition should be obeyed i.e. velocity of sound should be less than that of light throughout the model. In addition to the above the velocity
of sound should be decreasing towards the surface i.e. and or i.e. the velocity of sound is
increasing with the increase of density and it should be decreasing outwards.
4) The adiabadtic index for realistic matter should be .
5) The red shift z should be positive, finite and monotonically decreasing in nature with the increase of r.
6) Electric field intensity E is positive and monotonically increasing from center to boundary and at the centre the electric intensity is zero.
7) The anisotropy factor should be zero at the center and increasing towards the surface.
8) For a stable anisotropic compact star, must be satisfied  .
4. A New Class of Solution
To solve Equation (2.9) we assume
where and are real and positive constants. We assume anisotropy and electric field intensity in such a way that Equation (2.9) is integrable and provides well behaved solutions. The solutions are physically acceptable when both anisotropy and electric field intensity increase as we move from the centre to the surface. On plugging Equation (4.1) into Equation (2.9) we get the following solution,
A is an arbitrary constant of integration.
The expressions for density and pressure are given by
The proper charge density will be given by the following expression
Differentiating Equations (4.3), (4.4) and (4.5) with respect to r, we get the radial and tangential pressure gradients as well as density gradient.
We have shown the variation of metric potentials with radial coordinate for the star PSR J0348 + 0432 in Figure 1. It is clear from figure that metric potentials are regular and free form any singularity inside the fluid sphere. It can be observed from Figure 2 and Figure 3 that radial, transverse pressures and density monotonically decrease with the increase of radius. Electric field intensity increases when we move from centre towards surface (see Figure 4). Moreover the electric field intensity is zero at the centre, which shows that our solution satisfies the required acceptability conditions. From the graphs plotted in Figures 5-8 we notice that pressure to energy density ratio, red shift, proper charge density, square of radial and transverse velocity decrease with the increase of radial coordinate respectively. In Figure 9 we have shown the variation of anisotropy factor with radial coordinate r and we observe that anisotropy factor
Figure 1. Variation of metric potentials with radial coordinate for PSR J0348 + 0432.
Figure 2. Variation of radial and transverse pressures with radial coordinate for PSR J0348 + 0432.
Figure 3. Variation of matter density with radius for PSR J0348 + 0432.
Figure 4. Variation of electric field intensity with radial coordinate for PSR J0348 + 0432.
is zero at the centre and increases in moving towards boundary. Figure 10 shows the outmarch of compression moduli with radial coordinate inside the
Figure 5. Variation of pressure to energy density ratio with radius for PSR J0348 + 0432.
Figure 6. Variation of red shift with radial coordinate for PSR J0348+0432.
Figure 7. Variation of proper charge density with radius for PSR J0348+0432.
fluid sphere and it decreases when we move towards surface. The Null Energy Condition (NEC), Weak Energy condition (WEC) and Strong Energy Condition
Figure 8. Variation of v2/c2 with radius for PSR J0348 + 0432.
Figure 9. Variation of anisotropy with radial coordinate for PSR J0348 + 0432.
Figure 10. Variation of compression moduli with radius for PSR J0348 + 0432.
Figure 11. Variation of energy conditions with radius for PSR J0348 + 0432.
Figure 12. Variation of adiabatic index with radial coordinate for PSR J0348 + 0432.
5. Properties of the New Solution
The values of pressure and density at centre are given by
At the centre of a star the proper charge density must have a finite value
The pressure and density also decrease from the centre to the surface
The square of speed of sound is given by
and must have values less than 1 to fulfill the causality condition and the stability factor must satisfy . The expression for gravitational red shift and adiabatic index are
The gravitational red shift must be non zero and must have positive finite value at the centre. The compression moduli is given by
we can judge the density of compact star with the help of compression moduli, it must have highest value at the centre and should decrease on moving outwards from the centre.
6. Boundary Conditions
The interior solution obtained for super dense star must be matched with Reissner-Nordstrom metric at the boundary, the metric is given by
Here M represents the mass and e the total charge of compact fluid sphere.
Continuity of , and q across the boundary gives the following equations
where and is the charge at . Also the pressure at the boundary
In view of Equation (6.4) the constant A can be computed as
In view of Equations (6.2) and (6.3) the other constants can be evaluated as
Mass of super dense star can be expressed as
The surface density is given by
In Table 1 we have shown the variation of radial and transverse pressures, density, pressure to density ratio, electric field intensity with respect to radial coordinate r and we observe that these parameters show well behaved nature. In Table 2 the profile of proper charge density, red shift, compression moduli, square of the speed of sound, stability factor and anisotropy have been shown. We see that all these variables have physically reasonable behavior. In Table 3 we have mentioned the well behaved values of the parameters used for two compact stars PSR J0348 + 0432 and EXO 0748-676.
Table 1. The variation of radial and transverse pressures, density, pressure to density ratio and electric field intensity with radial coordinate corresponding to k = 2.8, X = 0.1, l = 0.1, m = 0.3, α = 1.6.
Table 2. Variation of proper charge density, red shift, compression moduli, square of the speed of sound, stability factor and anisotropy with radial coordinate with values of k = 2.8, X = 0.1, = 0.1, m = 0.3, α = 1.6.
Table 3. The well behaved values of the parameters used for two compact stars.
7. Results and Discussions
In the present article we have proposed a new model of anisotropic charged compact star which satisfies all the physical reality conditions and is free from central singularity with the exterior space time being the Reissner-Nordstrom metric. We electrified the uncharged fluid sphere e.g. Durgapal solution  with anisotropy taken into account. The charged solution is well behaved with positive values of charge parameter k.
The regularity of metric potentials throughout the interior of the star can be observed from Figure 1 by showing that the solution is free from geometrical or physical singularity. It has been observed from the Figures 2-7 that the physical parameters are positive inside the stellar model and they are monotonically decreasing function of radial coordinate. At the boundary of the star the radial pressure vanishes while the matter density and transverse pressure are non negative. The electric field intensity and adiabatic index also increase with radial coordinate by which we conclude that this is a physically viable configuration Figure 4 and Figure 12. Further the adiabatic index is more than 4/3 throughout the interior of the star, which reconfirms the stability of our model. We also observe that the anisotropic factor and increases when we move from centre towards boundary as observed from Figure 9. Moreover at the centre of a star vanishes which is an essential feature of a realistic star. The profile of velocity of sound has been shown in Figure 8. For our present model both radial and transverse velocity are less than 1 and monotonically decreasing towards the boundary which shows that causality condition holds good inside the star. From Figure 10 it is also observed that the compression moduli of star decreases with the increase of radius which indicates that star is more dense at the centre than outer surface. Furthermore our new model for anisotropic fluid distribution also satisfies the null energy condition , weak energy condition and and the strong energy condition Figure 11. We have shown the variation of required physical parameters with radial coordinate corresponding to , , , , in Table 1 and Table 2. Finally we have constructed a neutron star model by assuming surface density corresponding to, , and with mass and radius 14.51 kms which is close to the observational data of the neutron star PSR J0348 + 0432  . The parameters with values , and generates the model for EXO 0748-676 with radius 13.73 kms and mass  . We have provided the well behaved parameters of these two well known compact stars which gives mass and radius compatible to observed values in Table 3.