Topological Nambu monopole in two Higgs doublet models

We show that a topological Nambu monopole exists as a regular solution for a large range of parameters in two Higgs doublet models, contrary to the standard model admitting only non-topological Nambu monopoles. We analyze a Higgs potential with a global $U(1)$ symmetry and a discrete symmetry $\mathbb{Z}_2$. The monopole is attached by two topological $Z$ strings ($Z$ flux tubes) from both sides. Despite of a trivial second homotopy group, the discrete symmetry $\mathbb{Z}_2$ together with a non-trivial first homotopy group for $Z$ strings topologically ensures the topological stability. After analytically constructing an asymptotic form of such a configuration, we explicitly construct a solution of the equation of motion based on a 3D numerical simulation, in which magnetic fluxes spherically emanating from the monopole at large distances are deformed in the vicinity of the monopole. Since the monopoles are expected to be abundant in the present universe, they might be observed in the current monopole searches.


I. INTRODUCTION
Magnetic monopoles have attracted great interests from both experimental and theoretical physicists since the seminal work by Dirac [1], which was motivated to improve the asymmetry between electric and magnetic charges in the Maxwell equations, providing an explanation for the electric charge quantization.In field theoretical models, they have been theoretically realized as a regular solution [2,3] and have played crucial roles to study non-perturbative aspects of (non-)supersymmetric field theories [4][5][6].However, except for condensed-matter analogues [7,8], such monopoles have never been found in reality; for instance, such magnetic monopoles are predicted by all grand unified theories (GUTs) [9][10][11][12], and their search have been extensively conducted.Nevertheless, no GUT monopoles have been found probably because the cosmological inflation diluted them.To avoid dilution, monopoles should be produced below the inflation scale such as the electroweak scale.In fact, a monopole configuration in the Standard Model (SM) was first considered by Nambu [13].However, it is unstable because it must be attached from one direction by a string and the tension of the string pulls the monopole to infinity, although they were suggested to produce primordial magnetic fields before disappearance [14,15].
In this paper, we predict an existence of a topologically stable 1 and magnetic monopole without singularity in a quite simple extension of the SM, two Higgs doublet model (2HDM), in which one more Higgs doublet is added to the SM (for reviews, see, e.g., Refs.[19,20]).They might be observed in the current monopole searches since they are expected to be abundant in the present universe.Apart from cosmological production, our monopoles might also be produced and detected by the MoEDAL experiment [21] at LHC with masses of the order of a few TeV.Their discovery would not only be a realization of Dirac's hypothesis but also yield a solid evidence of new physics beyond the SM, since there are no stable and regular magnetic monopoles in the SM.
The reason of the non-existence of stable monopoles in the SM is its trivial topology, that is, the vacuum manifold is S 3 having a trivial second homotopy group π 2 , as well as trivial π 0 for domain walls and π 1 for cosmic strings.Nevertheless, (non-topological) electroweak Z-strings (or magnetic Z-fluxes) [22][23][24][25][26] have been studied, but they were shown to be unstable in realistic parameter region [27,28].Nambu monopoles are the end points of these strings [13].
2HDM is one of the most popular extensions of the SM with a potential to solve problems that are unanswered by the SM.It has four additional scalar degrees of freedom in addition to 125 GeV Higgs boson (h), which are charged Higgs bosons (H ± ), CP-even Higgs boson (H) and CP-odd Higgs boson (A).These additional scalars can be directly produced at LHC, though there is no signal so far today, placing lower bounds on masses of those additional scalar bosons.Those lower bounds highly depend on parameter choices of 2HDM as well as how SM fermions couple to the two doublets.For more detailed phenomenological studies, see, e.g., Refs.[29][30][31][32][33] and references therein.Moreover, 2HDM has a much richer vacuum structure than the SM, therefore allowing a variety of topologically stable solitons, in addition to non-topological solitons [34][35][36][37][38][39][40][41][42] analogous to the SM; domain walls [43][44][45][46], membranes [47,48], and cosmic strings such as topological Z strings [45,46,49,50] (see also Ref. [51]).However, magnetic monopoles were not examined because of a trivial second homotopy group π 2 as in the SM.Instead, the stability of our monopole is topologically protected by a combination of the following two symmetries of the Lagrangian; One is a global U (1) symmetry that ensures the stability of the topological Z strings.The other is a discrete symmetry Z 2 exchanging the topological strings.Consequently, our monopole is attached by two topological Z strings on both sides, where a Z flux is confined on each string and the string tensions pulling the monopole are balanced due to the Z 2 symmetry.We explicitly construct such a solution of the equation of motion based on a 3D full numerical simulation, in which magnetic fluxes emanate from the monopole.

