An Aronson-Bénilan / Li-Yau estimate in the JKO scheme in small dimension

Resumen Técnico: Estimación de Aronson-Bénilan / Li-Yau en el Esquema JKO en Dimensiones Pequeñas

1. Planteamiento del Problema

El artículo se centra en el análisis del esquema JKO (Jordan-Kinderlehrer-Otto), una aproximación variacional implícita (tipo Euler) para ecuaciones de difusión en el espacio de medidas de probabilidad con segundo momento finito. El objetivo es recuperar propiedades cualitativas y cuantitativas conocidas para el flujo continuo en el nivel discreto, uniformemente en el paso de tiempo $\tau$.

Las ecuaciones estudiadas son:
$$\partial_t \rho = \Delta \rho^m$$
donde $m$ define tres regímenes:

• $m > 1$: Ecuación de medio poroso (difusión lenta, fronteras libres).
• $m = 1$: Ecuación del calor (difusión lineal).
• $0 < m < 1$: Régimen de difusión rápida (soluciones estrictamente positivas).

El problema principal abordado es la derivación de una estimación de Aronson-Bénilan / Li-Yau en el esquema discreto. En el caso continuo, esta estimación establece un límite inferior sub-harmónico para la variable de presión $p$ (definida como $p = \frac{m}{m-1}\rho^{m-1}$ para $m \neq 1$ y $p = \log \rho$ para $m=1$):
$$\Delta p \geq -\frac{\alpha_{d,m}}{t}, \quad \text{donde } \alpha_{d,m} = \frac{d}{d(m-1)+2}$$
Esta desigualdad es fundamental para probar efectos de regularización $L^1-L^\infty$ y estudiar la regularidad de las soluciones y sus fronteras libres.

Brecha en la literatura: Aunque el caso superlineal ($m \geq 1$) está bien entendido, el caso de difusión rápida ($m < 1$) presenta dificultades debido al comportamiento no superlineal de la entropía. Además, la extensión de estas estimaciones al esquema JKO en dominios generales y dimensiones altas ha sido un desafío abierto.

2. Metodología

La estrategia de prueba se basa en un principio de máximo aplicado al determinante del Hessiano de los potenciales de Brenier, iterado paso a paso a lo largo del esquema.

• Variables Discretas: Se define la variable de presión discreta $p^\tau_k$ y una función auxiliar convexa $u^\tau_k = \tau p^\tau_k + \frac{1}{2}|\cdot|^2$.
• Análisis del Hessiano: El núcleo de la demostración es demostrar que el determinante del Hessiano de $u^\tau_k$ satisface una cota inferior que mejora en cada paso del esquema.
• Restricción Dimensional: La metodología se limita a dimensiones $d=1$ y $d=2$. Esto es crucial porque en estas dimensiones es posible controlar el determinante del Hessiano utilizando tanto la traza (Laplaciano) como la traza de la inversa del Hessiano mediante desigualdades algebraicas específicas (relacionadas con la desigualdad AM-GM para matrices). En dimensiones superiores, esta relación algebraica no se sostiene de la misma manera.
• Dominios Simples: Se trabaja en dominios convexos simples: toro, cubos, cuartos de espacio, semiespacios y el espacio completo. La restricción a dominios simples (especialmente cubos) permite manejar los puntos de máximo en la frontera mediante un análisis detallado del comportamiento del mapa de transporte óptimo en las caras y esquinas del dominio.

3. Contribuciones Clave

1. Tratamiento Riguroso del Caso $m < 1$:

   • El artículo llena un vacío en la literatura al probar la existencia, unicidad y condiciones de optimalidad para el problema de un paso del esquema JKO en el régimen de difusión rápida ($m < 1$), incluso en dominios no acotados (bajo condiciones de regularidad de volumen).
   • Se demuestra que los minimizadores son absolutamente continuos y se derivan las condiciones de optimalidad precisas, llenando la brecha teórica sobre el comportamiento de la entropía no superlineal.

2. Estimación de Aronson-Bénilan Discreta:

   • Se demuestra que la función $u^\tau_k$ es convexa y que su determinante del Hessiano satisface:
     $$\det(D^2 u^\tau_k)^{1/d} \geq 1 - X_k$$
     donde $(X_k)_{k \geq 1}$ es una secuencia universal en $[0, 1]$ que depende solo de $m$ y $d$.
   • Se establece el comportamiento asintótico de esta secuencia: $X_k \sim \frac{1}{(d(m-1)+2)k}$.

3. Convergencia al Límite Continuo:

   • Al tomar $\tau \to 0$ y $k\tau \approx t$, la estimación discreta recupera exactamente la constante de la estimación de Aronson-Bénilan continua:
     $$\Delta p \geq -\frac{d}{d(m-1)+2} \frac{1}{t}$$
   • Se demuestra que la estimación discreta es ligeramente más fuerte que la versión continua esperada debido a la relación entre el determinante y el Laplaciano en dimensiones bajas.

4. Acotación Uniforme $L^\infty$:

   • Como consecuencia directa, se obtienen cotas locales uniformes en $L^\infty$ para la densidad $\rho^\tau_k$, uniformes en el paso de tiempo $\tau$. Esto confirma la robustez del esquema JKO para recuperar los efectos de regularización $L^1-L^\infty$ clásicos.

4. Resultados Principales

• Teorema 1.1 (Estimación en JKO): Para $d=1, 2$ y dominios simples, la función $u^\tau_k$ es convexa y su determinante del Hessiano está acotado inferiormente por $1-X_k$. La secuencia $X_k$ decae como $O(1/k)$.
• Corolario 1.2 (Acotación $L^\infty$): Existe una constante $M$ tal que $\|\rho^\tau_k\|_{L^\infty(B_r)} \leq M$ para todo $\tau$ suficientemente pequeño y $t \geq t_0$. Esto coincide con los resultados clásicos de regularización para las ecuaciones de medio poroso, calor y difusión rápida.
• Análisis de Frontera: Se demuestra que, en dominios tipo cubo, el mapa de transporte óptimo preserva la estructura de las caras y esquinas, lo que permite aplicar condiciones de optimalidad en la frontera y asegurar que el gradiente de la función objetivo se anula o cumple condiciones de signo necesarias para el principio de máximo.

5. Significado e Impacto

• Fundamentos Teóricos: Este trabajo proporciona una justificación rigurosa de por qué el esquema JKO preserva las propiedades de regularidad más finas de las ecuaciones de difusión no lineales, incluso en regímenes donde la entropía no es estrictamente convexa ($m < 1$).
• Herramienta para Análisis Numérico: Las cotas uniformes en $L^\infty$ y la estimación de Aronson-Bénilan son herramientas esenciales para probar tasas de convergencia más rápidas y robustez en simulaciones numéricas basadas en el esquema JKO.
• Limitaciones y Futuro: La restricción a dimensiones $d \leq 2$ y dominios simples es una limitación técnica derivada de la necesidad de controlar algebraicamente el determinante del Hessiano. Extender estos resultados a dimensiones superiores o dominios con fronteras más complejas requeriría nuevas ideas, posiblemente relacionadas con la teoría de regularidad de ecuaciones de Monge-Ampère en dominios no lisos.

En resumen, el artículo logra un avance significativo al conectar la teoría de flujos de gradiente en espacios de Wasserstein con estimaciones de regularidad clásicas de EDPs, demostrando que el esquema JKO no solo converge, sino que preserva la estructura geométrica y analítica de la solución continua en el nivel discreto, al menos en dimensiones bajas.