Dyon in the $SU(2)$ Yang--Mills theory with a gauge-invariant gluon mass toward quark confinement

In the previous paper, we have shown the existence of magnetic monopoles in the pure $SU(2)$ Yang--Mills theory with a gauge-invariant mass term for the gluon field being introduced. In this paper, we extend our previous construction of magnetic monopoles to obtain dyons with both magnetic and electric charges. In fact, we solve under the static and spherically symmetric ansatz the field equations of the $SU(2)$"complementary"gauge-scalar model, which is the $SU(2)$ Yang--Mills theory coupled to a single adjoint scalar field whose radial degree of freedom is eliminated. We show that the novel dyon solution can be identified with the gauge field configuration of a dyon with a minimum magnetic charge in the massive Yang--Mills theory. Moreover, we compare the dyon of the massive Yang--Mills theory obtained in this way with the Julia--Zee dyon in the Georgi--Glashow gauge-Higgs scalar model and the dyonic extension of the Wu--Yang magnetic monopole in the pure Yang--Mills theory. Finally, we identify the novel dyon solution found in this paper with a dyon configuration on $S^1 \times \mathbb{R}^3$ space with nontrivial holonomy and propose to use it to understand the confinement/deconfinement phase transition in the Yang--Mills theory at finite temperature, instead of using the dyons constituting the Kraan--van Baal--Lee--Lu caloron.


I. INTRODUCTION
Quark confinement is a long-standing problem to be solved in the framework of quantum chromo-dynamics (QCD). One of the most promising scenarios for quark confinement is the dual superconductivity picture [1] for QCD vacuum. For this hypothesis to be realized, the existence of the relevant magnetic objects and their condensations are indispensable. In the gauge-scalar model, it is indeed well known that the 't Hooft-Polyakov magnetic monopole [2] exists as a topological soliton solution of the field equations of the Georgi-Glashow gauge-Higgs scalar model in which a single scalar field belongs to the adjoint representation of the gauge group SU (2). However, such a scalar field is not introduced in the original action of QCD. Thus, we are faced with the problem of showing the existence of magnetic monopoles in the Yang-Mills theory in absence of the scalar field. See e.g., [3] for a review.
In a previous paper [4], nevertheless, we have succeeded to obtain the magnetic monopole configuration with a nontrivial magnetic charge in the pure SU (2) Yang-Mills theory when a gauge-invariant mass term for the gauge field [5] is introduced without any scalar field. We call this theory the massive Yang-Mills theory and call the resulting magnetic monopole the Yang-Mills magnetic monopole. This result follows from the recent proposal for obtaining the gauge field configurations in the pure Yang-Mills theory from solutions of field equations in the "complementary" gauge-scalar model [5] in which the radial degree of freedom of a single adjoint scalar field is frozen. The gaugeinvariant mass term is obtained through change of variables and a gauge-independent description [5,6] of the Brout-Englert-Higgs (BEH) mechanism [7], which neither relies on the spontaneous breaking of gauge symmetry nor on the assumption of a nonvanishing vacuum expectation value of the scalar field.
As first shown by Julia and Zee [9], it is possible to provide magnetic monopoles with electric charges, in fact, a 't Hooft-Polyakov magnetic monopole [2] can have the electric charge in addition to the magnetic charge, which is called a dyon [10]. We show in this paper that the dyon configuration exists in the pure SU (2) Yang-Mills theory with a gauge-invariant mass term of the gauge field, i.e., the massive Yang-Mills theory. We call this dyon the Yang-Mills dyon. In fact, we solve under the static and spherically symmetric ansatz the field equations of the SU (2) "complementary" gauge-scalar model and obtain a gauge field configuration for a dyon with a minimum magnetic charge in the massive Yang-Mills theory. In particular, we compare the dyon obtained in this way in the massive Yang-Mills theory with the Julia-Zee dyon in the Georgi-Glashow gauge-Higgs model and the dyonic extension of the Wu-Yang magnetic monopole [11] in the pure Yang-Mills theory.
It is well known that topological solitons play very important roles in non-perturbative investigations of the Yang-Mills theory. For example, the Kraan-van Baal-Lee-Lu (KvBLL) caloron (instanton) [12] has been extensively used to reproduce the confinement/deconfinement phase transition in the Yang-Mills theory at finite temperature [13]. The KvBLL caloron is a topological soliton solution of the (anti-)self-dual equation of the SU (2) Yang-Mills theory on S 1 × R 3 space with a nontrivial instanton charge as the topological invariant, which consists of Bogomol'nyi-Prasad-Sommerfield (BPS) dyons having both electric and magnetic charges with a nontrivial holonomy at spatial infinity. In contrast, our dyon solution is a non-BPS solution of the (non-self-dual) field equation of the "complementary" gauge-scalar model. Our dyon has the nonvanishing asymptotic value as the nontrivial holonomy at spatial infinity, which is a common property to be comparable with (anti-)self-dual dyons as the constituents of the KvBLL calorons. The dyon solution of the "complementary" gauge-scalar model is identified with the dyon configuration of the massive Yang-Mills theory with a gauge-invariant mass term of the gauge field without any scalar field, which is regarded as the low-energy effective model of the Yang-Mills theory with a mass gap. Thus, we can propose another scenario for investigating the confinement/deconfinement phase transition in the Yang-Mills theory at finite temperature based on the novel non-self-dual dyon solution found in this paper.
This paper is organized as follows. In section II, we review the procedure [5] for obtaining the massive SU (2) Yang-Mills theory from the "complementary" SU (2) gauge-adjoint scalar model. In section III, we give a brief review of the Julia-Zee dyon in the (3 + 1)-dimensional Minkowski spacetime. In section IV, we give a Yang-Mills dyon solution of an SU (2) gauge-scalar model in which the radial degree of freedom of a single adjoint scalar field is frozen. In section V, we give the decomposition of the gauge field for the Yang-Mills dyon solution. In section VI, we give the magnetic and electric fields of the Yang-Mills dyon solution. In section VII, we give the energy density and the static mass of a Yang-Mills dyon. In section VIII, we discuss how the Yang-Mills dyon can be responsible for confinement/deconfinement phase transition of the Yang-Mills theory at finite temperature. The final section IX is devoted to conclusion and discussion.