II. THE MODEL
We introduce two SU (2) doublets, Φ 1 and Φ 2 , both with the hypercharge Y = 1.The Lagrangian which describes the electroweak and the Higgs sectors is written as Here, Y µν and W a µν describe field strength tensors of the hypercharge and the weak gauge interactions with µ (ν) and a being Lorentz and weak iso-spin indices, respectively.D µ represents the covariant derivative acting on the Higgs fields, and the index i runs i = 1, 2. V (Φ 1 , Φ 2 ) is the potential for the two Higgs doublets.In this paper, we assume that both Higgs fields develop real vacuum expectation values (VEVs) as Φ 1 = (0, v 1 ) T , Φ 2 = (0, v 2 ) T .Then the electroweak scale, v EW ( 246 GeV), can be expressed by these VEVs as v2 EW = 2v 2 1 + 2v 2 2 .For later use, we introduce the Higgs fields in a two-by-two matrix form [52], H, defined by , and therefore the covariant derivative on H is expressed as , and the potential can be written by using H as follows: where |H| 2 ≡ H † H and we have imposed a (softly-broken) Z 2 symmetry, H → Hσ 3 (or Φ 1 → +Φ 1 , Φ 2 → −Φ 2 ), in order to suppress Higgs-mediated flavor-changing neutral current processes.In this paper, to make the discussion simpler, we take the five parameters m 2 , m 3 , α 3 , α 4 , α 5 to 0. The potential thus reduces to the following simple form: . The custodial transformation acting on the matrix H [52,53] is defined as the following global SU (2) transformation: In addition, the SU (2) W gauge field also transforms as an adjoint representation.The potential V (H) has a symmetry under this transformation, which we call as the custodial symmetry.However, the presence of the U (1) Y gauge field explicitly breaks the symmetry down to U (1) EM (Z 2 ) C , where (Z 2 ) C transforms H to iσ 1 H(iσ 1 ) † and " " denotes the semidirect product because (Z 2 ) C acts on U (1) EM .Note that the symmetry (Z 2 ) C follows from m 2 = α 4 = 0, and yields tan In addition, since m 3 = α 5 = 0, the Lagrangian (not only the potential) is invariant under a global U (1) a transformation, which rotates the relative phase of the two doublets: After H gets a VEV, this U (1) a symmetry is spontaneously broken and the corresponding Nambu-Goldstone boson appears.

