A completely deformed anisotropic class one solution for charged compact star: a gravitational decoupling approach

In this article, we have investigated a new completely deformed embedding class one solution for the compact star in the framework of charged anisotropic matter distribution. For determining of this new solution, we deformed both gravitational potentials as ν↦ξ+αh(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu ~\mapsto ~\xi +\alpha \, h(r)$$\end{document} and e-λ↦e-μ+αf(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 } \mapsto ~e^{-{\mu }} + \alpha \,f(r)$$\end{document} by using Ovalle (Phys Lett B 788:213, 2019) approach. The gravitational deformation divides the original coupled system into two individual systems which are called the Einstein’s system and Maxwell-system (known as quasi-Einstein system), respectively. The Einstein’s system is solved by using embedding class one condition in the context of anisotropic matter distribution while the solution of Maxwell-system is determined by solving of corresponding conservation equation via assuming a well-defined ansatz for deformation function h(r). In this way, we obtain the expression for the electric field and another deformation function f(r). Moreover, we also discussed the physical validity of the solution for the coupled system by performing several physical tests. This investigation shows that the gravitational decoupling approach is a powerful methodology to generate a well-behaved solution for the compact object.


Introduction
It is well known that the astrophysical compact objects are formed due to the gravitational collapse. These compact stars belong to three different groups like neutron stars, white dwarfs, and black holes. The arrangement of these stars is constructed according to their internal structure. The black hole is the densest stars in the universe. Till now, it is an open problem in astrophysics to know about the behavior, nature and exact configuration of these astrophysical compact objects which are more compact than the normal compact objects, therefore in recent days, it is still an active field of research. Based on the hypothesis of strange matter, the a e-mail: sunil@unizwa.edu.om strange quark matter could be more stable than nuclear matter. Therefore, the neutron stars must essentially be composed of pure quark matter. The current expansion in a cosmological investigation has explained gradually the origin and distribution of substance and growth of compact objects in the Universe. Few of their assets, such as rotation frequencies, masses, and emission of radiation are computable. The measurements of the essential parameters decide the nature of compact objects which is still an observational challenge because these properties are not directly associated with the observations, such as the internal configuration of masses and radii that need the growth of hypothetical models. In view of theoretical computation, the mass and radii are calculated by solving the hydrostatic equilibrium equation that gives an equilibrium condition between the hydrostatic and gravitational force. In the case of Einstein general relativity, the Tolman-Oppenheimer-Volkoff (TOV) equations can represent the equilibrium condition for a spherical compact object. However, equation of state is necessary for describing the complete structure of these spherical compact objects. In the modern days, the experimental data shows the existence of such astrophysical objects which are observed at very high densities [1], and few of the astrophysical objects like Her X-1 (X-ray pulsar), PSR 0943+10, 4U 1820-30 (X-ray burster), 4U 1728-34, RX J185635-3754 (X-ray sources), and SAX J 1808.4-3658 (millisecond pulsar) extremely favor for the probability that these compact objects could essentially be strange stars. In further arguments, there is no any resilient confirmation to establish the mechanism for the understanding of astrophysical compact objects. According to the theoretical investigations, such compact objects are composed of a perfect fluid [2]. Mostly, the bag equation of state (EOS) and polytropic EOS ( p = k ρ γ ) have been extensively applied to describe less compact star and a white dwarf [3][4][5]. In this connection Herrera and Barreto [6] conveyed a general study on polytropic relativistic stars with anisotropic pressure while polytropic quark star models have been studied by Lai and Xu [7]. From a theoretical point of view, the impact of anisotropic compact objects was first introduced by Bowers and Liang [8]) and another study headed by Ruderman [9] revealed that nuclear matter may be composed of anisotropic pressure at a very high density which is the order of 10 15 g/cm 3 . In this scenario the several objects have been discovered by the authors (see Refs. [10][11][12][13][14][15][16][17][18][19][20][21]). It has also been investigated that the pressure anisotropy can disturb the critical mass, compactness, and steadiness of the extremely dense relativistic compact objects. Recently Chaichian et al. [22], Ferrer et al. [23] have proposed that the pressure anisotropy can be produced by the action of the magnetic field on a Fermi gas. Recently, Maurya et al. [24][25][26][27][28][29] have discovered several charged and uncharged solution for the compact stars. The effect and role of anisotropy with and without equation of state (EOS) were proposed by Harko and Mak [30], Maurya et al. [31][32][33][34][35], Varela et al. [36]. On the other hand, Sharma and Maharaj [37] have obtained an analytical extended solution by taking a specific choice of mass function. In the current days, the embedding class one solution has much attraction towards the researchers. Many mathematicians and physicists have explained the concept of embedding that the curved spacetime can be embedded into higher dimensional flat spacetime. The Riemann proposed the innovative idea of Riemannian geometry by giving an exhaustive study of the geometric objects in 1850. Sudden after, some mathematicians proposed some concept that the Riemannian manifold can be embedded into Euclidian space of higher dimensional. This concept was known as the isometric embedding problem which was unproved for many years. The first time in 1871, the Schlaefli guessed the prospect of embedding the Riemannian manifold locally and isometrically into the higher-dimensional Euclidean space under some constraint. Later on this conjecture was proved and founded by Janet [38], Cartan [39] and Burstin [40]. They have proved that any l dimensional Riemannian space can always be locally embedded into any pseudo-Euclidean space of dimension L ≥ l(l +1)/2. In the case of symmetric, we need at least L = (l + 1)-dimension of pseudo-Euclidian flat space for the embedding of l dimensional Riemannian geometry isometrically. This embedding theory has been also well explained by Stephani [41] and his collaborators in their pioneering work. In this concerns, the most popular work on embedding class one have been done by several authors [42][43][44][45][46][47][48][49][50][51][52][53][54][55].
In recent days, the technique of gravitational decoupling by minimal geometric deformation (MGD) has a great interest among the researchers. This is a powerful tool to extend any previously available solutions into a more complex form or provide more general new solutions. This effective minimal geometric deformation (MGD) mechanism and its extended case have been developed by Ovalle and his collaborators [56][57][58]. The MGD solutions in the context of Randall-Sundrum gravity have been discovered by Ovalle et al. [59][60][61]. Later on, Ovalle and Casadio et al. [62] have determined a new black hole solution by using this minimal geometric deformation (MGD) approach as well as they also explained that how the MGD introduce the anisotropy in the self-gravitating system and how spherically symmetric matter distribution changes the Schwarzschild vacuum solution. Using this MGD approach, several authors have obtained the solutions in a more complex forms, see Refs. [63][64][65][66][67][68][69][70][71][72][73][74][75][76]. Very recently, Singh and Maurya [77] have determined the first solution of embedding class one spacetime using the minimal geometric deformation (MGD) approach for anisotropic matter distribution. Now we are going to discuss how the MGD works in solving the Einstein field equation (EFE) which are the following: (I) First we choose the energy-momentum tensor T μν corresponding to already known solution. Then we can extend this known solutions of Einstein field equation into more complex form by joining of extra new stress tensor which will be introduced by MGD. Then the obtained new solution can be described by coupling of the T μν for the seed solution and new stress tensor through a non-dimensional coupling constant α as: similarly we can extend the new energy momentum tensor T (1) μν as, (2) and continue the same procedure as n times. In this process we can extend our simple initial solution of Einstein field equation connected with the source T μν into more generalised form which is related to the source T μν =T (n) μν . It is important to mention that each different component for the source T n μν is individually conserved, (II) On the other hand, this methodology can also work very well in reverse order to obtain the solution for compact objects. In order to proceed with this reverse methodology, first, we separate a more complex energy-momentum tensor T μν into many different simpler energy tensor components like T n μν . After these separations, we will solve the Einstein field equation corresponding to each component of T n μν for the compact object. In this process we will get multiple solutions that are connected with the above source T n μν . Finally the complete solution of the EFE for the original stress tensor T n μν can be obtained by joining of above each distinct solution. Moreover, we can describe this methodology by taking a particular example which can be given as: let us consider two gravitational sources namely S 1 and S 2 where S 1 is a source connected standard Einstein's equations while second source S 2 corresponding to the quasi-Einstein equations. After specifying both sources, we solve both systems of equations individually. Finally, we join both solutions and find a solution for the complete system which is given as S 1 ∪ S 2 . Moreover, it is well known that Einstein's field equations are highly nonlinear in nature. Therefore, the minimal geometric deformation (MGD) approach is a powerful methodology to solve the EFE.
The article is organised as follows: In Sect. 2, we defined the field equations for the charged matter distribution via gravitational decoupling (GD) and corresponding interior space-time. For this purpose, we divide this Sect. 2 into three Sects. 2.1, 2.2 and 2.3. In Sects. 2.1 and 2.2, we defined the action for two sources and corresponding field equations for both sources, respectively. The embedding class one condition is given in Sect. 2.3. We have obtained new solution for charged anisotropic matter distribution via GD in Sect. 3. In Sect. 4, we have performed the matching condition to determined the arbitrary constants and necessary parameters. The physical analysis has been presented in Sect. 4 which consist three subsections namely Sect. 5.1. Regularity of gravitational decoupling model, Sect. 5.2. Equilibrium condition for GD model, and Sect. 5.3. Stability of the stellar model via cracking. The mass-radius ratio and surface redshift are given in Sect. 6 where we discussed the gravitational and effective mass. The final section has been made for discussion and conclusion.
2 The field equations for the charged matter distribution via gravitational decoupling (GD) and corresponding interior space-time

