Feedback Control and Local Convexification of Wasserstein Gradient Flows

1. Planteamiento del Problema

El trabajo aborda el comportamiento dinámico de sistemas de partículas interactuantes modelados mediante el límite de campo medio, descrito por ecuaciones de Fokker-Planck no locales (dinámicas de McKean-Vlasov). Estas dinámicas se interpretan como flujos gradiente de Wasserstein en el espacio de medidas de probabilidad, minimizando una energía libre $F(\mu)$ de la forma:
$$F(\mu) = E(\mu) + \sigma \int_\Omega \mu \log \mu \, dx$$
donde $E(\mu)$ es una energía interna (que puede ser no convexa) y el término de entropía introduce difusión.

El problema central:
En muchos escenarios prácticos (como medios granulares o redes neuronales), el paisaje de energía libre es no convexo, presentando múltiples mínimos locales y puntos de equilibrio inestables (puntos de silla o máximos locales).

• En estos puntos inestables, el Hessiano de Wasserstein de la energía tiene autovalores negativos, lo que impide la convergencia local del flujo gradiente no perturbado.
• Incluso en casos convexos, la tasa de convergencia puede ser extremadamente lenta debido a un "hueco espectral" (spectral gap) pequeño.
• La teoría actual describe este comportamiento metastable pero carece de un mecanismo constructivo para forzar la convexidad local y estabilizar el sistema en puntos de equilibrio inestables o mejorar su tasa de relajación.

Objetivo: Diseñar un control por retroalimentación (feedback) de rango finito que modifique la dinámica para:

1. Desplazar el espectro del Hessiano de Wasserstein en un estado estacionario $\bar{\mu}$ por encima de un umbral positivo $\delta > 0$.
2. Garantizar la estabilización exponencial local del sistema no lineal hacia $\bar{\mu}$.
3. Lograr la convexificación local del paisaje de energía libre en coordenadas de transporte.

---

2. Metodología y Marco Teórico

Los autores desarrollan un marco teórico de operadores que conecta la geometría del espacio de Wasserstein con la teoría de control óptimo.

A. Realización del Hessiano de Wasserstein

El primer paso clave es identificar el Hessiano de Wasserstein de la energía libre $F$ en un punto estacionario $\bar{\mu}$ con un operador autoadjunto bien definido.

• Se trabaja en el espacio de variables potenciales $X = H^1_{zm}(\bar{\mu})$ (funciones con media cero respecto a $\bar{\mu}$).
• Se demuestra que el Hessiano corresponde a un operador autoadjunto $A$ en $X$ con resolvente compacto.
• El espectro de $A$ consiste en autovalores reales $\lambda_k \to +\infty$. Los autovalores negativos indican modos inestables.

B. Conexión con la Dinámica Linealizada

Se establece una relación de conjugación entre el Hessiano y el generador de la dinámica linealizada en variables de densidad:
$$L = -I_{\bar{\mu}} A I_{\bar{\mu}}^{-1}$$
donde $I_{\bar{\mu}}$ es la isometría de Riesz inducida por la estructura $H^1$ ponderada. Esto implica que controlar la dinámica linealizada es equivalente a modificar el espectro del operador Hessiano.

C. Diseño del Control por Retroalimentación

Para estabilizar el sistema, se propone un control de rango finito que actúa sobre los modos inestables (autovalores $\lambda_k < \delta$):

1. Perturbación: Se introduce un campo de fuerza controlado de la forma $\sum u_j(t) \nabla \cdot (\mu \nabla \alpha_j)$, donde $\alpha_j$ son perfiles espaciales (eigenfunciones del Hessiano).
2. Síntesis Óptima: Se formula un problema de control lineal cuadrático (LQR) con un funcional de costo ponderado en el tiempo.
3. Ecuación Riccati: Se resuelve una Ecuación Algebraica de Riccati (ARE) para obtener una ganancia de retroalimentación $\Pi$.
   $$u(t) = B^* \Pi \xi(t)$$
   Esto modifica el operador original $A$ a un operador en lazo cerrado $A_\Pi = A + BB^*\Pi$.

