# Anisotropic Charged Fluid Sphere in Isotropic Coordinates

Authors: Neeraj Pant; N. Pradhan; Ksh. Newton Singh
Journal: Journal of Gravity
 (2014)

---

## Abstract

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 intensityE 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 0<u≤0.2086, for all values of charge parameter K lying in the range 0.04≤K≤0.111, and for all values of anisotropy parameter α lying in the range 0.016≥α≥0. 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 ρb=2.7×1014g/cm3 and strange quark stars ρb=4.6888×1014g/cm3 . For K=0.06 and α=0.01, the maximum mass of neutron star is observed as M=1.53M⊙ and radius R=11.48km. Further for strange quark stars M=1.16M⊙ and R=8.71km are obtained.


---

## Body

## 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 withReissner-Nordstromsolution. 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ρb=2.7×1014 g/cm3 (Astashenok [10]) and mass 1.4M⊙-2.9M⊙. However, Astashenok [10] established that f(R) models with realistic equation of state of neutron star have upper limit mass 2M⊙ and minimal radius close to 9 km.(b)
A strange quark star has surface densityρb=4.6888×1014 g/cm3 (Fatema and Murad [11]; Zdunik [12]) and possible maximum mass 2M⊙. 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 2.4M⊙.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(1)ds2=-eω[dr2+r2(dθ2+sin2⁡θdϕ2)]+eνc2dt2,
where ω and ν are functions of r.The Einstein-Maxwell field equations for a nonempty space-time are(2)Rji-12δjiR=-8πGc4Tji=-8πGc4[(-FimFjm+14δjiFmnFnm)(p⊥+ρc2)vivj-p⊥δjihihhhh+(pr-p⊥)χjχihihhhh+14π(-FimFjm+14δjiFmnFnm)],
where Rji is Ricci tensor, Tji is energy-momentum tensor, R is the scalar curvature, Fjm is the electromagnetic field tensor, pr denotes the radial pressure, p⊥ is the transversal pressure, ρ is the density distribution, χi is the unit space-like vector in the radial direction, and vi is the velocity vector, satisfying the relation
(3)gijvivj=1.
Since the field is static,
(4)v1=v2=v3=0,v4=1g44.
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:
(5)8πGc4pr=e-ω[ω′24+ω′r+ω′ν′2+ν′r]+q2r4(6)8πGc4p⊥=e-ω[ω′′2+ν′′2+ν′24+ω′2r+ν′2r]-q2r4(7)8πGc2ρ=-e-ω[ω′′+ω′24+2ω′r]-q2r4,
where prime (′) denotes differentiation with respect to r. From (5) and (6) we obtain the following differential equation in ω and υ:
(8)e-ω[ω′′2+ν′′2+(υ′)24-(ω′)24-ω′·ν′2-ω′2r-υ′2r]-2q2r4=8·π·Gc4(p⊥-pr).
Our task is to explore the solutions of (8) and to obtain the fluid parameters pr, p⊥, 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 E of the following form:
(9)E2C=q2Cr4=KCr2(1+Cr2)-10/6B2,
where K 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, E remains regular and positive throughout the sphere. In addition, E vanishes at the center of the star and increases towards the boundary.We also take(10)8·π·Gc4(p⊥-pr)=Δ=2·α·C2r2(1+Cr2)-10/6B2,
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 ofeω and eυ.(ii)
The radial pressurepr must be vanishing, but the tangential pressure p⊥ may not vanish at the boundary r=rb 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 pr, p⊥ should be positive inside the star.(iv)
(dpr/dr)r=0=0 and (d2pr/dr2)r=0<0 so that pressure gradient dpr/dr is negative for 0≤r≤rb.(v)
(dp⊥/dr)r=0=0 and (d2p⊥/dr2)r=0<0 so that pressure gradient dp⊥/dr is negative for 0≤r≤rb.(vi)
(dρ/dr)r=0=0 and (d2ρ/dr2)r=0<0 so that density gradient dρ/dr is negative for 0≤r≤rb.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,0≤dpr/c2dρ<1 and 0≤dp⊥/c2dρ<1. In addition to the above the velocity of sound should be decreasing towards the surface; that is, (d/dr)(dpr/dρ)<0 or (d2pr/dρ2)>0 and (d/dr)(dp⊥/dρ)<0 or (d2p⊥/dρ2)>0 for 0≤r≤rb; 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ρ≥pr+2p⊥ and ρ+pr+2p⊥≥0.(ix)
The central red shiftZ0 and the surface red shift Zb should be positive and finite; that is, Z0=[e-ν/2-1]r=0>0 and Zb=[e-νb/2-1]>0 and both should be bounded.(x)
Electric intensityE, such that Er=0=0 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 E and the anisotropy factor Δ in such a manner that the solution can be obtained and physically viable. Thus we have
(11)eω/2=B(1+Cr2)-1/6,x=Cr2,y=dυdxq2=KB2Cx3(1+x)-10/6,8·π·Gc4(p⊥-pr)=Δ=2·α·C·x·(1+x)-10/6B2.
On substituting the above in (8), we get the following Riccati-differential equation in y:
(12)dydx+y3(1+x)+y22=36K+36α-1036·(1+x)2
which yields the following solution:
(13)eν/2={1+A(1+Cr2)2S/12}(1+Cr2)(4-S)/12B2,
where A, B, C, and K are arbitrary constants and
(14)S=72(K+α)-4,
where S is real for K+α≥1/18.The expressions for density, radial pressure, and transversal pressure are given by(15)8πGρc2=C36B2f20(72+20Cr2-36KCr2)(16)8πGprc4=C36B2f20{1+Af2S}×[Cr2[Af2S(36K+12+8S)hhhhhh+(36K+12-8S)f2S]+12[Af2S(2+S)+(2-S)]](17)8πGp⊥c4=C36B2f20{1+Af2S}×[Cr2[Af2S(36K+72α+12+8S)hhhhhh+(36K+72α+12-8S)f2S]+12[Af2S(2+S)+(2-S)]],
where
(18)f=(1+Cr2)1/12.


