A priori error analysis of a mass-lumped midpoint finite element method with a structure-preserving solver for the Landau--Lifshitz--Gilbert equation with Dzyaloshinskii--Moriya interaction

Resumen Técnico: Análisis de Error A Priori de un Método de Elementos Finitos de Punto Medio con Masa Agrupada para la Ecuación LLG con DMI

1. Planteamiento del Problema
El artículo aborda la aproximación numérica de la ecuación de Landau–Lifshitz–Gilbert (LLG) que gobierna la dinámica de la magnetización en materiales ferromagnéticos a temperaturas por debajo del punto de Curie. Específicamente, el estudio se centra en la inclusión de la interacción Dzyaloshinskii–Moriya (DMI) de volumen, la cual introduce un término de primer orden al campo efectivo e induce una condición de contorno quiral natural. Los principales desafíos matemáticos son la restricción de longitud unitaria puntual ($|m|=1$), la no convexidad del espacio de la solución y la estructura fuertemente no lineal de la ecuación. Si bien el modelo de solo intercambio ha sido ampliamente estudiado, la inclusión de la DMI modifica la estructura variacional y el comportamiento de contorno, lo que hace que la construcción de esquemas de elementos finitos estables, precisos y que preserven la estructura sea particularmente delicada.

2. Metodología
Los autores proponen y analizan un esquema de elementos finitos de punto medio con masa agrupada (mass-lumped) totalmente discreto. La metodología involucra los siguientes componentes:

• Discretización: El dominio espacial se discretiza utilizando elementos finitos de Lagrange conformes de bajo orden ($P_1$). La integración temporal emplea la regla del punto medio implícito.
• Agrupación de Masa (Mass-Lumping): Se utiliza un producto interno reducido de integración (masa agrupada). Esta técnica es crucial, ya que permite que la restricción de longitud unitaria nodal se cumpla exactamente en cada nodo, una propiedad que los productos internos $L^2$ estándar a menudo no logran preservar en esquemas no lineales.
• Operador de Campo Efectivo: El operador de campo efectivo discreto incorpora tanto la rigidez de intercambio como la DMI de volumen. El término DMI se maneja mediante una forma bilineal discreta que respeta naturalmente la condición de contorno quiral. Para facilitar el análisis, se introduce un operador de campo efectivo desplazado para asegurar la coercitividad uniforme de la forma bilineal discreta.
• Estrategia de Solución: Dado que el esquema de punto medio resulta en un sistema algebraico no lineal en cada paso de tiempo, el artículo propone una iteración de punto fijo que preserva la estructura. A diferencia de los resolvedores genéricos, esta iteración está diseñada para preservar la restricción de longitud unitaria nodal en cada iteración.
• Marco de Análisis de Error: El análisis compara la solución numérica con la proyección elíptica de la solución exacta. La demostración se basa en las cancelaciones geométricas inherentes al método de punto medio, la estimación cuidadosa de la contribución de primer orden de la DMI y las estimaciones de estabilidad derivadas de una ley de energía discreta.

3. Contribuciones Clave
El artículo realiza tres contribuciones primarias a la literatura de micromagnetismo computacional:

1. Formulación y Estabilidad: Los autores formulan un esquema de punto medio con masa agrupada totalmente discreto para la ecuación LLG con intercambio y DMI de volumen. Demuestran que este esquema preserva exactamente la restricción de longitud unitaria nodal y satisface una identidad de energía discreta, proporcionando una base para la estabilidad sin requerir restricciones de malla geométricas (como ángulos agudos) o condiciones estrictas de tipo CFL que se encuentran a menudo en otros métodos que preservan la estructura.
2. Estimaciones de Error A Priori Óptimas (Esquema Ideal): Bajo supuestos de regularidad adecuados de la solución exacta, el artículo demuestra una estimación de convergencia de orden óptimo para el esquema ideal (asumiendo la solución exacta del sistema no lineal). El error converge con orden $O(h + k^2)$ en la norma de energía ($H^1$), representando una precisión de primer orden en el espacio y de segundo orden en el tiempo. Hasta donde saben los autores, este es el primer estimación de error a priori de este tipo para el método de punto medio con masa agrupada, incluso en ausencia de DMI.
3. Análisis de Resolvedores Inexactos: Reconociendo que el sistema no lineal debe resolverse iterativamente en la práctica, el artículo analiza el impacto de los errores del resolvedor. Demuestra que si la iteración de punto fijo se termina basándose en un criterio de parada de residuo específico, el esquema inexacto resultante mantiene la misma tasa de convergencia que el esquema ideal, siempre que la tolerancia del resolvedor sea suficientemente pequeña en relación con los parámetros de discretización.

4. Resultados

• Límites Teóricos: El principal resultado teórico (Teorema 3.7 y Teorema 4.8) establece que el error en la norma de energía está acotado por $C(h + k^2)$ más un término dependiente de la tolerancia del resolvedor. Específicamente, para una tolerancia de resolvedor uniforme $\varepsilon_{fp}$, la condición $\varepsilon_{fp} = O(h)$ es suficiente para mantener la tasa de convergencia óptima.
• Experimentos Numéricos:
  • Tasas de Conververgencia: Pruebas numéricas sobre una solución suave confirman las tasas de convergencia teóricas: aproximadamente primer orden en $H^1$ y segundo orden en $L^2$ para el refinamiento espacial, y segundo orden para el refinamiento temporal.
  • Tolerancia del Resolvedor: Los experimentos variando la tolerancia del punto fijo demuestran que el error escala linealmente con $\varepsilon_{fp}/h$, validando la predicción teórica respecto a los errores inducidos por el resolvedor.
  • Preservación de la Estructura: Simulaciones que involucran un skyrmion de tipo Bloch y un experimento de "burbujeo" de alto gradiente confirman que el método preserva la restricción de longitud unitaria nodal hasta la precisión de la máquina y captura correctamente el comportamiento de disipación de energía.
  • Efectos de la DMI: Los experimentos ilustran la influencia asimétrica del signo de la DMI ($\kappa_d$) en la dinámica. Un $\kappa_d$ negativo acelera la concentración de gradientes (formación de burbujas), mientras que un $\kappa_d$ positivo tiende a suprimir la concentración y promover la relajación, consistente con la naturaleza quiral de la interacción.

5. Significado y Reivindicaciones
El artículo afirma proporcionar el primer análisis riguroso de error a priori para el método de elementos finitos de punto medio con masa agrupada aplicado a la ecuación LLG con DMI de volumen. Su importancia radica en cerrar la brecha entre las propiedades de estabilidad conocidas de este método (preservación exacta de la restricción y disipación de energía) y las tasas de convergencia cuantitativas.

Los autores enfatizan que su análisis no depende de condiciones de malla restrictivas (por ejemplo, condiciones de ángulo) o de restricciones de paso de tiempo estrictas que suelen afectar a otros esquemas que preservan la estructura. Además, al demostrar que el resolvedor inexacto práctico no degrada la tasa de convergencia, este trabajo proporciona una justificación rigurosa para el uso de este algoritmo computacionalmente eficiente y preservador de la estructura en la simulación de fenómenos magnéticos quirales complejos, tales como los skyrmions. Los resultados respaldan la robustez del método para manejar la delicada interacción entre el intercambio, la DMI y las condiciones de contorno quirales.