II. THE MASSIVE YANG-MILLS THEORY "COMPLEMENTARY" TO THE GAUGE-ADJOINT SCALAR MODEL
In this section, we review the procedure [5] for obtaining the massive SU (2) Yang-Mills theory from the "complementary" SU (2) gauge-adjoint scalar model. For this purpose, we introduce the two products for the Lie-algebra valued fields P := P A T A and Q = Q A T A (A = 1, 2, 3): where T A are the generators of the Lie algebra su(2) of the group SU (2). We choose the Hermitian basis T A by using the Pauli matrices σ A (A = 1, 2, 3) as We introduce the Hermitian SU (2) gauge field A µ (x) and the Hermitian scalar field φ(x) by A µ (x) := A A µ (x)T A and φ(x) = φ A (x)T A as the Lie algebra su(2) valued fields. Then, we introduce the SU (2) gauge-adjoint scalar model by the Lagrangian density where u is the Lagrange multiplier field to incorporate the radially fixing constraint Here F µν denotes the field strength of the SU (2) gauge field A µ and D µ [A ]φ is the covariant derivative of the scalar field φ(x) defined by First of all, we construct a composite vector boson field X µ (x) from A µ (x) andφ(x) as by introducing the normalized scalar fieldφ which can be identified with the color direction field n(x) in the gauge-covariant field decomposition of the gauge field, see [3,8]. Notice that X µ (x) transforms according to the adjoint representation under the gauge transformation U (x) ∈ SU (2): Then the kinetic term of the scalar field is identical to the mass term of the vector field X µ (x): as long as the radial degree of freedom of the scalar field is fixed [5]. It is obvious that the obtained mass term of X µ (x) is gauge-invariant by observing (10). Therefore, X µ (x) can be identified with the massive component without breaking the original gauge symmetry. This gives a gauge-independent definition of the massive modes of the gauge field in the operator level. It should be emphasized that we do not need to choose a specific vacuum of φ(x) and hence no spontaneous symmetry breaking of the gauge symmetry occurs. By using the definition of the massive vector field X µ (x), the original gauge field A µ (x) is separated into two pieces [3,5]: where the field V µ (x) can be written in terms of A µ (x) andφ(x): Here V µ (x) is called the restricted (or residual) part which is expected to give the dominant contribution to quark confinement, while X µ (x) is called the remaining (or broken) part which is identified with the massive mode which is expected to decouple in the low-energy or long-distance region. Then, we regard a set of field variables {c µ (x), X µ (x),φ(x)} as obtained from {A µ (x),φ(x)} based on a change of variables: and identify c µ (x), X µ (x) andφ(x) with the fundamental field variables for describing the massive Yang-Mills theory anew, which means that we should perform the quantization with respect to the variables {c µ (x), X µ (x),φ(x)} appearing in the path-integral measure.
In the gauge-scalar model, A µ (x) andφ(x) are independent field variables. However, the Yang-Mills theory should be described by A µ (x) alone. Hence the scalar field φ(x) must be supplied by the gauge field A µ (x) due to the strong interactions, or in other words, φ(x) should be given as a functional of the gauge field A µ (x).
Moreover, the independent degrees of freedom of the original gauge field A A µ (x) in the pure SU (2) Yang-Mills Here, we have omitted the infinite degrees of freedom of the space-time points. On the other hand, the new field variables have independent degrees of freedom: where the massive vector field X µ (x) obeys the condition: We can therefore observe that the theory with the new field variables has two extra degrees of freedom if we wish to obtain the (pure) Yang-Mills theory from the "complementary" gauge-scalar model. These extra degrees of freedom are eliminated by imposing the two constraints which we call the reduction condition. We choose e.g., the reduction condition The reduction condition indeed eliminates the two extra degrees of freedom introduced by the radially fixed scalar field into the Yang-Mills theory, since Following the Faddeev-Popov procedure, we insert the unity to the functional integral to incorporate the reduction condition: where χ θ := χ[A , φ θ ] is the reduction condition written in terms of A µ (x) and φ θ which is the local rotation of φ(x) by θ = θ(x) = θ A (x)T A and ∆ red := det δχ θ δθ denotes the Faddeev-Popov determinant associated with the reduction condition χ = 0. Then, we write the vacuum-to-vacuum amplitude of the gague-scalar model subject to the reduction condition is translated into the massive Yang-Mills theory with the gauge-invariant mass term of the field X as where the Jacobian J associated with the change of variables is equal to one, J = 1 [3]. Therefore, we obtain the massive Yang-Mills theory which keeps the original gauge symmetry: The obtained massive Yang-Mills theory indeed has the same degrees of freedom as the usual Yang-Mills theory because the massive vector boson X µ (x) is constructed by combining the original gauge field A µ (x) and the normalized scalar fieldφ(x) whereφ(x) is now a (complicated) functional of A µ (x) obtained by solving the reduction condition (15). It should be remarked that the solutions of the field equations of the gauge-scalar model satisfy the reduction condition automatically, although the converse is not true [5]. The field equations besides the constraint equation (5) are obtained as We take the inner product of (21) and φ(x) and use (5) to obtain which is used to eliminate the Lagrange multiplier field u in (21). Indeed, the field equations (20) and (21) are rewritten in terms of A µ (x) andφ(x) into By applying the covariant derivative D ν [A ] to the equation (23), the reduction condition is naturally induced: Moreover, by taking the exterior product of (24) andφ(x), the reduction condition is induced again: Hence, the simultaneous solutions of the coupled field equations (23) and (24) automatically satisfy the reduction condition (15). From this relation, we find that the solutions of the coupled field equations of the gauge-scalar model (23) and (24) can become the field configurations satisfying the reduction condition (15), which gives the field configuration to be taken into account in constructing the massive Yang-Mills theory through the path-integral (18).

