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

Resumo Técnico: Integradores Baseados em Ressonância para Equações de Schrödinger Estocásticas

Enunciado do Problema
Este trabalho aborda a aproximação numérica da equação de Schrödinger estocástica (SSE) impulsionada por ruído de Wiener $Q$-aditivo. O estudo foca em dois casos específicos: a equação linear com potencial irregular ($F(\xi) = V(x)\xi$) e a equação não linear cúbica ($F(\xi) = |\xi|^2\xi$). Um desafio primário neste domínio é a baixa regularidade temporal das soluções induzida por convoluções estocásticas, o que tipicamente limita a ordem de convergência de esquemas numéricos padrão. Além disso, análises teóricas existentes para tais equações são amplamente restritas a escalas de tempo curtas. O artigo visa desenvolver e analisar integradores de baixa regularidade baseados em ressonância que alcancem ordens de convergência mais altas sob suposições reduzidas de regularidade espacial e, crucialmente, estabeleçam limites de erro uniformes de longo prazo em regimes de potenciais e ruídos pequenos.

Metodologia
A metodologia proposta baseia-se em análise de Fourier e na construção de integradores numéricos baseados em ressonância. A estratégia central envolve:

1. Variáveis Torcidas: Transformar a solução $\xi(t)$ em uma variável torcida $v(t) = e^{-it\Delta}\xi(t)$ para isolar os termos oscilatórios dominantes.
2. Integração Exata de Termos Ressonantes: Os integradores são construídos integrando exatamente os termos dominantes que surgem na expansão de Fourier das interações não lineares ou de potencial. Isso envolve aproximar integrais temporais de funções oscilatórias (por exemplo, $\int_0^\tau e^{is(\ell_1^2 + 2\ell_1\ell_2)} ds$) mantendo a integral exata para casos ressonantes (onde a fase se anula) e aproximando casos não ressonantes.
3. Técnica de Compensação de Oscilação por Regularidade (RCO): Uma extensão novel da técnica RCO, anteriormente utilizada em contextos determinísticos, é adaptada ao quadro estocástico. Esta técnica é essencial para superar a perda de regularidade temporal (tipicamente $O(\sqrt{\tau})$) causada por convoluções estocásticas. Ao explorar a cancelação de oscilações em termos não ressonantes e utilizar um parâmetro de corte para altas frequências, os autores derivam estimativas de erro melhoradas.
4. Estratégia de Truncamento: Para o caso não linear cúbico, onde a convergência forte é impedida pela falta de integrabilidade exponencial, os autores empregam uma estratégia baseada em truncamento combinada com técnicas de baixa regularidade para derivar resultados de convergência pontual (pathwise).

Principais Contribuições e Resultados

• Equação Linear com Potencial Irregular:

  • Convergência: Um integrador de baixa regularidade ressonante (SLR1) é construído. Ele alcança convergência forte de ordem $\gamma \in (0, 1]$ em $L^2(\Omega; H^\sigma)$ para soluções em $H^{\sigma+\gamma}. Notavelmente, para convergência de primeira ordem em$H^\sigma$, o método requer apenas regularidade $H^{\sigma+1}$ para os dados iniciais e o potencial, melhorando sobre o requisito clássico $H^{\sigma+2}$.
  • Limites de Longo Prazo: Em um regime com potencial $O(\varepsilon^2)$ e ruído $O(\varepsilon)$, mostra-se que a solução exata possui limites de momentos uniformes até tempos de $O(\varepsilon^{-2})$. Um integrador não ressonante (SNRLR1) é projetado para este regime, exibindo limites de erro uniformes de longo prazo de tamanho $O(\varepsilon^2 \tau)$ em $L^2(\Omega; L^2)$, uma melhoria significativa sobre métodos clássicos.

• Equação Não Linear Cúbica:

  • Convergência: Tanto integradores ressonantes (SLR2) quanto não ressonantes (SNRLR2) são introduzidos. Resultados de convergência pontual são derivados em baixa regularidade, alcançando ordem $\gamma \in (0, 1)$ em $H^\sigma$ para soluções em $H^{\sigma+\gamma}$.
  • Limites de Longo Prazo: No regime fracamente não linear (não linearidade $O(\varepsilon^2)$, ruído $O(\varepsilon)$), prova-se que o integrador não ressonante (SNRLR2) satisfaz limites de erro pontuais uniformes de longo prazo de tamanho $O(\varepsilon^2 \tau^\delta)$ para qualquer $\delta < 1$ até tempos de $O(\varepsilon^{-2})$.

• Inovação Teórica:

  • O artigo apresenta a primeira extensão da técnica RCO para equações dispersivas estocásticas. Esta extensão permite aos autores recuperar estimativas de erro melhoradas que superam a limitação padrão de $\sqrt{\tau}$ imposta pelo ruído estocástico, reduzindo os erros de longo prazo por um fator de $\varepsilon^2$ em comparação com outros integradores de baixa regularidade.

Significado e Afirmações
Os autores afirmam que este trabalho representa o primeiro estabelecimento de limites de erro de longo prazo para integradores de baixa regularidade aplicados a equações dispersivas estocásticas. O significado reside em:

1. Requisitos Reduzidos de Regularidade: Os métodos exigem menos regularidade espacial ($H^{\sigma+1}$ versus $H^{\sigma+2}$) para alcançar convergência de primeira ordem, tornando-os aplicáveis a potenciais e dados iniciais mais irregulares.
2. Estabilidade de Longo Prazo: A derivação de limites de erro uniformes em escalas de tempo de $O(\varepsilon^{-2})$ aborda uma lacuna na literatura, onde resultados anteriores estavam confinados a escalas de tempo curtas.
3. Extensão Estocástica da RCO: A adaptação bem-sucedida da técnica RCO ao cenário estocástico fornece uma nova ferramenta analítica para lidar com a irregularidade temporal inerente às convoluções estocásticas.

Verificação Numérica
Os resultados teóricos são confirmados através de experimentos numéricos. Estes experimentos demonstram:

• Taxas de convergência forte e pontual de ordem um para o caso linear sob regularidade $H^1$.
• A limitação dos erros para o integrador não ressonante (SNRLR1) ao longo de longos intervalos de tempo ($O(\varepsilon^{-2})$), contrastando com o crescimento de erro observado nos métodos ressonantes (SLR1) e exponenciais clássicos.
• A redução de $O(\varepsilon^2)$ na magnitude do erro para o esquema não ressonante no regime de potencial/ruído pequeno.
• Comportamento de convergência e longo prazo similar para o caso não linear cúbico utilizando os integradores propostos.