Linear statistics in the Coulomb gas: cumulants beyond rotational invariance

Resumo Técnico: Estatísticas Lineares no Gás de Coulomb: Cumulantes Além da Invariância Rotacional

Enunciado do Problema
O artigo investiga as flutuações de estatísticas lineares em um gás de Coulomb clássico de $N$ partículas em dimensão $d \ge 2$, confinado por um potencial externo $U$ geral. As partículas interagem via repulsão logarítmica ($d=2$) ou Newtoniana ($d=3$). O objeto principal de estudo é a estatística linear $L_N = \sum_{i=1}^N f(x_i)$, onde $f$ é uma função de teste suave.

Embora a variância de ordem principal de $L_N$ seja bem compreendida e governada por um teorema do limite central com flutuações Gaussianas (escalonando como $O(1)$), o comportamento de cumulantes de ordem superior para potenciais não simétricos e funções de teste gerais permanece amplamente inexplorado. Resultados anteriores para cumulantes superiores (por exemplo, em [29]) foram restritos a configurações de invariância rotacional onde a gota de equilíbrio (o suporte da densidade de partículas) permanece um disco sob perturbação. Este trabalho visa derivar fórmulas explícitas para os principais cumulantes de terceira e quarta ordem para gotas $D_0$ suaves e gerais, sem assumir invariância rotacional tanto para o potencial $U$ quanto para a função de teste $f$.

Metodologia
Os autores empregam uma combinação de cálculo de forma, teoria do potencial e técnicas de mapeamento conformal:

1. Potencial Inclinado e Função Geradora de Cumulantes: Os cumulantes são derivados da função geradora de cumulantes (CGF) $\chi(s, N) = \log \langle e^{-Ns L_N} \rangle_U$. Ao introduzir um potencial inclinado $U_s = U + s f / \beta$, o problema é reduzido à análise da função de partição $Z_{N, U_s}$. O $q$-ésimo cumulante é obtido diferenciando $\chi(s, N)$ $q$ vezes em relação a $s$ em $s=0$.
2. Problema do Obstáculo e Velocidade de Fronteira: A essência da derivação baseia-se na conexão entre a medida de equilíbrio e o problema do obstáculo. À medida que o parâmetro $s$ varia, o suporte da densidade de equilíbrio (a gota $D_s$) se deforma. Os autores utilizam o resultado conhecido de que a velocidade normal $v_s$ da fronteira da gota $\partial D_s$ é determinada pelo salto no derivada normal da extensão harmônica da função de teste. Especificamente, $v_s = [\vec{\nabla} f_{D_s}] \cdot n_s / (\beta \Delta U_s)$.
3. Derivadas de Forma: Os autores calculam as derivadas da energia de Dirichlet associada à extensão harmônica de $f$ em relação à deformação do domínio. Isso envolve o cálculo de integrais de superfície sobre a fronteira móvel $\partial D_s$, utilizando a fórmula de Hadamard e as propriedades do operador Dirichlet-para-Neumann exterior ($\Lambda_+$).
4. Mapeamento Conformal ($d=2$): Para o caso bidimensional, uma derivação independente é fornecida usando análise complexa. Esta abordagem utiliza o mapa de Riemann exterior, a fórmula de variação de Julia para a evolução do mapeamento conformal sob movimento de fronteira, e a transformada de Hilbert para relacionar a velocidade de fronteira com a função de teste.
5. Justificativa Rigorosa ($d=2, \beta=2$): A Seção 4 fornece uma prova rigorosa das fórmulas assintóticas para os principais cumulantes de terceira e quarta ordem no caso determinantal ($d=2, \beta=2$) com dados real-analíticos. Isso se baseia em assíntotas de polinômios ortogonais e na expansão de Bergman no bulk para justificar a troca de limites ($N \to \infty$ e $s \to 0$).

Contribuições Principais e Resultados

