Extended gravitational decoupling (GD) solution for charged compact star model

In the present article, we have obtained a new solution for the charged compact star model through the gravitational decoupling (GD) by using a complete geometric deformation (CGD) approach (Ovalle, Phys Lett B 788:213, 2019). In this approach, the initial decoupled system is separated into two subsystems namely Einstein–Maxwell’s system and quasi-Einstein system. We solve Einstein–Maxwell’s system by taking well known Tolman–Kuchowicz spacetime geometry in the context of the perfect fluid matter distribution. On the other hand, the second system introduce the anisotropy inside the matter distribution which is solved by taking an EOS in θ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta $$\end{document} components. The boundary conditions have been derived to determine the constants parameter. To support the mathematical and physical analysis of the present GD solution, we have plotted all the graphs for the compact objects PSR J1614-2230, 4U1608-52 and Cen X-3 corresponding to the constant α=0.001\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =0.001$$\end{document}, 0.0012 and 0.0014, respectively. Moreover, we also studied the equilibrium and stability of the solution. The present study shows that the GD technique is a very significant tool to generalize the solution in a more complex form or one matter distribution to another matter distribution.


Introduction
It is a great challenge to find the exact solutions of Einstein's field equation (EFE) due to its non-linearity. Last 100 years, many exact solutions for EFE have been obtained in the context of perfect fluid matter distribution [1]. But only few of them are well behaved that describe astrophysical compact objects. Therefore, still, the researchers are busy to find the exact solution of EFE in different contexts. Usually, there are two methodologies that have great interests among all the researchers for solving the field equations. Fist procedure is that consider one metric function or decreasing energy a e-mail: sunil@unizwa.edu.om (corresponding author) density in order to solve the system. This approach leads an equation with other unknowns which is determined by using the pressure isotropy condition. But this process does not provide always an admissible exact solution. In a second way, we can solve the EFE by specifying an equation of state (EOS) that relates the pressure and density. In this way, we also did not get an exact solution easily for perfect fluid matter distribution due to its complicated integral from. On the other hand, it is not necessary that the matter distribution should be perfect fluid. This was proved by a theoretical investigation that the stellar configuration may contain the anisotropic pressure ( p r = p t ) when the matter density is very high, in particular, if it is more than the nuclear density [2][3][4][5][6][7]. This anistropic pressure introduce an anisotropy inside stellar configuration. Then the solution of Einstein's field equation in the presence of anisotropy can be determined easily as compared to perfect fluid matter distribution. In this connection, Many pioneering works have been done by several authors in different context . The several authors have been also obtained anisotropic and charged solution using embedding class one condition which is a much popular methodology from last few years . The presence of the electric charge in the static charged fluid compact object creates a repulsive force that averts the gravitational collapse [61,62]. Therefore, we can avoid the singularities by introducing the electric field inside the matter distribution. This repulsive force, due to the electric charge, balances the equilibrium for the dust distribution against the gravitational collapse. Recently, several authors have discovered the solution of the Einstein-Maxwell equations in a different framework which describe strange quark stars, charged black hole, and other astrophysical compact objects [63][64][65][66][67][68][69][70][71]. Varela et al. [72] studied self-gravitating, charged, anisotropic fluid in solving Einstein-Maxwell equations. In doing so they considered Krori and Barua's [73] metric potentials and linear or nonlinear equation of state. Arbanil et al. [74] studied a class of charged compact fluid spheres obeying the polytropic equation of state. They studied the Oppenheimer-Volkoff limit, Buchdahl limit for the charged polytropic spheres, and they found that the limit is extremal as it is a quasi black hole. So, as a result, we can say that the study of the charged fluid sphere is very interesting. That's why many researchers obtained exact solutions of Einstein-Maxwell field equations by electrifying some well-known uncharged fluid spheres e.g. Heintzmann's [75] solution by Pant et al. [76], Durgapal IV [77] by Pant and Rajasekhara [78], Durgapal V [77] by Gupta and Maurya [79], Pant I [80] solution by Maurya and Gupta [81]. The positive values of the charge parameter K describe completely the interior of the charged super-dense astrophysical object. One may consult the following works of literature [82][83][84][85][86][87][88] to study the relativistic compact stellar system more precisely with an electric charge. In this connection, recently the anisotropic charged and uncharged models in modified f (R, T ) theory gravity have been studied [89][90][91][92]. On the other hand, last few years the gravitational decoupling using minimal geometric deformation (MGD) approach has attracted the researchers. In this directions, few well recognize solutions of the Einstein field equations have been investigated in the framework of charged, uncharged and anisotropic matter distribution by using MGD . This MGD technique was developed to deform Schwarzschild space-time [114,115] into the Randall-Sundrum framework [116,117]. Ovelle et al. [118] and Contreras and Bargueno [119] have extended the black hole in 3 + 1 (Schwarzschild outer space-time) and BTZ black hole in 2 + 1 dimensions using the MGD, respectively. The technique for determining an anisotropic solution from any isotropic solution and inverse problem in the black hole context using the MGD is investigated by Contreras and his collaborators [120,121]. Recently, Ovelle [122] has proposed the extended case of MGD to solve the Einstein-Maxwell equations. In this connection, the extended gravitational decoupling in 2 + 1 dimensional space-time with cosmological context was discovered in [128]. This MGD approach has been also applied in other areas like Klein-Gordon scalar fields as an extra matter content [125], in cloud of strings [124], an extension of isotropic coordinates [126], ultra-compact Schwarzschild star or gravastar [127], as well as extended Durgapal isotropic model in the context of class one space-time for charged anisotropic matter distribution [129]. More recently MGD approach was used to discover the higher dimensional compact objects [130] and extent in the context of Lovelock [131] and modified f (R, T ) [91] gravity theories as well as in cosmological problems [132]. Moreover, through this GD approach, we can extend any well behaved isotropic solution (as discussed above) of the Einstein field equation into anisotropic or charged, as well as both anisotropic-charged domains by adding extra source θ i j in the original energy-momentum tensor T i j or by defining the action for total energy tensor. The presence of this extra source θ i j in the system can create an anisotropy or electric charge inside the matter distribution.
In the present article, we have applied the CGD approach by defining the modified action for total energy momentum tensor for the charged matter distribution which includes the extra source θ i j . In this situation, the field equations for complete or original system will contain ten unknown functions for charged perfect fluid matter distribution. By taking this into account, in the present article we solve these field equations using gravitational decoupling via complete geometric deformation (CGD) approach by transforming the metric potentials as: ξ → ν = ξ + α h(r ) and μ → e −λ = μ + α f (r ). In this way, we arrive two systems of equations namely Einstein-Maxwell's system for perfect fluid and quasi-Einstein system corresponding to the source T i j and θ i j , respectively. However, we suppose this θ i j will produce anisotropy inside the stellar structure. Then, we need to solve the Einstein-Maxwell's system of equations for anisotropic matter distribution. To do this, first, we solve Einstein-Maxwell's system by using well known Tolman-Kuchowicz gravitational potentials in the context of perfect fluid matter distribution. Then for finding the solution of a second system (quasi-Einstein system) we use the an equation of state (EOS), that relates the θ components, and specific physically motivated ansatz for deformation function f (r ). In this way, we obtain a well-behaved expression for the deformation function h(r ). After solving both systems individually, we combine both solutions and discussed several physical analysis in order to get the viability of this methodology.
The structure of the article is organized as follows: The Sect. 2 is the Einstein's-Maxwell field equations for the gravitationally decoupled system. This section is divided into three sections namely as: (2.1) the action for gravitational decoupling system, (2.2) basic stellar equation for the decoupled system, and (2.3) the gravitational decoupling approach and stellar field equations for T i j and θ i j . The matching conditions are performed to determine all the constants in Sect. 3. In the Sect. 4, we obtain a gravitational decoupling solution for the charged anisotropic compact star model by using CGD approach. The physical analysis has been done in Sect. 5 where we discussed regularity condition and causality condition in order to put the bound on the constants. The Sect. 6 contains the detailed discussion on the compactness and redshift. We have performed the equilibrium and stability analysis of the solution in Sect. 7. The physical feasibilty of the energy conditions have been verified in Sect. 8. The final Sect. 9 has been made for discussion and conclusion.
2 The Einstein-Maxwell field equations for gravitationally decoupled system

