[Paper Review] Higher Dimensional Coulomb Gases and Renormalized Energy Functionals

Nicolas Rougerie, Sylvia Serfaty|arXiv (Cornell University)|Jul 10, 2013

Random Matrices and Applications89 references17 citations

TL;DR

This paper introduces a renormalized energy functional to characterize next-to-leading order corrections in the ground state energy of higher-dimensional Coulomb gases under mean-field scaling. By splitting the Hamiltonian and using a smeared jellium model, it rigorously derives that fluctuations are governed by a universal renormalized energy, extending Sandier and Serfaty's 2D results to arbitrary dimensions $ d \geq 2 $, with applications to free energy asymptotics, Gibbs measures, and charge fluctuations.

ABSTRACT

We consider a classical system of n charged particles in an external confining potential, in any dimension d larger than 2. The particles interact via pairwise repulsive Coulomb forces and the coupling parameter scales like the inverse of n (mean-field scaling). By a suitable splitting of the Hamiltonian, we extract the next to leading order term in the ground state energy, beyond the mean-field limit. We show that this next order term, which characterizes the fluctuations of the system, is governed by a new "renormalized energy" functional providing a way to compute the total Coulomb energy of a jellium (i.e. an infinite set of point charges screened by a uniform neutralizing background), in any dimension. The renormalization that cuts out the infinite part of the energy is achieved by smearing out the point charges at a small scale, as in Onsager's lemma. We obtain consequences for the statistical mechanics of the Coulomb gas: next to leading order asymptotic expansion of the free energy or partition function, characterizations of the Gibbs measures, estimates on the local charge fluctuations and factorization estimates for reduced densities. This extends results of Sandier and Serfaty to dimension higher than two by an alternative approach.

Motivation & Objective

• To characterize the next-to-leading order correction in the ground state energy of classical Coulomb gases in dimensions $ d \geq 2 $, beyond the mean-field approximation.
• To define and rigorously analyze a new renormalized energy functional that captures the effective interaction energy of a jellium system with point charges screened by a uniform background.
• To extend previous 2D results of Sandier and Serfaty to higher dimensions using a novel splitting and smearing approach.
• To derive consequences for statistical mechanics, including asymptotic expansions of the partition function, properties of Gibbs measures, and local charge fluctuations.

Proposed method

• A splitting of the total Hamiltonian into mean-field and fluctuation terms, isolating the next-order energy correction.
• Introduction of a renormalized energy functional via smearing of point charges at a small scale, inspired by Onsager’s lemma, to remove divergences.
• Use of periodic boundary conditions and Green’s formula to compute the renormalized energy in terms of the gradient of a solution to a screened Poisson equation.
• Application of Fourier series to express the Green’s function on the torus, leading to Eisenstein series and connections to Epstein zeta functions in 2D.
• Establishment of equivalence between the renormalized energy $ \mathcal{W} $ and the energy $ W $ of a periodic configuration, enabling minimization over lattices.
• Use of the first Kronecker limit formula in 2D to relate the renormalized energy to the zeta function of the dual lattice, enabling minimization over lattices.

Experimental results

Research questions

• RQ1 How can the next-to-leading order correction in the ground state energy of a $ d $-dimensional Coulomb gas be characterized beyond the mean-field limit?
• RQ2 What is the correct renormalized energy functional that governs the effective energy of a jellium system in arbitrary dimension $ d \geq 2 $?
• RQ3 How does the renormalized energy relate to known mathematical objects such as Eisenstein series and Epstein zeta functions in 2D?
• RQ4 Can the statistical mechanics of the Coulomb gas—such as the partition function, Gibbs measures, and charge fluctuations—be rigorously analyzed using this new functional?
• RQ5 What is the role of lattice geometry in minimizing the renormalized energy, particularly in dimensions $ d = 2, 8, 24 $?

Key findings

• The next-to-leading order correction in the ground state energy is governed by a universal renormalized energy functional $ \mathcal{W} $, defined via the gradient of a solution to a screened Poisson equation.
• The renormalized energy $ \mathcal{W} $ is finite and well-defined even for singular point charges, achieved by smearing them at a small scale, consistent with Onsager’s regularization.
• In 2D, the renormalized energy for a lattice is given by $ \mathcal{W}(\Lambda) = c_d^2 \lim_{x\to 0} \left( E_\Lambda(x) - \frac{w(x)}{c_d} \right) $, where $ E_\Lambda $ is an Eisenstein series.
• The minimization of $ \mathcal{W} $ over lattices of unit volume in 2D corresponds to minimizing the Epstein zeta function, with the triangular lattice as the unique minimizer.
• The asymptotic expansion of the partition function is derived, showing that the free energy correction is governed by $ \mathcal{W} $, with precise control on fluctuations.
• Local charge fluctuations and reduced density factorization estimates are obtained, showing that the system exhibits universal behavior described by $ \mathcal{W} $.