Linear statistics in the Coulomb gas: cumulants beyond rotational invariance

Resumen Técnico: Estadísticas Lineales en el Gas de Coulomb: Cumulantes Más Allá de la Invarianza Rotacional

Planteamiento del Problema
El artículo investiga las fluctuaciones de las estadísticas lineales en un gas de Coulomb clásico de $N$ partículas en dimensión $d \ge 2$, confinado por un potencial externo general $U$. Las partículas interactúan mediante una repulsión logarítmica ($d=2$) o newtoniana ($d=3$). El objeto principal de estudio es la estadística lineal $L_N = \sum_{i=1}^N f(x_i)$, donde $f$ es una función de prueba suave.

Aunque la varianza de orden principal de $L_N$ está bien comprendida y está gobernada por un teorema del límite central con fluctuaciones gaussianas (escalando como $O(1)$), el comportamiento de los cumulantes de orden superior para potenciales generales no simétricos y funciones de prueba sigue siendo ampliamente inexplorado. Los resultados previos para cumulantes superiores (por ejemplo, en [29]) estaban restringidos a entornos con invarianza rotacional donde la gota de equilibrio (el soporte de la densidad de partículas) permanece como un balón bajo perturbación. Este trabajo tiene como objetivo derivar fórmulas explícitas para los principales cumulantes de tercer y cuarto orden para gotas suaves generales $D_0$ sin asumir invarianza rotacional ni para el potencial $U$ ni para la función de prueba $f$.

Metodología
Los autores emplean una combinación de cálculo de formas, teoría de potenciales y técnicas de mapeo conforme:

1. Potencial Inclinado y Función Generadora de Cumulantes: Los cumulantes se derivan de la función generadora de cumulantes (CGF) $\chi(s, N) = \log \langle e^{-Ns L_N} \rangle_U$. Al introducir un potencial inclinado $U_s = U + s f / \beta$, el problema se reduce al análisis de la función de partición $Z_{N, U_s}$. El $q$-ésimo cumulante se obtiene diferenciando $\chi(s, N)$ $q$ veces con respecto a $s$ en $s=0$.
2. Problema del Obstáculo y Velocidad de la Frontera: La esencia de la derivación reside en la conexión entre la medida de equilibrio y el problema del obstáculo. A medida que el parámetro $s$ varía, el soporte de la densidad de equilibrio (la gota $D_s$) se deforma. Los autores utilizan el resultado conocido de que la velocidad normal $v_s$ de la frontera de la gota $\partial D_s$ está determinada por el salto en la derivada normal de la extensión armónica de la función de prueba. Específicamente, $v_s = [\vec{\nabla} f_{D_s}] \cdot n_s / (\beta \Delta U_s)$.
3. Derivadas de Forma: Los autores calculan las derivadas de la energía de Dirichlet asociada a la extensión armónica de $f$ con respecto a la deformación del dominio. Esto implica calcular la variación de las integrales de superficie sobre la frontera móvil $\partial D_s$, utilizando la fórmula de Hadamard y las propiedades del operador de Dirichlet-al-Neumann exterior ($\Lambda_+$).
4. Mapeo Conforme ($d=2$): Para el caso bidimensional, se proporciona una derivación independiente mediante análisis complejo. Este enfoque utiliza el mapa de Riemann exterior, la fórmula de variación de Julia para la evolución del mapa conforme bajo el movimiento de la frontera, y la transformada de Hilbert para relacionar la velocidad de la frontera con la función de prueba.
5. Justificación Rigurosa ($d=2, \beta=2$): La Sección 4 proporciona una prueba rigurosa de las fórmulas asintóticas para los terceros y cuartos cumulantes en el caso determinantal ($d=2, \beta=2$) con datos real-analíticos. Esto se basa en la asintótica de polinomios ortogonales y la expansión de Bergman en el interior (bulk) para justificar el intercambio de límites ($N \to \infty$ y $s \to 0$).

Contribuciones Clave y Resultados

