A family of charged compact objects with anisotropic pressure

Utilizing an ansatz developed by Maurya and co-workers we present a class of exact solutions of the Einstein-Maxwell field equations describing a spherically symmetric compact object. A detailed physical analysis of these solutions in terms of stability, compactness and regularity indicate that these solutions may be used to model strange star candidates. In particular, we model the strange star candidate Her X-1 and show that our solution conforms to observational data to an excellent degree of accuracy. An interesting and novel phenomenon which arises in this model is the fact that the relative difference between the electromagnetic force and the force due to pressure anisotropy changes sign within the stellar interior. This may be a an additional mechanism required for stability against cracking of the stellar object.


Introduction
General relativity (GR) has proved to be extremely fruitful in describing the physical universe and the structures contained therein. Over and above in contributing to our understanding of astrophysical bodies and the evolution of the Universe, GR is quintessential in understanding the nature of gravity and it's behaviour in the presence of extremely dense sources as well as in higher dimensions [1,2]. With the discovery and observations of ultra-compact objects such as pulsars, neutron stars and black holes, the search for exact solutions of the Einstein field equations has moved away from mere mathematical excursions into the realm of modeling physical objects based on observational data [3]. It's just over a century since the first solution of the Einstein field equations describing a self-gravitating, bounded object was first obtained by Karl Schwarzschild [4]. The Schwarzschild interior solution describes a uniform density sphere and is a first approximation in describing the gravitational field of a static, spherically symmetric object. This solution is a e-mail: sunil@unizwa.edu.om b e-mail: megandhreng@dut.ac.za arXiv:1705.04292v2 [gr-qc] 5 Jul 2017 highly idealised and propagation speeds within the object exceed the speed of light thus rendering the model noncausal. Since the pioneering effort of Schwarzschild in obtaining exact solutions of the Einstein field equations describing the gravitational field of a bounded mass, there has been a massive drive to obtain realistic solutions describing self-gravitating objects. An indepth study of the available exact solutions in the literature indicate that many of them fall short of describing physically realizable stellar structures [5]. Many of these solutions are only valid in some region of the object, other solutions display unphysical behaviour in the density and pressure profiles while many stellar models are unstable against radial perturbations. The wide body of exact solutions which currently exist were obtained through various assumptions on the spacetime geometry, matter content or both [6]. Spherical symmetry is the natural assumption to make when modeling static stars. There is however, more freedom in choosing the matter content of the stellar fluid. In the past researchers have worked with perfect fluids, charged interiors, pressure anisotropy, bulk viscosity and scalar fields. More recently, motivated by developments in cosmology, modelling of stellar structures have included dark energy, dark matter and phantom energy [7,10].
Departure from spherical symmetry has been utilised in modeling stars. The Vaidya-Tikekar (VT) superdense stellar model incorporates a spheroidal geometry for the interior of the star [8]. The VT model has been shown to approximate the behaviour of neutron stars to a very good approximation [9]. Tikekar and co-workers have successfully modeled stars with paraboidal symmetry. The spheroidal parameter which appears in the gravitational potential measures the deviation from spherical symmetry. Work on rotating stars utilise an axis-symmetric metric to describe the stellar interior [11,12]. Recently there has been a surge in obtaining exact solutions of the Einstein field equations via embedding [13,14,15,16,17,18,19,20,21]. In 1947 Karmarkar obtained a restriction which is a necessary condition for embedding a spherically symmetric spacetime in four dimensions into a flat five dimensional spacetime. In general, an n-dimensional Riemannian spacetime is said to be of class p if it can be embedded into a flat space of dimension n + p [22]. The Karmarkar condition relates to class 1 spacetimes. Pandey and Sharma later showed that the Karmarakar condition is only a necessary condition for a spacetime to be of class 1 [23]. A further requirement has to be imposed for sufficiency of the Karmarkar condition. The derivation of the Karmarkar condition is purely geometric in nature which gives a relationship between the two gravitational potentials. This is useful because in order to obtain a complete description of the gravitational behaviour of the model one needs to specify one of the metric functions and the other is obtained via the Karmarakar condition. It is also interesting to note that the Karmarkar condition together with the assumption of pressure isotropy picks out the interior Schwarzschild solution as the only bounded matter configuration with vanishing pressure anisotropy. It follows that if the interior metric of a bounded sphere is of class 1, then the matter content is necessarily anisotropic or charged, with the Schwarzschild interior solution being the only exception.
In this paper we present a model of a spherically symmetric, charged object obtained by embedding a spherically symmetric static metric in Schwarzschild coordinates into a five dimensional flat space. The pressure within the fluid distribution is anisotropic. The resulting condition arising from the embedding reduces the problem of finding an exact solution to the Einstein-Maxwell equations to specifying one of the gravitational potentials and the behaviour of the electric field intensity.

