Minimally deformed anisotropic model of class one space-time by gravitational decoupling

In this article, we have presented a static anisotropic solution of stellar compact objects for self-gravitating system by using minimal geometric deformation techniques in the framework of embedding class one space-time. For solving of this coupling system, we deform this system into two separate system through the geometric deformation of radial components for the source function λ(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda (r)$$\end{document} by mapping: e-λ(r)→e-λ~(r)+βg(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{-\lambda (r)}\rightarrow e^{-\tilde{\lambda }(r)}+\beta \,g(r)$$\end{document}, where g(r) is deformation function. The first system corresponds to Einstein’s system which is solved by taking a particular generalized form for source function λ~(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{\lambda }(r)$$\end{document} while another system is solved by choosing well-behaved deformation function g(r). To test the physical viability of this solution, we find complete thermodynamical observable as pressure, density, velocity, and equilibrium condition via. TOV equation etc. In addition to the above, we have also obtained the moment of inertia (I), Kepler frequency (v), compression modulus (Ke\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_e$$\end{document}) and stability for this coupling system. The M–R curve has been presented for obtaining the maximum mass and corresponding radius of the compact objects.


Introduction
General Relativity (GR) theory of gravitation has been established by Einstein in 1916 in which gravitational properties spread with the speed of light and the law of physics articulated to be invariant with respect to accelerated observers [1][2][3][4][5]. The main assumptions of Einstein's theory of gravitation are based on the (i) all events in the universe as a a e-mail: ntnphy@gmail.com b e-mail: sunil@unizwa.edu.om c e-mail: mahmoodkhalid@unizwa.edu.om d e-mail: rahaman@associates.iucaa.in 4-dimensional Riemannian manifold, which is called spacetime, (ii) the curvature related with the metric is related to the matter by Einstein's field equations (EFE). There will be various fields on the space-time which describe the matter content the space-time through energy-momentum tensor T μν . The field equation of GR is non-linear 2nd order partial differential equation of hyperbolic type, which permits clear freedom of change of coordinates. The first solution of EFE describing a self-gravitating, bounded object was obtained by Karl Schwarzschild [6]. This interior solution represents a constant density model with the outer space-time being empty. However, the velocity of sound within the sphere exceeds the velocity of light thus such a model is not realistic. Therefore, this encouraged us to search for physically viable solutions to the Einstein field equations which describe the realistic models. The space-time of Schwarzschild's exterior solution [7] was obtained in 1916 which describes the gravitational field outside a spherical mass by imposing spherical symmetry on the space-time manifold. The theory of spherical symmetric space-time has been investigated by Takeno (1966) from the point of view of invariant classification, group of motion, conformal transformation and embedding classes, etc.
The method of gravitational decoupling by Minimal Geometric Deformation (MGD) is a great and powerful technique that extends known solutions into more difficult situations. By using this MGD technique, In Ref. [8], Gabbanelli et al. have extended isotropic Durgapal-Fuloria solution in the anisotropic domain while Ovalle and his collaborators [9] have shown that how a spherically symmetric fluid modifies the Schwarzschild vacuum solution and necessity of anisotropy in the fluid. In this connection, several other authors have used the MGD approach to discover the more complex solution which can be seen in the following Refs. [10][11][12][13][14][15][16][17][18][19][20][21][22][23]. The GD was developed by Ovalle [25,26] as a consequence of the Minimal Geometric Deformation (MGD) [24,27] (see also Ref. [28]) in the framework of Randall-Sundrum gravity [29,30]. In Ref. [28], the author extended the Minimal Geometric Deformation approach to investigate a new black hole solution.
The key features of this approach for new solutions to EFE are available in literature as [31]: I. Considering the energy-momentum tensorT μν for known metric as a source and extend known solutions of EFE into more complex situations. The new source is coupled with theT μν associated with the seed solution through a nondimensional coupling constant β which can be written as: T μν →T (1) μν =T μν + β (1) T (1) μν (1) and then continue to same the process with more sources, likê T (1) μν →T (2) μν =T (1) μν + β (2) T (2) μν (2) and so on. In this approach, we can spread direct solutions of the Einstein equations related with the simplest gravitational sourceT μν into the province of more complex forms of gravitational sourcesT μν =T (n) μν . II. It is noted that we can also use the reverse of the above methodology in order to find an exact solution to Einstein's field equation. In this procedure, we can split a more difficult energy-momentum tensorT μν into simpler components, saỹ T (i) μν , and then Einstein's equations have to be solved for each one of these components. In this situation, there will be many solutions corresponding to each component ofT (i) μν associated with the original energy-momentum tensor. At last, we can find the solution of Einstein's equations for the original energy-momentum tensorT μν by combining all the above individual solution. However, we would like to mention that this procedure works very well as each source satisfies the conservation equation identically, which can be written as Now we explain the procedure to explore MGD-decoupling methodology which is as follows: Suppose we have two gravitational sources namely R and S where we will first solve the standard Einstein's equations corresponding to the source R and then the other set of quasi-Einstein equations are solved for the source S. At last, we combine these two solutions to determine the complete solution for the total system of R ∪ S. As we know that Einstein's field equations are non-linear, therefore the above procedure leads a powerful technique to find the solutions and their analysis, especially during the situations that away form trivial cases. Many analytical solutions were created by Tolman [32], which describe the structure of the interior stellar geometry for the perfect fluid models. However, the anisotropic models where the tangential and radial pressures are unequal, allows a better understanding of the highly-dense matter. Ruderman [33] and Bowers and Liang [34] have been studied the anisotropic fluid distribution that has explored the most updated research. In this continuation Mak and Harko [80] has suggested that anisotropy plays an important role to understand the variation of properties of the dense nuclear matter for a strange star. On the other hand, the presence of anisotropy could be identified through the existence of a solid core or type 3A superfluid [81], pion condensed phase [82], and different kinds of phase transitions [83]. The positive anisotropy inside star provides a realistic star which have been studied by several authors. If we have the anisotropy factor Δ(r ) = p t − p r > 0, then anisotropic force inside the stellar system outward-directed which improve the stability and equilibrium criteria, and if the anisotropy factor Δ(r ) = p t − p r < 0, then anisotropic force inside the stellar system is directed inward that introduce instability in the system. Now the finding of the solutions to EFE is a great challenge to meet the requirements of physical acceptability [35][36][37]. Recently Jasim et al. [38] have constructed an anisotropic fluid sphere model by supposing a specific form of the potential metric functions e λ and e ν to EFE with MIT bag EOS in presence of cosmological constant. This model yields a realistic fluid sphere such as PSR 1937 +21. In this connection, an extensive study has been conducted by several authors to understand the role of anisotropy of the interior of stellar objects . The purpose of this article is to the study of minimally deformed solution for class one space-time by using gravitational decoupling method that gives a generalised solution for anisotropic compact star models. The paper is organized as follows. In Sects. 2 and 3, MGD have been used to scrutinize the interior space-time and field equation to find new embedding class one solution. The non-singularity of the solution has been discussed in Sect. 4 while, in Sect. 5, the boundary condition and determination of constraints have been analyzed. The Sect. 6 is devoted for analysis of the slow rotation approximation, moment of inertia and Kepler frequency. The elastic property of compact stars in Sect. 7 has been studied to focus and determine K e Âȃon compression modulus K e while, in Sect. 8 the energy conditions have been considered and confirmed at all points in the interior of a star. In Sect. 9, we analyzed the physical features as well as stability of the resulting solution with the help of graphical illustrations considering the equilibrium under various forces, causality, adiabatic index, and Harrison-Zeldovich-Novikov static. We also compared the solution for different values of β and analyzed for a wellbehaved solution. Finally, we summarize the results in the last section.

Einstein equations for two sources
The action of this modified matter distribution for MGD defined as [26] provided L M and L 1 are matter fields and additional Lagrangian density due to the extra source respectively. However R denotes the Ricci scalar, g is the determinant of the metric tensor g μν and β is a coupling constant. Now we define the energy-momentum tensor for both Lagrangian matter which are given bỹ By varying the action, (4) with respect to the metric tensor g μν , we get the general equations of motion where, The symbols used in above equations have their usual meanings. Let us consider the inside of spherical body is filled of an anisotropic fluid matter, therefore in the current situation the stress-energy tensorT μν takes the following form where the covariant component u ν denote the 4-velocity, fulfilling u μ u μ = −1 and u ν ∇ μ u μ = 0. Here,ρ,p r andp t represent the matter density and pressures (radial and tangential) for anisotropic matter. It is noted that the presence of this extra source Θ μν in Eq. (8) produce an anisotropies in self gravitating system that can be a scalar, vector or tensor field. We would like to mention here that the Einstein tensor is always divergence free, therefore the stress tensor T μν in Eq. (7) must satisfy the conservation law, To describe the space-time geometry inside the body for complete system we assume a spherically symmetric line element of the form, where ν and λ are functions of the radial coordinate 'r ' only. Then the Einstein field equations (7) together with the Eqs. (8), (9) and line element (11) give the following set of the equations, where the effective density (ρ), effective radial pressure ( p r ), effective tangential pressure ( p t ) and effective anisotropy corresponding to the energy momentum tensor T μν can be defined as, Further, we would like to mention that the linear combination of Eqs. (12)-(14) satisfy the conservation equation for the energy-momentum tensor T 2.2 Gravitational decoupling by MGD approach Since the system of Eqs. (12)-(14) contains seven unknown functions which are namely p(r ), ρ(r ), ν(r ), λ(r ) and three independent components of Θ. Therefore, this system has infinitely many solutions. Now we will apply the MGD approach to solve this system of equations. In this approach, the system will be converted in the such way that the field equations connected with the source Θ μν will satisfy "effective quasi-Einstein" [namely Eqs. (29)- (31)]. Now we can start by taking a solution of the Eq. (7) for the stress-tensor T μν which correspond to GR perfect fluid solution [that will be same as the Eqs. (12)- (14) when β → 0] with the line element, where the gravitational potential eλ (r ) can be defined as, here the m(r ) represents the Misner-Sharp mass function for the standard general relativity. The influence of the extra source Θ μν on the energy-momentum tensorT μν can be determined by the geometric deformation via. perfect fluid geometry {ν(r ),λ(r )} in Eq. (20) as where f (r ) and g(r ) are the deformation functions associated with the temporal and radial components of line elements, respectively. It is noted that these deformation functions depend only on radial coordinate while constant β is a free parameter. The considered MGD method allows to set g = 0 or f = 0, then for this situation the deformation will be performed only on the radial component and other temporal one unaltered (it corresponds to f = 0). By setting f = 0 we get This is called as the Minimal Geometric Deformation (MGD) along the radial component of the line element. After plugging the Eq. (24) into field equations (12)- (14), we get two sets of equations as first set corresponding to the standard Einstein field equations for an energy-momentum tensorT μν which is given as, along with the conservation equation, while second set of equations for the source Θ μν , called as quasi-Einstein equations, is given as The corresponding conservation equation ∇ ν Θ μν = 0 gives, The above expression is a linear combination of the quasi-Einstein equations. At this stage it is noted that both sources T μν and Θ μν are individually conserved, which implies that both systems interact only gravitationally.

Embedding class one condition associated with line element (20)
If the space-time (20) satisfies the Karmarkar [84] condition then the two metric functionsλ andν can be link viã The above condition implies that the four dimensional space-time (20) is embedded into five dimensional pseudo-Euclidean space i.e. embedding class one solutions. The Karmarkar condition must satisfy the Pandey and Sharma condition [85] R 2323 = 0 to describe a class one solution. On integrating (34) we get where A and B are constants of integration.
Using the definition of anisotropy in Eqs. (26) and (27) together with Eq. (35) we express the anisotropic factor,Δ(r), corresponding to energy-momentum tensorT μν by some manipulation [64] as Here some comments are in order: (i) in this paper, we want to determine a new solution for the field equations (25)- (27) using embedding class one condition. However, it is well known that we can achieve only two kinds of perfect fluid solutions, namely Schwarzschild interior solution (by vanishing of first factor in Eq. (36)) and Kohlar-Chao solution (by vanishing of second factor in Eq. (36)), for embedding class one space-time which have been already discussed in the literature. Therefore, we choose anisotropic matter distribution corresponding to energy momentum tensorT μν for determining a new solution of embedding class one spacetime. (ii) After solving the system (25)- (27), we will solve the quasi-Einstein equations (29)-(31) by taking a suitable form of the deformation function g(r ) which we are going to discuss in next section.

New embedding class one solution by MGD
Here the field equations (25)- (27) corresponding to energy momentum tensorT i j depends upon two unknown source functions, namelyλ andν. Once these source functions are determined then immediately we can obtain the thermodynamical observable likeρ,p r andp t to describe the complete structure of the proposed model. In this connection, the physical validity of this source functionν(r ) has been proposed by Lake, in which the function should be monotonic increasing function with a regular minimum at r = 0 that gives a physically viable static spherically symmetric perfect fluid solution of Einstein's equations which is regular at r = 0 [86]. On the other hand, the form of another source functioñ λ(r ) must ensure that eλ (r )) = 1 + O(r 2 ). This form ofλ gives a sufficient condition for a static perfect fluid solution to be regular at the centre. Therefore in view of above points we consider a new source functionλ(r ), to generate the physical viable solution, of the form as, where c and a are arbitrary parameters with the dimension of length −2 while n is a positive constant. Now form Eq. (37) we observe that eλ → 1 + O(r 2 ) and regular at centre as eλ (0) = 1. On inserting Eq. (37) into Eq. (35) and integrate we obtain, We observe thatν(r ) is regular at centre r = 0 and positive increasing throughout within the stellar compact object which provides a realistic compact star model. Therefore, the solution of the field equations (25)- (27) can be given by following line element, Now our next task to find all the components of Θ μ ν to describe the complete structure of the model. As we see that all these components depend on the deformation function g(r ) that requires some restrictions to lead the well-behaved solutions i.e. free of undesired mathematical and physical singularities, and non-decreasing nature, etc. These choices of deformation functions g(r ) have been widely considered by authors [87][88][89]. However, we can also be considered other options like as radial pressure associated with Θ r r to the mimic radial pressure i.e. Θ r r = p r and density associated with Θ t t to mimic energy density i.e. Θ t t = ρ, or relate only Θ-sector components through a polytropic, barotropic, or linear equation of state. It is worth mentioning that both later cases are too complicated for determining the deformation function f (r ). Therefore, we adopt the first procedure to construct the physically acceptable model then deformation function g(r ) has the following form, This form was inspired by the property of e λ that it tends to unity if r → 0 and increasing outward to describe the wellbehaved solutions. Then the explicit form of the complete space-time associated with the energy momentum tensor T μν can be written as (1 + cr 2 e nar 2 )(cr 2 + 1) where, e λ(r ) = (1 + cr 2 e nar 2 )(cr 2 + 1) The variation of both gravitational functions are shown in Fig. 1. Now the density, pressures (radial and transverse), and anisotropy for complete system are given as, 8πρ(r ) = 1 r 2 2anr 2 + 1 cr 2 e anr 2 + 1 − 2 anr 2 + 1 cr 2 e anr 2 + 1 2 − βcnr 2 cr 2 + 3 8π p r (r ) = c cr 2 + 1 −1 cr 2 e anr 2 + 1 −1 a An cr 2 e anr 2 + Bcre anr 2 ar An ce anr 2 e anr 2 cr 2 (βn − 1) − 1 +βn + Bre anr 2 cn 2ar 2 (βn + 1) + β +ce anr 2 cr 2 βn 2anr 2 a An cr 2 e anr 2 + Bcre anr 2 where, The variations of density, pressure and anisotropy are shown in Figs. 2, 3 and 4.
As we see that the central density is decreasing when n moves from 0.5 to 3.5 while the surface density increases for same n. This implies that the core will be more denser if n increases. However, it will have reverse situation for surface density. On the other hand, the pressure has totally opposite behavior than density. From Fig. 3, it can be observed that central pressure increases with increasing value of n. In this connection, we would like to mention that local anisotropies play an important role in the study of the compact objects. Since positive anisotropy leads the repulsive force which allows to more compact objects while negative anisotropy gives inward force that encourage the compact in collapsing direction. The Fig. 4 shows that the local anisotropy is increasing towards the boundary for each value of n.
The mass, compactness parameter, equation of state parameter and red-shift can be evaluated from For a realistic equation of state, the equation state parameters must be less than unity. The variations of equation state parameter and red-shift are shown in Figs. 5 and 6.

Non-singularity of the solution
The physical cogency of the solution confirm that the central values of pressure and density must be finite i.e. ρ c = 3c(βn − 1) 8π > 0, ∀ βn < 1 and c < 0, Also, it needs to ensure that any physical fluid fulfills the Zeldovich's criterion i.e. p rc /ρ c ≤ 1 which implies Now a limitations on B/A arises due to (56) and (57) as

Boundary conditions and determination of constants
The line element (2) which describes the interior of the star should join continuously with the exterior Schwar zschild metric, and can be written as with the radial coordinate r must be greater than 2m. At the pressure free interface (r = R), which needs the equality of corresponding potential functions e ν and e λ across the boundary (r = R), to get the following equations [90] 1 Whereas the extrinsic curvature or the 2nd fundamental form of stars K μν = ∇ μ r ν , where the unit radial vector r μ is normal to any surface of radius r which is also continuous at the interface (r = R). This can be stated in terms of Einstein's tensorial form as [91] [G μν r ν ] Σ = lim By using the field equations and (61) we get which implies On using the boundary conditions (60) and (63) we get +βn] cR e an R 2 e an R 2 cR 2 (βn − 1) − 1 + βn The observed values of compact stars are providing us the values of M and R, where a, n, β as free constraints.

Slow rotation approximation, moment of inertia and Kepler frequency
The moment of inertia for a uniformly rotating star with angular velocityω is assumed by [92] where, the Hartle's equation has been satisfied for the rotational dragω [93] d dr  [94] as where x = (M N R /R N R ) · km/M . The solution so obtained have been plotted mass vs I in Fig. 7 that demonstrated as n increases, the mass and moment of inertia are increasing till up to convinced value of mass and then decreases. From M-R diagram and by comparing Figs. 7 and 8, we have noticed that the mass corresponding to I max and M max are not equal. Actually, the mass corresponding to I max is lesser by ∼ 1.46% from the M max . This occurs to the EoSs due to hyperonization or phase transition to an unusual state without any strong high-density softening [95]. Using this graph we can estimate the maximum moment of inertia for a particular Due to centrifugal force, the radius at the equator increases as some factor as compare to the static one. Cheng and Harko [96] find out the approximate radius formulas for static and rotating stars as R R /R N R ≈ 1.626. Assuming the compact star is rotating in Kepler frequency Ω K = (G M N R /R 3 N R ) 1/2 and on using the Cheng-Harko formula we have plotted the M-R for rotating and non-rotating (Fig. 9). The corresponding frequency of rotating can be determined as [97,98] ν ≈ 1.22 The variation of frequency with mass is shown in Fig. 10.

Energy conditions
The energy conditions null energy condition (NEC), dominant energy condition (DEC), weak energy condition (WEC) and strong energy condition (SEC) have to be confirmed at all points in the interior of a star. Therefore, if the following inequalities hold, then the energy conditions will be satisfied simultaneously: where T μν t μ ∈ nonspace-like vector.
where i ≡ (radial r, transverse t), t μ and l μ are time-like vector and null vector respectively. With the help of graphical illustrations, the energy conditions have been checked. In Fig. 11, the L.H.S of the above inequalities have been plotted which confirms that all the energy conditions are fulfilled at the interior of stellar object.
Now the overall TOV-equation becomes −ν 2 (ρ +p r ) − dp r dr The components for different effective forces due to MGD gravitational decoupling namely gravitational force (F g ), hydrostatic force (F h ) and anisotropic force (F a ) can be defined as, The profile of three different forces are plotted in Fig. 12. From this figure we can observe that the system is in equilibrium position. Moreover, the gravitational force F g is balanced the system by joint action of anisotropic force F a and hydrostatic force F h . However, the parameter n plays an important effects on different forces as gravitational force and anisotropic force decreases in magnitude when n → 0.5 to 3.5 while hydrostatic force F h is increasing when n moves from 0.5 to 3.5.

Causality and stability condition
With the help of Causality condition, the stability situation have been analyzed. The causality condition will occurs when the sound velocities (radial (v 2 r ) and transverse (v 2 t ) are greater than zero and less than 1. The radial velocity and transverse velocity of sound can be achieved as v 2 r = dp r dρ , v 2 t = dp r dρ . Figure 13, shows the profile of radial and transverse velocities of sound, which indicates that our model fulfills the  Fig. 14 shows the stability condition that proposed by Abreu [99,100]

Adiabatic index and stability condition
The adiabatic index syndicates the basic features of the EoS on the randomness formulae and consequently contains the link between the relativistic structure of the anisotropic spheres and the EoS of the interior fluid. The stability is linked to the adiabatic index Γ , which can be written as [101], The stability of a Newtonian sphere condition is Γ r > 4/3 while, for Γ = 4/3 is the condition for a neutral equilibrium [102]. Due to the regenerative effect of pressure, this condition changes for a relativistic isotropic sphere, which is unstable. For the anisotropic fluid sphere, if the stability depends on the type of anisotropy then the situation becomes more complicated [100,101]. A recent work by Moustakidis [103] reveals that the critical value of adiabatic index strongly depends on the M/R. The critical value was found to be Figure 15, confirms that the model under consideration is stable, due to the adiabatic index is greater than 4/3.

Elastic property of compact stars
Assuming neutron stars exhibits isotropic bcc polycrystal structure one can defined the elastic properties via equation where K e and μ represents compression and shear modulus respectively. The stress tensor is given by In this work we will focus on compression modulus K e which can be determine as K e = n b (∂ p r /∂n b ) = Γ r p r [98]. The variation of compression modulus w.r.t. radius and mass are shown in Figs. 16 and 17 respectively. Harrison et al. [104] and [105] have revealed that the adiabatic index are the same of pulsating star and in slowly deformed matter, which leads to be stable if the mass of the star is increasing with central density i.e. ∂m/∂ρ c > 0 and unstable if ∂m/∂ρ c < 0.
Therefore, the mass can be furnished as a function of central density and can be defined as Its can be verified using the M − ρ c graph in Fig. 18.

Effect of coupling parameter β on the models
To complete the analysis in detail, we need to observe the effects of β on the nature of the solution. To proceed, we are needed to assume a particular value of n and then see the behavior by changing β. For n = 3 we have analyzed thoroughly and found that the physically acceptable range of β is limited with a range between 0 and 0.7. As β increases, the central density and pressure decrease while the anisotropy changes very little. However, the central values of the adiabatic index and sound speed increase with an increase of β. Therefore, the corresponding equation of state gains its stiffness along with β i.e. as the MGD+GTR coupling gets stronger we can obtain a very stiff equation of states, which may explain the current observations of very massive neutron stars (i.e. M > 2M ). Although, the range of β for a physically acceptable solution depends on the assumed values of n.

Discussion and conclusion
In this article, we have successfully incorporated the concepts of embedding class one in the gravitational decoupling formalism for the first time. This method makes a simple way of exploring new solutions in MGD, which leads to the new window of re-investigating all the existed embedding class one solutions in MGD formalism and their responses due to the additional source. The paper represents a new embedding class one solution which is deformed minimally by the gravitation decoupling technique. The MGD methodology demonstrated its adaptability in this area, making it an important and acceptable solution for EFE. It is reproduced through the graphical analysis, where the variation of the metric functions with the radial coordinate r , (see Fig. 1) for M = 1.74M and R = 9.528 km considering a = 0.001 km −2 and β = 0.01, and the deformation function g(r ) as in Eq. (24). The deformation function will be vanished at r = 0, while g(r ) = n as r approaches to infinity. Thus, the metric potential functions are well-behaved at the center and finite and regular throughout of stars. Therefore they are proper to produce new models for anisotropic compact stars. Figures 2 and 3, shows the behavior of density, pressures(radial and transverse) with respect to M = 1.74 M and R = 9.528 km considering a = 0.001 km −2 and β = 0.01, which detected that the model is non-singular, furthermore the model is positively finite, and monotone decreasing functions throughout the interior of the star, and achieve their maximum value at the center. Also, radial pressure vanishes at the boundary. The anisotropy factor, which is given in (Fig. 4) with radial coordinate r . However, Δ(0) = 0 at the center and it is positively increasing away from the center. From Fig. 5, it is clear that the EoS is characterized by the parameters ω r and ω t relating to radial coordinate r , in which, the EoS factors of the model are less than unity i.e. within the region and demonstrated as a well-behaved model. Figure 6 shows the variations of surface red-shift with radial coordinate r . Thus, the surface of redshift z(r ) → 0 as r → 0 and subsequently monotonically increasing onto the boundary. For n = 2 yields larger moment of inertia M max and K e (Figs. 7, 8). Also, we have noticed that from Figs. 7 and 8 the M − I graph is more sensitive and/or sharp in the stiffness of equation of states than M-R graph. By We are using the same color notation in all the graphs from Figs. 19, 20 and 21 as β = 0 (black), 0.014 (brown), 0.028 (blue), 0.042 (green), 0.056 (red) and 0.07 (cyan) using the concepts of Ghosh and Cheng-Harko i.e. Eq. (70) and R R /R N R ≡ 1.626 one can compare the M-R graphs of rotating and non-rotating limits in one frame (Fig. 9), while the variation of frequency with mass for different values of n is shown in Fig. 10. The EC with radial coordinate r for 4U 1608-52 (M = 1.74M , R = 9.528 km) by taking a = 0.001 km 2 and β = 0.01 are shown in Fig. 11, which confirms that all EC are fulfilled at the interior of the stellar object. While Fig. 12, shows the profile of three different forces to observe that the system is equilibrium i.e variation of forces in TOV-equation with radial coordinate for 4U 1608-52 (M = 1.74 M , R = 9.528 km) by taking a = 0.001 km −2 and β = 0.01. In the Table 1, we have presented the values of physical parameters for different values of n. The profile of radial and transverse velocities of sound has been motivated in Fig. 13, which indicates that our model fulfills the causality condition. While (Fig. 14) shows the stability condition proposed by Abreu [99].
The satisfactions of static stability criterion, modified TOV-equation and Herrera's cracking method also ensures that the solution is static, equilibrium and stable. Also, we noticed that the system is stable due to the adiabatic index Γ r is greater than 4/3 as shown in (Fig. 15), and is also increasing monotonically outward. The K e − M graphs (Fig. 16) signifies that as the mass of compact star increases, the compression modulus decreases while the K e with radial coordinate r graph implies an increasing trend of K e as r approaches the center i.e .the parameter n increases the stiffness of the corresponding equation of states increases. While in Fig. 17 a variation of K e with radial coordinate r for 4U 1608-52 (M = 1.74M , R = 9.528 km) by taking a = 0.001 km −2 Fig. 20 Variation of anisotropy, TOV-equation (F a −dashed, F g −small dashed, F h −solid), sound speed (v r −solid, v t −dashed) and stability factor with radial coordinate r for 4U 1608-52 (M = 1.74M , R = 9.528 km) by taking a = 0.001 km −2 and n = 3. We are using the same color notation in all the graphs from Figs. 19, 20 and 21 as β =0 (black), 0.014 (brown), 0.028 (blue), 0.042 (green), 0.056 (red) and 0.07 (cyan) Fig. 21 Variation of Γ r , energy condition (EC) ( p r + ρ−solid, p t + ρ−small dashed, p r + 2 p t + ρ−dashed), redshift and compression modulus with radial coordinate r for 4U 1608-52 (M = 1.74M , R = 9.528 km) by taking a = 0.001 km −2 and n = 3. We are using the same color notation in all the graphs from Figs. 19, 20 and 21 as β = 0 (black), 0.014 (brown), 0.028 (blue), 0.042 (green), 0.056 (red) and 0.07 (cyan) and β = 0.01 are motivated for different n. Also, the central values of these physical parameters can be confirmed by the Fig. 18 for n = 2 than n = 0.5. The mass function and compactness effects are able to illustrate the mass growth to choose an appropriate specific mass M. So, the profiles of understanding of the compactness and mass function are shown in Fig. 18. The m(r ) → 0 as r → 0 and monotonically increasing toward the boundaries. The former tendency can be explained as the mass increases the central density also increases which may leads to generation of many interesting particles that unstiffen the equation of state and the compression modulus. The later one is due to the central density is highly dense as compare to the surface that leads to more compression modulus at the core than its surface. As the maximum mass M max increases when n increases, the spinning rate ν max also increases. We have also compared the nature and behavior of the solution by assuming a particular value of n with different values of β. For this case, we have found that the central values of density and pressure decreases with increase in β (Fig. 19). However, other physical behaviors of solution for different β are given in Figs. 20 and 21. On the other hand, the stiffness increases with increase in β as the adiabatic index increases and the speed of sound approaches the speed of light. Although, the anisotropy changes in a very small amount when changing the coupling constant β. The acceptable range of β depends upon the chosen value of n. For n = 3 the possible range lies in 0 ≤ β ≤ 0.7. If β > 0.7, the trend of the density start increase slightly and decreasing near the surface, and the solution start violating causality condition.
Summing up, we can conclude that our models are physically acceptable to describe minimally deformed class one space-time by gravitational decoupling based on the results so obtained.