The action for two sources
We define the action for the modified matter distribution by using the extra source as [58] where L M and L E represent corresponding the matter fields and additional Lagrangian matter density due to presence of the extra source. It is well defined notation as R represents a Ricci scalar while g is the determinant of the metric tensor g μν . Then the energy-momentum tensor corresponding to both Lagrangian matter distributions can be written as, After varying the action (4) with respect to the metric tensor g μν , we obtain the general equations of motion for the extra stress tensor Θ μν as Let us consider the spacetime being static and spherically symmetric that describes the interior of the compact object. This spacetime can be written in the following form as, Now using Eqs. (7) and (8), we obtain Einstein's equations for static spherically symmetric spacetime for the energy tensor T eff μν as , Then the linear combination of Eqs. (9)-(11) satisfy the following conservation equation Our next aim to write the field equations corresponding to both stress tensor and solve them separately using the minimal geometric deformation approach.

The field equations for T μν and Θ μν via MGD approach
In this approach, we will look that the system will be converted in such a way that the equation of motions associated with the stress tensor Θ μν satisfy the effective "quasi-Einstein system". However, the impact of the extra source Θ μν on the energy tensor T μν can be determined by the geometric deformations as, where f (r ) and h(r ) denote the geometric deformation functions undergone by the radial and temporal metric component, respectively while the constant α is a free parameter. These deformations are known as the complete geometric deformation (GD) or extended geometric deformation along the radial and temporal components of the line element. Now let us substitute the deformed metric function from Eqs. (13) and (14) in the field Eqs. (9)-(11) which separate the system into two sets namely: (a) first set as the standard Einstein field equations corresponding to energy-momentum tensor T μν which is described by line element (8), and (b) The second set gives the equation of motions for the source Θ μν , known as "quasi-Einstein equations", which are given as, where F 1 and F 2 can be defined as, Along with the conservation equations for both systems, Now by motivation of Ovalle [58], we define Maxwell energy-momentum tensor for a charged self-gravitating distribution which corresponds to the extra source Θ μν in Eq. (7). Then electromagnetic field tensor Θ μν can be written as, It is well known that the anti-symmetric electromagnetic field tensor F μν in Eq. (23) satisfies the Maxwell field equations, where the electromagnetic four current vector j i is defined as with charge density σ = e ν/2 J 0 (r ). It is noted that only the non-vanishing component of the four-current is J 4 for static matter distribution. Since in the case of spherical symmetry, the four-current component is only a function along radial direction 'r'. Then corresponding non-zero components for electromagnetic field tensor are F 01 and F 10 that is related as F 01 = −F 10 . These components describe the electric field along the radial vector r . From Eqs. (24)-(26), the expression of electric field component for spherically symmetric line element (8) can be written in the form of electric charge as, where q(r ) describes the electric charge which is contained within the sphere of radius r. Then q(r ) can be defined by using relativistic Gauss law and corresponding electric field E as, In order to define T μν , we consider the fluid matter inside of spherical body is completely filled by anisotropic fluid. In this case the energy momentum tensor T μν can be given as, where the covariant component u ν denote the 4-velocity, fulfilling u μ u μ = −1 and ν μ u μ = 0. Here, ρ, p r and p t represent the matter density and pressures (radial and tangential) for anisotropic matter distribution. Then the components for T μ ν and Θ μ ν can be expressed as, By plugging the Eqs. (31) and (32) in the field equations (15)- (17) and (18)-(20) we obtain, and Then the Eqs. (21) and (22) yield the following conservation equations, where Ψ = q 2 r 4 .

