Compact star models in class I spacetime

In the present article, we have presented completely new exact, finite and regular class I solutions of Einstein’s field equations i.e. the solutions satisfy the Karmarkar condition. For this purpose needfully we have introduced a completely new suitable grr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{rr}$$\end{document} metric potential to generate the model. We have investigated the various physical aspects for our model such as energy density, pressure, anisotropy, energy conditions, equilibrium, stability, mass, surface and gravitational red-shifts, compactness parameter and their graphical representations. All these physical aspects have ensured that our proposed solutions are well-behaved and hence represent physically acceptable models for anisotropic fluid spheres. The models have satisfied causality and energy conditions. The presented models are also stable by satisfying Bondi condition and Abreu et al. condition, in equilibrium position and static by satisfying TOV equation, Harrison–Zeldovich–Novikov condition, respectively. For the parameters chosen in the paper are matching in modeling Vela X-1, Cen X-3, EXO 1785-248 and LMC X-4. The M–R graph generated from the solutions is matching the ranges of masses and radii for the considered compact stars. This work also estimated the approximate moment of inertia for the mentioned compact stars.


Introduction
After the pioneering work on the anisotropic relativistic compact star by Bowers and Liang [1], many researchers have investigated the possible origin of anisotropy in compact stars. Ruderman [2] proposed that at high density ∼ 10 15 g/cc the matter starts interacting relativistically that arises the anisotropy in pressure. It is also suggested that due to the a e-mail: nayan.mathju@gmail.com b e-mail: ntnphy@gmail.com c e-mail: susmita.mathju@gmail.com d e-mail: rahaman@iucaa.ernet.in formation of super-fluid neutrons inside neutron stars, pressure anisotropy may also arise [3]. As a results of research by many authors pressure anisotropy can be trigger by types of phase transition [4], pion-condensation [5], slow rotation [6], strong magnetic field [7] etc. Letelier and his co-author [8][9][10] have shown that such anisotropic matters can be considered as the composition of two perfect fluids, or a perfect fluid and a null fluid, or two null fluids.
Many researchers are also interested to investigate the properties of compact stars in higher dimensions. The concepts of extra dimensions were first proposed by Kaluza [11] and Klein [12] independently when they unify the gravity and EM-force. Liddle et al. [13] have analyzed the effect of extra dimensions in the maximum mass of neutron stars (NS). They have assumed an equation of state for non-interacting cold neutrons and found that the maximum mass of NS was reduces by the presence of extra dimensions. Chattopadhyay and Paul [14] have considered the Vaidya-Tikekar spacetime in n-dimensional Einstein's field equations to analyze the properties of compact stars. Bhar et al. [15] incorporated the conformal Killing equations in n-dimensional Einstein's field equation and discussed anisotropic compact stars.
Various extended theories of gravity have been used to investigate the physical properties of compact stars. Pani et al. [16] have used general class of alternative theories which includes scalar-tensor theories, a scalar field coupled to quadratic curvature invariant and indirectly f (R) theories as special cases. In their work, they have developed a systematic tool to rule out physical theories that are incompatible with observations. The concept of embedding 4-dimensional spacetime in 5-dimensional hyperspace was used by Castro et al. [17]. Here they have embedded four dimensional spacetime into five dimensional braneworld. It is shown that on choosing equation of state for hadronic, hybrid and quark stars, the maximum mass is controlled by brane tension λ which lies in the range 3.89 × 10 36 (≡ 1.44M ) < λ < 10 38 dyne/cm 2 . It is important to remind that construct-ing exact interior solutions representing non-uniform stellar distributions is nearly impossible in the context of the braneworld. This is because the nonlocality and non-closure of the braneworld equations, produced by the projection of the bulk Weyl tensor on the brane, lead to a very complicated system of equations which make the study of non-uniform distributions very hard [18]. It is yet to discover the criteria about what restriction should be imposed on braneworld equations to obtain a closed system [19]. To solve this problem, it is necessary to understand the bulk geometry and how a 4D spacetime cab be embedded. On the other hand, Karmarkar [20] embedded 4-dimensional spacetime into 5dimensional Euclidean space known as class I. This method implies an equation that links the metric coefficients g tt and g rr thereby simplifying to solve the field equations. A similar concept was also used in string theory on embedding branes e.g. in the Randall-Sundrum model [21].
Assuming anisotropic fluid distribution, Mak and Harko [22] derived the condition for lower limit of mass for any anisotropic compact stars which was strongly depends on degree of anisotropy. Dev and Gleiser [23] derived the critical condition for compactness parameter 2M/R for anisotropic relativistic stars and was also found strongly depends on nature and degree of anisotropy. They have also shown that for anisotropic compact stars the red-shift can go arbitrary large. In there second paper, Dev and Gleiser [24] stressed on stability of anisotropic compact stars. Here they have found that, depends on the nature of anisotropy there can exist stable anisotropic star even at < 4/3 while its isotropic counterpart isn't. Herrera et al. [25] obtained an algorithm to generate all spherically symmetric anisotropic fluid distributions from two generators. One of the generators was linked with anisotropy and other with re-shift function. On the other hand, Lake [26] used Newtonian hydrostatic equation for isotropic fluid distribution and generate infinite number of anisotropic solutions simply by assuming density.
Recently, Ivanov [27] derived a condition which is similar to Karmarkar condition, for conformally flat spacetime. However, the solutions resulting from the two theories are completely different. In Karmarkar spacetime there are no physical solutions describing isotropic fluid distributions, however, physical solution exist if electric charge or anisotropy or both are incorporated [28][29][30][31][32][33][34][35][36][37][38][39][40][41][42]. In this article also we are exploring new physical solutions satisfying field equations under class I category and discuss the solutions to model compact stars.
The present article has been designed as follows: In Sect. 2, we have written the Einstein field equations for static and spherically symmetric matter distribution. The Embedding class I spacetime satisfying the Karmarkar condition has described in Sect. 3. Section 4 contains our new class I solutions along with mass, compactness parameter, surface redshift and gravitational redshift. We have explored the accep- Metric potential functions (above) and energy density (below) vs the radial coordinate r corresponding to numerical values of constants given in Table 2 for four well-known compact stars tance of our solutions in Sect. 5. The determination of constants by using the matching condition and analysis of all energy conditions have been done in Sects. 6 and 7, respectively. The equilibrium and stability are analyzed in the subsections of Sect. 8: Sect. 8.1 displayed the equilibrium situation and the stability analysis has done in Sect. 8.2. The generating functions for our model have calculated in Sect. 9. Moment of inertia and Mass relationship are analyzed in Sect. 10. Finally, The result and discussions have been done in Sect. 11.