III. ELECTROWEAK STRINGS
In Refs.[45,46,49,50], it is pointed out that, unlike in the SM case, 2HDM allows topologically stable strings to exist thanks to the global U (1) a symmetry.First, consider topological strings with the Z flux (topological Z strings).
w a b S 9 1 m r W O l t u l 8 e P o + t u / q g r N E g d f q j 8 9 S 5 S x 6 H s V 5 N 3 2 m f Y t 9 I 6 + 1 j h 5 W V 9 a S 3 i z 7 J w 9 k / 8 z 1 m L 3 d A O z 9 q p f 5 v j a K c L 0 A a m f z 9 0 N N t w a b S 9 1 m r W O l t u l 8 e P o + t u / q g r N E g d f q j 8 9 S 5 S x 6 H s V 5 N 3 2 m f Y t 9 I 6 + 1 j h 5 W V 9 a S 3 i z 7 J w 9 k / 8 z 1 m L 3 d A O z 9 q p f 5 v j a K c L 0 A a m f z 9 0 N N t w a b S 9 1 m r W O l t u l 8 e P o + t u / q g r N E g d f q j 8 9 S 5 S x 6 H s V 5 N 3 2 m f Y t 9 I 6 + 1 j h 5 W V 9 a S 3 i z 7 J w 9 k / 8 z 1 m L 3 d A O z 9 q p f 5 v j a K c L 0 A a m f z 9 0 N N t l a t e x i t > ⇣ < l a t e x i t s h a 1 _ b a s e 6 4 = " m v q a t a p 8 S C Y o T H K 0 n p m 6 Y x T P y 7 I = " > A A A C a H i c h V F N L w N B G H 6 6 v q o + u j g g L q I h T s 2 s S I i T x M U R t d o E k d 0 1 m H S / s j t t 0 j b + g I s j 4 k Q i I n 6 G i z / g 0 J + A I 4 m L g 3 e 3 m w i C d z I z z z z z P u 8 8 M 2 P 6 t g g l Y 8 2 U 0 t b e 0 d m V 7 s 7 0 e p J 0 V V h p X R V q q S S j R D + B L K x A e 6 L o u k < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " m v q a t a p 8 S C Y o T H K 0 n p m 6 Y x T P y 7 I = " > A A A C a H i c h V F N L w N B G H 6 6 v q o + u j g g L q I h T s 2 s S I i T x M U R t d o E k d 0 1 m H S / s j t t 0 j b + g I s j 4 k Q i I n 6 G i z / g 0 J + A I 4 m L g 3 e 3 m w i C d z I z z z z z P u 8 8 M 2 P 6 t g g l Y 8 2 U 0 t b e 0 d m V 7 s 7 0 FIG.1: Tension of the one-parameter family of the strings made with U (ζ). ζ = 0, π/2, π correspond to the north pole, the equatorial point and the south pole on the moduli space S 2 , respectively.Although this value itself depends on the infra-red cut off, the difference between the maximum and minimum values does not.
There are two types of topological Z strings corresponding to which one of the two Higgs doublets is to be wound.To see that, let us take W ± µ = A µ = 0.Here we have defined as The solution called a (1, 0)-string3 is given by where and ϕ is the rotation angle around the z-axis.The boundary conditions imposed on the profile functions are On the other hand, the solution called a (0, 1)-string is given by . Both the (1, 0)-and (0, 1)-strings have winding number 1/2 for the global U (1) a and thus they are topological vortex strings.Note that both of them have logarithmically divergent tension due to the kinetic term of the Higgs field: which is a quarter of that for a global U (1) a integer vortex because of the half winding number for U (1) a [45].On the other hand, the contribution from the half winding in the gauge orbit is canceled by the gauge fields Z (0,1) i or Z (1,0) i and exponentially suppressed, and as a consequence, the Z flux is squeezed into the flux tube, which exponentially vanishes as a usual Abrikosov-Nielsen-Olesen vortex [54,55].In addition, the Z fluxes flowing inside them are ±2π cos θ W /g along the z-axis, respectively, which are half of that of a non-topological Z string in the SM.
Our potential has the custodial SU (2) C symmetry.When sin θ W = 0, this symmetry is the actual symmetry of the Lagrangian.However, the presence of the topological string solutions spontaneously break it down to U (1) C , and thus they have S 2 ( SU (2) C /U (1) C ) moduli.Each point on the S 2 moduli space corresponds to a physically different string solution having a common winding number 1/2 for the global U (1) a .We parametrize the S 2 moduli space by two parameters 0 ≤ ζ ≤ π, 0 ≤ χ < 2π, where ζ and χ correspond to the zenith and azimuth angles, respectively.We identify the (1,0)-string, Eq.( 4), as the one associated with the south pole of the S 2 moduli space, ζ = π.On the other hand, the (0,1)-string corresponds to the north pole, ζ = 0. String solutions on a generic point of the S 2 moduli space can be obtained by acting an SU (2) C transformation on the (1, 0)-string.In particular, we can construct a one-parameter family on S 2 connecting a (0, 1)-string (north pole) and a (1, 0)-string (south pole) using the σ 2 -axis rotation which corresponds to the one-parameter path on S 2 : 0 ≤ ζ ≤ π, χ = 0.This fact will be important in the construction of the magnetic monopole solution later.
We should note that in our case, sin θ W = 0 and the custodial symmetry is not exact.As a consequence, almost all the points of the S 2 moduli space are energetically lifted.As studied in Refs.[45,46], the two Z strings, (1, 0)-string and (0, 1)-string, are the most stable with degenerate among the topological strings.On the other hand, solutions on the equatorial points of the moduli space, which contain a W flux and are called as W strings, are the most unstable.Fig. 1 shows a plot of the string tension of the one-parameter family U (ζ)H (1,0) U (ζ) † .As we stated above, the (Z 2 ) C symmetry remains, which corresponds to flipping the upper and lower half spheres of the moduli space S 2 (or ζ → π−ζ) 4 .Thus, the presence of one of the Z strings spontaneously breaks the (Z 2 ) C symmetry.