Embedding class one condition associated with the geometry {T μν , ξ , μ}
If the space-time for T μν satisfies the Karmarkar [78] condition, then where the Karmarkar condition must satisfy the Pandey and Sharma condition [79] R 2323 = 0 to describe a class one solution. Then Eq. (41) provides an differential equation which relates both metric functions ξ and μ as, If the space-time satisfy the above condition (41) then corresponding four dimensional space-time can be embedded into five dimensional pseudo-Euclidean space which describe an embedding class one space-time.
On integration of the Eq. (42) which yields where A and B are constants of integration.
Using the definition of ansitropy in Eqs. (34) and (35) together with Eq. (43) we express the anisotropic factor, Δ(r ) corresponding to energy-momentum tensor T μν by some manipulation [43] as From Eq. (44) we see that if Δ = 0 then we can recover only two kinds of isotropic fluid solutions, namely Schwarzschild interior solution which can be obtained by taking the first factor in Eq. (44) to be zero while vanishing of the second factor in Eq. (44) gives a Kohlar-Chao solution. Both solutions are already discussed in the literature. Therefore, we choose energy-momentum tensor T μν for the anisotropic matter distribution to find a new solution for embedding class one space-time. Before proceeding to the next section we will mention some important comments here: (a) Initially we will solve the field Eqs. (33)-(35) via embedding class one condition to obtain a new solution for the anisotropic star.
(b) Then we will solve of Maxwell Eqs. (36)-(38) (known as quasi-Einstein equations) together with conservation Eq. (40). After solving both systems we will determine the complete thermodynamic observable for the charged anisotropic model.

