Exploring wormhole solutions with global monopole charge in the context of $f(Q)$ gravity

This study explores the potential existence of traversable wormholes influenced by a global monopole charge within the $f(Q)$ gravity framework. To elucidate the characteristics of these wormholes, we conducted a comprehensive analysis of wormhole solutions employing three different forms of redshift function under a linear $f(Q)$ model. Wormhole shape functions were derived for barotropic, anisotropic, and isotropic Equations of State (EoS) cases. However, in the isotropic EoS case, the calculated shape function failed to satisfy the asymptotic flatness condition. Additionally, we observed that our obtained shape functions adhered to the flaring-out conditions under an asymptotic background for the remaining EoS cases. Furthermore, we examined the energy conditions at the wormhole throat with a radius $r_0$. We noted the influences of the global monopole's parameter $\eta$, the EoS parameter $\omega$, and $n$ in violating energy conditions, particularly the null energy conditions. Finally, we conducted a stability analysis utilizing the Tolman-Oppenheimer-Volkov (TOV) equation and found that our obtained wormhole solution is stable.


I. Introduction
The wormhole concept was first proposed by L. Flamm [1] in 1916.Subsequently, Einstein and Rosen further explored the notion of a wormhole, introducing a hypothetical structure known as the Einstein-Rosen bridge [2].Wormholes are theoretical constructs with topological properties that provide a conduit for connecting distinct regions of Space-time.These structures are envisioned as tube-like formations, asymptotically flat at both ends and connected by a throat.Wormholes are categorized based on the nature of their throats, distinguishing between static and nonstatic variants.A wormhole with a constant throat radius is termed a static wormhole, while a non-static wormhole exhibits a variable throat radius.Fuller and Wheeler demonstrated that traversing through the Einstein-Rosen bridge, even for a photon, would lead to its instantaneous collapse upon formation [3].Subsequently, Morris et al. [4] proposed that exotic forms of matter threaded through a wormhole might stabilize it; however, the feasibility of such requirements remains uncertain.Finally, Morris and Thorne [5] presented static traversable wormholes, which offer a theoretical framework for interstellar travel and constitute exact solutions within the framework of General Relativity (GR).The authors demonstrated that a wormhole could be traversable if it possesses exotic matter with a minimal surface area that satisfies the flaring-out condition.This exotic matter is essential for constructing a wormhole in GR; indeed, in classical GR, wormhole solutions violate all energy conditions, especially null energy conditions [6].However, this type of hypothetical matter exhibits unusual properties within the framework of GR.In contrast, from the perspective of quantum gravity, it can arise as a natural consequence of fluctuations in the Space-time topology [7].Thus, minimizing the violation of the energy conditions or reducing the quantity of exotic matter at the throat is essential.
Kanti et al. made important advances in the research of traversable wormholes by proposing wormhole constructs within the framework of quadratic gravitational theories [8,9].This unique technique provides an intriguing concept in which gravity plays an important part in keeping the wormhole throat open, removing the necessity for exotic matter.Kuhfittig [10] has studied rotating axially symmetric wormholes by introducing time-dependent angular velocity to generalize static and spherically symmetric traversable wormholes.B öhmer et al. [11] investigated wormhole solutions by assuming a linear relationship between energy density and pressure.Barros and Lobo [12] developed wormhole solutions using three-form fields and explored their compatibility with weak and null energy conditions.Also, stability of thin shell wormholes with barotropic fluids present has been done in [13].Wormhole geometries have garnered considerable attention not only in modified gravity theories but also in higher-dimensional gravitational theories [14][15][16][17][18] and Kaluza-Klein gravity [19][20][21].These theoretical frameworks offer distinct advantages, such as eliminating the need for nonstandard fluids, which has been a significant driving force behind extensive research in modified gravity theories.Moreover, modifications to Einstein's gravity introduced additional degrees of freedom within the gravitational sector, presenting new avenues for addressing challenges such as dark energy and dark matter.Additionally, Falco et al. [22] concentrated on f (R) metric, f (Q) symmetric teleparallel, and f (T) teleparallel models to explore spherically symmetric static wormhole solutions.Recent literature [23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38] provides an insightful investigation of wormhole geometry under several modified gravity theories.These works contribute substantially to the ongoing debate about the theoretical underpinnings and observational implications associated with wormholes in alternative gravitational contexts.Jimenez et al. [39] proposed f (Q) gravity, which relies on the nonmetricity Q to explain gravitational interactions.Over the past few years, f (Q) gravity has been under strong observational investigation, with Lazkoz et al. [40] suggesting significant restrictions on its validity.Their study examined the plausibility of f (Q) gravity using a variety of observable datasets, including the Type Ia Supernovae, expansion rate, quasars, baryon acoustic oscillations, gamma-ray bursts, and cosmic microwave background distance.Mandal et al. [41] established the validity of f (Q) gravity models under energy conditions.They developed an embedding approach that incorporates the non-trivial contributions of the nonmetricity function into energy restrictions, confirming the consistency of f (Q) gravity with fundamental physical laws.These efforts contribute substantially to understanding the potential of f (Q) gravity as a convincing alternative to traditional gravitational theories, clearing the way for further investigation into the basic concept of gravity and its significance for cosmology.Also, Gadbail et al. [42] presented f (Q) gravity reconstruction from Friedmann-Laimatre-Robertson-Walker (FLRW) evolution, demonstrating exact ΛCDM expansion and showing how generic functions of Q require extra degrees of freedom in the matter part.Furthermore, the impact of the Generalized Uncertainty Principle (GUP) on the Casimir wormhole space-time in f (Q) gravity [43] has been analyzed systematically.In their investigation, the authors considered two distinct forms of f (Q) and meticulously derived both analytic and numerical solutions.These solutions were obtained by incorporating the influence of the GUP.Besides that, they investigated the ADM mass and used volume integral quantifiers to calculate the required amount of exotic matter at the throat of the wormhole.Also, Hassan et al. [44] studied how the Casimir effect affects wormhole geometry under f (Q) gravity.To improve experimental feasibility, they looked at three Casimir effect systems: (i) two parallel plates, (ii) two parallel cylindrical plates, and (iii) two spheres separated by a significant distance and discovered that certain arbitrary quantities violated classical energy conditions at the wormhole throat.Importantly, various intriguing astrophysical experiments have been carried out using the f (Q) gravity theory and its extended gravity.References [45][46][47][48][49][50][51][52] provide valuable resources for further investigation.
Monopoles were first proposed by Dirac [53] when he had shown by geometric arguments that if magnetic monopoles exist, then the charge of monopoles would be quantized.Even though the electromagnetic field is a rather simple abelian (U(1)) gauge field, it was later shown by G 't Hooft [54] and Polykov [55] that more complicated non-abelian gauge fields would also lead to similar types of monopoles.As we know, the hot Big Bang model of the universe in the Plank era had all the forces in equilibrium, and then the four forces decoupled; this gives a rather compelling argument that, indeed, we can find the monopoles from the early universe, which was produced by the spontaneous symmetry breaking.In this article, we rather focused on a simple Barriola valentine [56] type monopole (which forms when SO(3) group spontaneous symmetry breaks into SO(2) group) for which the effect of the gravitational field is well known.Moreover, Zloshchastiev [57] gave a comparative description of monopole in the quantum mechanical context.Also, De-Chang Da et al. [58] have studied the phenomenological aspects of the wormhole under the influence of global monopole.Further, Ahmed [59] studied a topologically charged fourdimensional wormhole and the energy conditions, including global monopole charge, and found that the energy-momentum tensor associated with this wormhole complies with the weak energy condition and the null energy condition.In reference [60], Shi et al. provided an exact solution for the nonlinear equation that describes the global monopole (GM) in flat space.Their study of the metric outside the GM revealed a repulsive gravitational field beyond the core and a solid angular deficit.Reference [61] demonstrated that the gravitational fields of GMs can cause matter clustering and anisotropies in the Cosmic Microwave Background (CMB).Additionally, reference [62] suggested that GMs might contribute to forming supermassive black holes.Studies [60,63] have shown that this topological defect induces a negative gravitational potential, leading to a repulsive gravitational field.Many researchers have highlighted the essential role of topological defects in the formation of cosmological structures and the evolution of the cosmos, as noted in references [64][65][66][67][68]. Recently, reference [69] investigated wormholes within the Milky Way galaxy that possess a global monopole charge.This article investigated whether such a monopole could cause the wormhole in modified f (Q) gravity.It is well known that wormholes can be sustained by global monopoles.Here, we have investigated whether the same holds in symmetric teleparallel gravity.While previous studies have examined wormholes in the context of monopoles [70][71][72][73], we provide a more comprehensive perspective by considering f (Q) gravity.
The organization of this paper is as follows: Section II outlines the criteria defining a traversable wormhole imbued with a Global Monopole and presents the formalism elucidating f (Q) gravity.In Section III, we derive the field equations within a linear model framework.Subsequently, in Section IV, three linear EoS cases are introduced, and the necessary conditions for the existence of a traversable wormhole are examined, along with a scrutiny of the energy conditions under three different choices of redshift function.Section V encompasses a stability analysis conducted utilizing the Tolman-Oppenheimer-Volkov (TOV) equation.Finally, Section VI consolidates the findings and engages in a comprehensive discussion concerning the implications of the results derived from this study.