IV. MAGNETIC MONOPOLE CONFIGURATION
Let us make a magnetic monopole configuration in 2HDM as follows.First, prepare a (0,1)-string and a (1,0)-string, cut both in half and connect them smoothly.In other words, make a configuration where the upper half (z > 0) is a (0, 1)-string, but the lower half (z < 0) is a (1, 0)-string.Such a smooth connection is realized by the one-parameter path on the moduli space made with U (ζ) (0 ≤ ζ ≤ π).Since (0,1)-string and (1,0)-string have the Z fluxes in the opposite directions, Z fluxes of 4π cos θ W /g in total flows from the connection.This is the same as the amount of a Z flux that a Nambu monopole discharges.Thus, the same argument for the conservation of a U (1) Y flux ensures that 4π sin θ W /g magnetic flux spreads from the connection, and thus it is a magnetic monopole.Interestingly, this can be regarded as one in which the non-topological Z string attached to a Nambu monopole is divided in two topological Z strings and these two fractions are pulled to the opposite directions to each other.
Let us discuss the stability of such a configuration.In the present case tan β = 1 and thus the tensions of the two Z strings are equal.Therefore, the monopole will not be pulled to one side, and such configuration is expected to be stable.This argument is intuitive and easy to understand, but we give a more rigorous discussion.The stability is ensured by the following two reasons.First, the whole of this configuration (both the strings and monopole) has a topological charge 1/2 for the global U (1) a symmetry in any cross sections with z = const..In other words, concentrating on the U (1) a charge, the configuration is an infinitely long global string.Therefore, it cannot be broken into pieces as usual global strings, and thus (1, 0)-and (0, 1)-strings must always be connected by a path on the lifted moduli space S 2 .Second, the string tensions on S 2 has a degenerate double-well structure (Fig. 1), where (1, 0)-and (0, 1)-strings are the most stable, and thus the above connection can be regarded as a topological kink interpolating between the two minima of the (Z 2 ) C symmetric potential.Therefore, the stability of the monopole is topologically ensured as in the case of a Z 2 kink.These arguments can be summarized in the statement that the stability is ensured by the U (1) a symmetry and the (Z 2 ) C symmetry.
We construct such a configuration concretely in the limit where we treat vortices as infinitesimally small (deltafunction like) topological defects.This limit corresponds to those far enough from vortices.The actual regular solution will be constructed numerically by the relaxation later.
The Higgs matrix H (0,1) of a (0,1)-string at large distance (ρ → ∞) is expressed as H (0,1) → H 0 ≡ v diag(1, e iϕ ).Let us introduce a function ζ(r, θ) and a unitary matrix U ( ζ) = e i 2 ζσ2 , and then acting a "local custodial transformation" on H 0 , we obtain Here, r and θ are the distance from the origin and the zenith angle from the z-axis, respectively, and ζ(r, θ) is assumed to be a monotonically increasing function with respect to θ and ζ(r, 0) = 0, ζ(r, π) = π.We should note that H ζ has a winding number 1/2 for U (1) a on any planes with z = const.and hence describes a global string along the z-axis.In fact, its kinetic energy behaves as tr|D i H ζ | 2 ∼ v 2 /ρ 2 for ρ → ∞ and hence the tension on z = const.is logarithmically divergent independently of ζ (or z), where the tension from the winding in the gauge orbit is canceled by ζ-dependent gauge fields and remains finite (See Appendix A).Furthermore, the string behaves as (0,1)-and (1,0)-strings around the positive and negative sides of the z-axis (θ ∼ 0, π), respectively, while it behaves locally as the W -string in a region in which ζ ∼ π/25 .The two Z-strings connect at the origin, and thus the Z flux flows upward and downward from the origin as we see later.In addition, the electromagnetic flux emanates from the connection and it behaves as a magnetic monopole.

