⚛️ nuclear theory

The hidden sunrise in the energy-energy correlator

This paper presents a complete analytic computation of the energy-energy correlator (EEC) in $\mathcal{N}=4$ super Yang-Mills theory at next-to-next-to-leading order by resolving the elliptic sector through its identification with the sunrise integral's modular curve and expressing the result in terms of iterated Eisenstein integrals.

Original authors: Matthew D. Schwartz, Xiaoyuan Zhang

Published 2026-10-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Matthew D. Schwartz, Xiaoyuan Zhang

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 subatomic world, particles do not simply bounce off one another like billiard balls; they explode into showers of new particles, spraying energy in every direction. To understand the rules governing these violent collisions, physicists measure how that energy is distributed. One particularly powerful tool is the energy-energy correlator, a measurement that asks a simple question: if you detect a particle moving in one direction, how likely is it that you will find another particle moving at a specific angle away from it? This measurement acts like a sensitive probe, revealing the hidden geometry of the forces that hold matter together. While some aspects of these collisions can be predicted with standard mathematical tools, the most complex interactions often require a more sophisticated language, one that involves shapes and patterns far beyond simple curves.

For decades, physicists have been able to calculate these energy patterns with increasing precision, but a specific, stubborn piece of the puzzle remained unsolved. In the theoretical framework known as N = 4 super Yang-Mills theory—a simplified model used to test the limits of our understanding of particle physics—the calculation for a high level of precision had stalled. The result was almost complete, consisting of a large collection of known mathematical functions, but it was missing a final, crucial component. This missing piece was a two-dimensional integral that resisted all attempts to be written down in a closed form. It was a barrier that prevented scientists from seeing the full picture of how energy flows in these idealized collisions.

In this work, the researchers finally broke through that barrier. They discovered that the stubborn, unsolved integral was not a random mathematical mess, but was deeply connected to a famous geometric shape known as the sunrise curve. This shape, which appears in the study of particle loops, is defined by a specific type of elliptic curve, a geometric object that is more complex than a circle or a parabola. By recognizing this connection, the authors were able to map the energy correlator onto this known curve. They found that the two shapes shared an identical mathematical fingerprint, allowing them to translate the difficult integral into a new language of modular forms. These are special functions that repeat in a structured way, much like a pattern on a wallpaper, but with a complexity that captures the behavior of the elliptic curve.

The team did not just identify the shape; they solved the problem completely. They expressed the final answer using a set of tools called iterated Eisenstein integrals. These are mathematical objects that can be evaluated with extreme speed and precision, turning a calculation that once seemed intractable into something that can be computed in seconds. This achievement is significant because it provides the first complete, analytic description of the energy-energy correlator at this level of precision. It means that for the first time, physicists have a precise formula that describes how energy is distributed in these collisions, including the most complex, elliptic parts that were previously hidden.

Beyond solving the equation, the researchers used this new result to test a different approach to physics called the bootstrap method. This method attempts to build physical predictions from scratch, using only the known rules of symmetry and the behavior of particles at extreme limits, without doing the heavy lifting of the original calculation. The team asked: if they had not known the answer beforehand, could they have deduced it just by knowing the rules and checking a few numbers? They found that the rules of physics alone were not enough to pin down the answer. While the physical constraints eliminated many wrong possibilities, they left thousands of unknown numbers still floating in the equation. To find the correct values, they had to rely on high-precision numerical data. This revealed a vital lesson for the future: while the geometric rules of the universe are powerful, they are not always sufficient on their own. To fully understand the complex machinery of particle collisions, we must combine the elegance of mathematical rules with the brute force of precise numerical measurement.

The implications of this work extend beyond this specific theoretical model. The techniques developed here, particularly the ability to handle elliptic curves and evaluate them quickly, offer a new path forward for understanding real-world particle collisions in quantum chromodynamics, the theory that describes the strong force holding atomic nuclei together. By mastering these complex mathematical structures in a simplified setting, the researchers have provided a blueprint for tackling the even more difficult problems in the real world. They have shown that the hidden sunrise of the energy correlator can be brought into the light, offering a clearer view of the fundamental forces that shape our universe.

Problem Statement
The Energy-Energy Correlator (EEC) is a collider observable that measures energy flow correlations and is infrared finite, making it a prime candidate for high-order perturbative calculations. While the EEC in $\mathcal{N}=4$ Super Yang-Mills (SYM) theory was previously computed to next-to-next-to-leading order (NNLO) in Ref. [1], the result was incomplete. It consisted of harmonic polylogarithms (HPLs) up to weight five plus a remaining unevaluated two-fold integral containing elliptic curves. Unlike polylogarithmic integrals, which can be systematically solved using canonical differential equations and expressed in terms of generalized polylogarithms, elliptic Feynman integrals pose significant challenges. There is no general algorithmic method to express such results in a universal basis, and numerical evaluation to high precision is difficult. This paper addresses the completion of the NNLO EEC calculation in $\mathcal{N}=4$ SYM by evaluating the remaining elliptic sector and expressing the full result in a form amenable to high-precision numerical evaluation.

Methodology
The authors employ a multi-step analytical and numerical strategy:

Partial Fraction Decomposition and Galois Symmetry: Starting from the unevaluated two-fold integral $F_R(\zeta)$ , the authors perform a partial fraction decomposition on the rational coefficients of the integrand. Guided by Galois symmetry, they separate the integral into a polylogarithmic piece ( $F_{red1}$ ) and two irreducible pieces ( $F_{K1}$ and $F_{R2}$ ) that carry square roots of quartic polynomials.
Identification of the Elliptic Curve: The authors demonstrate that the two irreducible pieces, $F_{K1}$ and $F_{R2}$ , live on the same elliptic curve. By comparing their $j$ -invariants and applying a Möbius transformation, they identify this curve as the equal-mass sunrise curve, a well-known two-loop Feynman integral. They show that the difference $\Delta(\zeta) = F_{R2} - F_{K1}$ is purely polylogarithmic.
Differential Equations and Modular Forms: The authors derive a system of coupled second-order differential equations for the elliptic integrals. They identify the Picard-Fuchs operator governing the system as that of the sunrise integral. By moving to a modular frame defined by the modular group $\Gamma_1(6)$ , they express the solution in terms of iterated Eisenstein integrals. The variable $\zeta$ is mapped to a modular parameter $q$ via a mirror map involving Dedekind eta functions.
Analytic Regression and Bootstrap: With the full analytic form derived, the authors initiate a "Landau bootstrap" study. They construct a function space based on the identified alphabets (polylogarithmic and elliptic) and propose an ansatz for the rational coefficients. They apply physical constraints (regularity, leading power behavior in collinear and back-to-back limits, and sum rules) to fix a fraction of the coefficients. The remaining undetermined parameters are fixed using analytic regression (lattice reduction) against high-precision numerical data generated from the analytic result itself.

Key Contributions and Results

Complete Analytic Result: The paper provides the first complete analytic expression for the NNLO EEC in $\mathcal{N}=4$ $N = 4$ SYM. The result is decomposed into a polylogarithmic sector ( $F_P$ $F_{P}$ ) and an elliptic sector ( $F_E$ $F_{E}$ ).
- The polylogarithmic sector is expressed in terms of HPLs of weight up to five.
- The elliptic sector is expressed as a linear combination of iterated Eisenstein integrals of $\Gamma_1(6)$ , multiplied by rational functions and the modular form $\psi_1$ .
Numerical Evaluation: The authors implement the evaluation of iterated Eisenstein integrals as truncated $q$ -series. This method allows the NNLO distribution to be evaluated to high precision (e.g., 50+ digits) within seconds, overcoming previous difficulties in evaluating elliptic functions.
Singular Expansions: The authors derive singular expansions for the EEC in both the collinear ( $\zeta \to 0$ ) and back-to-back ( $\zeta \to 1$ ) limits to high orders. These expansions reveal the structure of large logarithms and constants (including $\log 2$ , polylogarithms, and specific elliptic periods like $S_0(1)$ ) required for resummation.
Bootstrap Analysis: The paper presents a first study of the perturbative bootstrap for the EEC. It demonstrates that while physical constraints (regularity, OPE limits, sum rules) fix a significant portion of the function space, they are insufficient to determine the rational coefficients uniquely. The authors show that analytic regression using high-precision numerical data can successfully reconstruct the full analytic form, fixing the remaining parameters.
Technical Implementation: The authors note that many technical steps in the derivation and the bootstrap regression were completed with the assistance of AI under human supervision, utilizing tools like the BootLoops repository.

Significance
The paper claims significance in several areas:

Completing the NNLO Calculation: It resolves the "hidden sunrise" in the EEC, providing a closed-form analytic result that was previously inaccessible.
Elliptic Function Space: It establishes that the elliptic sector of the EEC is governed by the geometry of the sunrise curve and the modular group $\Gamma_1(6)$ . This provides a concrete handle on the elliptic sector of the EEC in QCD, suggesting that similar geometric analyses could simplify future QCD calculations.
Bootstrap Methodology: The study serves as a proof-of-concept for applying bootstrap methods to collider observables with elliptic sectors. It highlights the necessity of combining physical constraints with analytic regression from numerical data to solve for rational coefficients in complex observables.
Phenomenological Utility: The ability to evaluate the result to high precision quickly via $q$ -series makes the NNLO EEC a viable tool for precision QCD studies, such as the measurement of the strong coupling constant, once the QCD result is obtained.

The authors remain modest regarding the immediate application to QCD, noting that while the function space is fixed by geometry, the specific integrand details for QCD remain to be analyzed. However, they posit that knowing the function space in advance can bypass intermediate steps in direct calculations and simplify results.

More like this