arxiv: v1 [hep-th] 10 Oct 2017

Similar documents
arxiv: v3 [hep-th] 11 Nov 2018

The integral equations of Yang-Mills and its gauge invariant conserved charges

arxiv: v2 [hep-th] 28 Nov 2018

Lecturer: Bengt E W Nilsson

Contact interactions in string theory and a reformulation of QED

Chern-Simons Theory and Its Applications. The 10 th Summer Institute for Theoretical Physics Ki-Myeong Lee

Introduction and Vectors Lecture 1

UNIVERSITI TEKNOLOGI MARA EXACT SOLUTIONS OF THE SU(2) YANG-MILLS-HIGGS THEORY

Kink Chains from Instantons on a Torus

Lecturer: Bengt E W Nilsson

Singular Monopoles and Instantons on Curved Backgrounds

arxiv: v2 [hep-th] 28 Aug 2014

Lecture 9: RR-sector and D-branes

8.324 Relativistic Quantum Field Theory II MIT OpenCourseWare Lecture Notes Hong Liu, Fall 2010 Lecture 3

Spinning strings and QED

arxiv: v2 [hep-th] 21 Aug 2009

Maxwell s equations for electrostatics

752 Final. April 16, Fadeev Popov Ghosts and Non-Abelian Gauge Fields. Tim Wendler BYU Physics and Astronomy. The standard model Lagrangian

arxiv:hep-th/ v2 14 Oct 1997

Interacting non-bps black holes

Quantum Nambu Geometry in String Theory

Coordinate/Field Duality in Gauge Theories: Emergence of Matrix Coordinates

Monopoles and Skyrmions. Ya. Shnir. Quarks-2010

BPS Solitons and Killing Spinors in Three Dimensional N =2Supergravity

Electric Dipole Moment of Magnetic Monopole

Some Quantum Aspects of D=3 Space-Time Massive Gravity.

Physics 557 Lecture 5

A Brief Introduction to AdS/CFT Correspondence

Instantons and Monopoles in Maximal Abelian Projection of SU(2) Gluodynamics

Fourth Aegean Summer School: Black Holes Mytilene, Island of Lesvos September 18, 2007

Lecture 24 Seiberg Witten Theory III

Curved spacetime and general covariance

Note 1: Some Fundamental Mathematical Properties of the Tetrad.

Dynamics of Multiple Kaluza-Klein Monopoles in M- and String Theory

Lecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)

Physics 236a assignment, Week 2:

Magnetostatics. Lecture 23: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay

Baby Skyrmions in AdS 3 and Extensions to (3 + 1) Dimensions

Continuity of the Deconfinement Transition in (Super) Yang Mills Theory

Introduction to Group Theory

Physics/Astronomy 226, Problem set 4, Due 2/10 Solutions. Solution: Our vectors consist of components and basis vectors:

PoS(LAT2005)324. D-branes and Topological Charge in QCD. H. B. Thacker University of Virginia

Two-monopole systems and the formation of non-abelian clouds

Half BPS solutions in type IIB and M-theory

Infrared Propagators and Confinement: a Perspective from Lattice Simulations

Abelian and non-abelian Hopfions in all odd dimensions

A Multimonopole Solution in String Theory

8.821 F2008 Lecture 18: Wilson Loops

arxiv: v2 [hep-th] 9 Feb 2012

Topological Solitons from Geometry

August 2008 Solitons and AdS string solutions 2

arxiv:hep-th/ v1 7 Nov 1998

Continuity of the Deconfinement Transition in (Super) Yang Mills Theory

Chapter 7 Curved Spacetime and General Covariance

Electrically and Magnetically Charged Solitons in Gauge Field Th

Solitons, Vortices, Instantons

free space (vacuum) permittivity [ F/m]

Lecture 4-1 Physics 219 Question 1 Aug Where (if any) is the net electric field due to the following two charges equal to zero?

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015

Introduction to Chern-Simons forms in Physics - I

The Role Of Magnetic Monopoles In Quark Confinement (Field Decomposition Approach)

Non-Abelian Berry phase and topological spin-currents

Lorentz-covariant spectrum of single-particle states and their field theory Physics 230A, Spring 2007, Hitoshi Murayama

The stability of the wall is due to the potential and the gradient balancing each other (Derrick s Theorem). Simple argument: = E grad

Higher-Spin Black Holes and Generalised FZZ Duality

SUSY Breaking in Gauge Theories

arxiv: v2 [hep-th] 26 Jan 2018

Instantons in string theory via F-theory

arxiv:hep-th/ v2 19 Sep 2001

Spectrum of Holographic Wilson Loops

HIGHER SPIN CORRECTIONS TO ENTANGLEMENT ENTROPY

Two-loop Remainder Functions in N = 4 SYM

Quantum Fields in Curved Spacetime

Tachyon Condensation in String Theory and Field Theory

All symmetric AdS n>2 solutions of type II supergravity

Symmetries and solutions of field equations of axion electrodynamics

t Hooft Loops and S-Duality

NTNU Trondheim, Institutt for fysikk

Relativistic Mechanics

Chern-Simons gauge theory The Chern-Simons (CS) gauge theory in three dimensions is defined by the action,

B r Solved Problems Magnetic Field of a Straight Wire

Newman-Penrose formalism in higher dimensions

INSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD

Overthrows a basic assumption of classical physics - that lengths and time intervals are absolute quantities, i.e., the same for all observes.

Dynamical compactification from higher dimensional de Sitter space

Vectors. January 13, 2013

Electromagnetic. G. A. Krafft Jefferson Lab Jefferson Lab Professor of Physics Old Dominion University TAADI Electromagnetic Theory

YANG-MILLS GAUGE INVARIANT THEORY FOR SPACE CURVED ELECTROMAGNETIC FIELD. Algirdas Antano Maknickas 1. September 3, 2014

arxiv: v2 [hep-th] 2 Apr 2015

Regularization Physics 230A, Spring 2007, Hitoshi Murayama

Cosmic Strings and Topological Defects

Dyon degeneracies from Mathieu moonshine

Integrable Spin Systems From Four Dimensions

QM and Angular Momentum

Donoghue, Golowich, Holstein Chapter 4, 6

Physics 342 Lecture 30. Solids. Lecture 30. Physics 342 Quantum Mechanics I

Simple Evaluation of the Chiral Jacobian with the Overlap Dirac Operator

Gauge Theory of Gravitation: Electro-Gravity Mixing

Show, for infinitesimal variations of nonabelian Yang Mills gauge fields:

Vacuum Polarization in the Presence of Magnetic Flux at Finite Temperature in the Cosmic String Background

Transcription:

A direct test of the integral Yang-Mills equations through SU) monopoles C. P. Constantinidis, L. A. Ferreira and G. Luchini arxiv:171.3359v1 hep-th] 1 Oct 17 Departamento de Física Universidade Federal do Espírito Santo UFES), CEP 975-9, Vitória-ES, Brazil Instituto de Física de São Carlos; IFSC/USP; Universidade de São Paulo Caixa Postal 369, CEP 1356-97, São Carlos-SP, Brazil Abstract We use the SU) t Hooft-Polyakov monopole configuration, and its BPS version, to test the integral equations of the Yang-Mills theory. Those integral equations involve two complex) parameters which do not appear in the differential Yang-Mills equations, and if they are considered to be arbitrary it then implies that non-abelian gauge theories but not abelian ones) possess an infinity of integral equations. For static monopole configurations only one of those parameters is relevant. We expand the integral Yang-Mills equation in a power series of that parameter and show that the t Hooft-Polyakov monopole and its BPS version satisfy the integral equations obtained in first and second order of that expansion. Our results points to the importance of exploring the physical consequences of such an infinity of integral equations on the global properties of the Yang-Mills theory.

