Charged strange star model in Tolman–Kuchowicz spacetime in the background of 5D Einstein–Maxwell–Gauss–Bonnet gravity

In this article, we provide a new model of static charged anisotropic fluid sphere made of a charged perfect fluid in the context of 5D Einstein–Maxwell–Gauss–Bonnet (EMGB) gravity theory. To generate exact solutions of the EMGB field equations, we utilize the well-behaved Tolman–Kuchowicz (TK) ansatz together with a linear equation of state (EoS) of the form pr=βρ-γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_r=\beta \rho -\gamma $$\end{document}, (where β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} and γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} are constants). Here the exterior space-time is described by the EGB Schwarzschild metric. The Gauss–Bonnet Lagrangian term LGB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {L}_{GB}$$\end{document} is coupled with the Einstein–Hilbert action through the coupling constant α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}. When α→0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \rightarrow 0$$\end{document}, we obtain the general relativity (GR) results. Here we present the solution for the compact star candidate EXO 1785-248 with mass=(1.3±0.2)M⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=(1.3 \pm 0.2)M_{\odot }$$\end{document}; radius =10-1+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$= 10_{-1}^{+1}$$\end{document} km. respectively. We analyze the effect of this coupling constant α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} on the principal characteristics of our model, such as energy density, pressure components, anisotropy factor, sound speed etc. We compare these results with corresponding GR results. Moreover, we studied the hydrostatic equilibrium of the stellar system by using a modified Tolman–Oppenheimer–Volkoff (TOV) equation and the dynamical stability through the critical value of the radial adiabatic index.The mass-radius relationship is also established to determine the compactness factor and surface redshift of our model. In this way, the stellar model obtained here is found to satisfy the elementary physical requirements necessary for a physically viable stellar object.


