Compact stars with a small electric charge: the limiting radius to mass relation and the maximum mass for incompressible matter

One of the stiffest equations of state for matter in a compact star is constant energy density and this generates the interior Schwarzschild radius to mass relation and the Misner maximum mass for relativistic compact stars. If dark matter populates the interior of stars, and this matter is supersymmetric or of some other type, some of it possessing a tiny electric charge, there is the possibility that highly compact stars can trap a small but non-negligible electric charge. In this case the radius to mass relation for such compact stars should get modifications. We use an analytical scheme to investigate the limiting radius to mass relation and the maximum mass of relativistic stars made of an incompressible fluid with a small electric charge. The investigation is carried out by using the hydrostatic equilibrium equation, i.e., the Tolman-Oppenheimer-Volkoff (TOV) equation, together with the other equations of structure, with the further hypothesis that the charge distribution is proportional to the energy density. The approach relies on Volkoff and Misner's method to solve the TOV equation. For zero charge one gets the interior Schwarzschild limit, and supposing incompressible boson or fermion matter with constituents with masses of the order of the neutron mass one gets that the maximum mass is the Misner mass. For a small electric charge, our analytical approximating scheme valid in first order in the star's electric charge, shows that the maximum mass increases relatively to the uncharged case, whereas the minimum possible radius decreases, an expected effect since the new field is repulsive aiding the pressure to sustain the star against gravitational collapse.

Abstract. One of the stiffest equations of state for matter in a compact star is constant energy density and this generates the interior Schwarzschild radius to mass relation and the Misner maximum mass for relativistic compact stars. If dark matter populates the interior of stars, and this matter is supersymmetric or of some other type, some of it possessing a tiny electric charge, there is the possibility that highly compact stars can trap a small but non-negligible electric charge. In this case the radius to mass relation for such compact stars should get modifications. We use an analytical scheme to investigate the limiting radius to mass relation and the maximum mass of relativistic stars made of an incompressible fluid with a small electric charge. The investigation is carried out by using the hydrostatic equilibrium equation, i.e., the Tolman-Oppenheimer-Volkoff (TOV) equation, together with the other equations of structure, with the further hypothesis that the charge distribution is proportional to the energy density. The approach relies on Volkoff and Misner's method to solve the TOV equation. For zero charge one gets the interior Schwarzschild limit, and supposing incompressible boson or fermion matter with constituents with masses of the order of the neutron mass one gets that the maximum mass is the Misner mass. For a small electric charge, our analytical approximating scheme valid in first order in the star's electric charge, shows that the maximum mass increases relatively to the uncharged case, whereas the minimum possible radius decreases, an expected effect since the new field is repulsive aiding the pressure to sustain the star against gravitational collapse.

