Coulomb gas and the Grunsky operator on a Jordan domain with corners

This paper establishes an exact expression for the planar Coulomb gas partition function on a Jordan domain in terms of the Grunsky operator and derives its large-$n$ asymptotic behavior, proving that the geometry of the boundary (specifically whether it is a Weil-Petersson quasicircle or has corners) determines the leading order corrections to the partition function via the Loewner energy and corner opening angles.

The Gist

Imagine you have a magical, elastic sheet of rubber (a 2D surface) and you sprinkle thousands of tiny, identical magnets onto it. These magnets repel each other fiercely—they hate being close together. However, you've placed them inside a specific shape drawn on the table, like a circle, a square, or a weirdly shaped blob with sharp corners. The magnets are also trapped by a "hard wall" (the boundary of the shape) that they cannot cross.

This setup is what mathematicians call a Coulomb gas. The paper by Kurt Johansson and Fredrik Viklund asks a very specific question: How does the shape of the container affect the total "energy" or "disorder" of these magnets as you keep adding more and more of them?

Here is the breakdown of their discovery using simple analogies:

1. The Perfect Circle vs. The Bumpy Shape

If your container is a perfect circle, the magnets arrange themselves in a very predictable, smooth way. The math describing this is clean and elegant.

But what if your container has corners? Imagine a star shape or a polygon. At the sharp points (the corners), the magnets get confused. They can't spread out smoothly; they get "cramped" or "stretched" depending on how sharp the angle is.

• The Analogy: Think of a crowd of people trying to exit a room. If the door is a smooth arch, everyone flows out evenly. If the door is a sharp, jagged V-shape, people get stuck or pile up differently at the tips. The paper calculates exactly how much "extra trouble" (mathematical energy) these sharp corners cause.

2. The "Grunsky Operator": The Shape's Fingerprint

To solve this, the authors use a complex mathematical tool called the Grunsky operator.

• The Analogy: Imagine you want to describe the shape of a cookie without showing the cookie itself. You could take a photo of its shadow, or measure how light bends around it. The Grunsky operator is like a "mathematical shadow" or a fingerprint of the shape. It translates the physical geometry of the boundary (the wall) into a giant grid of numbers (a matrix).
• The paper proves that if you look at this grid of numbers, you can tell if the shape is "smooth" (mathematically perfect) or if it has "rough edges" (corners).

3. The Main Discovery: The "Corner Penalty"

The authors found a precise formula for what happens when the shape has corners.

• The Smooth Case: If the shape is perfectly smooth (no corners), the energy grows in a standard, predictable way.
• The Corner Case: If the shape has corners, there is an extra cost to the energy. This cost doesn't just add a little bit; it grows logarithmically (like the slow, steady ticking of a clock) as you add more magnets.

The formula they derived looks like this:
$$\text{Extra Energy} \approx \text{Number of Corners} \times \text{Function of the Angle}$$

The "Function of the Angle" is the star of the show. It tells you that:

• A very sharp corner (like a needle point) causes a huge amount of "friction" or energy cost.
• A wide, blunt corner causes less cost.
• Interestingly, the formula is universal. It doesn't matter if the shape is a star, a house, or a random polygon; the math only cares about the angle of the corner. It's like saying, "Every time you turn a corner at a 30-degree angle, you pay this specific tax."

4. Why Does This Matter? (The "Why Should I Care?")

You might wonder, "Who cares about magnets in a box?"

• Physics: This helps explain how electrons behave in tiny quantum devices or how materials change state (like ice melting).
• Random Matrices: In computer science and statistics, we often deal with huge tables of random numbers. The patterns in these tables often look like these magnet arrangements. Understanding the "corners" helps us predict how these systems behave when they get huge.
• Universality: The most exciting part is that the math they found appears in completely different fields, like the study of heat spreading through a metal plate or the behavior of light. It suggests that nature has a "standard fee" for sharp corners, regardless of what the system actually is.

Summary

The paper is a detective story where the authors:

