Mapas de transporte óptimo inestables

Información Básica

• ID del Artículo: 2510.13265
• Título: Mapas de transporte óptimo inestables
• Autor: Cyril Letrouit (Université Paris-Saclay, CNRS)
• Clasificación: math.OC (Optimización Matemática y Control)
• Fecha de Publicación: 15 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2510.13265

Resumen

La estabilidad de los mapas de transporte óptimo respecto a perturbaciones de las distribuciones marginales es un problema importante que involucra la validez del marco de transporte óptimo linealizado, análisis numérico y estadística. Bajo diversas hipótesis sobre la medida fuente, se sabe que la variación de los mapas de transporte óptimo respecto a cambios en la medida objetivo es estable. Este artículo se enfoca en el mecanismo opuesto—los factores que conducen a la inestabilidad. Los autores identifican dos mecanismos e ilustran la inestabilidad e incluso alta inestabilidad de los mapas de transporte óptimo mediante ejemplos de medidas fuente absolutamente continuas $\rho$ en $\mathbb{R}^d$. Primero se demuestra que la no acotación de la densidad puede conducir a inestabilidad: se construye una densidad fuente en la bola unitaria de $\mathbb{R}^d$ que explota de manera superpolinomial en dos puntos de la frontera, y el mapa de transporte óptimo correspondiente es altamente inestable. Luego se demuestra que incluso para densidades uniformes en conjuntos abiertos acotados, el mapa de transporte óptimo puede ser bastante inestable cuando se aproxima a configuraciones donde se pierde la unicidad del plan óptimo.

Antecedentes y Motivación de la Investigación

Descripción del Problema

Un problema central en la teoría del transporte óptimo es estudiar la estabilidad del mapa de transporte óptimo $T_\mu$ (desde una medida fuente fija $\rho$ a una medida objetivo $\mu$) respecto a perturbaciones de la medida objetivo. Específicamente, se investiga si existen constantes $C, \alpha > 0$ y $p \geq 1$ tales que para cualesquiera medidas de probabilidad $\mu, \nu$ en un conjunto compacto $Y \subset \mathbb{R}^d$, se cumpla: $\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha$

Importancia

La importancia de este problema se manifiesta en múltiples aspectos:

1. Validez del marco de transporte óptimo linealizado: El mapa $\mu \mapsto T_\mu$ proporciona una inmersión del espacio de Wasserstein $(P(Y), W_2)$ en el espacio de Hilbert $L^2(\rho, \mathbb{R}^d)$
2. Análisis numérico: Cuando la medida objetivo $\mu$ solo se puede obtener mediante una aproximación $\hat{\mu}$, es necesario garantizar que el $T_{\hat{\mu}}$ calculado esté cerca del verdadero $T_\mu$
3. Estimación estadística: En transporte óptimo estadístico, se necesita estimar el mapa de transporte óptimo a partir de muestras

Limitaciones de la Investigación Existente

Aunque se han establecido desigualdades de estabilidad cuantificadas bajo diversas hipótesis sobre la medida fuente $\rho$, nunca se ha demostrado que los mapas de transporte óptimo puedan ser inestables. Los resultados existentes se enfocaban principalmente en resultados positivos de estabilidad, careciendo de una comprensión profunda de los mecanismos de inestabilidad.

Contribuciones Principales

1. Primera construcción de mapas de transporte óptimo altamente inestables: Se demuestra la existencia de medidas fuente absolutamente continuas para las cuales ninguna forma de desigualdad de estabilidad cuantificada se cumple
2. Identificación de dos mecanismos de inestabilidad:
   • Explosión superpolinomial de la densidad
   • Aproximación a configuraciones donde se pierde la unicidad del plan óptimo
3. Provisión de construcciones concretas de contraejemplos: Mediante construcciones matemáticas precisas se exhiben fenómenos de inestabilidad
4. Revelación de limitaciones de regularidad del mapa $\mu \mapsto T_\mu$: Se demuestra que en ciertos casos este mapa no es $\frac{1}{2}$-Hölder continuo

Explicación Detallada de Métodos

Definición de la Tarea