Einstein's field equations
To describe the interior of a static and spherically symmetric fluid sphere we consider the line element in Schwar zschild coordinate system (t, r, θ, φ) as: where e ν(r ) and e λ(r ) are unknown functions of the radial coordinate r only and called metric potential functions.
As we are going to investigate the solutions of Einstein's field equations for anisotropic compact stars, so we consider the energy momentum tensor for the anisotropic fluid sphere, which is of the following form: (2) where ρ(r ), p r (r ) and p t (r ) are representing the energy density, radial pressure and transverse pressure of the fluid configuration, respectively. χ α and U α are the unit space-like vector and four velocity, respectively, satisfying Therefore, the Einstein field equations for the line element (1) and energy momentum tensor (2) are as follows: where = d dr , represents the derivative with respect to the radial coordinate r .
The anisotropic factor is defined as (r ) = p t (r )− p r (r ). Therefore, with the help of Eqs. (4) and (5) we can get the general expression of anisotropic factor for static and spherically symmetric fluid configuration in the following form:

The Karmarkar condition
The K armarkar condition is an important tool to study the stellar fluid spheres because of its special character. In general theory of relativity, it is recognized that an n-dimensional spacetime is called of class p if it is embedded in (n + p)dimensional Pseudo Euclidean flat space. In the year 1921, Kasner [43] investigated that the 4-dimensional spacetime of spherically symmetric object can always be embedded in 6-dimensional Pseudo Euclidean space and later Gupta and Goyel [44] have shown the same result with respect to another coordinate transformation. In 1924, Eddington [45] found that an n-dimensional spacetime can always be embedded in m-dimensional Pseudo Euclidean space with  Table 2 for four well-known compact stars m = n(n + 1)/2 and the required minimum extra dimension to embedded is less than or equal to the number (m − n) or same as n(n −1)/2. Therefore, the 4-dimensional spherically symmetric line element (1) is of embedding class II. In the literature, there are some special class spacetimes such as the Schwar zschild interior and exterior solutions are of class I and class II, respectively, Friedman-Robertson-Lemaitre [46][47][48] spacetime is of class I and the Kerr metric is of class V [49]. In the year 1948, K armarkar derived a condition [20] in terms of the components of Riemannian curvature tensor as: This condition is known as the K armarkar condition. Any 4-dimensional spacetime can be embedded in 5-dimensional flat space i.e. becomes an embedding class I whenever it satisfies the K armarkar condition. The K armarkar condition is only the necessary condition to become a class I, Pandey and Sharma [50] provided the sufficient condition as R 2323 = 0. Now, all the non-zero components of Riemannian curvature tensor R αβγ δ for the metric (1) are : Taking into account all these non-vanishing components of Riemannian curvature tensor, the K armarkar condition (7) yields the following differential equation:

Fig. 3
Anisotropic factor (above) and adiabatic index (below) vs the radial coordinate r corresponding to numerical values of constants given in Table 2 for four well-known compact stars After solving the above differential Eq. (8) we get a relationship between two metric potential functions where A and B are non-zero unknown integration constants, which shall be determined by applying some boundary conditions at the surface of compact stars. The condition (9) is a peculiar characteristic of the class I spacetime, where two metric potential functions are related to each other and hence researches can easily generate class I models for anisotropic fluid spheres with the help of a suitable form of one of the metric potential functions. Therefore, the expression of anisotropic factor given in Eq. (6) becomes [51]: by using Eqs. (8) and (9). From Eq. (10) one can see that the pressure anisotropy (r ) is zero throughout the fluid sphere if either first or sec-  Table 2 for four well-known compact stars ond or both the factors on right-hand side of Eq. (10) are zero. The vanishing of first factor on the right side of Eq. (10) yields the Kohlar-Chao solution [52] whereas if the second factor is zero then the corresponding solution will be the Schwar zschild interior solution [53]. Finally, we arrive at a point to find out class I solutions but we have four independent equations (3)- (5) and (9) with five unknowns, namely λ(r ), ν(r ), ρ(r ), p r (r ) and p t (r ). Therefore, it is impossible to find the exact solutions, for this purpose we shall consider a new metric potential function to generate our models.