The Einstein-Maxwell Field Equations for Charged Anisotropic Matter Distribution
The line element describing the interior of a static, spherically symmetric matter distribution is given in Schwarzschild coordinates [25,24] x i = (r, θ, φ, t) as follows: where we seek the radial dependence of the potentials λ and ν.
In this work we study charged compact objects within the framework of classical general relativity. The Einstein-Maxwell field equations relating the spacetime geometry to the matter content are where κ = 8π is the Einstein coupling constant. We use geometrized units in which G = 1 = c with G and c being the Newtonian gravitational constant and speed of photons in vacuum, respectively. We assume that the radial and tangential stresses within the interior matter distribution are unequal thus implying that the matter within the star is locally anisotropic. The energy-momentum tensor of fluid the distribution and electromagnetic field are defined respectively as [29] T i where v i is the four-velocity, v i = e ν(r)/2 δ i 4 , θ i is a unit space-like vector in the radial direction, θ i = e λ(r)/2 δ i 1 , ρ is the energy density, pr is the radial pressure and and p t is the tangential pressure. The components for T i j and E i j are defined respectively as: Since we are employing 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 describes the radial component of the electric field. 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 Eq. (7), we obtain For the spherically symmetric metric (1), the Einstein-Maxwell field equations may be expressed as the following system of ordinary differential equations [29] ν r e −λ − λ r where the prime denotes differentiation with respect to r. If we define the anisotropy parameter as ∆ = p t − pr, then from equations (9) and (10) we obtain [30] We note that when ∆ = 0 the pressure is isotropic at each interior point of the matter distribution. The term 2(p t − pr)/r appears in the conservation equations T i j;i = 0 and represents a force due to the anisotropic nature of the fluid. When p t > pr the force associated with anisotropy is directed outwards and when p t < pr (∆ < 0), the force is directed inwards. For ∆ > 0 it has been shown that it is possible to construct more compact objects compared to their neutral counterparts (∆ = 0) [26]. Various approaches have been adopted in finding solutions of the above system of equations. The system (9) -(11) consists of six unknowns ρ, pr, p t , ν, λ and E 2 = q 2 /r 4 , the electric field intensity. In attempting to find exact solutions of the system describing anisotropic charged compact objects one could specify an equation of state of the form p = p(ρ), choose the gravitational potentials ν and λ based on physical grounds or prescribe the behaviour of the anisotropy parameter, ∆.
By using Eqs.(9)-(11), the expression of pressure gradient in terms of anisotropy, metric functions and charge read as dpr dr where the above Eq. (13) represents the charged generalization of the well-known Tolman-Oppenheimer-Volkoff (TOV) equation of hydrostatic for anisotropic stellar structure [24,25]. If the mass function for electrically charged fluid sphere is denoted by m(r) we may then write 2m(r) Various studies of anisotropic charged compact objects have specified the behaviour of the mass function to obtain the gravitational potential λ(r). This is then fed back into the Einstein-Maxwell field equations to obtain the complete gravitational and thermodynamical behavior of the model.