New embedding class one solution by gravitational decoupling
In this section, we will present the solutions for both systems. Firstly, we solve the field Eqs. (33)-(35) corresponding to energy-momentum tensor T i j which involves two unknown source functions, namely μ and ξ . Once both metric functions are specified then immediately we can obtain the complete thermodynamical observable ρ, p r and p t for stress tensor T i j . It is well known that according to the Lake algorithms, the source function ξ(r ) should be monotonic increasing and regular minimum at r = 0 which provides a physically viable static spherically symmetric perfect fluid solution of Einstein's equations that is regular at r = 0 [80]. Therefore, we chose the source function e ξ same as Durgapal-IV [81] solution, where a and C are positive constants and the dimension of a is length −2 . The metric function ξ(r ) in Eq. (45) is regular at centre r = 0 and positive increasing throughout within the stellar compact object that gives a realistic compact star model. Now using the embedding class one condition by plugging the Eq. (45) in Eq. (43), we get a following expression for e μ of the form, where D = 16 a C/B. From Eq. (46) we note that the metric function e μ can be expressed in the form of e μ → 1+ O(r 2 ). This form provides a necessary condition to be a regular model at the origin. Therefore, the solution of the Einstein field Eqs. (33)-(35) can be given by following line element, Now we will determine electric field E and deformation functions f (r ) and g(r ) by solving of Maxwell Eqs.  42) is a first order ordinary differential equation which can be written in form of integral as, Substituting the value of ξ and μ from Eqs. (45) and (46) On inserting the value of h (r ) in Eq. (49) and the integrate w.r.t r we get, Now substituting the Eq. (51) in Eq. (36) we obtain, On integrating of Eqs. (50) and (52) which yield the deformation functions h(r ) and f (r ) as, where, h 1 (r ) = 3a(1 + 2D)r 2 + 2a 2 D(6 + D)r 4 , h 2 (r ) = (210 + 336ar 2 + 280a 2 r 4 + 120a 3 r 6 + 21a 4 r 8 ).
Then the complete deformed metric potentials are given as, and corresponding space-time can be expressed as, where h(r ) and f (r ) are given by Eqs. (53) and (54). This deformed space-time (57) will provide a complete structure of the matter distribution for effective energy momentum tensor T eff μν . Now the effective pressures and effective density can be written as, Then the corresponding expressions for ρ eff , p eff r and p eff t can be given as, a 2 r 2 16+D(4+27ar 2 +36a 2 r 4 +14a 3 r 6 ) . Using Eqs. (62) and (63), we obtain effective anisotropy Δ eff as a 2 r 2 16+D(4+27ar 2 +36a 2 r 4 +14a 3 r 6 ) .