I. INTRODUCTION
Owing to the difficulties encountered by the general theory of relativity (GTR) in explaining the anomalous behavior of gravitational events such as the accelerated expansion of the cosmos in a late-time [1,2], alternative or extended gravity theories have suddenly gained considerable importance.Conjecturing the presence of exotic matter fields, including quintessence fields (QFs), ghost fields (GFs), dark energy (DE), and dark matter (DM), to name a few, is one approach to solving this problem.There is currently no empirical evidence for these conjectures, but a variety of experiments are being carried out.In this concern, de Rham [3] proposes that the graviton is not massless but actually bears a small mass to explain the dark sector.This has several implications for physics, which have already been addressed previously in the literature.Reexamining the geometrical side of the field equations offers an alternative approach, meanwhile higher curvature impacts may have a role to play.Specifically, the Einstein-Gauss-Bonnet (EGB) theory has shown to be promising in this aspect and is hence widely investigated.It should be noted that the EGB is part of a more generic category of theories named Lovelock's polynomial Lagrangians which are the most comprehensive tensor theory yielding at most 2nd-order motion equations.The most common theory is owed to Horndeski [4] if one allows the Lagrangian to involve both tensor and scalar fields.The reality that the Gauss-Bonnet Lagrangian naturally manifests in the applicability of heterotic string theory at the low energy limit [5] provides another compelling argument in favor of the EGB theory.For inhomogeneous distributions of dust [6,7] and null dust [8], the causal structure of the singularities deviates from general relativity.The challenge we are studying seems to be whether gravitational effects with higher curvature play an important role in stellar evolution.In fact, Einstein's GTR forms the basis for the majority of our knowledge of discoveries in relativistic astrophysics.So it is therefore quite natural to wonder what consequence the higher curvature contributions have on stellar configurations if the GTR is to be replaced by a higher curvature theory that maintains 2nd order motion equations and shrinks to GTR in the solar system scale constraint.
Although 5-dimensional celestial bodies are not physically obtainable, their occurrence has not been excluded.Kaluza [11] and Klein [12] made the first advances in higher dimensional gravity by analyzing a 5-dimensional manifold and explaining the behavioral patterns of the electrodynamic field by four components of the metric tensor, ten by the typical 4-dimensional space-time manifold, and one additional dimension by a scalar field.Following that, modern works on brane-world cosmologies requiring higher dimensions have received significant coverage, despite the fact that its emphasis has since waned.The additional dimensions are usually explained by the fact that they are topologically coiled and extremely tiny in size.Keep in mind that the LHC experiment probed for additional dimensions on a large scale but was unsuccessful in finding any.Nevertheless, this does not exclude the possibility of additional dimensions at the microscopic level.Indeed, quantum field theory actually requires space-time dimensions of the size of Ten and Eleven.As we demonstrate in this paper, despite their small size, they may have a significant impact on certain features of the gravitational field.
For a long time, the existence of spherically symmetric static black hole (BH) solutions in the EGB theory paradigm has been studied extensively [13][14][15].Subsequently, several other factors have been explored for the GB solution in de Sitter (dS) and anti-de Sitter (AdS) space [16,17], including thermophysical properties related to the cosmological horizon and the BH horizon.Additionally, some authors have recently examined specific solutions relating to BHs in great detail (see e.eg.Refs.[18][19][20][21]).Considerable efforts have been made to study the geodesic motion of a test particle [22], the radius of photon spheres [23], the Hawking evaporation of AdS BHs [24], the phase transition of RN-AdS BHs [25], solutions of regular BHs [26], solutions of wormholes fulfilling energy conditions [27,28], and the gravitational collapse of an incoherent spherical dust cloud was proposed in [29][30][31][32].
The study of EGB gravity may also be important with respect to the prospect of overcoming some issues in strong field regimes.The study of cosmic structures may also offer important restrictions on the alternative theories of gravity that are being considered in this context.A significant number of mass-radius relations are accessible at the moment through electromagnetic measurements of extreme events like short gamma-ray bursts (SGRBs), and more recently through gravitational waves (GWs).It is worth noting that the theoretical building of the neutron star (NS) equation of state (EoS) is strongly constrained by the observational data for NSs with masses around 2 M ⊙ [33,34].Until now, the composition and structural properties of relativistic compact bodies are not well recognized in detail.Many solutions to Einstein's gravitational field equations defining the internal construction of relativistic compact bodies have now been derived, all assuming static and spherical symmetry.Due to the complexity of the Einstein field equations, however, it is typically difficult to find exact solutions in GR that are physically realizable.In light of modified theories of gravity, the situation gets more challenging.So, in order to arrive at exact solutions, researchers employ a number of mathematical methodologies.Many such concepts have been discussed, in particular, an algorithm that produces all regular static, spherically symmetric perfect-fluid solutions of Einstein's equations by selecting a single monotone function [68] and its enlargement to locally anisotropic fluids in [36].With compact objects, however, there is flexibility in selecting the appropriate interior solutions.Schwarzschild found the first exact solution to Einstein's field equations in 1916 [37], which improved our ability to predict many physical events.By using spherically symmetric perfect fluid solutions of Einstein's equations, Tolman [38] initially derived the simplest model to represent stellar interiors.Following that, several works-both theoretical and observational-on various aspects of modified gravity theories have been reported in recent literature .
The foregoing arguments encouraged us to make a thorough study of the stability and appropriateness of anisotropic solutions for relativistic compact stars.In this article, we are mainly interested in exploring the salient features of the static-charged anisotropic fluid sphere consisting of a charged perfect fluid in the framework of 5D EMGB gravity theory.We are using the well-behaved TK ansatz together with a linear EoS of the form p r = βρ − γ, (where β and γ are constants) to generate exact solutions of the EMGB field equations and analyze the effect of the Gauss-Bonnet Lagrangian term L GB which is coupled with the Einstein-Hilbert action through the coupling constant α on the principal physical characteristics of the stellar model.We would also like to point out that the characteristics of the matter configuration and higher dimensional effects are directly related.
The manuscript is organized in the following way: In Sects.II, III and IV, we are providing the relevant EGB gravity formalism as well as the EMGB field equations and their solutions, respectively.The discussion of matching conditions for constant evaluations can be found in Sect.V.In Sect.VI, we have been discussing significant results for the actual charged model and its stability in Sect.VI F. Finally, we are presenting our concluding remarks in Sect.VII.Throughout the paper, we have employed the essentially positive signature (−, +, +, +, +).

