SciELO - Scientific Electronic Library Online

vol.65 issue4Heat transfer augmentation in water-based TiO2 nanoparticles through a converging/diverging channel by considering Darcy-Forchheimer porosityA model for low mass compact objects author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand




Related links

  • Have no similar articlesSimilars in SciELO


Revista mexicana de física

Print version ISSN 0035-001X

Rev. mex. fis. vol.65 n.4 México Jul./Aug. 2019  Epub May 06, 2020 


A charged perfect fluid model with high compactness

G. Estevez-Delgadoa 

J. Estevez-Delgadob 

M. Pineda Duranc 

N. Montelongo Garcíad 

J.M. Paulin-Fuentese 

aFacultad de Químico Farmacobiología de la Universidad Michoacana de San Nicolás de Hidalgo, Tzintzuntzan No. 173, Col. Matamoros, 58240, Morelia Michoacán, México. *e-mail:

bFacultad de Ciencias Físico Matemáticas de la Universidad Michoacana de San Nicolás de Hidalgo, Edificio B, Ciudad Universitaria, 58030, Morelia Michoacán, México. *e-mail:

cInstituto Tecnológico Superior de Tacámbaro, Av. Tecnológico No 201, Zona el Gigante, C: P. 61651, Tacámbaro Michoacán, México. e-mail:

dInstituto Tecnológico de Morelia, Av. Tecnológico 1500, Lomas de Santiaguito 58120, Morelia Michoacán, México. Universidad Tecnológica de Morelia, Av. Vicepresidente Pino Suárez 750, 4ta Etapa de Ciudad Industrial 58200, Morelia Michoacán, México. e-mail:

eDivisión Académica de Ciencias Básicas de la Universidad Juárez Autónoma de Tabasco, Carretera Cunduacán-Jalpa KM. 1, Col. La Esmeralda, 86690, Cunduacán, Tabasco, México. e-mail:


A relativistic, static and spherically symmetrical stellar model is presented, constituted by a perfect charged fluid. This represents a generalization to the case of a perfect neutral fluid, whose construction is made through the solution to the Einstein-Maxwell equations proposing a form for the gravitational potential gtt and the electric field. The choice of electric field implies that this model supports values of compactness u = GM / c2R ≤ 0.5337972212, wich are higher than the case without electric charge (u = 0.3581350065), being this feature of relevance to represent compact stars. In addition, density and pressure are positive functions, bounded and decreasing monotones while the electric field is a monotonously increasing function as well as satisfying the condition of causality, so the model is physically acceptable. Additionally, the internal behavior of the hydrostatic functions and their values are obtained taking as data the corresponding to a star of 1 Mʘ for different values of the charge parameter, obtaining an interval for the central density ρc ≈ (7.9545, 2.7279)1019 Kg/m3 characteristic of compact stars.

Keywords: Exact solutions; perfect fluid; stars solutions

PACS: 04.40.Nr; 04.20.Jb; 04.20.Dw

1. Introduction

The compact objects: neutron stars, pulsars, strange stars have been approached from different perspectives, mainly through particle physics and gravitation, which have allowed us to describe their behavior at the microscopic and macroscopic levels respectively [1,2]. A consequence of the investigations carried out since the last century leads us to the fact that compact stars are constituted by electrons, neutrons and possibly quarks and that their equilibrium is in part due to the degeneration pressure among the particles. In the case of white dwarf stars this pressure of degeneration is due to the electrons and this consideration led to determine the so-called Chandrasekhar limit for the mass of this type of stars that corresponds to 1.44 solar masses [3,4], while the degenerating pressure of the neutrons implies that the mass of the neutron stars is approximately Mʘ. Stellar objects even more dense than these are constituted by quarks, although in reality these stars are hybrid, constituted predominantly by some of these particles. In the astrophysical aspect when we consider the microscopic behavior, particle interaction, we can introduce some equation of state associated to MIT bag model for the quark matter P = (1/3)(c2ρ - Bg) [5] or a polytropic state equation P = μρ1+1/n [6] if we refer to a neutron star or a white dwarf. Of course, there is a possibility that the interior of a star is decried by some other state equation P = P(ρ) depending on density orders [7-9] even if the state equation is not expressed implicitly, that is, F(P, ρ, S) = 0, where S represents some parameters or, in the case of electrically charged objects, include an electric charge relation. This may be the approach, although it is not indicated, present in some research reports in which solutions are found to Einstein’s equations with some source of matter that allows modeling compact stellar objects. Sometimes it is possible to determine a state equation, in the static and spherically symmetric case, since the pressure P = P(r) and density ρ = ρ(r) could get to express one in function of the other P = P(ρ) = 0 by inversion of the radial coordinate [10], although in general this does not happen due to the complexity of the solution.

