A new model for spherically symmetric anisotropic compact star

In this article we obtain a new anisotropic solution for Einstein’s field equations of embedding class one metric. The solution represents realistic objects such as Her X-1 and RXJ 1856-37. We perform a detailed investigation of both objects by solving numerically the Einstein field equations with anisotropic pressure. The physical features of the parameters depend on the anisotropic factor i.e. if the anisotropy is zero everywhere inside the star then the density and pressures will become zero and the metric turns out to be flat. We report our results and compare with the above mentioned two compact objects as regards a number of key aspects: the central density, the surface density onset and the critical scaling behaviour, the effective mass and radius ratio, the anisotropization with isotropic initial conditions, adiabatic index and red shift. Along with this we have also made a comparison between the classical limit and theoretical model treatment of the compact objects. Finally we discuss the implications of our findings for the stability condition in a relativistic compact star.


Introduction
Recent development in a cosmological deep survey has clarified progressively the origin and distribution of matter and evolution of compact objects in the Universe. Some of their properties, such as masses, rotation frequencies and emission of radiation are measurable, whereas measurements of important parameters which determine the nature of compact stars still represent an observational challenge. The properties that are not directly linked to observations, such as the internal composition or masses and radii, require the development of theoretical models. On the theoretical side, the a e-mail: sunil@unizwa.edu.om b e-mail: kumar001947@gmail.com c e-mail: baiju@unizwa.edu.om d e-mail: saibal@associates.iucaa.in mass and the radius are determined by solving the hydrostatic equilibrium equation which expresses the equilibrium between gravitational and pressure forces. In the framework of general relativity the equilibrium of a spherical object is described by the Tolman-Oppenheimer-Volkoff (TOV) [1,2] equations, and, for completeness, the equation of state is required. Very recently other theoretical advances in modelling of densely neutral gravitating objects in strong gravitational fields have generated much interest in last couple of decades. This is because of its importance in describing relativistic astrophysical objects such as neutron stars, quark stars, hybrid proto-neutron stars, bare quark stars etc.
The main theoretical routes have been used to study of stellar structure and evolution assuming that the interior of a star can be modelled as perfect fluid. The perfect fluid model necessarily requires the pressure in the interior of a star to be isotropic. Spacetime fuelled by a rotating anisotropic fluid has been used to model the interior of the star. One common source is a fluid with anisotropy in pressure. In particular at very high densities conventional celestial bodies are not composed purely of perfect fluids so that radial pressures are different from tangential pressures. The model of Bowers and Liang [3] is conceptually different from isotropic matter, but it possesses anisotropic matter in the study of general relativity. Mak and Harko [4] and Sharma et al. [5] suggest that anisotropy is a sufficient condition in the study of dense nuclear matter with a strange star. Some argument against the existence of anisotropy could be verifiable through the existence of a solid core or the presence of a type 3A fluid; see the work of Kippenhahn and Weigert [6]. On the other hand, this can arise from different kinds of phase transition and pion condensation as pointed out by Sokolov and Sawyer [7]. The structure of compact objects in general relativity will depend on several parameters, including fluid and magnetic stresses, entropy gradients, composition, heat flow and neutrino emission. However, we restrict our attention to the case of an anisotropic perfect fluid with equilibrium composition.
The theoretical investigation of compact objects has been done by several workers by using both analytical and numerical methods. However, emphasis has always been given on the importance of the local pressure anisotropy. This seems to be very reasonable explaining the matter distribution under a variety of circumstances. Also, this has been proved to be very useful to explore characteristics of relativistic compact objects [8][9][10][11][12][13][14][15]. In recent years, many exact solutions to the Einstein field equations have been generated by different approaches [16][17][18][19][20][21][22]. Therefore, the Einstein-Maxwell spacetime geometry for a compact object having a local anisotropic effect has attracted considerable attention in various physical investigations. However, the physical acceptability of the solution depends on the number of criteria which include the fulfillment of various energy conditions of general relativity.
Against this background we would like to mention that in one of the earlier works Maurya et al. [23] have proposed an algorithm of a charged anisotropic compact star while in the later works [24,25] they have given a new approach for finding an anisotropic solution of Einstein's field equations by using the metric potentials function. The present work is a sequel of the work done by Maurya et al. [26,27], in which the authors have obtained the charged compact star and the structure of a relativistic electromagnetic mass model under the condition of a class one metric. It would be desirable to do a systematic stability analysis of our model based on anisotropic spacetime. In this work, we check the massradius relation, stability and surface redshift of our models and find their behaviour is well behaved.
The present article is organized as follows: Sect. 2 contains the spherically symmetric metric and the Einstein field equations. Also we find the metric function λ in terms of the metric function ν by applying the class one condition. In Sect. 3, we obtain the interior structure of the anisotropic star under the class one condition. The values of the arbitrary constants and total mass of the compact star of radius R are obtained by using the boundary conditions in Sect. 4. In Sect. 5, we discuss several required physical conditions for anisotropic models along with the stability analysis which is vital one. Section 6 contains some concluding remarks on the anisotropic models.