II. Basic Criteria of a traversable wormhole with Global
Monopole and formalism of f (Q) gravity The wormhole metric in the Schwarzschild coordinates (t, r, θ, Φ), as defined by [5,6], is expressed as: (1) Here, b(r) denotes the shape function defining the geometry of the wormhole.The function ϕ(r) shows the redshift function associated with gravitational redshift.To render a wormhole traversable, the shape function b(r) must satisfy the flaring-out condition, as specified by (b − b ′ r)/b 2 > 0 [5].At the wormhole throat b(r 0 ) = r 0 , the condition b ′ (r 0 ) < 1 is imposed, where r 0 denotes the throat radius.Additionally, the asymptotic flatness condition demands that as r → ∞, the ratio b(r) r tends to 0. Furthermore, to avoid an event horizon, ϕ(r) must remain finite at all points.Adherence to these conditions ensures the possibility of exotic matter existence at the wormhole throat within the framework of Einstein's GR.Now, the formulation of the action for symmetric teleparallel gravity with a global monopole charge, derived from a four-dimensional action without considering the cosmological constant and minimally coupled triplet scalar field (c = G = 1), can be expressed as follows: (2) In the given context, f (Q) denotes an arbitrary function of Q, where Q represents the non-metricity scalar.L m signifies the matter Lagrangian density, L denotes the Lagrangian density of the monopole, and g stands for the determinant of the metric tensor g µν .The Lagrangian density of a self-coupling scalar triplet ϕ a is expressed as: In this expression, a = 1, 2, 3, η, and λ represent the gauge-symmetry breaking and self-interaction term scales, respectively.The field configuration corresponding to the monopole is described as follows: where the variable x a is defined as (rsinθcosϕ, rsinθsinϕ, rcosθ), ensuring that ∑ a x a x 2 = r 2 .The Lagrangian density can be expressed in terms of F(r) by utilizing the field configuration as (5) and for the field F(r), the Euler-Lagrangian equation is The energy-momentum tensor is constructed from Eq.
(3) as Hence, all four elements of the Energy-Momentum tensor can be found using Eq. ( 7): As demonstrated in Eq. ( 6), obtaining a precise analytical solution can be challenging.Therefore, it is often advantageous to approximate the region outside the wormhole to simplify the analysis and derive results.As a result, the components of the reduced energymomentum can be expressed as follows: Furthermore, the non-metricity tensor can be expressed by the following equation [39].
Additionally, the superpotential, also known as the nonmetricity conjugate, has the following formal definition: ) Also, traces of the non-metricity tensor can be provided by The non-metricity scalar is given by [39] The disformation tensor, denoted as L β µν , is defined as follows: The equations of motion for gravity can be obtained by varying the action concerning the metric tensor g µν .This equation is expressed as follows: where f Q = ∂ f ∂Q and T µν is the sum of the energymomentum tensor of anisotropic fluid part and matter field part.Therefore, it can be written as The elements of the energy-momentum for the anisotropic fluid are where ρ is the energy density, p r is the radial pressure, and p t is tangential pressure.The non-metricity scalar Q, associated with metric (1), is defined in accordance with the reference [74] as Therefore, the field equations governing f (Q) gravity for the wormhole, taking into account the Global Monopole Charge, can be expressed as: Utilizing these precise field equations, a comprehensive exploration of different wormhole solutions within the framework of f (Q) gravity models becomes feasible.