However, the construction of stellar models, although it does not allow us to obtain any equation of state, guides us in a better understanding of the inner behavior of the stars through the use of gravitational theories such as the theory of Einstein’s General Relativity. For this theory, the vast majority of models or analyzes proposed near the interior of the star consider that it is described by a perfect fluid [11,12], a perfect fluid charged [13,14], an anisotropic fluid [13,14] or a charged anisotropic charged perfect fluid [13,18]. The class of solutions presented in each case attends to the type of objects that one wishes to describe or understand. The perfect fluid models are used for less compact stars, while the anisotropic or charged models allow to represent more compact stellar objects. Regarding the last case, there is a minimum value [19] for the ratio of mass and radius which generalizes to Buchdal [20]. Through a graphical analysis it has been shown that the compactness is even greater than the upper limit for a perfect fluid [21]. In this work we present a charged perfect fluid interior solution in the context of the Einstein- Maxwell theory, which is physically acceptable, i.e., satisfied with this feature.

The construction of a static and symmetrically spherical space-time is proposed from the assumption of a specific form of gravitational potential gtt [22] and the magnitude of the electric field. The solution satisfies conditions that make it physically acceptable and allows this to represent astrophysical compact stars with a compactness u ≤ 0.53379772212. In the Sec. sec2 we present the Einstein - Maxwell field equations for a charged perfect fluid charged fluid and the conditions are given for the solution to be physically acceptable. In the Sec. sec3 the solution of the field equations is presented and the validation intervals of the parameters are determined considering the physical conditions that must be satisfied, while in the Sec. 4 we determine the type of astrophysical objects that can be modeled and shown by graphical representations of the internal behavior, as well as an analysis of the possible values of density and pressure inside according to a stellar object with a solar mass. We finished the research work by presenting some conclusions in Sec. 5.

2. The system

Einstein-Maxwell field equations of gravitation for a perfect fluid charged are given for Gμν=Tμν(fp)+Tμν(E) having the contribution of two parts, where [23]

TμvE=k c2ρ+Puμuv+Pgμv,

TμvE= 14πFμαFvα-14FαβFαβgμv (1)

with P representing the pressure distribution, ρ the density distribution measured by an observer with velocity vector uμ, Fμν the Maxwell tensor and k = 8πG / c4. We consider the static and spherically symmetric metric of the interior in Schwartzschild coordinates [24]

ds2=-y(r)2dt2+dr2B(r)+r2(dθ2+sin2θdϕ2). (2)

We consider that the contribution of the Maxwell tensor is due to an electric field, so the non-zero components of the electric part are:

TEt= TErr=-TEθθ=TEϕϕ=B2y2Ftr2=Er2



with q(r)the total charged inside a sphere of radius r. The field equations are write as [25]:

kc2ρ+E2r=-B'r+1-Br2 (4)

kP-E2r=2By'ry-1-Br2 (5)

kP+E2r=ry''+y'Bry-ry'+yB'2ry (6)

where ' denotes the derivative with respect to the coordinate r. This is the system of equations that we must solve for the construction of a stellar model given by a perfect charged fluid. The system has three Eqs. (4)-(6) for five functions (ρ, P, E, y, B) so two restrictions are required, in our case these are the shape of the gravitational potentiall gtt = -y2 and the function of the magnitude of the field E2, the choice of the first of these functions is motivated by previous work in which a perfect fluid neutral model has been developed, while the magnitude of the electric field is proposed so that we can represent objects with greater compactness than in the neutral case and with greater diversity of behavior. The solution to the Eqs. (4)-(6) requires to satisfy conditions that make it physically acceptable, which will be described in the next section.

2.1. Physicals conditions

