In this review, we discuss the present status of the description of confining flux tubes in SU(N) pure Yang–Mills theory in terms of ensembles of percolating center vortices. This is based on three main pillars: modeling in the continuum the ensemble components detected in the lattice, the derivation of effective field representations, and contrasting the associated properties with Monte Carlo lattice results. The integration of the present knowledge about these points is essential to get closer to a unified physical picture for confinement. Here, we shall emphasize the last advances, which point to the importance of including the non-oriented center-vortex component and non-Abelian degrees of freedom when modeling the center-vortex ensemble measure. These inputs are responsible for the emergence of topological solitons and the possibility of accommodating the asymptotic scaling properties of the confining string tension.

Recently, based on a new procedure to quantize the theory in the continuum, it was argued that Singer’s theorem points toward the existence of a Yang-Mills ensemble. In the new approach, the gauge fields are mapped into an auxiliary field space used to initially determine sectors labeled by center vortices, and then separately fix the gauge on them. In this work, we study this procedure in more detail. We provide examples of configurations belonging to these sectors and discuss the existence of non-Abelian degrees of freedom. Then, we discuss the importance of the mapping injectivity, and show that this property holds infinitesimally for typical configurations of the vortex-free sector and for the simplest example in the one-vortex sector. Finally, we show that these examples are free from Gribov copies.

In this work, we propose a $3D$ ensemble measure for center-vortex worldlines and chains equipped with non-Abelian degrees of freedom. We derive an effective field description for the center-element average where the vortices get represented by $N$ flavors of effective Higgs fields transforming in the fundamental representation. This field content is required to accommodate fusion rules where $N$ vortices can be created out of the vacuum. The inclusion of the chain sector, formed by center-vortex worldlines attached to pointlike defects, leads to a discrete set of $Z(N)$ vacua. This type of SSB pattern supports the formation of a stable domain wall between quarks, thus accommodating not only a linear potential but also the L\"uscher term. Moreover, after a detailed analysis of the associated field equations, the asymptotic string tension turns out to scale with the quadratic Casimir of the antisymmetric quark representation. These behaviors reproduce those derived from Monte Carlo simulations in $SU(N)$ $3D$ Yang-Mills theory, which lacked understanding in the framework of confinement as due to percolating magnetic defects.