Wormhole in the Milky Way galaxy with global monopole charge

Wormholes are tunnels or short-cuts in space-time, and their existence is very important for human civilization to express the vastness of space and time. So, it is necessary to analyze our own Milky Way galaxy if it can harbour any wormhole. This work is dedicated to the existence of wormhole geometry(at least theoretically) in the bulge and halo of the Milky Way Galaxy. The structure and existence of wormholes are verified in both the bulge and the halo region of the Milky Way galaxy (MWG). Different dark matter profiles like pseudo-isothermal, NFW and Universal Rotational Curve (URC) are analyzed to harbour these cosmic tunnels. Three kinds of redshift functions are used for each dark matter profile with the global monopole charge to cover all the possibilities of MWG supporting wormhole geometry.


Introduction
In the year 1915, Albert Einstein has published his most important work by equating gravity with geometry. In the later years, it was experimentally/observationally verified [1] and is being used by researchers all over the globe as a stepping stone to understand the Universe in a better way. This theory gave rise to bizarre events like gravitational lensing, black holes, gravitational waves, wormholes etc. Among them, wormhole geometry is another solution of Einstein's Field Equations (EFE) apart from black hole solution that connects two distinct space-time of the same or different universes widely separated. The first wormhole solution, the Einstein-Rosen bridge, was proposed by Albert Einstein and Nathan Rosen [2]. However, it was not paid much attention because the wormhole was non-traversable and considered only a theoretical solution. In 1973, Ellis showed a wormhole having spherical symmetry with the help of a ghost massless scalar field [3]. Later Morris and Thorne developed a theoreta e-mail: das.priyam888@gmail.com (corresponding author) b e-mail: kalam@associates.iucaa.in ical traversable wormhole which, unlike the Einstein-Rosen bridge, does not contain any singularity or event horizons [4]. Matt Visser and Maccone have shown that inter universe travel is possible with the help of traversable stable wormholes [5,6].
Recent discoveries about the Universe has revealed that all the planets, stars and visible matter comprise around 5% of the total matter and energy, where the rest is dominated by dark matter(DM) and dark energy, among which 27 % is dark matter. Although this DM does not interact with photons like normal baryonic matter, but only has gravitational effects. Astronomer Zwicky in the year 1933, first determined the existence of DM in a distant galaxy cluster [7,8]. The presence of DM in the Milky Way galaxy was first speculated by Oort [16,17] and later confirmed observationally by Diemand and Springel [9,10]. In the later years, further experiments regarding the rotational curve of outer regions of galaxies and gravitational microlensing were conducted to understand the DM distribution [11]. Dark matter is an important entity for the fabrication of wormholes. The DM density profile proposed by researchers over the year plays a significant role in theoretically developing a fully stable and traversable wormhole.
This research aims to investigate the existence of wormhole geometry only in the bulge and the halo region of MWG. The bulge is the core of the Milky Way galaxy formed with compact stars, also seen in other spiral galaxies. The core of MWG contains a supermassive black hole around which a large numbers of old star clusters are there. This galactic centre is known as the bulge, and it roughly extends to 0.5 kpc [12]. Studies conducted by Mamon et al. [13] of the mass and shape of the dark galactic halo revealed that the core of the galaxy is being dominated by baryonic matter rather than dark matter. In 2011, studies of Martinez-Valpuesta showed that the galactic centre must have separate flat long bar components that are twisted according to the barred bulge of the galaxy [14]. The halo of the Milky Way Galaxy (MWG) comprises old stars and globular star clusters. Although at the disk of the MWG star formations take place, the halo has less gas to offer. Thus no star formation takes place [15]. As DM may be an important entity for the formation of wormholes, the halo of MWG is analyzed along with the Universal Rotational Curve. Studies revealed DM being dominant in the halo region of the MWG, thus it may be the best candidate to harbour wormhole structures.
Global monopole charge is a topological defect that came into existence due to the phase transition mechanism in the early universe [18]. It is the result due to global O(3) symmetry breaking of the self-coupling triplet of scalar fields φ a that belongs to objects like cosmic strings, domain walls etc. It is an important topic in modern physics and has caught the attention of researchers to investigate its space-time geometry [19][20][21][22]. Studies have shown that this kind of topological defect contains negative gravitational potential [23]. The global monopole charge has also been used in cosmology [24][25][26]. This paper deals with a new kind of Morris Thorne wormhole solution with a global monopole charge in 3 + 1 gravity minimally coupled to a triplet of scalar fields. Sarkar et al. and Jusufi analyzed the stability of wormholes with global monopole charge with MacMillan DM profile in the bulge of MWG [27,28]. Whereas we have used pseudo-isothermal, Novarro-Frank-White(NFW) and Universal Rotational Curve (URC) DM profile in the bulge and halo region of the MWG to check the existence of wormholes under this DM profiles along with global monopole charge.

