arxiv: v4 [hep-th] 17 Aug 2017

Transcription

1 Microscopic structure of the QCD vacuum D.G. Pak, 1, 2, Bum-Hoon Lee, 2, Youngman Kim, 5 Takuya Tsukioka, 6 and P.M. Zhang 1 1 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 70000, China 2 Sogang University, Seoul , Korea Chern Institute of Mathematics, Nankai University, Tianjin 00071, China Asia Pacific Center of Theoretical Physics, Pohang, 790-0, Korea 5 Rare Isotope Science Project, Institute for Basic Science, Daejeon , Korea 6 School of Education, Bukkyo University, Kyoto , Japan arxiv: v [hep-th] 17 Aug 2017 We propose a new class of regular stationary axially symmetric solutions in a pure QCD which correspond to monopole-antimonopole pairs at macroscopic scale. The solutions represent vacuum field configurations which are locally stable against quantum gluon fluctuations in any small spacetime vicinity. This implies that the monopole-antimonopole pair can serve as a structural element in microscopic description of QCD vacuum formation through the monopole pair condensation. PACS numbers: q, 1.20.Dh, t, m Keywords: monopole, confinement, vacuum I. INTRODUCTION One of most attractive mechanisms of confinement is based on idea that QCD vacuum represents a dual superconductor which is formed due to condensation of color magnetic monopoles [1 5]. Realization of this idea in the framework of a rigorous theory has not been possible due to two principal difficulties: the first one, the color monopoles as physical particles have not been found in experiment, and so far there is no strict theoretical evidence of their existence in the real QCD. The second obstacle represents a long-standing problem of vacuum stability since late 70s when it was shown that Savvidy- Nielsen-Olesen vacuum is unstable [6, 7]. There have been numerous attempts to resolve this issue. Unfortunately, most of vacuum models lack the local vacuum stability at microscopic space-time scale. Besides, the known approaches to vacuum problem are based on assumption that vacuum field is static. Such an assumption is not consistent with the quantum mechanical principle which implies that elementary vacuum fields (vortices, monopoles etc.) should vibrate at microscopic level [8, 9]. In the present Letter we propose a novel approach for a microscopic description of QCD vacuum based on stationary vacuum field configurations. We construct a new class of regular stationary axially symmetric solutions which describe monopole-antimonopole pairs at macroscopic scale. The solutions represent vacuum fields which are locally stable against quantum gluon fluctuations in any small vicinity of each space-time point and can be served as structural elements of the QCD vacuum formed through the monopole pair condensation. We start with a standard Lagrangian of a pure SU(N) Yang-Mills theory and corresponding equations of motion L 0 = 1 F a µνf aµν, (D µ F µν ) a = 0, where (a = 1, 2,.., N), and (µ, ν = r, θ, ϕ, t) are the space-time indices in the spherical coordinates. Let us first consider stationary solutions in a pure SU(2) QCD. One can generalize a static axially symmetric Dashen- Hasslacher-Neveu (DHN) ansatz [10] as follows A 2 r = K 1 (r, θ, t), A 2 θ = K 2 (r, θ, t), A ϕ = K (r, θ, t), A 1 ϕ = K (r, θ, t), A 2 t = K 5 (r, θ, t). (2) The ansatz leads to a system of five differential equations which are invariant under residual U gauge transformations with a gauge parameter λ(r, θ, t) [11 1] K 1 = K 1 r λ, K 2 = K 2 + θ λ, K 5 = K 5 + t λ, K = K cos λ + K sin λ, K = K cos λ K sin λ. () To fix the residual symmetry we choose a Lorenz gauge by introducing the gauge fixing terms L g.f. = 1 2 ( rk r 2 θk 2 t K 5 ) 2. () We apply a method used in solving equations for the sphaleron solution [12, 1]. First, we decompose the functions K i in the Fourier series K i=1,2,, (r, θ, t) = δ i C + K 5 (r, θ, t) = n=1,2,... n=1,2,... i (nmr, θ) cos(nmt), 5 (nmr, θ) sin(nmt), (5) where C is an arbitrary number, M is a mass scale parameter which defines a class of conformally equivalent field configurations due to the presence of scaling invariance in a pure QCD, (r Mr, t Mt). The structure of the series expansion (5) implies that all time averaged color electric components of the field strength Fµν a vanish identically. Substituting the series decomposition truncated at a finite order n f into the classical action, and performing integration over the time period, one results in a reduced action which implies Euler equations for the