A. Energy conditions
The classical energy conditions, stemming from the Raychaudhuri equations, constitute a fundamental framework for exploring physically feasible matter configurations.These conditions, which encompass weak, null, strong, and dominant categories, play a crucial role in understanding the properties of matter within the context of GR.Among them, the null energy condition assumes particular significance, especially in the context of GR's treatment of wormhole solutions.The null energy condition holds pivotal importance due to its direct relationship with the energy density necessary to maintain the integrity of a wormhole's throat.Any departure from the null energy condition near the throat of a wormhole implies the existence of exotic matter characterized by negative energy density-an attribute not typically associated with conventional matter sources.Energy conditions function as constraints on the stressenergy tensor, delineating the distribution of matter and energy in space-time.They are represented as follows: • The weak energy condition (WEC) : ρ ≥ 0, ρ + p t ≥ 0, and ρ + p r ≥ 0.
• The strong energy condition (SEC) : ρ + p t ≥ 0, ρ + p r ≥ 0, and ρ + p r + 2p t ≥ 0. In conclusion, energy conditions provide crucial constraints on the behavior of matter in the Universe and play a crucial role in our exploration of wormholes.

III. Linear f (Q) model
In the present study, we adopt a linear functional form of f (Q) gravity, represented as: Here, α ̸ = 0 represents a model parameter.This value can easily be modified to regain GR when α = 1 and β = 0. Using this linear model, Solanki et al. [75] explored the late-time cosmic acceleration without including any dark energy component in the matter part.Avik De and Tee-How Loo's work [76] shows that this model maintains energy conservation.This model is used in [77] to analyze traversable wormholes inspired by non-commutative geometries with conformal symmetry, indicating the possibility of wormhole solutions with viable physical features.Moreover, this symmetric teleparallel model is also used in Casimir wormhole [43,44].Also, this model successfully explains strange stars, which is consistent with observation [78].Further, S. K. Maurya et al. [79] attempted to find an anisotropic solution for a compact star generated by this linear form of the f (Q) model.One can look into the following literature [80,81] for more extensive studies conducted in the context of symmetric teleparallel gravity.Therefore, the reduced field equations with an arbitrary redshift function ϕ(r) are stated as follows: (27) This study introduces three different forms of redshift function ϕ(r) as a component in deriving the wormhole solution.We introduce additional cases and study the energy conditions to obtain analytical solutions for the resulting field equations.These will be elaborated in the next Section IV.

