⚛️ high-energy theory

Chaos in Yang-Mills

This paper investigates the chaotic dynamics of pure-gauge SU(N) Yang-Mills theory by analyzing Wilson loop out-of-time-order correlators, finding that chaos is restricted to non-abelian theories in dimensions greater than two and proposing a critical loop size that determines whether scrambling occurs due to infrared sensitivity at the magnetic scale.

Original authors: Barel Skuratovsky

Published 2026-08-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Barel Skuratovsky

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast landscape of modern physics, there is a fundamental question about how complex systems lose their order and become unpredictable. This process, known as scrambling, is central to understanding everything from the behavior of hot plasma to the mysterious nature of black holes. Scientists study this by watching how tiny disturbances spread through a system over time. If a system is chaotic, a small change at the beginning will quickly ripple out, making the entire system look completely different and impossible to predict. To measure this, researchers look at specific mathematical quantities that grow exponentially when chaos is present. The speed of this growth is called the Lyapunov exponent, a number that tells us how fast a system forgets its initial state.

The paper at hand investigates this phenomenon within pure Yang–Mills theory, a mathematical framework that describes the strong force holding atomic nuclei together. This theory involves fields that carry a property called color charge, which comes in different types. The researchers focused on a specific object called a Wilson loop, which is essentially a closed path drawn through space. In the language of this theory, such a loop acts as a probe, measuring the properties of the force field along its path. The team wanted to know how this loop changes as time passes and whether it exhibits the rapid, chaotic growth characteristic of scrambling. They were particularly interested in whether the size of the loop mattered, and under what conditions the theory becomes chaotic at all.

The researchers began by asking a simple but profound question: how does a Wilson loop evolve in time? They discovered that the answer depends entirely on two factors: the symmetry group of the theory and the number of spatial dimensions. They found that for a loop to grow and scramble, the underlying theory must be non-abelian, meaning the order in which operations are performed matters, and the universe must have more than two spatial dimensions. If the theory were abelian, like the simpler theory of electromagnetism, or if it existed in only two dimensions, the loop would simply slide along without changing its fundamental structure. In those cases, the system does not scramble. The author concludes that pure gauge theories are chaotic only when they are non-abelian and exist in more than two dimensions.

Moving to the conditions of weak coupling, where the interactions between particles are relatively feeble, the team performed detailed calculations to see how fast this scrambling happens. They treated the loop as a collection of gluons, the particles that carry the strong force, and tracked how they interacted over time. Their calculations revealed a surprising sensitivity to the size of the loop. They found that the rate of chaos is not a fixed number but depends on the loop's diameter relative to a specific scale set by the surrounding plasma. This scale is related to how far the plasma can screen electric charges, a distance known as the Debye length.

The most striking result emerged when they compared the loop size to this screening distance. They found a critical threshold. If the loop is smaller than this critical diameter, it scrambles rapidly, behaving chaotically. However, if the loop is larger than this threshold, the scrambling stops, and the system decays instead. At the exact point where the loop size matches this critical diameter, the time it takes for the system to scramble becomes infinite. This suggests a continuous transition where the behavior of the system changes smoothly from chaotic to non-chaotic as the loop grows. The researchers propose that this transition is driven by the fact that a large loop cannot "see" or interact with very soft, low-energy fluctuations in the plasma that are smaller than its own size.

While the calculations were performed in a regime where the interactions are weak, the author also considered what might happen when the interactions are strong, a condition more typical of the real world inside atomic nuclei. In this strong-coupling regime, the particles are no longer the primary actors; instead, the theory is better described by strings or surfaces. The author suggests that the evolution of the loop can be visualized as a surface stretching between the loop's position at the start and its position at the end. If the loop splits and rejoins, the surface develops a handle, much like a tube turning into a pair of pants. They propose that the probability of such events is governed by the topology of these surfaces, but they note that calculating the exact rate of chaos in this strong regime remains a difficult challenge, as the necessary mathematical tools are not yet fully developed.

The study concludes by highlighting that the chaotic behavior of these loops is deeply tied to the infrared, or long-distance, properties of the theory. Even in the weak-coupling limit, the rate of chaos is influenced by the non-perturbative magnetic scale, a region where the theory becomes strongly coupled and standard calculation methods break down. By using the loop's own size as a physical cutoff, the researchers were able to map out a transition that links the microscopic chaos of the theory to the macroscopic size of the probe. This work provides a concrete framework for understanding how chaos emerges in gauge theories and suggests that the size of an observer is a crucial factor in determining whether a system appears chaotic or stable. The findings offer a new perspective on the conditions required for chaos, pointing to a delicate balance between the geometry of the observer and the screening properties of the environment.

Technical Summary: Chaos in Yang–Mills

Problem Statement
This work investigates the chaotic dynamics of pure-gauge $SU(N)$ Yang–Mills theory by analyzing the time evolution of spatial Wilson loops. The primary diagnostic tool is the out-of-time-order correlator (OTOC), defined as the thermal expectation value of the squared commutator of a Wilson loop with itself at two different times. While OTOCs have been established as measures of chaos in many-body systems and calculated in weakly coupled large- $N$ field theories, their application to non-local, gauge-invariant operators like Wilson loops in pure gauge theory remains underexplored. The paper seeks to determine the conditions under which these operators grow (scramble) and to compute the associated Lyapunov exponent, $\lambda_L$ , particularly addressing the infrared (IR) sensitivity inherent in thermal gauge theories.

Methodology
The analysis is divided into weak and strong coupling regimes, with a primary focus on the weak coupling limit ( $\lambda = g^2N \ll 1$ ) in $3+1$ dimensions at large $N$ .

