Universal Central Limit Theorem for non-exchangeable interacting diffusions

Resumo Técnico: Teorema do Limite Central Universal para Difusões Interagentes Não-Excambiáveis

Enunciado do Problema
O artigo investiga o comportamento de grande-$n$ de um sistema de $n$ difusões interagentes $(X^1, \dots, X^n)$ em $\mathbb{R}^d$, governadas pela equação diferencial estocástica:
$$dX^i_t = \left( b_0(t, X^i_t) + \sum_{j=1}^n \xi_{ij} b(t, X^i_t, X^j_t) \right) dt + \sigma dB^i_t,$$
onde $\xi = (\xi_{ij})$ é uma matriz de interação determinística com entradas não-negativas e diagonal zero. Embora o regime da "Lei dos Grandes Números" (propagação do caos) para tais sistemas não-excambiáveis tenha sido extensivamente estudado, particularmente a convergência da medida empírica para um limite de McKean-Vlasov, o comportamento das flutuações em torno desse limite permanece menos compreendido para estruturas de interação heterogêneas.

A questão central abordada é se o campo de flutuação global, definido como $\eta^n_t := \sqrt{n}(\mu^n_t - \mu_t)$ (onde $\mu^n_t$ é a medida empírica e $\mu_t$ o limite determinístico), converge para o mesmo limite Gaussiano que no caso clássico de campo médio excambiável (onde $\xi_{ij} = 1/(n-1)$). Especificamente, os autores buscam determinar as condições sobre a matriz de interação $\xi$ sob as quais $\eta^n$ converge em distribuição para a solução única de uma EDP estocástica (SPDE) linear.

Metodologia
A prova segue uma estratégia de "aperto-limite-unicidade" (tightness-limit-uniqueness), conduzida primordialmente em espaços de Sobolev negativos $H^{-k}$ para lidar com a natureza distributiva do campo de flutuação.

1. Aperto (Tightness): Os autores estabelecem o aperto da sequência $(\eta^n)_{n \in \mathbb{N}}$ no espaço $C([0, T]; H^{-k})$. Isso é alcançado através de:

   • Prova de aperto dos marginais temporais usando análise de Fourier e limites $H^1$ ponderados uniformemente.
   • Verificação do critério de Aldous para o aperto no espaço de caminhos.
   • Crucialmente, os autores utilizam estimativas quantitativas aguçadas de propagação do caos (de Shkolnikov [58]) combinadas com desigualdades funcionais (Poincaré e desigualdades de transporte) para controlar a variância do campo de flutuação. Eles decompõem o segundo momento do emparelhamento da flutuação em um termo de variância (controlado por desigualdades de Poincaré) e um termo de viés (controlado pela entropia relativa entre as leis das partículas e o limite de McKean-Vlasov).

2. Convergência (Identificação do Limite): Para identificar qualquer limite subsequencial como uma solução da SPDE alvo, os autores analisam a decomposição semimartingal de $\langle \eta^n_t, \phi \rangle$ para funções de teste $\phi$.

   • A parte de martingala converge para um termo de ruído Gaussiano caracterizado pela variação quadrática.
   • O termo de deriva (drift) envolve o operador de McKean-Vlasov linearizado. Um desafio técnico significativo surge porque os coeficientes na deriva são limitados, mas não necessariamente decaem no infinito, enquanto a convergência de $\eta^n$ ocorre apenas no espaço não-ponderado $H^{-k}$.
   • Para resolver isso, os autores empregam um argumento de corte (cutoff): eles localizam funções de teste em conjuntos compactos onde a convergência não-ponderada ocorre e controlam as caudas usando estimativas de momentos ponderados uniformes derivadas das propriedades do núcleo de Bessel e desigualdades de transporte.
   • Eles mostram que os termos de erro decorrentes da heterogeneidade de $\xi$ desaparecem no limite sob as condições de densidade especificadas.

