Charged Vaidya-Tikekar model for super compact star

In this work, we explore a class of compact charged spheres that have been tested against experimental and observational constraints with some known compact stars candidates. The study is performed by considering the self-gravitating, charged, isotropic fluids which is more pliability in solving the Einstein-Maxwell equations. In order to determine the interior geometry, we utilize the Vaidya-Tikekar geometry for the metric potential with Riessner-Nordstrom metric as an exterior solution. In this models, we determine constants after selecting some particular values of M and R, for the compact objects SAX J1808.4-3658, Her X-1 and 4U 1538-52. The most striking consequence is that hydrostatic equilibrium is maintained for different forces, and the situation is clarified by using the generalized Tolman-Oppenheimer-Volkoff (TOV) equation. In addition to this, we also present the energy conditions, speeds of sound and compactness of stars that are very much compatible to that for a physically acceptable stellar model. Arising solutions are also compared with graphical representations that provide strong evidences for more realistic and viable models, both at theoretical and astrophysical scale.


I. INTRODUCTION
After the discovery of general relativity in 1915 by Einstein, it's became an important tool for understanding and explaining the gravitational system. In particular, obtaining a singularity free interior solution for compact astrophysical objects is an important issue in relativistic astrophysics for the past two and a half decades. At some points, we usually refer compact objects as a collectively of different types of high density objects like white dwarfs, neutron stars and quark stars, that form at the end of their stellar evolution. Therefore, in order to study the structure of such stars form microscopic composition and properties of dense matter in an extreme condition is one of the most fundamental problems in modern astrophysics. However, in spite of the fact that such extreme densities nuclear matter may consist not only of leptons and nucleons but also consists of some mesons, hyperons, baryon resonances as well as strange quark matter (SQM), in their different forms and phases at the time of stellar evaluation. Indeed, works that reviewed by different authors demonstrate the difficulties associated with constructing solutions for obtaining a comprehensive description of compact objects.
From the above description one can easily understand that obtaining a reliable description of dense matter compact object is not an easy task, though there are several theoretical investigations, laboratory experiments as well as observations test have been performed during the last few decades. But the observational data form compact stars may soon provide some information about largest uncertainties in nuclear physics that rely heavily on the equations of state (EoS) at nuclear and supranuclear densities. In general, this could be achieved by estimating their mass and radius [2][3][4] which depends on EoS [5][6][7][8]. The motivation to undertake such a task was initiated by the discovery of the pulsar PSR J1614-2230 that produced by hot spots on the surface of compact stars [9] led to several interpretative problems such as neutron stars (NSs) and quark stars (QSs). However, this pulsar profile does not depend only on the stellar mass and radius, but also depends on several other features like moment of inertia, quadrupole moment and higher multi-pole moments etc. According to recent observation there are many compact objects, namely, X-ray burster 4U 1820-30, X-ray pulsar Her X-1, X-ray sources 4U 1728-34, PSR 0943+10 and RX J185635-3754, whose masses and radii are not compatible with the standard neutron star models. The problem is that we still have lack of information about the nuclear matter density, so theoretical studies hints that pressure is likely to be anisotropic within the stellar radii, i.e., the radial pressure and the tangential pressure. The search for anisotropic superdense stars were initiated by Bowers and Liang [10] and then a number of articles appeared related with this in [11][12][13][14][15][16][17][18][19].
The search for an exact solutions of Einstein field equations for static isotropic and anisotropic astrophysical objects are excellent testbeds with growing interest to mathematician as well as physicists. However, most of the exact interior solutions for both isotropic and anisotropic cases do not satisfy the general physical required conditions of the stellar systems. Therefore, exact solutions of Einstein-Maxwell field equations are also important in relativistic astrophysics. The conjecture is the obtained solution may be utilized to model for a charged relativistic star which match to the Reissner-Nordström exterior spacetime at the boundary. In this context one can avoid the gravitational collapse of a spherically symmetric distribution of matter, because of gravitational attraction is counterbalanced by the repulsive Coulombian force in addition to the pressure gradient. In this connection, Ivanov [46] and Sharma et al. [60] showed that due presence of the electric field affects the nature of luminosities, redshifts and maximum masses of relativistic star. Later on, Takisa and Maharaj [59] have proposed an exact solutions for charged anisotropic polytropic spheres. Imposing different equations of state some of these charged solutions have studied in [20][21][22][23][24][25].
The simplest known procedure that can be added for constructing a static charged perfect fluid interior solution for compact objects are either to prescribe a metric ansatz or an equation of state relating pressure and density. In the present article we consider a well known ansatz for one of the metric functions, namely, Vaidya and Tikekar [1], who prescribed an ansatz for the geometry of the t = constant hypersurface. The main important feature is that the solution characterizes a class of static spherically symmetric perfect fluid configuration and provides an exact solution of Einsteins equations. Following this technique, a large number of solutions have been studied in [26][27][28][29][30]. In a recent treatment Naveen and Bijalwan [31,32] have obtained a charged perfect fluid model with generalized electric intensity for all K except for 0 < K < 1 and extending this point of view Kumar and Gupta [33,34] obtained another solution for 0 < K < 1. In the present problem we focused on a charged fluid sphere starting with Vaidya and Tikekar [1] metric potential and tested our model with some standard observed mass and radius of the compact stars candidates as proposed in [35].
The paper is organized as follows: following a brief introduction in Sec. I, then we give the general relativistic formulation of Einstein-Maxwell system of equations for a relativistic stellar model in Sec. II. Paying particular attention to solve the system of equations analytically, we assume a particular form of metric potential, namely, Vaidya-Tikekar to generate exact solutions and obtain the expression for energy density and pressure in the same section. Then, we match the interior charged fluid to the exterior Reissner-Nordström line element in Sec. III. Next, in Sec. IV, we have discussed briefly some physical features of the proposed model maintaining the regularity conditions and obtained results are compared with observational data. Finally in Sec. V, we give a brief discussion.