2 2 coefficient functions (n n f ) i (r, θ) [12, 1]. The obtained reduced equations for the modes (n n f ) i (r, θ) represent a well-defined system of elliptic partial differential equations which can be solved by applying standard numeric recipes. To solve the equations for (n n f ) i (r, θ) we impose the following Dirichlet boundary conditions i (r, θ) r=0 = 0, i (r, θ) θ=0,π = 0. (6) Asymptotic behavior at (r ) of the functions is determined by the equations of motion as follows i 1 a 1 (θ)sin(nmr) r 2, 5 a 5 (θ)cos(nmr) r 2, 2,, a 2,, (θ) cos(nmr), (7) where a i (θ) are periodic angle functions satisfying the finite energy density constraint. We choose the lowest angle modes for a i (θ), in particular, in the leading order one has a i 2 (θ) = c i sin 2 θ, a 2 (θ) = c 2 sin(2θ). With this one can solve numerically the equations up to the sixth order of series decomposition in the numeric domain (0 r L, 0 θ π). The obtained numeric solution (n=2k) implies that even order coefficient functions 1,2,,5 and (n=2k 1) odd order functions vanish identically. A typical solution in the leading and subleading order approximation is presented in Fig. 1. One has fast convergence for the modes i, and the solution profile functions of the third and fourth order, (n=,) i, provide corrections less than 2% by a norm with respect to the solution obtained in the leading and subleading order. The fifth and sixth order corrections are less than 0.1%. The energy density has an absolute maximum at the origin and a local maximum along the torus center line in the plane θ = π/2. For large values of the radial size L of the numeric domain, the energy 1 density decreases as. The solution provides an absolute maximum of the classical action defined in the L2 whole infinite space-time. This implies classical stability of the solution in the space of general axially symmetric fields. In the leading order one has non-vanishing 1 time averaged color magnetic fields F F 1 θϕ 2 and rϕ 1, which create magnetic fluxes corresponding to a pair of non-topological monopole and antimonopole located at the origin. Non-topological monopoles have non-integer magnetic charges which are not gauge invariant quantities. The color magnetic flux through a sphere surrounding the monopole is not conserved and depends on a radius of the sphere. Note that we have not found a regular stationary solution describing a pair of separated Wu-Yang monopole and antimonopole [1]. It is unexpected that a regular monopole-antimonopole pair solution in the limit of zero distance between the monopoles does exist in the (c) (e) FIG. 1: Solution profile functions in the leading and subleading order: 1 ; 2 ; (c) ; (d) 5 ; (e) ; (f) the energy density plot (C = 2, g = 1, M = 1). (2) real QCD. A non-trivial singular field configuration describing a pair of topological monopole and antimonopole located at one point was known before [15]. Such a field is defined by the gauge potential A a µ = δ ϕ µ δ a sin 2 θ corresponding to the mapping S 2 S 2 with double covering the sphere which belongs to a trivial homotopy class. An example of a static regular monopole-antimonopole solution of this type is given in a special CP 1 model [15]. Let us now consider stationary solutions in SU() QCD. We are looking for essentially SU() field configurations which do not reduce to embedded SU(2) solutions. An axially symmetric ansatz contains the following non-vanishing components of the gauge potential corresponding to three I, U, V -type SU(2) subgroups of SU() (d) A 2 1 = K 1, A 2 2 = K 2, A 1 = K, A 2 = K 5, A 5 1 = Q 1, A 5 2 = Q 2, A = Q, A 5 = Q 5, A 7 1 = S 1, A 7 2 = S 2, A 6 = S, A 7 = S 5, A = K, A 8 = K 8, (8) where the fields K i, Q i, S i depend on space-time coordinates (r, θ, t). We employ the same method as in the case of SU(2) QCD using series decomposition for K i, Q i, S i (f)

