Strong convergence of finite element schemes for the stochastic Landau--Lifshitz--Bloch equation

1. Problema e Contexto

O artigo aborda a dinâmica de magnetização em ferromagnetos a altas temperaturas (acima da temperatura de Curie). Nesse regime, as flutuações térmicas e transições induzidas por ruído tornam-se significativas, exigindo uma modelagem estocástica. O modelo matemático central é a Equação Estocástica de Landau–Lifshitz–Bloch (sLLB), uma equação diferencial parcial estocástica (EDPE) vetorial e quase-linear.

A equação (1.1) no domínio $D \subset \mathbb{R}^d$ ($d=1, 2$) é dada por:
$$du = \left( \kappa_1\Delta u + \gamma u \times \Delta u - \kappa_2(1 + \mu |u|^2)u \right) dt + (\kappa_1 g + \gamma u \times g) \circ dW$$
onde $u$ é o vetor de magnetização, $W$ é um processo de Wiener, e os termos não lineares (produto vetorial e termo cúbico) apresentam desafios significativos para a análise numérica, pois não são globalmente Lipschitzianos nem monotônicos.

O objetivo principal é desenvolver e analisar esquemas de Elementos Finitos (EF) totalmente discretos (semi-implícitos e implícitos) que garantam a convergência forte em $L^2(\Omega)$ para a solução exata, com taxas de convergência explícitas, superando limitações de trabalhos anteriores que apenas garantiam convergência em probabilidade ou convergência forte apenas até tempos de parada discretos.

2. Metodologia

A análise matemática e numérica baseia-se em várias técnicas avançadas:

• Estimativas de Momentos Exponenciais: Um dos pilares do trabalho é a derivação de novas limites de momentos exponenciais para a solução exata (Teorema 3.5). Essas estimativas são cruciais para controlar o crescimento da solução e permitir a análise de erro em todo o espaço de amostra, não apenas em subconjuntos de probabilidade alta.
• Técnica de Localização (Localization): Para lidar com a não-linearidade não globalmente Lipschitziana, o autor utiliza a técnica de localização. Define-se um conjunto de "bons" eventos $\Omega_{\kappa, m}$ onde a norma $H^1$ da solução exata é limitada por uma constante $\kappa$. A análise de erro é realizada condicionada a esses eventos.
• Decomposição de Erro: O erro numérico é decomposto em duas partes: o erro de aproximação do operador de projeção (erro de interpolação) e o erro discreto da solução numérica.
• Esquemas Numéricos Propostos:
  • Para $d=1$: Um esquema semi-implícito (baseado em elementos $C^0$) e um novo esquema totalmente implícito. O esquema implícito é projetado para obter melhores propriedades de estabilidade.
  • Para $d=2$: Um esquema semi-implícito utilizando elementos finitos conformes $C^1$ (necessário devido à regularidade exigida pela regularização da equação, Eq. 1.2, que introduz um termo de quarta ordem $\Delta^2$).
• Estabilidade e Unicidade: O trabalho prova a estabilidade exponencial em média quadrática e a unicidade da medida invariante para o caso unidimensional sob a suposição de ruído suficientemente pequeno.

3. Principais Contribuições e Resultados

A. Limites de Momentos Exponenciais (Teorema 3.5)

O autor estabelece limites exponenciais para as normas $L^2$ e $H^1$ da solução exata. Isso é um resultado novo para a equação sLLB e é fundamental para provar a convergência forte global. Como corolário, demonstra-se a estabilidade exponencial em média quadrática e a unicidade da medida invariante para $d=1$ com ruído pequeno.

B. Taxas de Convergência Aprimoradas

O trabalho refina os resultados existentes (como os de [14]) fornecendo taxas de convergência mais agudas:

• Esquema Semi-implícito ($d=1$): Convergência forte em $L^2(\Omega)$ com taxa $O(h + k^{1/2-\delta})$ e convergência em probabilidade com taxa $O(h + k^{1/2-\delta})$.
• Esquema Semi-implícito ($d=2$, regularizado): Convergência forte com taxa $O(h^2 + k^{1/2-\delta})$ e em probabilidade com a mesma taxa.
• Esquema Implícito ($d=1$): O autor propõe um novo esquema implícito que atinge uma taxa de convergência ótima em probabilidade na norma $H^1$, especificamente $O(h + k^{1/2-\delta})$.

C. Convergência Forte Global

Diferente de trabalhos anteriores que limitavam a convergência forte a tempos de parada, este artigo prova a convergência forte em $L^2(\Omega)$ sobre todo o espaço de probabilidade para todos os esquemas analisados. As taxas de convergência forte são subótimas (dependem de logaritmos ou raízes de logaritmos), mas garantem que o erro médio quadrático tende a zero globalmente.

D. Resultados Numéricos

O artigo apresenta experimentos numéricos realizados no ambiente FEniCS que validam as previsões teóricas:

• Simulação 1 ($d=1$, Esquema Implícito): Confirma a convergência de ordem 1 no espaço (norma $H^1$) e ordem 1/2 no tempo, consistente com a teoria.
• Simulação 2 ($d=2$, Esquema Semi-implícito): Embora a análise teórica rigorosa para $d=2$ no esquema não regularizado seja complexa, os experimentos mostram o comportamento esperado de convergência.
• Simulação 3 (Estabilidade): Valida a Proposição 3.6, mostrando que, para ruído pequeno, a diferença entre duas soluções com condições iniciais distintas decai exponencialmente. Para ruído grande, essa estabilidade é perdida, ilustrando a dependência da estabilidade da magnitude do ruído.

4. Significado e Impacto

Este trabalho é significativo por várias razões:

1. Rigor Matemático: Preenche lacunas na análise de convergência forte para EDPEs quase-lineares com não-linearidades não monotônicas, um problema aberto na literatura numérica de micromagnetismo estocástico.
2. Novas Ferramentas Analíticas: A derivação de limites de momentos exponenciais específicos para a sLLB abre caminho para análises de estabilidade e ergodicidade mais robustas.
3. Aplicabilidade Prática: Ao fornecer taxas de convergência explícitas e esquemas estáveis, o trabalho oferece uma base teórica sólida para simulações computacionais precisas de dispositivos de armazenamento magnético (como HAMR - Gravação Magnética Assistida por Calor) operando em altas temperaturas.
4. Validação Numérica: A confirmação experimental das taxas de convergência e do comportamento de estabilidade reforça a utilidade prática dos métodos propostos.

Em resumo, o artigo avança o estado da arte na simulação numérica de fenômenos magnéticos estocásticos, provando teoremas de convergência forte que eram anteriormente inatingíveis e fornecendo esquemas numéricos eficientes e validados.