The action for gravitational decoupled system
The modified action for gravitational decoupled system can be defined by introducing an extra Lagrangian density for an extra source as [122], where L M and L e denote the Lagrangian for matter field and electromagnetic field, respectively, while L θ is a Lagrangian density corresponding to additional source. However, the symbols R and g has their usual meanings. Now we define T i j as a energy-momentum tensor for Lagrangian matter field L M which can be written as, As we know that matter Lagrangian L M depends on only the components of metric tensor g i j and not on their derivatives, so we obtain On the other hand, we denoted E i j and θ i j as electromagnetic field tensor and extra source corresponding to Lagrangian density L e and L θ , respectively, which can be given by, Now by varying the action (1) with respect to the metric tensor g i j we obtain the general equations of motion for the decoupled charged system as, where, Here, G i j denotes the Einstein tensor and the relativistic units are to be taken G = c = 1. We define T i j (corresponding to perfect fluid matter distribution) and E i j as, Here ρ and p represent the matter density and pressure for charged matter distribution while u ν is a covariant component for the 4-velocity which fulfil u μ u μ = 1 and u μ ∇ μ u μ = 0.
On the other hand, the anti-symmetric electromagnetic field tensor F μν satisfies the Maxwell's field equations, where j i , denotes the electromagnetic four current vector, is given by, where σ denotes the charge density and defined as σ = e ν/2 J 0 (r ). Moreover, the non-vanishing component for the four-current static fluid matter distribution is J 4 . Therefore, the four-current component acts only along radial direction due to the spherical symmetry. Then the corresponding nonzero components for electromagnetic field tensor are F 01 and F 10 that are related as F 01 = −F 10 . Then the electric field along the radial direction r can be defined by these components F 01 and F 10 only.