II. EINSTEIN-GAUSS-BONNET GRAVITY
To construct the field equations in EGB gravity, we need a modified action that differs from the Einstein scenario.In this study, we will be working on five dimensions.In five-dimensional spacetime with a matter field, the modified action for EGB gravity can be written as, where g signifies the determinant of the metric g ij , R displays the five-dimensional Ricci scalar, and Λ stands for the cosmological constant, whereas S matter denotes the action related to the matter field.The GB coupling constant α is correlated to the string tension in string theory and has length squared dimensions, representing ultraviolet corrections to Einstein's theory.It should be noted that according to Maedaa [30], the coupling constant α must be nonnegative for Minkowski spacetime to be stable.Consequently, we will constrain our study to the case where the GB coupling constant α is non-negative.The Ricci scalar, Ricci tensor, and Riemann curvatures are specifically combined in the GB term L GB , which is stated by We also incorporate the matter action S matter , which triggers the stress-energy tensor for the matter field, T ij .We now get the following equations of motion by varying the action (1) in relation to the metric g ij , Here, the Einstein tensor and the tensor containing the contributions from the GB component are denoted by the symbols G ij and H ij , respectively.These both tensors are explicitly expressed as The GB term does not exist for n ≤ 4 but actively contributes to the curvature of the spacetime for n ≥ 5 since it is a topological invariant in the 4-dimensional spacetime.In writing the Einstein-Maxwell-Gauss-Bonnet (EMGB) field equations, here we are using geometrized units, and thus we have taken G = c = 1 throughout the discussion.

III. FIELD EQUATIONS
The 5-dimensional line element for a static spherically symmetric spacetime has the standard form in coordinates (x i = t, r, θ, φ, ψ).By considering the comoving fluid velocity as u a = e −ν δ a 0 , the Einstein-Maxwell-Gauss-Bonnet (EMGB) field equation (3) yields the following set of independent equations in view of the metric (5), where ρ, p r and p t respectively denote the matter density, radial and transverse pressure of the fluid.Here E is the electric field intensity, σ is the proper charge density and κ = 8π.Note that a ′ denotes the differentiation with respect to the radial coordinate r.If q(r) represents the total charge contained within the 5− D sphere of radius r, then it can be defined by the relativistic Gauss law as

IV. SOLUTION OF THE FIELD EQUATIONS
To solve the above field equations ( 6)-( 8) we utilize the following ansatz where a, b, B, C are constants.These metric potentials conform to the well-known Tolman-Kuchowicz [9,10] spacetime.Using this Tolman-Kuchowicz ansatz given in Equations ( 11), the field equations ( 6)-( 8) can be rewritten as, The anisotropic factor ∆ = p t − p r assumes the following form where, Ψ = (1 + ar 2 + br 4 ), Along with TK-metric we also assume that the radial pressure p r is related to the matter density ρ by a linear equation of state, where β and γ are positive constants.Adding Eqs. ( 12)-( 13) we get, Solving Eqn.(17) with the help of Eqn. ( 16), the matter density(ρ) and radial pressure(p r ) are obtained as, By using the expression of ρ in Eqn. ( 6) the expression for the electric field intensity is obtained as, Consequently the expression for the transverse pressure p t is obtained as, V. EXTERIOR SPACETIME AND BOUNDARY CONDITION In order to determine the constant parameters i.e., a, b, B, and C for our proposed model it is necessary to match the interior space-time solution in a smooth way with an appropriate static and spherically symmetric exterior vacuum Schwarzschild solution.Since here we are investigating the interior anisotropic solution in 5-D Einstein-Gauss-Bonnet gravity, the most suitable and appropriate static exterior vacuum Schwarzschild solution is proposed by Boulware-Deser [67].This Einstein-Gauss-Bonnet-Schwarzschild solution is expressed by the following line element: where . The quantities M is the gravitational mass, Q represents the charge of the fluid as determined by an isolated observer at spatial infinity and K is an arbitrary constant.
Using continuity of the metric functions and their derivatives, namely g rr , g tt and ∂gtt ∂r across the boundary r = R we get, where, Now by solving Eqn.(23), we obtain the expression for b as a function of α and a as, where, Now from the condition p r (r = R) = 0 we get, where, Ψ R = (1 + aR 2 + bR 4 ) At the center of the star the electric field is obtained as, We know from regularity of electric charge that the electric field vanishes at the center of the fluid configuration, E 2 (r = 0) = 0. Thus by simultaneously solving the Eqns.p r (r = R) = 0 and E 2 (r = 0) = 0, we obtain the expressions for B and γ as functions of α, β, a and b as, Now by solving Eqn.(24), we obtain the expression for C in terms of B as, Here we have used the approximated mass and radius of the compact star candidate EXO 1785-248 to determine the values of b, C, B and γ.To design a realistic stellar model, we fix the values of β and a as β = 0.25 and a = 0.031 with K = 0.29 so as to obtain physical validity.Using these determined constants, we will now discuss the physical characteristics of our obtained solution.