1. Set the scene: Magnets repelling each other inside a container.
2. Found the clue: Used a "fingerprint" tool (Grunsky operator) to analyze the container's shape.
3. Solved the mystery: Proved that sharp corners create a specific, predictable "energy tax" that grows slowly as the system gets bigger.
4. The takeaway: Geometry isn't just about looks; sharp corners leave a permanent, calculable mark on the physics of the system, and that mark follows a beautiful, universal rule.

In short: If you have a shape with corners, the universe charges you a specific "corner fee" in energy, and the authors finally wrote down the exact receipt.

Technical Summary

1. Problem Statement and Motivation

The paper investigates the asymptotic behavior of the partition function $Z_n(D)$ for a planar Coulomb gas confined to a bounded Jordan domain $D$ with a hard wall boundary $\eta = \partial D$. The system consists of $n$ charged particles interacting via a logarithmic potential at inverse temperature $\beta=2$ (the determinantal case).

The central question is how the geometry of the boundary $\eta$, specifically the presence of corner singularities, influences the large $n$ expansion of the free energy $\log Z_n(D)$.

• Context: For smooth boundaries, the free energy is closely related to the Loewner energy ($I_L$) and the Weil-Petersson class of quasicircles. If the boundary is smooth (or a Weil-Petersson quasicircle), the normalized free energy converges to a finite limit proportional to $I_L$.
• The Challenge: If the boundary has corners (non-smooth points), the Loewner energy diverges ($I_L = \infty$). The authors aim to quantify this divergence and determine the precise asymptotic correction terms introduced by the corners.
• Specific Goal: To derive an asymptotic formula for $\log Z_n(D)$ when $\eta$ is a piecewise analytic curve with $m$ corners having interior opening angles $\alpha_p \pi$.

2. Methodology

The authors employ a sophisticated blend of complex analysis, spectral theory, and asymptotic analysis.

A. The Grunsky Operator Connection

The core of the methodology relies on an exact identity relating the partition function to the Grunsky operator $B$.

• Let $g: D^* \to D^*$ be the exterior conformal map normalized at infinity. The Grunsky coefficients $a_{k\ell}$ are derived from the expansion of $\log \frac{g(\zeta)-g(z)}{\zeta-z}$.
• The normalized partition function $\bar{Z}_n(D)$ satisfies:
  $$\log \bar{Z}_n(D) = \log \det(I - P_n B B^* P_n)$$
  where $P_n$ is the projection onto the first $n$ coordinates of $\ell^2(\mathbb{Z}^+)$.
• Thus, the problem reduces to analyzing the asymptotic behavior of the Fredholm determinant of the truncated Grunsky operator as $n \to \infty$.

B. Asymptotics of Grunsky Coefficients

The main technical breakthrough is the derivation of precise asymptotics for the Grunsky coefficients $b_{k\ell} = \sqrt{k\ell}a_{k\ell}$ when the domain has corners.

• Contour Integration: The authors use a double contour integral representation for the coefficients. They deform the integration contours to pass near the pre-images of the corners on the unit circle.
• Local Analysis: Near a corner with exterior angle $\gamma_p \pi$ (where $\gamma_p = 2 - \alpha_p$), the conformal map behaves like $z \mapsto z^{\alpha_p}$. This singularity dictates the leading order behavior of the coefficients.
• Main Asymptotic Result (Theorem 1.7): The authors prove that for large $k, \ell$:
  $$b_{k\ell} \approx K(k, \ell) + \text{error}$$
  where $K(k, \ell)$ is a sum of "Hardy kernels" $K_p(k, \ell)$ associated with each corner. The error term decays sufficiently fast ($O((k+\ell)^{-1-\rho})$) to be negligible in the leading logarithmic asymptotics.

C. Trace Expansion and Integral Approximation

To evaluate the determinant $\log \det(I - P_n B B^* P_n)$, the authors expand the logarithm into a series of traces:
$$\log \det(I - C_n) = -\sum_{i=1}^\infty \frac{1}{i} \text{tr}(C_n^i)$$