D. Estabilización No Lineal

El teorema principal demuestra que el desplazamiento espectral (convexificación infinitesimal) se traduce en estabilización exponencial local para la ecuación de Fokker-Planck completa no lineal. Esto se logra mediante:

• Estimaciones de regularidad en el resto no lineal (en coordenadas de transporte).
• Argumentos de "bootstrap" que utilizan la contracción del mapa proximal en un entorno pequeño de la variedad de transporte.

---

3. Contribuciones Clave

1. Identificación Espectral del Hessiano: Demostración rigurosa de que el Hessiano de Wasserstein en un punto estacionario es un operador autoadjunto con resolvente compacto en un espacio de Hilbert de potencias, permitiendo el uso de teoría espectral estándar.
2. Convexificación Constructiva: Propuesta de un mecanismo de control por retroalimentación que desplaza activamente el espectro del Hessiano, convirtiendo un punto de silla o inestable en un mínimo local estrictamente convexo con una tasa de convergencia predefinida $\delta$.
3. Estabilidad No Lineal Local: Prueba de que la estabilización lineal (espectral) implica la estabilidad exponencial global para la PDE no lineal completa, bajo supuestos de regularidad adecuados.
4. Generalización: Extensión del marco a:
   • Sistemas de múltiples especies (acoplados).
   • Ecuaciones de Fokker-Planck con restricciones de momentos lineales.
   • Variedades Riemannianas cerradas (incorporando la curvatura de Ricci y el tensor de Bakry-Émery).

---

4. Resultados Principales

• Teorema de Desplazamiento Espectral (Teorema C): Para cualquier $\delta > 0$, existe un operador de retroalimentación de rango finito tal que el Hessiano modificado satisface $Hess_{W^2} F_{cl}(\bar{\mu}) \geq \delta I$. Esto garantiza que la energía libre modificada es localmente fuertemente convexa.
• Teorema de Estabilización Exponencial (Teorema D): Bajo condiciones de regularidad, la ley de control produce una solución única global que converge exponencialmente a $\bar{\mu}$ con tasa al menos $\delta$:
  $$\|\mu_t - \bar{\mu}\| \leq C e^{-\delta t} \|\mu_0 - \bar{\mu}\|$$
• Validación Numérica:
  • Potencial de doble pozo: Se muestra cómo el control elimina la lentitud del túnel cuántico/metastable, acelerando la relajación.
  • Modelo de Kuramoto: Estabilización de un estado incoherente inestable en un sistema de osciladores acoplados.
  • Variedad Esférica ($S^2$): Estabilización en una esfera con un potencial no convexo donde la curvatura de Bakry-Émery es negativa en el ecuador, demostrando que el control puede superar la falta de convexidad geométrica intrínseca.

---

5. Significado e Impacto

Este trabajo es significativo por varias razones:

• Puente entre Geometría y Control: Conecta la geometría de Wasserstein (convexidad de desplazamiento) con la teoría de control de sistemas distribuidos (PDEs), ofreciendo una interpretación física de la convexificación como un control activo.
• Solución a la Metastabilidad: Proporciona una herramienta teórica y práctica para evitar que los sistemas de partículas queden atrapados en mínimos locales o estados inestables, un problema común en inferencia variacional, aprendizaje automático y física estadística.
• Marco Unificado: La metodología basada en operadores y ecuaciones de Riccati es lo suficientemente general para aplicarse a sistemas complejos (multiespecie, restricciones, geometrías curvas), superando las limitaciones de los métodos de control tradicionales que a menudo se restringen a espacios euclidianos o dinámicas lineales.
• Aplicaciones Potenciales: Tiene implicaciones directas en el diseño de algoritmos de muestreo más rápidos (Langevin), optimización de redes neuronales (evitando mínimos locales) y control de sistemas de agentes colectivos.

En resumen, el artículo establece que la convexidad local de un paisaje de energía no convexo puede ser forzada artificialmente mediante un control de retroalimentación óptimo, garantizando así una convergencia rápida y robusta hacia el equilibrio deseado.