Linear statistics in the Coulomb gas: cumulants beyond rotational invariance

Résumé technique : Statistiques linéaires dans le gaz de Coulomb : les cumulants au-delà de l'invariance par rotation

Énoncé du problème
L'article étudie les fluctuations des statistiques linéaires d'un gaz de Coulomb classique de $N$ particules en dimension $d \ge 2$, confiné par un potentiel externe général $U$. Les particules interagissent via une répulsion logarithmique ($d=2$) ou newtonienne ($d=3$). L'objet principal d'étude est la statistique linéaire $L_N = \sum_{i=1}^N f(x_i)$, où $f$ est une fonction test lisse.

Bien que la variance d'ordre un des statistiques linéaires soit bien comprise et régie par un théorème de la limite centrale avec des fluctuations gaussiennes (de l'ordre de $O(1)$), le comportement des cumulants d'ordre supérieur pour des potentiels généraux non symétriques et des fonctions tests demeure largement inexploré. Les résultats précédents pour les cumulants supérieurs (par exemple dans [29]) étaient restreints aux contextes d'invariance par rotation où la goutte d'équilibre (le support de la densité de particules) reste une boule sous perturbation. Ce travail vise à dériver des formules explicites pour les troisième et quatrième cumulants principaux pour des gouttelettes lisses générales $D_0$, sans supposer l'invariance par rotation pour le potentiel $U$ ou la fonction test $f$.

Méthodologie
Les auteurs emploient une combinaison de calcul de forme, de théorie du potentiel et de techniques de mise en correspondance conforme :

1. Potentiel incliné et fonction génératrice des cumulants : Les cumulants sont dérivés de la fonction génératrice des cumulants (CGF) $\chi(s, N) = \log \langle e^{-Ns L_N} \rangle_U$. En introduisant un potentiel incliné $U_s = U + s f / \beta$, le problème est réduit à l'analyse de la fonction de partition $Z_{N, U_s}$. Le $q$-ième cumulant est obtenu en dérivant $\chi(s, N)$ $q$ fois par rapport à $s$ en $s=0$.
2. Problème de l'obstacle et vitesse de la frontière : La dérivation centrale repose sur la connexion entre la mesure d'équilibre et le problème de l'obstacle. À mesure que le paramètre $s$ varie, le support de la densité d'équilibre (la gouttelette $D_s$) se déforme. Les auteurs utilisent le résultat connu selon lequel la vitesse normale $v_s$ de la frontière de la gouttelette $\partial D_s$ est déterminée par le saut de la dérivée normale de l'extension harmonique de la fonction test. Spécifiquement, $v_s = [\vec{\nabla} f_{D_s}] \cdot n_s / (\beta \Delta U_s)$.
3. Dérivées de forme : Les auteurs calculent les dérivées de l'énergie de Dirichlet associée à l'extension harmonique de $f$ par rapport à la déformation du domaine. Cela implique le calcul de variations d'intégrales de surface sur la frontière mobile $\partial D_s$, en utilisant la formule de Hadamard et les propriétés de l'opérateur de Dirichlet-vers-Neumann extérieur ($\Lambda_+$).
4. Mise en correspondance conforme ($d=2$) : Pour le cas bidimensionnel, une dérivation indépendante est fournie en utilisant l'analyse complexe. Cette approche utilise la carte de Riemann extérieure, la formule de variation de Julia pour l'évolution de la carte conforme sous le mouvement de la frontière, et la transformée de Hilbert pour relier la vitesse de la frontière à la fonction test.
5. Justification rigoureuse ($d=2, \beta=2$) : La section 4 fournit une preuve rigoureuse des formules asymptotiques pour les troisième et quatrième cumulants dans le cas déterministe ($d=2, \beta=2$) avec des données réal-analytiques. Cela repose sur les asymptotiques des polynômes orthogonaux et l'expansion de Bergman dans le volume pour justifier l'interversion des limites ($N \to \infty$ et $s \to 0$).

Contributions clés et résultats

1. Formule générale du troisième cumulant :
   L'article dérive une formule explicite pour le troisième cumulant principal $\langle L_N^3 \rangle_c$ pour toute dimension $d \ge 2$ et toute gouttelette lisse $D_0$. Le résultat est exprimé comme une intégrale de bord intrinsèque sur la gouttelette d'équilibre :
   $$\langle L_N^3 \rangle_c \simeq -\frac{1}{N \beta^2 \Omega_d} \int_{\partial D_0} dS \frac{\left( [\vec{\nabla} f_{D_0}] \cdot n \right)^3}{\Delta U}$$
   où $[\vec{\nabla} f_{D_0}] \cdot n$ est le saut de la dérivée normale de l'extension harmonique de $f$ à travers la frontière, $n$ est la normale unitaire sortante, et $\Delta U$ est le laplacien du potentiel externe (proportionnel à la densité de fond). La formule dépend explicitement du cube de ce saut.