IV. Linear equation of states
In this particular section, we consider three different EoS such as barotropic (p r = ωρ), anisotropic (p t = np r ), and isotropic (p t = p r ) to find the exact analytical solution (b(r)).These cases are as follows:

A. Barotropic EoS
In this section, we are delving into the construction of wormhole solutions by utilizing a linear barotropic Equation of State (EoS) as presented in [82]: Here, ω represents the EoS parameter.In late-time cosmic acceleration, the description of dark energy within the ΛCDM framework is distinguished by an EoS parameter of ω = −1, which has exhibited remarkable success thus far.Another extensively discussed model for time-dependent dark energy is the quintessence model, characterized by an EoS parameter ω > −1.Conversely, the least theoretically understood form of dark energy is phantom energy, identified by an EoS parameter ω < −1.Reference [83] explores asymptotically flat phantom wormhole solutions.Additionally, references [84,85] point out the considerable challenge of discovering wormhole solutions with a linear barotropic EoS within the context of non-linear models in teleparallel and symmetric teleparallel gravity.Also, Solanki et al. [86] successfully derived the precise wormhole solution within the framework of the non-linear f (R, L m ) model while maintaining a linear barotropic Equation of State (EoS).Nevertheless, we manage to derive an exact wormhole solution within the framework of our linear f (Q) model while adhering to a linear barotropic EoS using redshift functions as discussed in below Subsubsections IV A 1-IV A 3.
The function ϕ(r) = c [84,87], where c represents a constant, is widely regarded as an appropriate selection for a redshift function owing to its consistent behavior with respect to the radial coordinate r.Notably, when c equals zero, this function is commonly referred to as the tidal force.Upon solving equations ( 25) and ( 26) while incorporating the relation (28), we arrive at the following first-order differential equation given by b Upon integrating the aforementioned equation with the throat condition b(r 0 ) = r 0 , we get the explicit form of shape function as Now, we will delve into the visual representation of the shape function and explore the crucial conditions for the existence of a wormhole.To accomplish this, we will meticulously select appropriate parameters.We will showcase the graphical depiction of the asymptotic behavior of the shape function in Contour plot 1 for the EoS parameter ω (in the left plot) and global monopole parameter η (in the right plot).It confirms that as the radial distance increases, the ratio b(r) r tends towards 0. In addition, we verify the satisfaction of the flaringout condition b ′ (r 0 ) < 1 at the wormhole throat r = r 0 through Contour plot 2 for the same parameters ω and η.Here, we designate the wormhole throat at r 0 = 1.Also, we present the embedding diagram and its comprehensive visualization in Fig. 3.By employing equation (30) in equations ( 25)-( 27), one can obtain the expressions for the energy density, radial pressure, and tangential pressure as follows Now, we present the graphical representation of the energy density using Contour plot 4, illustrating its positively decreasing behavior throughout the Space-time for the parameters ω and η concerning the radial coordinate r.Additionally, NEC at the throat of wormhole r = r 0 is given by ) Moreover, we present the graphical illustration of this condition in Contour plots 5 and 6, providing a visual understanding of the radial and tangential pressures at the wormhole throat.Through Contour plot 5, it is clear that NEC for radial pressure shows negatively increasing behavior as radial distance increases.Also, the Contour plot 6 shows the positive behavior of NEC for tangential pressure.Overall, NEC is violated in this barotropic case, which may show the presence of exotic matter that maintains stability at the wormhole's throat.In addition, we extend our analysis toward SEC, where we observe satisfaction near the throat for both parameters and can be visualized through Contour plot 7.