III. JULIA-ZEE DYON SOLUTION IN THE GEORGI-GLASHOW MODEL
In this section, we give a brief review of the Julia-Zee dyon [9] in the (3 + 1)-dimensional Minkowski spacetime R 1,3 . The Georgi-Glashow model is introduced by the Lagrangian density where g, λ and v > 0 are respectively the gauge coupling constant, the scalar coupling constant and the value of the magnitude |φ(x)| of the adjoint scalar field φ(x) at the vacuum which is to be realized at infinity |x| = ∞. By varying the action with respect to the fields A µ (x) and φ(x), the field equations are obtained as The Julia-Zee ansatz with a unit magnetic charge is given by where Roman indices j, k run from 1 to 3 and r is the radius r := x 2 + y 2 + z 2 in the three-dimensional space with the Cartesian coordinates (x, y, z). Note that the electric charge cannot be specified at this stage in contrast to the magnetic charge specified by the form of the ansatz (31) for A A j (x) which has the same form as that used for obtaining the 't Hooft-Polyakov magnetic monopole with a unit magnetic charge [2]. This is because the electric charge depends on the asymptotic value of the profile function a(r) in A A 0 (x), to be obtained by solving the coupled field equations simultaneously for the other unknown functions f (r) and h(r), as will be performed below.
The field equations (29) and (30) are rewritten in terms of the profile functions a, f , and h as The scaled dimensionless variable ρ, and the scaled functions a, f , and h of ρ defined by are introduced to make the field equations (32)-(34) dimensionless: where the prime denotes the derivative with respect to ρ hereafter. In order to determine the boundary conditions, we consider the static energy E given by where we have defined the mass scale M X by and the energy density e(ρ) by For the energy (39) to be finite, the regularity of the fields at ρ = 0 is required: For large ρ, by the same reason, the scalar field must go to its vacuum expectation value at infinity: As the gauge field A A µ (x) goes to the pure gauge form at ρ → ∞, the profile function f (ρ) should take In order to solve the field equations (36)-(38), however, the asymptotic value a ∞ of a(ρ) must be specified: Notice that a ∞ is not completely arbitrary. As will be shown below, indeed, for f (ρ) not to oscillate at large ρ so that the spatial components A A j (x) of the gauge field become the pure gauge form at ρ ≈ ∞, the constant a ∞ should take the value |a ∞ | < 1. Notice that, if a(ρ) is a solution of equations (36)-(38), then −a(ρ) is also a solution of them. Therefore, a ∞ is restricted to take the nonnegative value 0 ≤ a ∞ < 1 without loosing the generality. The solution a(ρ) ≡ 0 with a vanishing A 0 component corresponds to the 't Hooft-Polyakov magnetic monopole, which leads to a ∞ = 0. Since the parameters g and v are factorized out, dyon solutions are distinguished by the value of a ∞ and the scalar coupling λ.
We further consider the asymptotic forms of the profile functions. For small ρ, ρ ≈ 0, we assume the power-series expansion in ρ By substituting the above power-series expansion into the field equations (36)-(38), the asymptotic forms for small ρ are obtained in agreement with the boundary conditions (42) for small ρ.
For large ρ, ρ → ∞, by introducing b(ρ) and k(ρ) by the field equations (36)-(38) reduce to the linear differential equations These equations can be solved independently under the boundary conditions b(∞) = 0, f (∞) = 0, and k(∞) = 0 as where C, F and H are arbitrary constants. The asymptotic solution (55) for large ρ indicates that for the profile function f (ρ) of the spatial components A j (x) of the gauge field not to oscillate at large ρ, the constant a ∞ should take the value |a ∞ | < 1. Therefore, a ∞ can be restricted to 0 ≤ a ∞ < 1 without loosing the generality. We define the chromomagnetic field B A j (x) and chromoelectric field E A j (x) by and the magnetic charge q m and electric charge q e by where we have introduced the normalized scalar field By using the asymptotic forms (54)-(55), the magnetic and electric charges q m and q e are calculated as where we assume that the coordinates x j are also dimensionless as well as ρ. We find that the magnetic charge q m is indeed nontrivial and has a minimal value in unit of 4π/g. We also find that the coefficient C of the next-leading term of a(ρ) in (54) is nothing but the ratio of the charges q e /q m : It should be noticed that, although some physical quantities such as the chromoelectric field E A j (x) and the electric charge q e do not depend on the asymptotic value a ∞ of a(ρ), they depend on the next-leading coefficient C of a(ρ), namely, the ratio of the charges. Therefore, in order to compare the solutions and the corresponding physical quantities with the same electric charge q e by varying the scalar coupling λ, we adopt the following boundary condition If we restrict the solution a(ρ) to non-negative one a(ρ) > 0 as mentioned in the above, then we have only a ∞ ≥ 0 and C = q e /q m ≥ 0.