1 Introduction The purpose of this paper is to perform a test of the integral equations of Yang-Mills theories, recently proposed in 1, ], using the SU) t Hooft-Polyakov monopole solution 3, 4] as well as its exact analytical BPS version 5, 6]. The main motivation for such a test is that these integral equations involve two complex parameters that are not present in the Yang-Mills partial differential equations. If those parameters are arbitrary, it means that contrary to abelian electromagnetism, Yang-Mills theories possess in fact an infinity of integral equations. Indeed, by expanding the Yang-Mills integral equations in power series of those parameters, we check that the SU) t Hooft-Polyakov monopole, and its BPS version, do satisfy the integral equations appearing in that expansion, up to second order in one of the parameters. The cancelations involved in such a check are highly non-trivial and give strong evidence on the arbitrariness of those parameters. As shown in 1, ] the integral Yang-Mills equations lead in a quite natural way to gauge invariant conserved charges. Such charges involve those two parameters in a way that if they are indeed arbitrary it would imply that in principle, the number of charges is infinite. However, due to some special properties of BPS multi-dyon solutions 7, 8], shown in 9], the higher charges are not really independent for such solutions, being in fact powers of the first ones the electric and magnetic charges). The same is true for the SU) t Hooft-Polyakov monopole. There remains to be investigated if other non-bps solutions also present such special properties or not, and so possess or not an infinity of charges. In order to discuss the role of such parameters in a more concrete way let us start by the theory of electromagnetism described by the Maxwell equations µ f µν = j ν µ f µν = 1.1) where f µν = µ a ν ν a µ, f µν = 1 ε µνρλf ρλ, j µ being the external four current, and a µ the electromagnetic four vector potential. The integral version of those equations is obtained through the abelian Stokes theorem for a rank-two antisymmetric tensor b µν on a space-time 3-volume Ω, as Ω b = Ω d b, where Ω is the border of Ω. Taking b µν as a linear combination of f µν and its Hodge dual, and using 1.1), one gets ˆ α f µν + β f ] ˆ µν dσ µν = β j µνρ dv µνρ 1.) Ω where j µνρ = ε µνρλ j λ is the Hodge dual of the external current and α and β are arbitrary parameters used in the liner combination. By considering α and β to be arbitrary, the integral equations 1.) correspond to the four usual integral equations of electromagnetic theory, which in fact preceded Maxwell differential equations. Indeed, taking α = and Ω to be a purely spatial 3-volume one gets the Gauss law. On the hand, taking β = and Ω to be a solid cylinder with its height in the time direction, and its base on a spatial plane, one gets the Faraday law, and so on. The role of the parameters α and β are not really important here because 1.) is linear in them. The situation becomes more complex in a non-abelian gauge theory. 1 Ω

The Yang-Mills theories were formulated à la Maxwell in terms of partial differential equations, the so-called Yang-Mills equations1] D µ F µν = J ν D µ F µν = 1.3) where F µν = µ A ν ν A µ +i e A µ, A ν ], with e being the gauge coupling constant, F µν = 1 ε µνρλf ρλ, J µ being the external matter current, and D µ = µ +i e A µ, ], and A µ being the non-abelian gauge field taking value on the Lie algebra of the gauge group G. In order to construct the integral form of Yang-Mills equations 1.3) one needs the non-abelian version of the Stokes theorem for a non-abelian) rank-two antisymmetric tensor B µν on a space-time 3-volume Ω. Even though the non-abelian Stokes theorem for a one-form connection on a -surface was known for some time, the same theorem for a two-form connection was constructed only more recently in 11, 1] using concepts on generalized loop spaces. Conceptually everything becomes more clear if one uses the two-form B µν, defined on space-time, to construct a one-form connection on the generalized loop space. Using such generalized non-abelian Stokes theorem, the integral form of Yang-Mills equations were constructed in 1, ]. The formulas involve path, surface and volume ordered integrals as follows. Consider a space-time 3-volume Ω, and choose a reference point x R on its border Ω. Scan Ω with closed -surfaces based on x R, labelled by a variable ζ, such that ζ = corresponds to the infinitesimal surface around x R, and ζ = ζ to the border Ω. Then scan each closed -surface with loops, starting and ending at x R, labelled by a variable τ. Each loop is parameterized by a variable σ. The integral form of Yang-Mills equations 1.3) is given by 1, ] V Ω) P e ie Ω dσw αf µν+β F µν)w xµ x ν σ τ = P 3 e Ω dζv J V U Ω), 1.4) where P and P 3 mean surface and volume ordered integration respectively, as explained above, and J = ˆ σf + e ˆ σ { deβ J µνλ W dx µ dx ν dσ dσ d xκ d x µ d σ d σ dx λ dζ α 1) F W κρ + β F W κρ ) σ ), d x ρ σ ) d τ d x ν σ) d ζ d xρ σ ) d ζ αfµν W + β F ) ] µν W σ) d x ν )} σ) d τ with J µνλ = ε µνλρ J ρ, being the Hodge dual of the external matter current. In order to simplify the formulas we have used the notation 1.5) X W W X W 1.6) with X standing for the field tensor, its Hodge dual, or the dual of the matter currents. The quantity W appearing above stands for the Wilson line, defined on a path parameterized by σ through the equation dw dσ + iea dx µ µ dσ W = 1.7)

and so W = 1 i e ˆ σ dσ A µ σ ) dx µ dσ + i e) ˆ σ dσ A µ σ ) dx µ dσ ˆ σ dσ A ν σ ) dx ν... 1.8) dσ The quantity V, called the Wilson surface, is defined on a surface parameterized by σ and τ, through the equation with dv V T τ) =, 1.9) ˆ σf T τ) = ie dσ W αf µν + β F ) µν W xµ x ν σ τ. 1.1) and the integration being on the closed loops used in the scanning of Ω, as explained above. The initial and final values of σ, denoted and σ f respectively, correspond to the initial and final points of the loop, which in fact are the same point since the loop is always closed. Therefore the solution of 1.9) is the surface-ordered series ˆ τ ˆ τ V τ) = 1 + T τ ) + ˆ τ T τ )T τ ) +... 1.11) The l.h.s. of 1.4) is obtained by integrating 1.9) on the -surface Ω, i.e. the border of Ω. On the other hand the r.h.s. of 1.4) is obtained by integrating the equation on the 3-volume Ω, and where du dζ K U = 1.1) ˆ K = Σ V J V 1.13) with Σ being the closed -surfaces scanning Ω, labelled by ζ, and J given by 1.5). The solution of 1.1) is given by the volume-ordered series ˆ ζ ˆ ζ ˆ ζ Uζ) = 1 + dζ Kζ ) + dζ dζ Kζ ) Kζ ) +... 1.14) Note that 1.4) does reduce to 1.) in the case that the gauge group G is U1). However, for non-abelian gauge groups the dependence of both sides of 1.4) on the parameters α and β are highly non-linear. Indeed, if such parameters are arbitrary one can expand both sides of 1.4) in a power series on them. The coefficient of each term of such series on the l.h.s. of 1.4) will have to equal the corresponding coefficient of the series on the r.h.s., leading to an infinity of integral equations. Consequently any solution of the Yang-Mills equations 1.3) will have to satisfy such an infinity of integral equations. It is this test that we want to perform with the t Hooft-Polyakov monopole, and its exact analytical BPS version 8]. We shall consider the 3-volume Ω to be purely spatial, and consequently only the spatial components of the field tensor and its dual, i.e. F ij and F ij, i, j = 1,, 3, will be present on both sides of 1.4). However, F ij is proportional to the electric field and so it vanishes for those static monopole solutions. In addition, only the component J 13 J 3