Introduction
Compact stars and their properties have been a theme of great relevance on several grounds. Chandrasekhar's celebrated work [1] on the maximum mass for white dwarfs advanced the way to the understanding of the nature and structure of compact stars. By using a cold equation of state in which the degeneracy electron pressure is the most relevant form of pressure for the support of a white dwarf against gravitational collapse, a radius-mass relation for these stars was deduced, from the non-relativistic electron regime in relatively large white dwarfs up to the relativistic electron regime in the most compact stars. He found that as the radius of the star approached zero the mass would go to a maximum value of 1.44 M ⊙ , where M ⊙ is the sun's mass. This is the Chandrasekhar limit. It uses Newtonian gravitation. Landau [2] through heuristic arguments found that the mass limit for white dwarfs could be written as M ∼ M 3 pl /m 2 n where m n is the neutron or the proton mass and M pl is the Planck mass, M pl = c/G, with being the Planck constant, G the Newton's constant of gravitation, and c the velocity of light, or, setting units such that G = 1 and c = 1, which we do from now on, one has M pl = √ . Putting in the numerical values for M pl and m n , the mass M of the star is about the Chandrasekhar mass limit M ∼ 1 M ⊙ . Landau further deduced that the stars should have a radius of about λ e M pl /m n , where λ e is the electron's Compton wavelength, λ e = /m e , m e being the electron's mass, giving a radius of the order of 5000 km. He also found that there was another regime in which the star is composed of neutrons, supported by the degeneracy pressure of these particles, and has a maximum mass of about M ∼ M 3 pl /m 2 n . These neutron stars are much more compact with a radius λ n M pl /m n , λ n being the neutron's Compton wavelength, giving about 10 km. For objects with a radius tending to zero one should use general relativity, rather than Newtonian gravitation. In general relativity a compact star can be defined neatly as a star that has a geometrical mass M (or, GM/c 2 if one restores G and c) somehow comparable to its radius R, i.e., R/M ∼ a, with a a number not much bigger than 1, say of order of 10 or less. Whereas for an extended star like the Sun R/M ⊙ ∼ 5 × 10 5 , one has for a white dwarf R/M ⊙ ∼ 3 × 10 3 , and for a neutron star R/M ⊙ ∼ 6, showing that the latter is really compact. On a general relativistic basis, Oppenheimer and Volkoff [3] worked out further the nature and structure of neutron stars. Setting up a stiff equation of state for stars matter made of neutrons they found roughly the results of Landau [2], namely, the mass limit is about 1 M ⊙ and R/M ⊙ ∼ 6. This limit is called the Landau-Oppenheimer-Volkoff limit. Improvements have been made on these limits. Using a cold equation of state valid in the full range of highly compressed matter, the Harrison-Wheeler equation of state, the full set of equilibria in an R-M relation were found, in particular, the two maxima masses corresponding to the Chandrasekhar and Landau-Oppenheimer-Volkoff limits appear naturally [4]. See also [5] for further discussion on compact stars and [6] for issues of their stability. These mass limits, as seen in the context of general relativity, appear because at these stages the energy associated to the pressure is so strong that its gravitating effect overwhelms the self support effect.
The properties of these compact stars get modified if either the constituent material is altered or an alternative theory of gravitation is used. The constituent material can be altered by introducing or adding new matter with different properties or even adding a new field to the star. It is also well known that the introduction of new matter fields can be mimicked by modifications of the gravitational field. One example, even in Newtonian gravitation, is that the effects of the dark matter can be mimicked by modifying the gravitational field, e.g., in the MOND theories [7,8,9]. The same type of choice holds true in tensor theories of gravitation, since one can pass the excess of the gravitational field present in the left-hand side of the Einstein equations to their right hand side giving an effective energy-momentum tensor in a form of a new field, e.g., see [10]. If this dark matter inhabits the core of stars it can imprint onto compact stars new properties [11,12]. On the other hand, the structure of stars like neutron stars in alternative theories of gravity has been analyzed in an Einstein-Dilaton-Gauss-Bonnet gravity [13], in Eddington-like theories [14] with the conclusion that more massive stars than in general relativity can form, in a quadratic gravity theory [15,16] and in f (R) theories [17,18]. Stars and compact stars in braneworlds have also been studied in [19,20,21,22].
In contrast to the fermion stars mentioned above, there are stars made of bosons, the boson stars, that can have a wide range of mass limits, namely, where m b is the mass of the boson that makes up the star [23,24,25], see [26,27] for reviews. These stars could have been formed in the beginning of the universe from the primordial gravitational collapse of the boson particles and have been proposed as alternatives to the usual compact objects [27], such as neutron stars and black holes, and also as part of the dark matter [28]. These stars can, in principle, be detected [29,30,31] and a way to tell the the difference between fermionic and bosonic dark matter has been proposed in [32]. Now, the first compact star ever displayed in its full structure was a general relativistic star with a very stiff equation of state, a star made of an incompressible perfect fluid, i.e., ρ(r) = constant, and isotropic pressure (where ρ(r) is the energy density at the radius r) [33]. This interior Schwarzschild star solution is spherically symmetric and has a vacuum exterior. An incompressible equation of state is interesting from various aspects, since one can extract clean results and it also provides compactness limits. Furthermore, this incompressible fluid applies to both fermion and boson particles, as long as the fluid is at an incompressible state. As a drawback for such an equation of state, one can mention that the speed of sound through such a medium is infinite, but generically the overall structure is not majorly changed. Schwarzschild found that there was a limit, when the central pressure p c goes to infinity and that the star's radius to the mass limit is R/M = 9/4 = 2.25, indeed a very compact star [33]. Volkoff [34] and Misner [35] rederived the Schwarzschild interior limit of R/M = 9/4 using the propitious Tolman-Oppenheimer-Volkoff (TOV) equation, a differential equation for the pressure profile as a function of the other quantities [3], see also [36]. In addition, Misner [35] even found a maximum mass for a given density of the incompressible fluid, the Misner mass. The Schwarzschild limit yields, for the same mass, a radius that is well below the radius of a neutron star, and also yields the most possible compact boson star. An incompressible star whose radius is below this limit will in principle collapse into a singularity leaving a black hole to the exterior.
One can ask if the Schwarzschild limit can be modified, allowing for instance a lower R/M relation. As mentioned above, one way is to have some kind of different matter or a new field in the star, or an alternative theory of gravitation. Rather than introducing an alternative theory of gravitation we here opt to study the case in which we add a matter field to the existent matter. We consider matter with a small electric charge, introducing thus an additional electric field in addition to the usual matter and gravitational fields. This addition of an electric charge and an electric field to the Schwarzschild incompressible matter configurations brings insight to the configurations overall structure in more complex situations and its study in stars mimics other fields and possible alterations in the gravitational field.
Spheres of electric charge tend to be unstable. Since like charges repel each other, if in one way or another a charged cloud forms it will blow out at once due to the electrical repulsion. However, gravitation is universal and attractive. For a sphere with a small mass the electric repulsion is stronger than the gravitational pull and the matter suffers an expansion. On the other hand, for a highly massive sphere with a small quantity of charge the gravitational pull on the charged matter can overcome the electric repulsion and thus the system collapses. In between the two situations, one can have an equilibrium situation when the gravitational pull and the electric repulsion balance each other and the sphere stays static. In general, there is also matter pressure which adds to the balance of forces. For an extended star the pressure acts to counterbalance the gravitational field, but for highly compact stars the pressure acts as an energy source that adds to the star's mass and energy increasing in an overwhelming fashion the gravitational pull against its own pressure push. Thus, in the balance between the gravitational, electric and pressure fields, in highly compact stars electric charge can remain trapped within the stars. In case the compact stars supports some electric charge, their structure and properties are modified, and in particular, the radius-mass relations for the corresponding stars should change.
One important quantity that gives a measure of how much electric charge a star can support is the ratio of the mass m to the charge q of the main fundamental constituents of the star. For normal matter the net electric matter in a star is utterly negligible as the ratio of the proton mass m p to the proton charge e is m p /e = 10 −18 , giving thus Q/M ≃ (m p /e) 2 ≃ 10 −36 , where Q is the star's total charge [5] (see also [37]). However, stars can contain some dark matter in their interior, and of the several dark matter fluid candidates some could be electrically charged. Indeed, natural candidates to compose the dark matter are supersymmetric particles. The lightest supersymmetric particles that make the bulk of dark matter should be neutral, one possible candidate is the neutralino [38], however, some of these particles could be electrically charged. The mass m to charge q ratio of these supersymmetric particles are much higher than the baryonic mass to charge ratio, indeed current supergravity theories indicate that some particles can have a ratio of one. For a m/q ∼ 0.1 − 0.3 one has Q/M ≃ 0.01 − 0.1, a small but non-negligible electric charge. Thus, if dark matter populates the interior of stars, and some of it is made of electrically charged particles there is the possibility that stars have some electric charge.
That electric charge can influence the structure of a compact star was proposed earlier by Bekenstein [39] who wrote the appropriate TOV equation. Some electric compact configurations with an incompressible equation of state for the matter were studied numerically and the corresponding generalized Schwarzschild limit, i.e., central pressure going to infinity in these configurations, was analyzed [40,41,42,43]. Other equations of state for electrically charged matter, such as polytropic equations were used in [44,45,46], where star configurations and their structure were studied and the Schwarzschild electric limit for the given equation of state and for a given charge was considered. In particular, in [44] it was argued that upon gravitational collapse a star could retain, significantly, part of its electric charge. This electric charge excess could remain trapped in the final configuration, be it a highly compact star or a black hole. Other interesting equations of state were proposed and studied in [47,48], and in [49] where electrically charged dust was studied. Electrically charged boson stars have been also studied and their properties analyzed [50,51]. Bounds on the radius to mass relation for uncharged and charged stars have been put forward in [52,53], respectively, see also [54,55,56]. It is also worth mentioning some work on charged Newtonian stars. That compact stars could exist was even noticed in the 18th century. A subset of these Newtonian compact stars are the dark stars of Mitchell, also mentioned later by Laplace, see [57]. The Chandrasekhar white dwarfs of very small radius, including the one with zero radius that gives the Chandrasekhar maximum mass, are also Newtonian compact stars. Indeed the stars that have very small radii, from the gravitational radius to zero radius, provide instances of the dark stars of Mitchell and Laplace. Of course these stars cannot exist in nature as for strong gravitational fields Newtonian gravitation is invalid. That compact Newtonian stars could be electrified was raised in [58] and further developed in [59]. Turning the table around, the real analogue of dark stars are the quasiblack holes considered in, e.g., [48].
In order to understand the effects of a small electric charge on a star, and in particular, on the interior Schwarzschild limit, we use an analytical scheme and investigate the limiting radius to mass relation and the maximum mass of relativistic compact stars made of an electrically charged incompressible fluid. The investigation is carried out using the hydrostatic equilibrium equation, i.e., the TOV equation, and the other structure equations, with the further hypothesis that the charge distribution is proportional to the energy density. The approach relies on Volkoff and Misner's method [34,35] to solve the TOV equation. For zero charge one gets the interior Schwarzschild limit and under certain assumptions one gets the Misner mass. Our analysis for stars with a small electric charge shows that the maximum mass increases relatively to the uncharged case, whereas the minimum possible radius decreases, an expected effect since the new field is repulsive aiding the pressure to sustain the star against gravitational collapse.
The paper is organized as follows. In Sec. 2 we give the general relativistic equations, the equations of structure for a static spherically symmetric configuration, the equations of state for energy density and charge density, and discuss the exterior spacetime and the boundary conditions. In Sec. 3 we review the pure, uncharged, interior Schwarzschild limit using the Volkoff and Misner's formalism to set the nomenclature. We also give the Misner mass. In Sec. 4 we study analytically the interior electric Schwarzschild limit using the Volkoff and Misner's formalism and give the electric correction to the Misner mass. In Sec. 5 we conclude. In the appendix we derive some equations necessary in our study.