New embedding class I solutions
According to earlier information, we consider a completely new metric potential function to obtain the closed form solutions of Einstein's field equation: where c > 0 (km −2 ), a = 0 (km −2 ) are undetermine constants, to be determine from boundary condition and erf[x] is known as the err or function, which is defined as We get another metric potential function by employing Eq. (11) in Eq. (9) as: Now, taking into account these two metric potential functions given in Eqs. (11) and (13), we obtain the exact expressions of ρ(r ), p r (r ), p t (r ) and (r ) as: whereas Therefore, the gradients of energy density ρ(r ) and radial pressure p r (r ) are obtain as whereas, Mass function (above) and compactness parameter (below) vs the radial coordinate r corresponding to numerical values of constants given in Table 2 for four well-known compact stars We get the following expressions for the mass function m(r ) and compactness parameter u(r ), respectively Also, the interior redshift Z g (r ) and surface redshift Z s (r ) are vs the radial coordinate r corresponding to numerical values of constants given in Table 2 for four well-known compact stars To explore the behaviors of the mass, compactness parameter, surface and gravitational redshifts we provide the graphical representations of them in Figs. 5 and 7. The mass and compactness parameter are zero at the centres of the fluid spheres, gradually increasing towards the surfaces of the fluid configurations and become maximum at the surfaces (see Fig. 5). The surface redshift and gravitational redshift have the opposite behavior throughout the fluid spheres, shown in Fig. 7. Moreover, at the surfaces (r = R) Z (R) = Z s (R) = Z g (R), obvious from Fig. 7. The the behavior of interior redshift explain the profile of interior density as well. If a photon comes out from center to surface, it has to travel longer path and much denser region (i.e. core). This leads to more dispersion resulting into loss of energy. Whereas a photon comes out from near the surface will travel shorter path and less denser region and therefore less dispersion and less energy loss. Hence, the interior redshift is maximum at the center and minimum at the surface. However, the surface redshift depends on overall mass and radius or in other word the surface gravity. As mass increases the radius will also increases slightly which will yields more surface gravity and more surface redshift. Thus, the trend of interior and surface redshifts are opposite.

The central values and physical analysis
For our model, we get the following results for the density and pressure at the centres of compact stars:  Table 2 for four well-known compact stars where g 0 = g(0) = erf [1]. The positive behavior of the central pressure ( p r c = p tc > 0) implies Moreover, according to the Zeldovich condition [54] p r c /ρ c ≤ 1. Therefore, from this condition we obtain the following inequality: Therefore, the inequalities (26) and (27) yield a boundary restriction on the constants A and B as: To demonstrate the physically acceptance of our models on anisotropic compact stars we draw the graphs for all the physical parameters, metric potential functions e −λ(r ) , e ν(r ) , energy density ρ(r ), radial pressure p r (r ), transverse pressure p t (r ), anisotropic factor (r ) and equation of state (EoS) parameters ω r (r ), ω t (r ) in Figs. 1, 2, 3 (above) and 4, respectively. Figure 1 (above) shows that e −λ(r ) and e ν(r ) are finite and non-singular everywhere within the fluid spheres and they coincide at the surfaces. The energy density is positive and as usual maximum at the centres and then decreasing towards the the boundaries, clear from Fig. 1 (below). The radial and transverse pressures both are positive and monotonically decreasing from the centre of the stars. Further, p r (r ) becomes zero at the boundaries of the stellar objects and p t (r ) is non-vanishing there (see Fig. 2). The anisotropic factor (r ) is positive and finite within the stars, maximum and zero at the boundaries and centres of the stars, respectively, obvious from Fig. 3 (above). The parameter η(r ) = p r (r ) p t (r ) can analyze the effect of anisotropy in the equilibrium position of fluid configurations in similar way of Newtonian treatment. The profile of η(r ) (see Fig. 10 (below)) discloses that the anisotropic force is increasing through out the fluid spheres i.e. behaves like outward force, similar to the result in Fig. 3 (above). Figure 4 indicates the behaviors of EoS parameters ω r (r ) and ω t (r ) and from the figure one can notice that both the parameters are with in the region 0 < ω r (r ), ω t (r ) < 1 i.e. our solutions are for the real feasible fluid distributions. Therefore, all this significant results confirm that our solutions are physically acceptable.

The matching condition and determination of constants
We have seen that there are some constants within our solutions, so to determine the values of these involving constants we are going to match our interior solution with the Schwar zschild vacuum solution at the surface of the star r = R. The Schwar zschild vacuum solution is given by the following metric where M is the total mass of the anisotropic fluid sphere contained within the sphere of radius R. By matching at the boundary of the compact star r = R (> 2M to avoid the singularity) we get: Also, the radial pressure p r (r ) at the boundary r = R vanishes i.e.
Using these boundary conditions (30)-(32), we obtain where a is a free parameter to obtain well-behaved solutions in all respects and mass M, radius R will be chosen accordingly different stars.