Basic stellar equation for decoupled system
In order to describe interior spacetime for the spherically symmetric static stellar system, we take the following line element as, where the metric potentials ν and λ are depend on the radial coordinate r only. Then the Einstein's field equations for the static spherically symmetric spacetime (13) corresponding to decoupling system (6) can be written as, − e −λ 8π where, Then the linear combination of equations (14) -(16) yields the following conservation equation, However for the spherically symmetric line element (13), the non-vanishing anti-symmetric electric field components F 01 and F 10 can be given in terms of the electric charge as, where q(r ) denotes the electric charge contained within the compact star of radius r . By using relativistic Gauss law, the electric charge q(r ) as well as electric field E can be defined as, Then the components for T i j and E i j can be expressed as, Then the conservation equation (20) will take the following form as, On the other hand, it is clearly noted that the inclusion of new source θ i j in the system introduced the anisotropy inside the matter distribution, if θ 2 2 = θ 1 1 . In this situation, the total energy tensor T tot , given in Eq. (7), can be described as, Hence, the anisotropy factor is given as, 2.3 The gravitational decoupling approach and stellar field equations corresponding to T i j (for charge matter distribution) and extra source θ i j In this section, we will see that how the extended gravitational decoupling converts the field equations (14)- (16) in two separate systems which are known as the Einstein-Maxwell system and quasi-Einstein system, respectively. As the previous discussion, we take the Einstein-Maxwell system with the perfect matter distribution for T i j while the quasi-Einstein system corresponding to extra source θ i j . Now we apply the Ovalle [122] transformation in metric potentials as, where h(r ) and f (r ) represents the geometric deformation functions corresponding to the temporal and radial metric component. The coupling constant α is a real number. Moreover, the above transformation is the extended case of minimal geometric deformation (MGD) which is called a complete geometric deformation (CGD) or extended geometric deformation along with both radial and temporal components of the line element. By substituting the deformed metric functions (31) and (32) in the field equations of decoupled system (14)- (16) with Eqs. (24) and (25), we arrive at the following set of equations: (I) The Einstein's equations of the charged perfect fluid matter distribution for energy momentum tensor T i j are given as (corresponds to α = 0), along the conservation equation, The solution of the field equations (33)-(35), satisfies conservation equation (36) that can describe the internal structure of the compact object for the charged perfect fluid matter distribution. The corresponding solution can be given by the following line element, (37) where the interior mass function m 0 is a mass function of the charged matter distribution for the standard GR expressions (33)- (35), which can be given as (II) Now let us go on the parameter α to see the effects of the extra source θ i j on the charged perfect fluid solution {ξ, μ ρ, p, }. For this purpose we write the field equations for the quasi-Einstein system associated with the source θ i j as , and the linear combination of quasi-Einstein equations (39)-(41) yields the following conservation equation, where, F 1 and F 2 are given as, It is clearly noted that the presence of θ -sector in the system will introduce the anisotropy inside the stellar structure.