VI. PHYSICAL ANALYSIS OF THE PRESENT CHARGED MODEL
The model must satisfy the following regularity and reality requirements within the stellar configuration in order to be a viable, regular, and stable charged anisotropic compact stellar model with a certain radius.The behavior of anisotropic solutions of our compact stellar model is investigated through graphical analysis of metric potentials, matter variables, energy conditions, causality conditions, Herrera's cracking approach, adiabatic index, compactness factor, and redshift parameters.

A. Regularity of metric co-efficient
In the survey of relativistic compact stars, the presence of singularities is recognized as a crucial problem.With this in mind, we are investigating the nature of the metric constituents in the core of these compact relativistic stars.For the solutions to be physically realizable [68], the metric constituents inside the compact stars must be regular and positive.From Eq. ( 11), we can see that This gives d dr [e ν ] r=0 = 0 and d 2 dr 2 [e ν ] r=0 > 0 at the center.We can also notice in Eq. ( 11) that the other metric element ensures the form e λ r =0 = 1 + O(r 2 ) near the center.Therefore, the metric components considered are non-singular at r = 0 and monotonically increasing.They thus satisfy all the conditions required to be physically acceptable which can be verified from Fig. (1).

B. Regularity of pressure and density
To reinforce the physical validity of our solution, the central values of the other physical quantities, namely matter density, and pressure must also be finite and non-singular to show their finiteness at the center.In this context, we found the central values of the physical quantities at the center r = 0. Firstly, the central density is obtained as, and secondly, the expression for central pressure is obtained as, The radial and tangential pressures are equal and always positive at the center, where r = 0.A graphical representation of the requirements of the star enables for testing if they are always satisfied.

C. Evolution of fluid parameters
Here we explore the viable properties of charged anisotropic compact stellar configuration and describe its evolution graphically.Furthermore, the comparison of our theoretical findings with empirical data could provide strong evidence in favor of the considered 5-D EMGB model.We start by analyzing some realistic characteristics of the stellar configurations such as the matter density, radial and lateral pressures along with their derivatives, the anisotropy factor, and the trace profile through graphs.We acknowledge that due to the extremely dense stellar profile, the matter contents should be higher in the stellar interior.In this respect, we have shown in Fig. 2 that the behavior of the physical parameters, namely ρ, p r and p t is positive and decreases monotonically towards the stellar surface, which shows a very compact profile of the considered stellar configuration.Then, we analyzed the density and pressure gradients viz., dρ/dr, dp r /dr, dp t /dr and demonstrated that their behaviors are negative and regular, which ensures the presence of an intensive stellar object in the 5-D EMGB model, as illustrated in Fig. 3. It's worth noting that the matter variables have high values compared to the GTR.
An interesting key factor when analyzing the matter configuration is the pressure anisotropy, whose direction relies on the pressure components.When ∆ = p t − p r < 0, the anisotropy is negative, indicating that the pressure is oriented internally, while when ∆ = p t − p r > 0, the anisotropy is positive, indicating that the pressure is oriented externally.The graphical visualization of the anisotropy factor ∆ is yielded in Fig. 4 (left panel), which shows that the current stellar model has a positive anisotropy factor ∆, where p t > p r then ∆ > 0. Therefore, the stellar structure is disclosed to a repulsive force that offsets the gravitational gradient, this phenomenon enables the production of a more compact stellar configuration.The radial and transverse pressures coincide at the center, i.e., ∆ vanishes at r = 0. Additionally, the radial and transverse pressure values diverge with increasing radius, which causes the anisotropy to keep moving in the direction of the surface of the spherical stellar configuration.Besides, Fig. 4 (right panel) demonstrates that the trace profile i.e., ( pr +2pt ρ ) is positive with increasing behavior throughout the interior of the object.This result advocates the hefty profile of stellar matter variables, thereby representing compact environments of the considered stellar configuration. ) are plotted against the radial distance r.