V. Z FLUX AND MAGNETIC FLUX
The gauge fields induced from the Higgs field H ζ are determined to minimize the kinetic energy : After minimizing, we obtain asymptotic forms of the field strengths of the Z field and the electromagnetic field (See Appendix A for the derivation).The former is given by and F Z 23 = F Z 31 = 0. From Eq.( 6), it can be seen that Z flux 2π cos θ W /g flows on the z-axis toward z = ±∞ from the origin.
For the electromagnetic field, we obtain where t [ij] ≡ t ij − t ji for any tensor t.From Eq. ( 7), it is clear that there is a magnetic flux in a region where ∂ i ζ = 0, and that it is coming out of the connection at the origin.The total magnetic flux Φ B can be obtained by integrating the flux density B i ≡ − 1 2 ijk F EM jk on an infinitely large sphere covering the system and using the Gauss's theorem: where we have used ζ(r, 0) = 0 and ζ(r, π) = π in the last equality.Therefore, there is a magnetic monopole at the origin with a magnetic charge 4π sin θ W /g.This monopole carries a U (1) modulus coming from that of the localized W -string. Interestingly, the above analysis does not rely on details in the form of ζ(r, θ).The existence of magnetic monopole is determined only by the information of the end points that ζ = 0, π for θ = 0, π, respectively.This property is similar to that of topological kinks.On the other hand, the "shape" of the magnetic flux spreading from the monopole depends on the functional form of ζ(r, θ).By minimizing the magnetic energy, it is determined to be ζ = θ and we obtain the magnetic flux density B i = (sin θ W /g) x i /r 3 , which means that the magnetic flux is distributed spherically at large distances from the monopole.

VI. NUMERICAL SIMULATION
Here we show a stable regular magnetic monopole solution numerically constructed by relaxation.The procedure is as follows.We smear out the infinitely small defects (singularities) in the configuration discussed above using some profile functions with a typical scale ∼ v −1 .The smeared configuration is regular anywhere and approaches at large distances to the asymptotic form constructed above.Then, we evolve it by the relaxation until it sufficiently approaches the solution of the equation of motion.We have taken the parameters g, g , m 1 , α 1 , α 2 so that the physical parameters are given by sin 2 θ W = 0.23, m W = 80 GeV, v EW = 246 GeV, m h = 125 GeV, m H = m H ± = 400 GeV, where m W and m h are masses of the W boson and the SM Higgs, respectively.In addition, m H and m H ± are the masses of the heavier CP-even Higgs and the charged Higgs bosons, respectively.
To carry out the numerical computation, we have used a length unit in which v is normalized to unity.We have computed the relaxation in a 3D box with the size L x = L y = 8, L z = 12 and adopted the Dirichlet boundary conditions such that the field values on the boundaries are fixed to the asymptotic ones constructed above.After the relaxation time t = 20, the variation of the energy density per time is O(10 −5 ), thus we have regarded the convergence achieved.We also have confirmed that the field values converged up to the same order.Fig. 2 shows plots of the energy density, magnetic flux and Z flux of the solution.The energy is localized in the form of a string.Also, the magnetic flux is rising from the center, and it is clear that the magnetic monopole exists.Note that the flux is spherically symmetric at large distances but not near the monopole.The Z flux flows upward and downward on the string, which indicates that the two Z strings are attached to the monopole.
We would like to emphasize that the energy density at any cross sections with z = const.includes the power-law tail leading to the logarithmic divergence, but it is only related to the global U (1) a winding.Therefore, neither the Z-fluxes at the both sides of the configuration, which decay exponentially fast, nor the electromagnetic monopole suffer from the logarithmic divergence.We define the total energy of the monopole (monopole mass) as the difference between the energy of the configuration and that of the Z-string without the monopole, which has no divergence and is well-defined because the energy density of the monopole (the difference of their energy densities) is exponentially localized.In our numerical computation, the monopole mass is typically O(1) TeV.