ϕ(r)
The function ϕ(r) = log 1 + r 0 r [87] serves as another viable option as a redshift function, effectively avoiding the event horizon beyond the wormhole throat.Upon solving equations ( 25) and ( 26) while incorporating the relation (28), we deduce the following first-order differential equation: By integrating the aforementioned equation with the condition b(r 0 ) = r 0 , we obtain an explicit expression for the shape function as where 2 F 1 is the hypergeometric function and it can be It is evident that the above expression is not asymptotically flat, which means that the ratio b(r) r does not tend to 0 as r → ∞ and graphically presented in the Contour plot 8.As a result, it is challenging to obtain asymptotically flat wormhole solutions for barotropic EoS under a particular logarithmic redshift function for the f (Q) model given in Eq. (24).
The function ϕ(r) = 1 r [46,88] presents another suitable alternative as a redshift function, given its consistent behavior for r > 0, effectively circumventing the event horizon beyond the wormhole throat.Now, by again solving equations ( 25) and ( 26) while incorporating the relation (28) for the particular redshift function ϕ(r) = 1 r , the first-order differential equation is given by : b By integrating the above first order differentiable equation with the throat condition, we get the shape function as below Further, visual representation of the above shape function is shown in the Contour plot 9 and it is clear that the shape function obtained for the redshift function ϕ(r) = 1 r case is not asymptotically flat.

B. Anisotropic EoS
In this specific segment, we focus on exploring the anisotropic energy-momentum tensor fluid relation within the context of asymptotically flat wormhole solutions.Unlike scenarios where p r equals p t , we delve into cases where a disparity exists between these radial and tangential pressures.To elucidate this relationship between p r and p t , we refer to the formulation presented in prior works [89,90] as In this context, the parameter n assumes significance, as it governs the characteristics of the anisotropic fluid under examination.It is imperative to note that in our study, n cannot hold the value of one, denoted as n ̸ = 1.This restriction is crucial because if n were equal to one, it would simplify into that of a perfect fluid.Moreover, we study the case where n = 1 in the next subsection IV C. Hence, our analysis exclusively per-tains to scenarios where n deviates from unity, thereby retaining the complexities associated with anisotropic fluid behavior by employing redshift functions, as discussed in Subsubsections IV B 1-IV B 3.