For an inner solution of the Einstein-Maxwell equations with perfect fluid associated with compact objects to describe a physically acceptable model, the following conditions must be fulfilled [12]:

  • The solution must not have singularities, geometric or physical variables, i.e., for 0 ≤ rR the curvature scalars must be regular and the metric functions (y2, B), the density and pressure must be bounded.

  • The pressure and density must be positive and monotonically decreasing functions as a function of radial distance, with its maximum value in the center, in, particular in the origin:



while for r ≠ 0 ρ' < 0 and P' < 0.

  • The energy conditions must be [26]

  • NEC null energy condition ρ + P ≥ 0

  • WEC weak energy condition ρ ≥ 0 and ρ + P ≥ 0

  • SEC strong energy condition ρ + 3P ≥ 0 and ρ + P ≥ 0

  • DEC dominant energy condition ρ ≥ 0 and ρ ± P ≥ 0

From the previous requirement for density and pressure we have that the only additional restriction corresponds to the DEC.

The causation condition must not be violated, i.e. the magnitude of the speed of sound must be less than the speed of light


and additionally we will impose that the speed of sound is a monotonous function decreasing towards the surface.

  • For the stability of the solution, in the relativist case, it is required that the adiabatic index [1,27]


  • There must be a region r = R, the surface of the star, where the pressure is zero P(R) = 0.

  • On the boundary r = R the interior solution should match continuously with an exterior Reissner-Nordstrom solution

ds2=-1-2GMc2r+Q2r2dt2+1-2GMc2r+Q2r2-1dr2+r2dθ2+sin2θdϕ2, rR (7)

This requires the continuity of y2(r), B(r) and q(r) = Q across the region r = R, where M and Q represent the total mass and charge inside the fluid sphere respectively.

  • Electric field intensity E is such that E(0) = 0 and taken to be monotonically increasing, i.e., dE / dr > 0 for 0 < r ≤.

These basic requirements allow to determine which interior solution can be useful as a model for the description of some compact object. As a result of the conditions mentioned above it is known that one characteristic of these models is that the radio mass ratio, as a generalization to the limit of Buchdahl [20] GM / c 2 R ≤ 4/9, has a maximum value

GMc2R3+c2R9+Q23R, (8)

this expression implies the possibility of having solutions with a greater compactness value than 1/2 due to the effect of the charge Q, property that is relevant to our model. In addition to this inequality there is also a lower bound for the radio mass ratio


both relations have a good limit for the neutral case Q = 0.

Another important quantity in the stellar models that is associated with the mass and the radius is the gravitational redshift on the surface zb=gtt-(1/2)R-1 that in the case of charged solutions also involves the charged. The continuity of the metric on the surface implies


3. The solution and their analysis

The system of equations described in the previous section supports a solution with

yr=S1+10ar21+ar23/2, (9)

in the case that the perfect fluid is neutral. In most previous investigations the form of this function has been proposed as y(r)=(1+ar2)n [12,29] which has allowed the construction of physically acceptable interior solutions. Nowadays it is known that this form guarantees the regularity of the geometry in the vicinity of the center, since for this it is only required that [30]




So other functions can be proposed with these properties. The proposed form for the function y(r) and other similar functions have allowed to show the relevance of the existence of anisotropic pressures to be able to have physically acceptable models [31] besides that they have been applied for the description of physically acceptable stellar models [30,32,33], some of which are characterized because the speed of sound is a decreasing monotonic function as a function of the radial coordinate [34].

The construction of interior solutions is useful to have a better clarity of the inner behavior of the stars [12], some solutions with perfect fluid sources without charge [29,35,36] these became generalized to charged case [37-39] to represent compact stars.

The proposal of this section is to obtain a charged model that generalizes to the previously constructed case [22]. Some of the advantages presented by charged models is that their compactness ratio becomes greater than their counterpart without charge as a result of the non-neutrality of the fluid. The behavior of the electric charge or equivalently of the electric field E = q(r)/r2 as already mentioned in the previous section, it must be zero in the center and it must be a growing monotonous function, so in our case we choose the electric field in the form

E2r=271+16ar2-30a2r4a2vr221+10ar221+10ar22 (10)

v ≥ 0represents the charge parameter that in the case v = 0 will allow us to recover the solution without electric charge [22]. The shape of the electric field magnitude is not unique, there are a variety of possible acceptable functions for the same function y(r) [37-39], the choice of this serves particular interests of the approach of the model to be proposed or of the research proposed, even this can be given through the charge density [40] or through a relationship with the state equation [41].

