Universal Central Limit Theorem for non-exchangeable interacting diffusions

Resumen Técnico: Teorema del Límite Central Universal para Difusiones Interactuantes No Intercambiables

Planteamiento del Problema
El artículo investiga el comportamiento de gran-$n$ de un sistema de $n$ difusiones interactuantes $(X^1, \dots, X^n)$ en $\mathbb{R}^d$, gobernadas por la ecuación 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,$$
donde $\xi = (\xi_{ij})$ es una matriz de interacción determinista con entradas no negativas y diagonal cero. Si bien el régimen de la "Ley de los Grandes Números" (propagación del caos) para tales sistemas no intercambiables ha sido estudiado extensamente, particularmente la convergencia de la medida empírica hacia un límite de McKean-Vlasov, el comportamiento de las fluctuaciones alrededor de este límite permanece menos comprendido para estructuras de interacción heterogéneas.

La cuestión central abordada es si el campo de fluctuación global, definido como $\eta^n_t := \sqrt{n}(\mu^n_t - \mu_t)$ (donde $\mu^n_t$ es la medida empírica y $\mu_t$ el límite determinista), converge al mismo límite gaussiano que en el caso clásico de campo medio intercambiable (donde $\xi_{ij} = 1/(n-1)$). Específicamente, los autores buscan determinar las condiciones sobre la matriz de interacción $\xi$ bajo las cuales $\eta^n$ converge en distribución al solución única de una EDP estocástica lineal.

Metodología
La demostración sigue una estrategia de "compacidad-límite-unicidad", realizada principalmente dentro de espacios de Sobolev negativos $H^{-k}$ para manejar la naturaleza distributiva del campo de fluctuación.

1. Compacidad: Los autores establecen la compacidad de la secuencia $(\eta^n)_{n \in \mathbb{N}}$ en el espacio $C([0, T]; H^{-k})$. Esto se logra mediante:

   • La demostración de la compacidad de los márgenes temporales utilizando análisis de Fourier y cotas de $H^1$ con pesos uniformes.
   • La verificación del criterio de Aldous para la compacidad en el espacio de trayectorias.
   • Crucialmente, los autores utilizan estimaciones cuantitativas agudas de propagación del caos (de Shkolnikov [58]) combinadas con desigualdades funcionales (Poincaré y de transporte) para controlar la varianza del campo de fluctuación. Deponen el segundo momento del emparejamiento de la fluctuación en un término de varianza (controlado por desigualdades de Poincaré) y un término de sesgo (controlado por la entropía relativa entre las leyes de las partículas y el límite de McKean-Vlasov).

2. Convergencia (Identificación del Límite): Para identificar cualquier subsecuencia límite como una solución de la EDP estocástica objetivo, los autores analizan la descomposición de semimartingala de $\langle \eta^n_t, \phi \rangle$ para funciones de prueba $\phi$.

   • La parte de martingala converge a un término de ruido gaussiano caracterizado por la variación cuadrática.
   • El término de deriva involucra el operador linealizado de McKean-Vlasov. Un desafío técnico significativo surge porque los coeficientes en la deriva son acotados pero no necesariamente decaen en el infinito, mientras que la convergencia de $\eta^n$ es solo en el espacio sin peso $H^{-k}$.
   • Para resolver esto, los autores emplean un argumento de corte: localizan las funciones de prueba en conjuntos compactos donde la convergencia sin peso es válida y controlan las colas utilizando estimaciones de momentos con pesos uniformes derivadas de las propiedades del núcleo de Bessel y las desigualdades de transporte.
   • Muestran que los términos de error que surgen de la heterogeneidad de $\xi$ desaparecen en el límite bajo las condiciones de densidad especificadas.

3. Unicidad: El artículo demuestra la unicidad por trayectoria para la EDP estocástica límite. Al considerar la diferencia $\delta_t = \eta_t - \eta'_t$ entre dos soluciones impulsadas por el mismo ruido, los términos de ruido se cancelan, dejando una ecuación lineal determinista. Los autores derivan una desigualdad de energía en una norma de Sobolev negativa con peso y aplican la desigualdad de Gronwall para mostrar que $\delta_t \equiv 0$.

Supuestos y Condiciones Clave
El resultado principal se sostiene para matrices de interacción $\xi$ que satisfacen la Suposición 2.1, la cual incluye:

• Sumas de Filas: $\sum_j \xi_{ij} = 1$ para todo $i$.
• Sumas de Columnas: Sumas de columnas acotadas, con la suma al cuadrado promedio convergiendo a 1.
• Umbral de Densidad: Las entradas de la matriz deben satisfacer una condición de decaimiento específica. La condición más crítica es que la entrada máxima de $\xi$ debe ser $o(n^{-1/2})$. Específicamente, se requiere la condición $\frac{1}{n^{3/2}} (\sum_{i,j} |\tilde{\xi}_{ij}|) \max_{i,j} \xi_{ij} \to 0$, donde $\tilde{\xi}$ captura la desviación respecto al caso de campo medio.
• Regularidad: Las funciones de interacción $b_0$ y $b$ deben tener derivadas acotadas hasta orden $k+1$.

Resultados Principales
El artículo establece un Teorema del Límite Central Universal (Teorema 2.4):

• Bajo las condiciones estructurales y de densidad establecidas para $\xi$, el campo de fluctuación $\eta^n$ converge en distribución en $C([0, T]; H^{-k})$ a un proceso gaussiano único $\eta$.
• Este límite $\eta$ es la solución única de la EDP estocástica lineal:
  $$\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 es idéntica al límite de fluctuación en el caso de campo medio intercambiable.
• El resultado se aplica a:
  • Matrices de adyacencia escaladas de grafos regulares $m_n$ donde $m_n / \sqrt{n} \to \infty$.
  • Grafos de Erdős–Rényi $G(n, p_n)$ donde $n p_n / \sqrt{n} \to \infty$ (en el sentido quenched).
  • Modelos de interacción espacial en el régimen subcrítico ($\alpha < 1/2$).

Significancia y Agudeza
El artículo afirma que el umbral de $n^{-1/2}$ para la entrada máxima de la matriz de interacción es agudo.

• En el contexto de los modelos de interacción espacial (donde la fuerza de interacción decae como $|x_i - x_j|^{-\alpha}$), la condición $\max \xi_{ij} \sim n^{\alpha-1}$ implica que el umbral de $n^{-1/2}$ corresponde a $\alpha < 1/2$.
• Los autores señalan que en el régimen supercrítico ($\alpha > 1/2$), las fluctuaciones, cuando se reescalan por $\sqrt{n}$, divergen, y un escalamiento diferente ($n^{1-\alpha}$) conduce a un límite determinista en lugar de uno gaussiano.
• La universalidad del límite gaussiano es significativa porque demuestra que, siempre que la gráfica de interacción sea suficientemente densa (específicamente, que la fuerza de interacción no esté demasiado concentrada en un solo par), la heterogeneidad microscópica de la red no altera las estadísticas de fluctuación macroscópica; el sistema se comporta como si fuera totalmente intercambiable.

Conclusión
Este trabajo extiende la teoría clásica de fluctuaciones de las difusiones interactuantes de campo medio a una amplia clase de sistemas no intercambiables. Al combinar la propagación del caos cuantitativa con desigualdades funcionales y un análisis cuidadoso en espacios de Sobolev con peso, los autores demuestran que la universalidad del límite de fluctuación gaussiano se mantiene siempre que la matriz de interacción satisfaga una condición de densidad específica, la cual se muestra óptima para los modelos espaciales.