IV. CONSTRUCTION OF THE YANG-MILLS DYON
Next, we discuss the dyon in the massive SU (2) Yang-Mills theory through the "complementary" SU (2) gaugescalar model. Taking the Julia-Zee ansatz (31), the field equations (20), (21), and the constraint (5) become Here the last equation is nothing but the radially fixing constraint and can be used to eliminate h(r) from the other equations to obtain Hence, the Lagrange multiplier field u = u(r) can be determined through (71) once the remaining two equations (69) and (70) for a(r) and f (r) are solved. In order to make the field equations dimensionless, we define the dimensionless variable ρ = gvr and the rescaled functions of ρ: a(r) = gva(ρ), f (r) = f (ρ). Then (69) and (70) read By repeating the same procedure for obtaining the boundary condition as the Julia-Zee dyon, we find that it is sufficient to impose the following boundary conditions for the Yang-Mills dyon: where C is an arbitrary constant. We find that these conditions are enough to guarantee the regularity of the fields at the origin ρ = 0 to obtain a finite energy. For large ρ, the equations (72) and (73) reduce to (51) and (52) respectively and therefore the profile functions behave like the Julia-Zee dyon: For small ρ, however, the asymptotic forms are much different from the Julia-Zee case. To realize this, we shall linearize the field equations by assuming f (ρ) = 1 + g(ρ) and |g(ρ)|, |a(ρ)| ≪ 1. Then, the field equations (72) and (73) become The first equation (77) is solved as The second equation (78) is the same as the Yang-Mills monopole case and can be solved as [4] f Here note that the equation cannot be satisfied by a simple power-series of ρ without the logarithmic term. The field equations (72) and (73) can be solved numerically, which is shown in Fig. 1 for e.g., C = 0.5. The Julia-Zee dyon solution with a large coupling λ ≫ 1 approaches the Yang-Mills dyon except for the neighborhood of the origin ρ ≈ 0. This is the similar situation to the Yang-Mills monopole. For a(ρ), the naive limit λ → ∞ of the Julia-Zee dyon completely agrees with the Yang-Mills dyon.
In particular, f (ρ) ≡ 0 is also the solution of (73) and hence the exact solution of (73) is given by These solutions, however, do not satisfy the boundary conditions (74) for ρ → 0 and a(ρ) diverges at the origin ρ = 0 for C = 0. In view of these, the dyon constructed by f (ρ) ≡ 0 and (81) has a diverging energy and is regarded as a dyonic extension of the Wu-Yang monopole. For C = 0, the solution leads a(ρ) ≡ 0, which means that this solution is nothing but the Wu-Yang monopole with a vanishing electric charge.

V. GAUGE FIELD DECOMPOSITION FOR A DYON
In what follows, we shall omit the tilde (˜) for the profile functions f and h to simplify the notation. In the present ansatz (31), the normalized scalar fieldφ(x) with h(r) = +1, leads to the decomposed fields which are explicitly written as Notice that the time component of X µ (x) is identically zero X 0 (x) = 0, according to (8) by taking into account the facts thatφ(x) is time-independent and that the time component A 0 (x) of the gauge field and the normalized scalar fieldφ(x) are parallel in the color space gA A 0 (x) = a(r)φ A (x): In what follows, we adopt the polar coordinate system (r, θ, ϕ) for the spatial coordinates. Then the decomposed field have the following components: and where we have defined Fig. 2 is the plot of the fields A 0 , A ,V , and X as functions of ρ = gvr. We find that the spatial component A(r) of the original gauge field A µ (x) is indeed regular at the origin: It should be remarked that although both decomposed fields V (x) and X(x) are singular at the origin, the singularities cancel between the two decomposed fields to yield the regular A(r). On the other hand, the time component A 0 (r) of the gauge field is regular at the origin and behaves around the origin although the time component X 0 (x) of the massive mode is absent, X 0 (x) = 0 and the time component A 0 (r) is identical to the restricted field V 0 (r). By choosing the gauge transformation matrix U (x) as the fields are transformed into where we have defined We find that the dyonic contribution appears in the Wu-Yang potential V (x), while there are no effects in the massive mode X (x).