Matching conditions
For the study of stellar interiors, it is a necessary and important point that the interior spacetime geometry must join with the exterior spacetime geometry in a smooth way. For doing this, we will introduce the well known Israel-Darmois junction conditions [82,83]. These junction conditions require the continuity of the gravitational potentials across the surface Σ boundary (defined by r = R) or equivalently The above expressions are called as the first fundamental form. On the other hand, It is also required for a vanishing radial pressure at the boundary of the star as, The above condition is known as the second fundamental form which determines the size of the compact object i.e the radius R. Moreover, since the compact object contains a charge, then the continuity of the electric charge q(r ) on the surface of the compact object must be added to the first fundamental form, where Q is the total charge of the compact object. On the other hand, if the electric charge is a continuous function of the radial coordinates across the boundary then the electric field E(r ) will also continuous across the boundary of the compact object. Then, it is observed that the Riemannian exterior manifold can be described by the well known Reissner-Nordström spacetime, which is given as, we obtained the explicit expressions by matching of the first fundamental form as, Here the parameter M represents the total mass containing within the compact star. By solving of the Eqs. (71) and (72) and second fundamental form given by (68) we get the expressions for the constant parameters D, M and C as, where, +3(48 + 9736α) a 4 R 8 +153α 2 a 4 R 8 + 24α(928 + 3α)a 5 R 10 , D 4 (R) = 4(1524 − 7α)αa 6 R 12 − 16α 2 a 7 R 14 + 4α 2 a 8 R 16 .

Regularity of the gravitational decoupling model
For the physical acceptable model, the solution must be free from the geometrical and physical singularity i.e. the deformed gravitational potentials must be regular at the centre and throughout within the compact object. From Eqs. (55) and (56) we can easily observe that e ν(r ) and e λ(r ) fulfil the requirement of regularity at centre as, Since the gravitational potentials must be positive within the compact object then the Eq. (76) provided that the constant C must be positive. The profile of gravitational potentials has been shown in Fig. 2 which is shows that the gravitational potentials are regular and positive throughout the fluid model. Moreover we have also determined the behavior of effective physical parameters such as density (ρ eff ), radial pressures ( p eff r ) and tangential pressure ( p eff t ) within the selfgravitating stellar object. For the regular stellar model, the above physical parameters must be positive and decreasing outward as well as the maximum value must attains at the centre of the compact object (Fig. 3). For this purpose we determine the central value of these physical parameters to fix the range of the arbitrary constants which are as follows, Now, the ratio of pressure and density should be less than unity that gives, In addition to above the causality condition must satisfy everywhere within the star i.e. v 2 r < 1 and v 2 t < 1. Then we obtain the values of (v 2 r ) r =0 and (v 2 t ) r =0 as, Since D is positive then from the Eqs. (80) and (81) we get, Using the Eqs. (77)- (79) and (82), (83), we get the following inequality as, It can be easily verified that If the coupling parameter α is very small as near to zero then the inequality (84) holds for D while the Eq. (85) will provide the restriction on D for the large values of α.

Equilibrium condition for GD model
In this section, we will discuss the equilibrium condition for the geometric deformation (GD) for decoupling source Θ μν via generalised Tolman-Oppenheimer-Volkoff (TOV) equation. It is well known that the equilibrium condition can be described from the conservation equation of T eff μν as follows, The explicit form of the Eq. (86) can be described as, After substituting the components of T μ ν and Θ μ ν from Eqs. (31) and (32) into Eq. (87) we get, Now the Eq. (88) shows that the system is balanced under different forces as, Gravitational force: Electric force: Anisotropic force: Coupling force: These above forces form the equilibrium of the system. The variation of these forces is intensely connected to the constant parameter α. Figure 4 shows that how we reached the equilibrium of the GD system under these different forces. Now we will focus on several points based on the different As we see that the electric force (F e ) and extra coupling force (F α ) have no effect to balance this mechanism however the gravitational force (F g ) is balanced by joint action of hydrostatic (F h ) and anisotropic (F a ) forces. On the other hand, the Figures top right, bottom left and bottom right have been plotted for α = 0.10, α = 0.25 and α = 0.50, respectively. From these figures, it is interesting to note that the electric force F e and coupling force F α have same values in magnitude within the star but they are acting in the opposite directions while F g can be balanced by combining of F a and F h . Therefore, the equilibrium is achieved for each different value of α.