2. Formule générale du quatrième cumulant :
   Le quatrième cumulant est dérivé comme une intégrale de bord impliquant des opérateurs géométriques et spectraux additionnels :
   $$\langle L_N^4 \rangle_c \simeq -\frac{1}{\beta^3 \Omega_d N^2} \int_{\partial D_0} dS \left[ 3 \frac{K^2}{\Delta U} \Lambda_+ \left( \frac{K^2}{\Delta U} \right) + 4 \frac{K^3 \Delta f}{(\Delta U)^2} + 2M \frac{K^4}{(\Delta U)^2} + \frac{K^4 \partial_n \Delta U}{(\Delta U)^3} \right]$$
   où $K = [\vec{\nabla} f_{D_0}] \cdot n$, $\Lambda_+$ est l'opérateur de Dirichlet-vers-Neumann extérieur, $M$ est la courbure moyenne de la frontière, et $\partial_n \Delta U$ est la dérivée normale de la densité de fond. Contrairement au troisième cumulant, le quatrième implique l'action de l'opérateur de Dirichlet-vers-Neumann sur le carré du saut, ainsi que des termes de courbure et de gradient de densité.

3. Géométries spéciales :
   Les formules générales sont spécialisées pour des géométries spécifiques, incluant :

   • Disque ($d=2$) : Exprimé en termes de transformée de Hilbert et de modes de Fourier.
   • Demi-plan ($d=2$) : Simplifié pour impliquer la transformée de Hilbert de la dérivée de $f$.
   • Sphère ($d=3$) : Exprimé en utilisant les harmoniques sphériques et les valeurs propres de l'opérateur de Dirichlet-vers-Neumann extérieur.
     Ces résultats spécialisés récupèrent les formules radiales connues lorsque $f$ est invariante par rotation.

4. Validation :
   Les auteurs vérifient leurs formules générales par rapport à des déformations non radiales exactement solubles :

   • Gouttelettes elliptiques traduites en $d=2$ (Appendice A).
   • Déformations sphéroïdales d'une boule en $d=3$ (Appendice B).
     Dans le cas $d=2, \beta=2$, il est montré que les résultats correspondent aux formules proposées par Wiegmann et Zabrodin [51], fournissant une preuve rigoureuse pour les troisième et quatrième cumulants dans ce cadre spécifique.

Signification et revendications
L'article affirme fournir les premières formules explicites pour les troisième et quatrième cumulants principaux des statistiques linéaires lisses dans les gaz de Coulomb pour des potentiels généraux non invariants par rotation et des fonctions tests en dimensions arbitraires $d \ge 2$.

• Au-delà des fluctuations gaussiennes : Bien que la variance soit bien connue pour être de l'ordre de $O(1)$ et gaussienne, ces cumulants d'ordre supérieur représentent les premières corrections non gaussiennes. Les formules dérivées montrent explicitement les facteurs d'échelle $N^{-1}$ et $N^{-2}$ dans les intégrandes, cohérents avec le comportement observé dans les cas invariants par rotation, mais la nouveauté principale réside dans la structure intégrale complexe impliquant les sauts de frontière, les opérateurs de Dirichlet-vers-Neumann et la courbure.
• Dépendance de frontière intrinsèque : Une conclusion clé est que ces cumulants d'ordre supérieur dépendent uniquement des valeurs de bord de la fonction test et de ses dérivées, ainsi que de la géométrie de la frontière de la gouttelette (courbure, dérivées normales de la densité). Cela étend l'universalité de bord observée pour la variance aux ordres supérieurs.
• Fondement rigoureux : Pour le cas $d=2, \beta=2$, l'article propose une dérivation rigoureuse qui valide les résultats heuristiques précédents de la littérature de physique (Wiegmann et Zabrodin).
• Limitations : Les auteurs notent que la dérivation générale suppose que les gouttelettes perturbées restent lisses et ne subissent pas de transitions topologiques (par exemple, division ou fusion). La preuve rigoureuse pour un $d$ général et un $\beta$ général repose sur une condition suffisante concernant les dérivées de la fonction de partition (Remarque 4.6), qui n'est pas entièrement prouvée pour tous les cas, mais qui est établie pour le cas déterministe.

Ce travail fait le pont entre la mécanique statistique, la théorie des matrices aléatoires et les problèmes de frontière libre (problèmes de l'obstacle), offrant un cadre unifié pour comprendre les fluctuations non gaussiennes dans les systèmes de particules en interaction au-delà des configurations symétriques.