appears on the r.h.s. of 1.4), and that vanishes because the solution is static and we shall work in the gauge where A =. Remember that the only contribution for the matter current for such a solution comes from the triplet Higgs field φ, and that is of the form J µ φ, D µ φ ]. Therefore, all terms involving the parameter β are not present on both sides of 1.4), for static monopoles when Ω is purely spatial, and so it reduces to with J ˆ = e ˆ σf P e i e α Ω dσw F ij W xi x j σ τ = P 3 e αα) Ω dζv J ˆ V, 1.15) dσ F W k,l σ ), ˆ σf ] x dσfij W k x i ) x l σ) σ σ τ σ ) xj xl σ) ζ ζ σ ) xj τ σ) 1.16) where i, j, k, l = 1,, 3, repeated indices are summed, and where we have denoted Fij W W F ij W. Note that we have explored the symmetry of J ˆ in σ and σ to replace σ f dσ σ dσ 1 σf dσ σ f dσ. Equation 1.15) is what we call generalized integral Bianchi identity. Note that one would expect the integral Bianchi identity to be 1.15) for α = 1, i.e. P e i e Ω dσw F ij W xi x j σ τ = 1l 1.17) Indeed, that is what leads to the quantization of the magnetic charge. From a physical point of view it is intriguing that by rescaling the field tensor magnetic field) as F ij α F ij, leads to the appearance of a term like the r.h.s of 1.15), making the magnetic flux through Ω to change drastically. However, the validity of 1.4), and so of 1.15), is guaranteed by the generalized nonabelian Stokes theorem for a two-form B µν and the partial differential Yang-Mills equations 1.3) as proved in 11, 1, 1, ]. The intriguing non-linear phenomenon that we want to directly check in this paper, is if one can expand both sides of 1.15) in powers of α, and if the SU) t Hoot-Polyakov monopole and its exact analytical BPS version, satisfy each one of the integral equations obtained through such an expansion. The paper is organized as follows. In section we perform the expansion of the generalized integral Bianchi identity 1.15) in powers of the parameter α, and we show that each term of the expansion can be expressed solely in terms of the Wilson line operator. In section 3 we calculate explicitly the Wilson line operator for the SU) t Hooft-Polyakov and BPS monopoles using a suitable scanning of surfaces and volumes. The result is quite simple and it is given in 3.14). In section 4 we check the validity of the integral equation in first order of the α-expansion, and in section 5 we do the same for the integral equation in second order of that same expansion. We present our conclusions in section 6 and in the appendix A we give the results of the numerical calculations of the integrals needed to perform the check of the integral equations. The expansion of the Yang-Mills integral equations Assuming that α and β are indeed arbitrary we expand both sides of 1.4) in power series in those parameters. As we have said the l.h.s. of 1.4) is obtained by integrating 1.9), and its r.h.s. 4

by integrating 1.1). By writing the expressions on the l.h.s and on the r.h.s of the integral equation 1.4) in terms of 1.11) and 1.14), and collecting the coefficients at first order in α and zeroth order in β, we get the integral equation at first order in α ˆ σf dσ F W µν x µ σ x ν ˆ ζ τ ζ=ζ = ie dζ d xκ d x µ d σ d σ ˆ σf d x ρ σ ) d τ ˆ σ dσ dσ F W κρ σ ), Fµν W σ) ] d x ν σ) d ζ d xρ σ ) d ζ d x ν ) σ) d τ where ζ is the value of ζ corresponding to the closed surface Ω, in the scanning of the 3-volume Ω, which is the border of Ω see explanation of the scanning in the paragraph above 1.4)). On the other hand, the integral equation appearing in order β and zeroth order in α, in the expansion of 1.4) is given by ˆ σf dσ F W µν x µ σ x ν τ ζ=ζ = ˆ ζ + i e dζ ˆ σ dσ d xκ d x µ d σ d σ ˆ σf { dσ J W dx µ dx ν dx λ µνλ dσ dζ Fκρ W σ ), F ] µν W σ) d x ρ σ ) d τ d x ν σ) d ζ d xρ σ ) d ζ d x ν )} σ) d τ Note that in the case where the gauge group G is the abelian group U1), the equation.1) corresponds to 1.) for α = 1 and β =. Equation.) corresponds to 1.) for α = and β = 1. Note in addition that in the case where the 3-volume Ω is purely spatial, the commutator term in.1) involving the field tensors can be interpreted as a density of non-abelian magnetic charge associated to the gauge field configuration inside Ω. The commutator term in.) involving the field tensor and its Hodge dual can be interpreted as a density of non-abelian electric charge associated to the gauge field configuration inside Ω. In the case where Ω has time components, those commutators will be associated to flows of non-abelian electric and magnetic charges. We have explored further those facts to obtain the integral form of the non-abelian Gauss, Faraday, etc., laws, and the physical implications of these new terms the commutator terms) should be further explored in some other opportunity. As one goes higher in the expansion, the integral equations become more and more complex. However, for the case we are considering in this paper, namely the static t Hooft-Polyakov monopole and its BPS version, and where the 3-volume Ω is purely spatial, there is an important simplification. As we have argued in the paragraph above 1.15), only the spatial components of the field tensor magnetic field) appear in the formulas, since the spatial components of its Hodge dual electric field) vanish. As explained in section of 11], or in the appendix of ], if one performs an infinitesimal variation, x µ σ) x µ σ) + δx µ σ), of the curve where the Wilson line W 1.7) is calculated, but keeping its end points fixed, the infinitesimal variation of the Wilson line operator is given by.1).) W σ f ) δw σ f ) = ie ˆ σf dσ W F µν W d xµ d σ δxν.3) 5