Stability of the stellar model via cracking
Normally after analyzing of equilibrium condition of the system, it is necessary to determine the stability of the model. For this purpose, it is essential to examine whether the perturbation introduced by the pressure anisotropy yields a stable equilibrium system or not. To analyze this condition, we use Abreu's criteria [84] which has been developed by Herrera's cracking criterion [85]. The Abreu's criteria for the stability analysis of the stellar compact objects can be described by the subliminal radial and tangential sound speeds as, Clearly it can be observed that both radial and tangential subliminal speed of sound must fulfill the causality condition i.e. v 2 r < 1 and v 2 t < 1 everywhere inside the compact object (by taking the speed of light c = 1 in relativistic geometrized units). From the Fig. 5, it is clear that the obtained stellar model satisfies the causality condition. But we would mention an interesting point that for increasing value of α the radial velocity behaves dual nature i.e. first it is decreasing and then start increasing when α ≥ 0.275 while tangential velocity decreases throughout for all values of α. Now let us discuss the stability of the model through the cracking concept. Form the Fig. 6, we observe that stability factor v 2 t − v 2 r lies between −1 and 0 and the factor v 2 r − v 2 t lies between 0 and 1. This shows that radial velocity is greater than the tangential velocity throughout the star as well as there is no cracking appears within the model. Therefore, the obtained self-gravitating model is stable.
6 Mass-radius ratio and surface redshift

Gravitational and effective mass
To describe the structure of the compact object it is necessary to study the upper limit of the mass-radius ratio. It is well known that the upper limit of the mass-radius ratio for the isotropic matter distribution is given by Buchdahl's [86]. This limit was founded in the context of perfect fluid which has a decreasing energy density outwards. The upper bound of mass-radius can be given as, where M is the total mass of the perfect fluid sphere and R the radius of the star in Schwarzschild coordinates, which is determined by the vanishing of pressure on the boundary of  Fig. 4 the star. The above inequality has a significant implication, which gives the maximum limit of the surface redshift z s . However, the Buchdahl's limit is changed in the presence of electric charge. In this situation, the upper and lower limits are given by Andreasson [87] and Bohmer and Harko [88] as, Here we would like to mention that the mass M appears in the Eq. (96) is not the same as the total mass in (95). More clearly, we can say that the total mass in (95) describes the mass for an isotropic fluid sphere which is exactly equal to the gravitational and effective mass in case of the perfect fluid matter distribution. Explicitly it can be written as, On the other hand, the mass M appears in expression (96) defines the gravitational mass for the charged compact object which is given as From the Eqs. (97) and (98), it can be observed that the gravitational mass has more value in the case of electrically charged compact structures as compared to uncharged perfect fluid matter distribution. Therefore, the effective mass for the charged compact object can be given as,

Surface and gravitational redshift
It was discussed in the above section that the mass-radius ratio plays an important role to decide the value of surface redshift Z s . Then the surface red-shift can be calculated as, while the gravitational redshift within the stellar models can be calculated by following formula, z = e −ν(r )/2 − 1.
The variation of gravitational redshift within the stellar model is given in Fig. 7. From this figure, we predict that the gravitational redshift is maximum at centre and decreasing away from centre. The numerical values of surface redshift is given in Table 2 for different values of coupling parameter α.

