⚛️ high-energy theory

Analytic structure of holographic thermal correlators from Fourier series

This paper computes holographic Euclidean thermal correlators for scalar operators using a distributional Fourier series on the thermal circle, yielding a manifestly periodic result that recovers all OPE coefficients (including the double-trace sector) and reveals bouncing singularities as non-perturbative transseries sectors with vanishing parameters.

Original authors: Paolo Arnaudo, Benjamin Withers

Published 2026-03-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Paolo Arnaudo, Benjamin Withers

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

Imagine you are trying to understand the weather inside a giant, invisible storm cloud (a black hole) by listening to the sound of rain hitting the ground outside. In the world of theoretical physics, this "storm" is a black hole, and the "rain" is the heat and energy of the universe around it.

This paper by Paolo Arnaudo and Benjamin Withers is about finding a new, clever way to listen to that rain to understand the storm's secrets, specifically how particles interact when the universe is hot.

Here is the breakdown of their work using simple analogies:

1. The Problem: Listening to a Storm

In physics, scientists use a tool called Holography. Think of it like a 3D movie projector. The "movie" (our 3D universe with gravity and black holes) is actually a projection of a "film strip" (a 2D quantum world without gravity).

The authors wanted to calculate how two particles "talk" to each other in a hot universe (a thermal state). Usually, to do this, they have to solve very difficult math equations (differential equations) that describe how waves move near a black hole. It's like trying to predict the exact path of every single raindrop in a hurricane. It's messy, hard, and often requires brute-force computer simulations.

2. The New Tool: The Fourier Series (The Musical Score)

Instead of tracking every raindrop, the authors decided to look at the music of the storm.

They used a mathematical tool called a Fourier Series. Imagine a complex sound, like a symphony. A Fourier Series breaks that sound down into a list of individual musical notes (frequencies).

The Twist: In this specific "hot" universe, the notes don't add up to a smooth, perfect melody. If you try to play them all at once, the sound is chaotic.
The Insight: The authors realized that while the notes don't make a smooth song, they do make sense if you treat them as a statistical pattern (a "distribution"). It's like realizing that while you can't predict exactly where one raindrop will land, you can perfectly predict the pattern of the puddle it makes.

3. The "Bouncing" Mystery

One of the biggest mysteries in black hole physics is the Singularity. This is the center of the black hole where the laws of physics break down.

The Old Theory: Scientists thought that if you tried to look at the "music" of the black hole from certain angles, you would see "ghost notes" or "echoes" caused by waves bouncing off the singularity. These are called bouncing singularities.
The Discovery: The authors found that in this specific "hot" (Euclidean) setup, those ghost notes are silent. The "bouncing" happens, but the math says the signal is zero. It's like a ball bouncing off a wall in a soundproof room; the bounce happens, but no sound comes out.
Why it matters: This confirms that the "music" of the universe is smooth and predictable in the safe zone between the start and the end of time, without any weird glitches from the singularity interfering.

4. The "Double-Counting" Trick (OPE Coefficients)

In physics, when particles interact, they often create "double" effects (like two people clapping creating a louder sound than one). This is called the Double-Trace sector.

The Old Way: Previously, scientists had to calculate the "single clap" (stress-tensor data) and then use a complicated guessing game (bootstrapping) to figure out the "double clap."
The New Way: Because the authors used the Fourier Series (the musical notes), the "double clap" information was already baked into the notes. They didn't need to guess; they just had to read the sheet music. This is a huge shortcut that saves a lot of time and effort.

5. The "Magic Glasses" (Padé Approximants)

Since the musical notes (the Fourier series) are chaotic and don't converge nicely, the authors used a mathematical trick called Padé Approximants.

The Analogy: Imagine you are looking at a blurry, pixelated image of a face. You can't see the details. But if you use a special filter (Padé approximants) that connects the dots intelligently, the image suddenly snaps into focus, revealing the face clearly.
The Result: This filter allowed them to see the "shape" of the black hole's influence in the complex mathematical space, confirming that everything is smooth and orderly in the region where it should be.

Summary: Why Should You Care?

This paper is a masterclass in changing the perspective.

Instead of fighting the messy math of black holes, they found a way to listen to the "notes" of the universe.
They discovered that the "ghosts" (singularities) are silent in this specific context, which simplifies our understanding of how the universe behaves at high temperatures.
They found a shortcut to calculate complex particle interactions that used to take days of supercomputer time.

In short, they took a chaotic, noisy storm and showed us that if you listen to the right frequency, the music is actually a perfect, harmonious song.

1. Problem Statement

The paper addresses the computation of the Euclidean thermal two-point correlation function for scalar operators in a strongly coupled Quantum Field Theory (QFT) at finite temperature $T = \beta^{-1}$ . In the context of the AdS/CFT correspondence (holography), this corresponds to studying scalar perturbations in an asymptotically AdS black hole spacetime (specifically the planar Schwarzschild-AdS $_5$ black brane).

The primary challenges addressed are:

Analytic Structure: Understanding the singularities of the correlator in the complex Euclidean time plane ( $\tau$ ), particularly those associated with "bouncing" geodesics off the black hole singularity.
Operator Product Expansion (OPE): Extracting both the stress-tensor sector and the double-trace sector ( $[\mathcal{O}\mathcal{O}]$ ) directly from the holographic calculation without relying on "thermal bootstrap" methods or image sums to recover the double-trace data from stress-tensor data.
Convergence: The standard Fourier series representation of the thermal correlator does not converge as a function but only as a distribution. The authors aim to rigorously handle this distributional nature to extract physical data.

