Exact solution of anisotropic compact stars via. mass function
Abstract
The interest in studying relativistic compact objects play an important role in modern astrophysics with an aim to understand several astrophysical issues. It is therefore natural to ask for internal structure and physical properties of specific classes of compact star for astrophysical observable, and we obtain a class of new relativistic solutions with anisotropic distribution of matter for compact stars. More specifically, stellar models, described by the anisotropic fluid, establish a relation between metric potentials and generating a specific form of mass function are explicitly constructed within the framework of General Relativity. New solutions that can be used to model compact objects which adequately describe compact strange star candidates like SMC X1, Her X1 and 4U 153852, with observational data taken from Gangopadhyay et al. (1). As a possible astrophysical application the obtained solution could explain the physics of selfgravitating objects, might be useful for strongfield regimes where data are currently inadequate.
I. Introduction
Over the past few years, there has been growing interest in static spherically symmetric compact object consistent with general relativity which are most often observed as pulsars, spinning stars with strong magnetic fields. However, in spite of the fact, it is important to know the exact composition and nature of particle interactions, which allow us to completely describe them in terms of their mass, spin angular momentum and charge. In particular, observations of compact stars are considered primary targets of the forthcoming field of gravitational wave astronomy. In the recent times we have experimental evidence that such objects do exist from observations with very high densities (2), and some of the compact objects like Xray pulsar Her X1, Xray burster 4U 182030, Xray sources 4U 172834, PSR 0943+10, millisecond pulsar SAX J 1808.43658, and RX J1856353754 strongly favour the possibility that they could actually be strange stars. In other words, there is no any strong evidence to make conclusion/ understood the mechanism about compact objects. In spite of this drawback, the internal composition and consequent geometry of such objects are still considered an open question in the scientific community. From a theoretical point of view, such compact objects are composed of a perfect fluid (34). Generally, polytropic equation of state (EOS) in the form P = k and bag model (4); (5) have widely used to describe a white dwarf or a less compact star (6). An example, Herrera and Barreto (7) carry out a general study on polytropic general relativistic stars with anisotropic pressure whereas Lai and Xu (8) have studied polytropic quark star model.
Though, interior of a star is an important astrophysical question and hence it is pertinent to construct relativistic models by assuming anisotropy fluid distribution. The theoretical study of the influence of anisotropic compact objects was first initiated by Bowers and Liang (9) and another investigation led be Ruderman (10) showed that nuclear matter may have anisotropic features at least in certain very high density ranges g/, where the nuclear interaction must be treated relativistically. Later on, several models (see for example, (11); (12); (13)) have been proposed in this direction. The procedure has been found to be interesting and useful on the physical grounds that anisotropy affects the critical mass, massradius relation and stability of highly compact relativistic bodies. It is also wellestablished fact that a magnetic field acting on a Fermi gas produces pressure anisotropy was discussed recently in (14); (15); (16). Recently, Maurya and Gupta in (17); (18) have studied charged anisotropic stars whereas without charged solution has been analysed in (19); (20). The role of anisotropy, with the linear equation of state was pursed by Mak and Harko (21), and extending the work by Sharma and Maharaj (22) obtained an exact analytical solutions assuming a particular form of mass function. In yet another paper, Victor Varela et al (23) studied Charged anisotropic matter with linear or nonlinear equation of state .
However, it has recently been proposed that an arbitrary d dimensional (pseudo) Riemannian space can always be locally embedded into any Riemann space of dimension N d(d + 1)/2. Riemann’s seminal work in 1868 (24) had inspired Schlfli that how one can locally embed such manifolds in Euclidean space (25). He discussed the local form of the embedding problem and conjectured that maximum number of extra dimensions that can embedded as a submanifold of a Euclidean space with N= or, equivalently, when the codimension r = N  n it gives r = n(n  1)/2. Furthermore, one motivation in the treatment of embedding has been devoted, namely, Randall Sundrum braneworld model is based on the assumption that four dimensional spacetime is a threebrane, embedded in a fivedimensional Einstein space (26). On the other hand Nash in 1956 (27) established the idea of global isometric embedding theorem of into Euclidean space .
In a sense one could say that all Riemannian manifold has a local and a global isometric embedding in an Euclidean space. This result opened new perspectives in embedding theorems with increasing degrees of generality and soon it became an powerful tool to construct and classify solutions of GR. Inspired by these advances, a popular approach has been emerged in embedding 4dimensional spacetime into 5dimensional flat spacetime by using the spherical coordinates transformation and known as embedding class one if it satisfies the Karmarkar condition (28). One of the primary motivations of such embedding is to establish a relationship between metric potentials and obtained exact solutions of Einstein’s equation to a singlegenerating function. In connection to this an exact anisotropic solution of embedding class one, has been developed in (29); (30); (31). They showed that in seeking solution for relativistic static fluid spheres one can utilized this technique successfully. Motivated by the above facts, we will see how this theorem can be utilized for obtaining an exact solution of Einstein’s field equation for compact star models. In this article we consider the static spherically symmetric spacetime metric with embedding class one conditions, which can be altered to fit with a set of astrophysical objects.
This paper is outlined in the following manner. In Sec. II, we present the structural equations for anisotropic fluid distributions of stellar models applying the embedding class one conditions. Specific models are then analyzed in a brief description by obtaining a particular form of mass function. In Section III, we match our interior solution to an exterior Schwarzschild vacuum solution at the boundary surface and then determine the values of constant parameters. In Section IV, we continue our discussion thorough geometrical analysis of the solution such as energy conditions, hydrostatic equilibrium and stability of the star by fixing of certain parameters. We end our discussion by concluding remarks in V.
II. Class one condition for Spherical symmetric metric and General relativistic equations
The simplest configuration for a star is the static and spherically symmetric geometry has the usual form
(1) 
where the coordinates (t, r, , ) are the spherical coordinates and the metric coefficient and are the functions of the radial coordinate , and yet to be determined by solving the Einstein equations.
As the metric (1) is time independent and spherically symmetric, we restrict ourselves that the space time is of emending class one, if it satisfies the Karmarker condition (see Ref. (28); (29) for more details discussion) and metric functions and are dependent on each other as:
(2) 
where , and . However = constant, leads to flat space time metric.
The energy momentum tensor associated to a spherical distribution of matter bounded by gravitation is locally anisotropic, that is = , where and are the radial and tangential pressures and is the energy density of the fluid, respectively. Thus, the Einstein field equation, = 8, where is the Einstein tensor then reduce to the following ordinary differential equations for the metric (since we use natural units where G = c = 1) as:
(3) 
where the primes () denote differentiation with respect to r. Then one may write the solution in a very familiar form of mass function of the compact object,
(4) 
It is worth noting that the mass is the density inside a proper volume element within a radius r. Using the above Eq. (4) and Eq. (2), we can write the metric component which is given by the equality
(5) 
Rewriting the above expression in terms of mass function as:
(6) 
Since and because , and which is already proved by Herrera et al. (38) and Maurya et al. (39). From Eq. (6), we obtain and . This solution clearly implies that is positive and monotonic increasing function of r.
It is possible to derive the mass function in terms of F(r). Aiming this the value of F(r) from Eq. (6), we have
(7) 
To solve this integral, we shall assume , that we claim should be zero at centre i.e. due to m, which indicating that the mass function is regular at the origin. Equation (7) can be rewritten as
(8) 
To make the above integral equations tractable and construct a physically viable model, we suppose that in a particular form, which gives
(9) 
where two constants a, b are positive and makes the Eq. (8) in a simplest form
(10) 
Solving the integral gives
(11) 
where A stands for integrating constant. Substituting Eq. (11) in Eq. (2), we obtain one of the metric coefficient in the form
(12) 
and the mass function is given by
(13) 
The motivation of choice for particular dimensionless function lies on the fact that the obtained mass function a monotonic decreasing energy density in the interior of the star. To construct a physically viable model this type of mass function is not new but similar works have been considered earlier by Matese & Whitman (35) and Finch & Skea (36) for isotropic fluid spheres, and Mak & Harko (37) for anisotropic fluid spheres.
The system of equations used to study for spherically symmetric configurations with anisotropic fluid distribution are
(14) 
(15) 
(16) 
where we use for our notational convention.
Firstly, we will present the anisotropic effect by a term , by taking into account Eqs. (1516), may be expressed in the following equivalent form
(17) 
which representing a force that is due to the anisotropic nature of the fluid. The anisotropy will be repulsive or directed outwards if , and attractive or directed inward when . The effects of anisotropy forces maintain the stability and equilibrium configurations of a stellar stricture, as we discuss later.
In more detail, it is important to emphasize the results. That’s why we extend our calculation for the first order differential equation, which are
(18) 
(19) 
(20) 
for our notational conventional we use
, , and
,
,
,
,
.
We exercised our results in various ways first by checking calculation and putting restriction on the physical parameter based on logarithmic principle and then by using graphical representation which are illustrated in Figs. (13), describing the metric functions, energy density, radial and transverse pressures and measure of anisotropic within the given radius. In the next section we will look for supplementary restrictions on the model to make it physically viable.
III. Boundary conditions
To do the matching properly, we start by joining an interior spacetime , to an exterior , Schwarzschild vacuum solution at the boundary surface r = R. At this boundary the metric should be continuous. The Schwarzschild solution is given by
where denote the total mass of the compact star. In deciding criterions for an anisotropic compact star the radial pressure must be finite and positive inside the stars, and it should vanishs at the boundary of the star (40).
In determining the value of constrains parameters we use the boundary condition , which can express as
(22) 
In order to match smoothly on the boundary surface, we must require the continuity of the first and the second fundamental forms across that surface. Then it follows the condition which gives the value of constant parameter B as
(23) 
and using the condition , we have
(24) 
This represents the total mass of the sphere as seen by an outside observer. Additionally, bounds on stellar structures is an important source of information and classification criterion for compact objects, to determine the massradius ratio which was proposed by Buchdahl (33). This bound has been considered for thermodynamically stable perfect fluid compact star with ratio 2M/R, must be less than 8/9. We have carried out the analysis for our compact star candidates SMC X1, Her X1 and 4U 153852, that are used to calculate the values of constants A and B. For this purpose, we present these results (central density, surface density, central pressure, massradius ratio) in table I & II, respectively.
Table I  Numerical values of the parameters or constants for different compact stars (1)  