Einstein's field equations with global monopole charge and wormhole metric
The formulation of the Einstein field equations with global monopole charge obtained from 3 + 1 dimensional action without considering the cosmological constant minimally coupled triplet scalar field (c = G = 1) is given by where k = 8π . Lagrangian density of a self-coupling scalar triplet φ a is given by where a = 1, 2, 3 and η, λ are scale of gauge-symmetry breaking, self interaction term respectively. Field configuration described by monopole is described as where the variable x a is defined as (rsinθ cosφ, rsinθ sinφ, rcosθ) so that a x a x 2 = r 2 .
Since our main interest lies in analysing traversable wormholes in the MWG, we consider a Morris-Thorne wormhole in four-dimensional spacetime. The wormhole metric in Schwarzschild coordinate is denoted by [29] where φ(r ) and b(r ) are the redshift and shape functions, respectively. Determination of these functions is very important to examine the wormhole geometry. The redshift function must be finite since the Morris-Thorne wormhole is devoid of any singularity and event horizons. In the equation, the minimum value of r is r th , also known as the throat radius of the wormhole. The shape function must obey the relation b(r th ) = r th . The redshift function and shape function must be well defined beyond r > r th for the wormhole to be traversable. The flare out condition must also be satisfied by the shape function b(r), given by, for r > r th in which b (r th ) < 1. By using the field configuration, the Lagrangian density can be written in terms of f (r ) as and the Euler-Lagrangian equation for the field f (r ) is From Eq. (2), the energy momentum tensor can be obtained as Now, from Eq. (8) all the four components of Energy-Momentum tensor are as follows: As seen in Eq. (7), it is very complex to find an exact analytical solution. So the approximation of f (r ) → 1 outside the wormhole is enough to simplify it to get the solution. Therefore, the reduced energy-momentum components are written as Now, we know that the Einstein field equation is where T μν is the sum of the energy-momentum tensor of anisotropic fluid part and matter field part. Therefore, it can be written as The components of the energy-momentum part of the anisotropic fluid are Solving the Einstein field equation for Morris-Thorne wormhole metric, the components we get are, Now, combining the Einstein field equation for the wormhole along with the Global Monopole Charge, the solutions can be written as where ρ is the energy density, P r is the radial pressure and P t is the transverse pressure. We have used c = G = 1 and f = f(r).

Different dark matter profiles
Sarkar et al. [27] has recently able to shown the wormhole geometry in the bulge of MWG using the Mc Millan DM profile. We have considered three possible DM profiles to generalize and maximize the probability of wormhole geometry supported by our galaxy-the NFW, pseudo-isothermal and URC profile. The NFW profile as studied by Novarro et al. [30]. represented the dark matter profile which is defined by where, ρ 0 and r 0 is the core density and the core radius respectively. This profile along with the Global Monopole Charge (GMC) will be used for both the bulge and the halo region of the Milky Way galaxy to investigate its effect in wormhole geometry. Now, we take another DM profile known as the Pseudo-Isothermal, represented by, Although this profile is not used in the determination of wormhole in the galactic halo because the wormhole does not satisfy the flare out condition and the necessary violation needed for the NEC. The URC or Burkert DM density profile along with the Global Monopole Charge is also considered for the investigation of wormhole geometry in both the bulge and halo of MWG. The density relation is given by, where r 0 and ρ 0 are the core radius and density respectively.