2. Methodology

The authors adopt a strategy based on the Fourier series representation of the correlator on the thermal circle $S^1_\beta$ , rather than solving partial differential equations (PDEs) numerically for the full spacetime or bootstrapping from OPE data.

A. Fourier Series Representation

The Euclidean correlator $G_E(\tau, k)$ at fixed spatial momentum $k$ is expressed as:
$G_E(\tau, k) = \frac{1}{\beta} \sum_{n=-\infty}^{\infty} \tilde{G}_E(\zeta_n, k) e^{i\zeta_n \tau}, \quad \zeta_n = \frac{2\pi n}{\beta}$
The coefficients $\tilde{G}_E(\zeta_n, k)$ are computed by solving the radial Ordinary Differential Equation (ODE) for the scalar field in the black hole background.

B. Holographic Computation of Coefficients

The radial ODE reduces to the Heun differential equation. The Fourier coefficients are determined by the ratio of connection coefficients between local solutions at the AdS boundary and the black hole horizon. The authors employ three complementary techniques to evaluate these coefficients:

Numerical Wronskians: Direct numerical evaluation of Heun functions and their derivatives to compute connection coefficients.
Instanton Expansion: An analytic approximation using the Seiberg-Witten geometry viewpoint, where connection coefficients are expanded in a parameter $X$ (related to the instanton counting parameter).
Asymptotic $1/n$ Expansion: An analytic expansion of the coefficients for large Matsubara frequency $n$ , derived using the method of [6].

C. Resurgence and Distributional Analysis

Resurgence Analysis: The authors perform a Borel resummation of the asymptotic $1/n$ expansion. They analyze the transseries structure to determine if non-perturbative sectors (associated with bouncing singularities) are active.
Regularization: Since the Fourier series diverges as a function (coefficients grow as $n^{2\Delta-4}$ ), the authors utilize Padé approximants and $i\epsilon$ prescriptions (Abel summation) to analytically continue the series into the complex $\tau$ plane and extract finite local data.

3. Key Contributions and Results

A. Direct Extraction of Double-Trace OPE Coefficients

A major result is the ability to extract all OPE coefficients, including the double-trace sector, directly from the asymptotic expansion of the Fourier coefficients.

The stress-tensor sector coefficients ( $a_m$ ) are determined by the leading asymptotic behavior.
The double-trace coefficients ( $b_q$ ) are determined by subleading terms in the same expansion.
Significance: This bypasses the need for previous "thermal bootstrap" techniques where double-trace data had to be inferred from stress-tensor data via periodicity constraints. The Fourier series representation is manifestly periodic, naturally containing both sectors.

B. Analytic Structure and Bouncing Singularities

Analyticity: The authors confirm that the correlator is analytic in the strip $0 < \text{Re}(\tau) < \beta$ .
Bouncing Singularities: In the complex $\tau$ plane, singularities associated with geodesics bouncing off the black hole singularity are located at $\tau_j = \frac{\beta}{\sqrt{2}} e^{i\pi/4(1+2j)}$ .
Transseries Analysis: Using resurgence theory, the authors show that while the transseries for the Fourier coefficients contains terms corresponding to these bouncing singularities (exponentially small terms $e^{-w_j n}$ $e^{- w_{j} n}$ ), the transseries parameters ( $\sigma_j$ ) are exactly zero for the Euclidean correlator.
- Contrast: This differs from the retarded Green's function, where these non-perturbative terms are active.
- Implication: The Borel resummation of the perturbative $1/n$ series is exact for the Euclidean case; no non-perturbative corrections are needed.

C. Distributional Nature of the Correlator

The paper rigorously treats the correlator as a distribution rather than a function.

The Fourier series converges in the sense of distributions (tested against smooth functions).
Using the instanton approximation for integer conformal dimensions $\Delta$ , the authors explicitly derive that the correlator consists of contact terms (derivatives of delta functions) at the coincident point $\tau = 0$ (and its images), consistent with QFT expectations.
Away from singularities, the regularized series reproduces the expected smooth behavior.

D. Numerical Validation

The authors validate their results by comparing:

OPE coefficients derived from the Borel-resummed asymptotic expansion.
OPE coefficients derived from Padé approximants of the numerically computed Fourier series.
The results show excellent numerical agreement, confirming the consistency of the three computational approaches.

4. Significance

Methodological Advancement: The work demonstrates that working directly with the Fourier series on the thermal circle is a powerful and efficient alternative to solving PDEs or using bootstrap constraints. It provides a unified framework for computing both stress-tensor and double-trace OPE data.
Insight into Holographic Singularities: The finding that non-perturbative "bouncing" sectors are turned off in the Euclidean thermal correlator (but would be active in the Lorentzian/retarded case) provides a precise understanding of how holographic singularities manifest in different regimes.
Distributional Framework: The paper reinforces the necessity of treating thermal correlators as distributions, particularly when dealing with coincident points and analytic continuation, aligning holographic calculations with rigorous QFT axioms.
Practical Utility: The authors provide explicit formulas (Eqs. 5.6, 5.7) for OPE coefficients and release a Mathematica notebook with 70 coefficients of the asymptotic expansion, facilitating further research in thermal holography.

In summary, this paper provides a complete, non-perturbative (in the sense of exact Borel resummation) and distributionally consistent calculation of holographic thermal correlators, successfully extracting the full OPE structure and clarifying the role of black hole singularities in the Euclidean regime.