The Wilson line operators W appearing in the Yang-Mills integral equations 1.4) are evaluated on the paths that scan the closed surfaces which in their turn scan the 3-volume Ω. Thefore, as we vary the parameter τ which labels the loops, we vary the loop along a given surface, and so δx µ = dxµ δτ. When we vary the parameter ζ which labels the surfaces, the loops vary perpendicular to that surface and so δx µ = dxµ dζ δζ. Consequently, from.3) we get the following two useful formulas ˆ σf ie ie ˆ σf dσ W F µν W dxµ dσ dσ W F µν W dxµ dσ dx ν dx ν dζ = W dw = W dw dζ..4) As we have shown, for the static t Hooft-Polyakov monopole and its BPS version, and a purely spatial 3-volume Ω, the integral Yang-Mills equation 1.4) becomes the generalized integral Bianchi identity 1.15). Therefore, from 1.11), 1.1) and.4), one gets that the l.h.s. of 1.15) is given by V Ω) P e i e α Ω dσw F ij W xi x j σ τ = 1 + α + α W dw τ) ζ=ζ.5) ˆ τ = 1 + α V 1) + α V ) +... dw W τ )W dw τ) ζ=ζ +... From.4) one gets that 1.16) becomes J ˆ = W dw, W dw dζ ]..6) Therefore, from 1.14) and.6), the r.h.s. of 1.15) becomes U Ω) P 3 e αα) ˆ Ω dζv J ˆ ζ V = 1 α α 1) dζ ˆ ζ ˆ ζ + α α 1)] dζ dζ V V dw dw W, W dζ = 1 + α U 1) + α U ) +... V W dw, W dw dζ W dw, W dw dζ ] V ] ) V τ, ζ) ] ) τ V, ζ ) +....7) where V in.7) is evaluated with the same expansion as in 1.11) with β =, and so an expansion similar to.5). Therefore, by equating.5) to.7) one gets that the term in first order in α leads to the integral equation V 1) = W dw ζ=ζ = ˆ ζ dζ W dw, W dw dζ ] = U 1)..8) 6

Similarly, the term in order α gives the following integral equation V ) = + + ˆ ζ ˆ ζ dζ dζ ˆ τ ˆ ζ dw W ˆ τ dζ W dw ζ=ζ = dw W, ˆ ζ dζ W dw, W dw dζ W dw, W dw dζ ] ] W dw, W dw dζ ] dw dw W, W dζ ].9) ] = U ). We are going to verify if the SU) t Hooft-Polyakov monopole and its BPS version 3, 4, 5, 6, 8], satisfy the integral equations.8) and.9). Note that the only quantity appearing in.8) and.9) is the Wilson loop W. In the next section we evaluate it for those monopole solutions. 3 The Wilson loop for t Hooft-Polyakov and BPS monopoles reads The spherically symmetric t Hooft-Polyakov ansatz 3, 4] for a SU) static magnetic monopole φ = 1 H ζ) ˆr T e r A = 3.1) A i = 1 e ɛ x j ijk r 1 Kζ)) T k with r = x 1 + x + x 3, ˆr i = x i /r, ζ = ear, a being the vacuum expectation value of the Higgs field in the triplet representation, and T i being the generators of the SU) Lie algebra T i, T j ] = i ε ijk T k 3.) The exact analytical BPS monopole solution corresponds to the functions 5, 6] Kζ) = ζ sinh ζ Hζ) = ζ coth ζ 1 3.3) For the t Hooft-Polyakov monopole the functions Kζ) and Hζ) are obtained numerically, but they have qualitatively the same behaviour as 3.3), i.e. we have that K) = 1, and then it decays monotonically exponentially) to zero as r, and H) =, and then it grows monotonically with r and for r, such grow is linear in r. The function Hζ) will not be important in our calculations because the Higgs field does not appear in our integral equations for the case of static solutions and for Ω being purely spatial see.8) and.9)). The important simplification we obtain in our calculations is due to the fact that Kζ) is a monotonic function of ζ, and so it admits an inverse function. We will then trade the parameter ζ by the function K, and our calculations will not depend upon the explicit form of the function Kζ). 7

We have chosen to evaluate both sides of the integral equations.8) and.9) on a purely spatial 3-volume Ω which is a ball centered at the origin of the Cartesian coordinate system x i, i = 1,, 3, used in the ansatz 3.1). We then scan that volume Ω with closed surfaces which are spheres also centered at the origin of the Cartesian coordinate system, with radii varying from zero to the radius of Ω. However, since the reference point x R have to lie on the border Ω of Ω, and since the surfaces scanning it have to be based at x R, we shall attach to the ball Ω a infinitesimally thin cylinder lying on the negative x 1 -axis, and locate the reference point x R at x 1, x, x 3 ) =,, ). The cylinder has a radius ε, which will be taken to zero at the end of the calculations. The surfaces scanning Ω will have the form depicted in Figure, i.e. an infinitesimally thin cylinder on the negative x 1 -axis and a sphere centered at the origin of the Cartesian coordinate system. With the attachment of the thin cylinder we can keep the surfaces based at x R, and centered at the origin. In addition, x R being at infinity allows us to have the volume Ω with any radius. We shall label the surfaces scanning Ω with the parameter ζ, which is the same as the one appearing in the ansatz 3.1). Then ζ = corresponds to the surface made of the thin cylinder and a sphere of radius zero attached to it, and ζ = ζ corresponds to the border Ω, made of the thin cylinder attached to a sphere of radius ζ, the same as the radius of Ω, centered at the origin. The loops will be labelled by a parameter τ, they start and end at the reference point x R, and there will be three types of loops, as follows: 1. Loops of type I), scanning the thin cylinder, as depicted in Figure 1. For such loops the parameter τ varies from to π, with τ = corresponding to the infinitesimal loop around x R, and τ = π corresponding to a straight line from x R to the border of the sphere, then encircling the joint of the cylinder with the sphere, and coming back to x R through the same straight line. The three parts of such loops will be denoted I.1), the first straight line, I.), the circle and I.3) the second straight line. We parameterize the loops with σ, such that the points on the loops have the following coordinates: I.1) x 1 = τ + σ ζ + π x = x 3 = ε σ ) I.) x 1 = τ ζ + π x = ε sin σ x 3 = ε cos σ σ π) I.3) x 1 = τ + π σ ζ + π x = x 3 = ε π σ ) with fixed ζ and τ π. x R " x 1 Figure 1: Scanning of type I). The gap between the straight lines is only a visual resource.. Loops of type II), scanning the thin sphere, as depicted in Figure. For such loops the parameter τ varies from π to π. A loop of this type is made of a straight line from x R to 8

the border of the sphere, then making a circle on the surface of the sphere, starting and ending at the junction of the cylinder with the sphere, and lying on a plane perpendicular to the plane x 1 x 3, that makes an angle τ with the plane x 1 x. Finally, it returns to x R through the same straight line. Again, the three parts of such loops will be denoted as II.1) for the first straight line, II.) for the circle and II.3) labels the second straight line. We parameterize the loops with σ, such that the points on the loops have the following coordinates: II.1) x 1 = σ ζ x = x 3 = ε σ ) II.) x 1 = ζ cos τ1 cos σ) 1 ) x = ζ cos τ sin σ σ π) x 3 = ζ cos τ sin τ1 cos σ) II.3) x 1 = σ + π ζ x = x 3 = ε π σ ). with fixed ζ and where in II.) the parameter τ varies from π to π. x x R... x 1 x 3 Figure : Scanning of type II). 3. Loops of type III), scanning the thin cylinder backwards, as depicted in Figure 3. For such loops the parameter τ varies from π to, and they are made of two straight lines. The first one starting at x R and ending on some point on the side of the cylinder with coordinates x 1, x, x 3 ) = x 1,, ε). The second part of the loop is the same straight line reversed) going back to x R. We shall denote the first straight line III.1) and the second III.). We parameterize the loops with σ, such that the points on the loops have the following coordinates: III.1) x 1 = π τ ζ + σ x = x 3 = ε σ ) III.) x 1 = π τ ζ σ x = x 3 = ε σ ) with fixed ζ, and where π τ. An important simplification is made by observing that the Wilson line is constant along loops scanning the thin cylinder. Indeed, we observe that on the segments I.1), I.3), II.1), II.3), III.1) and III.), the coordinate x 1 is linear in σ, and x and x 3 are independent of it. Therefore, using 9