VI. CHROMOELECTRIC AND CHROMOMAGNETIC FIELDS OF A DYON
In this section, all the expressions are given for the Julia-Zee dyon, i.e., radially variable case. The expressions for the Yang-Mills (radially fixed) dyon can be easily obtained by setting h(ρ) = 1.
In the similar way to the Yang-Mills monopole [4], we examine the magnetic charge q m and electric charge q e obtained by the chromomagnetic field B A j (x) and chromoelectric field E A j (x): The magnetic charge q m and its density ρ m (r) are defined by Similarly, the electric charge q e and its density ρ e (r) are defined by The charge densities ρ m (r) and ρ e (r) can be written in terms of the profile functions: ρ e (r) =r 2 d dr a(r) d dr h(r) + 2 a(r)h(r)f 2 (r). Fig. 3 is the plots of the charge densities ρ m and ρ e as functions of ρ = gvr at C = 0.5. We also illustrate the a ∞ -dependence of the ratio C of the charges, C = q e /q m , which is shown in Fig. 4. The a∞ dependence of the ratio qe/qm between two charges qe, qm, namely C. Throughout in this paper, we fix the magnetic charge qm to the unit 4π/g, which means that C is nothing but the electric charge qe in unit of 4π/g.
Next, we investigate the behavior of the chromomagnetic field B A j (x) and chromoelectric field E A j (x), especially around r ≈ 0. To do so, we return to the polar coordinate representation: In order to define the gauge-invariant field strength, we take the inner product between (103) or (104) andφ(x), and the other components are zero. Fig. 5 is the plot of the gauge-invariant chromomagnetic and chromoelectric field strengths as functions of ρ = gvr at C = 0.5. The chromoelectric field E r (x) is regular at the origin even for the Yang-Mills dyon. However, the chromomagnetic field B r (x) of the Yang-Mills dyon diverges logarithmically at the origin, just like the Yang-Mills monopole.

VII. ENERGY DENSITY AND STATIC MASS OF A DYON
We define the energy integral I as a function of a ∞ and λ by integrating the energy density e(ρ) defined by (41) so that the energy E takes the form: Fig. 6 is the plot of the energy density e(ρ) as a function of ρ obtained from the solution a(ρ), f (ρ) and h(ρ) at C = 0.5, which should be compared with the Julia-Zee solution. We find that at the origin the energy density e(ρ) of the Yang-Mills dyon takes the value e(0) = 1, while the Julia-Zee dyons behave e(0) = 0 for any values of 0 ≤ λ < ∞. This difference is caused by the radially fixing condition h(ρ) = 1. The sixth term h 2 (ρ)f 2 (ρ) in (41) survives at the origin in the Yang-Mills dyon due to h(0) = 1, while for the Julia-Zee dyon it vanishes since h(0) = 0.
In the BPS limit λ = 0 of the 't Hooft-Polyakov monopole (a ∞ = 0), the integral I takes the value one:  which leads to For a fixed value of the scalar coupling λ, the energy E is monotonically increasing in a ∞ as seen in the left panel of Fig. 7 or in C = q e /q m as seen in the right panel of Fig. 7. The right panel of Fig. 7 has been obtained in [14] for the Julia-Zee dyon with finite values of λ and we have added the Yang-Mills dyon, which partially corresponds to the limit of λ → ∞ of the Julia-Zee dyon. For the Yang-Mills dyon and the Julia-Zee dyon for a sufficiently large coupling λ ≫ 1 above a critical value λ ′ of λ, the electric charge q e cannot reach the maximal limit q e /q m = 1: the maximal value of q e /q m is obtained in a numerical way as q e q m a∞→1,λ→∞ = 0.828.
By using the maximal value of the energy integral for the Yang-Mills dyon, the maximal value of the static mass of the Yang-Mills dyon can be estimated as where we have used the value for the off-diagonal gluon mass M X = 1.2GeV obtained by the preceding studies on a lattice [15] and the typical value of the running coupling constant α s (p) := g 2 (p)/4π ≈ 2.3 ± 0.1 at p ≃ M X ≈ 1.2GeV obtained in [16]. The Yang-Mills dyon mass, 1.18GeV, obtained in this paper is equal to the heaviest one in the family of Julia-Zee dyons in the Georgi-Glashow model, since the energy integral (106) is monotonically increasing in the coupling constant λ and the asymptotic value a ∞ of the time component of the gauge field. The Yang-Mills dyon mass, 1.18GeV, is 27% larger than the Yang-Mills monopole one: 0.93GeV found in [4], it still remains in the same order of the off-diagonal gluon mass: M X = 1.2GeV. In view of these, the existence of the Yang-Mills dyon with a reasonable mass tells us that the Yang-Mills dyons can play the role of the quark confiner instead of the Yang-Mills monopoles if their condensations occur according to the dual superconductor picture.