ϕ(r) = c
In this specific subsubsection, we use Eqs.( 26) and ( 27) with consideration of the relationship (40), we obtain the subsequent first-order differential equation: Upon integrating the aforementioned equation with the throat condition b(r 0 ) = r 0 , we derive the explicit form Now, we transition into the visual analysis of the shape function, where we explore the pivotal conditions governing the existence of a wormhole.To facilitate this examination, we meticulously select appropriate parameters.We begin by presenting Fig. 10, showcasing the asymptotic behavior of the shape function with respect to the various values of parameter ω (depicted in the left plot) and η (displayed in the right plot).This graphical representation confirms the trend where, as the radial distance increases, the ratio b(r) r tends towards 0, thus showing the satisfaction of the asymptotic flatness condition.Furthermore, we ensure the fulfillment of the crucial flaring-out condition b ′ (r 0 ) < 1 at the wormhole throat r = r 0 , as demonstrated in Fig. 11, utilizing the various values of the same parameters ω and η.Ad-ditionally, we provide a comprehensive visualization of the embedding diagram in Fig. 12.By utilizing Eq. ( 42) within Eqs. ( 25)- (27), one can derive the expressions for the energy density, radial pressure, and tangential pressure as follows: Moving forward, we showcase the energy density's behavior using Fig. 13, depicting its gradual decreasing behavior throughout Space-time concerning the parameters n and η across the radial coordinate r.Moreover, it is imperative to account for the NEC concerning the radial and tangential pressures at the wormhole throat, which can be expressed as: Further, Figs. 14 and 15 visually elucidate the NEC's radial and tangential pressure aspects at the throat.Fig. 14 reveals the NEC's negative increase trend in radial pressure with increasing radial distance, while Fig. 15 highlights the NEC's positive trend in tangential pressure.The violation of NEC in this barotropic scenario suggests the likely presence of exotic matter required for throat stability.Extending our analysis to the SEC, Fig. 16 demonstrates its satisfaction near the throat for both parameters, providing a comprehensive visualization of the SEC's behavior.26) and ( 27) with the relation (40), we obtain the subsequent first-order differential equation as b  Further, we get the explicit form of the shape function by integrating the above equation as where, "E" is an exponential integral function and defined by t n dt.Moreover, a visual representation of the above shape function is shown in Fig. 17, and clearly, the asymptotic flatness condition is not satisfied here for this redshift choice.Also, from the above expression, it is clear that the shape function is not defined for n = −1, 1  2 .Hence, the plot for n = −1 is not visible in the Fig.

ϕ(r) = 1 r
By employing this particular form of the redshift function ϕ(r) = 1 r alongside Eqs. ( 26) and ( 27), incorporating the relation (40), we derive the following subsequent first-order differential equation: Furthermore, the above-mentioned differential equation does not have a solution in closed form.Even after extensive attempts using computational tools like Mathematica, no closed-form solution has been attained, indicating that it cannot be expressed in terms of elementary functions.Consequently, we are unable to make any conclusive remarks regarding the condition of asymptotic flatness.

C. Isotropic EoS
Typically, the energy-momentum tensor of asymptotically flat wormholes in GR exhibits anisotropy, where the radial pressure p r differs from the tangential pressure p t .We have previously delved into wormhole solutions for this anisotropic EoS (IV B) scenario.In this section, we shift our focus to investigating wormhole solutions under the following isotropic relation: Now, by utilizing Eqs. ( 26) and ( 27) while incorporating the relationship (51), we arrive at the following firstorder differential equation: By integrating the above-mentioned equation with the condition b(r 0 ) = r 0 , we derive an explicit expression of the shape function as It is evident that the above expression is not asymptotically flat, which means that the ratio b(r) r does not tend to 0 as r → ∞.As a result, it is challenging to obtain asymptotically flat wormhole solutions for isotropic pressure under a constant redshift function for the f (Q) model given in Eq. (24).This type of result can be found in the reference [86].Moreover, our findings are consistent with those discussed in Subsubsections (IV A 2-IV A 3) and (IV B 2-IV B 3) when examining the redshift functions ϕ(r) = log 1 + r 0 r and ϕ(r) = 1 r .Specifically, we determined that the shape functions derived from these redshift functions do not meet the criteria for asymptotic flatness.