x R " x1 Figure 3: Scanning of type III). 3.1), we have that A i dx i dσ straight lines= ±A 1 = 1 e ε r 1 K) T as ε 3.4) with the upper signs valid for the segments I.1), II.1) and III.1), and the lower signs valid for I.3), II.3) and III.). On the segment I.), on the other hand, we get that dx i A i dσ I.) = ε cos σ A + sin σ A 3 ] 3.5) = ε 1 K) e r ε T 1 + τ ζ + π ) ] cos σ T 3 + sin σ T ) as ε The only non-vanishing contribution comes from the segment II.) which gives A i dx i dσ II.) = 1 e 1 K) cos τ cos τ sin τ1 cos σ)t 1 + sin τ sin σ T + sin τ1 cos σ) 1 ) T 3 ] = 1 e 1 K) cos τ ei τ T e i σ T 3 e i τ T T 3 e i τ T e i σ T 3 e i τ T 3.6) Therefore, integrating 1.7) one gets that the Wilson lines on the loops of type I and III are trivial, i.e. W I) = W III) = 1l. On the loops of type II one gets that W II) = W 3 W W 1, where W a are the Wilson lines obtained by integrating 1.7) on the segments II.a), a= 1,, 3. Due to 3.4) we have that W 1 = W 3 = 1l. Under a gauge transformation A i Āi = g A i g + i e ig g, with g = e i τ T e i σ T 3 e i τ T, one gets that W W = g f W g i = e i τ T e i π T 3 e i τ T W 3.7) where g i and g f are the values of g at the initial and final points of the loop II.), and so g i = 1l, and g f = e i τ T e i π T 3 e i τ T. Therefore, one gets that and the equation 1.7) for W becomes dx i A i dσ dx i II.) Āi dσ II.)= 1 e K cos τ T 3 sin τ T 1 ] 3.8) d W d σ + i K cos τ T 3 sin τ T 1 ] W =. 3.9) Since the connection term K ζ) cos τ T 3 sin τ T 1 ] is independent of t follows that the path ordering is unimportant and the integration on the loops II.) gives W = e i π K cos τ T 3 sin τ T 1 ] 3.1) 1

Using the fact that e i π T 3 = ±1l, depending if the representation used is of integer +) or half-integer ) spin, we get that g f = ± 1l, and so W = W II) = W = ± e i π Kζ) cos τ T 3 sin τ T 1 ] 3.11) where we have equated W to W II), because, as shown above W I) = W III) = 1l. Therefore, in 3.11) we have τ varying from π to π. The calculations concerning the Wilson line can be simplified defining γ as with Then 3.11) can be written as from which we get cos γ = K cos τ ; sin γ = sin τ F F 3.1) F ζ, τ) = K ζ) cos τ + sin τ. 3.13) W = ± e i γ T e i π F T 3 e i γ T, 3.14) W W = i e i γ T { π F T 3 + γ cos π F ) 1) T + sin π F ) T 1 ]} e i γ T = i { π F sin γ + γ sin π F ) cos γ] T 1 + γ cos π F ) 1] T + π F cos γ + γ sin π F ) sin γ] T 3 }. 3.15) We then have with W d W d τ = i cos τ N j K, τ) T j 3.16) N 1 K, τ) = 1 F π 1 K ) ] sin τ + K sin π F ) F N K, τ) = K F 1 cos π F )] 3.17) cos τ N 3 K, τ) = K sin τ F π 1 K ) ] cos sin π F ) τ +. cos τ F In addition W d W d τ, W d W d ζ ] d F = π i d τ d γ d ζ d γ ) d F d τ d ζ e i γ T 1 cos π F )) T 1 + sin π F ) T ] e i γ T = i π K cos τ M j K, τ) T j 3.18) with M 1 K, τ) = K cos τ F 1 cos π F )] M K, τ) = sin π F ) F M 3 K, τ) = sin τ F 1 cos π F )] 11 3.19)

where K stands for dk dζ, and where we have used the formulas d F d τ = sin τ cos τ F 1 K ), d F d ζ = K cos τ K, F d γ d τ = K F, d γ d ζ = sin τ cos τ F K. 3.) With these expressions we are ready to perform the calculations of section for the SU) monopoles. 4 Check of first order integral equations for SU) monopoles The integral equation for a purely spatial volume Ω, in first order in α, for the SU) monopoles t Hooft-Polyakov or BPS) is given by expression.8). However, since the Wilson line is unit for the loops of type I and III see section 3) we get that.8) is only non-trivial for loops of type II, where τ varies from π to π, and so.8) becomes V 1) = W dw ˆ ζ = dζ ζ=ζ π π W dw, W dw dζ ] = U 1) 4.1) where the l.h.s is an integration on a closed surface of radius ζ and the r.h.s is an integration in the volume contained inside that surface. Our goal is to evaluate both sides of this equation using the results obtained in the expressions 3.16) and 3.18). In order to perform the integration of the l.h.s term, a better choice of variables is the following: y = sin τ ; y 1 ; z = Kζ) cos τ ; z 1 4.) with ζ, Kζ) 1, and π τ π. Note that we are using the fact that Kζ) is monotonically decreasing function of ζ for both, the t Hooft-Polyakov and BPS monopole solutions. The explicit form of the function Kζ) is not important here. In these variables we get F = y + z = K + 1 K )y 4.3) and so using 3.16) and 3.17) we get that the l.h.s. of 4.1) becomes with K Kζ ), and V 1) = i dy N j K, y) T j 4.4) N 1 K, y) = π { } F y 1 K ) + K sinπf ) πf N K, y) = K {1 cos πf } 1 y F 4.5) N 3 K, y) = { } K y π F sin πf ) 1 +. 1 y F πf Note that N 3 is an odd function of y and so integrating we get V 1) K ) = i J 1 K )T 1 + i J K )T, 4.6) 1

with J 1 K ) = π dy J K ) = dy 1 K K + 1 K )y 1 K )y sin π ) K + 1 K )y + π K + 1 K )y )} K {1 1 y K + 1 K )y) cos π K + 1 K )y. Note that as ζ varies from to ζ, one has that K varies from 1 to K K ζ ) < 1. Therefore, the integration domain on the r.h.s. of 4.1) is a truncated semi-disc shown in Figure 4. The absolute value of the Jacobian of the variable transformation ζ, τ) y, z), given in 4.), is K cos τ = K cos τ, since K is strictly negative. In addition, it is more appropriate to perform a further change of variables to evaluate the integration on the r.h.s. of 4.1). We define the polar type coordinates s, θ) as y = s cos θ ; z = s sin θ ; S K, θ) s 1 ; θ π 4.7) with K SK, θ) p 1 cos θ1 K ) 4.8) Therefore one has that ˆ ζ ˆˆ dζ K cos τ = truncated semi-disc dzdy = π ˆ dθ ds s 4.9) SK,θ) 1.8.6 z = p 1 y z.4. z = K p 1 y 1.5.5 1 y Figure 4: The integration domain in the new polar coordinates. Each value of K fixes a new domain by shortening the area of the disk from below. 13