VII. CONCLUSION AND OUTLOOK
We have constructed a stable magnetic monopole solution in 2HDM under a potential with the U (1) a symmetry and the (Z 2 ) C symmetry, which yields tan β = 1, based on full 3D numerical simulation.In the solution, the magnetic flux (a) energy density a / q k r u 6 e 3 j 5 / f 2 B g c G h 4 J D g 6 l n L M s q 3 y p G r q p p 1 R Z I f r m s G T Q h M 6 z 1 g 2 l 0 u K z t P K 0 U Z 7 P 1 3 h t q O Z x q 6 o W X y / J B c N r a C p s i A q H 5 z I C V 4 V j a g 6 n w t F 9 i K h g l 4π sin θ W /g emanates from the monopole, and the Z strings are attached to its both sides.The stability is ensured by the fact that the monopole can be regarded as a topological (Z 2 ) C kink interpolating between the two degenerate vacua on the lifted moduli space S 2 .
We should note that our monopole is essentially different from monopoles confined in vortices, which are put in the Higgs phase [56][57][58][59][60][61][62][63][64] 6 .In those cases, the magnetic flux and the physical degrees of freedom are confined in the string.On the other hand, our monopole is not in the Higgs phase since U (1) EM is not broken.Consequently, while the Z flux is confined, the magnetic flux is not confined but emanates spherically at large distances.
In this paper, we have assumed the U (1) a symmetry, and thus the corresponding NG boson appears after H gets a VEV.Since such massless modes are severely constrained by experiments, we should introduce explicit U (1) a breaking terms by switching on m 3 and α 5 , which give a mass to the NG boson (the CP-odd Higgs boson).This effect makes the topological strings attached by domain walls [45,46].The same is expected to happen with the string-monopole complex that we have considered above.As a result, our monopole (and string) are pulled by the wall and cannot be static.There is another option to avoid the NG boson, which is gauging the U (1) a symmetry [67].In such models, the monopoles are not attached by domain walls because the U (1) a symmetry is exact.Therefore, their evolution is easier to study.
Let us comment on phenomenological properties of our monopole.Interestingly, they did not be diluted by the cosmological inflation because it is produced at the electroweak symmetry breaking.This raises a question of whether our monopole causes the so-called cosmological monopole problem [68].A naive answer is no, because our monopole has a much lighter mass O(1) TeV than that of GUT monopoles, and thus it would not dominate the energy density of the universe at any epoch assuming that a few monopoles per horizon volume were produced at the symmetry breaking.Another famous bound on monopoles is the Parker bound [69,70].Again, our monopole would not conflict the bound because of its lightness.Of course, these arguments are quite naive.Therefore, it would be interesting to investigate the evolution of our monopoles (and strings) in the early universe and the abundance in the present universe.Such studies provide predictions for the current and nearly future monopole searches, and enable us to impose bounds on 2HDM.To this end, it will also be important to consider whether the monopoles can be generated, not only by the Kibble-Zurek mechanism [71,72], but also by the reconnection between the Z strings.
Another phenomenological application of our result may be baryogenesis.In Ref. [73], a sphaleron-like configuration is made from a Nambu monopole-antimonopole pair (dumbbell) in the SM.It is interesting to find out whether the same argument is applicable to our monopole.In addition, it is also interesting to consider the direct detection of the monopoles by collider experiments.
Finally, let us briefly comment on the case without the (Z 2 ) C symmetry, in which the Z fluxes of the two Z strings are not equal, and a difference occurs in the finite parts of those tensions.Therefore, the monopole is pulled to one side, resulting in the instability.A quantitative study of the instability will be done elsewhere.
g 4 0 a q m C 4 0 I I d m A i l r S M P g i / c B p r C B Y J U t M 9 o I S n a m m S x Z J j C b s t Y k d V 6 z L q y b t c M I 7 U t p z j S A 1 F m M E b X d E a 3 9 I f O 6 S 8 9 v F m r G d V o e 2 n I b H W 0 7 J d T + 8 N L 9 / 9 V V W X W + P W k e t e z x h a m I 6 9 K v P s R 0 7 6 F 3 d H X d w 9 u l 2 Y X x 5 r j d E T / x P 9 v u q E r u Y F b v 7 O P F 3 j x E E n 5 g P z L 5 3 4 N V n 7 k 8 o I X J r O F X P w V v R j F V 0 z I e 0 + h g D n M o y j n 7 u E E F 7 g 0 E s Y 3 I 2 d Q J 9 V I x J o h P A t j 5 h F l z 5 O J < / l a t e x i t > (b) magnetic flux < l a t e x i t s h a 1 _ b a s e 6 4 = " r G 8I v A h c Y X 6 K M T b 7 F 2 r M 7 k a Y U Z U = " > A A A C f n i c h V H L S s N A F D 2 N r 1 p f V T e K m 2 J R d G G d i m B x V X D j U l t b B S s l i d M 6 m C Y h m R a 1 F P f + g A t X C k V E0 I 9 w 4 w + 4 6 C e I y w p u X H i T B k R F v S E z Z 8 6 9 5 8 6 Z G c 0 2 h C s Z a 4 W U r u 6 e 3 r 5 w f 2 R g c G h 4 J D o 6 l n e t q q P z n G 4 Z l r O j q S 4 3 h M l z U k i D 7 9 g O V y u a w b e 1 w z U v v 1 3 j j i s s c 0 s e 2 3 y vo p Z N U R K 6 K o k q R i c L k h / J + p w 2 X 4 h 5 K S 6 F H i s Z 1 a N G M R p n C e Z H 7 C d I B i C O I D a s 6 D U K 2 I c F H V V U w G F C E j a g w q V v F 0 k w 2 M T t o U 6 c Q 0 j 4 e Y 4 G I q S t U h W n C p X Y Q x r L t N o N W J P W X k / X V + u 0 i 0 G / Q 8 o Y Z t g T u 2 F t 9 s h u 2 T N 7 / 7 V X 3 e / h e T m m W e t o u V 0 c O Z v I v v 2 r q t A s c f C p + t O z R A k p 3 6 s g 7 7 b P e K f Q O / r a y X k 7 u 5 q Z q c + y K / Z C / i 9 Z i z 3 Q C c z a q 9 7 c 5 J k L R O g B k t + v + y f I L y W S h D e X 4 + l E 8 B R h T G E a c 3 T f K 0 h j H R v I 0 b 6 n a O I O 9 w q U W W V B W e y U K q F A M 4 4 v o a Q + A D K 1 k v 0 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " r G 8 I v A h c Y X 6 K M T b 7 F 2 r M 7 k a Y U Z U = " > A A A C f n i c h V H L S s N A F D 2 N r 1 p f V T e K m 2 J R d G G d i m B x V X D j U l t b B S s l i d M 6 m C Y h m R a 1 F P f + g A t X C k V E0 I 9 w 4 w + 4 6 C e I y w p u X H i T B k R F v S E z Z 8 6 9 5 8 6 Z G c 0 2 h C s Z a 4 W U r u 6 e 3 r 5 w f 2 R g c G h 4 J D o 6 l n e t q q P z n G 4 Z l r O j q S 4 3 h M l z U k i D 7 9 g O V y u a w b e 1 w z U v v 1 3 j j i s s c 0 s e 2 3 y vo p Z N U R K 6 K o k q R i c L k h / J + p w 2 X 4 h 5 K S 6 F H i s Z 1 a N G M R p n C e Z H 7 C d I B i C O I D a s 6 D U K 2 I c F H V V U w G F C E j a g w q V v F 0 k w 2 M T t o U 6 c Q 0 j 4 e Y 4 G I q S t U h W n C p X Y Q x r L t N o N W J P W X k / X V + u 0 i 0 G / Q 8 o Y Z t g T u 2 F t 9 s h u 2 T N 7 / 7 V X 3 e / h e T m m W e t o u V 0 c O Z v I v v 2 r q t A s c f C p + t O z R A k p 3 6 s g 7 7 b P e K f Q O / r a y X k 7 u 5 q Z q c + y K / Z C / i 9 Z i z 3 Q C c z a q 9 7 c 5 J k L R O g B k t + v + y f I L y W S h D e X 4 + l E 8 B R h T G E a c 3 T f K 0 h j H R v I 0 b 6 n a O I O 9 w q U W W V B W e y U K q F A M 4 4 v o a Q + A D K 1 k v 0 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " r G 8 I v A h c Y X 6 K M T b 7 F 2 r M 7 k a Y U Z U = " > A A A C f n i c h V H L S s N A F D 2 N r 1 p f V T e K m 2 J R d G G d i m B x V X D j U l t b B S s l i d M 6 m C Y h m R a 1 F P f + g A t X C k V E 0 I 9 w 4 w + 4 6 C e I y w p u X H i T B k R F v S E z Z 8 6 9 5 8 6 Z G c 0 2 h C s Z a 4 W U r u 6 e 3 r 5 w f 2 R g c G h 4 J D o 6 l n e t q q P z n G 4 Z l r O j q S 4 3 h M l z U k i D 7 9 g O V y u a w b e 1 w z U v v 1 3 j j i s s c 0 s e 2 3 y v o p Z N U R K 6 K o k q R i c L k h / J + p w 2 X 4 h 5 K S 6 F H i s Z 1 a N G M R p n C e Z H 7 C d I B i C O I D a s 6 D U K 2 I c F H V V U w G F C E j a g w q V v F 0 k w 2 M T t o U 6 c Q 0 j 4 e Y 4 G I q S t U h W n C p X Y Q x r L t N o N W J P W X k / X V + u 0 i 0 G / Q 8 o Y Z t g T u 2 F t 9 s h u 2 T N 7 / 7 V X 3 e / h e T m m W e t o u V 0 c O Z v I v v 2 r q t A s c f C p + t O z R A k p 3 6 s g 7 7 b P e K f Q O / r a y X k 7 u 5 q Z q c + y K / Z C / i 9 Z i z 3 Q C c z a q 9 7 c 5 J k L R O g B k t + v + y f I L y W S h D e X 4 + l E 8 B R h T G E a c 3 T f K 0 h j H R v I 0 b 6 n a O I O 9 w q U W W V B W e y U K q F A M 4 4 v o a Q + A D K 1 k v 0 = < /l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " r G 8 I v A h c Y X 6 K M T b 7 F 2 r M 7 k a Y U Z U = " > A A A C f n i c h V H L S s N A F D 2 N r 1 p f V T e K m 2 J R d G G d i m B x V X D j U l t b B S s l i d M 6 m C Y h m R a 1 F P f + g A t X C k V E 0 I 9 w 4 w + 4 6 C e I y w p u X H i T B k R F v S E z Z 8 6 9 5 8 6 Z G c 0 2 h C s Z a 4 W U r u 6 e 3 r 5 w f 2 R g c G h 4 J D o 6 l n e t q q P z n G 4 Z l r O j q S 4 3 h M l z U k i D 7 9 g O V y u a w b e 1 w z U v v 1 3 j j i s s c 0 s e 2 3 y v o p Z N l a t e x i t > (c) Z flux < l a t e x i t s h a 1 _ b a s e 6 4 = " e i d 5 D i h a a w C y S M E t 7 n G l B Y / S o 4 G W P b I b 1 m I P 7 I 4 9 s f d f a z X c G m 0 v N Z q V j p Z b + Z G T y c T r v 6 o S z Q K H n 6 o / P Q s U s O p 6 1 c i 7 5 T L t W 6 g d f a V + 2 k q s 7 c w 2 5 t g l e y b / F 6 z J 7 u k G R u V F v d r m O 2 c I 0 A f E v z / 3 T 5 B a j M U J b y + H 1 2 P e V / g x h R l E 6 b 1 X s I 5 N b C F J 5 9 Z x j h v c + t 6 k a S k q L X R S J Z + n G c e X k J Y + A C U 9 k S U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e i d 5 D i h a a w C y S M E t 7 n G l B Y / S o 4 G W P b I b 1 m I P 7 I 4 9 s f d f a z X c G m 0 v N Z q V j p Z b + Z G T y c T r v 6 o S z Q K H n 6 o / P Q s U s O p 6 1 c i 7 5 T L t W 6 g d f a V + 2 k q s 7 c w 2 5 t g l e y b / F 6 z J 7 u k G R u V F v d r m O 2 c I 0 A f E v z / 3 T 5 B a j M U J b y + H 1 2 P e V / g x h R l E 6 b 1 X s I 5 N b C F J 5 9 Z x j h v c + t 6 k a S k q L X R S J Z + n G c e X k J Y + A C U 9 k S U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e i d 5 D i h a a w C y S M E t 7 n G l B Y / S o 4 G W P b I b 1 m I P 7 I 4 9 s f d f a z X c G m 0 v N Z q V j p Z b + Z G T y c T r v 6 o S z Q K H n 6 o / P Q s U s O p 6 1 c i 7 5 T L t W 6 g d f a V + 2 k q s 7 c w 2 5 t g l e y b / F 6 z J 7 u k G R u V F v d r m O 2 c I 0 A f E v z / 3 T 5 B a j M U J b y + H 1 2 P e V / g x h R l E 6 b 1 X s I 5 N b C F J 5 9 Z x j h v c + t 6 k a S k q L X R S J Z + n G c e X k J Y + A C U 9 k S U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e i d 5 D i h a a w C y S M E t 7 n G l B Y / S o 4 FIG. 2: Plots for the numerical solution constructed by the relaxation.In all plots, v is normalized to unity.(a): Energy density.The color represents its value, where red is the largest and blue is the smallest.There is a string-like object that contains the energy density along the z axis.(b): Magnetic flux.The direction of the arrows indicates that of the flux.Also, the color and size of the arrows indicate the flux density, where red is the strongest, blue is the weakest.We can see the existence of a magnetic monopole at the center from which the magnetic flux emanates.The flux is spherical at large distances, but not in the vicinity of the monopole.(c): Z flux.The direction, color and size of the arrows are the same as those for the magnetic flux.The Z flux flows upward and downward along the string from the monopole.