Class one condition for spherical symmetric metric:
Let us consider a 5-dimensional flat line element where we suppose the coordinates z 1 , z 2 , z 3 , z4 and z 5 assume the following forms: In the above transformations K is a positive constant. We can then write the differential forms of the above components as where prime denotes differentiation with respect to the radial coordinate r. Substituting the above expressions dz 1 , dz 2 , dz 3 , dz 4 and dz 5 into the metric (15), we obtain: A direct comparison of metric (21) and metric (1) yields The Eq.(22) provides the embedding class condition. The pressure anisotropy factor, (p t − pr = ∆), is readily obtained from Eqs. (9-10) together with relation (22) as It is clear from (23) that when ∆ = 0, we have p t = pr at each interior point of the fluid distribution. At this point we should highlight the fact that the embedding class condition (22) together with pressure isotropy (∆ = 0) yields only two exact solutions for uncharged fluids, (i) the interior Schwarzschild solution and (ii) the Kohler-Chao solution [27]. The Kohler-Chao solution cannot be used to model a bounded configuration such as a star since there is no surface at which the radial pressure vanishes. Such a surface would define the boundary of the star. The Schwarzschild interior solution describes the interior gravitational field of a uniform density sphere and suffers various pathologies such as the prediction of superluminal propagation velocities within the fluid as well it being unstable against radial perturbations. In a recent paper Maurya and Govender [28] modeled charged compact objects with isotropic pressure via embedding. The pressure isotropy condition becomes a definition for the electric field intensity or the charge distribution. Just as in our approach here the embedding relates the two metric functions ν(r) and λ(r). Recently the modeling of compact objects such as neutron stars, pulsars and strange stars has attracted huge attention amongst researchers. This is mainly due to the fact that a large number of data sets are available in the literature against which the strengths and merits of the various theoretical models can be tested. The role of pressure anisotropy within the stellar core has been highlighted in many of these models.

Generalized charged anisotropic solution for compact star
We can recast Eqs. (9), (10) and (11) in terms of mass function as follows: In this paper we would like to construct a generalized model by adopting a single generic function ν(r). The invariance of the Ricci tensor requires that the energy density ρ(r), radial pressure pr(r) and tangential pressure p t (r) should be finite at the origin. The regularity of Weyl invariants require that mass m(r) and electric charge q(r) should attain minimum values at the centre r = 0 of the configuration (m(0) = q(0) = 0)and attain maximum values at the surface of the star. i.e. m(0) = 0, m (r) > 0 and q(0) = 0, q (r) > 0. In the modeling of charged anisotropic compact stars Maurya et al. [31] have shown that the metric function ν(0) = is finite constant, q(0) = 0, ν (0) = 0 and ν (0) > 0. Since energy density and radial pressure are positive finite and continuous it follows that r > 2m(r) [36,37]. Form pr ≥ 0 with r > 2m(r) we have ν (r) = 0. This shows that the generic function ν(r) is regular minimum at the centre and monotone increasing function of r. Bearing in mind these observations we suppose the generic function ν(r) has the following form: where we have two cases: Here B is positive constant. We observe that ν(0) = ln B, ν = 2 n Ar (1+Ar 2 ) and ν = 2 nA(1−Ar 2 ) (1+Ar 2 ) 2 . It follows that ν(0) > 0, ν (0) = 0, ν (0) = 2 nA > 0 and ν(r) = 0 with r = 0 for both cases (i) and (ii). This implies that this generic source function ν(r) is monotone increasing function of r with regular minimum at r = 0 (Fig.1). Substituting the value of ν into Eq. (22) we obtain This form of the metric function is well-motivated and has been utilised by numerous authors to model compact stars arising from the Karmarkar condition. The parameter n plays a pivot role in the structure and stability of the compact object. Table 5 provides an overview of the class one solutions using the ansatz (27) for the metric function ν(r). We observe that there is a strong connection the range of n which admits physically viable models and nature of the matter content of the star. Anisotropy and electric charge or the absence thereof dictates the admissibility of the range of n. It is evident from Eqn. (27) that in the case of vanishing n the spacetime is rendered flat. In this study we will consider solutions for both n > 0 and n < 0. This approach will allow us to investigate the impact of the 'switch', n on the various thermodynamical properties of the model. The solution is not well behaved in range −3 < n < 2.7. For −7.5 < n ≤ −3 and 2.7 ≤ n < 4. we will get star with low mass. For n ≤ −7.5 and n ≥ 4. we will get stellar models describing compact objects such as Her X-1.  In order to determine the mass function for electrically charged compact star we suppose the electric charge function q(r) is of the form q(r) = E 0 A 2 r 6 (1 + Ar 2 ) n (29) where E 0 is positive constant. We note that the electric field (E = q/r 2 ) vanishes at the center of the configuration. We note that in the case of isotropic pressure the condition of pressure isotropy can be treated as a definition for the charge as a function of the radial coordinate. This approach does not guarantee that q(r) will have physical desirable properties such as the function defined in (29). Utilising Eq. (14) and (29) we readily obtain the mass function m(r) as We observe from Eqs. (29) and (30) that q(0) = 0 and m(0) = 0. However both q (r) and m (r) are positive for r > 0 in both cases (i) and (ii). This indicates that q(r) and m(r) are increasing monotonically away from centre and attains regular minimum at r = 0. By plugging the values of the Eqs (30) and (29) into Eqs. (24), (25) and (26) we get [by assuming f = (1 + Ar 2 ), ψ = Ar 2 ), Relation Eq.(12) provides the anisotropic factor ∆ which is given by where, The gradients of the pr, p t and ρ assume the following forms dpr dr where, Dp r1 = 2 n + E 0 f n−1 [1 + 3 ψ 3 + ψ (5 + 2 C f n ) + 7 ψ 2 + C (2 + n) ψ 2 f n ], The physical viability of our model will be pursued in the next section. < 0. For this purpose we have calculated: Case(i). if n < 0 and A < 0 : Then the Eqs. (38) and (39) provide C > 2n and C < 0. However the Eqs. (40)(41)(42) give 0 ≤ E 0 < n 2 .