3 as in (5) (the Abelian potentials K,8 have similar decomposition). After integration over time in the classical action one can derive the Euler equations for the field modes i, Q i, S i. It is remarkable that the obtained system of Euler equations can be simplified further by applying the following reduction ansatz Q 1,2,5 = S 1,2,5 = 1,2,5,,8 = 2 Q = ( = ( S ),, ), (9) where (n = 1,, 5,...) and all even modes vanish. The ansatz reduces the number of equations to four equations for independent functions in each order. Note 1,2,,5 that the Abelian potentials,8 are equal to each other as it takes place in the case of Abelian Weyl symmetric homogeneous magnetic fields providing an absolute minimum of the quantum effective potential [16, 17]. We demonstrate existence of another type of monopole pair solution which is different from the SU(2) stationary monopole pair considered above. To find such a solution i we apply asymptotic conditions (7) for the functions containing even angle modes a 2, (θ) and odd angle modes a 1,5 (θ). Imposing vanishing Dirichlet boundary conditions at the origin and along the boundaries (θ = 0, π), we solve the equations for i in the fifth order approximation. The obtained solution profile functions i and energy density plots are presented in Fig. 2. Consider the magnetic flux structure of the obtained solution when the observation time is much larger than the period of time oscillations. The time averaged radial components of the Abelian field strength in the leading order include only non-linear terms F θϕ = 2, 8 F θϕ = + 2, (10) where the Abelian potential contains the angle dependent factor sin 2 θ, and 2 is an even function with respect to the reflection symmetry θ π θ. The,8 magnetic fields F θϕ create opposite magnetic fluxes through a sphere of radius R with a center at the origin which correspond to non-topological antimonopole and monopole located at one point. Careful numeric analysis shows that stationary solutions corresponding to minimal energy in a chosen finite numeric domain are classified by only two parameters, the conformal parameter M and amplitude c 0 of the oscillating Abelian potential. A detailed structure of the numeric solution up to the fifth order series decomposition is given in [1]. The most important property of the proposed stationary solutions is that they possess quantum stability against small vacuum gluon fluctuations. To verify (c) (e) FIG. 2: Profile functions for the monopole pair solution in the leading order: 1 ; 2 ; (c) ; (d) 5 ; (e) the energy density plot; (f) the energy density contour plot in Cartesian coordinates (g = 1, M = 1). the stability of the solutions it is suitable to apply the quantum effective action formalism. Consider one-loop quantum effective action expressed in terms of functional operators [18] S 1loop eff = 1 Tr ln[kab 2 µν] + Tr ln[mfp], ab K ab µν = δ ab δ µν 2 t δ µν (D ρ D ρ ) ab 2f acb F c µν, (11) M ab FP = (D ρ D ρ ) ab, where the Wick rotation t it has been performed to provide the causal structure of the operators, and D µ, Fµν c are defined with a classical background field Bµ. a The operators Kµν ab and MFP ab correspond to gluon and Faddeev-Popov ghosts. Note that the expression (11) is valid for arbitrary background field and does not depend on a chosen gauge for the background field due to use of a gauge covariant background formalism [19]. Effective action describes the vacuum-vacuum amplitude, and the presence of an imaginary part of the action implies vacuum instability. Therefore, if the operator Kµν ab is not positively defined then an unstable mode will appear as (d) (f)