Operator Dynamics (Section 2):
The author derives the time evolution of a spatial Wilson loop using the Heisenberg equation. By expanding the loop in nested commutators with the Hamiltonian, they identify two distinct modes:
- Displacement: A rigid translation of the contour driven by the electric field ( $F_{0i}$ ).
- Splitting: A source term driven by the magnetic field ( $D_j F_{ji}$ ) that inserts field strength at a point on the contour, effectively cutting the loop.
  They demonstrate that the splitting mechanism, which is necessary for operator growth, relies on the non-abelian structure constants ( $f^{abc}$ ) and the existence of at least two spatial dimensions ( $d > 2$ ). Consequently, pure gauge theories in $d=2$ or abelian groups exhibit no operator growth.
Weak Coupling Computation (Section 3):
The OTOC is expanded to leading order, reducing the problem to a gluon four-point function. The author identifies 36 surviving Wick contractions of the eight gauge fields involved, which fall into three topological classes:
- Class (i): Each commutator contracts with itself (four "rails").
- Class (ii): One line per commutator spans the time interval (two rails), with the remaining fields forming equal-time "caps."
- Class (iii): Crossed topology (discarded as it lacks a fixed ladder skeleton).
The author constructs a Bethe-Salpeter equation (BSE) by resumming ladder diagrams formed by "rungs" (exchanged gluons) between the rails.
- Color Algebra: In Class (i), a single exchanged gluon between commutators vanishes due to the trace of a single adjoint generator ( $\text{Tr}(F^e)=0$ ), forcing the leading rung to be $O(\lambda^2)$ . In Class (ii), a single gluon exchange is allowed, entering at $O(\lambda)$ , competing directly with the damping of the rails.
- Numerical Solution: The author solves the BSE for Class (ii) numerically. They place the "free blocks" (the initial four-point function) on-shell and treat the rungs as operators acting on the eigenfunctions.
- IR Cutoff: Recognizing that the Lyapunov exponent is IR-sensitive due to the non-perturbative magnetic scale in the gluon self-energy, the author proposes a physical cutoff. They argue that the Wilson loop, being a color singlet of diameter $2R$ , is blind to momentum exchanges softer than $1/2R$ . This scale is substituted for the magnetic scale in the damping rate calculation.
Strong Coupling (Section 4):
At strong coupling, where gluons are not the fundamental excitations, the author proposes a "worldsheet" description. The evolution is modeled as a sum over cobordisms (surfaces bounded by the initial and final loops) graded by genus. The OTOC is interpreted as a $2 \to 2$ closed-string scattering amplitude. However, the author notes that a controlled calculation of the scattering amplitude in a confining background is currently unavailable.

Key Results

Necessity of Non-Abelian $d>2$ : The paper establishes that operator growth in pure gauge theory requires a non-abelian group and $d > 2$ . The author conjectures that pure gauge theories are chaotic only when non-abelian at $d > 2$ , as the splitting mechanism vanishes in $d=2$ or for abelian groups where operators merely displace without scrambling.
Lyapunov Exponent and IR Sensitivity: Even at weak coupling, the Lyapunov exponent $\lambda_L$ is sensitive to the IR magnetic scale. Without an IR regulator, the exponent cannot be computed perturbatively.
Loop-Size Transition: By imposing the loop diameter $2R$ as the IR cutoff (replacing the magnetic scale in the damping rate), the author finds a critical transition in loop size.
- For loops smaller than a critical diameter $R^*$ , the exponent is positive ( $\lambda_L > 0$ ), indicating scrambling.
- For loops larger than $R^*$ , the exponent becomes negative, indicating decay dominated by the damping of constituent gluons.
- The critical diameter scales as $2R^* m_D \simeq 0.56 \lambda^{-1/2}$ , where $m_D$ is the Debye mass. This implies the critical size tracks the magnetic screening length rather than the Debye length.
- At the transition point, the scrambling time diverges, characteristic of a continuous dynamical transition.
Scaling of the OTOC: The author calculates the prefactor for the OTOC of a normalized Wilson loop, finding it scales as $O(N^{-4})$ rather than the generic $O(N^{-2})$ for single-trace operators, due to the specific color counting of the loop expansion.

Significance and Claims
The paper claims to provide the first framework for describing chaos in pure-gauge Yang–Mills via Wilson loops. Its primary contributions are:

Theoretical Constraint: It rigorously proves that operator growth in pure gauge theory requires non-abelian symmetry and $d>2$ , and conjectures that this condition is necessary for chaos in such theories, linking the existence of propagating gluons and the splitting mechanism.
Size-Dependent Chaos: It proposes that chaos in gauge theories is not a universal property of the theory alone but depends on the scale of the probe. There exists a "scrambling threshold" determined by the loop size relative to the plasma's screening lengths.
Strong Coupling Proposal: It offers a geometric interpretation of the OTOC at strong coupling as a sum over cobordisms, identifying the obstruction to calculation as the need for a regularized scattering amplitude above the effective string cutoff.

Limitations and Open Questions
The author explicitly states several limitations:

The calculation is restricted to Class (ii) topologies; the dominant Class (i) channel (by color factors) was not solved due to complexity.
The IR cutoff via loop size is a physical ansatz; a first-principles derivation of the thermal decay width in the presence of the loop is required to confirm the transition is physical.
The results are derived in the deconfined phase; the behavior of the exponent across the confinement-deconfinement transition remains an open question.
The "butterfly velocity" and momentum dependence of the exponent were not computed (zero injected momentum was assumed).

The paper concludes that while the transition appears robust within the proposed framework, establishing its physical reality requires further investigation into the thermal decay width and the non-perturbative magnetic sector.

More like this