D. Electric field
It is generally recognized that the regularity of the electric charge must be effectively zero at the center and positive for r > 0, which aids in avoiding gravitational collapse.To confirm this regularity, we have presented in Fig. 5 (left panel), the trend in the electric field.In this context, we notice that the electric field vanishes at the stellar core and is positive throughout the stellar configuration, but it takes a maximum value within the stellar configuration at r ≈ 5.5 km instead of the stellar surface.It is well-known that the stellar core can become unstable under the effect of strong electric fields.On the other hand, quasi-static equilibrium states can be generated by the presence of charges up to 10 20 coulombs.Due to pair production within the stellar configuration caused by very high charge densities and extremely strong electric fields, the core of the stellar configuration becomes unstable.Next, Fig. 5 (right panel) illustrates how the charge density behaved for the hypothetical stellar configuration.The charge density is consequently seen to be non-singular at the stellar configuration's center, decreases monotonically to the stellar configuration's surface layer, and attains a non-zero value at the stellar frontier.

E. Energy conditions
An ordinary matter can be distinguished from an exotic fluid on the basis of the energy conditions (ECs) on matter variables.These restrictions give detailed information about the physical viability of the solutions.Strong, weak, null, and dominant energy conditions (SEC, WEC, NEC, and DEC) for charged stellar configuration are given by the following inequalities, Fig. 6 illustrates the graphical behavior of these energy conditions which shows that our system is consistent with all these ECs.This guarantees that the charged matter that composes our stellar model is realistic and physically viable.

F. Stability of the stellar model
In this section, we will investigate the stability of our stellar model by using the two methods given by (i) Causality condition via Herrera's cracking concept and (ii) determining the relativistic adiabatic index.

Causality condition via Herrera's cracking concept
A very important condition to be a physically acceptable stellar model is the causality condition, which states that the square of the radial sound velocity (v 2 r = dpr dρ ) and that for tangential sound velocity (v 2 t = dpt dρ ) for anisotropic compact object should be less than unity everywhere [72], i.e., In other words, these conditions state that the velocity components of sound must be less than the speed of light.
Since we have assumed a linear equation of state for radial pressure p r , then from Eqn. ( 16), we obtain v 2 r = β.Therefore, the radial velocity component is constant throughout the stellar model and totally depends on parameter β, while tangential velocity component varies with parameter α.The expression of v 2 t in EMGB gravity is highly non-linear and complicated that has to be analyzed numerically.From Fig. 7, it is very clear that the values of all velocity components lie within [0, 1].
On the other hand, Herrera et al [69][70][71] proposed the cracking concept for self-gravitating compact stellar model by considering perfect fluid and anisotropic matter distributions.Following this concept, we have checked the cracking condition |v 2 t − v 2 r | ≤ 1 and thus from Fig. (7) we find that the present compact stellar model is potentially stable by satisfying both the causality conditions and the Herrera's cracking concept.

Relativistic adiabatic index
The adiabatic index, which measures the stiffness of the matter, is crucial for understanding the dynamical stability of a relativistic compact object.Chandrasekhar first proposed the dynamical stability of stellar systems against infinitesimal radial adiabatic perturbation [75].In his pioneering work, he predicted that for stability of a stellar system the value of adiabatic index Γ should be greater that 4/3.Several researchers have tested the predictions of Chandrasekhar for both isotropic and anisotropic stellar objects [76][77][78][79].The adiabatic index is defined as, We have successfully plotted the behaviour of the adiabatic index Γ in Fig. (8) which shows that Γ is greater than 4/3 for our proposed model.This confirms the stability of this stellar system.