Line element for class one metric and Einstein's field equations
We consider the static spherically symmetric metric for describing the spacetime of the compact stellar configuration The energy-momentum tensor of interior matter for a strange star may be expressed in the following standard form: where ρ, p r and p t correspond to the energy density, radial and tangential pressures, respectively, of the matter distribution.
The Einstein field equations can be written as Here G = c = 1 in geometrized relativistic units.
In view of the metric (1), Eq. (3) yields the following differential equations [28]: The metric (1) may represent the spacetime of the embedding class one, if it satisfies the condition of Karmarker [29]. This condition gives the following relation between the metric potentials ν and λ [24]: Here C is a positive constant quantity.

New anisotropic models for compact star
To find an interior solution of the anisotropic compact star in class one, we consider the pressure isotropy condition as given through the expression of the anisotropic factor as follows: For finding the non-zero expression for the anisotropic factor, we assume the metric potential to be in the form where A and B are positive constants. Hence from Eqs. (7) and (9), we get where From Eqs. (9) and (10), we observe that e λ(0) = 1 and e ν(0) = B at the centre, r = 0. This shows that metric potentials are singularity free and positive at the centre. Also both are monotonically increasing functions, which shows that these metric potential are physically valid [30]. These features can be observed from Fig. 1.
By plugging Eqs. (9) and (10) into Eq. (8), we get We note from Fig. 2 that the anisotropy is zero at the centre r = 0 and is monotonically increasing with the increase of r . Also from Eq. (12), we observe that the anisotropic factor vanishes everywhere inside the compact star if and only if A = 0.
Equations (4), (5) and (6) give the expressions for the radial pressure p r , the tangential pressure p t and the energy density ρ: The radial and tangential pressures at the centre, r = 0, can be given by p r = A (4−D)/8 π and p t = A (4−D)/8 π .
Since A and D are positive and the pressure should be positive at the centre, this implies that D < 4. In a similar way, we can find the density at the centre, r=0, as ρ 0 = (3 A D/8 π). Since the density should be positive at the centre, D is positive due to positivity of A. As D, A, B all are positive, C is also a positive quantity. The behaviour of p r and p t are shown in Figs. 3, 4. We suppose that the radial and tangential pressures of the star are related to the matter density by the parameters ω r and ω t as p r = ω r ρ and p t = ω t ρ.
Then the expressions for the parameters ω r and ω t are given by From Fig. 5 it is clear that the ratios ω r = p r /ρ and ω t = p t /ρ are less than 1. This implies that the density dominates over the pressures throughout inside the star. However, this also implies that the underlying fluid distribution is nonexotic in nature [39].

Matching condition
For any physically acceptable anisotropic solution, the following boundary conditions must be satisfied: (i) At the surface of the compact star, the interior of metric (1) for anisotropic matter distribution match with the exterior of Schwarzschild solution [31], which is given by the metric where M is a constant representing the total mass of the compact star at r = R. (ii) The radial pressure p r must be finite and positive at the centre r = 0 and it must vanish at the surface r = R of the star [32]. The condition p r (R) = 0 gives This readily yields the radius R of the compact star as Using the continuity of metric coefficients e ν , e λ and ∂g tt ∂r across the boundary of the star gives the following equations: The following equations with Eq. (19) give the value of the unknowns M, B and C as follows: On the other hand, the value of the constant A can be determined by assuming the density at the surface of the star i.e. ρ s at r = R.

Sound speed
The speed of sound should monotonically decrease throughout from the centre to the boundary of the star and it must be within the range 0 ≤ V i = √ d p i /dρ < 1. It is argued by Canuto [33] that the sound speed should decrease outwards for the EOS with an ultra-high distribution of matter. From Fig. 6, it is clear that the speed of sound is monotonically decreasing outwards.

Energy conditions
The anisotropic fluid must satisfy the following energy conditions: the Null energy condition (NEC), the Weak energy condition (WEC) and the Strong energy condition (SEC). Therefore, the following inequalities should hold simultaneously at each point inside the compact star corresponding to the above conditions (Fig. 7):

Tolman-Oppenheimer-Volkoff equation (TOV)
The generalized TOV equation for the anisotropic fluid distribution is given by [1,2] We can write the above TOV equations as follows: where M G is the effective gravitational mass and it can be given by Equation (32) describes the equilibrium condition for an anisotropic fluid distribution subject to the gravitational (F g ), the hydrostatic (F h ) and the anisotropic stress (F a ) so that where its components can be defined as The explicit form of the above forces can be expressed as (Fig. 8)