Matching conditions for the stellar structure
The matching condition is a critical part in the study of the stellar distributions at the surface of the star (r = R) between interior (r < R) and exterior (r > R) spacetime geometries.
In the present study, the interior stellar geometry is given by the extended geometric deformation line element, wherem(r ) is internal mass of the stellar structure for total energy tensor T tot i j . Now the inner metric (43) should be smoothly matched with an exterior spacetime geometry. In his regard, Ovalle [122] has proposed that the exterior spacetime for the extended gravitational decoupling corresponding to "Maxwell version" of the vacuum solution T i j = 0 can be given by the well-known Reissner-Nordstrom solution as, wherem(R) =M and q(R) = Q is the total mass and total charge for the compact object of the radius R, respectively. For smooth joining we apply the Israel-Darmois junction conditions procedure [133,134] to match the inner manifold M − with the external one M + at the boundary .
The procedure of joining both space-time at the boundary is known as the continuity of the first and second fundamental forms across the surface . The first fundamental form says that the intrinsic geometry described by the metric tensor g i j induced by M − and M + on the interface meets By writing the first fundamental in the explicit form as, While the continuity of the second fundamental can be described as, here r j is a unit vector. Using Eqs. (6) and (49) we can find which gives, Now we will explain by simple way, for the standard GR case the T i j = 0 is the vacuum solution for the region r > R.
Here r = R is the boundary of the stellar structure. But the exterior spacetime (r > R) may not be a vacuum anymore due to presence of new fields which is coming from the θ sector. This extra field introduce the anisotropy inside the self-gravitating system. Then the above boundary condition (51) will take the following final form, where p r (R) = p − r (R). The condition given by Eq. (52) is called the general expression for the second fundamental form connected with the Einstein field equations described by Eq. (6). Now we substitute the value of (θ 1 1 ) − (R) for the interior geometry from Eq. (40) into the condition (52). Then the second fundamental form (52) can be written as, where, f * (R) and h * (R) are geometric deformation functions for the outer Reissner-Nordstrom solution (44) due to the extra source θ μν . It is a important remark that if the exterior solution is given by the Reissner-Nordstrom solution (44) then we must have f * (r ) = 0 and h * (r ) = 0. Now from Eq. (54) we get, which can also read as, The conditions (47), (48) and (56) are the necessary and sufficient conditions for determining the arbitrary constants involve in the system.