V. Equilibrium Conditions
This section utilizes the generalized Tolman-Oppenheimer-Volkov (TOV) equation [91][92][93] to assess the stability of our wormhole solutions.The generalized TOV equation can be expressed as where ϖ = 2ϕ(r).Depending on the anisotropic matter distribution, the hydrostatic, gravitational, and anisotropic forces are outlined as follows: ) To achieve equilibrium in the wormhole solutions, it is essential for the equation F H + F G + F A = 0 to hold.In this study, we have assumed a constant redshift function ϕ(r), which implies that the gravitational contribution F G vanishes in the equilibrium equation.Thus, the equilibrium equation simplifies to: Utilizing Eqs. ( 31), (32), and (55), we obtain the following equations for the hydrostatic and anisotropic forces for the linear model ( 24) using barotropic EoS case, i.e., p r = ωρ as Also, for the anisotropic EoS case, i.e., p t = np r , considering Eqs. ( 44), (45), and (55), the equations for the hydrostatic and anisotropic forces are given as The hydrostatic and anisotropic forces are visually represented in Figs.18 and 19 for both the cases.A notable observation from the graphs is that despite their opposite directions, both forces exhibit similar patterns.This symmetrical balance between the forces suggests the stability of the wormhole solutions obtained in our study.

VI. Conclusion
Wormholes serve as hypothetical tunnels connecting two distinct spatial regions separated by a spacelike interval within the same Space-time manifold.Despite being a fascinating theoretical concept, they have yet to be observed.In GR, the wormhole solutions of interest, particularly those that are stable and traversable, encounter a known issue where energy conditions are violated.The necessity for exotic matter, as dictated by standard model extensions within classical GR, poses a  significant hurdle.To address this challenge, researchers have turned to modified versions of gravity.Exploring wormholes within modified theories of gravity offers a promising avenue to overcome the issue of energy condition violations.This approach enables the derivation of stable solutions without compromising the fundamental energy conditions.By delving into modified gravity frameworks, researchers aim to refine our understanding of wormholes and potentially pave the way for their realization in the Universe.In our current investigation, we have thoroughly explored wormhole solutions within the framework of f (Q) gravity under the influence of a global monopole charge.This extensive study involved the utilization of three distinct linear EoS juxtaposed against the computed energy density within the f (Q) gravity framework.Our analysis has resulted in a collection of viable and internally consistent solutions to the field equations, specifically tailored to describe a Morris and Thorne wormhole (Eq.( 1)) with linear EoS.These linear EoS comprise the barotropic (Eq.( 28)), anisotropic (Eq.( 40)), and isotropic (Eq.( 51)) cases.The functional forms of the shape functions have been derived, showcasing re-markable accuracy within the specified conditions.The key highlights of our investigation are succinctly outlined as follows: 1.In this investigation, we have focused on the linear form of f (Q) (see Eq. ( 24)) model with three different forms of redshift function.
2. It is found that the choice of redshift functions like ϕ(r) = log 1 + r 0 r and 1 r do not satisfy the asymptotic criteria of the shape function and hence they are not suitable for this study when working under the influence of global monopole charge (see Subsubsections IV A 2-IV A 3 and IV B 2-IV B 3).Thus, we discuss the below conclusion on the basis of the constant redshift function only.   5.These observations underscore the consistency and applicability of the calculated shape functions within these theoretical frameworks, indicating the viability of our solutions within the context of wormhole physics.
6. Embedding diagrams for the calculated shape functions given by Eqs.(30)  9. Notably, in the isotropic EoS case, the calculated shape function fails to satisfy the asymptotic flatness condition, i.e., the ratio b(r) r does not tend to 0 as r → ∞.Consequently, obtaining asymptotically flat wormhole solutions for isotropic pressure under a constant redshift function for the f (Q) model given in Eq. ( 24) becomes challenging.