Dados $\rho, \mu \in P_2(\mathbb{R}^d)$, donde $\rho$ es absolutamente continua, el teorema de Brenier garantiza la existencia de un único mapa de transporte óptimo $T_\mu \in L^2(\rho)$ desde $\rho$ a $\mu$. Se estudia la estabilidad de $T_\mu$ respecto a perturbaciones de $\mu$.

Método de Construcción del Teorema 1.1

Construcción de la Medida Fuente

Sea $A = (1, 0, \ldots, 0)$, $A' = (-1, 0, \ldots, 0)$, $E = \{A, A'\}$. Se define la función: $f(r) = r^{-d} \min(1, (\log r)^{-2})$

Se construye la densidad: $\rho(x) = c_0 f(\text{dist}(x, E))$ donde $c_0 > 0$ hace que $\rho$ sea una densidad de probabilidad en la bola unitaria $B_{\mathbb{R}^d}(0,1)$.

Secuencia de Medidas Objetivo

Se consideran $B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0)$ y $B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0)$, definiéndose: $\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})$

Estimaciones Clave

Mediante análisis geométrico se demuestra: $\|\|T_{\mu_\theta} - T_{\mu_0}\|\|^2_{L^2(\rho)} \geq \frac{c_0 c_1 \sigma_{d-1} R^2}{|\log(\theta/4)|}$

Mientras que $W_p(\mu_0, \mu_\theta) \sim R\theta/2$, por lo tanto la desigualdad de estabilidad falla para cualquier $\alpha > 0$.

Método de Construcción del Teorema 1.2

Construcción del Conjunto Soporte

Se construye la unión de infinitos pares de paralelepípedos rectangulares: $X = \bigcup_{i=1}^{+\infty} S_i$ donde $S_i = T^+(A_i^+, \ell_i, r_i) \cup T^-(A_i^-, \ell_i, r_i)$.

Selección de Parámetros

Se eligen secuencias satisfaciendo:

• $r_i = c_0 2^{-i}$ (decaimiento superpolinomial)
• $\ell_i = w_i = c_0 c_1 i^{-2}$
• Distancia entre celdas mucho mayor que la distancia dentro de celdas

Demostración de Inestabilidad

Mediante análisis de cambios en el mapa de transporte óptimo bajo perturbaciones específicas, se demuestra: $\frac{\|\|T_\mu - T_{\nu_i}\|\|^2_{L^2(\rho)}}{W_p(\mu, \nu_i)^{2\alpha}} \geq C w_i^{3-\frac{2\alpha}{p}} r_i^{1-\frac{2\alpha}{p}-2\alpha}$

Cuando $\alpha > \frac{p}{2(p+1)}$, el lado derecho tiende a infinito.

Configuración Experimental

Verificación Teórica

Este es un trabajo puramente teórico, verificando resultados mediante demostraciones matemáticas rigurosas:

1. Análisis geométrico: Análisis detallado de la estructura geométrica del mapa de transporte óptimo
2. Estimaciones de medidas: Cálculo exacto de medidas en varias regiones
3. Cálculo de distancias: Estimaciones precisas de distancias de Wasserstein y distancias $L^2$

Verificación de Construcciones

Se verifica la corrección de las construcciones mediante los siguientes pasos:

1. Demostración de que la densidad construida es efectivamente una densidad de probabilidad
2. Verificación de expresiones explícitas del mapa de transporte óptimo
3. Confirmación de que la configuración geométrica satisface las condiciones de separación requeridas

Resultados Experimentales

Resultados Principales

Resultados del Teorema 1.1

Para la densidad fuente construida $\rho$ (que explota de manera superpolinomial en dos puntos de la frontera de la bola unitaria), se demuestra que para cualquier bola $Y = B_{\mathbb{R}^d}(0,R)$, cualesquiera $C, \alpha > 0$ y $p \geq 1$, la desigualdad: $\forall \mu, \nu \in P(Y), \|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha$ no se cumple.

Resultados del Teorema 1.2

Para la medida de probabilidad uniforme $\rho$ en un conjunto abierto acotado, se demuestra que para cualesquiera $C > 0$, $p \geq 1$ y $\alpha > \frac{p}{2(p+1)}$ (en particular, $\alpha = \frac{1}{2}$), la desigualdad de estabilidad falla.