• They replace the operator $B$ with the kernel operator $K$ derived from the asymptotics.
• They show that the traces of powers of $K$ can be approximated by integrals involving the Fourier transform of the kernel functions $H_p(t)$.
• A key step involves evaluating the integral:
  $$\int_{\mathbb{R}} \log\left(1 - \frac{\sinh^2(\beta x)}{\sinh^2 x}\right) dx$$
  which yields the specific geometric factor associated with the corner angles.

3. Key Contributions and Results

A. Main Theorem: Asymptotics for Domains with Corners

Theorem 1.3 provides the central result. For a domain $D$ with $m$ corners of interior angles $\alpha_p \pi$:
$$\lim_{n \to \infty} \frac{1}{\log n} \log \frac{\bar{Z}_n(D)}{\bar{Z}_n(\mathbb{D})} = -\frac{1}{6} \sum_{p=1}^m \left( \frac{\alpha_p + 1}{\alpha_p} - 2 \right)$$

• Significance: This confirms physics predictions regarding the universality of logarithmic divergences in free energy expansions for systems with conical singularities. The term $\left( \frac{\alpha_p + 1}{\alpha_p} - 2 \right)$ is the "corner anomaly."
• Interpretation: The divergence is of order $\log n$, distinct from the $O(n^2)$ and $O(n)$ terms found in smooth domains.

B. Loewner Energy and Equipotentials

Theorem 1.4 analyzes the Loewner energy $I_L(\eta_r)$ of a sequence of analytic curves $\eta_r$ (equipotentials) approximating the cornered boundary $\eta$ from the exterior.
$$\lim_{r \to 1^+} \frac{I_L(\eta_r)}{\log(\frac{1}{r-1})} = 2 \sum_{p=1}^m \left( \frac{\alpha_p + 1}{\alpha_p} - 2 \right)$$
This establishes a direct link between the divergence rate of the Loewner energy and the corner geometry, showing that the "reduced" Loewner energy is finite only for smooth curves.

C. Fekete-Pommerenke Energy

Theorem 1.6 extends the analysis to the Fekete-Pommerenke energy $I_F$ (related to the maximal discriminant/Fekete points). For exterior equipotentials $\eta_r$:
$$\lim_{r \to 1^+} \frac{I_F(\eta_r)}{\log(\frac{1}{r-1})} = 2 \sum_{p=1}^m \left( \frac{\gamma_p + 1}{\gamma_p} - 2 \right)$$
where $\gamma_p = 2 - \alpha_p$ are the exterior angles. This highlights a duality between interior and exterior corner contributions.

D. Characterization of Weil-Petersson Quasicircles

The paper reconfirms that a Jordan curve is a Weil-Petersson quasicircle if and only if the Loewner energy is finite. The authors show that if the boundary has corners, the free energy diverges logarithmically, providing a quantitative measure of how far the curve is from being Weil-Petersson.

4. Significance and Impact

1. Rigorous Validation of Physics Predictions: The paper provides a rigorous mathematical derivation of the "corner anomaly" term in the free energy expansion, which had previously been predicted by physicists studying critical systems and conformal field theory (CFT).
2. Bridging Geometric Function Theory and Statistical Mechanics: It establishes a deep connection between the spectral properties of the Grunsky operator (a classical object in complex analysis) and the thermodynamic properties of Coulomb gases.
3. Universality: The specific function of the corner angles appearing in the limit is shown to be universal, appearing in diverse contexts such as the heat trace expansion (Mellin transform of the heat kernel) and the Loewner energy.
4. New Analytical Tools: The detailed asymptotic analysis of Grunsky coefficients for piecewise analytic domains (Theorem 1.7) is a significant technical achievement. It provides a template for studying other spectral problems on domains with singularities.
5. Open Problems: The authors formulate conjectures (Conjectures 1.8, 1.10, 1.11, 1.12) regarding the existence of "reduced" energies and the behavior of partition functions for general $\beta$ and less regular curves, setting a roadmap for future research.

In summary, this paper successfully quantifies the impact of geometric singularities (corners) on the statistical mechanics of planar Coulomb gases, using the Grunsky operator as the primary analytical tool to bridge complex geometry and free energy asymptotics.