VIII. CONFINEMENT/DECONFINEMENT PHASE TRANSITION AND THE YANG-MILLS DYONS
In this section, we consider the Yang-Mills dyons on S 1 × R 3 space. We introduce the coordinates x = (τ, x) ∈ S 1 × R 3 with the "time" coordinate τ = x 4 for S 1 and spatial coordinates x = (x 1 , x 2 , x 3 ) for R 3 . Suppose we have performed the Wick rotation to the coordinates and the field in the Minkowski spacetime to obtain the Euclidean counterparts: τ = ix 0 and A 4 (τ, x) = −iA 0 (x 0 , x). Then the Euclidean action S mYM E is obtained as where Greek indices µ and ν run from 1 to 4. We have introduced the period T −1 of the "time" direction τ , which is regarded as the inverse of temperature. Note that the fields on S 1 × R 3 must be periodic in τ with the period T −1 on S 1 , that is to say, the fields must satisfy the periodic boundary condition A. The self-dual dyon and the KvBLL caloron in the massless Yang-Mills theory Before obtaining the Yang-Mills dyon solution, let us review the conventional dyon solution in the pure massless SU (2) Yang-Mills theory on S 1 × R 3 space. The action S YM E is given by removing the scalar field φ from (113): This action is nonnegative and has a lower bound: where we have introduced the Hodge dual ⋆ F µν of the field strength tensor F µν by The equality of (116) holds if and only if the equation is satisfied. The equation (118) is called the self-dual equation for the plus sign on the right hand side, and the anti-self-dual equation for the minus sign. In what follows, we will concentrate the self-dual equation. First, we adopt the ansatz for S 1 × R 3 space where Roman indices j and k run from 1 to 3 and r is the radius in R 3 , i.e., r = √ x j x j . In fact, the ansatz (140) is τ -independent and hence it trivially satisfies the periodic boundary condition (114). The nontrivial components of the self-dual equation (118) with the ansatz (119) is given by which is written in terms of the profile functions a(r) and f (r) as The solution of the self-dual equations (121) is exactly obtained as where V > 0 is an arbitrary parameter with a dimension of mass, which is related to the asymptotic holonomy. For later convenience, we set By introducing the dimensionless variable ρ and functions a and f in the same way as (35) and using the relation (124), the solutions (122) and (123) are cast into Notice that the solutions (125) and (126) of the self-dual equations (121) also satisfy the second-order field equations which is obtained by substituting the ansatz (119) to the Yang-Mills field equation as suggested from (36) and (37) by the replacement a → −ia and h → 0. By recalling the asymptotic behavior (76) of the profile function a(ρ), the ratio C of the charges for the self-dual dyon is fixed for any values of a ∞ . The action S YM E is also rewritten as where we have defined the dimensionless energy I by The dimensionless energy I is calculated by using the solutions (125) and (126) as and hence the action S YM The gauge field A µ (x) for the self-dual dyon is asymptotically written as by performing the gauge transformation to the unitary gauge (or the stringy gauge) in which the A 4 component is constant and diagonal at spatial infinity. Second, the pure massless SU (2) Yang-Mills theory has an another topological soliton solution. The KvBLL calorons are the solution of the (anti-)self-dual equation (118) of the pure massless Yang-Mills theory in S 1 × R 3 space with a nontrivial holonomy. Indeed, the gauge field A µ (x) of a KvBLL caloron is given by [12] gA µ (τ, x) = V δ µ4 T 3 +η 3 µν ∂ ν log x 3 with three functionsψ, ψ, and χ defined bŷ whereη A µν is the 't Hooft symbol defined byη A µν = ǫ Aµν − δ Aµ δ ν4 + δ Aν δ µ4 andV := 2πT − V . This solution indeed has the periodicity T −1 , i.e., A µ (τ + T −1 , x) = A µ (τ, x). Let  [17] where the direction is chosen to be along the third spatial direction. The core sizes of the two dyons M, L are respectively given by V −1 andV −1 . See Fig. 8.
The nontrivial holonomy of the KvBLL caloron originates from the first term of (135) for V = 0. If we take V = 0, the KvBLL caloron reduces to the Harrington-Shepard (HS) caloron solution with a trivial holonomy [18]. In the zero temperature limit T → 0, the KvBLL caloron reduces to the BPST instanton [17] of size ̺.
The constituent dyon can be identified if it is in the vicinity of one of its constituent dyons and far away from the other, namely, at large separations. For instance, near the M dyon center and far away from the L dyon (r ≫ 1/V ), the KvBLL caloron solution (135) exhibits the asymptotic behavior of the M dyon by performing the gauge transformation to the unitary gauge in which the A 4 component is constant and diagonal at spatial infinity. This is nothing but the self-dual dyon (134).
B. The Yang-Mills dyon on S 1 × R 3 space Let us get back to the massive Yang-Mills theory (113). We adopt the ansatz for S 1 × R 3 space This ansatz has the same form as the "static" Julia-Zee ansatz on the (3 + 1)-dimensional Minkowski spacetime R 1,3 , which is obtained by the replacement a(r) → −i a(r) according to the Wick rotation. Therefore, we can use the arguments in the previous sections. By introducing the dimensionless variable ρ and functions a and f in the same way as (35), the equations (72) and (73) are rewritten as We set the boundary conditions of a(ρ) and f (ρ) as Notice that in the Euclidean space there is no restriction for the asymptotic value a ∞ of a(ρ), since the asymptotic form of f (ρ) for large ρ satisfies the boundary condition f (∞) = 0 and exhibits no oscillating behavior for any values of a ∞ . This differs from the dyons in the Minkowski spacetime. We find that if a(ρ) is a solution of the equations (141) and (142), then −a(ρ) is also a solution of them. Therefore, a ∞ is restricted to take the nonnegative value a ∞ ≥ 0 without loosing the generality. The last term −f (ρ) on the left hand side of (142) originates from the kinetic term of the radially fixed scalar field φ(x), or equivalently, the gauge-invariant gluon mass term. If the term −f (ρ) is absent from the equation (142), which means the absence of the gauge-invariant gluon mass term, the system (113) reproduces the pure massless Yang-Mills theory (115). If a ∞ is sufficiently large, a ∞ ≫ 1, the term −f (ρ) in (142) can be negligible and hence the field equation (142) can be approximated by (128). This means that the Yang-Mills dyon behaves as the self-dual dyon (125) and (126) for large a ∞ except the neighborhood of the origin ρ ≈ 0. From (129), we observe that the ratio C of the charges for the Yang-Mills dyon is equal to 1 for large a ∞ : In the Yang-Mills dyon, the asymptotic value a ∞ of the profile function a(ρ) can be regarded as a function of the ratio C of the charges q e and q m , namely, C = q e /q m . Fig. 10 is a plot of the ratio C of the charges as a function of a ∞ for the fixed magnetic charge q m = 4π/g. We find that the electric charge q e of the Yang-Mills dyon depends on the asymptotic value a ∞ of the profile function a(ρ), while the electric charge of the self-dual dyon is fixed q e = q m . We confirm numerically that the Yang-Mills dyon differs from the self-dual dyon for any finite a ∞ , but approaches the self-dual dyon in the limit a ∞ → ∞. This means that the upper bound of the electric charge q e exists and is given by The action (113) for a single Yang-Mills dyon is given by using (106) and (107) where we have defined the dimensionless energy integral I(a ∞ ) by which should be compared with (131) for the self-dual dyon. The function I(a ∞ ) of a ∞ is monotonically increasing in a ∞ with the lower bound I(a ∞ ) ≥ I(a ∞ = 0) = 1.787 based on the numerical calculations as given in the left panel of Fig. 11. The action for a KvBLL caloron is given by S = 8π 2 /g 2 , which is T -independent. This is not the case for both the self-dual dyon (133) and the Yang-Mills dyon (148). The a ∞ -dependence of I(a ∞ ) is obtained as follows. By taking the derivative of (149) with respect to a ∞ , we have where we have integrated by parts. The first term of (150) vanishes, since the boundary conditions (144) of f (ρ) is a ∞ -independent. The third term of (150) also vanishes due to the field equations (141) and (142). Thus, we obtain We find that from (146) and (151), I(a ∞ ) behaves as and hence I(a ∞ ) linearly diverges in a ∞ : This is consistent with the numerical calculation shown in the left panel of Fig. 11. Notice that the gauge field A µ of the Yang-Mills dyon has the asymptotic behavior dominated by V µ , A µ ≈ V µ , at long distance, e.g., in the unitary gauge, which is the same as the asymptotic field (139) of a constituent dyon of the KvBLL caloron and hence the asymptotic field (134) of the self-dual dyon. Note that the self-dual dyon and Yang-Mills dyon no longer have τ -dependence and hence they are trivially periodic, unlike the KvBLL caloron.
C. The Yang-Mills dyon versus the self-dual dyon from the KvBLL caloron toward confinement/deconfinement phase transition In order to discuss the confinement/deconfinement phase transition in the Yang-Mills theory at finite temperature, we define the Polyakov loop operator L(x) as L(x) := 1 tr (1) tr P exp ig where P denotes the path-ordering prescription. The asymptotic holonomy P ∞ is defined by the Polyakov loop operator at the spatial infinity By performing the gauge transformation to the unitary gaugeφ A (x) = δ A3 , so that the "time" component A 4 (x) of the gauge field becomes diagonal the asymptotic holonomy can be calculated as where we have used (35): a(∞) = gva(∞) = gva ∞ . Fig. 10 implies that the asymptotic holonomy P ∞ depends on the electric charge q e through a ∞ : since we have fixed the magnetic charge q m to the unit q m = 4π/g. In the limit of vanishing electric charge q e → 0, the Yang-Mills dyon reduces to the Yang-Mills monopole and the asymptotic holonomy P ∞ becomes trivial P ∞ → 1 according to a ∞ → 0. In other words, the asymptotic holonomy P ∞ becomes nontrivial as long as the Yang-Mills dyon has a nonzero electric charge. Note that this is not the case of the self-dual dyon. Since the electric charge q e of the self-dual dyon is fixed (129), the asymptotic holonomy P ∞ of the self-dual dyon does not depend on the electric charge The KvBLL calorons are extensively used to reproduce the confinement/deconfinement phase transition in the Yang-Mills theory at finite temperature [13] by using the dyon gas model derived from the Yang-Mills theory where the interactions among dyons are calculated from quantum fluctuations around the dyon solution [19]. It should be remarked that the essential degrees of freedom responsible for confinement/deconfinement are not the KvBLL calorons themselves, but the constituent dyons characterized by the asymptotic behavior (139) with nontrivial holonomy at spatial infinity.
The (non-self-dual) Yang-Mills dyon obtained in this paper has also the nontrivial holonomy and therefore can be used to explain the confinement/deconfinement phase transition in the Yang-Mills theory at finite temperature, instead of using the KvBLL calorons or the self-dual dyons. In fact, it is possible to calculate the effective potential for the Polyakov loop in the same framework of the massive Yang-Mills model and show the existence of confinement/deconfinement transition at a certain critical temperature T c , which is obtained as a definite ratio to the gluon mass M , see [20]. The evaluation assumes a non-vanishing uniform background field for the component A 4 and takes into account the quantum fluctuations to one-loop order. This procedure is regarded as the first approximation for the non-uniform gauge field A 4 originating from the Yang-Mills dyon solution. Therefore, the existence of the Yang-Mills dyon with nontrivial holonomy (160) justify the the calculation of the effective potential of the Polyakov loop operator in a constant background [20] based on the same framework.
An advantage of using the Yang-Mills dyon is to give a successful explanation for quark confinement at zero temperature as the zero temperature limit of the finite temperature case. In the zero temperature limit T → 0, the Yang-Mills dyon reduces to the Yang-Mills magnetic monopole [4]. It has been already shown that such Yang-Mills magnetic monopoles successfully explain quark confinement at zero temperature from the viewpoint of dual superconductivity, see e.g. [3] for a review. In contrast, the KvBLL caloron reduces to the BPST instanton in the zero temperature limit T → 0. To the best of the author knowledge, however, the BPST instantons have not yet succeeded to explain quark confinement at zero temperature from the first principles without assuming additional inputs, see e.g., [21].

