[論文レビュー] A second-order numerical method for Landau-Lifshitz-Gilbert equation with large damping parameters

Yongyong Cai, Jingrun Chen+2 more|arXiv (Cornell University)|May 8, 2021

Magnetic properties of thin films参考文献 28被引用数 11

ひとこと要約

本稿では、大きな減衰パラメータを有するランダウ＝リフシッツ＝ギルバート（LLG）方程式に対して、2次精度で無条件安定な数値解法を提案する。減衰項を調和写像フローに再定式化し、非線形項を完全に明示的に取り扱うことで、各時間ステップで定数係数を有する対称的・正定値線形系のみを解くことができ、高速フーリエ変換（FFT）を用いたスムーズなスolverが可能になる。主な貢献は、大規模な減衰に対しても安定性と物理的精度を維持する計算効率の高い2次精度スキームの構築であり、1次元および3次元のシミュレーション（ドメインウォールの動的挙動を含む）によって検証されている。

ABSTRACT

A second order accurate numerical scheme is proposed and implemented for the Landau-Lifshitz-Gilbert equation, which models magnetization dynamics in ferromagnetic materials, with large damping parameters. The main advantages of this method are associated with the following features: (1) It only solves linear systems of equations with constant coefficients where fast solvers are available, so that the numerical efficiency has been greatly improved, in comparison with the existing Gauss-Seidel project method. (2) The second-order accuracy in time is achieved, and it is unconditionally stable for large damping parameters. Moreover, both the second-order accuracy and the great efficiency improvement will be verified by several numerical examples in the 1D and 3D simulations. In the presence of large damping parameters, it is observed that this method is unconditionally stable and finds physically reasonable structures while many existing methods have failed. For the domain wall dynamics, the linear dependence of wall velocity with respect to the damping parameter and the external magnetic field will be obtained through the reported simulations.

研究の動機と目的

• 現代の磁性材料に一般的に見られるが、既存の数値スキームでは十分に検討されていない大規模な減衰パラメータを有するランダウ＝リフシッツ＝ギルバート（LLG）方程式を解く、数値的に効率的かつ安定な手法の開発。
• 各時間ステップで非対称な線形系を解く必要がある従来の半陽解法の計算ボトルネックを克服し、GMRESなどの反復解法に依存する低効率性を解消すること。
• 時間方向1次精度であるガウス＝セイデル射影法（GSPM）を上回るため、同じ計算複雑度を維持しつつ時間方向2次精度を達成し、高速スケーラブルな解法を可能にすること。
• ドメインウォールの運動やエネルギー散逸といった物理的に意味のある動的挙動を、他の手法が失敗または不安定化する大規模な減衰条件下でも的確に捉える、堅牢な数値フレームワークの提供。

提案手法

• 減衰項を調和写像フローとして表現することで、拡散的挙動のより構造的な取り扱いが可能になるようにLLG方程式を再定式化する。
• 時間微分に2次の後退差分公式（BDF2）を適用し、時間方向に2次精度を確保する。
• 定数係数を有するラプラシアン（拡散）項をBDF2スキームにより陰的に取り扱い、対称的・正定値線形系を形成する。
• すべての非線形項（特にジルバーモーメント項および調和写像フローの非線形部分）を2次補間式を用いて明示的に取り扱い、安定性と効率性を維持する。
• 各時間ステップにおける得られる線形系は対称的かつ正定値であるため、FFTに基づく高速スケーラブルな解法が利用可能となり、計算効率が著しく向上する。
• 調和写像フローによる安定化効果のおかげで、数値実験により示されるように、大規模な減衰パラメータに対しても本手法は無条件安定である。

実験結果

リサーチクエスチョン

• RQ1 大規模な減衰パラメータを有するLLG方程式に対して、2次精度で無条件安定な数値スキームを開発可能か？既存手法はしばしば失敗または不安定化するが、そのような条件下でも有効か？
• RQ2 ベクトル的かつ非線形的性質を持つLLG方程式に対しても、FFTを用いた高速スケーラブルな解法が可能であるという点で、本手法が高い計算効率を維持できるか？
• RQ3 ドメインウォールの速度は減衰パラメータαおよび外部磁場heにどのように依存するか？本スキームはその依存関係を正確に捉えることができるか？
• RQ4 特に3次元シミュレーションにおける大規模な減衰条件下で、本手法はガウス＝セイデル射影法（GSPM）および半陽解法（SIPM）よりも効率的かつ正確であるか？

主な発見

• 本手法は1次元および3次元の両方のシミュレーションで時間方向に2次精度を達成し、テストした全誤差ノルムにおいて収束次数が約2.0に一致する。
• 本手法は大規模な減衰パラメータ（α = 40まで）に対しても無条件安定である一方、GSPMおよびSIPMは同様の条件下で不安定または物理的に意味のない挙動を示す。
• ドメインウォール速度は減衰パラメータαおよび外部磁場heに対して線形に依存し、最小 squares 擬合により得られた勾配は、αあたり0.91～1.02 m/s、heあたり0.91～0.95 m/sの範囲に収まる。
• 本手法はGSPMおよびSIPMに比べて著しく効率的である：1次元ではSIPMより少ないCPU時間を要し、GSPMに比べてはるかに少ない。3次元ではGSPMと同等の性能を示すが、SIPMに比べて同じ精度で著しく高速である。
• エネルギーの時間発展曲線は、α = 2, 5, 8の条件下では3つの手法で一貫した散逸パターンを示す。α = 10ではSIPMにわずかな違いが見られ、本手法の極端な減衰条件下での頑健性が裏付けられる。
• 他の手法が失敗するα = 40でも、安定したドメインウォールや均一な磁化状態といった物理的に妥当な磁化構造を的確に捉えることに成功している。