Basic general relativistic equations
We are interested in analyzing highly compacted charged spheres as described by the Einstein-Maxwell equations with charged matter. We set G = 1 and c = 1. The field equations are where Greek indices are spacetime indices running from 0 to 3, with 0 being a time index. The Einstein tensor G µν is defined as G µν = R µν − 1 2 g µν R, where R µν is the Ricci tensor R µν , g µν is the metric tensor, and R the Ricci scalar. The Faraday-Maxwell tensor F µν is defined in terms of an electromagnetic four-potential (1) is the Einstein equation, stating the relation between the Einstein tensor and the energy-momentum tensor T µν . T µν is written here as a sum of two terms, E µν is the electromagnetic energy-momentum tensor, which is given in terms of the Faraday-Maxwell tensor F µν by the relation M µν represents the matter energy-momentum tensor and we assume to be the energymomentum tensor of a perfect fluid, namely, with ρ and p being the energy density and the pressure of the fluid, respectively, and U µ is the fluid four-velocity. Equation (2) is the Maxwell equation, stating the proportionality between the covariant derivative ∇ ν of the Faraday-Maxwell tensor F µν and the electromagnetic four-current J µ . For a charged fluid, this current is given in terms of the electric charge density ρ e by The other Maxwell equation ∇ [α F βγ] = 0, where [...] means antisymmetrization, is automatically satisfied.

