Skip to content
DraconiticResearch
Search Sources API
203,234 studies31 sourcesUpdated 13 min ago
PreprintCC BY-NC-ND 4.0arXiv

An unconditionally convergent BDF2-type finite element integrator for the inertial Landau-Lifshitz-Gilbert equation

Abstract

We propose and analyze a fully discrete numerical scheme to approximate solutions to the inertial Landau-Lifshitz-Gilbert (iLLG) equation, a state-of-the-art model to describe the dynamics of the magnetization at subpicosecond time scales. The method employs first-order finite elements for the spatial discretization and the two-step backward differentiation formula (BDF2) for the time discretization, where the inertial term is treated by a convex combination of BDF2 and the backward Euler method. We prove well-posedness and unconditional stability of the scheme, and obtain explicit bounds for the error in the realization of the iLLG-inherent unit length constraint satisfied by the magnetization. Moreover, we derive the discrete energy law satisfied by the approximations, which resembles the inherent dissipative energy law of the model. Finally, with a constructive approach, we prove unconditional convergence of the approximations toward a weak solution of the problem. Numerical experiments assess the performance of the scheme and validate the theoretical findings.

Related

  1. BDF2-type integrator for Landau-Lifshitz-Gilbert equation in micromagnetics: IMEX time-stepping and full effective field
    Michele Aldé, Dirk Praetorius, Agus L. Soenjaya
    arXivPreprintCC BY-NC-ND 4.0Full textNumerical analysis & scientific computingarXiv
  2. An explicit third-order low-regularity integrator for the cubic nonlinear Schrödinger equation
    Shengxin Luo, Bin Wang, Yan Wang
    arXivPreprintarXiv licenceFull textNumerical analysis & scientific computingarXiv
  3. Asymptotic stability of Bloch wall for dispersive Landau-Lifshitz equation
    Slim Ibrahim, Christof Melcher, Ikkei Shimizu
    arXivPreprintarXiv licenceFull textAnalysis & PDEsarXiv
  4. Finite-energy travelling waves for the easy-plane Landau-Lifshitz equation
    Rowan Juneman, Ze Li, Juncheng Wei
    arXivPreprintarXiv licenceFull textAnalysis & PDEsarXiv
  5. The stochastic Landau-Lifshitz-Baryakhtar equation in critical spaces
    Federico Butori, Foivos Evangelopoulos-Ntemiris, Lorenzo Marino +1
    arXivPreprintCC BY 4.0Full textProbabilityarXiv
BibTeX
@misc{alde2026unconditionally,
  title = {{An unconditionally convergent BDF2-type finite element integrator for the inertial Landau-Lifshitz-Gilbert equation}},
  author = {Michele Aldé and Hywel Normington and Dirk Praetorius and Michele Ruggeri},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.35382},
  eprint = {2610.35382},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.35382},
  note = {CC BY-NC-ND 4.0},
}
RIS
TY  - UNPB
TI  - An unconditionally convergent BDF2-type finite element integrator for the inertial Landau-Lifshitz-Gilbert equation
AU  - Michele Aldé
AU  - Hywel Normington
AU  - Dirk Praetorius
AU  - Michele Ruggeri
PY  - 2026
DA  - 2026/10/29
PB  - arXiv
DO  - 10.48550/arxiv.2610.35382
UR  - https://arxiv.org/abs/2610.35382
AB  - We propose and analyze a fully discrete numerical scheme to approximate solutions to the inertial Landau-Lifshitz-Gilbert (iLLG) equation, a state-of-the-art model to describe the dynamics of the magnetization at subpicosecond time scales. The method employs first-order finite elements for the spatial discretization and the two-step backward differentiation formula (BDF2) for the time discretization, where the inertial term is treated by a convex combination of BDF2 and the backward Euler method. We prove well-posedness and unconditional stability of the scheme, and obtain explicit bounds for the error in the realization of the iLLG-inherent unit length constraint satisfied by the magnetization. Moreover, we derive the discrete energy law satisfied by the approximations, which resembles the inherent dissipative energy law of the model. Finally, with a constructive approach, we prove unconditional convergence of the approximations toward a weak solution of the problem. Numerical experiments assess the performance of the scheme and validate the theoretical findings.
ER  -