Herrera cracking concept
We know that for physically acceptable anisotropic models, the radial and transverse speed of sound should lie between 0 and 1, i.e., 0 ≤ V r < 1 and 0 ≤ V t < 1. We observe from this inequality that the parameters also should satisfy the inequality 0 ≤ V 2 r < 1 and 0 ≤ V 2 t < 1. Now we define the expression for the square of the velocity of sound as From Fig. 9, we conclude that the square of the radial and transverse speeds of sound are within the range everywhere inside the stars. Therefore, 0 ≤ |V 2 t − V 2 r | < 1. In order to examine the stability of the local anisotropic fluid distribution, we follow the cracking concept of Herrera and Aberu et al. [34,35] which states that the region is potentially stable where the radial speed of sound is greater than the transverse speed of sound. This implies that there is no change in sign V 2 r − V 2 t and V 2 t − V 2 r . So we calculate the difference between the radial and transverse speeds of sound: D e 2 Ar 2 [D 2 Ar 2 e 4 Ar 2 − 2 (5 + 4 Ar 2 ) + D e 2 Ar 2 (5 + 6Ar 2 + 8A 2 r 4 ) .
We note from Fig. 10 that the radial speed of sound is always greater than the transverse speed of sound and also 0 ≤ |V 2 t −V 2 r | < 1 everywhere inside the star. These features represent that the proposed physical models are stable.

Adiabatic index
In order to determine an equilibrium configuration, the matter must be stable against the collapse of local regions. This also requires Le Chatelier's principle (known as the local or microscopic stability condition), stating that the radial pressure p r must be a monotonically non-decreasing function of ρ such that d p r dρ > 0 [36]. Heintzmann and Hillebrandt [37] also proposed that a neutron star with an anisotropic equation of state is stable for γ (= p r +ρ p r d p r dρ ) > 4/3. From Fig. 11, it is clear that the adiabatic index (γ ) is higher than 4/3 everywhere inside the star.

Effective mass-radius ratio
This section contains the maximum allowable mass-radius ratio for the above proposed anisotropic fluid models. As Buchdahl [38] has already discussed, the maximum limit of mass-radius ratio for a static spherically symmetric perfect fluid star should satisfy the upper bound 2M/R < 8/9. Also Mak and Harko [4] have given the generalized expression for the same mass-radius ratio.
The effective mass of the anisotropic compact star is defined as Fig. 9 Variation of square of radial speed of sound and transverse speed of sound with radial coordinate r/R for Her X-1 and RXJ 1856-37 However, the compactness u of the star can be expressed as (45)

Surface redshift
The surface redshift (Z ) corresponding to the above compactness (u) is given by the expression (Fig. 12)

Conclusions
In the present article, we have obtained new anisotropic compact star models of the embedding class one metric. Our models satisfy all the physical reality conditions. Some of the special features of the present model are as follows: 1. We used the boundary conditions by joining the Schwarzschild metric with a class one metric at the boundary of the star r = R. Subsequently we obtained the arbitrary constants A, B, C along with the total mass   of the compact star and the corresponding numerical values are provided in Table 1. All these values match the observed data of real compact stars. 2. The metric potentials are free from any singularity at the centre, and positive and finite inside the star (Fig. 1). Also the ρ, p r and p t are positive, finite and monotonically decreasing away from the centre. However, the parameters ω r and ω t are within the range between 0 and 1 (Fig. 5). 3. The model is in static equilibrium. We observe from Fig. 8 that the gravitational force (F g ) is dominating over the hydrostatic force (F h ) and is counter balanced by the joint action of the hydrostatic force and the anisotropic stress. 4. The model has a density of the order 10 15 gm/cm 3 . The corresponding values for Her X-1 and RXJ 1856-37 are as follows: (i) at the centre ρ 0 = 1.8664 × 10 15 gm/cm 3 and ρ 0 = 2.3968 × 10 15 gm/cm 3 , (ii) at the surface ρ s = 1.3273 × 10 15 gm/cm 3 and ρ s = 1.6924 × 10 15 gm/cm 3 ( Table 2).
This density profile shows that our models may represent a realistic anisotropic objects. (5) The redshift is monotonically decreasing and attains its maximum value at the centre of the compact star.
The numerical values corresponding to the Her X-1 and RXJ 1856-37 are: (i) at the centre Z 0 = 0.2669 and Z 0 = 0.2796, (ii) at the surface Z s = 0.0474 and Z s = 0.0480. As a final comment, an interesting and puzzling point about the anisotropic compact model is that its stability depends on the unavoidable anisotropic pressure and the TOV equations used to place a constraint on the anisotropic parameters. It would be interesting to propose a richer model in which consideration of the pressure anisotropy on the compact relativistic objects could lead to a more realistic model of the anisotropization mechanism as regards compact relativistic objects.