Equations of structure
The line element for a static spherically symmetric spacetime is of the form where t, r, θ e φ are the usual Schwarzschild-like coordinates, and the metric potentials A(r) and B(r) are functions of the radial coordinate r only. The assumed spherical symmetry of the spacetime implies that the only nonzero components of a purely electrical Faraday-Maxwell tensor F µν are F tr and F rt with F tr = −F rt and where F tr is a function of the radial coordinate r alone, F tr = F tr (r). The other components of F µν are identically zero. It is advantageous to define the total electric charge q(r) inside a spherical surface labeled by the radial coordinate whose value is r by I.e., one swaps F tr for q(r). It is also opportune to define a new quantity m(r) in such a way that 1 I.e., one swaps A(r) for m(r). The new function m(r) represents the gravitational mass inside the sphere of radial coordinate r.
One of the Einstein equations can be substituted by the contracted Bianchi identity ∇ µ T µν = 0, which gives dB(r) dr = B(r) p(r) + ρ(r) a differential equation for B, q, and p. Einstein equations also give a differential equation for B(r) alone, i.e., (11) Now, we are ready to write the other three equations in a form we want to use. One finds that another of Einstein equations gives a differential equation for m(r), i.e., dm(r) dr = 4πρ(r)r 2 + q(r) r dq(r) dr .
Since m(r) represents the gravitational mass inside the sphere of radial coordinate r, Eq. (12) represents then the energy conservation as measured in the star's frame. The only non-vanishing component of the Maxwell equation (2) is given by Finally, replacing Eq. (13) and the conservation equation (10) into Eq. (11) it yields dp dr = − (p + ρ) (4πpr + m/r 2 − q 2 /r 3 ) (1 − 2m/r + q 2 /r 2 ) + ρ e q/r 2 where to simplify the notation we have dropped the functional dependence, i.e., m(r) = m, q(r) = q, ρ(r) = ρ, p(r) = p, and ρ e (r) = ρ e . Eq. (14) is the TOV equation modified by the inclusion of electric charge [39] (see also [46]). The system of equations (11)-(14) is the system we were looking for. We need now to specify the equation of state and the equation for the charge density profile.

