Free energy of the Coulomb gas in the determinantal case on Riemann surfaces

This paper derives the asymptotic expansion of the partition function for a Coulomb gas system on compact Riemann surfaces of any genus by employing a bosonization formula to relate analytic torsion and geometric quantities, thereby proving the geometric version of the Zabrodin-Wiegmann conjecture in the determinantal case.

The Gist

Imagine a crowded dance floor on a curved surface, like the surface of a sphere, a donut, or a pretzel with many holes. This is the setting for the "Coulomb gas" described in Lucas Bourgoin's paper.

Here is the story of what the paper does, broken down into simple concepts:

1. The Dance Floor and the Dancers

Imagine $N$ tiny, charged dancers (particles) on a closed, curved stage (a Riemann surface).

• The Interaction: These dancers repel each other. They want to get as far apart as possible, but they are stuck on the stage. This repulsion is like the "Coulomb force" (think of how two magnets with the same pole push away from each other).
• The Goal: The paper asks a very specific question: If we have a huge number of dancers (approaching infinity), what is the total "energy cost" or "free energy" of this chaotic dance?

In physics, this "free energy" is calculated using something called a Partition Function (let's call it $Z$). It's a giant mathematical recipe that sums up every possible way the dancers could arrange themselves.

2. The "Determinantal" Case: A Perfectly Organized Chaos

The paper focuses on a special scenario called the "determinantal case."

• The Analogy: Usually, if you have a crowd of people, they move randomly. But in this specific case, the dancers are like a perfectly choreographed troupe. Their movements are linked in a way that prevents them from ever bumping into each other.
• The Math: This "perfect organization" allows mathematicians to use a special tool called a determinant (a specific type of calculation used in linear algebra) to describe the system. It turns a messy, chaotic problem into a structured one that can be solved.

3. The Map and the Compass (Metrics and Green Functions)

To calculate the energy, the author needs a way to measure distances and forces on these curved surfaces.

• The Green Function: Think of this as a "force map." It tells you how strongly one dancer pushes another based on their distance.
• The Metrics: The paper uses two specific "rulers" to measure the surface:
  1. The Canonical Metric: A standard, natural way to measure the shape of the surface.
  2. The Arakelov Metric: A more complex, specialized ruler used in advanced geometry.
• The Trick: The author switches between these rulers to make the math easier, much like a cartographer switching between a flat map and a globe to measure a route.

4. The Magic Spell: Bosonization

This is the paper's main "magic trick."

• The Problem: Calculating the energy of $N$ interacting particles is incredibly hard.
• The Solution: The author uses a formula called the Bosonization Formula.
• The Analogy: Imagine trying to count the noise of a thousand people shouting. Instead of listening to every voice, the Bosonization formula is like a translator that converts the "shouting" (the particles) into a "symphony" (a single, elegant wave of sound).
• What it connects: It links the messy world of the dancing particles to the clean, quiet world of Analytic Torsion (a way of measuring the "vibration" or "shape" of the surface itself). It essentially says: "The energy of the crowd is directly related to the shape of the stage."

5. The Big Discovery: The Final Formula

After doing a massive amount of complex math, the author derives a final formula that predicts the energy as the number of dancers ($N$) gets huge.

The formula looks like this:
$$\text{Energy} \approx (\text{Big Number}) \times N^2 + (\text{Medium Number}) \times N \ln(N) + \dots + (\text{The Secret Constant})$$

• The Big Terms: The first few terms ($N^2$, $N \ln N$) describe the obvious, bulk behavior of the crowd.
• The Secret Constant ($b_0$): This is the most important part of the paper. The author proves that the final, constant term in the formula contains the logarithm of the determinant of the Laplacian.
  • What is the Laplacian? Think of it as a machine that measures how "curvy" or "wiggly" the surface is. Its "determinant" is a single number that summarizes the entire geometry of the stage.
  • Why it matters: The paper confirms a famous guess (the Zabrodin-Wiegmann conjecture). It proves that the "shape" of the universe (the Riemann surface) leaves a permanent fingerprint on the energy of the particles, even when there are infinite of them.

6. The "Fluctuations" (The Wiggles)

The paper also looks at what happens if the dancers don't follow the perfect choreography exactly.

• The Analogy: If the perfect dance is a straight line, the "fluctuations" are the tiny, random wiggles the dancers make around that line.
• The Result: The author proves that these wiggles follow a Normal Distribution (the famous "Bell Curve"). This means that while the dancers move randomly, their average behavior is predictable and follows a standard statistical pattern.

Summary

In simple terms, Lucas Bourgoin solved a puzzle about how a massive crowd of repelling particles behaves on a curved, multi-holed surface. By using a mathematical "translator" (Bosonization) to turn the crowd's behavior into a question about the shape of the surface itself, he proved that the surface's geometry is written into the final energy calculation. This confirms a long-standing prediction about how geometry and physics are deeply intertwined in these systems.

Technical Summary

Technical Summary: Free energy of the Coulomb gas in the determinantal case on Riemann surfaces

Problem Statement
The paper investigates the asymptotic behavior of the partition function $Z$ for a system of $N$ charged particles (a Coulomb gas) interacting via a Coulomb potential on a compact Riemann surface $M$ of genus $g$. Specifically, the study focuses on the "determinantal case" where the inverse temperature parameter $\beta = 1$. The primary objective is to derive an explicit asymptotic expansion of the free energy $\ln Z$ as $N \to \infty$.

The partition function is defined as:
$$Z_\beta = \frac{1}{N!} \int_{M^N} \mu^N \exp\left( 2\beta \sum_{1\le i<j\le N} G(z_i, z_j) + \beta(N+g-1) \sum_{i=1}^N V(z_i) \right)$$
where $G$ is the Green function of the scalar Laplacian, $V$ is a quasi-subharmonic potential, and $\mu$ is a volume form. The authors aim to resolve the "geometric Zabrodin-Wiegmann conjecture," which predicts the structure of this expansion, particularly the constant term $b_0$, which is expected to involve the regularized determinant of the scalar Laplacian.

Methodology
The derivation relies on a combination of complex geometry, spectral theory, and statistical mechanics techniques:

1. Bosonization Formula: The core tool is the bosonization formula, which relates the determinant of the magnetic Laplacian (Kodaira Laplacian) on a Hermitian line bundle $L$ of degree $k = N + g - 1$ to geometric quantities, including the Arakelov-Green function and the determinant of the scalar Laplacian.
2. Metric Selection: The computation is performed using specific metrics to simplify the application of the bosonization formula:
   • The canonical metric $\rho_{can}$ is used to define the Green function.
   • The Arakelov metric $\rho_{Ar}$ is used for the integration measure and the definition of the line bundle's admissible metric.
3. Asymptotic Analysis of Determinants: The authors utilize recent results by Finski regarding the asymptotic expansion of the determinant of the magnetic Laplacian $\det \square_{L, \rho, h}$ as the degree $k \to \infty$.
4. Modified Partition Function: To handle the theta function appearing in the bosonization formula, the authors first compute the asymptotic expansion of a "modified partition function" $Z^{(\theta)}$ which includes the theta function factor. They then recover the standard partition function $Z$ by integrating out the theta function over the Jacobian torus.
5. Variational Calculus: The transformation of functionals (Liouville, Mabuchi, Aubin-Yau) and determinants under changes of metric and potential is used to relate the system with potential $V$ to the system with $V=0$.

Key Contributions and Results

• Asymptotic Expansion: The paper establishes the following asymptotic expansion for the logarithm of the partition function as $N \to \infty$:
  $$\ln Z = b_2 N^2 + a_1 N \ln N + b_1 N + a_0 \ln N + b_0 + O(N^{-1} \ln N)$$
  The coefficients $a_i$ and $b_i$ are explicitly computed in terms of the genus $g$, the potential $V$, and geometric functionals of the Riemann surface.

• Resolution of the Geometric Zabrodin-Wiegmann Conjecture: The main result confirms the conjecture in the determinantal case. Specifically, the constant term $b_0$ is shown to contain the logarithm of the regularized zeta determinant of the scalar Laplacian, $\ln \det \Delta_{\rho_{Ar}}$. This term appears alongside other geometric functionals such as the Polyakov functional, the Liouville functional, and the Mabuchi functional.

• Explicit Coefficients:

  • The $N^2$ term ($b_2$) depends on the equilibrium measure defined by the potential $V$.
  • The $N \ln N$ term ($a_1$) is found to be $-1/2$.
  • The constant term ($b_0$) is expressed using the Faltings' delta invariant and the determinant of the Laplacian, confirming the presence of the Gaussian free field structure in the constant term.

• Fluctuations: As a corollary, the paper proves that the fluctuations of linear statistics, defined as $\text{Fluct}_N f = \sum f(x_i) - N \int f d\mu_V$, converge in distribution to a normal random variable. The variance is given by a Dirichlet energy term, and the mean involves the scalar curvature and the potential.

• Exact Computations for Low Genus: The paper provides explicit checks for genus $g=0$ (sphere) and $g=1$ (torus), where the results align with direct computations and known formulas involving the Dedekind eta function and theta functions.

Significance
The paper provides a rigorous derivation of the free energy expansion for Coulomb gases on compact Riemann surfaces in the determinantal case. By explicitly identifying the scalar Laplacian determinant within the constant term of the expansion, the work confirms the geometric version of the Zabrodin-Wiegmann conjecture. This establishes a concrete link between the statistical mechanics of Coulomb gases, the geometry of Riemann surfaces (specifically Arakelov geometry), and the spectral theory of Laplacians. The results are relevant to random matrix theory (specifically the Ginibre ensemble on curved surfaces) and the study of Quantum Hall states on compact manifolds. The methodology demonstrates the utility of the bosonization formula in connecting analytic torsion to physical partition functions.