Null energy condition and redshift function
In analyzing the formation of traversable wormhole geometry, the Null energy condition (NEC) is of utmost importance. The violation of NEC must be satisfied by the anisotropic matter content in the wormhole [31][32][33]. The NEC can be expressed in terms of general relativity as T μν v μ v ν ≥ 0, where n μ is a null vector [34] represented by ρ(r ) + P r (r ) ≥ 0. Thus for the violation of NEC v μ ≤ 0, and verification of this condition, one must know the P r (r ), and from Eq. (20), it is evident that the knowledge of the redshift function is necessary to compute the radial pressure and thus determine the violation of NEC. We have used three redshift functions for more generalization of our model. For the redshift function is called the tidal force [35]. Since the function is finite, it does not have any event horizons or singularity. By taking φ = 0, we have shown that this redshift function satisfies the violation of NEC along with Global Monopole Charge (η) ranges from 0 ≤ η ≤ 0.006. Thus the DM content in the bulge and halo of the Milky Way galaxy is well enough to sustain wormhole geometry in the two regions. The second redshift function is taken as where τ is an non zero constant and remains finite for r > 0. This kind of redshift function also does not have an event horizon. Moreover, substituting the redshift function in Eq. (20), we have obtained total NEC violation for all the dark matter profiles in both the bulge and halo of the Milky Way galaxy.
The final redshift function is obtained from the flat rotational curve of galaxies. In the studies of Chandrasekhar and Landau [36,37] the tangential velocity for stable circular geodesic motion can be represented by (v Φ ) 2 = r φ (r ). Rahaman et al. [38] proposed a theoretical function to describe the flat rotational curve which is where r is the distance in kpc and v Φ is the velocity in Km/s. Now by using the above equation, the redshift function can be easily obtained as where E[x,y] is exponential integral function and D is the integration constant and the value of the parameters to be k 1 = k 2 = 1.α = 0.8,β = 0.001. Similar to the previous redshift function, we used it to find the radial pressure and it is seen to violate the NEC condition for all the dark matter profiles we have considered. This will prove that the DM profiles considered in our articles are appropriate to sustain wormhole geometry in both the bulge and halo of the Milky Way Galaxy. Moreover, the Global Monopole Charge is studied in the recent paper of Rahaman. et al. shows that the tendency of the NEC violation decreases with the increase of η.

Embedding surface of the wormhole
The embedding surface is defined by z(r ) function and is described by the following equation The above function is not a conversing function at b(r th ) = r th instead it becomes vertical at the throat of the wormhole.
In this paper, the embedding function is plotted in a 3D graph to get the full visualization of the wormhole for each case.

Existence of wormhole geometries in the bulge of Milky Way galaxy
5.1 Wormhole using NFW dark matter profile in the galactic bulge Core radius and core density is taken as ρ 0 = 0.004 kpc −2 and r 0 = 2 kpc. We can get the throat radius, r th = 1.7 kpc from the linear graph b − r with r where b − r cuts the r axis. The shape function for the wormhole geometry is obtained as shown in Fig. 1a.
where A is int. constant and is determined by using the boundary condition b(r th ) = r th . In Figs. 2a and 3a, the flare out condition is satisfied by taking 0 ≤ η ≤ 0.006 as done by S. Sarkar et. al. in the article [27]. We can see b r < 1 and b−r.b b 2 < 1 are satisfied. In the above graphs, if the GMC η is taken to be 0, it will represent the Einstein gravity without any topological defect. Now using Eq. (25) in Eq. (20) to find the P r of the matter of wormhole for investigating the NEC is given by Fig. 1 The shape function b(r) with respect to r (kpc) and η is plotted by using ρ 0 = 0.004 kpc −2 , r 0 = 2 kpc for a NFW profile : r th = 1.7 kpc b Pseudo-Isothermal profile : r th = 1 kpc (c)URC profile : r th = 1.1 kpc Fig. 2 The ratio b(r ) r with respect to r (kpc) and η is plotted by using ρ 0 = 0.004 kpc −2 , r 0 = 2 kpc for a NFW profile : r th = 1.7 kpc, b Pseudo-Isothermal profile : r th = 1 kpc, c URC profile : r th = 1.1 kpc with respect to r (kpc)and η is plotted by using ρ 0 = 0.004 kpc −2 , r 0 = 2 kpc for a NFW profile : r th = 1.7 kpc, b Pseudo-Isothermal profile : r th = 1 kpc, c URC profile : r th = 1.1 kpc Fig. 4 NEC with respect to r (kpc) and η is plotted by using ρ 0 = 0.004 kpc −2 , r 0 = 2 kpc, r th = 1.7 kpc, k 1 = k 2 = 1, α = 0.8, β = 0.001 and τ = 2.5 for NFW profile Thus, the NEC is violated as can be confirmed by the plot in and the NEC is also violated as can be inferred from Fig. 4b. The final expression for the radial pressure is obtained by using Eq. (28) in Eq. (20).
Similarly the violation of null energy is shown in Fig. 4c.  Figure 7a shows the embedding diagram of the wormhole and Fig. 8a shows the full visualization of the wormhole structure.