IX. CONCLUSION AND DISCUSSION
In this paper, we have constructed the dyon configurations in the pure SU (2) Yang-Mills theory both in the (3 + 1)dimensional Minkowski spacetime R 1,3 and in S 1 × R 3 space by incorporating a gauge-invariant gluon mass term   [23] contains the HS caloron in a minimal topological charge sector. Moreover, in the limit V → 0, the massive Yang-Mills theory reduces to the ordinary massless Yang-Mills theory, since the scalar field φ decouples.
even in the absence of the scalar field. Such a gauge-invariant mass term is obtained through a gauge-independent description of the BEH mechanism proposed in [5]. The procedure for obtaining the relevant dyon is guided by the "complementarity" between the SU (2) gauge-adjoint scalar model with a single radially fixed scalar field and the massive SU (2) Yang-Mills theory. In fact, we have obtained the static and spherically symmetric dyon configuration in the SU (2) massive Yang-Mills theory by solving the field equations of the "complementary" SU (2) gauge-adjoint scalar model with a single radially fixed scalar field. We have found that the static energy or the rest mass of the obtained Yang-Mills dyon is finite and proportional to the mass M X of the Yang-Mills gauge field A representing the existence of the massive component X .
In the long-distance region, we observed that the Yang-Mills dyon configuration A reduce to the restricted field V , which agrees with the dyonic extension of the Wu-Yang magnetic monopole as a consequence of the suppression of the massive modes X in the long-distance region. This feature is similar to the usual Julia-Zee dyons. In the short-distance region, on the other hand, the Wu-Yang magnetic monopole becomes singular, while the Julia-Zee dyon remains non-singular even at the origin. In the Yang-Mills dyon, we found that the massive components X play the very important role of canceling the singularity of V in the short-distance region such that the original gauge field A remains non-singular at the origin. This regularity of the Yang-Mills dyon is guaranteed by the logarithmic behavior of the gauge field itself without the aid of the scalar field, which vanishes at the origin as seen in Julia-Zee dyons. This behavior renders the energy of the Yang-Mills dyon finite even if the magnitude of the scalar field is fixed. It should be remarked that the chromomagnetic field B is divergent at the origin due to the logarithmic behavior of the solution f (ρ), which is, however, unessential for obtaining finite physical quantities such as energy, magnetic and electric charge density, and magnetic flux. Moreover, in the Yang-Mills dyon configuration, the time-component of the gauge field A 0 is regular, whose regularity is supported by the absence of the time-component of the high-energy massive mode: X 0 ≡ 0. Furthermore, we estimated the static mass of the Yang-Mills dyon by using the values of the previous studies [15,16]. We found that the heaviest static mass of the Yang-Mills dyon M dyon ≈ 1.18GeV is around the off-diagonal gluon mass M X ≈ 1.2GeV. This is a quite reasonable result for quark confinement to be realized due to condensation of the relevant Yang-Mills monopoles according to the dual superconductor picture. We need, however, more careful investigations to conclude whether or not the interactions among monopoles are indeed sufficient for realizing the monopole condensations, as examined by Polyakov [22] in the three dimensional case.
We observed that the Yang-Mills dyon cannot acquire the electric charge which is equal to the magnetic one. This is caused by a gauge-invariant mass term. In the contexts of instantons, the electric charge is equal to the magnetic one by definition, i.e., the (anti-)self-dual condition. However, there do not exist such (anti-)self-dual objects in our theory due to the mass term. We found that the Yang-Mills dyon in S 1 ×R 3 space has a nontrivial holonomy. This implies that our (non-self-dual) dyon with nontrivial holonomy P ∞ can be used to explain the confinement/deconfinement phase transition in the Yang-Mills theory at finite temperature based on the dual superconductor picture for confinement, instead of using the traditional KvBLL calorons or the self-dual dyon.
Finally, we give a conjecture that there will exist a caloron in the massive Yang-Mills theory so that the (non-selfdual) Yang-Mills dyon found in this paper could be identified with a constituent of the caloron. We call such a caloron the Yang-Mills caloron. In [23], the non-self-dual calorons which have the axial symmetry were constructed in the pure massless SU (2) Yang-Mills theory in the four-dimensional Euclidean space. Such axially symmetric solutions were also found in the Yang-Mills-Higgs model, i.e., the radially variable model by adopting the Kleihaus-Kunz ansatz [24]. We are therefore led to consider the axially symmetric Yang-Mills calorons in the massive Yang-Mills theory. In the zero temperature limit T → 0, the Yang-Mills dyon reduces to the Yang-Mills magnetic monopole [4]. This property is expected to hold in the Yang-Mills calorons, which will reduce to the Yang-Mills monopole-antimonopole chains. See Table I for the properties of the KvBLL caloron and the conjectured Yang-Mills caloron. This issue will be explored in near future.