II. EINSTEIN FIELD EQUATIONS
In this work, we consider the static spherically symmetric spacetime for seeking solution of a compact stellar object, that can be written in Schwarzschild coordinates as where the unknowns ν(r) and λ(r) are both metric functions in terms of radial coordinate, which yet to be determined by solving the field equations.
The most important question that arise is the matter distribution inside a compact star. Here, we consider the case of a charged gravitating object with isotropic pressures. Then the energy momentum tensor T i j will include the terms from the Maxwell's equation E i j and the complete form of energy-momentum tensor is where u j , and δ i j stand for 4-velocity of the fluid, and the metric tensor, respectively. Generally first component standards for energy momentum tensor; e.g. for a perfect fluid we have (c 2 ρ + p)u i u j − pδ i j , where ρ is the matter density and p is the pressure of the fluid. The basic argument to assume the perfect fluid implies that the flow of matter is adiabatic, no heat flow, radiation, or viscosity is present [36]. The second term associate with electromagnetic stress-energy tensor from the Maxwell's field equation and hence they will follow the relation √ −gF ij ,j = 4πj i √ −g, where F ij denote the skew symmetric electromagnetic field tensor and j i is the four-current density. Now, accomplishing the effects due to the electric field and pressure isotropy, and using the metric (1) with stress tensor given in Eq. (2), the Einstein field equation, Rδ i j , provides the following relationships with κ = 8π (geometrized units G = c = 1) and prime denotes the differentiation with respect to the radial coordinate. The total charge inside a radius r is given by Also Note that Eqs. (3)-(5) provide three independent equations, for five unknown quantities i.e. ν, λ, ρ(r), p(r) and q(r) which we have to solve simultaneously to get our results. As, obtaining an explicit solutions to the Einstein field equations is a difficult tusk due to highly nonlinearity of the equations. Thus, we will reduce the number of unknown functions by assuming a well known form of metric potential [1] where C = −K/R 2 , and K & R are two parameters which characterize the geometry of the star. Our choice of Vaidya-Tikekar ansatz is physically well motivated and has been studied for uncharged superdense stars by Tikekar [37] and Maharaj and Leach [38]. This facilitated the model in an interesting geometric meaning as deviation from sphericity of 3-space geometry. It may also be noted that metric potential restricts the geometry of the 3-dimensional hypersurfaces t = const. to be spheroidal and when K = 0 the hypersurfaces t = const. become spherical. The solution of EFEs (3)-(5), is in a different but equivalent form if we introduce Cr 2 = x and e ν = Z 2 as a new variables. After a little bit algebraic calculation we have where prime denote differentiation with respect to the variable x. Furthermore, to transform the field equations to a more convenient form we introduce another variables defined by and plugging the values of Y and Z into the equation (9), we get where for notational simplicity we use In order to solve the second order differential equation (12) more easily, we have chosen where θ is a positive constant, and comparing with Eq. (13) leads to defining the total charge of the system as Moreover, using the value of χ from (13) into the Eq. (12), which yields Now, comparing Eq. (15) with a known standard differential equation which leads to the following relations P 0 = Y 2 θ 2 +b, P 1 = 0, P 2 = −2θ 2 , R = 0. Now, it is clear that the differential equation is exact and rearranging the terms one can write this in a conventional form P 2 − dP 2 dY + d 2 P 0 dY 2 = 0. Hence the primitive of the given equation is where A and b stands for arbitrary constants. Finally, the differential equation (17) leads to the following expression Now, using the relation (11) for Z, which yields where In addition, taking into account Eqs. (14) and (19) into equations (8)-(10), it is straightforward to achieve the following nonzero components of field equations where, N 1 = 1  (Table-II (Table-III). .
With the purpose of determining the non-zero components, we are in a position to determine the internal structure of compact stars. For compact star described by electrically charged fluid (within a certain radius) should comply with the following requirements throughout the interior radius: • The energy density and pressure should be positive and regular within the radius.
The above three conditions imply that pressure and density should be maximum at the center and monotonically decreasing towards the surface, which is clear from Fig. (1). We are concerned here with charged isotropic case and according to Thirukkanesh and Mahara [61] charge distributions are singular at the origin, where the electric field does not vanish. That means, at the origin all our sources have vanishing electric field and finite proper charge density. From Fig. (1), it is clear that the charge will increase with the increasing charge fraction.

III. BOUNDARY CONDITIONS
At this stage the interior solution is smoothly connected to the vacuum exterior Reissner-Nordström metric at the junction surface with radius r = R. It is important to note that these three constants A, b and θ are fixed by suitable junction conditions imposed on the internal and external metrics at the hyper-surface. The exterior metric is given by where M is the total gravitational mass of the fluid distribution and is defined by with the definition ζ(R) = κ 2 R 0 ρr 2 dr, ξ(R) = κ 2 R 0 rσqe λ/2 dr and Q = q(R) represents, the mass within the sphere, the mass equivalence of the electromagnetic energy of distribution and Q is the total charge inside the sphere as suggested by [39]. To do the matching properly, at the boundary surface r = R, we start by imposing the junction condition that the metric should be continuous [40]. Thus, by joining the interior metric function g rr = e λ and g tt = e ν with the metric coefficient of the exterior Reissner-Nordström spacetime given in (22), we obtain the following conditions by using the continuity p(R) = 0, and q(R) = Q.
Now, using the conditions (24) and (25) Table -4. We will demonstrate here the interior and surface gravitational redshift z S of these compact sources by using the definition z S = ∆λ/λ e = λ0−λe λe , where λ e is the emitted wavelength at the surface of a nonrotating star and λ 0 is the observed wavelength received at radial coordinate r. Thus, one of these quantities is defined according to where g tt (r) = e ν(R) = 1 − 2M R + Q 2 R 2 is the metric function. Note that interior redshift should decrease with the increase of radius and less than the universal bounds, found when different energy conditions holds. In the isotropic case gravitational redshift for perfect fluid spheres is given by z s < 2 [41,42]. For an anisotropic star, this value admits higher redshifts, Z s = 3.84, as given in Ref. [43,44], but in the presence of cosmological constant provides a significant increase up to z s ≤ 5 [45], which is consistent with the bound z s ≤ 5.211 obtained by Ivanov [46]. In Table-I, II and III, we tabulated the calculated values of the center and surface redshift for the compact objects SAX J1808.4-3658, 4U 1538-52 and Her X-1 by taking the same values which we have used for graphical representation in Fig. 1.

IV. PHYSICAL FEATURES AND STABILITY ANALYSIS OF COMPACT OBJECTS
To gain some insight into the electrically charged fluids stellar model, we perform some analytical calculations and studied physical properties of the interior of the fluid sphere. The structure of charged spheres are analyzed by plotting several figures and studied equilibrium conditions under different forces. The obtained solution in this paper is used to study relativistic compact stellar objects within specified observational constraint.

A. Tolman-Oppenheimer-Volkoff (TOV) equations
To address the question, how gravitational and other different of fluid forces counteract with increasing electrostatic repulsion when the pressure gradients tend to vanish towards the boundary and charged fluid becomes more diluted.This situation may be illustrated by considering the hydrostatic equilibrium under different forces. This can be easily achieved by adapting generalized Tolman-Oppenheimer-Volkoff (TOV) equation [47,48] in the presence of charge, as prescribed by Ponce de Leon [49] is where M G = M G (r) is the effective gravitational mass inside a sphere of radius r and q = q(r) is given by (14). The expression for the effective gravitational mass is given by Now, plugging the value of M G (r) in Eq. (27), we get We have plotted different forces, namely, gravitational (Fg), hydrostatic (F h ) and electric forces (Fe), respectively, to describe the equilibrium condition for charged fluid sphere. The result shows that gravitational force is dominated by hydrostatic and electric forces to maintain the equilibrium condition. We follow the same procedure for finding solutions as given in Fig. 1.
Eq. (29), gives the information about the stellar equilibrium configuration for charged relativistic fluid, subject to the gravitational (F g ), hydrostatic (F h ) and electric forces (F e ), respectively, which are defined as: F e (r) = σ q r 2 e λ 2 = 1 8 π r 4 In Fig. 2, the behavior of these forces for the onset of hydrostatic equilibrium are shown for the compact star candidates SAX J1808.4-3658, 4U 1538-52 and Her X-1. It is worth mentioning that, the whole system is counterbalanced by the components of gravitational force (F g ), hydrostatic force(F h ) and electric force (F e ) and the system attains a static equilibrium.

B. Stability Analysis
Now, we are interested in analysing, the speed of sound propagation v 2 s , which is given by the expression v 2 s = dp/dρ. In natural the velocity of sound is less than the velocity of light. Here one can consider the speed of light is c = 1, so the sound speed is always less than unity. At this stage we investigate the sound speed for charged fluid matter and for stable equilibrium configurations this should belongs to the interval 0 < v 2 = dp dρ < 1, as in ref. [50] for a subluminal sound speed. We take Eq. (20) and (21) for obtaining an explicit solution, which is dp In spite of these expression complexity, we use the graphical representation to represent it more conveniently. Using the expressions for all the terms in this formula, we have plotted Fig. 3. From Fig. 3, it is clear that the velocity of sound lies within the proposed interval and therefore our model maintains stability. Our investigation show that our proposed model for charged perfect fluid star satisfies both energy and stability conditions.

C. Energy conditions
To clarify the question of whether the present model satisfy all the energy conditions within the framework of general relativity or not. Clearly, such configurations depend on the relationship between matter density and pressure obeying certain restrictions. In connection with that, energy conditions are essential tools to understand many theorems of classical general relativity -such as the singularity theorems of stellar collapse. Basic definitions are given for instance in [51,52]. Among them, we will be particularly interested in the (i) the Null energy condition (NEC), (ii) Weak energy condition (WEC) and (iii) Strong energy condition (SEC). Such conditions have the following inequalities ρ + 3 p + q 2 4πr 4 ≥ 0.
(34d) Inequalities (34) hold automatically for the sources considered here. The weak energy condition imposes the requirement of a positive energy density as measured by a distant observer. Using this inequalities one can easily justify the nature of energy conditions for the specific stellar configuration as shown in Fig. 4, that are satisfied for our proposed model.

D. Relativistic adiabatic index and Stability
We have so far considered stars in both hydrostatic equilibrium and speed of sound. But an important question that remains to be answered is whether these centenarians are enough for stable analysis. As another application is the adiabatic index. Using this one can incorporate all the basic characteristics of an equation of state on the instability formulae. In the case of an EOS of neutron star matter, the adiabatic index γ, varies from 2 to 4 and for an anisotropic compact star will be stable if γ > 4/3 [53]. For an adiabatic perturbation, the adiabatic index γ, can be expressed in the form [54].
Starting with Eqs. (20), (21) and (35), one can easily justify the adiabatic index and the stable configuration for the specific stellar configuration SAX J1808.4-3658, 4U 1538-52 and Her X-1. In order to clarify further we plot in Fig. 5, the dependence of the averaged adiabatic index, which clearly indicates that the configurations developed in this paper is stable. Finally, it should be noted that the adiabatic index which we define in Eq. (35), is a local characteristic of a specific EOS and depends on interior fluid density. For further reference we refer our reader to follow [55].

V. CONCLUSION
The main point we wish to make is that a general solution of the Vaidya-Tikekar model for a spherically symmetric superdense star has been carried out, and it is found a physically valid solution. The spheres contain a charged perfect fluid matter and the exterior spacetime is represented by the Reissner-Nordström metric. We also intend to study more closely the physical features of compact stars, generating an exact solution, and compated with number of specific astronomical objects. We summarize the obtained results for the interior of the stellar configuration as follows: (i) Regularity conditions: The energy density and pressure is positive inside the star. The central density ρ 0 = (3C(K − 1)) /K > 0 and central pressure p r (r = 0) > 0, which shows that the density and pressure are positive inside the star. (ii) Generalized TOV equation: We consider the generalized TOV equation for desribing the equilibrium condition subject to gravitational (F g ), hydrostatic (F h ) and electric forces (F e ), respectively and we observe from Fig. 2, that gravitational force is balanced by the joint action of hydrostatic and electric forces to attain the required stability of the model. However, the effect of electric force is less than the hydrostatic force. (iii) Causality conditions: In addition, we have shown that inside the charged fluid sphere the speed of sound is less than the speed of light i.e. 0 < (dp/c 2 dρ) < 1 for different values of 0 < K < 1 (see Table I-III for more detail results). However, we can see from Fig. 3 that the sound speed is increasing monotonically towards the surface of the star. (iv) Energy conditions: The configurations of the objects satisfy the null energy condition (NEC), weak energy condition (WEC) and strong energy condition (SEC), simultaneously inside the star as displayed in Fig. 4.
(v) Surface Redshift: We have also studied another impotent features is the surface redshift of stellar structure. A general feature of the behavior is that it attend maximum value at the center and monotonically decreasing towards the boundary, which we can see from figure Fig. 5 (left panel). The numerical values corresponding to center and surface redshift for the SAX J1808.4-3658, 4U 1538-52 and Her X-1 are enlisted in Table I-III. (vi) Stability conditions: We have also considered the adiabatic index to gain some useful information at an extreme conditions and for stable configuration of a compact object. In this case it is possible to approximation the adiabatic index γ > 4/3. We demonstrate that in Tables I-III and graphical form as displayed in Fig. 5 (right panel). (vii) Mass-radius ratio: As an another application of the obtained upper bound of mass-radius relation together with the Buchdahl limit for perfect fluids satisfying the inequality R ≤ (9/8) R S = (9/4) G M/c 2 [41]. For a compact configuration with electric charge to the system, the Buchdahl limit is certainly increase. It has been found that Buchdahl limits for charged stars in [56][57][58]. In our case the obtained mass-to-radius ratio fall within the bounds as proposed for charged spheres (see Table -4 for more).
In conclusion, we will investigate other forms of metric potentials that could exhibit more general behavior and to make observationally distinguish between these compact objects.