In our case this form of the electric field was chosen so that the class of objects that can be represented can have a greater compactness, and this will be shown in the next section. Replacing this form of electric field and potential y in the Eqs. (5) y (6), after subtracting them, we come to the differential equation:


here we have defined H0(r) = 1 + 16ar2 - 30a2r4. The solution to this equation implies:

B(r)=1+a(ν+2)r2+(1+10ν)a2r41+10ar2+100Car2(1+ar2)5H-3313-2(1+10ar2)H+3313+2, (11)

Where C is the constant of integration and


This determines the solution, i.e., the geometry and the hydrostatic functions, so that when replacing (9), (10) y (11) in the Eqs. (4) and (5) we obtain

kc2ρr=300Ca1+31ar2+S11+ar24H_3-313H+3313+31+10ar22-32+53ar2+426a3r4+170a3r6+200a4r8av21+10ar221+ar22+3a8+25ar2-10a2r41+10ar22 (12)

kPr=100Ca1+65ar2+10a2r41+ar22H_2-3313H+2+33131+10ar22+2+179ar2+190a2r4+710a3r6+200a4r8av21+10ar221+ar22+46-25ar2+10a2r4a1+10ar22 (13)

where in the expression of density, Eq. (12), S1 = 150a2r4 + 1330a3r6 +10a4r8. Given the conditions that must be satisfied the solution to be physically acceptable, it is necessary to determine the speed of sound


due to its extension it is not written in the text, although it will be considered in the analysis of the conditions. The metric functions as well as the density, pressing and electric field are regular in the center, however, the explicit form of these is required in order to obtain the intervals of the system constants (a, C, v) as a result of the physical conditions described in the previous section. Then, from equations (12) and (13) as well as how the speed of sound is, we have:

kc2ρ0= -3H-03313H+0-3313Ca+38-va (14)

kP0=H-03313H+0-3313Ca+46+va (15)

v20c2=2C-5H-0-3313H+0331314-v10C-H-0-3313H+0331330+v (16)

According to the required conditions, each of these relationships must be positive, in addition to the speed of sound that must be smaller than the speed of light. On the other hand as ρ'(0) = 0, P'(0) = 0, density and pressure are decreasing monotone functions with their maximum in the center, so it is required that their second derivatives

kc2ρ''0=2710H-03313H+03313C-30+va2 (17)

kP''0=27 10H-03313H+03313C-514+va2 (18)

be negative, while the electric field satisfies E(0) = 0, E'(0) = 0 y E2''(0) = 27va2; then v ≤ 0 as had been assumed in the Eq. (10). To the set of inequalities generated from the conditions required in the center for the solution to be physically acceptable expressed by (14)-(18), it is convenient to complement them with the relation obtained from the existence of the surface of the star, identified by P(R) = 0. Considering (13) in r = R we can get the constant C:

C=-s3wv-246-25w+10w21+w2H-0+10w3313-2H+0+10w-3313-2S2w (19)

where w = aR2 with S2(w) = 200(1+65w+10w2)(1+w)6 and S3(w) = 2+179w+1902w 2+170w3+200w4. The set of inequalities is expressed only in terms of the dimensionless parameters (v,w). However, these conditions are only the content of the behavior in the object center. The parameter w is restricted by the condition that the electric field must be a growing monotonous function for all r ∈ [0,R]. This happens if inside the electric field has no maximum and the limit case occurs when at the border the electric field has a maximum, this is E2(R)' = 0 which implies that

1+21w-120w2-650w3+300w4=0, (20)

since w > 0 then w is constrained by the smallest positive root of this equation, which is w = 0.13828732, this value of w is greater than the maximum value allowed when you have the perfect fluid without load whose value is w(q=0) = 0.1073273425. The maximum value of the parameter v it is determined by the condition that the speed of sound must be positive. The behavior of this function, for the solution we are analyzing, requires that the speed of sound at the border be positive and this happens if vvmax where