We then have that the M i s, defined in 3.19), become M 1 = z sin θ 1 cos π F )] = F s M = sin π F ) F = sin π s) s M 3 = y cos θ 1 cos π F )] = F s 1 cos π s)] 1 cos π s)] Using 3.18) we get that in these coordinates the r.h.s. of 4.1), denoted by U 1), reads 4.1) U 1) = iπ dθ SK,θ) ds {sin θ1 cosπs))t 1 + sinπs)t + cos θ1 cosπs))t 3 } 4.11) from where we can easily perform the integration in s, obtaining where I 1 K ) = I K ) = dθ U 1) K ) = ii 1 K )T 1 + ii K )T 4.1) { )} πk dθ sin θ π 1 cos θ1 K ) + sin πk 1 cos θ1 K ) { )} πk 1 cos 1 cos θ1 K ) 4.13) The integral along the T 3 -direction in 4.11) vanishes since the integrand is odd, under reflection around θ = π, in the interval θ π note that S θ) is even in that interval). Therefore, in order to check the validity of the integral equation at first order in α, given in 4.1), we have to verify the equalities of the coefficients of T i in 4.6) and in 4.1). We have performed the numerical integration of the quantities I i K ) and J i K ) for several values of K, covering the range 1 K, corresponding to ζ. Note that the actual value of K for a given value of ζ is different for the t Hooft-Polyakov and BPS monopoles. However, the fact that K ζ) is a monotonically decreasing function of ζ, for both solutions, allowed us to trade the coordinate ζ by K, and perform one check that is valid for the two monopole solutions. In section A.1 we give the results of the numerical integrations of the quantities I i K ) and J i K ), i = 1,. As one observes in those tables the values of I i K ) and J i K ) are remarkably identical, differing in the worst case around the eighth decimal place, due to the numerical approximation. This indicates that the t Hooft-Polyakov and BPS SU) monopoles are indeed solutions of the first order integral equation.1), or equivalently 4.1), appearing in the expansion in α of the integral non-abelian Gauss law in.5) and.7). 5 Check of second order integral equations for SU) monopoles The integral equation for a purely spatial volume Ω, in second order in α, for the SU) monopoles t Hooft-Polyakov or BPS) is given by expression.9). However, since the Wilson line is unit for 14

the loops of type I and III see section 3) we get that.9) is only non-trivial for loops of type II, where τ varies from π to π, and so.9) becomes V ) = + + π ˆ ζ ˆ ζ dζ dζ ˆ τ π π ˆ ζ dw W W dw ˆ τ π dζ π = ζ=ζ ˆ ζ dζ dw W, W dw, W dw dζ U 1) + G + G 3 = U ) π W dw, W dw dζ π ]] W dw, W dw dζ ] ] dw dw W, W dζ where we have denoted G and G 3 the terms appearing on the second and third lines respectively of 5.1). In addition, we have used the fact that the first term on r.h.s. of the first line of 5.1) is the same up to a minus sign) as U 1) given on the r.h.s. of 4.1). We start by evaluating the l.h.s. of 5.1), using 3.16), and 4.5) to get ˆ τ dw W W dw ˆ y 3 = dy dy N i K, y ) N j K, y) T i T j ζ=ζ π π i,j=1 ] 5.1) = 1 3 ] dy N i K, y) T i=1 i dy ˆ y dy 3 i j=1 N i K, y ) N j K, y) T i T j 5.) where in the first term on the r.h.s. of 5.) we have used the symmetry of the integrand in y and y to transform the integral on the triangle y 1 and y y, to the integral on the square y, y 1. We now use the fact that N i K, y) = ε i N i K, y), with ε i = 1 for i = 1, and ε 3 = see 4.5)). Then we can write dy ˆ y dy N i K, y ) N j K, y) = 1 + ε i ε j dy ˆ y dy dy N i K, y ) N j K, y) y dy N i K, y ) N j K, y) 5.3) Therefore for the case where ε i ε j = 1, one can write further that dy ˆ y dy N i K, y ) N j K, y) = 1 dy dy N i K, y ) N j K, y) ; ε i ε j = 1 For the case ε i ε j =, we do not use 5.3), but instead write T i T j = 1 {T i, T j } + 1 T i, T j ] = 1 {T i, T j } + i ε ijk T k 5.4) Note that we are dealing here with products, and not only commutators, of the SU) Lie algebra generators. We have therefore to work with a basis in the enveloping algebra of SU), which in the 15

case of quadratic terms we shall take to be the 9 quantities T i, and the anti-commutators {T i, T j }, i, j = 1,, 3. If one works with the spinor representation given by the Pauli matrices, i = 1,, 3, then one has σ j = i ε ijk σ k + δ ij 1l, and non-diagonal terms vanish, i.e. {, σ j } =, for i j. However, if one works with the triplet or higher representations one has {T i, T j } = even for i j. So, we have to consider the coefficients of all the 9 elements of the basis of the enveloping algebra to be independent. Therefore, using 5.4) one gets that ˆ τ V ) = dw W W dw = N 1 K )T1 + N K )T + N 1 K ) {T 1, T } ζ=ζ where π π + N + 13 K ) {T 1, T 3 } + N + 3 K ) {T, T 3 } i N 13 K ) T + i N 3 K ) T 1 ] N i K ) = 1 N 1 K ) = 1 N 13 ± K ) = 1 N 3 ± K ) = 1 dy N i K, y)) i = 1, ) dy N 1 K, y) dy dy ˆ y ˆ y dy N K, y ) ) ˆ y ) N 1 K, y )N 3 K, y) ± dy dy N 3 K, y )N 1 K, y) ˆ y dy N K, y )N 3 K, y) ± dy dy N 3 K, y )N K, y) with the N i s defined in 4.5), and where we have dropped the term proportional to T 3 because N 3 is an odd function of y, and so its integral on the interval y 1, vanishes. Using 3.16), 3.18) and 4.), the term on the second line of 5.1), denoted G, becomes ˆ K G = i π ε ijk T k dk 1 dy 1 y M j K, y) ˆ y dy N i K, y ) 5.5) 5.6) i 4 π R k K ) T k 5.7) Using 3.18) and 4.1) the term on the third line of 5.1), denoted G 3, becomes ) G 3 = 4 π dθ SK,θ) ds s dθ SK,θ ) ds s 3 M i s, θ) M j s, θ ) T i T j 5.8) i,j=1 with K K, and so ζ ζ. Note that in the θ, s )-integration, K has to be taken as a function of θ and s. From 4.) and 4.7) one gets that K =. Note that the 1 s cos θ θ, s )-integration is the same as the one performed in 4.11), with K replaced by K. Therefore, similar to what happened, there will be no terms in the direction of T j for j = 3, since M 3 s, θ ) is odd under reflection of θ around θ = π see 4.1)). Since K and S K, θ) are even under the reflection of θ around θ = π, there will be no terms in the direction of T i for i = 3, since M 3 s, θ) is odd under that reflection. Using 5.4) one gets that G 3 = 4π S 1 K )T1 + S K )T ] + S 1 {T 1, T } + is 3 K )T 3 16 s sin θ