Gravitational decoupling solution
As we see that system of equations (33)- (35) and (39)-(41) contains ten unknowns. Therefore, in order to solve the first system of equations (33)- (35) we consider the seed spacetime (37) corresponding to well-known Tolman-Kuchowicz solution as, here, a, b, B are the parameters of dimension l −2 , l −4 , and l −2 , respectively while C is a dimensionless constant. As we observe that both μ and ξ are non-singular and well behaved throughout within the stellar model. These forms of the gravitational potential yield a physically viable solution. From the Eqs. (33) and (34) together with pressure isotropy condition (using in Eqs. (34) and (35)), we obtain the expressions for ρ, p and q 2 /r 4 as, where, Now we focus on the second system of equations (39)- (41). As we see that this system involes five unknown function namely θ 0 0 , θ 1 1 , θ 2 2 , and two deformation functions f (r ) and h(r ). For solving of this system, we propose a linear equation of state (EOS) in θ along with the one deformation function f (r ). But It is important to note that the choice of f (r ) should depend upon the following points: 1. For positive α together with positive and increasing function f (r ), the growth of μ(r ) must be faster than the the deformation function α. f (r ) in order to preserve e λ(r ) and mass function to be positive and increasing. 2. For positive α, and negative decreasing function f (r ), the metric function e λ(r ) and mass function will increase throughout automtically.
Based on above discussions we take, Now using (39), (40) together with EOS (62) and (63), we obtain other deformation function h(r ) as, with Now we have completely determined the deformation functions f (r ) and h(r ), which are physically viable (Fig. 1). Then the gravitational decoupling solution for the total energy momentum tensor T tot i j can be given by following line element, Hence, the explicit form of the deformed gravitational potentials read as, e ν(r ) = e B r 2 +C e α h(r ) .
where the interior mass functionm(r ) for decoupled system (14)- (16) can be given by, Then from Eqs. (27), (38) and (68), we get the following relation, It is noted that the massesm(r ) and m 0 will be equal and describe the mass function for the standard GR case when α = 0. Now the expressions for the components of θ -sector are determined as, where the used coefficients in θ 2 2 (r ) are given in appendix due to long expressions.
By substituting the Eqs.  29), we find the expressions for the total density (ρ tot ) and total radial pressure ( p tot r ) and total tangential pressure ( p tot t ) as, and anisotropy factor can be calculated by the formula:

Physical analysis
To be a regular model, it requires that the spacetime should be free from any mathematical and geometrical as well as physical singularity. For this purpose, we need to check the variation of the deformed metric potentials within the stellar interior. From Eqs. (66) and (67) we obtain: e λ(0) = 1 and e ν(0) = e C . The variation of the metric potentials is presented in Fig. 1 (bottom right) which shows metric potentials are regular and positive throughout the model.
On the other hand the central values of total pressure ( p tot 0 ) and total density (ρ tot 0 ) can be given , Now, the solution should satisfy the Zeldovich condition i.e. the ratio of pressure-density should be less than unity, then In addition to above, the causality condition must satisfy everywhere within the star i.e. 0 < v 2 r < 1 and 0 < v 2 t < 1. Then we obtain the central values (v 2 r ) r =0 and (v 2 t ) r =0 and together with causality condition, we obtain Then the inequality (80) gives, while inequality (81) yields, Where 1 and 2 are given in the appendix. Now combine all the inequalities (77), (79), (82) and (83), we obtain an inequality that restricts B as, Also, the variation of total density (ρ tot ) and total radial and tangential pressures ( p tot r and p tot t ) are shown in Fig. 2, which shows that those are positive and decreasing away from centre. However the maximum values attain at the centre of the