Discussion and conclusions
Our study is dedicated to a new class of completely deformed solution for the spherically symmetric charged anisotropic compact stellar model. This completely deformed solution has been obtained by using the gravitational decoupling approach in the framework of embedding class one spacetime. In our investigations, we confirm that the 'gravitational decoupling' technique is a powerful mechanism to determine a more general solution of the Einstein field equations in the scenario of embedding class one spacetime. To be more specific, in this study first we define the action for the modified matter distribution via introducing another extra source. Then we write the general equation of motion by  Fig. 7 The behavior of gravitational redshift (Z ) verses radial coordinate r/R. We have used the same data values for plotting of this figure as given in Fig. 1 varying this action along with the metric tensor g μν . Next, we obtain the Einstein field equations of spherically symmetric metric corresponding to the above general equation of motion for effective energy-momentum tensor T eff μν , which is the combination of two sources namely T μν and Θ μν . By using gravitational decoupling ξ → ν = ξ + α h(r ) and e −μ → e −λ = e −μ + α f (r ), we separated the effective system into two separate systems as first corresponds to T μν (known as Einstein field equations) and second for Θ μν (Known as quasi-Einstein field equations). Now our aim is to solve both systems individually. As we see that the first system [Eqs. (33)- (35)] depend on two unknowns ξ and μ. To solve this system we chose a particular form of gravitational potential, ξ = C(1 + ar 2 ) 4 , same as Durgapal IV solution which satisfies all the physical and mathematical requirements for a regular and well-behaved solution [89]. Then we apply the embedding class one condition to determine the solution for this system for anisotropic matter distribution. Now we will focus on second system which is given by the Eqs. (18)- (20). Before solving these equations, first we converted the quasi-Einstein Eqs. (18)-(20) into a Maxwell field Eqs. (36)- (38) by applying the Ovalle [58] approach. Now we see that field Eqs. (36)-(38) depend on mainly three unknowns as electric charge q, deformation functions h(r ) and f (r ). For determining the solution of this Maxwell system, we need to solve the conservation Eq. (40) which is a first-order ordinary differential (ODE) equation containing two unknown functions Ψ and h(r ), where (Ψ = q 2 /r 4 ). Therefore, we solve this ODE by taking a specific form of deformation function h(r ) [which is non-singular at the center and increasing away from the centre, see Eqs. (50) and (53)], and then obtained a closed-form expression for electric field (Eq. 51). Then after we obtain an another deformation function f (r ) [see Eq. (54)] by solving of the Eq. (36) together with Eq. (51). In this way, we have solved both systems completely and the obtained solution for this com- plete system is given by line element (57). Moreover, we also determined the effective quantities density, pressures (radial and tangential) and anisotropy factor to describe the structure of the charged compact object. To perform the graphical analysis we used the matching conditions and obtain all the arbitrary constants. The physical analysis of the solution for the charged compact object is given as follows: The behavior of electric charge q is shown in Fig. 1 for different values of α with respect to specific mass-radius ratio M/R = 0.2. From this figure we see that electric charge is zero at centre and increasing towards the boundary of star. We have also calculated the amount of charge on the boundary in Coulomb unit as: (i) Q = 0.5593 × 10 20 C for α = 0.10, (ii) Q = 0.8748×10 20 C for α = 0.25, (iii) Q = 1.0278×10 20 C for α = 0.35, (iv) Q = 1.2156 × 10 20 C for α = 0.50. This amount of charge shows that it is increasing as α increases.
In Fig. 2, the variation of deformations function h(r ) and f (r ) are given by top (left and right) figures while the bottom (left and right) figures are corresponding to gravitational potentials, viz. e λ and e ν for different α and fixed M/R = 0.2. It can be clearly seen that for each different values of α, the both deformation functions and gravitational potential are free from singularity and increasing throughout within the compact stellar model. The such behavior of h(r ), f (r ), e λ and e ν gives a physically well-behaved stellar model. The variation of effective matter distribution p eff r and p eff t against the radial coordinate r/R are given shown in top (left and right) of Fig. 3 while ρ eff and Δ eff are featured in bottom (left and right) of this Fig. 3. From this figure, we observe that effective pressures and density are maximum at the centre and decreasing outward. Also, the effec-tive anisotropy is zero at the centre (i.e. radial and tangential pressure are equal at center) and positive throughout within in stellar model. This feature of anisotropy shows that the force due to anisotropy is directed outward which allows the construction of a more compact object. Moreover, the numerical values for central pressure, central and surface density are given in Table 1. The surface density of the present model lies between 3.131ρ s − 3.112ρ s which is higher than the normal nuclear density ρ s [9,90,91]. The above range of surface density confirm that the obtained stellar structure is an appropriate candidate for the hypothetical ultradense compact object.
To check the equilibrium position of all forces, we have studied the generalized TOV equation for a gravitationally decoupled system. In Fig. 4, we have shown the distributions of all forces F h , F a , F g , F e and F α for the coupling parameter α = 0.0 (top left), α = 0.10 (top right), α = 0.25 (bottom left) and α = 0.50 (bottom right). From these individual figures, we see that all forces achieve the equilibrium condition, which shows the stability of the system. In this gravitational decoupled system, a new kind of force F α is coming due to the geometric deformation in the gravitational potential source. For all values of α > 0, this extra force F α acts along inward direction and behaves an attractive nature which is balanced by electric force F e such that F e + F α = 0, this yields an equilibrium condition for the extra source Θ μν . On the other hand, the gravitational force F g is balanced by the joint action of hydrostatic force F h and anisotropic force F a such that F g + F h + F a = 0, which provides the stable equilibrium for the matter distribution T μν . By combining of both equilibrium conditions, we get F g + F h + F a + F e + F α = 0.