An exact solution of planar QCD by a complex matrix ensemble

Abstract: We construct an exact solution to the planar QCD loop equation in four-dimensional Euclidean space using a novel matrix-valued momentum loop formalism. Central to this construction is a new loop calculus, where functional derivatives act on the loop velocity $\dot C(\theta)$. This framework yields finite, well-defined expressions for point and area derivatives in loop space and highlights the role of the loop-space Bianchi identity in ensuring consistency with gauge theory. The Wilson loop is represented as an average over matrix-valued momentum loops $P(\theta)$ tracing closed paths in a compact complex manifold constrained by self-duality and boundary conditions. These random walks are governed by a measure supported on such manifolds. The self-duality conditions ensure the vanishing of the classical Yang--Mills term in the loop equation, while additional constraints eliminate the contact terms in the planar ($N \to \infty$) limit. A key structural feature of the solution is its invariance under reflection--conjugation symmetry: reversing loop orientation and taking complex conjugation maps a solution to another. This holds for arbitrary compact Lie groups used in the momentum ensemble and is independent of the original SU($N$) gauge group. We conjecture that, in the classical limit, the matrix ensemble defines a generalized minimal surface embedded not in $\mathbb{R}^4$ but in a higher-dimensional matrix space. This (hypothetical) extended surface has self-dual area elements and satisfies the loop equation via geometry and the Bianchi identity.