Mass-radius ratio and surface redshift
It is required to discuss the maximum limit of the massradius ratio to describe the compactness of the stellar model. In the case of isotropic matter distribution, the maximum limit of the mass-radius ratio was proposed by Buchdahl's [139] in the framework of perfect fluid having decreasing energy density towards to boundary. This maximum limit for mass-radius is given as, where m(R) = M describes the total mass of the object in perfect fluid matter distribution, while the R represents the radius of the model which is obtained by taking the pressure to be zero on the surface. Moreover, the presence of an electric charge in the solution modifies Buchdahl's limit. In this case, Andreasson [140], and Bohmer and Harko [141] have provided the maximum and minimum limit for the mass-radius ratio, respectively as, where m 0 (R) = M 0 is the total mass of the compact object for the charged perfect fluid matter distribution. It is noted that the mass M 0 present in the Eq. (86) is not equal as the total mass appears in (85). This can be defined as,  Table 1 The numerical values of mass, radius, central pressure, central density, surface density, electric charge at surface and constants B and C of the different compact objects, namely PSRJ1614-2230 [135,136], 4U1608-52 [136,137], Cen X-3 [136,138] for fixed value of a = 0.00509km −2 , b = 0.00002km −4 , β = 1. 8 On the other hand, the gravitational mass (appears in Eq. (85)) and effective mass both will be same in the context of perfect fluid or anisotropic fluid matter distribution, which can written as, But the above situation is not same for the charged matter distribution. In this case, the effective mass for charged matter distribution can be given as, Now from the equations (87), (88) and (89), we observed the followings: (i) the total mass M 0 for charged stellar object will be more than the total mass M of the compact object corresponding to the prefect or anisotropic fluid matter distribution, (ii) the effective mass [M 0 ] eff of the charged compact stellar model will be same as the total mass M in context of prefect fluid or anisotropic fluid matter distribution. Of course, the gravitational mass for the electrically charged stellar object has more value than the uncharged perfect fluid stellar model. Moreover, the same scenario will happen for the present gravitational decoupling models 2 . Now our aim is to see that whether the compactness i.e. mass-radius ratio in the presence of gravitational decoupling will take more value than the without gravitational decoupling, and also it will go beyond to the above standard bounds or not?. Then from the Table 2, we see that the mass-radius ratio (M R ) for GD model is also lying within the range given in Eq. (86). But It has more value than the mass-radius ratio ( M 0 R ) in absence of anisotropy (means when α = 0). Then it can be concluded that the presence of θ -sector i.e. anistropy, in the system can produce more compact objects. On the other hand, we would like to mention here that the effective mass plays an important role to define the upper bound of the surface redshift z s for the compact object. As we have already discussed that the mass-radius ratio gives a significant bound to observe the 2 we refer the following notations which have been used inside the text to represent the mass function for different matter distribution: surface redshift z s value. This bound can be determined by the following formula, On the other hand the gravitational redshift inside the compact object for GD solution can be obtained as, We have shown the behavior of the gravitational redshift inside the star in Fig. 3. From this figure, we can see that the gravitational redshift is maximum at the centre and decreasing outward.
The explicit form of this conservation can be expressed as, which is same as Eq. (26). Also, the above Eq. (92) can be divided in different forces whose linear combination will balance the system and achieve a stable equilibrium of the solution. These forces can be written as: (i). gravitational force: . electric force: F e = 2 q q 8 π r 4 , and (iv). anisotropic force: F a = 2 r (θ 2 2 − θ 1 1 ). To find the variations of these forces we plot the Fig. 4 for modified TOV equation (92) for different values of α. From Fig. 4 we observe that gravitational force (solid lines) and hydrostatic force (long dashed-dot lines), and like anisotropic force (small dashed-dot lines) are increasing and achieve its maximum value at the point within the stellar interior and then start decreases in each case. while other forces like electric force (long dashed lines) are monotonically increasing for throughout within the stellar interior. Another important point is that the electric force (F e ) plays a major role to balance the system near to surface of the objects while ansitropic force introduces a less impact on the system. As we see that the gravitational forces can be balanced by joint action of all other forces like hydrostatic force, electric force and anisotropic force such that F g + F h + F e + F a = 0. In this way, we achieved the stable equilibrium of the obtained each stellar model.
After analyzing the hydrostatic equilibrium under different forces, it is also required to check the stability analysis of the stellar model. To check this analysis, we use Abreu's criteria [38] which has been initiated by Herrera's cracking concept [142]. According to the Abreu's criterion, the stellar model will be stable if the subliminal radial and tangential sound speeds satisfy the following inequalities (Fig. 5), The left panel has been plotted for radial speed (v 2 r ) and right panel for the tangential speed (v 2 t ) verses radial coordinate r/R of different compact objects as PSR J1614-2230 (black curve for α = 0.001), 4U1608-52 (green curve for α = 0.0012), and Cen X-3 (red curve for α = 0.0014). The numerical values of the free parameters a, b, β and γ are same as used in Fig. 4 0 In order to verify these above inequalities, first, we need to check whether the stellar model is satisfying causality condition or not?. Then from Fig. 5, we see that the radial and tangential subliminal speed of sounds are satisfying the causality conditions i.e. v 2 r < 1 (top figure) and v 2 t < 1 (bottom figure) throughout within stellar object (where the speed of light is taken as unity i.e. c = 1). But we mention an interesting point that both radial and tangential velocities of sound are decreasing monotonically and radial velocity (v 2 r ) is always greater than the tangential velocity (v 2 r ) throughout within the stellar compact object for each chosen value of α, which can be analyzed from Fig. (5). After verifying the causality condition, we check the stability condition of the stellar compact model through the Eq. (93). Form Fig. 6, we observe that stability factors v 2 t − v 2 r (dashed lines) and v 2 t − v 2 r (solid lines) are lying within the intervals [−1, 0] and [0, 1]. This implies that the radial velocity dominates the tangential velocity everywhere inside the stellar interior. Also, there is no cracking within the star. Therefore, we conclude that the obtained self-gravitating charged compact star models are stable.