## 5. Properties of the New Solution

The central values of pressure and density are given by(19)(8πGprc4)r=0=(8πGp⊥c4)r=0=C3B2(1+A)[2(1+A)+S(A-1)](20)(8πGρc2)r=0=2CB2.
The central values of pressure and density will be nonzero positive definite if the following conditions are satisfied:
(21)A>(S-2)(S+2),C>0.
The pressure to density ratios are given by
(22)prρc2=(Cr2[f2SA(36K+12+8S)+(36K+12-8S)]+12·[f2SA(2+S)+(2-S)])×((1+Af2S)[72+20Cr2-36KCr2])-1(23)p⊥ρc2=(Cr2[f2SA(36K+72α+12+8S)hhhhh+(36K+72α+12-8S)f2S]+12·[f2SA(2+S)+(2-S)])×((1+Af2S)[72+20Cr2-36KCr2])-1.
Subjecting to the condition that the ratio of pressure to density should be positive and less than 1 at the centre, that is, p0/ρ0c2≤1, the following inequality has to be satisfied.                         (24)(prρc2)r=0=(p⊥ρc2)r=0=2(1+A)+SA-S6(1+A)=16-S(1-A)6(1+A)≤1.
All the values of A which satisfy (21) will also lead to the condition p0/ρ0c2≤1.Differentiating (16) with respect to r, we get
(25)8πGc4dprdr=C2r216B2{1+Af2S}2f32×[-8·Cr2×[A2f4S{36K+12+8S}hhhhhhhhhhh+{36K+12-8S}hhhhhhhhhhh+2A·f2S(36K+12-2S2)]+12·[A2f4S{36K-28-12S}hhhhhhhh+(36K-28+12S)hhhhhhhh+2A·f2S(36K-28+2S2)]Cr2].
Thus it is found that maxima of pr occur at the centre; that is,
(26)pr′=0⟹r=0,8πGc4(pr′′)r=0=-ve.
Thus the expression of the right-hand side of (26) is negative for all values of A satisfying condition (21), showing thereby that the pressure (pr) is maximum at the centre and monotonically decreasing.Differentiating (17) with respect to r, we get
(27)8πGc4dp⊥dr=C2r216B2{1+Af2S}2f32×[-8·Cr2×[A2f4S{36K+72α+12+8S}hhhhhhhhhhhhhh+{36K+72α+12-8S}hhhhhhhhhhhhhh+2A·f2S(36K+72α+12-2S2)]+12·[A2f4S{36K+72α-28-12S}hhhhhhhhhh+(36K+72α-28+12S)hhhhhhhhhh+2A·f2S(36K+72α-28+2S2)]Cr2].
Thus it is found that maxima of p⊥ occur at the centre; that is,
(28)p⊥′=0⟹r=0,8πGc4(p⊥′′)r=0=-ve.
Thus the expression of the right-hand side of (28) is negative for all values of A satisfying condition (21), showing thereby that the transversal pressure is maximum at the centre and monotonically decreasing.Now differentiating equation (15) with respect to r we get
(29)8πGc2dρdr=C2r216B2f32[-1200-432K-160·Cr2+288·KCr2].
Thus the maxima of ρ occur at the centre; that is,
(30)ρ′=0⟹r=0,8πGc2(ρ′′)r=0=(-1200-432K)C2216B2.
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,(1/c2)(dp/dρ)r=0, is given by
(31)1c2(dprdρ)r=0=(36K-28)(1+A)2+12S+4AS2-12SA2(-100-36K)(1+A)2<1andpositive(32)1c2(dp⊥dρ)r=0=(36K+72α-28)(1+A)2+12S+4AS2-12SA2(-100-36K)(1+A)2<1andpositive.
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 coefficientsgij and their first derivatives gij,k in interior solution (I) as well as in exterior solution (E) are continuous up to and on the boundary B. The continuity of metric coefficients gij of I and B on the boundary is the known first fundamental form. The continuity of derivatives of metric coefficients gij of I and B 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(33)ds2=(1-2GMc2R+q2R2)c2dt2-(1-2GMc2R+q2R2)-1dR2-R2(dθ2+sin2θdϕ2),
where M is the mass of the ball as determined by the external observer and R 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:
(34)eυb=[1-2GMc2Rb+qb2Rb2](35)Rb=rb·eωb/2,q(r=rb)=qb(36)12(ω′+2r)brb=(1-2GMc2Rb+qb2Rb2)1/2(37)12(ν′)brb=(GMc2Rb-qb2Rb2)(1-2GMc2Rb+qb2Rb2)-1/2.
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 “d,” Schwarzschild parameters u=GM/c2Rb, and radius of the star Rb;
(38)C=6(1-d)(6d-4)rb2>0foru≤1036+qb22Rb2(39)A=(6(u-qb2Rb2)·(1+Crb2)(4-S)/12-(4-S)·Crb2·d·(1+Crb2)(-8-S)/12(u-qb2Rb2))×((u-qb2Rb2)(4+S)·d·Crb2·(1+Crb2)(-8+S)/12-6(u-qb2Rb2)(1+Crb2)(4+S)/12)-1(40)B=(1+Crb2)(4-S)/12+A(1+Crb2)(4+S)/12d,
where we define a new parameter called Reissner-Nordstrom parameter “d” given by
(41)d=[1-2u+qb2Rb2]1/2
whose value lies within 2/3<d<1 for Crb2>0.Surface density is given by(42)8πGρc2Rb2=Crb236(1+Crb2)2[72+20Crb2-36KCrb2].Central red shift is given by(43)Z0=[B21+A-1].The surface red shift is given by(44)Zb=[e-νb/2-1].