4 an eigenfunction corresponding to a negative eigenvalue of the following Schrödinger equation K ab µ Ψ b ν = λψ a µ, (12) where the wave functions Ψ a µ(t, r, θ, ϕ) describe gluon fluctuations. Note that the ghost operator MFP ab is positively defined and does not produce instability [7]. The potential in the operator Kµ ab does not depend on the asimuthal angle. Due to this one can separate a corresponding angle dependent part from the function Ψ a µ and solve the eigenvalue equation in a three-dimensional domain (t, r, θ). Substituting interpolation functions for the stationary magnetic solutions in the leading order, one can solve the eigenvalue equation (12). In the case of SU(2) stationary solution a full eigenvalue spectrum is divided into four sub-spectra corresponding to four decoupled systems of equations: (I) Ψ 2, (II) Ψ 1, Ψ, (III) Ψ 2 1, Ψ 2 2, Ψ 1, Ψ, (IV) Ψ 1 1, Ψ 1, Ψ 1 2, Ψ 2, Ψ 2. The lowest eigenvalue is positive, and it is reached by a solution satisfying the system of equations (II), the corresponding eigenfunctions are plotted in Fig. (g = 1, M = 1). Other systems of equations, (I,III,IV), have a similar FIG. : Non-vanishing eigenfunctions corresponding to the lowest eigenvalue λ = : Ψ 1 ; Ψ. structure of the eigenvalue spectrum. In the case of SU() stationary monopole pair solution the equations in (12) are not factorized, and one has to solve a full set of thirty two differential equations. Numerical results of solving the Schrödinger eigenvalue equation with the stationary monopole-antimonopole background field show that the eigenvalue spectrum is positively defined. The obtained dependence of the lowest eigenvalue on the size of the chosen numeric domain for large values of L confirms the positiveness of the eigenvalue spectrum, Fig.. This proves the quantum stability of stationary monopole-antimonopole pair solutions in QCD. The quantum stability of SU(2) and SU() stationary background fields has been checked for solutions with amplitude values of the Abelian potential K (r, θ, t) in the interval (0 c 0 2) and with conformal parameter values (0 M 2). For large values of M, c 0 unstable modes appears which destabilize the vacuum. Let us consider symmetry properties of essentially SU() stationary solutions. The classical Yang-Mills Lagrangian can be rewritten in terms of Weyl symmetric Λ FIG. : Dependence of the ground state eigenvalue λ(l) on the size L of the numeric domain in the cases of SU() (red curve) and SU(2) (blue curve) stationary monopoleantimonopole pair solutions (g = 1, M = 1). fields [17] L 0 = p=1,2, { 1 6 (G p µν) Dp µw p ν D p νw p µ 2 igg p µνw p µ W p ν } L () int [W ], (1) where Gµν p are Abelian field strengths containing the gauge potentials A,8 µ, the complex fields Wµ p represent off-diagonal gluons, and the index p counts the Weyl symmetric gauge potentials. Note that the Lagrangian is not Weyl symmetric under permutation of I, U, V subgroups SU(2) since the quartic interaction term L () int [W ] is not factorized into a sum of separate parts corresponding to I, U, V sectors. It is remarkable that the Lagrangian on the space of essentially SU() stationary solutions possesses high symmetric structure. First of all, substituting the reduction ansatz (9) for general functions K i, Q i, S i into the Lagrangian, one can verify that L () int obtains an explicit Weyl symmetric form L () int [W ] = 9 8 p=1,2, ( (W p µ W p µ) 2 (W µ) 2 (W ν ) 2) (1) Secondly, each I, U, V sector in the Lagrangian (1) contains cubic interaction terms corresponding to the anomaly magnetic moment interaction which is precisely the source of the Nielsen-Olesen vacuum instability [7]. It is surprising, within the framework of the ansatz (9) one has complete mutual cancellation of all cubic interaction terms. This implies that on the space of Weyl symmetric fields the classical action describes a generalized λφ theory. In conclusion, simple consideration shows that our approach to QCD vacuum problem based on stationary Weyl symmetric monopole pair solutions opens a new perspective towards construction of a microscopic theory of vacuum and vacuum phase transitions. An important feature of the monopole pair solution is that it represents

