Anisotropic fluid spheres of embedding class one using Karmarkar condition

We obtain a new anisotropic solution for spherically symmetric spacetimes by analyzing the Karmarkar embedding condition. For this purpose we construct a suitable form of one of the gravitational potentials to obtain a closed form solution. This form of the remaining gravitational potential allows us to solve the embedding equation and integrate the field equations. The resulting new anisotropic solution is well behaved, which can be utilized to construct realistic static fluid spheres. Also we estimated the masses and radii of fluid spheres for LMC X-4, EXO 1785-248, PSR J1903+327 and 4U 1820-30 by using observational data set values. The masses and radii obtained show that our anisotropic solution can represent fluid spheres to a very good degree of accuracy. The physical validity of the solution depends on the parameter values of a, b and c. The solution is well behaved for the wide range of parameters values 0.00393≤a≤0.0055\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.00393\le a \le 0.0055$$\end{document}, 0.0002≤b≤0.0025\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.0002 \le b \le 0.0025$$\end{document} and 0.0107≤c≤0.0155\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.0107 \le c \le 0.0155$$\end{document}. The range of corresponding physical parameters for the different compact stars are 0.3266≤vr0≤0.3708\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.3266\le v_{r0} \le 0.3708$$\end{document}, 0.1583≤vt0≤0.2558\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.1583\le v_{t0} \le 0.2558$$\end{document}, 0.3256≤zs≤0.4450\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.3256\le z_{s} \le 0.4450$$\end{document} and 4.3587≤Γ0≤5.6462\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.3587\le \Gamma _{0} \le 5.6462$$\end{document}.


Introduction
The presence of nonzero anisotropy, in which the radial pressure differs from the tangential pressure, is an important component in modelling relativistic stellar systems in the absence of an electric field. The pioneering paper by Bowers and Liang [1] introduced anisotropic spheres in general relativity. Subsequently there has been extensive research in studying the physics related to anisotropic pressures. It was shown by Dev and Gleiser [2,3] that pressure anisotropy influences a e-mail: sunil@unizwa.edu.om b e-mail: maharaj@ukzn.ac.za the mass, structure and physical properties of highly compact spheres. It is important to observe that the mass of the object and the redshift both change with different values of the magnitude of the anisotropy. In recent times there has been made a considerable effort in modelling observed astrophysical objects in the presence of anisotropy. Some recent research papers addressing this physical feature include the work of Sharma and Ratanpal [4], Ngubelanga et al. [5,6], Sunzu et al. [7,8], Murad and Fatema [9,10] and Murad [11], and the references therein. The physical analyses contained in these treatments confirm the importance of including nonzero anisotropy in modelling astrophysical objects. It should be noted that the presence of anisotropic pressures leads to values of observed compactness parameters for several astrophysical bodies including Her X-1, PSR 0943+10, 4U 1820-30, SAX J 1808.4-3658, and 4U 1728-34. It is therefore important to generate analytical models of the Einstein field equations, with a matter tensor containing anisotropy, which are consistent with physical requirements for astrophysical applications.
To generate a physically consistent model we need to find an analytical solution of the general relativistic field equations. This is achieved by restricting the spacetime geometry, the matter content or specifying an equation of state. A rather different approach is to use the embedding of a fourdimensional differentiable manifold into a higher dimensional Euclidean space. Embedding of curved spacetimes into spacetimes of higher flat dimensions has assisted in generating several new exact models in cosmology and relativistic astrophysics [12]. The embedding condition gives an additional differential equation, the so-called Karmarkar condition [13], in static spherical spacetimes relating the two gravitational potentials. A solution of the Karmarkar condition then helps to solve the Einstein field equations. This method has proved to be a fruitful mechanism to find new relativistic astrophysical models in recent investigations. Bhar et al. [14][15][16], Maurya et al. [17][18][19][20] and Newton Singh et al. [21][22][23] have generated different families of physically acceptable Karmarkar spacetimes that describe the interior regions of acceptable stars. In these analyses a particular form of one of the potentials is made, which solves the Karmarkar condition, giving the second potential, eventually leading to an astrophysical model. In our treatment we show that a very general form of the chosen potential, including hyperbolic functions, leads to an astrophysical model with desirable physical features. This shows that the Karmarkar condition allows for more complicated (and acceptable) physical behaviour than the earlier simpler forms used for the gravitational potentials.
In this paper we present a new interior anisotropic model for astrophysical compact stars by solving the embedding condition in static spherical spacetimes. We show that the resulting exact solutions can be used to describe a physically reasonable astrophysical matter distribution. The exterior region is characterized by the Schwarzschild vacuum metric. We discuss the Einstein field equations in Sect. 2 and present the Karmakar embedding condition. In Sect. 3, we generate an exact solution to the embedding condition and show how this leads to an exact solution of the field equations. In Sect. 4 we present the matching conditions between interior and exterior spacetimes regions; we also demonstrate that the parameters arising are bounded. A detailed physical analysis is undertaken in Sect. 5. In particular the stability, cracking and energy conditions of the relativistic sphere are studied. We also investigate the physical features of the model in connection with the objects LMC X-4, EXO 1785-248, PSR J1903+327 and 4U 1820-30; the results are presented in the form of tables and graphs. We discuss the significance of the results obtained in this paper in Sect. 6.