1. Fórmula General del Tercer Cumulante:
   El artículo deriva una fórmula explícita para el tercer cumulante principal $\langle L_N^3 \rangle_c$ para cualquier dimensión $d \ge 2$ y una gota $D_0$ suave general. El resultado se expresa como una integral de frontera intrínseca sobre la gota de equilibrio:
   $$\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}$$
   donde $[\vec{\nabla} f_{D_0}] \cdot n$ es el salto de la derivada normal de la extensión armónica de $f$ a través de la frontera, $n$ es la normal unitaria exterior y $\Delta U$ es el Laplaciano del potencial externo (proporcional a la densidad de fondo). La fórmula depende explícitamente del cubo de este salto.

2. Fórmula General del Cuarto Cumulante:
   El cuarto cumulante se deriva como una integral de frontera que involucra operadores geométricos y espectrales adicionales:
   $$\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]$$
   donde $K = [\vec{\nabla} f_{D_0}] \cdot n$, $\Lambda_+$ es el operador de Dirichlet-al-Neumann exterior, $M$ es la curvatura media de la frontera y $\partial_n \Delta U$ es la derivada normal de la densidad de fondo. A diferencia del tercer cumulante, el cuarto involucra la acción del operador de Dirichlet-al-Neumann sobre el cuadrado del salto, así como términos de curvatura y gradiente de densidad.

3. Geometrías Especiales:
   Las fórmulas generales se especializan para geometrias específicas, incluyendo:

   • Disco ($d=2$): Expresado en términos de la transformada de Hilbert y modos de Fourier.
   • Semiplano ($d=2$): Simplificado para involucrar la transformada de Hilbert de la derivada de $f$.
   • Esfera ($d=3$): Expresado utilizando armónicos esféricos y los autovalores del operador de Dirichlet-al-Neumann exterior.
     Estos resultados especializados recuperan las fórmulas radiales conocidas cuando $f$ es rotacionalmente invariante.

4. Validación:
   Los autores verifican sus fórmulas generales contra deformaciones no radiales exactamente solubles:

   • Gotas elípticas trasladas en $d=2$ (Apéndice A).
   • Deformaciones esferoidales de un balón en $d=3$ (Apéndice B).
     En el caso $d=2, \beta=2$, se muestra que los resultados coinciden con las fórmulas propuestas por Wiegmann y Zabrodin [51], proporcionando una prueba rigurosa para los terceros y cuartos cumulantes en ese entorno específico.

Significancia y Reivindicaciones
El artículo afirma proporcionar las primeras fórmulas explícitas para los principales terceros y cuartos cumulantes de estadísticas lineales suaves en gases de Coulomb para potenciales no rotacionalmente invariantes generales y funciones de prueba en dimensiones arbitrarias $d \ge 2$.

• Más allá de las Fluctuaciones Gaussianas: Mientras que la varianza es bien conocida por ser $O(1)$ y gaussiana, estos cumulantes superiores representan las primeras correcciones no gaussianas. Las fórmulas derivadas muestran explícitamente los factores de escala $N^{-1}$ y $N^{-2}$ en los integrandos, consistentes con el comportamiento observado en casos de invarianza rotacional, pero la novedad principal reside en la compleja estructura integral que involucra saltos de frontera, operadores de Dirichlet-al-Neumann y curvatura.
• Dependencia Intrínseca de la Frontera: Un hallazgo clave es que estos cumulantes superiores dependen únicamente de los valores de frontera de la función de prueba y sus derivadas, así como de la geometría de la frontera de la gota (curvatura, derivadas normales de la densidad). Esto extiende la "universalidad del borde" observada en la varianza a órdenes superiores.
• Fundamento Riguroso: Para el caso $d=2, \beta=2$, el artículo ofrece una derivación rigurosa que valida resultados heurísticos previos de la literatura de física (Wiegmann y Zabrodin).
• Limitaciones: Los autores señalan que la derivación general asume que las gotas perturbadas permanecen suaves y no experimentan transiciones topológicas (por ejemplo, división o fusión). La prueba rigurosa para $d$ general y $\beta$ general depende de una condición suficiente respecto a las derivadas de la función de partición (Observación 4.6), la cual no está plenamente probada para todos los casos, pero se establece para el caso determinantal.

Este trabajo une la mecánica estadística, la teoría de matrices aleatorias y los problemas de frontera libre (problemas del obstáculo), ofreciendo un marco unificado para comprender las fluctuaciones no gaussianas en sistemas de partículas más allá de las configuraciones simétricas.