Hallazgos Técnicos

1. Carácter crítico de la explosión de densidad: La explosión polinomial $r^{-d+\delta}$ mantiene estabilidad, pero la explosión superpolinomial conduce a inestabilidad completa
2. Limitaciones de regularidad Hölder: El mapa $\mu \mapsto T_\mu$ en ciertos casos no es $\frac{1}{2}$-Hölder continuo, e incluso no es $\frac{1}{3}$-Hölder continuo
3. Presagio de no unicidad: La pérdida de unicidad del plan óptimo es precedida por deterioro de estabilidad

Trabajo Relacionado

Resultados Positivos de Estabilidad

• Gigli (2011): Teoría de límites superiores de $\frac{1}{2}$-continuidad Hölder del mapa $\mu \mapsto T_\mu$
• Letrouit & Mérigot (2024): Desigualdades de estabilidad para densidades acotadas en dominios de John
• Otros trabajos: Resultados de estabilidad para densidades log-cóncavas, densidades con explosión polinomial, etc.

Investigación de Inestabilidad

Anteriormente solo existían ejemplos de inestabilidad del potencial de Kantorovich; este artículo construye por primera vez mapas de transporte inestables.

Conclusiones y Discusión

Conclusiones Principales

1. Los mapas de transporte óptimo pueden ser altamente inestables, y cualquier desigualdad de estabilidad cuantificada puede fallar
2. La explosión superpolinomial de la densidad es una condición suficiente para inestabilidad
3. Incluso para densidades uniformes, la inestabilidad aparece cuando se aproxima a configuraciones de no unicidad

Limitaciones

1. El Teorema 1.1 requiere densidad no acotada, lo cual puede no ser realista en aplicaciones prácticas
2. La construcción del Teorema 1.2 requiere infinitas componentes conexas, con estructura geométrica compleja
3. No se ha determinado si densidades acotadas pueden lograr la fuerte inestabilidad del Teorema 1.1

Direcciones Futuras

1. Investigación de la inestabilidad más fuerte para el caso de densidades acotadas
2. Desarrollo de métodos de estimación estadística para medidas fuente inestables
3. Exploración de conexiones profundas entre inestabilidad y estructura geométrica

Evaluación Profunda

Fortalezas

1. Avance teórico: Primera demostración de inestabilidad de mapas de transporte óptimo, llenando un vacío teórico
2. Construcciones ingeniosas: Mediante construcciones geométricas y selección de parámetros astutos se logran objetivos teóricos
3. Análisis riguroso: Demostraciones matemáticas completas con manejo adecuado de detalles técnicos
4. Perspectivas profundas: Revelación de conexiones intrínsecas entre estabilidad, unicidad y regularidad de densidad

Deficiencias

1. Limitaciones de aplicabilidad práctica: Los ejemplos construidos pueden no ser comunes en aplicaciones reales
2. Problemas abiertos: Ciertos problemas naturales (como la fuerte inestabilidad para densidades acotadas) permanecen sin resolver
3. Ausencia de verificación numérica: Como trabajo teórico, carece de experimentos numéricos de verificación

Impacto

1. Valor teórico: Proporciona contraejemplos importantes e perspectivas para la teoría del transporte óptimo
2. Advertencia para aplicaciones: Alerta a los practicantes sobre riesgos de inestabilidad en ciertas configuraciones
3. Inspiración para investigación: Proporciona nuevas direcciones y herramientas para investigación posterior

Escenarios Aplicables

1. Investigación profunda de la teoría del transporte óptimo
2. Análisis teórico del marco de transporte óptimo linealizado
3. Análisis de estabilidad en transporte óptimo estadístico
4. Evaluación de estabilidad de algoritmos numéricos de transporte óptimo

Referencias

El artículo cita literatura importante en el campo del transporte óptimo, incluyendo trabajos fundacionales de Brenier, resultados de regularidad de Gigli, e investigaciones recientes sobre estabilidad cuantificada, proporcionando una base teórica sólida para la investigación.