and S4 = 129600w468420w5 + 16200w6 + 1000w7. So the maximum value of v depends on w. The minimum value of v also depends on the specific value of w and is determined by the condition that the speed of sound is smaller than the speed of light, the value of r where this happens changes depending on the value of w. In particular for the maximum value of w = 0.13828732 the speed of sound matches the speed of light for v = 0.274043033 in r = 0.572156468R, so the range of v ∈ (0.27404 3033. 2.84837494). In general, the interval for v ∈ (vmin; vmax) and the value of vmin are not being determined so that the speed of sound is smaller than the speed of light.

4. Graphic analysis

The continuity of the metric and electric field allows us to determine the relationship which determines the compaction of the stellar objects that can be represented with the solution construed. The external geometry is given by the Reisner-Nordstrom solution

ds2= -frdt2+fr-1dr2+r2dθ2sin2θdϕ2

fr=1-2GMc2r+Q2r2, Er=Q2r2,  rR, (21)

where M and Q are the mass and the electric charge of the star respectively. Evaluating the interior metrics in r = R and matching terms we get the reason for compactness

u=GMc2R=271+16w-30w21+38w+10w2w2ν21+10w21+w21+65w+10w2+27w1+65w+10w2. (22)

The value of compactness u = u(w,v) is a monotonically increasing function of both variables so the maximum value of compactness occurs for the maximum value of v and w. Evaluating this relation for the maximum value of v = vmax we obtain


where S5 = 7100w5 + 1000w6, while evaluating this expression for the maximum permitted value of w w = 0.13828732 we obtain that the maximum value of possible compactness is umax = 0.5337972212, which compared to the neutral case that has a compactness value of u = 0.3581350065, the effect of charge allows the representation of more compact objects compared to the case without charge. Models not electrically charged allow a compactness value accepted by the limit of Buchdahl, however for the charged case the radius of the star could be smaller than the radius of Schwarzschild u > 1/2 as it happens in our model, the reason is associated with the effect of the charge. In the neutral case u = 1/2 indicates that the radius of the star coincides with the radius of Schwarzschild rs = 2GM / c2 which would imply the existence of a black hole. However when the Einstein-Maxwell equations are considered in a static space-time, spherically symmetric and asymmetrically flat, the respective black hole is given by the geometry of Reissner - Nordstrom represented by the metric (7) and in this case the external event horizon r+=GM/c2+G2M2/c4-Q2. In particular for Q = GM / c2 is obtained r+ = GM / c2, so the event horizon happens for u = 1, value that is greater than the neutral case, although not necessarily the maximum compactness value of charged case occurs for u = 1. Its maximum value of compactness is limited by inequality (8), from which it can be observed that u could be greater than 1/2 since in this the effect of the charge appears adding up as a quadratic quantity. Now we will present the results graphically for the maximum value of w = 0.13828732 with different values of v which generate the respective compactness values. From the Table I is observed that the maximum value of compactness is greater than 0.5 and it happens when the charge parameter v is maximum in conjunction with w maximum. As it happens for this maximum value of w in which the compactness decreases with the electric charge parameter, this behavior is shown for other values of w < wmax . Figure 1 corresponding to the speed of sound shows that for w = wmax there are regions of increasing monotony behavior near the center and decreasing monotony near the surface. The maximum value of compactness is determined because the speed of sound is zero at the border (blue solid line), while the minimum value (black solid line) is obtained because inside the speed the sound is equal to the speed of the light in x = 0.5721564648. In general, the behavior of this function allows to determine the ranges of the parameter v and consequently the value of compactness for a fixed value of w.

Table I Values of the compactness for different values of the charged parameter with w = 0:13828732. 

v 2.8483749 2.204792 1.5612094 0.9176260 0.2740430
u 0.5337971 0.4960593 0.4583215 0.4205838 0.3828460

Figure 1 Behavior of the speed of sound for different values of compactness. 

To graphically represent the behavior of density we define the dimensionless variable ρ → kc 2 R 2ρ as a function of x = r/R. On graph 2 its regular behavior is shown, bounded and monotone decreasing for different values of the parameter v with w = wmax. The central density is greater for the larger values of the charge parameter v while on the surface the opposite occurs. The graphic representation, Fig. 3, of the dimensionless function associated with pressure PkR2P in terms of the function of x = r/R shows its monotonically decreasing behavior as well as bounded and positive. By the effect of the electric charge manifested through the parameter v the pressure associated with the perfect fluid is lower for higher values of this which allows denser objects. As can be seen from the comparison of Figs. 2 and 3 for greater compactness the density is greater and the pressure is lower. Figure 4 shows that the solution satisfies the condition required for the adiabatic index, γ > 4/3, which guarantees its stability.