3. Unicidade: O artigo prova a unicidade por trajetória (pathwise uniqueness) para a SPDE limite. Ao considerar a diferença $\delta_t = \eta_t - \eta'_t$ entre duas soluções impulsionadas pelo mesmo ruído, os termos de ruído se cancelam, restando uma equação linear determinística. Os autores derivam uma desigualdade de energia em uma norma de Sobolev negativa ponderada e aplicam a desigualdade de Gronwall para mostrar que $\delta_t \equiv 0$.

Premissas e Condições Principais
O resultado principal é válido para matrizes de interação $\xi$ que satisfazem a Premissa 2.1, que inclui:

• Somas das Linhas: $\sum_j \xi_{ij} = 1$ para todo $i$.
• Somas das Colunas: Somas de colunas limitadas, com a média do quadrado das somas das colunas convergindo para 1.
• Limiar de Densidade: As entradas da matriz devem satisfazer uma condição de decaimento específica. A condição mais crítica é que a entrada máxima de $\xi$ deve ser $o(n^{-1/2})$. Especificamente, a condição $\frac{1}{n^{3/2}} (\sum_{i,j} |\tilde{\xi}_{ij}|) \max_{i,j} \xi_{ij} \to 0$ é necessária, onde $\tilde{\xi}$ captura o desvio do caso de campo médio.
• Regularidade: As funções de interação $b_0$ e $b$ devem ter derivadas limitadas até ordem $k+1$.

Principais Resultados
O artigo estabelece um Teorema do Limite Central Universal (Teorema 2.4):

• Sob as condições estruturais e de densidade declaradas para $\xi$, o campo de flutuação $\eta^n$ converge em distribuição em $C([0, T]; H^{-k})$ para um processo Gaussiano único $\eta$.
• Este limite $\eta$ é a solução única da SPDE linear:
  $$\partial_t \eta_t = \frac{\sigma^2}{2} \Delta \eta_t - \nabla \cdot (b_0 \eta_t) - \nabla \cdot \left( \eta_t \int b(\cdot, y) \mu_t(dy) \right) - \nabla \cdot \left( \mu_t \int b(\cdot, y) \eta_t(dy) \right) - \sigma \nabla \cdot (\sqrt{\mu_t} \dot{W}),$$
  que é idêntica ao limite de flutuação no caso excambiável de campo médio.
• O resultado aplica-se a:
  • Matrizes de adjacência escalonadas de grafos $m_n$-regulares onde $m_n / \sqrt{n} \to \infty$.
  • Grafos Erdős–Rényi $G(n, p_n)$ onde $n p_n / \sqrt{n} \to \infty$ (no sentido quenched).
  • Modelos de interação espacial no regime subcrítico ($\alpha < 1/2$).

Significância e Precisão (Sharpness)
O artigo afirma que o limiar de $n^{-1/2}$ para a entrada máxima da matriz de interação é preciso (sharp).

• No contexto de modelos de interação espacial (onde a força de interação decai como $|x_i - x_j|^{-\alpha}$), a condição $\max \xi_{ij} \sim n^{\alpha-1}$ implica que o limiar de $n^{-1/2}$ corresponde a $\alpha < 1/2$.
• Os autores observam que no regime supercrítico ($\alpha > 1/2$), as flutuações, quando reescalonadas por $\sqrt{n}$, divergem, e um escalonamento diferente ($n^{1-\alpha}$) leva a um limite determinístico em vez de um Gaussiano.
• A universalidade do limite Gaussiano é significativa porque demonstra que, desde que a rede de interação seja suficientemente densa (especificamente, a força de interação não esteja muito concentrada em um único par), a heterogeneidade microscópica da rede não altera as estatísticas de flutuação macroscópicas; o sistema comporta-se como se fosse totalmente excambiável.

Conclusão
Este trabalho estende a teoria clássica de flutuações de difusões interagentes de campo médio para uma ampla classe de sistemas não-excambiáveis. Ao combinar a propagação do caos quantitativa com desigualdades funcionais e uma análise cuidadosa em espaços de Sobolev ponderados, os autores provam que a universalidade do limite de flutuação Gaussiana mantém-se desde que a matriz de interação satisfaça uma condição de densidade específica, a qual é mostrada como ótima para modelos espaciais.