Junction conditions
In order to generate a model of a physically realizable bounded object we need to ensure that the interior spacetime M − must match smoothly to the exterior spacetime M + . Since the exterior spacetime is empty, M + is taken to be the Reissner-Nordstrom solution.
The boundary of the star is that surface for which the radial pressure vanishes, pr = 0 at r = R (Misner and Sharp [44]). For our model we obtain The constant B can be determined by using the condition e ν(R) = e −λ(R) , which yields: However the constant A can be determined using the surface density ρs of the star.

Energy conditions
The charged anisotropic fluid sphere should satisfy the following three energy conditions, viz., (i) null energy condition (NEC), (ii) weak energy condition (WEC) and (iii) strong energy condition (SEC). For satisfying the above energy conditions, the following inequalities must be hold simultaneously inside the charged fluid sphere:

Equilibrium condition
The Tolman-Oppenheimer-Volkoff (TOV) equation [24,25] in the presence of charge is given from Eq. (13) as dpr dr The above equation can be expressed into four different components gravitational force Fg = − ν (λ +ν ) 2 r e λ , hydrostatic force F h = −dpr/dr , electric force Fe = q 4 πr 4 dq dr and anisotropic force Fa = 2∆ r which are defined as: Observations of the various panels in Fig. 7 show that the force due to anisotropy dominates the electromagnetic force for small |n|. As |n| increases the difference in the magnitudes of Fa and Fe decrease until they are equal for a particular value of |n|. A further increase in |n| shows that Fe dominates Fa with the relative difference being more marked for large positive values of n.

velocity of sound & stability analysis:
The causality 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. d dr dpr dρ < 0 or d 2 pr dρ 2 > 0 and d dr dpt dρ < 0 or d 2 pt dρ 2 > 0 for 0 ≤ r ≤ r b i.e. the velocity of sound is increasing with the increase of density and it should be decreasing outwards. From Fig. 8 we observe that the speed of sound decreases monotonically from the center of star (high density region) towards the surface of the star (low density region). The sound speed is less than unity thus indicating that causality is preserved within the stellar core (Fig. 9). In their study of the stability of relativistic spheres, Herrera et al. adopted a perturbative scheme in which the energy density and the anisotropy were perturbed and the effect of these perturbations on the fluid elements were studied. They were able to show that different parts of the star respond differently to various degrees of anisotropy which may lead to cracking or overturning within the core. Abreu et al. utilised a different approach to studying cracking in static spheres. In their approach the difference in the tangential and radial sound speeds served as an indicator of potentially unstable regions. They further showed that stable regions within the stellar fluid are characterised by the stability factor, r | which has to be less than unity for a potentially stable configuration. Fig.  10 clearly indicates that our model is stable for a large range of |n|.