Energy conditions
It is very important to check the feasibility of some inequalities corresponding to the stress-energy tensor. For this purpose, we study these inequalities so-called energy conditions for representing physically realistic matter configuration. The respective energy conditions viz., the null energy condition (NEC), strong energy condition (SEC) and weak energy condition (WEC) are defined as where T μν l μ ∈ nonspace-like vector where i ≡ (radial r, transver se t), l μ and t μ are timelike vector and null vector respectively. In order to verify the viable feasibility of the stress-energy tensor, we need to study whether the compact stellar structure is consistent with the inequalities (94)- (96) or not? For this purpose, we plot the Fig. 7 for the above energy conditions and we observe that our stellar compact objects are consistent with all the energy conditions and hence ratifies that the physical acceptability of gravitational decoupling solution for compact objects.  Fig. 7 The behavior of energy condition verses radial coordinate r/R. The description of the figures as follows: the solid line for ρ tot , dotted line for ρ tot + p tot r , small dashed lines for ρ tot + p tot t , and long dasheddotted lines for ρ tot + p tot t + 2 p tot t for different compact objects PSR J1614-2230 (black color), 4U1608-52 (green color), and Cen X-3 (red color)

Discussion and conclusions
In this article, a completely deformed anisotropic charged fluid solution for the compact star model has been investigated by applying gravitational decoupling by means of a geometric deformation approach. To find the gravitational decoupling solution for the compact star model, first, we write the modified action for the gravitational decoupling system. Then we define the equation of motion by varying this action along with the metric tensor g i j . In this way, we arrive the Einstein's field equation for the coupled system which corresponds to the total energy-momentum tensor T tot i j in the context of spherically symmetric spacetime. The corresponding Einstein's field equations for coupled system will be solved by using the gravitational decoupling (GD) using a complete geometric deformation approach (CGD) in which both gravitational potentials have bee deformed as ξ → ν = ξ + α h(r ) and μ → e −λ = μ + α f (r ), as proposed by Ovalle [122]. This GD approach transforms the original Einstein's field equations for coupled system into two individual subsystems corresponding to the new sources T i j and θ i j . Here we would like to mention that the energy tensor T i j is considered for the charged perfect fluid matter distribution while θ i j introduce the anisotropy in the system (quasi-Einstein system). To solve the system for T i j , we use well-defined spacetime given by Tolman-Kuchowicz and find the solution for the first system, which provides the gravitational potentials μ and ξ and electric field E 2 . Now still five unknowns namely θ 0 0 , θ 1 1 and θ 2 2 , h(r ) and f (r ) are remaining in order to describe the complete structure of the stellar structure. Since it is clear that we have three independent equations for determining these five unknown which implies that we have two degree of freedom. To find this, we solve the quasi-Einstein system of Eqs. (39)- (41) cor-responding to θ -sector by specifying the linear equation of state (EOS), θ 0 0 = β θ 1 1 + γ , and well motivated ansatz for f (r ). After solving of this EOS, we obtain another deformation function h(r ). In this process, we achieved a complete deformed charged anisotropic solution for the compact star model. Moreover, we use the well known Isrial-Dormois junction condition for determining the constant parameters. The physical properties of the solution are described through the graphical analysis. We would like to mention that all the plots have been made for the same values of a = 0.00509, b = 0.00002, β = 1.8 and γ = 0.0001 with different α. In this way we obtained three different compact stars: (i) PSR J1614-2230 for α = 0.001, R = 11.0143km, (ii) 4U1608-52 for α = 0.0012 and R = 10.514 km, and (iii) Cen X-3 with α = 0.0014 and R = 9.922km. Here, we would like to mention that Demorest et al. [135] have predicted the radius corresponding to the differential neutron star EOS. They found that the radius of the compact objects with 1.97 solar mass lies between 11 and 15 km. But for strange quark matter (SS EOS), the radius will be less for this compact object with the same 1.97 solar mass. Later on, Gangopadhyay and his collaborators [136] have predicted the radii corresponding to 12 different stars and fitted the refined mass measurement of 12 pulsars using strange star equation of state (SS EOS) and they found that the star PSR J1614-2230 with solar mass 1.97 solar has the radius, approx. 9.69 km. In the present paper, By varying α we have predicted the radius of the realistic compact objects, as above, for fixed mass value. Then we observe the following points: if α increases then we will get less compact object as well lower masses stars. But for higher values of α ≈ 0.15 either causality will violates or cracking will appears in the system. The Fig. 1 shows the behavior of the deformation functions h(r ), | f (r )|, electric charge q(r ) and gravitational potentials e ν(r ) and e λ(r ) . From this figure, we observe that the functions | f (r )| (as f (r ) is negative and monotonic decreasing throughout the star) and h(r ) are zero at the centre and increasing towards the boundary of the stellar object and the same scenario happens in q(r ) also. The amount of charge on the surface of the star in Coulomb unit as: (i) Q = 2.8578 × 10 20 C for PSR J1614-2230 with α = 0.0001, (ii) Q = 2.4628 × 10 20 C for 4U1608-52 with α = 0.0012, and (iii) Q = 2.0361 × 10 20 C for Cen X-3 with α = 0.0014. However, the e ν and e λ both are positive and increasing monotonically for all values of α and free from singularity. The above well-defined regular gravitational potentials describe a well behaved physically viable gravitational decoupling solution for compact objects. On the other hand the Fig. 2 also shows the behavior of total radial and tangential pressures, ( p tot r and ( p tot t ), total density (ρ tot ), and the anisotropy ( ) against the radial coordinate r/R. As we can see in this Fig. 2 that the total pressures and density are maximum at centre and decreasing monotonically away from centre, while the anisotropy is zero at centre and increasing monotonically towards the boundary. The radial pressure is vanishes of the boundary of star which decide the size of the compact model i.e. radius R while tangential pressure is not. We determined the total mass (M), and surface redshift (z s ) of each obtained compact object which are the most important physical features of the model. It was proposed that the massM has larger value as compared to the total mass M 0 , which means that anisotropy introduce more massive objects. We have presented the numerical values of the compactness  [63] and Bowers and Liang [8]. Moreover, the variation of gravitational redshift within the charged ansitropic compact stellar models are shown in Fig. 4. The gravitational redshift is maximum at centre and decreasing towards the surface boundary, and attains minimum at surface. On the other hand, we verified the stable equilibrium of the gravitational decoupling solution via different forces. To do so, we need to study the modified Tolman-Oppenheimer-Volkoff (TOV) equation for anisotropic charged matter distribution. We have presented distributions of all forces F g , F h , F a , and F e in Fig. 4. From this Fig. 4, we found that the anisotroic force F α and electric force F e acts along an outward direction. However, the combined impact of the hydrostatic force F h , anisotropic force F a and electric force F e will balance the gravitational force F g such that F g + F h + F a + F e = 0, which yields the stable equilibrium for the matter distribution T tot i j . Therefore, our obtained GD solution is in the equilibrium stage. Finally, we discussed the causality and stability of the GD solution. The Fig. 5 shows that the speed of sounds is decreasing throughout the stellar interior and less than the speed of light. Moreover, the radial velocity is always greater than the tangential one which can be predicted from Fig. 5. On the other hand, from Fig. 6, we also see that the values of the stability factors v 2 t −v 2 r (dashed lines) and v 2 r −v 2 t (solid lines) belong to the intervals [−1, 0] and [0, 1], respectively. Also, The anisotropic model satisfies all the energy conditions (see Fig. 7). Finally, we would like to mention that the obtained charged anisotropic solution satisfies all the mathematical and physical requirements which shows that the extended gravitational decoupling by means of a complete geometric deformation approach is very effective and significant tool for generalizing or finding the new solution of the Einstein's equations.
Acknowledgements S. K. Maurya acknowledges that this work is carried out under TRC project-BFP/RGP/CBS/19/099 of the Sultanate of Oman.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There is no external data associated with this manuscript.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .