Energy-energy correlations at next-to-next-to-leading order

TL;DR

This paper presents a novel approach to calculating energy-energy correlations in $ ext{N}=4$ super Yang-Mills theory at NNLO, avoiding infrared divergences and revealing connections to elliptic functions.

Contribution

It introduces a simple formula relating EEC to a triple discontinuity of a four-point correlation function, enabling straightforward numerical evaluation at NNLO.

Findings

Derived a new formula for EEC in $ ext{N}=4$ sYM

Computed EEC at NNLO using a two-fold integral representation

Identified elliptic functions in the integral kernels

Abstract

We develop further an approach to computing energy-energy correlations (EEC) directly from finite correlation functions. In this way, one completely avoids infrared divergences. In maximally supersymmetric Yang-Mills theory ($\mathcal{N}=4$ sYM), we derive a new, extremely simple formula relating the EEC to a triple discontinuity of a four-point correlation function. We use this formula to compute the EEC in $\mathcal{N}=4$ sYM at next-to-next-to-leading order in perturbation theory. Our result is given by a two-fold integral representation that is straightforwardly evaluated numerically. We find that some of the integration kernels are equivalent to those appearing in sunrise Feynman integrals, which evaluate to elliptic functions. Finally, we use the new formula to provide the expansion of the EEC in the back-to-back and collinear limits.

Equations

$$\mathcal{E}(\vec{n})=$$

$$\mathrm{EEC}(\chi )=$$

$$\times \frac{\int d^{4}x{e}^{iqx}⟨0|O^{†}(x)\mathcal{E}({\vec{n}}_1)\mathcal{E}({\vec{n}}_2)O(0)|0⟩}{(q^0{)}^2\int d^{4}x{e}^{iqx}⟨0|O^{†}(x)O(0)|0⟩}.$$

$$\zeta =\frac{q^2(n_1\cdot n_2)}{2(q\cdot n_1)(q\cdot n_2)},$$

$$⟨0|O(x_1)O(x_2)O(x_3)O(x_4)|0⟩=\frac{\mathcal{G}(z,\bar{z})}{(x_{13}^2{x}_{24}^2{)}^2},$$

$$u=z\bar{z}\equiv \frac{x_{12}^2{x}_{34}^2}{x_{13}^2{x}_{24}^2},v=(1-z)(1-\bar{z})\equiv \frac{x_{14}^2{x}_{23}^2}{x_{13}^2{x}_{24}^2}.$$

$${\text{dDisc}}_{w=w_0}g(w)=g(w)-\frac{1}{2}g(w^{\circlearrowleft })-\frac{1}{2}g(w^{\circlearrowright }),$$

$$\mathrm{EEC}(\zeta )$$

$$\times {\text{dDisc}}_{\bar{z}=1}{\text{Disc}}_{z=0}\left[(z-\bar{z}){\mathcal{G}}_{ϵ}(z,\bar{z})\right],$$

$$z=\frac{\zeta t(t-\bar{z})}{t(\zeta -\bar{z})+(1-\zeta )\bar{z}}.$$

$${\mathcal{G}}_{ϵ}(z,\bar{z})\equiv (1-\bar{z}{)}^{ϵ}\mathcal{G}(z,\bar{z}).$$

$$z\bar{z}\Phi \equiv (1-z)(1-\bar{z})(z-\bar{z})\mathcal{G},$$

$$F(\zeta )\equiv 4{\zeta }^2(1-\zeta )\mathrm{EEC}(\zeta ).$$

$$F(\zeta ;a)=a{F}_{\mathrm{LO}}+a^2{F}_{\mathrm{NLO}}+a^3{F}_{\mathrm{NNLO}}+\mathcal{O}(a^4).$$

$${\Phi }^{(1)}=\frac{1}{2}{\mathrm{Li}}_2(\bar{z})-\frac{1}{2}{\mathrm{Li}}_2(z)+\frac{1}{4}\ln \left(\frac{1-\bar{z}}{1-z}\right)\ln \left(z\bar{z}\right).$$

$${\text{Disc}}_{z=0}\left[\frac{{\Phi }^{(1)}}{(1-\bar{z}{)}^{1-ϵ}}\right]=\frac{\pi }{4}\frac{\ln (1-\bar{z})-\ln (1-z)}{(1-\bar{z}{)}^{1-ϵ}}.$$

$${\text{dDisc}}_{w=0}{w}^{-1+ϵ}$$

$$\text{dDisc}[w^{-1+ϵ}{\ln }^{m}w]{\ln }^{n}w$$

$$\displaystyle=\partial^{m}_{\epsilon}\left[2\sin^{2}(\pi\epsilon)\,\partial^{n}_{\epsilon}\left(\frac{\delta(w)}{\epsilon}+\sum_{k=0}\frac{\epsilon^{k}}{k!}\big{[}w^{-1}\ln^{k}w\big{]}_{+}\right)\right]\,.$$

$$F_{\mathrm{LO}}(\zeta )=-\ln (1-\zeta ).$$

$$\left\{\zeta ,1-\zeta ,\frac{1-\sqrt{\zeta }}{1+\sqrt{\zeta }}\right\}.$$

$$F_{\mathrm{NNLO}}(\zeta )$$

$$\times \left[R_1(z,\bar{z})P_1(z,\bar{z})+R_2(z,\bar{z})P_2(z,\bar{z})\right],$$

$$R_1=\frac{z\bar{z}}{1-z-\bar{z}},R_2=\frac{z^2\bar{z}}{(1-z{)}^2(1-z\bar{z})}$$

$$K_1$$

$$K_2$$

$$F_{\text{NNLO}}(\zeta )\stackrel{\zeta \to 1}{\sim }$$

$$-\frac{61}{720}{\pi }^4\ln y-\frac{{\pi }^2}{3}{\zeta }_3-\frac{7}{2}{\zeta }_5.$$

$$F(\zeta )\sim \frac{1}{2}H(a){\int }_0^{\infty }dbb{J}_0(b)e^{-\frac{1}{2}{\Gamma }_{\text{cusp}}(a)L^2-\Gamma (a)L}.$$

$$F_{\mathrm{NNLO}}(\zeta )\stackrel{\zeta \to 0}{\sim }$$

$$+17-\frac{4{\pi }^2}{3}-{\zeta }_3+\frac{5{\pi }^4}{144}+\frac{3}{2}{\zeta }_5].$$

Full text

**1. Introduction. ** The energy-energy correlation (EEC) Basham et al. (1978) measures the energy flow through a pair of detectors separated by an angle $\chi$, see Fig. 1. The EEC has several nice properties, and may help to understand better the nature of jets in quantum field theory. It is an infrared-safe observable Kinoshita (1962); Lee and Nauenberg (1964) that can be computed perturbatively. Moreover, it has simple factorization properties in the back-to-back ($\chi\to\pi$) Collins and Soper (1981); Kodaira and Trentadue (1982); de Florian and Grazzini (2005); Tulipánt et al. (2017); Moult and Zhu (2018) and collinear ($\chi\to 0$) Konishi et al. (1978, 1979) limit. This knowledge can be used to match fixed-order predictions to resummation calculations de Florian and Grazzini (2005); Tulipánt et al. (2017); Moult and Zhu (2018). On the other hand, the EEC is experimentally measurable, and in particular has been used for precision tests of QCD and measurement of the strong coupling constant $\alpha_{s}$ Aaboud et al. (2017); Kardos et al. (2018); Verbytskyi et al. (2019).

The EEC at leading order (LO) is known since Basham et al. (1978), while the next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) results were obtained numerically in refs. Richards et al. (1982, 1983); Ali and Barreiro (1982); Falck and Kramer (1989); Kunszt et al. (1989); Glover and Sutton (1995); Clay and Ellis (1995); Kramer and Spiesberger (1997) and Del Duca et al. (2016a, b), respectively. Only very recently, the NLO result was computed analytically Dixon et al. (2018).

Although the EEC is infrared finite, the standard approach to computing it involves infrared divergent scattering amplitudes Gituliar and Moch (2017); Dixon et al. (2018). On the other hand, it can be defined starting from correlation functions, which are infrared finite Sveshnikov and Tkachov (1996); Korchemsky et al. (1997); Korchemsky and Sterman (1999); Belitsky et al. (2001); Hofman and Maldacena (2008). For example, for $e^{+}+e^{-}\rightarrow\gamma^{*}\rightarrow X$, the main ingredient is a four-point correlation function of two energy-momentum tensors (representing the two detectors), and two electromagnetic currents, which create the electron-positron pair from the vacuum.

To the best of our knowledge, this approach has not yet been implemented in QCD. On the other hand, these ideas were applied in $\mathcal{N}=4$ super Yang-Mills (sYM) Hofman and Maldacena (2008); Belitsky et al. (2014a, b), culminating in the first analytic calculation of and EEC at NLO Belitsky et al. (2014c). The structure of this result, and in particular the types of polylogarithmic functions appearing in it, foreshadowed the structures later found in QCD Dixon et al. (2018).

We show that for an analog of the electromagnetic current in $\mathcal{N}=4$ sYM, the EEC is computed by a new, extremely simple formula, given by a two-fold integral of a particular triple discontinuity of the four-point correlation function, see eq. ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) below.

We validate the new formula ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) by reproducing in an efficient way the known LO and NLO results. We then compute, for the first time, the EEC at NNLO. The result is given in terms of an analytic two-fold integral representation. We present plots of the numerically integrated answer.

We also use the integral formula ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) in order to obtain limits of the energy correlator, namely the small angle and the back-to-back limits.

2. EEC from correlation functions. The detector operator that measures the energy flux in the direction $\vec{n}$ is given by an integrated stress-energy tensor $T_{\mu\nu}$ Sveshnikov and Tkachov (1996); Korchemsky et al. (1997); Korchemsky and Sterman (1999); Belitsky et al. (2001); Hofman and Maldacena (2008),

[TABLE]

The standard definition for the EEC as a differential cross-section can then be recast as

[TABLE]

Here the operators $O$ (source) and $O^{\dagger}$ (sink) create the final state, whose particles are detected by the two calorimeters. Note that the operators are ordered as written, i.e. eq. ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex1)) expresses the EEC in terms of an integrated Wightman four-point correlation function. The choice of the local operator $O$ depends on the physical problem. For $e^{+}e^{-}$ annihilation, $O$ is given by an electromagnetic current.

It is convenient to introduce the variable

[TABLE]

where the null four-vector $n=(1,\vec{n})$ characterizes the position of the detector. In the rest frame of the source $q^{\mu}=(q^{0},\vec{0})$, we have $\zeta=\sin^{2}(\chi/2)$. In the same frame the averaging over the angles in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex1)) becomes trivial for a scalar source.

In ${\cal N}=4$ sYM, we may choose the source and sink to be scalar operators that are the bottom component of the supermultiplet of conserved currents. As such, they are natural analogs of the electromagnetic current, and have fixed conformal weight two. Moreover, the correlation function $\langle{0|OT_{\mu\nu}T_{\rho\sigma}O|0}\rangle$ can be obtained from the all-scalar correlation function $\langle{0|OOOO|0}\rangle$ by superconformal transformations. In this way, we can express the EEC in terms of the four-point correlation function of scalar operators Belitsky et al. (2014b, 2016); Korchemsky and Sokatchev (2015).

Let us describe this relation in detail. We define the Euclidean correlation function

[TABLE]

depending on the conformal cross-ratios

[TABLE]

The analytic continuation from the Euclidean to the Wightman function is easily obtained, for example, by using the Mellin transform of ${\cal G}(u,v)$ Mack (2009); Belitsky et al. (2014a, b, c). Alternatively, we can convert the detector time integration in (1) to integration of certain discontinuities of $\mathcal{G}(u,v)$ Belitsky et al. (2014a). This approach is further improved by exploiting a different double discontinuity, first proposed in Caron-Huot (2017), see also Simmons-Duffin et al. (2018); Kravchuk and Simmons-Duffin (2018) . The latter is defined by

[TABLE]

where $w^{\circlearrowleft}$ and $w^{\circlearrowright}$ refer to the mappings of the branch point $w=w_{0}$ to the two adjacent Riemann sheets. For example, $\text{dDisc}_{w=0}\,w^{\alpha}=2\sin^{2}(\pi\alpha)\,w^{\alpha}$. In addition, in four dimensions and for a source of conformal weight two, the Fourier transform in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex1)) is equivalent to another discontinuity Belitsky et al. (2014a). Putting all this together (for a detailed derivation see Henn et al. (2019)), we arrive at the triple discontinuity formula for the EEC

[TABLE]

where $0<\zeta<1$ 111One expects also contact terms at the endpoint, which are required e.g. to verify the relation $\int_{-1}^{1}\rm EEC(\chi)d\cos(\chi)=1$ Basham et al. (1978). In this letter, we assume the operators to be neither coincident nor back-to-back.. After taking the discontinuities, the variable $z$ is treated as a function of the two integration variables,

[TABLE]

The following comments are in order. In order to apply dDisc to power singularities, we employ an analytic regulator for the poles in $1-\bar{z}$,

[TABLE]

Instead of first computing the integrals in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) and taking the regulator $\epsilon\to 0^{+}$ at the end, it is much more efficient to treat the integrand as a singular distribution Gelfand and Shilov (1964) of the ‘plus’ type, $w^{-1+\epsilon}_{+}=\frac{1}{\epsilon}\delta(w)+w_{+}^{-1}+\epsilon[w^{-1}\ln w]_{+}+O(\epsilon^{2})$. This allows us to drop many irrelevant terms, as illustrated in the one-loop example below. Moreover, it allows us to automate the calculation at higher loops.

For convenience, we define

[TABLE]

and

[TABLE]

The correlation function ${\cal G}$ has the perturbative expansion ${\cal G}=\sum_{k\geq 0}a^{k}{\cal G}^{(k)}$ in the ‘t Hooft coupling $a=N_{c}g^{2}_{\rm{YM}}/(4\pi^{2})$. It was computed to two loops in Eden et al. (2000); Bianchi et al. (2000), and at three loops in Eden et al. (2012); Drummond et al. (2013). Note that to obtain the EEC at NkLO, one needs ${\cal G}$ at $(k+1)$ loops. For any $\zeta\in(0,1)$, we write the expansion

[TABLE]

3. EEC at LO. At Born level $\Phi^{(0)}\sim z-\bar{z},$ and our triple discontinuity formula ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) gives 0, which is the expected answer (up to contact terms). Let us now apply eq. ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) to the leading order (or one-loop) EEC, for which Usyukina and Davydychev (1993)

[TABLE]

The discontinuity at $z=0$ is

[TABLE]

The following operation $\text{dDisc}_{\bar{z}=1}$ is carried out keeping in mind the singular character of the distribution $w^{-1+\epsilon}$,

[TABLE]

This expression vanishes in the limit $\epsilon\to 0$. A non-trivial result is obtained due to the logarithms, with the help of the identity

[TABLE]

One such logarithm is already present in eq. (14). The $t-$integration in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) gives $A\ln(1-\bar{z})+B$ with coefficients $A,B$ analytic at $\bar{z}=1$. At LO we apply (16) with $m,n\leq 1$. The operation dDisc produces $\delta(1-\bar{z})$, i.e. we observe a transcendental weight drop of three. The $\bar{z}$ integral in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) is removed and we obtain, in accord with Engelund and Roiban (2013),

[TABLE]

4. EEC at NLO and NNLO. The example above shows the way to higher perturbative orders.

At NLO, $\Phi^{(2)}$ is given by so-called two-loop ladder and one-loop box squared integrals Eden et al. (2000); Bianchi et al. (2000), expressed in terms of weight-four harmonic polylogarithms (HPL) Remiddi and Vermaseren (2000) in $z$ and $\bar{z}$. The HPL shuffle algebra allows us to extract all singular terms as $z\to 0$ and $\bar{z}\to 1$. This is conveniently implemented in the Mathematica package HPL Maitre (2006). In this way, we obtain powers of $\ln z$ and $\ln(1-\bar{z})$, whose coefficients are HPL that are analytic at $z=0$ and $\bar{z}=1$. The triple discontinuity of the logarithms lowers the transcendental weight by three. We are left with a two-fold integration in $t$ and $\bar{z}$. We find that the integrand is linearly reducible Brown (2009a, b); Bogner and Lüders (2015); Panzer (2015a). (In the case of the one-loop box squared, we first change variables to achieve this.) As a consequence, we can carry out all integrals algorithmically Bogner and Brown (2015); Panzer (2015b). This is conveniently done using the Maple program HyperInt Panzer (2015b). The final result for the EEC can be converted to classical polylogarithms and fully agrees with Belitsky et al. (2014c). The symbol alphabet Goncharov et al. (2010) is

[TABLE]

At NNLO, the integral formula ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) exhibits some new algebraic features, due to the leading singularities (algebraic prefactors) of the three-loop integrals appearing in $\Phi^{(3)}$ and their symbol alphabet Drummond et al. (2013).

The set of integrals are three-loop ladders, products of one- and two-loop boxes, as well as so-called easy and hard integrals that come in three different orientations. A new feature compared to two loops is that the hard integral is given in terms of a class of Goncharov polylogarithms (GPL) Goncharov et al. (2010), whose symbol contains an entry $z-\bar{z}$. We implement the GPL shuffle relations to extract the singularities at $z=0$ and $\bar{z}=1$, and take the triple discontinuity. We find that all except two terms are linearly reducible and can be integrated, giving HPL with transcendental weight ranging from two to five. The symbol alphabet of these terms is the same as at NLO.

For certain orientations, the easy and hard integrals have leading singularities that introduce a rational factor $1/(1-u)$ or $u/(u-v)$. These terms do not admit a change of integration variables that allows the two-fold integration in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) to be automated in HyperInt. We present our final three-loop result as a sum of polylogarithms $f_{\rm HPL}$, plus a two-fold finite integral,

[TABLE]

where

[TABLE]

and $z$ is given by eq. (8). $P_{1},P_{2}$ are HPL in $z,\bar{z}$ of weight three. We provide their explicit expressions, as well as the polylogarithm function $f_{\rm HPL}(\zeta)$ with symbol alphabet (18), in the ancillary files.

After partial fractioning, $R_{1},R_{2}$ decompose into a sum of linearly reducible terms and irreducible kernels $K_{1},K_{2}$,

[TABLE]

Remarkably, upon the replacements $\zeta\rightarrow\frac{m^{2}}{m^{2}-S}$, $t\rightarrow x_{1}+1,\,\bar{z}\rightarrow x_{2}$, the denominator of $K_{1}$ becomes the Symanzik polynomial of the equal mass sunrise integral Broadhurst (1990); Laporta and Remiddi (2005); Bloch and Vanhove (2015); Adams et al. (2016a). The features of the integrand that we observe here is evidence that the integral is elliptic.

We evaluate numerically the two-fold integrals in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)) in Mathematica for 50 different values of $\zeta$ with precision $O(10^{-7})$. In Fig. 2 we plot the EEC shape at the first three perturbative orders.

5.  Back-to-back and collinear expansion. The EEC is logarithmically enhanced in the extreme kinematic situations where the two detectors are back-to-back ($\zeta\rightarrow 1$) or are close to each other ($\zeta\rightarrow 0$). We extract analytically the contribution of the two-fold integral in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)) in both limits.

In the back-to-back limit, we observe that as $\zeta\rightarrow 1$ the integrand as a function of $t$ and $\bar{z}$ is suppressed everywhere in the integration domain, hence the integral in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)) only contributes subleading powers of $y=1-\zeta$. The leading asymptotic behavior of $F_{\text{NNLO}}(\zeta)$ is determined by $f_{\rm HPL}(\zeta)$. At leading power in $y$,

[TABLE]

On the other hand, the $\zeta\rightarrow 1$ behavior of the EEC in $\mathcal{N}=4$ sYM is analogous to that in QCD. The Sudakov logarithms can be resummed to all loop orders Collins and Soper (1981),

[TABLE]

Here $L=\ln(e^{2\gamma_{E}}b^{2}/4y)$; $J_{0}$ is a Bessel function; $\Gamma_{\text{cusp}}(a)$ and $\Gamma(a)$ are the cusp and collinear anomalous dimensions, whose three-loop values can be found in Bern et al. (2005). The function $H(a)$ describes hard emissions in the EEC.

We verify that our expression ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex5)) agrees with the resummation formula (24) expanded to three-loop order. The coefficient of the single logarithm term allows us to determine the hard function up to NNLO, $H(a)=1-\zeta_{2}\,a+5\zeta_{4}\,a^{2}$.222 The hard function $H(a)$ was independently obtained at two loops in Korchemsky (2019).

Remarkably, eq. ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex5)) has homogeneous transcendental weight, and so do the ingredients of the resummation formula. In fact, in the back-to-back limit, we anticipate that $\mathcal{N}=4$ sYM describes the maximally transcendental part of the QCD asymptotics Belitsky et al. (2014c). This is confirmed by comparing with the recent QCD result Moult and Zhu (2018) upon the substitution $C_{F},C_{A}\rightarrow N_{c}$.

In the collinear limit, the elliptic integral in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)) is not suppressed. We extract its leading contribution by expanding in powers of $\zeta$ at the integrand level and integrating the leading term, which comes solely from the $R_{1}\times P_{1}$ piece. Combining with $f_{\rm HPL}(\zeta)$ in the $\zeta\rightarrow 0$ limit, we obtain the three-loop small-angle asymptotics

[TABLE]

A very nontrivial confirmation of our result is the agreement with the prediction of the light-ray OPE approach Koloğlu et al. (2019).

In Fig. 3, we show the comparison between the leading power asymptotics and the numerical evaluation of ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)).

6. Outlook. Our work opens the way for several applications. Our main formula ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) expresses the EEC in terms of a four-point function, a central object in any conformal field theory (CFT). The derivation of ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)), to be given in Henn et al. (2019), relied on the fact that the operator ${O}$ has fixed scaling dimension two. It would be interesting to see if a similar triple discontinuity formula exists in a generic CFT, where the source and sink are replaced by a conserved current. It would also be very interesting to explore how our formalism works in QCD.

We extracted the leading power asymptotics in the collinear and back-to-back limits at NNLO. We wish to emphasize that our result ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)) contains information about subleading powers as well. The latter constitute useful data for resummations of large logarithms at subleading power Laenen et al. (2008, 2011); Bonocore et al. (2016); Boughezal et al. (2017); Moult et al. (2017a); Feige et al. (2017); Moult et al. (2017b, 2018a, 2018b); Balitsky and Tarasov (2017); Beneke et al. (2018a, b). Moreover, thanks to eq. ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) it may be possible to understand these limits, at arbitrary coupling, in terms of properties of the four-point correlation function.

In $\mathcal{N}=4$ sYM, the integrand of the four-point function ${\cal G}$ is known to ten loops Eden et al. (2012); Bourjaily et al. (2016). Once the integrated four-loop result becomes available Drummond et al. (2007); Drummond (2013); Schnetz (2014); Golz et al. (2017); Caron-Huot and Henn (2014); Eden and Smirnov (2016), our formula ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex2)) can be used to obtain the EEC at N3LO.

We provided evidence that the EEC at NNLO contains elliptic functions. It suggests that the same type of function might also appear in QCD. It would be interesting to fully work out the remaining integrals in ([

$\scriptstyle\rm CERN-TH-2019-026,\ LAPTH-014/19,\ MPP-2019-54$

Energy-energy correlations at next-to-next-to-leading order](#S0.Ex4)), using multiple elliptic polylogarithms Adams et al. (2016b); Adams and Weinzierl (2018); Remiddi and Tancredi (2017); Ablinger et al. (2018); Broedel et al. (2018, 2019a, 2019b).

Finally, it would be interesting to apply our infrared finite approach to energy correlators with more than two detectors (cf. Moult et al. (2016); Komiske et al. (2018); Chivukula et al. (2018) for related generalizations of energy correlations). This may be important for the theoretical understanding of jet observables Coleman et al. (2018); Larkoski et al. (2017); Asquith et al. (2018); Komiske et al. (2019); Marzani et al. (2019); Datta et al. (2019).

Acknowledgements. We would like to thank D. Chicherin for his generous help with the calculations and T. Hahn, M. D. Schwartz and G. Zanderighi for helpful conversations. We are grateful to the authors of Ref. Drummond et al. (2013), as well as to Ian Moult and Tong-Zhi Yang for useful correspondence. We acknowledge many enlightening discussions with G. Korchemsky. This work received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme, Novel structures in scattering amplitudes (grant agreement No 725110). J. M. H., E. S.  and K. Y. also wish to thank the Galileo Galilei Institute for hospitality during the workshop “Amplitudes in the LHC era”. A. Z. is grateful to M. Kologlu, P. Kravchuk, D. Simmons-Duffin for many discussions on the subject.