Equation of state and the charge density profile
In the present model there are six unknown functions: B(r), m(r), q(r), ρ(r), p(r), and ρ e (r); and just four equations: Eqs. (11), (12), (13), and (14). Additional relations are obtained from a model for the cold fluid, which furnishes relations among the pressure and the energy density. For an electrically charged fluid, a relation defining the electric charge distribution is also needed.
Here we assume an incompressible fluid, i.e., So the energy density is constant along the whole star. Following [44] (see also [46]), we assume a charge density proportional to the energy density, where, in geometric units, α is a dimensionless constant which we call the charge fraction. The charge density along the whole star is thus constant as well. Other equations for the charge distribution could be considered, as more charge concentration on the core, or more charge on the outer layers, see , e.g., [40,41,42].
We have now four equations: Eqs. (11), (12), (13), and (14); and four unknowns: B(r), m(r), q(r), and p(r), as ρ and ρ e are given in (15) and (16), respectively. The resulting set of equations constitute the complete set of structure equations which, with some appropriate boundary conditions, can be solved simultaneously. We are not going to solve it (see [43]). Here we use this system of equations to find the Schwarzschild interior limit for the small charge case.

The exterior vacuum region to the star and the boundary conditions
The conditions at the center of the star are that m(r = 0) = 0, q(r = 0) = 0, and A(r = 0) = 1 to avoid any type of singularities, and that p(r = 0) = p c , ρ(r = 0) = ρ c , and ρ e (r = 0) = ρ ec , where p c is the central pressure, ρ c is the central energy density, and ρ ec is the central charge distribution, the two latter having the same constant values throughout the star (see Eqs. (15) and 16).
The interior solution is matched at the surface to the exterior Reissner-Nordström spacetime, with metric given by where An important quantity for the exterior metric is the gravitational or horizon radius r + of the configuration. The Reissner-Nordström metric, given through Eqs. (17)-(18), has as the solution for its own gravitational radius.