## 7. Discussions and Conclusions

From Figures1 and 2 it has been observed that the physical quantities pr, p⊥, pr/ρc2, p⊥/ρc2, dpr/c2dρ, dp⊥/c2dρ, and z are positive at the centre and within the limit of realistic state equation and monotonically decreasing while the quantities γ,Δ,E are increasing for all values of K, α, and u lying in the ranges 0.04≤K≤0.111, 0<α≤0.016, and 0≤u≤0.2086, respectively. For K+α>0.111 the pressure is negative. With the increase in the value of K from 0.04 to 0.111 the Schwarzschild parameter “u” increases; hence the mass increases, but the value of α has to be decreased to 0. With α=0 we recover the isotropic model.Figure 1
The variation ofp, ρ, p/ρc2, Z, γ, and E from centre to surface for K=0.06 and α=0.01.Figure 2
The variation of(1/c2)(dp/dρ) and anisotropy from centre to surface for K=0.06 and α=0.01.By increasingK 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 u satisfying the inequality 0≤u≤0.2086 for K 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.Figure 3
The variation of mass with radius of the superdense stars.In Tables1 and 2 we present a model of superdense neutron star and quark star based on the particular solution discussed above. By assuming surface density, ρb=4.6888×1014 g/cm3 corresponding to α=0.01, K=0.06 for which u=0.198, the resulting well-behaved solution has a maximum mass M=1.16 M⊙ and radius R=8.71 km (for quark star) and by assuming the surface density ρb=2.7×1014 g/cm3 the obtained maximum mass is M=1.53 M⊙ and radius R=11.48 (for neutron star).Table 1
Variation of maximum mass of neutron star and corresponding radiusRb with u for K = 0.06 and α = 0.01 by assuming the surface density ρb=2.7×1014 g/cm3.
u
(
8
π
G
/
C
2
)
ρ
r
b
2
R
b (km)
M
/
M
⊙
0.010
0.08605
3.41575
0.02299
0.020
0.16661
4.77451
0.06413
0.040
0.31173
6.59067
0.17622
0.060
0.44198
7.92405
0.32153
0.080
0.54874
8.91444
0.48155
0.100
0.63497
9.68137
0.65060
0.120
0.70596
10.31179
0.83478
0.140
0.75730
10.78723
1.01872
0.160
0.79117
11.13602
1.20163
0.180
0.80834
11.36756
1.38018
0.198
0.81007
11.47792
1.53364Table 2
Showing the variation of maximum mass of Quark star (ρb=4.6888×1014 g/cm3), its radius, surface red shift, surface electric intensity, and surface anisotropy with u for K = 0.06 and α = 0.01.
u
R
B (km)
M
/
M
⊙
Z
b
E
b
r
b
Δ
b
r
b
2
0.010
2.59201
0.01745
0.010101
0.008901
0.000026
0.020
3.62309
0.04866
0.020408
0.017722
0.000105
0.040
5.00127
0.13373
0.041667
0.035122
0.000411
0.060
6.01310
0.24399
0.064963
0.053045
0.000938
0.080
6.76465
0.36542
0.088850
0.070282
0.001646
0.100
7.34663
0.49370
0.113586
0.087017
0.002524
0.120
7.82502
0.63347
0.140901
0.104300
0.003626
0.140
8.18579
0.77305
0.169317
0.121076
0.004886
0.160
8.45047
0.91185
0.199472
0.137677
0.006318
0.180
8.62618
1.04734
0.231679
0.154199
0.007926
0.198
8.70992
1.16379
0.262626
0.169054
0.009526



---

*Source: 380320-2014-09-02.xml*