Effective theory of the D = 3 center vortex ensemble

By means of lattice calculations, center vortices have been established as the infrared dominant gauge field configurations of Yang–Mills theory. In this work, we investigate an ensemble of center vortices in D = 3 Euclidean space-time dimension where they form closed flux loops. To account for the properties of center vortices detected on the lattice, they are equipped with tension, stiffness and a repulsive contact interaction. The ensemble of oriented center vortices is then mapped onto an effective theory of a complex scalar field with a U(1) symmetry. For a positive tension, small vortex loops are favoured and the Wilson loop displays a perimeter law while for a negative tension, large loops dominate the ensemble. In this case the U(1) symmetry of the effective scalar field theory is spontaneously broken and the Wilson loop shows an area law. To account for the large quantum fluctuations of the corresponding Goldstone modes, we use a lattice representation, which results in an XY model with frustration, for which we also study the Villain approximation.


Introduction
The essential features of the QCD vacuum are confinement and the spontaneous breaking of chiral symmetry. A thorough understanding of these two phenomena, and of the infrared sector of QCD in general, is still lacking. However, substantial progress has been achieved during the last twenty years in identifying the relevant infrared degrees of freedom of QCD (or at least of Yang-Mills theory). From these studies consistent confinement pictures have emerged: The center vortex picture [1][2][3][4][5][6], the dual Meissner effect [7][8][9][10] and the Gribov-Zwanziger picture [11,12]. These different scenarios do not contradict each other but turn out to be just different manifestations of the same phenomena in different gauges. Center vortices detected on the lattice by the method a e-mail: oxman@if.uff.br of center projection [13] show the proper scaling behaviour towards the continuum limit [14,15] only in the so-called maximal center gauge. Analogously, magnetic monopoles are detected after Abelian projection and show proper scaling only in the maximal Abelian gauge [16][17][18][19]. (The so-called "indirect maximal center gauge" is done on top of the maximal Abelian gauge). Finally, the Gribov-Zwanziger picture [11,12] has been established in Coulomb gauge [20][21][22].
Magnetic monopoles are attached to center vortices [23] and change the direction of the flux of center vortices [24]. 1 Therefore, condensation of center vortices in the confined phase implies also the condensation of magnetic monopoles and thus the dual Meissner effect. Center vortices as well as magnetic monopoles live on the Gribov horizon in both Coulomb and Landau gauge [25]. Configurations on the Gribov horizon give rise to an infrared diverging ghost form factor, a necessary condition for the Gribov-Zwanziger picture to be realized. When center vortices are eliminated from the ensemble of gauge field configurations contributing to the Yang-Mills functional integral, the ghost form factor becomes infrared finite and confinement is lost [26,27]. Furthermore, the Gribov-Zwanziger mechanism in Coulomb gauge presupposes the dual Meissner effect [28].
Center vortices are localized field configurations for which the Wilson loop operator becomes a center element of the gauge group if the latter is non-trivially linked to the center vortex, see Eq. (1) below. Lattice calculations performed in the maximal center gauge have provided mounting evidence that the center vortices are the dominant infrared degrees of freedom of Yang-Mills theory: When center vortices are removed from the gauge field ensemble of the lattice functional integral [29] the string tension is lost, chiral symmetry is restored [30] and topological charge is lost [31]. The emergence of the string tension can be easily understood in a random center vortex model [32]. Furthermore, center vortices give also a simple explanation of the deconfinement phase transition [32].
In an idealized picture, which is realized, in particular, after center projection on the lattice, center vortices are pointlike objects in D = 2, closed strings in D = 3 and closed surfaces in D = 4. They are closed by the Bianchi identity and live on the dual lattice. The gross features of center projected Yang-Mills theory, like the emergence of the string tension or the deconfinement phase transition, can be reproduced in a center vortex model with an action given by the vortex area plus a penalty for the curvature of the vortices [33][34][35]. The latter accounts for the stiffness of the vortices. In addition, Casimir scaling can be reproduced by the vortex thickness [36]. In D = 4 the model has to be numerically simulated on the lattice since a (continuum) string theory cannot be treated analytically. Since D = 3 Yang-Mills theory has the same infrared properties as D = 4 it is useful to investigate the center vortex model in D = 3, where vortices are closed loops.
In this paper we study the ensemble of closed center vortices in D = 3 exploiting the fact that the partition function of a gas of one dimensional objects can be represented by a complex scalar quantum field theory [37][38][39]. Within this theory, we calculate the Wilson loop. To keep the soft Goldstone modes of the complex scalar field, we resort to a lattice representation which leads to the 3D XY model with frustration for which we also study the Villain approximation. Since we are interested here in the Wilson loop, we consider only oriented center vortices as the orientation (of the flux) of center vortices is irrelevant for the confining properties. Non-oriented center vortices arise in the presence of magnetic monopoles [24]. Such vortices were considered in D = 3 in Ref. [40]. Let us also mention that an ensemble of closed center vortices in 3D, generated by Monte Carlo methods applied to lines with stiffness that can grow, shrink and reconnect, was recently considered in Ref. [41]. Furthermore, Monte Carlo simulations to explore the statistical properties of the 3D XY model using a disorder parameter that creates flux vorticity were carried out in Ref. [42].
The organization of the paper is as follows: In the next section we consider the partition function of an ensemble of ideal center vortices in D = 3 Euclidean space time dimensions in the presence of a Wilson loop. The center vortex loops are endowed with tension, stiffness and a binary repulsive interaction. The latter is linearized by means of a real scalar field. Then the partition function of the vortex ensemble is reduced to the quantum transition amplitude of a single center vortex in an external scalar field. In Sect. 3 this amplitude is expanded in leading order in spherical harmonics and transformed into an effective theory of a complex scalar field. This field develops a non-trivial vacuum expectation value. To study the quantum fluctuations of this field, in Sect. 4 the theory is translated to a lattice where it results in an XY model, which is also investigated in the Villain approximation. Finally the Wilson loop is calculated at low and high temperatures in Sect. 5. A short summary and our conclusions are given in Sect. 6.