Compact Star  R (km)  M  A  
SMC X1  8.301  1.04  0.0000270  0.006648  0.439167064  452.08447 
Her X1  8.1  0.85  0.0000185  0.005620  0.481332147  618.82624 
4U 153852  7.866  0.87  0.0000247  0.006300  0.481382726  490.02078 
Let us now focus on surface gravitational redshift, which gives a wealth of information about compact objects, and defined by Z = = , where is the emitted wavelength at the surface of a nonrotating star and the observed wavelength . Thus one can estimate the gravitational redshift, defined by , from the surface of the star as measured by a distant observer by the following relation
(25) 
where = = . Measurement of the gravitational redshift for a static perfect fluid sphere is not larger than = 2 (33), whereas for an anisotropic fluid sphere this value may be increase up to = 3.84, as given in Ref. (34). We are trying to estimate the surface redshift given table II for different compact star candidates.
Table II  Estimated physical values based on the observational data and theory  

Compact star  
SMC X1  1.0710  7.3060  1.0052  0.36959 8/9  0.25948 
Her X1 
9.0538  6.6539  6.3783  0.30957 8/9  0.20348 
4U 153852  1.0149  7.3928  7.9036  0.32628 8/9  0.218346 
IV. Physical features and Comparative study of the physical parameters for compact star model
To proceed further discussion based on the obtained solution that must satisfy some general physical requirements. In order to simplify the analysis and make the solution more viable we explore some physical features of the compact star and carry out a comparative study between the data of the model parameters with a set of astrophysical objects in connection to direct comparison of some strange/compact star candidates.
.1 Energy Conditions
Let us first discuss a very simple but important features of a stellar model. The study of energy conditions within the framework of GR is an essential part for studying the compact objects. Here we examine the energy conditions, namely : (i) Null energy condition (NEC), (ii) Weak energy condition (WEC) and (iii) Strong energy condition (SEC), at all points in an interior of a star holds simultaneously, by the following inequalities
(26)  
(27)  
(28) 
From the above inequalities we provide a graphical representation that one can easily justify the nature of energy condition for three different compact objects in Fig. 3. Due to the complexity of the expression, we only able to write down the above inequalities and plotted the graphs against the above criterions. In fact, all the energy conditions are not violated and well behaved in the stellar interior.
.2 Generalized TolmanOppenheimerVolkov Equation
On the other hand, for a given compact star, it is possible to test for hydrostatic equilibrium under the different forces namely gravitational, hydrostatic and anisotropic forces. Nevertheless, in our investigations we need to apply generalized TolmanOppenheimerVolkov (TOV) equation (42); (43), which is commonly used for an anisotropic fluid distribution is given by
(29) 
where the effective gravitational mass is defined by
(30) 
Then equation (29) may be rewritten as
(31) 
Let us now attempt to explain the Eq. (31) from an equilibrium point of view, which was first shown by Tolman (41) and Oppenheimer and then Volkoff (42), where they predicted the stable equilibrium condition for the compact star as a sum of three different forces, viz. gravitational force (), hydrostatics force () and anisotropic force (). Thus the above condition assumed the following form, namely
(32) 
The components of forces can be expressed in explicit form as:
(33) 
(34) 
(35) 
At this point in the derivation, a stable configuration of the system is counter balance by the components of different forces. To illustrate this in more detail let us introduce an graphical representation as evidenced in Fig. 4.
Here, the combine forces of hydrostatic () and anisotropic (), dominate the gravitational force (), while the anisotropic stress has a less role to the action of equilibrium condition. Hence, following the reason outlined above and making use of Eq. (32), we established the stable configuration model as shown in Fig. 4.
.3 Stability Analysis
In this section we analyze the speed of sound propagation , is given by the relation . For a physically interesting stellar geometry should require that the sound speed does not exceed the speed of light, i.e., the physically relevant region is always less than unity. Here we will investigate the speed of sound for anisotropic fluid distribution and propagating along radial as well as transverse direction, should satisfy the bounds and , (32).
We first carry out an analysis of sound velocity with graphical representation. Our result for obtained velocity of sound for compact star with anisotropic matter is presented in Fig. 5, for strange star candidates SMC X1, Her X1 and 4U 153852. As the resulting expressions are very elaborated, so we only plot the results for different compact stars. It is interesting to note that both and monotonic decreasing function, which all show a behavior similar to that found for other compact objects. This is a sufficient condition for the solution to be causal. A notable characteristic of the sound velocity is the estimation of the potentially stable and unstable eras, by considering the expression , for stable potential (46). According to Fig. 5 (extreme right) our solution gives us the stable star configuration.
Our investigation shows that obtained mass function by using the Karmarker condition for strange compact star matter, satisfies both energy and stability conditions
I Concluding Remarks
To summarize, we have studied compact objects supported by anisotropic fluid distribution, which plays an important role in preventing gravitationally collapse. The motivation behind such a construction is that, we established a relation between metric potentials by imposing an embedding theorem, known as embedding class one within the framework of GR. As a result of this approach, we generate a mass function which we have used to study the interior of stellar objects with density decreasing outwards. In next, we have started by deriving the basic equations of Einsteinâs field equations that describing the structure of compact objects.
In section (II), using the structural equations we derived the energy density, radial and transverse pressures and measure the nature of anisotropy. The complicated expressions given by equations (1417) are plotted as a function of the radius for our compact stars candidates SMC X1, Her X1 and 4U 153852. An important feature of this model is the energy density might not vanish at the boundary r =R, though the radial pressure vanishes for all parameter values in according to the boundary condition. Additionally, at the boundary the interior spacetime have been matched by a Schwarzschild metric, and determine the values of arbitrary constants A and B for our compact stars candidates SMC X1, Her X1 and 4U 153852 in section (III). A comparative study of our results with that of the compact star candidates are provided in Table I and II.
Once the mass function is specified, in order to close the system based on physical requirements, we further proceed by investigating the energy conditions, hydrostatic equilibrium under the different forces, and velocity of sound. All the physical properties are well behaved within the stellar radius. We also showed that upper bound of the massradius ratio must be less than 8/9 as proposed by Buchdahl (33), for different compact star candidates which we have used for our model.
This indicates that the approach adapted in this paper is likely to produce other meaningful models with specific mass function that also greatly help in understanding the properties of other different static compact configurations. In our further study it would be interesting to investigate other forms of mass function that exhibit more general behaviour and thereby describe strange stars related with observational details.
Acknowledgments
AB is thankful to the authority of InterUniversity Centre for Astronomy and Astrophysics, Pune, India for providing research facilities.
References
 T. Gangopadhyay et al.: Mon. Not. R. Astron. Soc., 431, 3216 (2013).
 J. Lattimer: (2010) http://stellarcollapse.org/nsmasses.
 B. V. Ivanov: Phys. Rev. D, 65, 104011 (2002).
 E. Witten: Phys. Rev. D, 30, 272 (1984).
 N. K. Glendenning, Ch. Kettner & F. Weber: Phys. Rev. Lett., 74, 3519 (1995).
 S. L. Shapiro & S. A. Teukolosky: Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects (Wiley, New York,1983).
 L. Herrera & W. Barreto: Phys. Rev. D., 88, 084022 (2013).
 X.Y. Lai & R.X. Xu: Astropart.Phys., 31, 128134 (2009).
 R.L. Bowers & E. P. T. Liang: Class. Astrophys. J., 188, 657 (1974).
 R. Ruderman: Rev. Astr. Astrophys., 10, 427 (1972).
 K. Dev & M. Gleiser: Gen. Rel. Grav., 34, 1793, (2002).
 K. Dev & M. Gleiser: Gen. Rel. Grav., 35, 1435, (2003).
 L. Herrera, J. Martin & J. Ospino: J. Math. Phys., 43, 4889, (2002).
 M. Chaichian et al.: Phys. Rev. Lett, 84, 5261 (2000).
 A. Perez Martinez, H. Perez Rojas & H. J. Mosquera Cuesta: Eur. Phys. J. C, 29, 111 (2003).
 E. J. Ferrer et al.: Phys. Rev. C, 82, 065802 (2010).
 S.K. Maurya et al.: Eur. Phys. J. C, 75, 389 (2015).
 S.K. Maurya et al.: Astrophys. Space Sci., 361, 163 (2016).
 Sk. Monowar Hossein et al.: Int.J.Mod.Phys. D, 21, 1250088 (2012).
 Mehedi Kalam et al.: Eur. Phys. J. C, 72, 2248 (2012).
 M.K. Mak & T. Harko: Chin. J. Astron. Astrophys.,, 2, 248, (2002).
 R. Sharma & S. D. Maharaj: Mon. Not. R. Astron. Soc., 375, 1265 (2007).
 Victor Varela et al.: Phys.Rev. D, 82, 044052 (2010).
 B. Riemann & Abh. Knigl: gesellsch., 13, 1 (1868).
 L. Schl: Ann. di Mat. srie 5, 170 (1871).
 L. Randall & R. Sundrum: Phys. Rev. Lett., 83, 3370 (1999); : Phys. Rev. Lett., 83, 4690 (1999).
 J. Nash: Ann. Math., 63, 20 (1956).
 K. R. Karmarkar: Proc. Ind. Acad. Sci. A, 27, 56 (1948).
 S. K. Maurya et al.: Eur. Phys. J. A, 52, 191 (2016).
 S. K. Maurya et al.: Eur. Phys. J. C, 77, 45 (2017).
 Ksh. Newton Singh et al.: Int.J.Mod.Phys. D, 25, 1650099 (2016).
 L. Herrera : Phys. Lett. A, 165, 206 (1992).