Table I contains details regarding the energy con-
ditions for different EoS cases with constant redshift function.It is assumed that the wormhole throat is situated at r 0 = 1.
11.In Section V, we utilized a Tolman-Oppenheimer-Volkov tool to assess the stability of the wormhole.A notable observation from graphs 18-19 is that hydrostatic and anisotropic forces exhibit similar patterns despite their opposite directions.This symmetrical balance between the forces suggests the stability of the wormhole solutions obtained in our study.12.In addition, we maintain other parameters at constant values, including α = −2, β = 0.02, and r 0 = 1.
In conclusion, our derived solutions for the barotropic and anisotropic EoS cases are deemed physically viable within the f (Q) gravity framework with a very small contribution of global monopole charge.As a potential avenue for future research, extending this analysis to explore other modified theories of gravity could offer valuable insights and contribute to a deeper understanding of gravitational dynamics in diverse astrophysical scenarios.Moreover, exploring wormhole solutions within the framework of f (Q) gravity, considering non-constant redshift functions and incorporating various matter sources, would be an intriguing avenue for future research.

8 FIG. 1 .FIG. 2 .
FIG.1.The contour plot illustrates the barotropic EoS, depicting changes in the asymptotic flatness condition for ω (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

FIG. 5 . 18 FIG. 6 .
FIG.5.The contour plot illustrates the barotropic EoS, depicting changes in the radial pressure for NEC for ω (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

16 FIG. 7 .
FIG.7.The contour plot illustrates the barotropic EoS, depicting changes in the SEC for ω (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

4 FIG. 8 .
FIG.8.The contour plot illustrates the barotropic EoS, depicting changes in the asymptotic flatness condition for ω (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r' with the redshift function ϕ(r) = log 1 + r 0 r .Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

FIG. 9 .FIG. 10 .
FIG.9.The contour plot illustrates the barotropic EoS, depicting changes in the asymptotic flatness condition for ω (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r' with the redshift function ϕ(r) = 1 r .Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

FIG. 11 .FIG. 12 .
FIG. 11.The plot illustrates the anisotropic EoS, depicting changes in the flaring-out condition for various values of n (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

FIG. 13 .FIG. 14 .
FIG.13.The plot illustrates the anisotropic EoS, depicting changes in the energy density for various values of n (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

FIG. 15 .FIG. 16 .
FIG.15.The plot illustrates the anisotropic EoS, depicting changes in the NEC for tangential pressure for various values of n (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

FIG. 17 .
FIG. 17.The contour plot illustrates the anisotropic EoS, depicting changes in the asymptotic flatness condition for n (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r' with the redshift function ϕ(r) = log 1 + r 0 r .Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

3 .
The shape function derived under the linear barotropic EoS case demonstrates adherence to the asymptotic flatness condition b(r) r → 0 as r → ∞ and the flaring-out condition b ′ (r 0 ) < 1 for the EoS parameter ω and Global Monopole parameter η, as depicted through contour plots 1-2.

40 FIG. 18 .
FIG.18.The contour plot illustrates the barotropic EoS (p r = ωρ), depicting changes in the hydrostatic and anisotropic forces for ω (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1. .

FIG. 19 .
FIG.19.The plot illustrates the anisotropic EoS (p t = np r ), depicting changes in the hydrostatic and anisotropic forces for various values of n (depicted on the left) and η (depicted on the right) concerning the radial coordinate 'r'.Additionally, the plot maintains constant values for other parameters, such as α = −2, β = 0.02, and r 0 = 1.

TABLE I .
Barotropic EoS (p r = ωρ)Anisotropic EoS (p t = np r ) Isotropic EoS (p t = p r ) Summary of energy conditions for results of linear EoS cases under constant redshift function and (42) under both EoS cases are presented in Figs. 3 and 12. 7. Analysis of the energy density, as depicted in Figs. 4 and 13 reveal a positively decreasing behavior throughout space-time under both linear EoS cases.8. Examination of Figures 5-6 and 14-15 indicates violation of the NEC for radial pressure while it is satisfied for tangential pressure.Consequently, both the NEC and WEC are violated overall, suggesting the potential presence of exotic matter at the wormhole throat to keep the throat open.Additionally, it is noteworthy that the SEC (Refer Figs.7 and 16) is fulfilled for barotropic and anisotropic EoS cases.