Field equations and the Karmarkar condition
We assume that the interior matter of the star is locally anisotropic which is given by the following line element (by taking c = 1) where e λ(r ) and e ν(r ) represent the gravitational potentials of stellar structure. The Einstein field equations corresponding to an anisotropic fluid distribution is given by (assume where T i j and R i j represent the energy-momentum tensor and contracted Ricci tensor, respectively, while R is the scalar curvature tensor. The energy tensor for the anisotropic matter distribution can be defined as where the contravariant quantity v i is the four-velocity vector and u i is the unit spacelike vector in the radial direction. Here p r , p t and ρ denote the radial pressure, tangential pressure and matter density for anisotropic matter. In view of line element (1), the Einstein field equations (2) provide the following differential equations for the anisotropic star: 2r Here primes denote the derivative with respect to the radial coordinate r . The value of the velocity of light (c) and the gravitational constant (G) are taken to be unity in the above coupled differential equations. Furthermore we obtain the anisotropic factor by using the pressure isotropy condition with Eqs. (5) and (6) as

Karmarkar condition
It is well known that the spherical symmetric line element (1) can always be embedded in six dimensional flat spacetime which implies that the spherical symmetric line element is of embedding class two in general. On the other hand we can also embed the spherical line element into five dimensional flat spacetime if it satisfies the Karmarkar condition [13]. Then it represents the spacetime of embedding class one. However, it is a necessary and sufficient condition for the spherically symmetric spacetime to be of class one. The Karmarkar condition is given in terms of the curvature components by By plugging these components of hi jk in Eq. (8) we get the following differential equation: On solving the differential equation (9) we obtain the potential, where C and D are nonzero arbitrary constants of integration.

Tolman-Oppenheimer-Volkoff equation
Now to derive the Tolman-Oppenheimer-Volkoff (TOV) equation we evaluate p r + ρ from Eqs. (4) and (5) to get The derivative of the radial pressure is Then using Eqs. (7), (11) and (12) we get If the gravitational mass within a compact star of radius r is denoted M G (r ), then it is given by the Tolman-Whittaker formula [25]: Then from Eqs. (13) and (14), we obtain The above equation (15) represents the well known generalized Tolman-Oppenheimer-Volkoff (TOV) equation which provides the equilibrium condition for anisotropic stellar matter distribution.

Anisotropic solution of embedding class one
Equations (4)-(5) have five unknowns, namely ν, λ, ρ, p r and p t . However, the Karmarkar condition provides a relation between ν and λ, which implies that we have four con-ditions (including three equations) to solve this system of equations. For this paper we consider a totally new expression for gravitational potential e λ , which has not been used before. We take where a, b and c are nonzero constants, The units of the constants a and c are length −2 . We need to check whether the given expression for λ is physically valid. For this purpose we obtain e λ at the centre and plot Fig. 1. We observe that it is increasing monotonically away from the centre and e λ = 1 at the centre. This behaviour of e λ indicates that it is physically acceptable. By plugging the value of λ into Eq. (10), we obtain 2 is finite and positive at the centre. Also we may observe from Fig. 1 that the gravitational potential e ν is increasing with r throughout the star. This implies that the above expression of ν may be suitable to obtain a physically valid anisotropic solution according to Lake [26].
The expressions for the matter density, radial and tangential pressures are obtained (by taking ψ = a r 2 + b, We plot the variation of matter density, radial and tangential pressures in Fig. 2. We can see the density is maximum at the centre and the minimum occurs on the boundary of the star. The radial and tangential pressures are both monotonically decreasing away from the centre. However, p r becomes zero at the boundary of the star, which gives the radius of the star. Note that p t is nonzero and positive. It is worth pointing out here that the central density is of the order of 10 15 g/cm 3 , which indicates that the nuclear matter is more appropriate for the anisotropic fluid distribution. The expression for the anisotropic factor is given by The pressure anisotropy of the matter plays an important role in determining the stability of the model. Since the force due to anisotropy can be defined as F = 2 r , which implies if is positive then the direction of the force is outward; however, the force will be directed inward if is negative. But the existence of a repulsive force (in the case in which tangential pressure is higher than radial pressure i.e. p t > p r ) allows the construction of more compact star models when using an anisotropic fluid than when using an isotropic fluid [27]. From Fig. 3, we can see that is positive and finite throughout inside the star. Also it is zero at the centre and attains a maximum at the boundary of the star.
The gradients of pressure and density are given by In the above we have set +a c r 2 (r ) cosh 2ψ + 2 a 2 r 2 cosh 3ψ ,

Bounds on the parameters
Since the radial pressure p r and the tangential pressure p t are positive and finite inside the star we obtain the upper bound of B as Also the fluid model must satisfy the Zeldovich condition i.e. p r /ρ < 1 and p t /ρ < 1 everywhere inside the star, which gives the lower bound of B:    (26), we get the following inequality: The behaviour of p r /ρ and p t /ρ is shown in Fig. 4. From this figure it is clear that both p r /ρ and p t /ρ are less than 1 everywhere inside the anisotropic star, which shows that our fluid model satisfies the Zeldovich condition.

Matching conditions
To find the arbitrary constants A and B, we must match our interior solution to the exterior Schwarzschild solution at the boundary of the star. The line element of the exterior Schwarzschild solution is given by where the constant mass M provides the total mass of the anisotropic star within r = R. By matching of e ν and e λ at the surface of the star (r = R) (continuity of first fundamental form) we get For fixing the arbitrary constants, the second fundamental form (∂g 44 /∂r ) also has to be matched at the boundary r = R. By matching of ∂g 44 /∂r at the surface of the star we get a zero radial pressure at the boundary [28]. Then p r = 0 at r = R provides the value of the arbitrary constant, Then using Eqs. (29) and (30) together with Eq. (31), we obtain the values of constant A and total mass M:

Well-behaved property of the solution
For a well-behaved nature of the solution, the velocity of sound must be less than the velocity of light, and it should decrease monotonically throughout the anisotropic star (Fig. 5). For this purpose we have to calculate the radial and tangential speed of sound as

Dominant energy conditions
For a physically reasonable anisotropic solution the energy momentum tensor has to obey the following dominant energy conditions: (i) The null energy condition (NEC) implies that local mass-energy density must not be negative: ρ ≥ 0, (ii) The weak dominant energy condition (WDEC) implies that the flow of energy inside star must not be faster than the velocity of light: ρ − p r ≥ 0, ρ − p t ≥ 0, (iii) The strong dominant energy condition (SDEC) implies that the flow of energy inside the star must not be faster than one-third of the light velocity: ρ − 3 p r ≥ 0, ρ − 3 p t ≥ 0.
From Fig. 6 it is clear that all energy conditions are satisfied within the compact star.

Mass-radius relation
For any physical valid star model according to Buchdahl [29] the mass to radius ratio for a perfect fluid model should be 2M/R < 8/9. Later on Mak and Harko [31] have proposed this relation in a more generalized form which can be written In this connection we would like to compare our proposed compact star model with the observed data of different realistic objects. For this purpose we have calculated the physical parameters for the models (see Tables 1, 2) by taking the mass of the compact stars LMC X-4, EXO 1785-248, PSR J1903+327 and 4U 1820-30. The obtained radii of the different compact stars are given in Table 1, which are in good agreement with the proposed values of Gangopadhyay et al. [30].

Surface redshift
As is well known the compactification parameter of the compact star is given by Then in connection with the above compactification parameter the surface redshift (z s ) can be obtained directly by From Eq. (38) we can see that the surface redshift cannot be arbitrarily large as it depends upon the compactness parameter u = m/R. The behaviour of the redshift inside the star can be seen in Fig. 7, which can be obtained by the formula z = e −ν/2 − 1.

Stability of anisotropic models via cracking
In our anisotropic fluid model, to verify stability we plot the radial (v r = √ d p r /dρ) and transverse (v t = √ d p t /dρ) sound speeds in Fig. 5. It can be observed that both velocities satisfy the inequalities 0 < v 2 r < 1 and 0 < v 2 t < 1 everywhere within the stellar object (Fig. 8), which obeys the anisotropic fluid models [32,33].
To check whether the local anisotropic matter distribution is stable or not, we use the proposal of Herrera [32], known as the criterion of cracking (or overturning) of the star. This indicates that the region is potentially stable where the radial velocity of sound is greater than the transverse velocity of sound. We can easily say that 0 < v 2 r − v 2 t < 1 (dotted line) and −1 < v 2 t − v 2 r < 0 (solid line) as can be seen in Fig. 9. Hence, we conclude that our compact star model provides a stable configuration.

Stability via adiabatic index
The stability of the relativistic anisotropic star also depends upon the adiabatic index . Heintzmann and Hillebrandt [34] proposed that neutron star models with anisotropic equation of state are stable if γ > 4/3. However, according to Newton's theory of gravitation the isotropic neutron star model has no upper mass limit for the adiabatic index γ > 4/3 [35]. The adiabatic index is defined by For a relativistic anisotropic fluid sphere the stability condition is given by where p r 0 , p t0 , and ρ 0 denote the initial radial pressure, tangential pressure and energy density, respectively, in the static equilibrium condition which satisfies the TOV equation (15). However, the second and last terms inside the square brackets correspond to the anisotropic and relativistic corrections (being positive quantities), which increase the instability range of the adiabatic index. For this purpose we have plotted vs. r/R. The behaviour of the adiabatic index is shown in Fig. 10, which shows that > 4 3 everywhere inside the compact star model.

Stable equilibrium condition via TOV equation
The Tolman-Oppenheimer-Volkoff (TOV) equation describes the interior structure of the compact star which is a relationship between two physical quantities, the radial pressure and the density. Using the TOV equation we want to examine whether our present model is in a stable equilibrium stage under the three forces, viz. anisotropic force (F a = 2( p t − p r )/r ), hydrostatics force (F h = −d p r /dr ) and gravitational force (F g = −ν (ρ + p r )/r ). This implies that the sum of three different forces becomes zero: The explicit form of the above three different forces for this solution is given by The behavior of above forces can see in Fig. 11.

Harrison-Zeldovich-Novikov stability criterion:
Recently, Singh et al. [40] have employed the Harrison-Zeldovich-Novikov [41,42] criterion for investigating the stability of the anisotropic compact star models (Fig. 12). This criterion is that the configuration is stable only if the mass of the compact star is increasing with central density i.e. dM/dρ 0 > 0; however, it is unstable if dM/dρ 0 < 0. Let us define the mass function of our static anisotropic fluid solution in terms of the central density by The derivative of Eq. (45) with respect to ρ 0 gives From Fig. 12, it is clear that our anisotropic models are stable according to the Harrison-Zeldovich-Novikov stability criterion. It is noted that the stability of our configurations is increased with increasing radii and there is no change after attaining a maximum value for the respective central densities. It important to note that the curve starts leaving off, i.e. dM/dρ 0 = 0, after reaching its maximum value, which indicates that the configuration is rendered unstable. 5.6 Herrera's generating functions for the present embedding class one solution: The algorithm for all possible anisotropic fluid solution via generating functions for the Einstein field equations is given by Herrera et al. [43], where E is arbitrary constants of integration; then the corresponding generating functions are Using the above algorithm the generating functions for the present embedding class one solution (using Eq. (10)): By plugging Eq. (16) into Eqs. (48) and (49) we get the following generating functions z(r ) and : It is clear that the present system can be completely determined by one generating function through e λ(r ) and an additional ansatz in the form of a class one condition.

Physical analysis and discussion
In the present paper we have investigated the nature of anisotropic fluid spheres, by utilizing the Karmarkar condition, which are also known as spacetimes of embedding class one. To outline this investigation we have considered the following assumption for the gravitational potential: (i) Firstly we have developed a relation between the gravitational potentials e ν and e λ for the spherically symmetric line element by using the Karmarkar condition. Due to this relation, we can convert all the differential equations in terms of one of the metric coefficients (the full details can be seen in the references by Maurya et al. [36,37]). For this purpose we have assumed a totally new metric potential e λ = [1 + 2c r 2 + cosh 2(ar 2 +b)]/[1 + cosh 2(ar 2 +b)] to find the anisotropic solution for realistic fluid spheres. The variation of e ν and e λ can be seen in Fig. 1.
(ii) The fluid spheres are purely anisotropic because embedding class one solutions can give only two types of perfect fluid solutions which are namely the Kohlar-Chao or Schwarzschild solutions; otherwise the metric turns out to be flat. As we can see the radial pressure p r is zero at the boundary but the tangential pressure p t does not vanish at r = R; however, both are equal at the centre of the fluid sphere (Fig. 2). Also the pressure anisotropy vanishes everywhere inside the fluid sphere if and only if c = 0. In this situation the pressures and density become zero and the metric turns out to be flat. As we can see in Fig. 3 the anisotropy is increasing throughout, while it is zero at the centre, which implies that p r and p t are not equal except at the centre of the fluid sphere. (iv) Since p r /ρ and p t /ρ lie between 0 and 1 everywhere within the sphere, our fluid sphere satisfies the Zeldovich condition. Moreover, p r /ρ and p t /ρ are monotonically decreasing throughout the inside the sphere (Fig. 4). (v) For the well behaved nature of the solution the velocity of sound should be decreasing throughout the fluid sphere and must be less than the velocity of light. From Fig. 5, it is clear that both v r and v t are decreasing and less than 1, which shows that our anisotropic solution is well behaved. Also all dominant energy conditions are satisfied, a necessary physical requirement everywhere inside the fluid sphere (Fig. 6). (vi) The surface redshift is also determined by using the compactness factor for the fluid sphere. For the fluid sphere PSR J1903+327, the surface redshift turns out to be z s = 0.44505, which is a maximum. The redshift without cosmological constant for isotropic fluid spheres is bound by z s ≤ 2 [29,38,39]. However, in the presence of the cosmological constant Bohmer and Harko [39] argued that the surface redshift must satisfy the restriction z s ≤ 5 for anisotropic fluid spheres. Therefore the value of the surface redshift for our anisotropic spheres seems to be compatible with realistic compact objects. Moreover, the surface redshift cannot be arbitrarily large because it depends on the compactness factor u = M/R. (vii) We have also discussed the stability of the fluid sphere (which is the most vital physical requirement) by using the following facts: (a) the cracking concept proposed by Herrera [32], (b) the variation of the adiabatic index inside the fluid spheres, (c) the stable equilibrium condition by using the TOV equation and (d) the Harrison-Zeldovich-Novikov stability criterion. The results are as follows. It can be observed in Fig. 8 that the velocity v 2 r and v 2 t lie in the ranges: (i) 0.053 ≤ v 2 r ≤ 0.107, 0.017 ≤ v 2