H. A. Buchdahl: Phys. Rev., 116, 1027 (1959);
H. A. Buchdahl: Astrophys. J., 146 275, (1966).  B. V. Ivanov: Phys. Rev. D, 65, 104011 (2002).
 J. J. Matese & P. G. Whitman: Phy. Rev. D, 11, 1270 (1980).
 M. R. Finch & J. E. F. Skea: Class. Quantum Grav., 6, 467 (1989).
 M. K. Mak & T. Harko: Proc. Roy. Soc. Lond., A, 459, 393 (2003).
 L. Herrera, J. Ospino & A. Di Parisco: Phys. Rev. D, 77, 027502 (2008).
 S. K. Maurya, Y.K. Gupta & S. Ray: arXiv: 1502.01915 [grqc] (2015).
 C. W. Misner & D. H. Sharp: Phys. Rev. B, 136, 571 (1964).
 R.C. Tolman: Phys. Rev., 55, 364 (1939).
 J.R. Oppenheimer & G.M. Volkoff: Phys. Rev., 55, 374 (1939).
 J. Ponce de León: Gen. Relativ. Gravit., 25, 1123 (1993).
 V. Varela: Phys. Rev. D, 82, 044052 (2010).
 J. Devitt & P.S. Florides: Gen. Relativ. Gravit., 21, 585 (1989).
 H. Andreasson: Commun. Math. Phys., 288, 715 (2009).