Equations
Before we treat the small charge case analytically, we consider the exact Schwarzschild interior solution as given by Volkoff [34] and displayed later in Misner's lectures [35]. For this we put q = 0 in Eqs. (11)- (14). Equation (11) is of no direct interest here, Eq. (12) gives Eq. (13) is trivially satisfied in this case, and finally, Eq. (14) simplifies to dp dr Since, by equation (15), the density is constant we can integrate equation (20) obtaining where R is the star radius and we have imposed that there is no point mass in the center. Defining a characteristic length R c as we can rewrite the mass function, Eq. (22), as Interchanging ρ and R c as necessary and noting that 2 ρ R 2 c = 3 4π we get from Eq. (21), dp (25)

The interior Schwarzschild limit: The R and M relation and the minimum radius
Equation (25) is separable and can be integrated as with the boundary condition that the surface of the star R has zero pressure, i.e., Defining a new radial coordinate χ by Eq. (26) can be put in the form subjected to the boundary condition where χ s is given through Integrating Eq. (29), subjected to the boundary condition (30), yields the pressure The central pressure, p c is the pressure computed at zero radius r = 0, i.e., χ = 0, so that This blows up, This is equivalent to Now, Eqs. (24) and (28) allow us to write where M ≡ m(R) is the star's total mass. Thus, Eqs. (35) and (36) yield Equation (37) is the Schwarzschild limit found in [33].

Misner mass bound
Following Misner [35] we can also display a mass bound. Equation (24) gives Eliminating R in Eqs. (37) and (38), and noting that p c ≤ ∞, one gets the mass bound To get a mass we have to have a density and thus an R c . One can make sense of a constant density if one takes it as the density at which matter is almost incompressible and the pressure throughout the star is very high. If the fluid is an ideal gas this happens when the particles have relativistic velocities of the order 1. For fermions this happens when the Fermi levels are near the rest mass m n of the fermions, neutrons say, while for bosons this means that the gas temperatures are of the order of the rest mass m b of the particles. This gives, for both fermions and bosons, a density of one particle per cubic Compton wavelength. I.e., for a particle with mass m and Compton wavelength λ given by λ = /m the density is ρ ∼ m 4 / 3 . In the case of a star composed of neutrons, Misner [35] obtains where M ⊙ is the Sun's mass. This bound is similar to the Chandrasekhar limit  (14), treating the charge q(r) as a small perturbation, thus assuming α small. To do so, we note that the solutions for the mass and the charge will be of the form where we are assuming that the non-perturbed charge is zero q 0 (r) = 0, m 0 (r) is the mass of the uncharged star given by Eq. (20), or Eq. (22), and q 1 (r) and m 1 (r) are the perturbed small charge and mass functions to be determined. The pressure is also assumed to be given by the expansion where p 0 is the pressure in the uncharged case, given by equation (25), or (32), and p 1 (r) is the perturbation induced in the pressure when a small charge is considered. Note that, while the boundary condition for the non-charged star was simply p(R) = p 0 (R) = 0, the boundary condition for the charged star becomes At this point, it will prove useful to introduce the dimensionless variable where R c is the characteristic length defined in Eq. (23). The expressions for the mass, charge, and pressure in this new variable are generically defined as From Eq. (45) we defined x s as the x at the surface, so that Accordingly, we define 4.1.2. Calculation of the perturbed charge distribution q 1 Expanding Eq. (13) for small α, we get in the x variable that up to first order in α. Solving the above equation subject to the condition q 1 (0) = 0 and expressing the solution in terms of the variable x, results in

Calculation of non-perturbed and perturbed masses
The unperturbed mass m 0 can now be expressed simply as One can also find an expression for m 1 , namely, Indeed, from Eqs. (12) and (42), it is clear that the equation for the perturbed mass m 1 is given by Inserting Eq. (50) into Eq. (53), we can integrate it using the boundary condition that the total mass at the center of the star is m(0) = 0, which implies that m 1 (0) = 0 since Eq. (22), or Eq. (51), satisfies m 0 (0) = 0. Doing this, we are led to Eq. (52).