with S a K ) = π S 1 K ) = π S 3 K ) = π dθ dθ SK,θ) SK,θ) ds s ds s dθ dθ SK,θ ) SK,θ ) ds s M a s, θ) M a s, θ ) ; a = 1, ds s M 1 s, θ) M s, θ ) + M s, θ) M 1 s, θ )] 5.9) dθ SK,θ) ds s dθ SK,θ ) ds s M 1 s, θ) M s, θ ) M s, θ) M 1 s, θ )] The s -integration can be performed analytically and so, using 4.1) and the fact that K = s sin θ 1 s, we get cos θ and SK,θ ) SK,θ ) ds s M 1 s, θ ) = sin θ 1 s sin θ s sin θ + 1 s ) sin θ )] 1 π sin π s sin θ s sin θ + 1 s ) sin θ )] ds s M s, θ ) = 1 π s sin θ + cos π s sin θ + 1 s ) sin θ 5.1) 5.11) Note that the above integrals are symmetric under the reflection of θ and θ around π. The quantities M 1 s, θ), M s, θ), S K, θ) and K s, θ) are also symmetric under the reflection of θ around π. Therefore, the integration in θ and θ can be performed in the interval from zero to π, by multiplying 17

the result by two. So, we then get that S 1 K ) = P 1 K ) = dθ ds SK,θ) dθ 1 cos πs)) sin θ sin θ )] πs sin θ π s sin θ + 1 s ) sin θ + sin πs sin θ s sin θ + 1 s ) sin θ P 1 K ) = dθ + cos π SK,θ) ds dθ 1 cos πs)) sin θ sin θ )] πs sin θ s sin θ + 1 s ) sin θ dθ ds dθ sin πs) S θ) πs sin θ s sin θ + 1 s ) sin θ + sin )] πs sin θ s sin θ + 1 s ) sin θ 5.1) 5.13) 5.14) S K ) = where we have introduced ˆ π dθ ds dθ sin πs) + cos SK,θ) )] πs sin θ s sin θ + 1 s ) sin θ S 1 = 1 P 1K ) + P 1 K )) S 3 = 1 P 1K ) P 1 K )), 5.15) Summarizing, we have obtained both sides of the integral equation in second order in α for a given K, given in 5.1). From 5.5) we have that V ) = in3 T 1 + in13 T N 1 T1 N T N 1 {T 1, T } N 13 + {T 1, T 3 } N 3 + {T, T 3 } 5.16) and, from 5.1), 4.1), 5.7) and 5.9) we have that U ) = ii 1 +4π R 1 )T 1 + ii +4π R )T i4π R 3 +4πS 3 )T 3 4πS 1 T1 4πS T 4πS 1 {T 1, T }. 5.17) We have to check the equality between the coefficients of each element of the basis of the SU) enveloping algebra on the expansion of V ) and U ). Those coefficients involve integrals which are calculated numerically for a set of values of K. The results are presented in the tables of section A. in the Appendix. The consistency is remarkable and with that check we can state clearly that the the t Hooft-Polyakov monopole and its BPS version satisfy the integral Yang-Mills equations up to second order in α. 18

6 Conclusions The integral Yang-Mills equations appeared from an attempt to understand integrability in higher dimensional space-times 1, ]. Through a loop space formulation 11, 1] one can construct a suitable generalization of the non-abelian Stokes theorem for two-form fields that can be used naturally to define conservation laws, thus mimicking the so-called zero curvature representation of integrable field theories in 1 + 1)-dimensions. That has led us to consider the applications of such non-abelian Stokes theorem to construct the integral equations for non-abelian gauge theories, generalizing the well known abelian version of such integral equations used to describe the laws of electrodynamics. That was indeed possible, as we have shown in 1, ], and the usual differential Yang-Mills equations are obtained from these integral equations when the appropriate limit is taken. The present paper shows that there is more to be explored. The integral Yang-Mills equations allow the introduction of two c-numbers as parameters which arise naturally in the construction of the equations, and as being non-linear, produce a quite non-trivial dependence on those parameters of the surface and volume ordered integrals appearing on both sides of the equation. We have tested the assumption that the integral Yang-Mills equations are in fact a collection of an infinite number of equations, each one corresponding to the coefficients of the above mentioned expansion in powers of those parameters. This was done by considering the fact that, by construction, a solution of the differential Yang-Mills equation is also a solution of the integral Yang-Mills equation. Thus, using the t Hooft-Polyakov and BPS monopoles as such configurations, we tested the validity of the equations arising at first and second order in the parameter expansion of the integral Yang- Mills equation. Despite the quite different structures of the terms resulting from the surface and volume ordered integrals, we have checked their equalities with a high numerical precision of at least one part in 1 7. In addition, much of the check has been done analytically, and we have obtained an exact expression for the Wilson line operator, on each loop scanning the surfaces and volumes, for the SU) t Hooft-Polyakov monopole solution and its BPS version see 3.14)). That result can certainly be useful in many other applications. The fact that those configurations are solutions of both of the highly non-trivial equations at each order of the expansion, indicates that the parameters could indeed be arbitrary. The arbitrariness of the parameters leads to a variety of important consequences which can now be considered, such as their role in the conserved charges that arise dynamically from the integral equations and the significance of having an infinite number of integral equations. Aknowkedgments: The authors are grateful to partial financial support by FAPES under contract number 447/15. LAF is partially supported by CNPq-Brazil. 19

A Numerical Results In this section we show the results of the numerical integrations related to the terms on the l.h.s and r.h.s of the expansion of the integral equation performed at first and second order in α. The coefficients of the generators of the algebra eventually, up to a common factor of i ) are compared for different values of K. For each integral estimative there is an associated upper bound on the error, which we represent by using the following notation: 1.37 ±.8 1.37 ± 8). A.1 Equation V 1) = U 1) Coefficients of T 1 K I 1 K ) J 1 K ).1 1.561418 ± ) 1.5614186.1 1.771874199 ± 6) 1.771874199. 1.771916 ± ) 1.7719168.3 8.68787584 ± 8) 8.6878758454.4 6.38863598 ± 4) 6.3886359858.5 4.1697936 4.1697936.6.385155685 ± 6).385155685.7 1.3849 ± 1) 1.384978.8.315114413 ± 9).3151144136.9.3915944383.3915944383.99.3798 ± 3).3798666 ± 1) Coefficients of T K I K ) J K ).1.191716844 ± 3).19171684 ± 1).1-1.85885115± 9) -1.8588511± ). -3.37694768± 1) -3.37694768± 4).3-4.97856758-4.978567 ± 6).4-4.54181669 ± 8) -4.5418166± ).5-4.499388345-4.499388± ).6-3.156196± ) -3.156196± ).7-1.97418488± 6) -1.9741848± 1).8 -.9381858851± 6) -.938185885± 8).9 -.391335± 3) -.391335± 9).99 -.33663984± 4) -.33663 ± 7) Table 1: Numerical verification of the validity of equation 4.1): the coefficients of T 1 and T in 4.6) and in 4.1) agree up to the eighth order.