5 5 a non-trivial essentially non-abelian field configuration which describes a pair of monopole and antimonopole located at one point. This implies that monopole and antimonopole, as well as monopole-antimonopole pairs with opposite color orientations, merge into a stable state with a finite energy density in the limit of zero distance between monopole and antimonopole. This is contrary to the case of Dirac and Wu-Yang pair of monopole and antimonopole which annihilate when they meet each other. We expect that QCD vacuum is formed due to condensation of monopole-antimonopole pairs, and the microscopic vacuum structure is characterized by few parameters: the conformal parameter M, the amplitude of oscillations c 0 of the Abelian gauge potential K in the asymptotic region, and the concentration of monopole pairs at zero temperature. Numeric analysis shows that with increasing temperature the internal energy of each monopole pair increases and at some critical values of the parameters (M > 2, c 0 > 2) the monopole-antimonopole pair becomes unstable. Note that in the confinement phase the vacuum averaging value of the gluon field operator 0 A a µ 0 vanishes since the size of the hadron is much larger than the characteristic length λ M = 2π/M of the vacuum monopole field oscillations. To describe dynamics of the vacuum structure at finite temperature one should apply the Euclidean functional integral formalism with time integration in the finite interval (0 t β = 1/kT ). It is clear that at high enough temperature the upper integration limit β will be less than λ M. This will lead to a non-vanishing vacuum averaging value of the gluon field operator, 0 A a µ 0, and transition to the deconfinement phase with spontaneous symmetry breaking where the gluon can be observed as a color object. Acknowledgments One of authors (DGP) thanks Prof. C.M. Bai for warm hospitality during his staying in Chern Institute of Mathematics and Dr. Ed. Tsoy for useful discussions of numeric aspects. The work is supported by NRF of Korea, grant MSIP No.201R1A2A1A (ERND), and by Brain Pool Program (KOFST). [1] S.J. Brodsky, G.F. de Teramond, and H.G. Dosch, Int. J. Mod. Phys. A29, 101 (201). [2] Y. Nambu, Phys. Rev. D10, 262 (197). [] S. Mandelstam, Phys. Rep. 2C, 25 (1976). [] A. Polyakov, Nucl. Phys. B120, 29 (1977). [5] G. t Hooft, Nucl. Phys. B190, 55 (1981). [6] G.K. Savvidy, Phys. Lett. B71, 1 (1977). [7] N.K. Nielsen and P. Olesen, Nucl. Phys. B1, 76 (1978). [8] N.K. Nielsen and P. Olesen, Nucl. Phys. B160, 80 (1979). [9] J. Ambjørn and P. Olesen, Nucl. Phys. B170, 60 (1980). [10] R. Dashen, B. Hasslacher and A. Neveu, Phys. Rev. D10 (197) 18. [11] N.S. Manton, Nucl. Phys. B15, 19 (1978). [12] C.Rebbi and P. Rossi, Phys. Rev. D 22, 2010 (1980). [1] J. Kunz, B. Kleihaus and Y. Brihaye, Phys. Rev. D 6, 587 (1992). [1] D.G. Pak, B.-H. Lee, Y. Kim, T. Tsukioka and P.M. Zhang, Microscopic structure of the QCD vacuum - Supplemental Material. [15] L.P. Zou, P.M. Zhang, D.G. Pak, Phys. Lett. B 728C, 650 (201). [16] H. Flyvbjerg, Nucl. Phys. B176, 79 (1980). [17] Y.M. Cho, J.H. Kim, D.G. Pak, Mod. Phys. Lett. A21, 2789 (2006). [18] Y.M. Cho and D.G. Pak, Phys. Rev. D65, (2002). [19] L.F. Abbott, Acta Phys. Polon. B1, (1982).