The Schrödinger problem on metric graphs

Resumen Técnico: El Problema de Schrödinger en Grafos Métricos

Planteamiento del Problema
Este artículo investiga el problema de Schrödinger en grafos métricos, un problema de interpolación estocástica que se ocupa de encontrar la evolución más probable de partículas de gas independientes entre dos distribuciones observadas, $\mu_0$ y $\mu_1$. El estudio se centra en establecer la relación entre este problema y el transporte óptimo (OT) dentro del entorno geométrico específico de los grafos métricos. A diferencia de las variedades de Riemann o los espacios $RCD^*(K, N)$ donde existe una teoría exhaustiva, los grafos métricos presentan desafíos únicos debido a su falta de suavidad y a la ausencia de fuertes cotas sintéticas de curvatura inferior. Los autores pretenden derivar rigurosamente las formulaciones estáticas y dinámicas del problema de Schrödinger en estos espacios, demostrar su equivalencia y analizar su convergencia hacia la distancia de Wasserstein al cuadrado a medida que el parámetro de regularización se anula.

Metodología
Los autores emplean una combinación de técnicas analíticas y métodos numéricos:

1. Análisis Estático y $\Gamma$-convergencia:

   • El problema estático de Schrödinger se reformula como un problema de transporte óptimo entrópico (EOT) con regularización. La medida de referencia se elige como el núcleo de calor $h_{\beta/2}$ sobre el grafo métrico.
   • Para probar la $\Gamma$-convergencia del problema EOT hacia el problema de transporte óptico de Kantorovich estático cuando $\beta \to 0$, los autores utilizan una técnica de "aproximación por bloques" adaptada de [CDPS17].
   • Crucialmente, aprovechan las cotas superiores e inferiores gaussianas para el núcleo de calor (derivadas de las propiedades del semigrupo y de las fórmulas de suma de caminos) para manejar la función de costo, la cual difiere de la distancia cuadrática estándar debido a la estructura del grafo.

2. Formulación Dinámica y Equivalencia:

   • Se deriva una versión dinámica del problema de Schrödinger (DSP), modelada como una minimización sobre soluciones débiles de la ecuación de continuidad. El funcional incluye un término de energía cinética y un término de información de Fisher (derivado de la regularización de la entropía).
   • La equivalencia entre las formulaciones estática (S$\beta$) y dinámica (DSP) se establece construyendo soluciones específicas utilizando el semigrupo del flujo de calor. Esto implica definir densidades dependientes del tiempo mediante el núcleo de calor y verificar que satisfacen la ecuación de continuidad y las condiciones de acoplamiento en los vértices del grafo.
   • Los autores abordan la falta de cotas de curvatura fuertes (que típicamente facilitan la construcción de secuencias de recuperación en entornos dinámicos) desplazando los requisitos de regularidad hacia la medida de referencia en lugar de hacia la propia secuencia de recuperación.

3. Existencia y Generalización:

   • La existencia de soluciones para el problema dinámico se demuestra mediante el método directo del Cálculo de Variaciones. Para manejar casos donde las densidades pueden no estar acotadas, los autores extienden el funcional para incluir las partes singulares de las medidas y utilizan un funcional de acción modificado basado en la función convexa $\Psi(u, v) = |v|^2/2u$.
   • La compacidad del conjunto de soluciones se establece mediante estimaciones de Hölder y propiedades del núcleo de calor en grafos métricos.

4. Investigación Numérica:

   • Se propone un algoritmo primal-dual basado en formulaciones de Lagrangiano aumentado para resolver el problema dinámico. El método emplea una formulación mixta que introduce una variable auxiliar para el gradiente de la densidad.
   • La discretización se realiza mediante cocientes de diferencias centrales para los puntos interiores y diferencias de un solo lado en los vértices, preservando la estructura de simetría temporal.

Contribuciones Clave y Resultados

• Equivalencia Estática y Convergencia: El artículo demuestra que el problema estático de Schrödinger en grafos métricos es equivalente a un problema de transporte óptimo entrópico con una función de costo modificada. Establece la $\Gamma$-convergencia de este problema entrópico estático al problema de Kantorovich estático (distancia de Wasserstein al cuadrado) cuando $\beta \to 0$.
• Equivalencia Dinámica: Los autores derivan rigurosamente una formulación dinámica de tipo Benamou-Brenier para el problema de Schrödinger en grafos métricos. Demuestran que, para datos iniciales y finales acotados, las formulaciones estática y dinámica son equivalentes, y que el problema dinámico admite un minimizador único.
• Convergencia de Mínimos: Como consecuencia de la equivalencia estática-dinámica y de la $\Gamma$-convergencia estática, los valores mínimos del problema dinámico de Schrödinger convergen a la distancia de Wasserstein al cuadrado, y los minimizadores convergen a geodésicas de Wasserstein.
• Manejo de Déficits de Curvatura: El artículo destaca una distincción teórica significativa: mientras que la $\Gamma$-convergencia se cumple para las formulaciones estáticas, los autores señalan que la $\Gamma$-convergencia para las formulaciones dinámicas no se sigue directamente en grafos métricos debido a la falta de fuertes cotas de curvatura inferior (específicamente, la ausencia de una estimación de contractividad de Wasserstein nítida sin una constante adicional). La convergencia de los mínimos se logra, en cambio, mediante la equivalencia estática-dinámica.
• Existencia para Datos Generales: La formulación dinámica se extiende a datos iniciales y finales más generales (eliminando el supuesto de densidad acotada), y se demuestra la existencia de soluciones en este contexto más amplio.
• Validación Numérica: Experimentos numéricos en grafos con forma de estrella demuestran la convergencia de las soluciones de Schrödinger dinámicas hacia la geodésica de Wasserstein a medida que $\beta \to 0$. Los resultados ilustran el efecto de suavizado del término de información de Fisher (continuidad en los vértices) para $\beta > 0$, contrastando con la naturaleza discontinua de las geodésicas de transporte óptimo estándar en grafos.

Significación
El artículo proporciona la primera imagen exhaustiva del problema de Schrödinger en grafos métricos, llenando un vacío en la literatura donde tales resultados estaban limitados previamente a los espacios $RCD^*(K, N)$. Demuestra que, a pesar del fallo de las condiciones de curvatura sintética fuerte, las relaciones fundamentales entre los problemas de Schrödinger estático y dinámico y su conexión con el transporte óptimo pueden establecerse rigurosamente. El trabajo clarifica el papel específico de la regularidad del núcleo de calor para compensar la falta de cotas de curvatura y ofrece un marco numérico robusto para resolver estos problemas en espacios métricos no suaves.