Figure 2 Density behavior with w = 0:13828732 for different values of compactness. 

Figure 3 Pressure for different values of compactness u with a fixed value of w = 0:13828732. 

Figure 4 Graphical representation of the magnitude of the electric field for w max

From Fig. 5 an increase of the electric field strenght is observed as the electric charge increases.

Figure 5 Adiabatic index. 

Another condition imposed for a model to be physically acceptable is that the energy conditions are satisfied since the density and the positive pressure are the same as the magnitude of the electric field E2. The only one that requires a verification is the dominate energy condition and this is shown graphically in the Fig. 6 where it is observed that this is satisfied.

Figure 6 Pressure for different values of the parameter y

The graphical representation of the hydrostatic functions presented in terms of dimensionless functions, is helpful to observe their behavior, however to obtain the ranges of the possible physical values of these functions, it is necessary to introduce the units to see that the values which the model are consistent with orders of magnitude associated with this type of objects. To determine the orders of magnitude of the density and the pressure below we consider a star with a mass equal to the mass of the sun and we obtain the values of the radius, pressure and central density as well as its density at the surface for the values of the compactness used for the graphics.

From the Table II, it is observed that the ranges of values of the central density ρc and the density on the surface ρb , are of the order of magnitude associated with compact stars. So also the central pressure Pc is consistent with the expected values. The values of the speed of sound in the center νc and on the surface νb are smaller than the speed of light, i.e., these do not violate the condition of causality and even in the case of greater compactness the speed of sound at the surface is zero. From the graphic behavior and the value of the adiabatic index γ in the center given in the table is obtained that the solution is stable, because γ > 4/3. Two other quantities given for the different values of compactness are the net charge and the value of gravitational redshift on the surface, both values are within the expected ranges [40].

Table II Values of the hydrostatic variables in the center and on the surface for different values of the compactness as well as M = Mʘ

u 0.533797 0.496059 0.458322 0.420584 0.382846
R(m) 2765.976 2976.398 3221.473 3510.524 3856.563
ρc (1019 Kg/m3) 7.954465 5.769167 3.909380 1.658865 2.727820
ρb (1018Kg/m3) 5.369969 5.425123 5.303410 5.032169 4.638761
Pc (1036P) 2.421304 2.298425 2.080038 1.875830 1.667535
v2c(c2) 0.402442 0.447095 0.506627 0.589964 0.714951
v2b (c2) 0 0.163272 0.351349 0.570351 0.828586
γ 1.929431 1.443294 1.441248 1.498013 1.643794
Q(1020C) 1.648824 1.560997 1.421715 1.187770 0.713078
zb 1.257805 1.171623 1.094611 1.025251 0.962353

Finally we apply our model for the data of SAX J1808.4-3658 (SS1) star of mass M = 1.435 Mʘ and radio R = 707 km [42], its compactness value is u = 0.2996795, so the value of v is determined from the Eq. (22)


From the analysis of the restrictions for the model to be physically acceptable, we obtain that the validity interval of w ∈ [0.030456, 0.040500]. The only difference in the behavior of the functions for this data and the data for greater compactness is that there is now a subinterval of w ∈ [0.030456, 0.32117] in which the speed of sound is a decreasing monotone function. The Fig. 7 shows the behavior of the speed when there is no charge, w = 0.0405, being this a monotonously increasing function. In the charged case there are different behaviors for 0.032117 < w < 0.0405. The function of the speed of sound presents a region where it is monotonously increasing and another where it is growing monotonous. For w ∈ [0.030456, 0.032117] the speed of sound is monotonically decreasing as a result of the effect of the charge.

Figure 7 Behavior of the speed of sound for different values of compactness. 

In the Table III, the values are shown with the magnitudes of the hydrostatic variables as well as the charge for the same parameter values w, in this it is observed that for the maximum value of the charge, the speed of sound on the surface approaches zero. In addition, the orders of magnitude of the densities are similar to the values reported with other models [40].

Table III Values of the hydrostatic variables in the center and on the surface SS1. 