The energy conditions
It is well-known in the literature that the physical mass distributions must satisfy all the energy conditions within its interiors. The energy conditions are: (2) null energy condition (NEC), (2) weak energy condition (WEC) and (3) strong energy condition (SEC), these conditions represent by following inequalities: Here i ≡ (r, t), r for radial and t for tranverse components.
For the verification of all energy conditions, we plot all the L.H.Ss of above inequalities in Fig. 6. From Fig. 6 it is evident that our solutions satisfy all the energy conditions within the interiors of fluid configurations.

The equilibrium and stability analysis
The equilibrium position and stable state are most important situations for the non-collapsing compact object within our Universe. In this section, we are going to check the equilibrium and stability of the fluid distributions represented by our solutions.

The equilibrium condition
The together effect of gravitational force, Hydrostatic force and Anisotropic force holds any anisotropic stellar object in the equilibrium position. At the equilibrium position, the balancing force equation is known as T OV -equation. The generalized T OV -equation for anisotropic fluid distribution can be written as [55,56] Fig. 8 Forces for the compact stars Vela X-1 (above left), Cen X-3 (above right), EXO 1785-248 (below left), LMC X-4 (below right) vs the radial coordinate r corresponding to numerical values of constants given in Table 2 − M g (r ) where M g (r ) represents the effective gravitational mass. The Tolman-Whittaker mass formula gives the exact form of the effective gravitational M g (r ) as: On using the energy-momentum tensor (2) and field equations (3) and (4), the Eq. (39) becomes Therefore, using the value of M g (r ), Eq. (38) reduces as: Moreover, the Eq. (41) can be written as: where F g (r ) = − 1 2 ν (r ){ρ(r ) + p r (r )}, F h (r ) = − dp r (r ) dr and F a (r ) = 2 (r ) r are called the gravitational force, hydrostatic force and anisotropic force, respectively. For our present solutions the three different forces are obtained in the following forms: whereas, S 13 = Bπ e 2h 2 r + 2Bcr 3 + Ar π e 2h 2 √ cg and S 1 , S 2 ,S 6 , S 7 , S 8 , S 9 , S 10 , S 11 , S 12 are given in Eqs. (16) and (19), respectively.
The behaviors of these three forces are shown in Fig. 8 and it is ascertained that our solutions representing fluid distributions are in equilibrium positions.

The stability condition
Here, we are willing to analysis the stability condition for our models with the help of (1) stability factor, (2) Adiabatic index and (3) Harrison-Zeldovich-Novikov criterion.

Causality condition
In General Theory of Relativity, the maximum velocity is the velocity of light, which is equal to 1 in the gravitational unit.
1. Causality condition: The causalit y condition states that whenever sound passes through the physical fluid distribution then its velocity must be less than the velocity of light, otherwise non-physical stellar configuration. The radial velocity v r (r ) and transverse velocity v t (r ) of sound inside the compact star can be determined using the following formulae: v r (r ) = dp r (r ) dρ(r ) and v t (r ) = dp t (r ) dρ(r ) . (46) Therefore, according to causalit y condition 0 ≤ v r (r ), v t (r ) < 1. Figure 9 reveals that our models satisfy the causalit y condition i.e. our models represent the physical stellar fluids.

Stability condition:
In the year 1992, Herrera [57] proposed the cracking method to study the stability of an anisotropic stellar fluid under the radial perturbations. Later, using the concept of cracking Abreu et al. [58] provided the conditions for anisotropic fluid model with respect to the stability factor {v t (r )} 2 − {v r (r )} 2 as: Our solutions satisfy the condition −1 < {v t (r )} 2 − {v r (r )} 2 < 0, clear from Fig. 10 (above) and hence our models of anisotropic compact stars are potentially stable.