Maximum allowable mass and redshift:
The well-known Buchdahl [32] limit for relativistic static spheres 2M/R ≤ 8/9, has been generalised for static charged spheres. Work by Andréasson [33] and Böhmer and Harko [34] showed that mass to radius ratio in the presence of charge was restricted to The compactness u(r) can be defined in terms of the effective mass, m ef f : where and the metric potential e −λ is given in (28). In their study of anisotropic static spheres, Bowers and Liang [35] showed that the surface redshift can be arbitrarily large. In the case of isotropic stars the surface redshift has an upper bound of Zs = 4.77. The relative magnitude of the radial and tangential stresses within the core plays an important role in determining the magnitude of Zs. As pointed out by Maurya et al. [28] when p t > pr the associated surface redshift is greater than its isotropic counterpart. The gravitational surface red-shift (Zs) can be calculated from: We note that the surface redshift depends on the compactness u which should in principle, be constrained the Buchdhal limit. Tables 3. and 4. show that the  surface redshift decreases with an increase in |n| (Fig. 11). For very large |n| the surface redshift is approximately constant.   Table 5 List of embedding class one solutions with well behaved nature of dp i /dρ for the ansatz e ν(r) = B(1 + A r 2 ) n n & A Electric charge Pressure Well behaved Reference function (E) anisotropy (∆) nature of dp i /dρ

Discussion of results
We have presented an exact static model of the Einstein-Maxwell equations which describes a spherically symmetric charged body arising from the requirement that the internal geometry is of embedding class I. The energy momentum tensor de-scribes an anisotropic fluid with an electromagnetic field. Fig. 1 displays the trend in the gravitational potentials as a function of the dimensionless ratio, r/R. The gravitational potentials are continuous and increase smoothly from the center of the star towards the surface. An increase/decrease in n has no appreciable effect in the magnitude or nature of the gravitational potentials. In Fig. 2 we present the trend in the charge (left panel) and the mass (right panel). Both the charge and the mass vanish at the centre of the configuration and increases monotonically towards the surface of the star. An increase in |n| is accompanied by an increase in both the charge and the mass. The divergence is greater towards the surface layers of the star. Fig. 3 illustrates the behaviour of the radial pressure (left panel) and the tangential pressure (right panel). It is clear that the radial and tangential pressures are monotonically decreasing functions of the radial coordinate. The radial pressure vanishes at some finite radius which defines the boundary of the star. We note that the tangential pressure is nonvanishing at the stellar surface. We also note that the radial and tangential pressures increase with an increase in |n|. This increase is noticeable closer to the inner layers of the star and is indistinguishable as the surface layers are approached. A very large increase in |n| of the order of 10 3 has very little effect on the relative magnitudes of both the radial and transverse stresses throughout the interior of the star. The trend of the density is profiled in fig. 4. We observe that the density is a monotonically decreasing function attaining a maximum value at the centre of the star. We observe an interesting trend in the density as |n| increases. An increase in |n| is accompanied by a decrease in the density at each interior point of the gravitating body. For very large values of |n| the density profile is approximately the same for each interior point. The decrease in the density is most noticeable closer to the center of the star. The anisotropy parameter is displayed in Fig. 5. As pointed earlier, the anisotropy parameter vanishes at the centre of the star and increases monotonically outwards towards the surface. An increase |n| is accompanied by an increase in ∆ with the relative differences being more marked towards the surface layers of the star. All the energy conditions are satisfied at each interior point of the configuration as displayed in fig. 6. The various forces operating within the stellar interior are plotted in fig. 7. It has been pointed out that a change in |n| effects changes in the force due to anisotropy and the electromagnetic force. For small values of |n| we pointed out that the anisotropic force dominates the electromagnetic force. This trend switches over for large |n|. Our model obeys the causality condition throughout the stellar interior (Fig.8).The stability of our model was studied by looking at the relative sound speeds squared (Fig. 9). The Abreu et al. stability analysis shows that there are no unstable regions within the stellar core indicating that the likelihood of cracking occurring within our model is remote (Fig. 10). Fig. 7 reveals a new phenomenon associated with these models. We note that for the first time that the relative difference between the electromagnetic force and the force due to anisotropy can change sign and this is directly related to an increase in |n|. Figure 11 shows the trend of the redshift inside the star. The details of the embedding class one solutions with well-behaved nature of dpi/dρ for the ansatz (27) ν(r) = B (1 + A r 2 ) n is given by Table 5.