w 0.040500 0.038309 0.035478 0.032117 0.030456
ν 0 1.475197 3.952775 8.085688 10.79803
Pc (1035P) 2.207465 1.936262 1.526330 1.059583 .8073013
ρc (1018kg/m3) 9.491643 7.652449 8.927432 8.523171 7.748338
ρb (1018 kg/m3) 2.908835 2.767600 2.572201 2.318696 2.183865
vc (c2)2 .5000328 .4219832 .3371526 .2531569 .2166285
vb (c2)2 0.652298 0.465718 0.264912 0.076179 0.000004
Q(1020C) 0. 1.223631 1.875167 2.459075 2.711817
γ 1.908003 2.109820 1.911706 2.236546 1.455289

5. Conclusions

Here we have obtained a charged perfect fluid solution that generalized a perfect fluid uncharged solution in a spacetime static and spherically symmetric. The solution depends of two parameters w and v the latter related to charge. The built solution considers a new form of gravitational potential gtt = (1+10ar3)3/(1+ar2)3 [22] and a particular form of the increasing monotonous charge, this is the one that generates different profiles of behavior of the speed of sound, differing from the case without charge where the speed is only monotone increasing. Another point in which the electric charge generates changes with respect to the case without electric charge is that the maximum umax = 0.5337972212 is much higher than at the maximum possible value for a perfect neutral fluid u < 4/9.

This remarkable characteristics of the model, as far as we know, have not been reported previously analytic solutions with such large compactness values, although in previous reports mention has been made of the possibility of having electrically charged models with greater compactness than 4/9, upper bound for the case of a perfect fluid neutral [21].

As a result of the analysis of the solution we have that the density and pressure, as radial functions, are monotonous, decrescent, bounded and regular. Additionally, considering the data of a star with mass Mʘ, we have it has that the orders of magnitude that the model generates with respect to density and pressure are consistent with those expected for neutron stars or quarks, with the advantage of having a range of possible values of the density, which, as expected, stars with the same radius and mass data may have different internal behavior. In our case what determines in a specific way the internal structure of the stars is the charge. On the other hand, the adiabatic index is a growing monotonous function and γ > 4/3 which guarantees stability. We also show that the model for smaller values of compactness, in particular for u = 0.2996795 can be associated with the data of the star J1808.4-3658(SS1). There is a parameter range w in which the speed of sound is a decreasing monotonous function, characteristic that some authors consider as desirable for some models.

Finally, it is worth mentioning that although the graphical analysis and physical hydrostatic values were presented for specific parameters, the behavior for other values of the parameters, is similar so it is concluded that the model can be used for other compact stars u < 0.5337972212. This work gives rise to some questions that could be developed in the future, such as what are the characteristics that the electric field function must have that allows to represent compact objects with a wide spectrum of possibilities in the representation of highly compact objects. Another question is the construction and analysis of an anisotropic charged model which generalizes the one presented here and the determination of the effect of the anisotropy and the electric charge and their comparison between both of them.


We appreciate the facilities provided by the Universidad Michoacana de San Nicolás de Hidalgo. The author JMPF would like to thank CONACyT México for the individual research grant received. We thank the reviewer of the work for his suggestions made to clarify some points.


1. S. L. Shapiro and S. A. Teukolsk, WILEY-VCH Verlag GmbH and Co. KGaA. (2004). [ Links ]

2. Norman K. Glendenning, 2nd Edition, (Springer 2000). [ Links ]

3. S. Chandrasekhar, Astrophys. J. 74 (1931) 81. [ Links ]

4. S. Chandrasekhar, Mon. Not. R. Astron. Soc. 95 (1935) 207. [ Links ]

5. E. Witten, Phys. Rev. D. 30 (1984) 272. [ Links ]

6. U. S. Nilsson and C. Uggla, Annals Phys. 286 (2001) 292. [ Links ]

7. A. Lavagno and G. Pagliara, Eur. Phys. J. A S01 27 (2006) 289. [ Links ]

8. J. W. Negele and D. Vautherin, Nucl. Phys. A 207 (1973) 298. [ Links ]

9. G. Baym, C. Pethick and P. Sutherland, Astrophys. J. 170 (1971) 299. [ Links ]