Adiabatic index
The relativistic adiabatic index plays an important role to analyze the stability of stellar fluid spheres. The relativistic adiabatic index r (r ) is defined as: Fig. 9 Radial sound speed and transverse sound speed vs the radial coordinate r corresponding to numerical values of constants given in Table 2 for four well-known compact stars For Newtonian limit, any stable configuration will alter its stability by initiating an adiabatic gravitational collapse if r (r ) ≤ 4/3 and catastrophic if < 4/3 [59]. According to Chan et al. [60] this condition changes for relativistic and/or anisotropic fluid which depends on the nature of anisotropy. Figure 3 (below) shows that the adiabatic index r (r ) > 4/3 for our solutions with positive anisotropy. More strict condition on adiabatic index for stable region was derived by Moustakidis [61] and found that the critical value of adiabatic index crit depends on ξ -parameter (amplitude of lagrangian displacement from equilibrium) and compactness parameter β = M/R. On assuming particular form of ξ -parameter he obtained the constraint as Any stable configuration should have ≥ crit . For the presented this condition is fulfilled (Fig. 3 below).

Harrison-Zeldovich-Novikov criterion
The Harrison-Zeldovich-Novikov [54,62] static stability criterion states that the mass should increase with the increase of central density ρ c for the stable state of compact stars, i.e. ∂ M(ρ c ) ∂ρ c > 0 for the stable state of compact stars. For our solutions, we obtain the mass as a function of the central density as: The obtained solutions hold static stability criterion and hence stable, clear from Fig. 11.  Fig. 10 Stability factor (above) and η(r ) (below) vs the radial coordinate r corresponding to numerical values of constants given in Table 2 for four well-known compact stars  Table 2 for four well-known compact stars

Generating functions
To obtain all possible anisotropic solutions of Einstein's field equations, Herrera [6] proposed an algorithm with the help of generating functions as:  Table 2 for four well-known compact stars +2Z dr dr (50) where, C is an arbitrary integration constant and the corresponding generating functions are: Now, adopting the class I condition (9), the generating functions in Eqs. (51) and (52) take the following forms: The generating functions for our present solutions are obtained as: where (r ) is given in Eq. (15).

Moment of inertia and mass relationship
Bejger and Haensel [67] adopted a method where a static solution can become rotating by using an approximate expression of moment of inertia I given as  Table 2 for four well-known compact stars Using the above expression we have plotted the trend of I w.r.t. mass M in Fig. 13. From this graph and using the observed ranges of masses in Table 1 we have predicted the possible approximate moment of inertia (see the last section). It is to be noted that I -M graph is more sensitive to the equation of state as compared to that of M-R graph.

