Numerical stationary states for nonlocal Fokker-Planck equations via fixed points of consistency maps

Resumen Técnico: Estados Estacionarios Numéricos para Ecuaciones de Fokker-Planck No Locales mediante Puntos Fijos de Mapas de Consistencia

Planteamiento del Problema
El artículo aborda la computación numérica de estados estacionarios para una clase de ecuaciones de tipo Fokker-Planck no locales de la forma:
$$\partial_t u = \nabla \cdot (u\nabla(H'[u] + \Phi[u](x) + V(x))),$$
donde $u(x,t)$ es una densidad, $H(u)$ representa la energía interna, $V(x)$ es un potencial externo y $\Phi[u]$ codifica interacciones no locales, típicamente no lineales. Estas ecuaciones surgen en diversos campos como la dinámica de fluidos, la separación de fases, la adhesión celular y la agregación biológica.

El desafío principal reside en la dinámica de largo tiempo de estos sistemas. Los resolvedores dependientes del tiempo estándar (por ejemplo, esquemas de volumen finito) están limitados a la detección de estados dinámicamente estables. A menudo fallan al identificar estados estacionarios inestables o ramas de bifurcación complejas, particularmente en sistemas que exhiben metaestabilidad o múltiples escalas de tiempo. Además, el software de continuación numérica existente suele tener dificultades con los requisitos estructurales específicos de la difusión no lineal general y los términos no locales presentes en esta clase de ecuaciones.

Metodología
Los autores proponen un marco numérico basado en puntos fijos que evita la discretización directa de los operadores diferenciales. En su lugar, el problema estacionario se reformula como una ecuación de punto fijo no lineal $T u - u = 0$, donde $T$ es un mapa derivado de la PDE original y sus términos de interacción no local.

1. Reformulación de Punto Fijo:
   Al establecer $\partial_t u = 0$ y utilizar la conservación de la masa y la positividad, se demuestra que el estado estacionario $u^*$ satisface una condición de consistencia en su soporte: $H'[u] + \Phi[u] + V = c_0$. Asumiendo que $H'$ es invertible, la solución estacionaria se expresa como un punto fijo del mapa:
   $$T u := (H')^{-1}[c_0(u) - \Phi[u](x) - V(x)].$$
   Para la difusión lineal ($H[u] = \sigma u(\log u - 1)$), este inverso es una función exponencial.

2. Estrategia de Solución Numérica:
   Los autores resuelven el sistema discretizado $F_n(u) = T_n(u) - u = 0$ utilizando un método de Newton-Krylov libre de Jacobiano (JFNK).

   • Bucle Externo: Una iteración de Newton actualiza la solución.
   • Bucle Interno: Un método de subespacio de Krylov (específicamente GMRES) resuelve el sistema lineal $J \delta u = -F(u)$.
   • Enfoque Libre de Matriz: El método no construye explícitamente la matriz Jacobiana. En su lugar, calcula la acción de la derivada de Fréchet $DT[u](\phi)$ sobre un vector $\phi$.
   • Aproximación de la Derivada: Los autores comparan dos enfoques para la derivada de Fréchet:
     • Analítica: Derivar la expresión de la forma cerrada exacta de la derivada del mapa $T$.
     • Diferencia Finita: Aproximar la derivada utilizando diferencias centrales.

3. Discretización:
   Los términos no locales (convoluciones) se manejan mediante Transformadas Rápidas de Fourier (FFT) en dominios periódicos o reglas de cuadratura apropiadas. Esto traslada la carga de precisión de los estenciles de diferenciación a la precisión de la convolución y las reglas de integración.

Contribuciones Clave

• Desarrollo de Marco: Los autores desarrollan y prueban su robusto marco de punto fijo para tres modelos representativos: la ecuación de McKean-Vlasov, un modelo de bandada de Cucker-Smale y una ecuación de Fokker-Planck neuronal.
• Derivadas Analíticas vs. Discretas: El artículo proporciona una comparación directa entre el uso de la derivada de Fréchet analítica y una aproximación de diferencia finita. Demuestra que el enfoque analítico ofrece una precisión y eficiencia computacional garantizadas, aunque el enfoque de diferencia finita sigue siendo viable cuando las formas analíticas son intratables.
• Análisis de Bifurcación: El método reproduce con éxito diagramas de bifurcación conocidos y detecta tanto estados estacionarios estables como inestables, incluyendo ramas supercríticas y subcríticas.
• Descubrimiento de Nuevo Comportamiento: El marco revela estructuras de bifurcación complejas no observadas previamente en estos problemas específicos, particularmente en casos de núcleos multimodales y dominios bidimensionales.

Resultados
Los autores validan su enfoque contra resultados analíticos existentes y teorías de bifurcación conocidas:

• Precisión: Para la ecuación de McKean-Vlasov con un núcleo de tipo Kuramoto, las soluciones numéricas convergen a la solución de referencia analítica con una tasa de decaimiento del error de $O(N_n^{-4})$, donde $N_n$ es el número de puntos de la malla. Los umbrales de bifurcación críticos se recuperan con alta precisión ($\sim 10^{-5}$).
• Núcleos Multimodales: Utilizando núcleos de interacción bimodal, el método identifica múltiples ramas de solución, incluyendo bifurcaciones subcríticas y ramas secundarias que emergen de los puntos de bifurcación primarios. En casos donde el operador lineal tiene un núcleo bidimensional (bifurcación degenerada), el método detecta un comportamiento de tipo transcritico que la teoría local estándar tiene dificultades para caracterizar.
• Estabilidad: Al combinar el resolvedor de punto fijo con un resolvedor dependiente del tiempo para el análisis de perturbaciones, los autores mapean la estabilidad de varias ramas. Observan que, mientras las ramas supercríticas son generalmente estables, las ramas subcríticas son inestables, y ocurren interacciones complejas (por ejemplo, intersecciones de ramas con perfiles distintos).
• Generalizabilidad: El enfoque se extiende con éxito a:
  • Dominios 2D: Demostrando aplicabilidad en el toro $\mathbb{T}^2$.
  • Diferentes Dominios: Manejando dominios no acotados (Cucker-Smale) y dominios semi-infinitos con condiciones de contorno sin flujo (Fokker-Planck neuronal).
  • Diversos Núcleos: Probando núcleos de tipo "top-hat", triangulares y de atracción-repulsión, revelando ricas estructuras de solución incluyendo bifurcaciones secundarias.

Significancia
El artículo sostiene que este enfoque de punto fijo ofrece una excelente herramienta para explorar las soluciones estacionarias de las ecuaciones de Fokker-Planck no locales. Su principal significancia radica en su agnosticismo respecto a la estabilidad dinámica, lo que permite la detección de estados inestables y ramas de bifurcación complejas que los resolvedores de evolución temporal pasan por alto. Al aprovechar las propiedades estructurales de la PDE para recategorizar el problema como un mapa de punto fijo, el método evita las dificultades asociadas con la discretización directa de los operadores diferenciales. Los autores señocan que, si bien el trabajo actual se centra en la difusión lineal, el marco está diseñado para ser extensible a la difusión no lineal (por ejemplo, medio poroso) y sistemas de poblaciones interactuantes, siempre que se cumplan las condiciones estructurales necesarias para la formulación de punto fijo.