Resonance-based integrators for stochastic Schrödinger equations. Convergence and long-time error bounds

Resumen Técnico: Integradores Basados en Resonancia para Ecuaciones de Schrödinger Estocásticas

Planteamiento del Problema
Este trabajo aborda la aproximación numérica de la ecuación de Schrödinger estocástica (SSE) impulsada por ruido de Wiener $Q$ aditivo. El estudio se centra en dos casos específicos: la ecuación lineal con un potencial irregular ($F(\xi) = V(x)\xi$) y la ecuación no lineal cúbica ($F(\xi) = |\xi|^2\xi$). Un desafío principal en este dominio es la baja regularidad temporal de las soluciones inducida por las convoluciones estocásticas, lo cual típicamente limita el orden de convergencia de los esquemas numéricos estándar. Además, los análisis teóricos existentes para tales ecuaciones están en gran medida restringidos a escalas de tiempo cortas. El artículo tiene como objetivo desarrollar y analizar integradores basados en resonancia de baja regularidad que logren órdenes de convergencia más altos bajo suposiciones reducidas de regularidad espacial y, crucialmente, establezcan cotas de error uniformes a largo plazo en regímenes de potenciales y ruido pequeños.

Metodología
La metodología propuesta se basa en el análisis de Fourier y la construcción de integradores numéricos basados en resonancia. La estrategia central implica:

1. Variables Enrolladas: Transformar la solución $\xi(t)$ en una variable enrollada $v(t) = e^{-it\Delta}\xi(t)$ para aislar los términos oscilatorios dominantes.
2. Integración Exacta de Términos Resonantes: Los integradores se construyen integrando exactamente los términos dominantes que surgen en la expansión de Fourier de las interacciones no lineales o del potencial. Esto implica aproximar integrales temporales de funciones oscilatorias (por ejemplo, $\int_0^\tau e^{is(\ell_1^2 + 2\ell_1\ell_2)} ds$) reteniendo la integral exacta para los casos resonantes (donde la fase se anula) y aproximando los casos no resonantes.
3. Técnica de Compensación de Oscilación por Regularidad (RCO): Se adapta una extensión novedosa de la técnica RCO, previamente utilizada en contextos deterministas, al marco estocástico. Esta técnica es esencial para superar la pérdida de regularidad temporal (típicamente $O(\sqrt{\tau})$) causada por las convoluciones estocásticas. Al explotar la cancelación de oscilaciones en términos no resonantes y utilizar un parámetro de corte para frecuencias altas, los autores derivan estimaciones de error mejoradas.
4. Estrategia de Truncamiento: Para el caso no lineal cúbico, donde la convergencia fuerte se ve obstaculizada por la falta de integrabilidad exponencial, los autores emplean una estrategia basada en truncamiento combinada con técnicas de baja regularidad para derivar resultados de convergencia casi segura.

Contribuciones y Resultados Clave

• Ecuación Lineal con Potencial Irregular:

  • Convergencia: Se construye un integrador de baja regularidad resonante (SLR1). Logra una convergencia fuerte de orden $\gamma \in (0, 1]$ en $L^2(\Omega; H^\sigma)$ para soluciones en $H^{\sigma+\gamma}$. Notablemente, para una convergencia de primer orden en $H^\sigma$, el método requiere solo regularidad $H^{\sigma+1}$ para los datos iniciales y el potencial, mejorando el requisito clásico de $H^{\sigma+2}$.
  • Cotas a Largo Plazo: En un régimen con potencial $O(\varepsilon^2)$ y ruido $O(\varepsilon)$, se demuestra que la solución exacta tiene cotas de momentos uniformes hasta tiempos de $O(\varepsilon^{-2})$. Se diseña un integrador no resonante (SNRLR1) para este régimen, que exhibe cotas de error uniformes a largo plazo de tamaño $O(\varepsilon^2 \tau)$ en $L^2(\Omega; L^2)$, una mejora significativa sobre los métodos clásicos.

• Ecuación No Lineal Cúbica:

  • Convergencia: Se introducen tanto integradores resonantes (SLR2) como no resonantes (SNRLR2). Se derivan resultados de convergencia casi segura a baja regularidad, logrando orden $\gamma \in (0, 1)$ en $H^\sigma$ para soluciones en $H^{\sigma+\gamma}$.
  • Cotas a Largo Plazo: En el régimen débilmente no lineal (no linealidad $O(\varepsilon^2)$, ruido $O(\varepsilon)$), se demuestra que el integrador no resonante (SNRLR2) satisface cotas de error uniformes a largo plazo casi seguras de tamaño $O(\varepsilon^2 \tau^\delta)$ para cualquier $\delta < 1$ hasta tiempos de $O(\varepsilon^{-2})$.

• Innovación Teórica:

  • El artículo presenta la primera extensión de la técnica RCO a ecuaciones dispersivas estocásticas. Esta extensión permite a los autores recuperar estimaciones de error mejoradas que superan la limitación estándar de $\sqrt{\tau}$ impuesta por el ruido estocástico, reduciendo los errores a largo plazo en un factor de $\varepsilon^2$ en comparación con otros integradores de baja regularidad.

Significado y Afirmaciones
Los autores afirman que este trabajo representa el primer establecimiento de cotas de error a largo plazo para integradores de baja regularidad aplicados a ecuaciones dispersivas estocásticas. El significado radica en:

1. Requisitos de Regularidad Reducidos: Los métodos requieren menos regularidad espacial ($H^{\sigma+1}$ frente a $H^{\sigma+2}$) para lograr una convergencia de primer orden, haciéndolos aplicables a potenciales y datos iniciales más irregulares.
2. Estabilidad a Largo Plazo: La derivación de cotas de error uniformes en escalas de tiempo de $O(\varepsilon^{-2})` aborda una brecha en la literatura, donde los resultados anteriores estaban confinados a escalas de tiempo cortas.
3. Extensión Estocástica de RCO: La adaptación exitosa de la técnica RCO al entorno estocástico proporciona una nueva herramienta analítica para manejar la irregularidad temporal inherente a las convoluciones estocásticas.

Verificación Numérica
Los resultados teóricos se confirman mediante experimentos numéricos. Estos experimentos demuestran:

• Tasas de convergencia fuerte y casi segura de orden uno para el caso lineal bajo regularidad $H^1$.
• La acotación de los errores para el integrador no resonante (SNRLR1) en intervalos de tiempo largos ($O(\varepsilon^{-2})$), en contraste con el crecimiento del error observado en los métodos resonantes (SLR1) y exponenciales clásicos.
• La reducción de $O(\varepsilon^2)$ en la magnitud del error para el esquema no resonante en el régimen de potencial/ruido pequeño.
• Un comportamiento de convergencia y a largo plazo similar para el caso no lineal cúbico utilizando los integradores propuestos.