Equations for the pressures, solution for the zeroth order pressure, and calculation of the perturbed pressure at the star's radius (i) Equations for the pressures
To find the equations for the pressures p 0 (x) and p 1 (x), we begin by expressing Eq. (14) for the total pressure in terms of the variable x given in Eq. (45), Now we can expanded the right side of the above equation in powers of α and retain the two lowest terms. By doing so, and using the expansion (43) on the left side of Eq. (54) and Eqs. (41)-(42) on the right hand side, we can equate the terms in equal powers of α, thus obtaining two differential equations. The first one, obtained from the 0th power in α is dp 0 ( which is simply the differential equation for the unperturbed pressure. The second differential equation, to first order in α 2 , is which is the differential equation for the perturbed pressure p 1 , where, again, to shorten equations we have dropped the dependence of variables p 1 , p 0 , m 1 , m 0 , and q 1 on x, and we have also defined the auxiliary functions f 1 = f 1 (x) and f 2 = f 2 (x) by and Ultimately, we want to obtain an equation for the radius R for which the central pressure blows up. From Eq. (43), the central pressure is p(0) = p 0 (0) + p 1 (0). In the Appendix we show that p 1 (0) is always finite. So we have to find a solution for the radius R at which p 0 (0) blows up.
(ii) Solution for the zeroth order pressure We start by obtaining a solution for p 0 . Since the boundary condition has changed relatively to the uncharged case, it is now given by Eq. (44), we cannot use a priori the form (32) for p 0 . We use the solution to Eq. (55) without specifying any boundary condition. In the variable x, this means where A is an integration constant. To find out what this constant is, we insert the above equation into Eq. (44) and solve the resulting equation with respect to A. After expanding for the small charge parameter α, and so for small p 1 , we are led to up to first order in the perturbed quantities. Then, the expression for p 0 , analogous to Eq. (32), becomes Since p 0 (x) depends on p 1 (x s ) in the denominator, we have to find p 1 (x s ), i.e., we have to calculate the perturbed pressure at the star's radius.

(iii) Calculation of the perturbed pressure at the star's radius
The equation for p 1 , Eq. (56), cannot be solved analytically for all x. However, we are only interested in the value of p 1 at the surface of the star. At this particular radius it is possible to obtain the exact value of the perturbed pressure without ever solving Eq. (56).
The reason for this is the fact that at the star's radius the pressure p(x s ) = p 0 (x s )+p 1 (x s ) is zero. Therefore, we can insert the boundary condition p(x s ) = 0 in the exact derivative of the pressure given by Eq. (54) and expand the resulting equation for small α, giving up to first order, and where the auxiliary functions f 1 (x) and f 2 (x) are given by Eqs. (57) and (58), respectively. Using the expansion (43) on the left side of Eq. (62), one can clearly see that there is a compatibility condition which must be physically required, namely that, at the star surface, the first term on the right side of Eq. (62) must be equal to Eq. (55) and the second term equal to Eq. (56). Hence, we arrive at the two equations 3 2 These equations give two different solutions, namely, p 0 (x s ) = 0 and p 1 (x s ) = 0, or p 0 (x s ) = − 4 3 and p 1 (x s ) = − 2 3 x −1 s f 2 (x s ). This latter solution does not satisfy the boundary condition (44) so the unique valid solution is given by and

Equation for the minimum radius
We see that the central pressure p 0 (0) given in Eq. (61) is divergent when the following condition holds, Expanding it in α 2 we arrive to valid in first order in α 2 . Using the expansion provided by Eq. (42), it can be shown that to first order in α 2 , we have Upon substituting Eq. (69) into Eq. (68) we conclude that Now, the minimum star radius R will not be just 8/9 R c but will have corrections of order α 2 . These corrections will induce changes of the order α 4 in Eq. (70). Thus, we can set x s = 8/9 in Eq. (70), i.e., up to order α 2 .
In converting from the variable x back to r, we use In addition, it can also be interesting to express α in terms of the total charge Q and mass M. The following relation valid for small q 1 , or small α, can be found up to order q 1 . We can then express α in terms of the ratio Q/M. Using Eqs. (50) and (51) in Eq. (74), and solving the resulting equation for α, we obtain

The electric interior Schwarzschild limit: The R, M and Q relation for small charge
We are now in a position to calculate the desired ratio (72) in terms of the quantities R, M and Q and find the appropriate relation. Inserting Eq. (73) into Eq. (72) we find which is one form of the interior Schwarzschild limit for small charge. Inserting Eq. (75) into Eq. (76) we get valid up to order Q 2 /M 2 . This is another form of the interior Schwarzschild limit for small charge. We can also express the limit in terms of the horizon radius, r + , for the Reissner-Nordström metric. The horizon radius is defined by Eq. (19), i.e., up to order Q 2 /M 2 one has, up to order Q 2 /M 2 . Equation (78) is yet another form of the interior Schwarzschild limit for small charge. The electric interior Schwarzschild limit for small charge presented in various forms in Eqs. (76), (77), and (78) is the main result of this work. All the three forms of the electric interior Schwarzschild limit for small charge show that, in comparison with the uncharged case Eq. (37), the star can be more compact. In particular, Eq. (78) shows that in the charged case the radius of the star can be a bit nearer its own horizon.
In [43] these compact stars were studied numerically. An R/M × Q/M relation was given numerically for 0 ≤ Q/M ≤ 1. For small charge, Q/M ≪ 1, one can extract from the numerical calculations in [43] This should be compared to our analytical calculation valid in first order of Q/M, given here in Eq. (77), i.e., R M = 9 4 − 0.568 Q 2 M 2 . It shows that the numerical code used in [43] is compatible with the analytical calculation. In that work [43] it was also shown numerically that in the other extreme, namely, Q/M = 1, one would obtain a star at its own gravitational radius, R/M = 1, i.e., an (extremal) quasiblack hole.
A related theme is the Buchdahl and the Buchdahl-Andréasson bounds. Buchdahl [52] by imposing a simple set of assumptions, namely, the spacetime is spherically symmetric, the star is made of a perfect fluid, and the density is a nonincreasing function of the radius, found that the radius to mass relation is R/M ≥ 9/4. Thus the Schwarzschild limit [33], i.e., the limiting R/M configuration when the central pressure goes to infinity, is an instance that saturates the Buchdahl bound. Following the line of reason of Buchdahl, Andréasson [53] obtained a bound for the minimum radius of a star using the following energy condition p + 2p T ≤ ρ where p T is the tangential pressure, p is the radial pressure and ρ is the energy density. This bound, the Buchdahl-Andréasson

Conclusions
In this work we have studied compact stars with a small electric charge and found the limiting radius to mass relation and the maximum mass through an analytical approach based on Volkoff and Misner's method to solve the TOV equation for incompressible matter. More specifically, we have analyzed the interior Schwarzschild limit of spherically symmetric star configurations composed of a fluid with constant energy density ρ and with a small electrical charge distribution ρ e proportional to ρ, ρ e = αρ with α ≪ 1. The exterior spacetime is described by the Reissner-Nordström metric. We have found through our analytical scheme that due to the electric charge distribution the limiting star configuration can have more mass and a smaller radius relatively to the limiting star with zero charge. This is expected since the electric charge distribution has a repulsive effect, adding to the pressure as a force that withstands the star. For stars containing some type of dark matter in their interior there is the possibility that they possess a small but non-negligible electric charge, in which case our analytical formula is opportune and can be confronted quantitatively with observational data.
Appendix A. Behaviour of the perturbed pressure for small radius Our goal is to compute the limit of infinite central pressure p(r = 0) in this charged case, so it is important to analyze the behaviour of p 1 for small radius x, x ∼ 0, which also means χ ∼ 0, and r ∼ 0. We start by obtaining the perturbed charge, q 1 , for small radius expanding equation (50), which gives The non-perturbed pressure, p 0 , given by equation (61) becomes Using equation (56) for the perturbed pressure, p 1 , together with equations (A.1), (A.2) and (61) it gives dp 1 dx ∼ −p 1 (3p 0 (x = 0) + 2) x + 1 2 α 2 x . So the unique way that the central pressure p(r) = p 0 (r) + p 1 (r) blows up is p 0 (r = 0) blowing up.