Stochastic Homogenization of Non-local Hamilton-Jacobi-Bellman equations

Resumen Técnico: Homogeneización Estocástica de Ecuaciones de Hamilton-Jacobi-Bellman No Locales

Planteamiento del Problema
Este artículo aborda la homogeneización estocástica de ecuaciones de Hamilton-Jacobi-Bellman (HJB) no locales en un medio aleatorio estacionario y ergódico. La ecuación específica bajo estudio es:
$$\begin{cases} \partial_t u^\epsilon(t, x) = \epsilon^{\alpha-1} \mathcal{L} u^\epsilon(t, x) + H\left(\nabla u^\epsilon(t, x), \frac{x}{\epsilon}, \omega\right), & (t, x) \in [0, \infty) \times \mathbb{R}^n, \\ u^\epsilon(0, x) = u_0(x), & x \in \mathbb{R}^n, \end{cases}$$
donde $\alpha \in (1, 2)$, $\omega$ representa un elemento de un espacio de probabilidad $(\Omega, \mathcal{F}, P)$, y $H$ es un hamiltoniano convexo con crecimiento superlineal con oscilaciones aleatorias rápidas. El operador $\mathcal{L}$ es un operador integro-diferencial no local definido por una medida de Lévy de la forma $K(\hat{y})|y|^{-n-\alpha}dy$, donde $K$ es una función suave, no negativa, par y uniformemente no degenerada en la esfera unidad. Este operador actúa como el generador infinitesimal de un proceso de difusión de saltos de tipo $\alpha$-estable.

El objetivo es demostrar que, a medida que $\epsilon \to 0$, $u^\epsilon$ converge uniformemente de forma local (casi seguramente) a la solución $u$ de una ecuación homogeneizada determinista:
$$\begin{cases} \partial_t u(t, x) = \bar{H}(\nabla u(t, x)), \\ u(0, x) = u_0(x), \end{cases}$$
donde $\bar{H}$ es un hamiltoniano efectivo determinista, convexo e independiente de la variable espacial $x$ y del parámetro aleatorio $\omega$.

Metodología
Los autores adaptan el método desarrollado por Kosygina, Rezakhanlou y Varadhan (CPAM 2006) para ecuaciones HJB locales con términos de Laplaciano evanescentes al entorno no local que involucra difusiones de saltos. La metodología se basa en tres pilares principales:

1. Representación de Control Óptimo Estocástico: La solución $u^\epsilon$ se representa mediante un problema de control estocástico impulsado por un proceso de difusión de saltos. Los autores establecen que el operador no local $\mathcal{L}$ puede reescribirse en "forma de divergencia" $\mathcal{L}u = \nabla \cdot (\mathcal{R}\nabla u)$, donde $\mathcal{R}$ es un operador integral regular. Esta caracterización estructural es crucial para manejar los términos no locales en el proceso de homogeneización.

2. Proceso Ambiental y Hamiltoniano Efectivo: Siguiendo el enfoque de [24], los autores elevan los procesos de difusión de saltos desde el espacio físico $\mathbb{R}^n$ al espacio de probabilidad $\Omega$ para definir un "proceso ambiental" $\theta_t$. Utilizando la ergodicidad de este proceso, definen el Lagrangiano efectivo $\bar{L}$ y, posteriormente, el Hamiltoniano efectivo $\bar{H}$ mediante una fórmula variacional que involucra medidas invariantes del proceso ambiental.

3. Super-correctores Aproximados: Para establecer la cota superior de la convergencia (la dirección difícil), los autores construyen super-soluciones aproximadas (super-correctores). Esto implica encontrar un campo vectorial aleatorio $v$ tal que una desigualdad específica que involucra el operador de divergencia no local $J$ y el hamiltoniano se cumpla:
   $$Jv(\omega) + H(p + v(\omega), \omega) \leq \bar{H}(p).$$
   La construcción de estos correctores es la principal novedad técnica, lo que requiere el manejo de la singularidad del operador no local y la falta de regularidad elíptica estándar.

Contribuciones Clave y Resultados

• Teorema Principal (Teorema 1.1): Bajo supuestos estándar de estacionariedad y ergodicidad, así como condiciones de convexidad/superlinealidad (Supuesto A), el artículo demuestra la convergencia casi segura uniformemente local de $u^\epsilon$ a $u$. El resultado se cumple bajo cualquiera de dos supuestos técnicos adicionales:

  • Supuesto (B): Una condición específica de convexidad-continuidad sobre el hamiltoniano $H$ (relacionada con la separación de variables).
  • Supuesto (C): Una condición que requiere que el exponente de crecimiento $\beta$ del hamiltoniano satisfaga $\beta > n$, lo que permite la aplicación de los teoremas de incrustación de Sobolev.

• Caracterización en Forma de Divergencia: El artículo caracteriza explícitamente el operador no local $\mathcal{L}$ como la divergencia de un operador integral regular que actúa sobre el gradiente (Sección 2.2). Esta representación, $\mathcal{L} = J \circ \nabla = \nabla \cdot \mathcal{R} \circ \nabla$, es esencial para definir los correctores en el entorno aleatorio.

• Estimaciones de Regularidad: Los autores derivan estimaciones de regularidad explícitas para la solución $u^\epsilon$ con respecto al espacio y al tiempo (Teorema 4.1). Estas estimaciones se utilizan para elevar la convergencia puntual a la convergencia uniforme local.

• Construcción de Super-correctores: El artículo detalla la construcción de super-correctores aproximados para el caso no local. A diferencia del caso local, donde los correctores suelen construirse mediante métodos de EDP, aquí los autores utilizan el análisis convexo en el espacio de variables aleatorias y técnicas de suavizado (mollification) para manejar el operador no local $J$.

Significado y Pretensiones
El artículo afirma ser el primer trabajo que aborda la homogeneización estocástica de ecuaciones HJB no locales con términos no locales evanescentes en un entorno aleatorio estacionario.

• Extensión de la Teoría Existente: El trabajo extiende los resultados seminales de Kosygina, Rezakhanlou y Varadhan [24] de entornos locales (difusión) a entornos no locales (difusión de saltos).
• Adaptación Metodológica: Demuestra que las técnicas de representación de control y de proceso ambiental, diseñadas originalmente para Laplacianos evanescentes, pueden adaptarse con éxito a los procesos de saltos. Los autores señalan que, si bien la teoría de soluciones viscosas para ecuaciones HJB no locales ha sido estudiada recientemente (por ejemplo, [6, 17, 13]), su enfoque se basa principalmente en representaciones de control estocástico y procesos ambientales abstractos en lugar de una teoría profunda de soluciones viscosas.
• Alcance y Limitaciones: Los autores mantienen un alcance modesto, centrándose en la homogeneización cualitativa. Establecen explícitamente que no abordan la tasa de convergencia (que se sabe es difícil incluso en el caso local) ni el caso de entornos aleatorios dinámicos (coeficientes dependientes del tiempo), aunque sugieren estos como direcciones naturales para trabajos futuros. El resultado está restringido al caso en el que el término no local ($\epsilon^{\alpha-1} \mathcal{L}$) desaparece cuando $\epsilon \to 0$.

En resumen, el artículo establece una base rigurosa para la homogeneización de ecuaciones HJB no locales impulsadas por difusiones de saltos, tendiendo un puente entre la teoría del control estocástico para procesos de salto y la teoría de la homogeneización en medios aleatorios.