Result and discussions
In the present article, we have developed few models for anisotropic compact stars in the framework of K armarkar condition by introducing a completely new metric potential function e λ(r ) . To demonstrate the physical acceptance of our proposed models, we have performed varies physical experiments with the help of different physical parameters. For clarity, we have provided the graphical presentations (Figs. 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 and 13) for the physical parameters involved in our models by matching our interior solution with the exterior Schwar zschild solution for four well-known compact stars. All the features of our model are following: • Metric potentials: The finite and nonsingular metric potential functions are necessary to generate the physical viable models of anisotropic compact stars. In our models, e λ(0) = 1 and e ν(0) = 1 4a 2 π 2 Aa √ π + B √ c e −1 + √ π g 0 2 = positive constant for nonzero a. i.e. the metric potential functions are well-behaved at the centre of the stars. The graphical representations of e −λ(r ) and e ν(r ) indicate that they are finite and regular throughout the radius of stars (see Fig. 1 (above)) and hence they are suitable to generate the models for anisotropic compact stars. • Energy density and pressures: The energy density ρ(r ), radial pressure p r (r ) and transverse pressure p t (r ) should remain positively finite throughout the interior of fluid spheres. The energy density and radial pressure are maximum at the centre and decreasing in nature towards the surface. Also, the radial pressure should vanish at the surface of fluid sphere. Figures 1 (below) and 2 approve that our obtained energy density and pressures (radial and transverse) are good in behavior as they have satisfied all those conditions. To compare with observational data we have calculated the numerical values of central, surface densities and central pressure for different compact stars given in Table 3. The central and surface densities both are of order 10 14 and the central pressure is of order 10 34 , which are almost same with observational data. • Equation of state parameters: For real matter distribution the equation of state(EoS) parameters ω r (r ) = p r (r ) ρ(r ) and ω t (r ) = p t (r ) ρ(r ) should lie in 0 < ω r (r ), ω t (r ) < 1 [63]. In our present models, both the EoS parameters are within the region 0 < ω r (r ), ω t (r ) < 1 (see Fig. 4), which is another important testimony of our well-behaved models.
• Anisotropy: The anisotropic factor (r ) = p t (r )− p r (r ) should vanishes at the stellar centre. If (r ) > 0 then the anisotropic force F a (r ) = 2 (r ) r is outward directed i.e. play as a impulsive force, which can support more compact construction. And if (r ) < 0 then the anisotropic Table 3 Estimated values of the central and surface densities, cental pressure, surface redshift at the boundary and twice of compactness parameter with buchdahl limit for four celestial compact stars corresponding to the values of constants given in Table 2 Stars ρ c (10 14 ) ρ s (10 14  F a (r ) is inward directed. In our models, (0) = 0 and positively increasing within the stellar interiors, clear from Fig. 3 (above) and consequently the anisotropic force is repulsive in nature, assist to construct more compact stars. The similar result is seen in the behavior of η(r ), it indicates the anisotropy as a outward force.  Fig. 5 for four compact stars. The m(r ) and u(r ) tends to zero when r tends to zero and monotonically increasing toward the surfaces. According to Buchdahl[] the mass to radius ration M R < 4 9 or equivalently 2u s = 2u(R) < 8 9 . We have computed the numerical values 2u s , provided in Table 3 and all these values indicate that our solutions satisfied the Buchdahl limit. (above) indicate that the sound velocities (radial v r (r ) and transverse v r (r )) are positive and less than 1 and the stability factor {v r (r )} 2 −{v r (t)} 2 is negative i.e. our solutions represent physical matter distributions, which are potentially stable. The profiles of adiabatic index (Fig. 3 (below)) and mass in terms of central density (Fig. 11) are another two evidences for stable configuration represented by our solutions. • M-RGraph: The profile of mass-surface radius relationship is shown in Fig. 12. It is found that for the ranges of masses given in Table 1 have the same corresponding radii provided on the same table. These means that our solutions yield the ranges of masses and radii very closed to that the observed values. Therefore, the solutions may represent the chosen compact stars.
• Surface redshift and gravitational redshift: The variations of surface redshift and gravitational redshift are shown in Fig. 7. From that figure, we can see that surface redshift Z s (r ) → 0 as r → 0 and thereafter monotonically increasing unto the surfaces of stellar spheres. The gravitational redshift Z g (r ) has opposite behavior, it is maximum at the centres and decreasing toward the surfaces of the fluid configurations. Moreover, Z s (r ) and Z g (r ) are coincided at the surfaces r = R, i.e. Z (R) = Z s (R) = Z g (R). Further, we have calculated the numerical values of surface redshift Z s (r ) at the boundaries of the celestial stars i.e. the maximum values, given in Table 3 and all these maximum values of Z s (r ) are within the range suggested by Ivanov [68].
Finally, with respect to all obtained significant results, we can conclude that our models are physically acceptable to describe celestial anisotropic compact stars.