Ensemble of thin center vortices in 3D
Lattice calculations in pure SU (N ) Yang-Mills (YM) theory have established center vortices as the essential degrees of freedom underlying confinement. While center vortex removal leads to a perimeter law, the projection to the center vortex ensemble gives an area law, with the N -ality properties observed in the complete Monte Carlo simulations. By the Bianchi identity, center vortices form closed manifolds, i.e. closed loops in D = 3, of color electric or magnetic flux.
The effect of a thin center vortex on a Wilson loop W C is topological. In a three dimensional Euclidean spacetime, when a closed center vortex worldline l links the Wilson loop C the latter gains a factor where z is an element of the center Z (N ) of the gauge group SU (N ) and L(l, C) ∈ Z is the linking number, which is a topological invariant that counts the number of times the loop l winds around the loop C, with a sign that depends on whether l has positive or negative orientation with respect to C. An explicit integral formula was given by Gauss, where x μ (s) andx μ (τ ) parametrize l and C, respectively.
The center elements of SU (N ) are given by the N th roots of unity where 1 denotes the N -dimensional unit matrix. Like all group elements, the center elements can be represented as exponentials of Lie algebra-valued vectors where ν(k) are the co-weights, which live in the Cartan subalgebra (with generators H a in the fundamental representation) and define the corners of the fundamental domain (Weyl alcove) of the SU (N ) algebra. From Eq. (1) it is clear that a center vortex is connected with a non-trivial center element z(k = 0). Since z(k) = z(1) k , there exist vortex branching for SU (N ≥ 3), which we will, however, not consider in the present paper. Therefore, in the following we will consider only vortices connected with z(1). This is sufficient for SU(2) gauge theory.
Our objective is to compute the large distance behavior of the ensemble average z(C) , summing over all possible numbers and shapes of closed vortex worldlines, after equipping them with appropriate physical properties.
Initially, we can rewrite (2) 2πν in terms of a gauge field Then the contribution of the center vortex loop l to the Wilson loop C in Eq. (1) becomes Using Stoke's theorem in Eq. (6) where the linking number can be expressed as the intersection number between l and a surface S(C) bordered by C. Note that A l μ (x) and B l μ (x) represent the gauge potential and the dual field strength of the center vortex localized on l. Indeed from Eq. (7) follows where we adopted s as the arc length parameter of the center vortex and defined Defining a vector field J C μ (x) localized on S(C) wherex(σ 1 , σ 2 ) is a parametrization of S(C), the contribution of a center vortex l belonging to the center element z = e i2πν = e i2π k/N can be written as Since the linking number L(l, C) is symmetric with respect to the interchange of the loops, L(l, C) = L(C, l), performing this change in Eq. (14) and comparing with Eq. (8), it follows that the vector fields A l μ (x) and J l μ (x) have to be equivalent. Indeed, A l μ (x) (7) can be gauge transformed to J l μ (13) [24,43]. (The former satisfies ∂ μ A l μ (x) = 0). Hence the vector field J C μ (x) (13) represents the gauge potential of a center vortex whose worldline is given by C (instead of l).
For a set of n center vortex wordlines l k , k = 1, . . . , n, parametrized by x (k) (s k ), s k ∈ [0, L k ], the contribution to the Wilson loop is, In order to identify observables in the center vortex ensemble with those in the effective field description to be derived below, it will be convenient to proceed with a general J μ . Furthermore, since the coweights ν = ν a (k)H a occur in the following equations only in the form e i2πν = e i2π k/N we can replace below ν by k/N , so that the field J μ (x) will no longer be algebra-valued. In addition, we will confine ourselves to a single vortex type k = 1, which is sufficient for the gauge group SU(2). The center vortex ensemble obtained in D = 4 Euclidean Yang-Mills theory after center projection can be modelled by vortices that are distributed according to an action which contains the vortex length and its curvature [33][34][35]. We therefore describe the intrinsic properties of the center vortices (i.e. n center vortex loops Here μ is the tension of the center vortices, i.e. the action cost per length, and 1/κ measures their stiffness. The larger κ the more flexible are the vortex loops since a smaller penalty is given for the curvatureu μuμ . As is well-known, a finite stiffness 1/κ is crucial to get a well-defined continuum limit when the worldlines are thought of as polymers, and are discretized in terms of monomers. Regarding the tension parameter, positive and negative μ favors small and very long (percolating) vortices in the ensemble, respectively. In Ref. [44] it was also shown that the center vortices in YMT do interact and their interaction scales properly in the continuum limit. Therefore we give the center vortex loops a binary interaction of the form resulting in the partition function Here the measure [Dl] n integrates over all the possible realizations of n center vortices and will be specified later, see Eq. (26) below. The average of the Wilson loop for the center vortex ensemble is obtained from with J C μ given by Eq. (13). We shall consider repulsive contact interactions which account for excluded volume effects, Then the interaction term of the action can be expressed by the scalar vortex density yielding It is convenient to linearize this term by means of a scalar field φ(x) Then, we arrive at where and the measure [Dl] n is given by, Here the variables x, u were collectively denoted by v, and [Dv(s)] L v,v integrates over center vortex worldlines of fixed length L starting and ending at the same position x with the same initial and final tangent vector u, which corresponds to smooth closed loops. In Eq. (24) the summation over the number of vortices can be carried out explicitly yielding where is the end-to-end probability for a worldline of length L to start at x 0 , with tangent u 0 , and end at x with u. duq(v, v, L) represents the partition function of a single closed vortex line of fixed length L. Although we will consider a general N , for N ≥ 3 center vortex ensembles also involve fusion rules, conserving the topological charge. This type of branching was not included in the discussion above. That is, in the Yang-Mills context, our vortex model will be particularly appropriate to describe vortex ensembles in SU (2) Yang-Mills theory. Further comments about this point will be given in Sect. 6.

Ensemble average
To perform the ensemble average, we shall closely follow the calculations given in Refs. [40,45], based on polymer techniques developed in Refs. [46,47], which are briefly reviewed here. Polymers are macromolecules formed by large linear arrays of small molecules (monomers). To understand systems of many interacting polymers, it is essential to initially characterize the behavior of a single polymer, which in turn depends on its microscopic structure. Among the commonly used descriptions, the continuous Gaussian model includes local stretching but no stiffness. On the other hand, the wormlike model is characterized by a fixed length polymer with an energy cost for local bending. This is indeed the case represented in Eq. (28). In the interacting wormlike model, the end-to-end probability satisfies a Schrödinger-like diffusion equation in real (x) and tangent (u) space. This has permitted to understand semiflexible polymers and the effect of solvent inhomogeneities and anisotropies in a controlled manner [46]. In the flexible limit, the formal solution to this equation is also the starting point to derive an effective field theory, when summing over all possible lengths L and numbers of monopole worldlines. In the functional integral (28), all the paths have fixed length L so they can be discretized in terms of M small segments ("monomers") of length L = L/M between the points x j and x j+1 , we have, Now, separating the integral over d 3 x M−1 d 2 u M−1 in Eq. (29) and where is the end-to-end probability for a line with initial condition x 0 , u 0 and length L − L, formed by M − 1 monomers, to end at x , u . Integrating over x and using the notation of the continuum, with an infinitesmial L, we find For a finite κ, the terms first order in L lead to, with the initial condition, The operatorL 2 u is the Laplacian on the sphere u μ ∈ S 2 . It arises from expanding q(. . . , u , . . .) in the integrand of Eq. (33) in powers of u − u, and computing the moments of the distribution ψ(u − u ). 2 As u μ carries angular momentum l = 1, in an expansion of q in terms of spherical harmonics, Eq. (34) couples the different l-sectors Using the completeness relation, the initial condition now reads, For small stiffness and large L, the correlations between the initial tangent direction u 0 and the final one u become small, thus favouring small angular momenta. In the flexible limit, which corresponds to a finite but small stiffness, the angular momenta l ≥ 2 can be disregarded in the expansion (34). In this case, the solution to Eq. (34) can be approximated by [46], Indeed, using Eqs. (41)- (43), it can be easily verified that where the second member involves an operator carrying angular momentum l = 2 (a traceless symmetric tensor), which was disregarded in the ansatz (41). Then, the equations close when restricted to the l = 0, 1 sectors. For a discussion involving the coupled equations for the whole tower of angular momenta, see Refs. [40,45]. Summarizing, keeping the dominant term Q 0 , and using Eq. (42), with the l = 0 initial condition in Eq. (40), the end-to-end probability turns out to be, Inserting this expression into Eq. (27), and using we find the following representation of the partition function for center vortices with small but non-zero stiffness 1/κ To exhibit the physical meaning of the complex field V in terms of the initial center vortex ensemble, we express Z [J μ ] in Eq. (18) using the total dual field strength B μ of the thin center vortices, [c.f. Eq. (11)] and include an external source J (x) to the scalar vortex density ρ(x), After linearizing the interaction term, this amounts to the substitution φ(x) → φ(x) + J (x) in Eq. (25). Then, following the same steps that led to Eq.
Taking functional derivatives in Eqs. (51) and (52), with respect to J (x 1 ), J (x 2 ), . . . and J μ 1 (x 1 ), J μ 2 (x 2 ), . . . at noncoinciding points, and setting the external sources to zero (including the sources J μ (x) induced by the Wilson loop), we find the following correspondence between correlation functions (in the absence of the loop C), where Since reversing the orientation of the vortex flux changes the sign of B μ (x), correlation functions with an odd number of B μ (x)'s vanish.

XY and Villain models with J μ (x)
In the previous sections, we have obtained the effective field representation for the average of center elements z(C) given by Eqs. (19) and (49). On the other hand, in the initial center vortex ensemble this average is represented by, [cf. Eqs. (18) and (50)]. Again, we have used the fact that vortex configurations come in pairs of opposite orientation, B μ and −B μ , to get an explicitly real expression, in accordance with the real integrand in Eq. (49). This field representation is valid for flexible vortices (small but nonzero values of 1/κ), a condition that has permitted us to keep only the smallest angular momenta in the tangent u-space and to obtain a quadratic kinetic term. If the initial Yang-Mills theory were coupled to a set of Higgs fields, such that the center vortices emerged as classical saddle points of the action, the parameter μ would be positive. This action cost would lead to an ensemble of small loops, and a perimeter law for large Wilson loops, see Sect. 5. Let us analyze this situation from the point of view of the effective field theory. When μ > 0, the functional integral over the complex field V in Eq. (49) can be computed perturbatively taking as reference the quadratic Lagrangian Then the partition function is dominated by the functional determinant of the inverse propagator of the massive field V , with squared mass 3κμ > 0, Recalling that the complex field V is minimally coupled to the "gauge field" J μ [cf. Eq. (34)], and that the effective action is gauge invariant, the partition function Z 0 [J μ ] can only depend on the "field strength" μνρ ∂ ν J ρ . That is, the average z(C) , where J μ = J C μ , depends on μνρ ∂ ν J C ρ . Now, from the definition of J C μ (x) (13), it follows that which is localized on the Wilson loop C. This, together with the nontrivial mass scale 3κμ > 0 implies a perimeter law in the Higgs phase. On the other hand, from lattice simulations, we know that center vortices percolate in the pure Yang-Mills vacuum and that, from the ensemble point of view, this leads to an area law for Wilson loops. Percolated vortices are necessary large and require μ < 0 in the vortex action. Let us now investigate how in this case the area law emerges in the effective field description. For μ < 0 we have for the potential term in the effective field theory (49) with v 2 = −μζ > 0. This potential breaks the underlying global U(1) symmetry spontaneously. Dropping the irrelevant constant in Eq. (60), Eq. (49) becomes, For sufficiently weak vortex interactions the potential in Eq. (61) tolerates only small fluctuations in the field h(x), while γ (x) is a Goldstone field, whose fluctuations are not restricted by the potential. Furthermore, γ (x) is a compact field defined modulo 2π . To have a well-defined description of the soft degrees of freedom v e iγ , and to keep their compact character, we switch to the lattice version of Eq. (61), where we have ignored the small amplitude fluctuation field h(x), putting ρ(x) to its vacuum value v. Note that the relevant dimensionful parameter here is not v 2 (scalar density) but √ η = v 2 /3κ, which also controls the vector vortex (current) density [cf. Eq. (54)]. Then, this parameter is expected to be related to the number of vortices intersecting a given surface per unit area, that for dimensional reasons should scale as ∼ η 2 .
For a cubic lattice with M sites x, spacing a, and oriented linksμ, the discretized covariant derivative is, When J C μ is smooth, α μ (x) = a J C μ (x). From the explicit form of J C μ (13) follows that α μ (x) = 2π N if the surface S(C) is crossed by the link (in the direction of the normal to S(C)), and zero otherwise. Therefore, we are led to the 3d XY model with frustration α μ (x), and partition function where The normalized average of the Wilson loop (19) is given by In fact, to make contact with the continuum, we are interested in the critical region where the correlation lengths become large with respect to the lattice spacing. In the literature, studies about the frustrated 3d XY model can be found for specific realizations of α μ (x). The fully frustrated case, with homogeneous frustration vector, has been extensively analyzed (see Ref. [49] and references therein). This vector has x, y and z components given by the circulation of α μ (x) along plaquettes on the yz, zx and x y-planes, respectively. Different frustration vectors have been studied, each one displaying its own critical properties. The case of a random phase shift has been discussed in Ref. [50]. To the best of our knowledge, there are no studies for the 3d XY model with phase shifts localized on a geometric region, as needed to compute the Wilson loop. However, in this case, the phase shift vanishes along the whole lattice but on those links that cross S(C). Moreover, the frustration vector is only nonzero on plaquettes that contain just one link with nontrivial shift, which are placed at the border of S(C). Then, to analyze z(C) latt in Eq. (66), we shall assume that the thermodynamic properties of the system 3 are those of the problem without frustration, that has a critical point at β c ≈ 0.454 (see Ref. [51] and references therein). In the critical regime, different models within the same universality class can be used. Outside this region, at very small (large) β, which means large (small) quantum fluctuations, the details are in general model dependent. Let as describe what happens when we go from very small β to β c . In this region, the XY model is in excellent agreement with the Villain model, which is in the same universality class. Let us summarize the main steps underlying this approximation following Ref. [51], where the XY and Villain approximations (without frustration) were extensively reviewed and studied for a superfluid. An expansion in powers of β leads to integrals of products of cosine functions. To organize the calculation, it is more convenient to use the Fourier decomposition, where I b (β) is the modified Bessel function of integer order b. Then, introducing integer valued variables b μ (x), integrating by parts on the lattice, and over the γ -variables yields The summation over {b μ (x)} runs over non-backtracking oriented closed loops of unit strength. On a given link (x,μ), if just one loop passes with the same (opposite) orientation aŝ μ, then b μ (x) takes the value +1 (−1). If n ≥ 2 loops pass on this link, all with the same orientation, then the variable takes the value b μ (x) = ±n depending on how the loops are oriented with respect toμ. These loops are analog to the fluxes of B μ through the plaquettes. The Villain approximation to Z XY (α μ ) in Eq. (66) is given by the replacement, with = ln(I 1 (β)/I 0 (β)).
This leads to Next, the Gaussian weights can be linearized with continuous real fields C μ (x). Then, the sum over n μ (x) can be carried out explicitly using the Poisson formula. This replaces the C μ (x) integrals by a sum over integers b μ (x), Finally, using Eq. (72), one finds which is to be compared with Eq. (70). That is, the Villain approximation amounts to the replacement, . (76)

Wilson loop β-behavior
For definiteness, let us consider a planar Wilson loop C and a planar surface S(C) whose normal points along the1-axis. This surface is placed between the sets of sites {z} and {z+1}, that is, it is crossed by the links that run from z to z+1. Then, we have, At very small β (high "temperatures"), the first contributions to Z XY (α μ ) and Z V (α μ ) coincide, and are of order β 4 . They correspond to loops of length 4, running on the sides of the plaquettes. If none of the loop sides is a link that crosses S(C), then b 1 (z) = 0. If two sides of the plaquette cross S(C), they have different orientations with respect to1, so they do not contribute to z b 1 (z). Only loops with just one side crossing S(C) give a nontrivial factor (77). If M P is the number of links running from z to z +1 and placed on the perimeter of S(C), then, This leads to the average, and where P = M P a is the perimeter of C. As β is increased (the "temperature" decreased), keeping away from β c , the expansion of the partition function will require higher orders in β.
More powers in β imply that larger loops and multiple smaller loops are produced. In the language of superfluids, loops of superflow are generated as the temperature is decreased toward the critical temperature. In our context, more and more center vortices are generated as we approach the continuum limit. Anyway, at any finite order in β a perimeter law (with a renormalized prefactor) is expected. Let us now analyze the situation close to β c . In what follows, we shall denote the total configuration space as {b μ (x)} = {0} ∪ {B 0 }, where {0} and {B 0 } represent the trivial, b μ (x) ≡ 0, and nontrivial configurations, respectively. This is the initial set we shall consider to perform a sequence of approximations. For any subset {A} ⊂ {B 0 }, we define We are interested in computing, The space {B 0 } can be partitioned into three disjoint subsets: where all the loops cross S(C) at least once, and the rest {R 0 }, where configurations contain at least one loop of each type. That is, Close to the transition point, it is well-known that smaller loops with higher fluxes ± 2, ± 3, · · · are irrelevant with respect to larger loops with unit flux, due to the difference in configurational entropy [51]. The replacement in Eq. (76) is exact for configurations {B 0 } characterized by b μ (x) = ±1, or 0 (I −1 (β) = I 1 (β)). This, together with the excellent agreement between the XY and Villain models, indicates that loop configurations that meet at a link are irrelevant with respect to those with single occupation. This refers to loops with the same orientation; those that meet with opposite orientations were already forbidden for non-backtracking loops. This is the way the initial properties of the ensemble are encoded in the statistical properties of the Villain model close to β c . Negative tension and positive stiffness is now related with the preference for larger non-backtracking loops to the detriment of smaller ones. In addition, the statistical irrelevance of multiple occupation or, in other words, excluded volume effects, can be traced back to the repulsive interactions. Then. the calculation can be approximated by Up to a factor a, x,μ |b μ (x)| adds the loop lengths, while This assumption, together with Eq. (84), gives We can proceed with a partition of {B 1 } into three subsets: and so on. To represent the Wilson loop in Eq. (82), each factor must be computed at α μ (with linking numbers) and then divided by the factor at α μ ≡ 0 (without linking numbers). Then, using Eq. (88), we have Next, the single self-avoiding loops in {A 1 } can be partitioned into subsets {z} labelled by the link z, z +1 where the loop intersects S(C) (b 1 (z) = ± 1). Other configurations in {A 1 } are combinations of 2, 3, . . . , Area/a 2 loops in different sets {z}. Then, we estimate, where we have used that for every loop there is a similar one with reversed orientation. Disregarding F {B 2 } , originated from single loops with multiple crossings, and approximating ξ z as z-independent quantities ≈ ξ , we arrive at an area law, where β c ≈ 0.454, M A is the number of links that cross S(C), which span an area A = M A a 2 , and we have used Eq. (67). Then, for SU (2) we get, Finally we notice that the result obtained in Eqs. (90) and (91) can be cast into the form where p z and 1 − 2 p z can be thought of as the probability for the link z, z + 1 to take the value e i 2π N (e −i 2π N ), and 1, respectively.
When N ≥ 3, center vortex ensembles involve branching and the effective model in Eq. (49) (based on a single field V ) is not expected to be applicable. In this respect, for quarks in a representation with N -ality k, the center element generated by center vortices with linking number 1 would be [z(1)] k . Then, if the average were performed along the same lines, with the same ensemble, the result would be, and for large N , and finite k, the string tension ratios would approach a squared sine law, However, this is not the expected behavior for Yang-Mills theory (see Ref. [52] and references therein).

Summary and conclusions
In this paper we have investigated the effective field representation of an ensemble of closed center vortex loops in D = 3 Euclidean spacetime dimensions, as they emerge in the continuum limit after center projection of lattice gauge theory. To account for the properties of center vortices extracted on the lattice, the vortex loops were equipped with tension (μ), stiffness (1/κ) and repulsive contact interactions. Such an ensemble of vortex loops, with definite orientation of the vortex flux, can be mapped onto an effective field theory of a complex scalar field, which has a global U(1) symmetry. We mainly focused on the evaluation of the expectation value of the Wilson loop, i.e. of center elements of the gauge group. For κ, μ > 0 small vortex loops are favoured and a perimeter law is obtained as expected. In this case, the modes of the effective field theory are massive. For κ > 0 and μ < 0, large center vortex loops are favoured and the corresponding effective field theory undergoes the spontaneous breaking of the U(1) symmetry. In this phase, the presence of Goldstone bosons introduces large quantum fluctuations. To deal with them, we formulated the problem on a lattice, using the XY model or its Villain approximation, which have been studied in the context of classical statistical mechanics and superfluids. At high temperatures, the calculation is a perturbative one involving combinatorics. However, to make contact with the continuum, it is important to consider the lattice models at the critical point, placed at β c ≈ 0.454, which has been accessed via Monte Carlo simulations. As this point is approached, the simulations show that more and more loops are generated. The phase transition is driven by single loops becoming infinitely long rather than by the proliferation of multiple smaller loops. In addition, the contribution of configurations with multiply occupied links is numerically irrelevant with respect to that originated from larger loops with single occupation due to the difference in entropy [51]. These properties led to an area law as an extensive property. An interesting point is to clarify at which stage the minimum area appears in the calculation. In this respect, the consideration of loops that intersect the surface S(C) only once seems to work better for the minimum area surface. This is an important ingredient we have used in the derivation. In addition, relevant effects can be incorporated by including non-oriented center vortices with (correlated) magnetic monopoles. No matter how small this component is, it will break the U (1) symmetry explicitly. In the Wilson loop calculation, the implied scale would enable a solitonic-like saddle point, localized on the minimum area surface, plus surface fluctuations that lead to the Lüscher term. Then, while the center vortex loops (oriented or non-oriented) are essential to provide a confining linear potential with N -ality, the correlated monopoles could be relevant to endow this potential with string-like features.