An electromagnetic extension of the Schwarzschild interior solution and the corresponding Buchdahl limit

We wish to construct a model for charged star as a generalization of the uniform density Schwarzschild interior solution. We employ the Vaidya and Tikekar ansatz [{\it Astrophys. Astron.} {\bf 3} (1982) 325] for one of the metric potentials and electric field is chosen in such a way that when it is switched off the metric reduces to the Schwarzschild. This relates charge distribution to the Vaidya-Tikekar parameter, $k$, indicating deviation form sphericity of three dimensional space when embedded into four dimensional Euclidean space. The model is examined against all the physical conditions required for a relativistic charged fluid sphere as an interior to a charged star. We also obtain and discuss charged analogue of the Buchdahl compactness bound.


I. INTRODUCTION
In 1916, almost immediately after the derivation of the unique exact solution describing the exterior gravitational field of a static spherically symmetric isolated object, Schwarzschild [1,2] obtained an interior solution for a uniform density fluid sphere of finite radius whose exterior was described by the former metric. Since then, a large number of exact solutions have been obtained for a more realistic description of relativistic compact stars, out of which only a few are shown to be physically viable, regular and well-behaved [3]. A natural extension of such models has been the inclusion of the electromagnetic field. The corresponding Einstein-Maxwell system, being highly non-linear, is extremely difficult to solve, and hence different simplifying techniques are often invoked to solve such a system. The choice of the Vaidya and Tikekar (VT) [4] metric ansatz, is one such approach which, ever since its inception, has found huge success for modelling of realistic astrophysical objects. The VT model was subsequently generalized by many investigators [5][6][7][8][9]. This ansatz is motivated by a geometric property that t = const. hypersurface of the associated spacetime, when embedded in a 4-Euclidean space is not spherical but spheroidal. The parameter k, which appears in the ansatz, indicates the departure from the sphericity of associated 3-space. The VT ansatz has been utilized to model a wide variety of compact stars (see e.g., Ref. [10][11][12][13][14][15][16][17]) and radiating stars (see e.g., Ref. [18][19][20][21]). The model has found its applications in higher dimension spacetimes as well (see, e.g., Ref. [22]). Recently an anisotropic stellar model has also been developed in the Buchdahl-Vaidya-Tikekar metric ansatz [23].
In this paper, we revisit the Vaidya-Tikekar metric ansatz to model a static spherically symmetric compact star with a charged fluid interior. Incorporation of the electromagnetic field in the modelling of relativistic astrophysical objects has a long history and is relevant for a wide variety of astrophysical systems at different stages of its evolution. Some of the pioneering works in this field include the investigations of Majumdar [24], De and Raychaudhari [25], Papapetrou [26], Cooperstock and Cruz [27] and Bekenstein [28], to name a few. Among many other factors, the observed high value of the electric field at the surface of ultra-compact stars [29] provides a sufficient ground to study the Einstein-Maxwell system. A large class of analytic solutions to the Einstein-Maxwell system corresponding to the exterior Reissner-Nordström metric has been compiled by Ivanov [30].
There has been considerable work in the literature on charged fluid model with the VT ansatz and with different prescriptions for the electric field, and models have been examined for their physical viability and acceptability [14,[31][32][33][34]. To get a solution of the Einstein-Maxwell field equations, one has always to prescribe some fall off behaviour for electric field and/or choose a specific value of the spheroidal parameter k. All these papers and the references given therein offer a good spectrum of different choices. In this paper, our main motivation is to superimpose electric field distribution on uniform density fluid distribution. We demand that when the electric field is switched off, the solution should reduce to the uniform density Schwarzschild solution. This limit would require the parameter k → 0, which means electric charge should be proportional to geometric parameter k. This, of course, amounts to prescribing a particular charge and electric field distribution. The novel and interesting feature of this model is clubbing of k with charge distribution. It is matched to R-N metric at the boundary defined by p = 0, and the solution is physically fully viable and satisfactory for modelling an astrophysical charged stellar object.
We obtain and discuss charged analogue of the Buchdahl limit [35] for the model. In relativistic astrophysics, measurement of the maximum mass to radius ratio of a stellar configuration has been a matter of prime importance ever since Buchdahl [35] derived the relation M/R ≤ 4/9 for an isotropic fluid sphere. In a very simple manner the Buchdahl bound could be found by demanding pressure at the centre being finite. This is what has motivated us to seek charged analogue of the uniform density Schwarzschild fluid sphere so that we could obtain the charged analogue of the Buchdahl limit. Theoretical developments in the analysis of Buchdahl limit for Schwarzschild as well as Reissner-Nordström background spacetimes is available in Ref. [36] and references therein. In this paper, by restricting the model parameters in such a manner that central pressure does not diverge, we obtain an upper bound on compactness of charged star which is a charged generalization of the Buchdahl limit.
The paper is organized as follows. In section II, for a static and spherically charged fluid distribution, we lay down the independent set of equations for the Einstein-Maxwell system in terms of a single generating function f (r). The idea of this generating function is such that by setting f (r) = 0, one gets back the Schwarzschild solution for a homogeneous fluid sphere. For f (r) = 0, we specify the charge distribution q(r) so as to have fallen back on the Schwarzschild solution. The most interesting choice we would explore is f (r) = kr 2 /C 2 , which is, in fact, the VT-metric ansatz, and we generate a new class of solutions where the electromagnetic field gets determined in terms of the parameter k. In section III, we match the interior solution to the exterior Reissner-Nordström metric at the boundary given by p = 0. Making use of the junction conditions, we determine the constants of the model in terms of total mass M , radius R and charge Q. In section IV, we discuss the Buchdahl limit, some physical properties and overall viability of the model. We end with a discussion in section V.

II. EINSTEIN-MAXWELL SYSTEM
We write the line-element describing the interior of a static spherically symmetric charged fluid sphere in the form in standard coordinates x i = (t, r, θ, φ) where ν(r) and µ(r) are the undetermined functions. The unknown functions can be determined by solving the Einstein-Maxwell field equations where, represent the energy-momentum tensor corresponding to matter and electromagnetic field, respectively. In Eqs. (4) and (5), ρ, p, σ and F ij denote the energy-density, pressure, charge-density and electromagnetic field tensor, respectively. u i is the 4-velocity of the fluid. Due to spherical symmetry, the only surviving independent component of the electromagnetic field tensor is F tr E. The Maxwell equations (3) yield where the total charge q(r) contained within the sphere of radius r is defined as From Eq. (6), we have The Einstein-Maxwell field equations (in system of units having G = c = 1) are then obtained as where, a prime ( ) denotes differentiation with respect to the radial parameter r. Combining equations (10) and (11), we obtain which becomes the definition of the electric field provided µ(r) and ν(r) are known. To solve the system, we assume the metric potential µ(r) in the Buchdahl-VT [4,35] form as where C is an arbitrary constant. It is not yet the Schwarzschild solution because ν is still undetermined. Note that when f (r) = 0, one can generate the Schwarzschid interior solution for an incompressible fluid sphere. A coordinate transformation x 2 = 1 − r 2 C 2 allows us to rewrite Eq. (13) in the form where f x represents derivative with respect to x. To solve equation (15), we introduce a new variable and rewrite equation (15) as Now let's try to superimpose charge distribution on the Schwarzschild solution. This suggests that ψ(r) should be a linear function and when charge distribution is switched off; i.e. f = 0. Thus, we demand d 2 ψ/dx 2 = 0 giving us and Note that the above solution is obtained for the particular choice (18) which in terms of the radial parameter r takes the form The choice of the charge distribution q(r) is, therefore, motivated by the fact that it provides a simple solution which may be treated as a generalization of the Schwarzschild interior solution for an incompressible fluid as discussed below. Moreover, the choice ensures that q(r) is well behaved at r = 0 as well as at all interior points of the star. Combining equations (16) and (19), the unknown metric function e ν is obtained as where, a, b, C, and k are constants to be determined by matching the solution to the exterior R-N metric at the boundary. Thus, the space-time metric of a static and spherically symmetric object in the presence of an electric field is obtained as Note that it reduces to the Schwarzschild solution when f (r) = 0, and it is f (r) that determines charge q(r) 2 in equation (20). The model would be fully determined when f (r) is prescribed. In the next section we shall now consider the particular choices of this function for building a model for a charged stellar object.
A. Case f (r) = 0 The electric field in this case vanishes and we regain the Schwarzschild interior solution for an (uncharged) incom- where R is the radius of the object and density and pressure are as given below That pressure should vanish at the boundary determines the constant R = C 1 − a 2 9b 2 which implies that we must have a < 3b. The central pressure takes the form which implies that we must have a > b. Combining the two conditions we get the bound, b < a < 3b.
We must choose the function f (r) such that it provides a regular, well-behaved and physically meaningful stellar model. Accordingly, we write f (r) = k r 2 C 2 which is the Vaidya-Tikekar (VT) [4] ansatz for the modelling of a relativistic compact star. The 3-hypersurface is spheroidal which becomes flat and sphere for k = −1, 0, respectively. The spacetime is well behaved for r < C and k > −1.
The energy-density, pressure and charge-density are then given by and at the centre they take the values as Obviously, the central density remains positive for k > −1 and pressure for a > b > a/3. For a = b, the central pressure diverges. The total mass within a sphere of radius r defined as integrates to give m(r) = 1 16 It is noteworthy that m(r = 0) = 0. Using Eq. (20), we have which also vanishes at the centre.

III. MATCHING CONDITIONS
The exterior spacetime of the static charged object is described by the Reissner-Nordström metric where M and Q represent the total mass and charge, respectively. The matching conditions are the continuity of e ν , e µ , and p = 0 at the boundary r = R. This means m(R) = M and q(R) = Q, and so we write The conditions explicitly take the following form at r = R, where, n = R 2 C 2 , u = M R and α 2 = Q 2 M 2 . Solving Eq. (40), we get so that where y = 1 − 2u + α 2 u 2 . Solving Eqs. (39) and (41), we obtain a = y[24 + 54kn + 4k 3 n 2 + k 2 n(9 + 25n)] Thus, all the constants are expressed in terms of k, M , R and Q. We rewrite Eq. (38) as and substitute the value of n to obtain where g(u, α) = 16 − 96α 2 + 288uα 2 + 96uα 4 − 456u 2 α 4 − 24u 2 α 6 + 232u 3 α 6 − 39u 4 α 8 .
Note that for α = 0, we have k = 0 while the converse is always true from Eq (39).

A. Application to astrophysical objects
A physically acceptable stellar interior solution should have the following features: (i) The density and pressure should be positive throughout the interior of the star i.e., ρ, p > 0; (ii) the pressure p should vanish at some finite radial distance i.e., p(r = R) = 0 and (iii) the causality condition should be satisfied throughout the star which implies that 0 ≤ dp dρ ≤ 1. To verify whether the above conditions are fulfilled in this model, we take the mass and radius of the pulsar 4U 1820 − 30 as input parameters. The estimated mass and radius of the star are M = 1.58 M and R = 9.1 km, respectively [37]. With these values the values of the constants, for different choices of the parameter k, are given in Table I It is noted that the parameter C, which goes inverse to energy density increases with increasing k; i.e. density decreases. Since k is directly related to charge, which means an increase in k means increase in charge. Subsequently, increase in repulsive component due to charge resulting in a decrease in fluid density. For these set of values, we show the behaviour of the physically interesting quantities in Figs. (1)-(7). The plots indicate that the model is regular and well-behaved at all interior points of the star. Figs (1) and (2) respectively show that density and pressure monotonically decrease with increasing radius. Note that central density is larger for larger k, in contrast central pressure in larger for lower k. That is central density and pressure are respectively largest and smallest for homogeneous fluid didtribution. The rate of fall is stronger for larger k, again reflecting the repulsive effect of charge. It is interesting to note that at r ∼ 6 km, all density curves cross the uniform density straight line downwards. This clearly indicates that homogeneous distribution has the largest density at the boundary. On the other hand, mass and charge as expected monotonically increase with the radius as shown in Figs. (3)-(4). Like density curves, mass curves also cross over at some r, and the rate of increase for larger k falls down the uniform density curve. On the other hand, charge always has a stronger rate of increase with increasing k. Fig. (5) shows that the electric field is zero at the centre, and it monotonically increases towards the boundary and the rate of growth is stronger for larger k. Radial variation of the charge density is shown in Fig. (6). Note that matter density decreases with r while charge density increases and attains maximum then slowly decreases. It is interesting to note that the charge density becomes maximum at a radial distance where the inhomogeneous density meets the uniform density profile. The sound speed for different k values are shown in Fig. (7) which indicates that even though for relatively higher values of k the causality condition is satisfied, as we approach the constant density case (k ∼ 0.1), sound speed becomes as expected greater than unity.

IV. BUCHDAHL LIMIT
Astrophysically it is of prime importance to find how compact a star could be, i.e. what is the upper bound on mass to radius ratio -the compactness ratio, M/R? From an intuitive perspective, the stiffest equation of state is of uniform density incompressible fluid which is uniquely described by the Schwarzschild interior solution. The compactness limit would be indicated by upper bound on pressure at the centre. As we have seen above, this would be the condition a ≥ b in Eq. (31) giving the compactness limit, M/R ≤ 4/9. Buchdahl was first to obtain this limit under very general conditions of density and pressure being positive, the former decreasing with radius outward and at the boundary, it is matched to Schwarzschild exterior solution [35]. The same limit is also obtained for anisotropic fluid by invoking the strong energy condition, ρ ≥ p r + 2p t [38,39].
For a charged object there exist more than one limits [40][41][42][43] obtained for different interior distributions and equation of state. One that is closest in spirit to the uniform density case [41] in which it is envisaged a uniform density distribution is enveloped by a thin charged shell and the Buchdahl analogue limit reads as (1 + 1 − 8α 2 9 ) , It reduces to the Buchdahl limit, M/R ≤ 4/9 when charge is switched off Q = 0, and M/R ≤ 2/3, 8/9 < 1 for α 2 = 1, 9/8, respectively. Interestingly it prescribes the upper bound on charge a star could have, α 2 ≤ 9/8 which is > 1. That is, a non-black hole charged object could indeed be overcharged relative to a charged black hole. Very recently an insightful and novel prescription [36] has been proposed for the compactness limit which is given by gravitational field energy being less than or equal to half of the non-gravitational matter-energy of the object. The remarkable feature of this definition is that it is entirely determined by the unique exterior R-N metric without any reference to interior distribution, may what that be! The limit that follows is the one given above.
In the following, we would like to obtain the compactness ratio M/R for our model. For a homogeneous sphere, the mass within a sphere of radius r is given by where ρ 0 = 3 8πC 2 is the homogeneous density. Defining the compactness factor as we write the metric potentials of the inhomogeneous sphere in the form We also express the physical quantities as Equation (56) shows that charge in our model is proportional to the constant density mass m 0 . The dimensionless parameter k, representing deviation from sphericity, also gets tagged into the expression of charge.
Since the departure from sphericity is expected to be small, expanding tan −1 √ kr in equation (34), we obtain the mass function as m(r) = m 0 32(1 + k) + (41 + 53k)kφ 0 − 6(11 + 7k)k 2 φ 2 which shows that for k = 0 we regain m(r) = m 0 . At the centre r = 0, we have e 2ν = (a − b) 2 , e 2µ = 1, q 2 = 0, The regularity of σ c demands that we must have k ≤ 2. We also note that p c → ∞ for a = b. In other words, for a stellar configuration, we must have a > b and consequently, the upper bound on the compactness vis-a-vis Buchdahl type limit can be obtained by setting a = b. Now, imposing the condition a ≥ b in Eqs. (44) and (45), we obtain In the absence of any simple solution of Eq. (58), let us first consider the case k = 0. Substituting k = 0, in Eq. (38), we obtain Q = 0 which implies that we must set α = 0 in this case. Note that in the uncharged case (k = 0), the exterior Reissner-Nordström metric should be replaced by the Schwarzschild solution by setting Q = 0. Therefore, substituting k = 0 = α in Eq. (58), we obtain which readily provides the Buchdahl limit u(= M/R) ≤ 4/9. For k = 0, substituting the value of k from Eq.(47) into Eq. (58), and solving it numerically, one can find the upper bound on u = M R for a given charge to mass ratio, α 2 = Q 2 M 2 . The results are shown in Table II and Fig. (8). Alongside Eq. (58) we also plot Eq. (48) clearly indicating how beautifully our model is coasting along with the exact Buchdahl bound for charged object.

A. Approximation method
To reduce the complexity of the equations, we set k = k , where 0 < << 1. This is a reasonable assumption as the departure from sphericity denoted by k is expected to be small. In this case, retaining terms up to O( ), we use (46) to obtain Eq. (60) ensures that we have α = 0 for k = 0. Inserting the value of α 2 in Eq. (58) and neglecting terms O( 2 ), we obtain Even though this is an approximate equation, a neat solution of the above is not available. Nevertheless, numerical solution of the above equation provides a generalization of the Buchdahl limit in the case of a charged sphere. It is not difficult to show that for k = 0, the equation yields the Buchdahl limit u ≤ 4/9. To find an analytic solution of (58) in the extreme case a = b (when the central pressure diverges), we make use of equation (47) and obtain a truncated equation up to the order O( ) as which readily yields Rationalizing the above, we finally obtain a charged analogue of the Buchdahl limit given by .
The above result provides an upper bound on α 2 ≤ 9/8 and u ≤ 8/9 < 1. Non black hole object would always have radius larger than the black hole. For α 2 = 0, we regain u ≤ 4/9. It is remarkable to note that by making use of a different technique and a specific model, we have been able to obtain the desired upper bound on the compactness of a charged object -a charged Buchdahl limit. By employing the Vaidya-Tikekar metric ansatz we have constructed a generalization of the homogeneous density distribution in which the parameter k gets coupled to charge distribution. When charge is set to zero, the solution goes over to the Schwarzschild uniform density fluid sphere. That means it is kind of charging the Schwarzschild uniform density solution.
This has facilitated computation of the charged analogue of the Buchdahl compactness limit in Eq. (58). In Fig. 8, the compactness ratio M/R is plotted for Eqs. (48) and (58), the former is the Buchdahl bound as found in Refs. [41] and [36] while the latter is computed for the present model. It is remarkable to see that how the one due to Eq. (refmaseq) coasts beautifully that due to Eq. (48).
Another worth noting feature is that where charge density (Fig. 6) attains the maximum value where energy density crosses the uniform density line (Fig. 1). This indicates that charge density increases while energy density decreases with radius until the latter crosses the uniform density line. Then the former attains maximum value and begins decreasing. This means for radius greater than the one where energy density becomes less than the uniform density value, charge density also begins decreasing. It is interesting to see how the behaviour of charge density is linked to the fact whether energy density is greater or less than the uniform density value.
It is remarkable that a charged object could have α 2 = Q 2 /M 2 ≤ 9/8 > 1; i.e. it could be overcharged relative to a charged black hole. The question then arises, is it possible to construct models with 1 ≤ α 2 ≤ 9/8 -an explicit example of overcharged object? It would be interesting to construct such a model and that's what we would like to take up next in a separate investigation.