arrow
Volume 16, Issue 1
A Second-Order Semi-Implicit Method for the Inertial Landau-Lifshitz-Gilbert Equation

Panchi Li, Lei Yang, Jin Lan, Rui Du & Jingrun Chen

Numer. Math. Theor. Meth. Appl., 16 (2023), pp. 182-203.

Published online: 2023-01

Export citation
  • Abstract

Electron spins in magnetic materials have preferred orientations collectively and generate the macroscopic magnetization. Its dynamics spans over a wide range of timescales from femtosecond to picosecond, and then to nanosecond. The Landau-Lifshitz-Gilbert (LLG) equation has been widely used in micromagnetics simulations over decades. Recent theoretical and experimental advances have shown that the inertia of magnetization emerges at sub-picosecond timescales and contributes significantly to the ultrafast magnetization dynamics, which cannot be captured intrinsically by the LLG equation. Therefore, as a generalization, the inertial LLG (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves as a hyperbolic system at sub-picosecond timescales, while behaves as a parabolic system at larger timescales spanning from picosecond to nanosecond. Such hybrid behaviors impose additional difficulties on designing efficient numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second-order temporal derivative of magnetization is approximated by the standard centered difference scheme, and the first-order temporal derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with second-order accuracy. At each time step, the unconditionally unique solvability of the unsymmetric linear system is proved with detailed discussions on the condition number. Numerically, the second-order accuracy of the proposed method in both time and space is verified. At sub-picosecond timescales, the inertial effect of ferromagnetics is observed in micromagnetics simulations, in consistency with the hyperbolic property of the iLLG model; at nanosecond timescales, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model.

  • AMS Subject Headings

35Q99, 65Z05, 65M06

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{NMTMA-16-182, author = {Li , PanchiYang , LeiLan , JinDu , Rui and Chen , Jingrun}, title = {A Second-Order Semi-Implicit Method for the Inertial Landau-Lifshitz-Gilbert Equation}, journal = {Numerical Mathematics: Theory, Methods and Applications}, year = {2023}, volume = {16}, number = {1}, pages = {182--203}, abstract = {

Electron spins in magnetic materials have preferred orientations collectively and generate the macroscopic magnetization. Its dynamics spans over a wide range of timescales from femtosecond to picosecond, and then to nanosecond. The Landau-Lifshitz-Gilbert (LLG) equation has been widely used in micromagnetics simulations over decades. Recent theoretical and experimental advances have shown that the inertia of magnetization emerges at sub-picosecond timescales and contributes significantly to the ultrafast magnetization dynamics, which cannot be captured intrinsically by the LLG equation. Therefore, as a generalization, the inertial LLG (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves as a hyperbolic system at sub-picosecond timescales, while behaves as a parabolic system at larger timescales spanning from picosecond to nanosecond. Such hybrid behaviors impose additional difficulties on designing efficient numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second-order temporal derivative of magnetization is approximated by the standard centered difference scheme, and the first-order temporal derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with second-order accuracy. At each time step, the unconditionally unique solvability of the unsymmetric linear system is proved with detailed discussions on the condition number. Numerically, the second-order accuracy of the proposed method in both time and space is verified. At sub-picosecond timescales, the inertial effect of ferromagnetics is observed in micromagnetics simulations, in consistency with the hyperbolic property of the iLLG model; at nanosecond timescales, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model.

}, issn = {2079-7338}, doi = {https://doi.org/10.4208/nmtma.OA-2022-0080}, url = {http://global-sci.org/intro/article_detail/nmtma/21348.html} }
TY - JOUR T1 - A Second-Order Semi-Implicit Method for the Inertial Landau-Lifshitz-Gilbert Equation AU - Li , Panchi AU - Yang , Lei AU - Lan , Jin AU - Du , Rui AU - Chen , Jingrun JO - Numerical Mathematics: Theory, Methods and Applications VL - 1 SP - 182 EP - 203 PY - 2023 DA - 2023/01 SN - 16 DO - http://doi.org/10.4208/nmtma.OA-2022-0080 UR - https://global-sci.org/intro/article_detail/nmtma/21348.html KW - Inertial Landau-Lifshitz-Gilbert equation, semi-implicit scheme, second-order accuracy, micromagnetics simulations. AB -

Electron spins in magnetic materials have preferred orientations collectively and generate the macroscopic magnetization. Its dynamics spans over a wide range of timescales from femtosecond to picosecond, and then to nanosecond. The Landau-Lifshitz-Gilbert (LLG) equation has been widely used in micromagnetics simulations over decades. Recent theoretical and experimental advances have shown that the inertia of magnetization emerges at sub-picosecond timescales and contributes significantly to the ultrafast magnetization dynamics, which cannot be captured intrinsically by the LLG equation. Therefore, as a generalization, the inertial LLG (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves as a hyperbolic system at sub-picosecond timescales, while behaves as a parabolic system at larger timescales spanning from picosecond to nanosecond. Such hybrid behaviors impose additional difficulties on designing efficient numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second-order temporal derivative of magnetization is approximated by the standard centered difference scheme, and the first-order temporal derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with second-order accuracy. At each time step, the unconditionally unique solvability of the unsymmetric linear system is proved with detailed discussions on the condition number. Numerically, the second-order accuracy of the proposed method in both time and space is verified. At sub-picosecond timescales, the inertial effect of ferromagnetics is observed in micromagnetics simulations, in consistency with the hyperbolic property of the iLLG model; at nanosecond timescales, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model.

Panchi Li, Lei Yang, Jin Lan, Rui Du & Jingrun Chen. (2023). A Second-Order Semi-Implicit Method for the Inertial Landau-Lifshitz-Gilbert Equation. Numerical Mathematics: Theory, Methods and Applications. 16 (1). 182-203. doi:10.4208/nmtma.OA-2022-0080
Copy to clipboard
The citation has been copied to your clipboard