10. R. C. Tolman, Phys. Rev. 55 (1939) 364. [ Links ]

11. J. R. Stone, Eur. Phys. J. A 52 (2016) 66. [ Links ]

12. M. S. R. Delgaty and K. Lake, Comput. Phys. Commun. 115 (1998) 395. [ Links ]

13. B. Dayanandan, S. K. Maurya and T. T. Smitha, Eur. Phys. J. A 53 (2017) 14. [ Links ]

14. N. Pant, N. Pradhan and M. H. Murad, Astrophys Space Sci 355 (2015) 137. [ Links ]

15. S. K. Maurya, Eur. Phys. J. A 53 (2017) 89. [ Links ]

16. S. K. Maurya, Y.K. Gupta, T.T. Smitha and F. Rahaman, Eur. Phys. J. A 52 (2016) 191. [ Links ]

17. S. K. Maurya and S. D. Maharaj, Eur. Phys. J. A 54 (2018) 68. [ Links ]

18. M. H. Murad and S. Fatema, Eur. Phys. J. C 75 (2015) 533. [ Links ]

19. H. Andréasson, Commun. Math. Phys. 288 (2009) 715. [ Links ]

20. H. A. Buchdahl, Phys. Rev. 116 (1959) 1027. [ Links ]

21. M. K. Mak, P. N. Dobson Jr. and T. Harko, Europhys. Lett. 55 (2001) 310. [ Links ]

22. G. Estevez-Delgado, J. Estevez-Delgado, N. Montelongo-García, M. Pineda Duran and J. M. Paulin-Fuentes, Submit to MPLA (2018). [ Links ]

23. D. D. Dionysiou, Astrophys Space Sci 85 (1982) 331. [ Links ]

24. C. W. Misner, K. S. Thorne and J. A Wheeler, Gravitation, (Freeman, San Francisco 1973). [ Links ]

25. J. Kumar and Y. K. Gupta, Astrophys Space Sci 345 (2013) 331. [ Links ]

26. M. H. Murad and S. Fatema, Eur. Phys. J. C 75 (2015) 533. [ Links ]

27. H. Heintzmann and W. Hillebrandt, Astron. Astrophys. 38 (1975) 51. [ Links ]

28. C. G. Böhmer and T. Harko, Gen Relativ Gravit 39 (2007) 757. [ Links ]

29. N. Pant, Astrophys Space Sci 331 (2011) 633. [ Links ]

30. G. Estevez-Delgado and J. Estevez-Delgado, Mod. Phys. Lett. A 33 (2018) 1850081. [ Links ]

31. G. Estevez-Delgado and J. Estevez-Delgado, J. Eur. Phys. J. C. 78 (2018) 673. [ Links ]

32. G. Estevez-Delgado, J. Estevez-Delgado N. Montelongo García and M. Pineda Duran, Mod. Phys. Lett. A 33 (2018) 1850237. [ Links ]

33. G. Estevez-Delgado, J. Estevez-Delgado N. Montelongo García and M. Pineda Duran, Can. J of Phys. Acepted 2018. [ Links ]

34. G. Estevez-Delgado, J. Estevez-Delgado, J. M. Paulin-Fuentes, M. Pineda Duran and N. Montelongo García, Rev Mex Fis. Acepted (2019). [ Links ]

35. M. C. Durgapal, J. Phys. A 15 (1982) 2637. [ Links ]

36. H. Heintzmann, Z. Phys. 228 (1969) 489. [ Links ]

37. Y. K. Gupta and S. K. Maurya, Astrophys Space Sci 332 (2011) 155. [ Links ]

38. S. K. Maurya and Y.K. Gupta, Astrophys Space Sci 332 (2011) 481. [ Links ]

39. N. Pant, R.N. Mehta and M. (Joshi) Pant Astrophys Space Sci 332 (2011) 473. [ Links ]

40. S. Ray, A.L. Espindola, M. Malheiro, J.P.S. Lemos, and V.T. Zanchin, Phys. Rev. D 68 (2003) 084004. [ Links ]

41. B. V. Ivanov, Phys. Rev. D 65 (2002) 104001. [ Links ]

42. M. Kalam et al., Int J Theor Phys 52 (2013) 3319. [ Links ]

Received: December 02, 2018; Accepted: February 04, 2019

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License