G. Hydrostatic Equilibrium via Modified TOV Equation
The four different forces, namely gravitational force, hydrostatics force, electric force and anisotropic force are acting in our present model.In this section, we will investigate the hydrostatic equilibrium of our proposed model via the modified Tolman-Oppenheimer-Volkoff (TOV) equation under different forces.To analyze, we write the general form of the modified TOV equation as follows: proposed by Tolman-Oppenheimer-Volkoff and named as TOV equation.
Where M G (r) is the gravitational mass within the radius r derived from the Tolman-Whittaker formula and the Einstein's field equations and is defined by Using the expression for M G (r) in equation ( 41), we get The above expression may also be written as where F g , F h , F a , and F e represent the gravitational, hydrostatic-gradient, anisotropic and electric force respectively given by the following expression.
H. Some physical implications of the solutions for a toy stellar model we can readily measure the mass function inside the radius R by computing the integral connected directly to the density level using the following definition, conversely employing the metric function suggested by On the other hand, using the formula suggested by Böhmer and Harko [73], it is easy to determine the lower and upper limits of the mass-radius ratio, In addition, one may get the effective mass for the distribution of charged matter by, We also signify compactification for the effective mass-to-radius ratio, u, which is given as follows, This serves as the basis for evaluating the surface redshift function z s , which is stated as, These four quantities' behaviors, namely, mass function, compactness factor, gravitational redshift, and surface redshift, are displayed in Fig. 10.The graph demonstrates that each physically significant variable meets the criteria of a realizable stellar configuration.It should be noted that the mass function at the center is guaranteed to be regular.We can also see that the three quantities, namely mass function, compactness factor, and surface redshift are monotonically increasing functions, but the quantity: gravitational redshift is a monotonically decreasing function as a function of the radius r and also all these quantities are positive inside the stellar body.In particular, the plot demonstrates that the mass-radius ratio and its lower and upper limits for the compact stellar candidate under consideration in this research comply completely with the criteria listed in Table I, which may be confirmed in [73].However, the authors [73,74] contend that the surface redshift of an anisotropic fluid sphere is hypothesized to be less than z s ≤ 5 or z s ≤ 5.211.Hence, under these restrictions, the surface redshift is fulfilled everywhere within the stellar configuration, which demonstrates the viability of our stellar model.Moreover, we analyzed whether the charged matter is realistic or exotic based on various constraints applicable to EoS components such as 0 ≤ w r < 1, & 0 < w t < 1.These matter-related constraints should be maintained for a realistically charged and uncharged matter composition.The EoS has the following mathematical forms, It is straightforward to establish that our model of the heavenly objects faithfully reflects the realistic nature of the composition of matter under the radial and transverse components of the EoS.These w r , & w t profiles' graphs are displayed in Fig. 11.The evolution of radial and transverse pressures versus density is also shown in Fig. 11 for different values of α.In this paper, we have extensively explored the viability and stability of a static charged anisotropic fluid sphere consisting of a charged perfect fluid in the framework of 5D EMGB gravity theory.What's more, we have employed the well-behaved TK ansatz together with a linear EoS of the form p r = βρ − γ, (where β and γ are constants) to generate exact solutions of the EMGB field equations and analyze the effect of the Gauss-Bonnet Lagrangian term L GB which is coupled with the Einstein-Hilbert action through the coupling constant α on the main physical properties of the stellar model.
The eventual model has been thoroughly analyzed to identify whether it conforms to strict regularity and stability conditions.In addition, the causality condition via Herrera's cracking concept and the relativistic adiabatic index was deemed to be satisfactory.The adherence of the principal features such as energy density, pressure components, anisotropy factor, sound velocity, hydrostatic equilibrium of the stellar system via a modified TOV equation, massradius relation, compactness factor, surface redshift with the recognized physical behavior for different values of coupling constant α has been confirmed by graphical plots using the stated parametric values.Furthermore, our analysis demonstrates that matter configuration attributes and higher dimensional effects are intimately associated.The interior spacetime and the exterior space-time described by the EGB Schwarzschild metric were also matched to settle all integration constants that arose along the way.
This study has shown the meaning and significance of the Gauss-Bonnet Lagrangian term L GB which is coupled with the Einstein-Hilbert action through the coupling constant α together with the TK ansatz and a linear EoS in the development of astrophysical models that harmonize with observed data.

12 FIG. 1 :
FIG.1: Variation of metric functions e λ and e ν with respect to r for EXO 1785-248.

FIG. 6 :
FIG.6: Variation of energy conditions with respect to the radial distance r.

12 FIG. 7 :
FIG. 7: Components of sound velocity: (Left) v 2 r is plotted against r, (Middle) v 2 t is shown against the radial distance r and (Right) the stability factor |v 2 t − v 2 r | is plotted against r.

12 FIG. 8 :
FIG.8: Behavior of the adiabatic index Γ as a function of the radial distance r.

FIG. 9 :
FIG.9: Variations of different forces, such as the gravitational force Fg, hydrostatic-gradient force F h , anisotropic force Fa and electric force Fe respectively versus radial coordinate r for different values of α.

12 FIG. 11 :
FIG.11:(Left) EoS parameter ωr is shown against the radial distance r, (Middle) EoS parameter ωt is shown against the radial distance r, (Right) Relation between pr and ρ, (Bottom) Relation between pt and ρ.

TABLE I :
Comparative study of lower bound, mass-radius ratio, compactness, upper bound, and surface red-shift of the stellar toy model for different values of α α Lower bound Mass-radius compactness (u) Upper bound Surface redshift Q 2 (18R 2 +Q 2 ) 2R 2 (12R 2 +Q 2 )