A. Equation V ) = U ) The tables below show the values of the coefficients of the algebra elements of 5.16) and 5.17). The fact that they agree implies on the validity of the equation obtained after expanding the Yang- Mills integral equation to second order in α and therefore, on the validity of the integral equation itself for any value of α, at least up to that order. Coefficients of T 1 K I 1 K ) + 4π R 1 K ) N 3 K ).1.16581± ).1658175± 3).1 1.139151158± 6) 1.13915115± 4). 3.47985977± ) 3.47985977± 1).3 6.616113± 8) 6.616114± 6).4 7.311678563± 4) 7.31167856± 1).5 6.7761777766± 5) 6.77617777± 9).6 4.8561878177± 6) 4.85618781± 4).7.56959344± 1).56959344117.8.87196± 9).8719± 9).9.1155473857.11554738± 3).99.113953± 3).11395339± 1) Coefficients of T K I K ) + 4π R K )) N 13 K ).1 -.46994457± 3) -.46994457± 1).1-4.3831934± 9) -4.3831934± 1). -7.5679398± 1) -7.5679398± 1).3-7.85686-7.856± 1).4-5.67178717± 8) -5.67178717± 6).5 -.51857335496± 9) -.518573354± ).6 -.9619543± ) -.96195434± 7).7.7148861± 6).7148861± 7).8.66759135± 6).6675913± 4).9.17786± 3).17786± 1).99.33498± 4).3349± 8) Table : Comparison between the coefficients of T 1 and T of equations 5.16) and 5.17). 1

Coefficients of T1 K 4πS 1 K ) N 1 K ).1 78.8943986± ) 78.894398654.1 7.997987± 8) 7.997988. 57.56846± ) 57.5684639.3 37.739593357 37.739593357.4.4733451± 3).473345139.5 8.69611198± ) 8.696111984.6.719937± 3).719937533.7.5387646± 1).5387646116.8.496479511± ).49647951165.9.7667316± 1).76673164.99 7.131853995 1 7.131853989± ) 1 Coefficients of T K 4πS K ) N K ).1.183776741± ).1837767± 1).1 1.7661178± 1) 1.7661178± ). 5.718878± 4) 5.718878± 7).3 9.357861936± 1) 9.357861± 1).4 1.143865± 4) 1.143865± 5).5 8.17899778± 4) 8.1789977± 4).6 4.84365889 4.843658± 3).7 1.9451734± ) 1.9451734± 1).8.44963774± 5).449637± 1).9.85636848378.856368± 1).99.798598954e-6.798± 8)e-6 Table 3: Comparison between the coefficients of T 1 and T of equations 5.16) and 5.17).

Coefficients of {T 1, T } K 4πS 1 K ) N 1 K ).1-1.4116884± ) -1.41168813.1-11.1488± 5) -11.148884. -18.6787384± 7) -18.67873847.3-18.669678-18.6696777.4-14.3877884± 1) -14.38778844.5-8.477816688± 6) -8.477816684.6-3.6157418± 4) -3.615741851.7-1.3689433± 4) -1.368943394.8 -.1478168413± 1) -.14781684139.9 -.46798144471± 3) -.467981444135.99-4.437451947 1 8-4.437477± 7) 1 8 Table 4: Comparison between the coefficients of {T 1, T } of equations 5.16) and 5.17). 3

Coefficients of T 3, {T, T 3 } and {T 1, T 3 } K 4π R 3 K ) + 4πS 3 K ) N 3 + K ) N 13 + K ).1 ± 1 8 ± 7 1 1 ± 9 1.1 ± 7 1 8 ± 1 9 ± 6 1 9. ± 1 1 7 ± 7 1 9 ± 1 1 8.3 ± 1 1 ± 1 1 8 ± 4 1 9.4 ± 7 1 8 ± 1 1 8 ± 6 1 9.5 ± 7 1 ± 9 1 9 ± 7 1 9.6 ± 5 1 9 ± 4 1 9 ± 7 1 9.7-4.445± ) 1 1 ± 6 1 1 ± 5 1 9.8-1± 7 1 1 ± 9 1 ± 3 1 9.9-3.7± 7) 1 3 ± 3 1 ± 5 1 9.99 9.599± 4) 1 6-1.6± ) 1 1 ± 8 1 9 Table 5: The coefficients above are the ones that should vanish in the equation obtained at second order in α; indeed, within a numerical precision, they are zero. 4

References 1] L. A. Ferreira and G. Luchini, Integral form of Yang-Mills equations and its gauge invariant conserved charges, Phys. Rev. D 86, 8539 1); arxiv:15.88 hep-th]]. ] L. A. Ferreira and G. Luchini, Gauge and Integrable Theories in Loop Spaces, Nucl. Phys. B 858, 336 1); arxiv:119.66 hep-th]]. 3] G. t Hooft, Magnetic Monopoles in Unified Gauge Theories, Nucl. Phys. B 79, 76 1974). doi:1.116/55-31374)9486-6 4] A. M. Polyakov, Particle Spectrum in the Quantum Field Theory, JETP Lett., 194 1974); Pisma Zh. Eksp. Teor. Fiz., 43 1974)]. 5] M. K. Prasad and C. M. Sommerfield, An Exact Classical Solution for the t Hooft Monopole and the Julia-Zee Dyon, Phys. Rev. Lett. 35, 76-76 1975). 6] E. B. Bogomolny, Stability of Classical Solutions, Sov. J. Nucl. Phys. 4, 449-454 1976) Yad. Fiz. 4, 861 1976)]. 7] N. S. Manton and P. Sutcliffe, Topological solitons, Cambridge Monographs on Mathematical Physics, 4, Cambridge University Press 8] E. J. Weinberg, Classical Solutions in Quantum Field Theory : Solitons and Instantons in High Energy Physics, Cambridge Monographs on Mathematical Physics, 1, Cambridge University Press 9] C. P. Constantinidis, L. A. Ferreira and G. Luchini, A remark on the asymptotic form of BPS multi-dyon solutions and their conserved charges, Journal of High Energy Physics, JHEP 151, 137 15); doi:1.17/jhep115)137; arxiv:158.349 hep-th]]. 1] C. N. Yang and R. L. Mills, Conservation of Isotopic Spin and Isotopic Gauge Invariance, Phys. Rev. 96 1954) 191. doi:1.113/physrev.96.191 11] O. Alvarez, L. A. Ferreira and J. Sanchez Guillen, A New approach to integrable theories in any dimension, Nucl. Phys. B 59, 689 1998); doi:1.116/s55-31398)4-3; hepth/971147]. 1] O. Alvarez, L. A. Ferreira and J. Sanchez-Guillen, Integrable theories and loop spaces: Fundamentals, applications and new developments, Int. J. Mod. Phys. A 4, 185 9); doi:1.114/s17751x943419; arxiv:91.1654 hep-th]]. 5