Wormhole using pseudo-isothermal dark matter profile in the galactic bulge
Core radius and core density is taken as ρ 0 = 0.004 kpc −2 and r 0 = 2 kpc. We can get the throat radius, r th = 1 kpc from the linear graph b − r with r where b − r cuts the r axis. Figure 1b shows the variation of the shape function for the wormhole.
The expression for the shape function as obtained for the pseudo-isothermal DM profile coupled with the global monopole charge is given by is being satisfied. And Fig. 2b it is shown that b(r ) r < 1. Moving over to the most important condition for the sustainability of the wormhole structure in any location of the universe, i.e. the NEC violation condition using the three redshift functions obtained in Sect. 4. Similarly, the radial pressure for all the three redshift functions is computed, and then the violation of NEC is checked. The P r (r ) for the function φ = 0 is given by Figure 5a represents the violation of the NEC for the radial pressure computed in Eq. (36). For the second redshift function φ = τ r , the new radial pressure P r (r ) is computed as P r (r ) = − 1 8πr 4 r 2 − 8η 2 πr 2 + 8.37758(r − 2τ ) 0.115 −0.1193r + η 2 (r 3 − 1) + 3r 2 0 rρ 0 −3r 3 0 ρ 0 T an −1 r r 0 (36) This radial pressure when added with the DM density profile to check for the violation of NEC shows clear result of violating the NEC condition as evident from Fig. 5b. Finally using the redshift function in Eq. (28) which has been obtained from the flat rotational curve, the radial pressure P r (r ) is expressed as, Using the above equation the NEC is satisfied for the final redshift function too. It is plotted in Fig. 5c. Figures 7b and 8b shows the embedding diagram and the full visualization of the wormhole in the galactic bulge respectively.