1. Fórmula Geral do Terceiro Cumulante:
   O artigo deriva uma fórmula explícita para o principal terceiro cumulante $\langle L_N^3 \rangle_c$ para qualquer dimensão $d \ge 2$ e gota $D_0$ suave e geral. O resultado é expresso como uma integral de fronteira intrínseca sobre a gota de equilíbrio:
   $$\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}$$
   onde $[\vec{\nabla} f_{D_0}] \cdot n$ é o salto da derivada normal da extensão harmônica de $f$ através da fronteira, $n$ é a normal unitária externa, e $\Delta U$ é o Laplaciano do potencial externo (proporcional à densidade de fundo). A fórmula depende explicitamente do cubo deste salto.

2. Fórmula Geral do Quarto Cumulante:
   O quarto cumulante é derivado como uma integral de fronteira envolvendo operadores geométricos e espectrais adicionais:
   $$\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]$$
   onde $K = [\vec{\nabla} f_{D_0}] \cdot n$, $\Lambda_+$ é o operador Dirichlet-para-Neumann exterior, $M$ é a curvatura média da fronteira, e $\partial_n \Delta U$ é a derivada normal da densidade de fundo. Diferente do terceiro cumulante, o quarto envolve a ação do operador Dirichlet-para-Neumann sobre o quadrado do salto, bem como termos de curvatura e gradiente de densidade.

3. Geometrias Especiais:
   As fórmulas gerais são especializadas para geometrias específicas, incluindo:

   • Disco ($d=2$): Expresso em termos da transformada de Hilbert e modos de Fourier.
   • Semiplano ($d=2$): Simplificado para envolver a transformada de Hilbert da derivada de $f$.
   • Esfera ($d=3$): Expresso usando harmônicos esféricos e os autovalores do operador Dirichlet-para-Neumann exterior.
     Estes resultados especializados recuperam as fórmulas radiais conhecidas quando $f$ é rotacionalmente invariante.

4. Validação:
   Os autores verificam suas fórmulas gerais contra deformações não radiais exatamente solúveis:

   • Gotas elípticas transladadas em $d=2$ (Apêndice A).
   • Deformações esferoidais de uma bola em $d=3$ (Apêndice B).
     No caso de $d=2, \beta=2$, mostra-se que os resultados coincidem com as fórmulas propostas por Wiegmann e Zabrodin [51], fornecendo uma prova rigorosa para os terceiros e quartos cumulantes neste cenário específico.

Significância e Alegações
O artigo afirma fornecer as primeiras fórmulas explícitas para os principais terceiros e quartos cumulantes de estatísticas lineares suaves em gases de Coulomb para potenciais gerais, não rotacionalmente invariantes, e funções de teste em dimensões arbitrárias $d \ge 2$.

• Além das Flutuações Gaussianas: Embora a variância seja bem conhecida como sendo $O(1)$ e Gaussiana, estes cumulantes de ordem superior representam as primeiras correções não Gaussianas. As fórmulas derivadas mostram explicitamente os fatores de escala $N^{-1}$ e $N^{-2}$ nos integrandos, consistentes com o comportamento observado em casos de invariância rotacional, mas a principal novidade reside na estrutura integral complexa envolvendo saltos de fronteira, operadores Dirichlet-para-Neumann e curvatura.
• Dependência Intrínseca da Fronteira: Uma descoberta fundamental é que estes cumulantes de ordem superior dependem apenas dos valores de fronteira da função de teste e de suas derivadas, bem como da geometria da fronteira da gota (curvatura, derivadas normais da densidade). Isso estende a "universalidade de borda" observada na variância para ordens superiores.
• Fundação Rigorosa: Para o caso $d=2, \beta=2$, o artigo oferece uma derivação rigorosa que valida resultados heurísticos anteriores da literatura de física (Wiegmann e Zabrodin).
• Limitações: Os autores observam que a derivação geral assume que as gotas perturbadas permanecem suaves e não passam por transições topológicas (por exemplo, divisão ou fusão). A prova rigorosa para $d$ geral e $\beta$ depende de uma condição suficiente em relação às derivadas da função de partição (Observação 4.6), que não é totalmente provada para todos os casos, mas está estabelecida para o caso determinantal.

Este trabalho une mecânica estatística, teoria de matrizes aleatórias e problemas de fronteira livre (problemas de obstáculo), oferecendo um arcabouço unificado para compreender flutuações não Gaussianas em sistemas de partículas além de configurações simétricas.