Wormhole using URC dark matter profile in the galactic bulge
Finally, the URC or Burkert dark matter profile is analysed to support the wormhole structure in the bulge of the Milky Way galaxy. The EFE are solved by coupling them with the Global Monopole Charge. The core radius and core density is taken as the same i.e. ρ 0 = 0.004 kpc −2 and r 0 = 2 kpc. We can get the throat radius, r th = 1.1 kpc from the linear graph b − r with r where b − r cuts the r axis. Figure 1c shows the variation of the shape function with the global monopole charge and distance.
The expression of shape function of the wormhole obtained from the URC dark matter profile is given as where, A is the integral constant that is derived by imposing the boundary condition of b(r th ) = r th . The flare out condition is satisfied as evident from the Figs. 3c and 2c shows b(r ) r ≤ 1. The value of η is taken to be from 0 ≤ η ≤ 0.006 Now, the other important condition, i.e. the NEC, is being checked for violation. To check the violation of NEC, the radial pressure P r (r ) is needed to be computed with the help of the three redshift functions obtained in Eqs. (25,26,28).  6 NEC with respect to r (kpc) and η is plotted by using ρ 0 = 0.004 kpc −2 , r 0 = 2 kpc, r th = 1 kpc, k 1 = k 2 = 1, α = 0.8, β = 0.001 and τ = 2.5 for URC profile Expression for P r (r ) obtained for φ = 0 is The violation of NEC can be observed using this radial function plotted in Fig. 6a. The P r (r ) expression as obtained from the second redshift function φ = τ r is represented as The NEC violation for this particular redshift function is demonstrated in the Fig. 6b. The P r (r ) expression as obtained from the redshift function, and φ obtained from flat rotational curve is represented as The NEC violation for this particular redshift function is demonstrated in the Fig. 6c. The embedding diagram and the full visualization of the wormhole in the galactic bulge for URC is shown in Figs. 7c and 8c.  In this section, the NFW dark matter density profile for the galactic halo is used for investigating the existence of wormhole geometry in the outer halo of the Milky Way galaxy. The core density ρ 0 = 0.05 kpc −2 and core radius r 0 = 1kpc are as per article by Rahaman et al. [38]. We can get the throat radius, r th = 1.7kpc from the linear graph b − r with r where b − r cuts the r axis and b(r)-r ≤ 0 for r > r th .
The shape function is same as that of central bulge with varying parameters. It is denoted by where S is the integration constant. The variation of shape function with the distance and the Global Monopole Charge is demonstrated in Fig. 9a. Now, in the Figs. 10a and 11a, the flare out conditions of the wormhole are satisfied. Now the violation of NEC is checked by using the expression of the P r (r ) obtained in Eqs. (39,40,41) and plotting the NEC in Fig. 12.   Fig. 9 The shape function b(r) with respect to r (kpc) and η is plotted by using ρ 0 = 0.05 kpc −2 , r 0 = 1 kpc with 0 ≤ η ≤ 0.006 for a NFW profile: r th = 1.7 kpc, b URC profile: r th = 1 kpc r with respect to r (kpc) and η is plotted by using ρ 0 = 0.05 kpc −2 , r 0 = 1 kpc with 0 ≤ η ≤ 0.006 for a NFW profile: r th = 1.7 kpc, b URC profile: r th = 1 kpc with respect to r (kpc)and η is plotted by using ρ 0 = 0.05 kpc −2 , r 0 = 1 kpc with 0 ≤ η ≤ 0.006 for a NFW profile : r th = 1.7 kpc, b URC profile : r th = 1 kpc Fig. 12 NEC with respect to r (kpc) and η is plotted by using ρ 0 = 0.05 kpc −2 , r 0 = 1 kpc, r th = 1.7 kpc, k 1 = k 2 = 1, α = 0.8, β = 0.001 and τ = 2.5 for NFW profile Fig. 13 NEC with respect to r (kpc) and η is plotted by using ρ 0 = 0.05 kpc −2 , r 0 = 1 kpc, r th = 1.7 kpc, k 1 = k 2 = 1, α = 0.8, β = 0.001 and τ = 2.5 for URC profile Thus from the above figures, it is evident that the NFW dark matter profile coupled with the Global Monopole Charge has the suitable fuel for the supporting and existing wormhole structures in the halo of the Milky Way galaxy. Fig. 14a shows the embedding graph for the Halo of Milky Way for NFW dark matter profile. Figure 15a shows full visualization of the wormhole.
6.2 Existence of wormhole in the galactic halo using URC dark matter profile The final dark matter profile that we are considering for our research for the existence of wormholes structure in our Milky Way galaxy is the URC profile. The core radius and core density is taken as ρ 0 = 0.05kpc −2 and r 0 = 1kpc. We can get the throat radius, r th = 1kpc from the linear graph b − r with r where b − r cuts the r axis.
The expression for the shape function is same as that of Eq. (38) with only the external parameters has been changed. The graphical representation of the shape function is given in Fig. 9b. Now, the shape function is tested for the flare out conditions and is verified from Figs. 10b, 11b, that it obeys the flare out condition.
Finally the NEC is checked for violation by plotting ρ + p against r and η. They are demonstrated in Fig. 13.
Embedding function for this particular wormhole is constructed as shown in Fig. 14b and the full visualization of the Wormhole is shown in Fig. 15b.

Wormhole in the expanding universe
It is well known now that our Universe is expanding with an acceleration due to the effect of dark energy. The main factor determining the behaviour of wormholes in an expanding Universe is that the violation of NEC i.e. (ρ + P) < 0. This indicates that the energy density increases with expanding cosmic time to an infinite value. Pedro has performed the detailed calculations. Gonzales Diaz and F. Pedro [40] has investigated the evolution of the circumference, throat and the embedding surface and has shown that the wormhole will inflate with time as the Universe expands. Inspiring with that we have also investigated similar circumstances and try to check how the global monopole charge affects the expansion of the wormhole's physical attributes. The equation describes the de-Sitter scale factor for the expanding Universe, [40] where χ = √ Λ/3 and Λ = 2.036 × 10 −35 s −2 is the cosmological constant. The state factor is given by p = ωρ and ω < −1. Thus, the Morris-Thorne wormhole metric in an expanding Universe changes to Throat radius of the wormhole is denoted by r th , so the proper circumference in the accelerating Universe is given by This explicitly shows that the expansion of the circumference is only depends on the scale factor of the expanding Universe with a quintessence field, and the GMC has no effect. The variation of the proper radial distance between two arbitrary large points has been computed. The de-Sitter scale factor g(t) and the GMC depends on time. The equation is described as The below diagram (Fig. 16) shows that the radial distance increases with time and the topological defect. Similarly, since the embedded function's evaluation is same as radial distance, the embedding function also increases with the expansion of the Universe. Moreover, after evaluating the wormholes for an expanding Universe, if we look at the first case where the wormhole has zero redshift function with no tidal force, the new Eq. (44) represents the de-Sitter metric for the expanding Universe.

Conclusion
This article investigates the existence of a wormhole in the bulge and halo of the Milky Way galaxy. Three different Dark matter profiles are taken for the investigation. Previous studies have shown the importance of Dark matter in constructing wormholes since it violates the null energy condition. Although NEC violation has not been artificially demonstrated in the laboratory, the DM is theorized to violate the NEC. The three Dark matter profiles taken are pseudo-Isothermal, Burkert and NFW. For generalization, the topological defect from the early universe, i.e. the global monopole charge, is considered to broaden the possibility of a wormhole in the Milky Way galaxy. Sarkar [39] has also shown the supporting of wormhole geometry in the halo of different massive spiral and dwarf galaxies. Thus, our study in this article further enforces the works on supporting wormhole geometry in our galaxy's halo and central bulge.
This article shows that the pseudo isothermal, Burkert and NFW DM profiles coupled with the global monopole charge accurately supports the wormhole structure in the central bulge region. The shape function gives the wormhole's embedded function and, upon its rotation along the z-axis, gives the full visualization of the wormhole structure. Interestingly, the throat radius for all the models is obtained between 1 to 2 kpc. The flare out conditions have also been satisfied by the above dark matter profiles, another important criterion for wormhole geometry. The NEC has been checked using three redshift functions for more generalization: tidal force, rotational redshift function, and redshift function obtained from the Universal rotational curve. All the redshift functions corresponding to e φ are finite. NEC has been checked individually for the redshift functions and is duly satisfied. However, the presence of a Global monopole charge is seen to minimize the violation of NEC, as shown in Fig. 16. Figure 8 shows the full visualization of the wormhole.
Contrary to the supporting of wormhole structure by different DM profiles coupled with the GMC in the bulge of the Milky Way galaxy, only two among these three in our study have been found to support wormhole geometry in the galactic halo region. The pseudo isothermal DM profile failed to violate the NEC when coupled with the global monopole charge. The existence of wormhole geometry has been successfully investigated for the NFW and URC dark matter profile. Figures 10 and 11 shows the satisfaction of the flare out conditions and Figs. 12 and 13 shows the violation of NEC by the three redshift functions individually. The throat radius of the wormholes is determined to be around 1 to 2 kpc. Finally, the embedding graph and the full visualization of the wormhole structures have been constructed in Figs. 14 and 15.
The wormholes are deemed traversable for a human being if the tidal gravitational forces acting are very small. Arguments have been made regarding some black hole solutions mimicking wormholes. Our study is intended to inspire researchers and scientists to investigate the existence of wormholes in the bulge and the halo of our galaxy, as it is entirely theoretically possible for them to exist. Scientist can search for several clues within our galaxy to confirm the presence of wormhole structure, which includes, first pinpointing areas of the galaxy that contains enough dark matter for NEC violation then searching for the iron K α line profiles that are emitted from accretion disks formed due to traversable wormholes [94], gravitational microlensing of background universe due to the presence of traversable wormholes, a gravitational anomaly in trajectories of nearby wormhole objects due to the gravitational effect of objects circulating near the wormhole in another space [95] and formulation of gravitational waves that will be generated upon objects traversing a wormhole. Since the analytical study about the existence of wormholes in our universe in this paper papers proves the existence of wormhole geometry with a wider generalized case.