[論文レビュー] Micromagnetic Monte Carlo method with variable magnetization length based on the Landau-Lifshitz-Bloch equation for computation of large-scale thermodynamic equilibrium states

Serban Lepadatu|arXiv (Cornell University)|Jun 10, 2021

Magnetic properties of thin films参考文献 56被引用数 10

ひとこと要約

本稿では、磁化の回転と長さの変化を許容する試行移動を可能にする、ランダウ＝リフシッツ＝ブロッハ方程式に基づくマイクロ磁気モンテカルロ（MMC）法を提案する。この方法により、大規模な熱力学的平衡状態の計算が効率的に行える。本手法は磁化長さのマクスウェル＝ボルツマン分布を正しく再現しており、GPU上で完全に並列化可能であり、確率的ランダウ＝リフシッツ＝ブロッハ動的アプローチよりも最大19.2倍の高速化を達成する。これにより、グラニュラー膜や人工スピンアイスを含む複雑な磁性系の有限温度シミュレーションに最適である。

ABSTRACT

An efficient method for computing thermodynamic equilibrium states at the micromagnetic length scale is introduced, using the Markov chain Monte Carlo method. Trial moves include not only rotations of vectors, but also a change in their magnetization length. The method is parameterized using the longitudinal susceptibility, reproduces the same Maxwell-Boltzmann distribution as the stochastic Landau-Lifshitz-Bloch equation, and is applicable both below and above the Curie temperature. The algorithm is fully parallel, can be executed on graphical processing units, and efficiently includes the long range dipolar interaction. This method is generally useful for computing finite-temperature relaxation states both for uniform and non-uniform temperature profiles, and can be considered as complementary to zero-temperature micromagnetic energy minimization solvers, with comparable computation time. Compared to the dynamic approach it is shown the micromagnetic Monte Carlo method is up to almost 20 times faster. Moreover, unlike quasi-zero temperature approaches which do not take into account the magnetization length distribution and stochasticity, the method is better suited for structures with unbroken symmetry around the applied field axis, granular films, and at higher temperatures and fields. In particular, applications to finite-temperature hysteresis loop modelling, chiral magnetic thin films, granular magnetic media, and artificial spin ices are discussed.

研究の動機と目的

• 全温度範囲、特にCurie転移温度を超える領域を含む、マイクロ磁気スケールにおける熱力学的平衡状態を効率的に計算する手法の開発。
• ゼロ温度でのエネルギー最小化法や動的sLLBアプローチの限界を克服し、磁化長さのフラクチュエーションを正しく捉えること。
• 磁化長さ分布を適切に扱うことで、有限温度におけるヒステリシスループ、キラル磁性薄膜、グラニュラー媒体、人工スピンアイスの正確なシミュレーションを可能にすること。
• 計算効率と物理的正確性の両立を図り、自由エネルギーの正しいマクスウェル＝ボルツマン分布を再現すること。

提案手法

• 試行移動にマイクロ磁気ベクトルの回転と長さ変更を含むマルコフ連鎖モンテカルロ法を採用する。
• 正しいマクスウェル＝ボルツマン分布への収束を保証するため、詳細つり合い条件を用いて受容確率を導出する。
• 縦磁化率を用いてパrameter化し、確率的ランダウ＝リフシッツ＝ブロッハ方程式の統計的挙動と直接的に関連付ける。
• 磁気分配場は1イテレーションに1回のみ更新され、長距離相互作用の相関が集合全体で緩やかであることに基づく。
• オープンソースのBorisフレームワークに実装されており、完全に並列化可能でGPU加速が可能である。
• 長距離ダイポール相互作用を、1イテレーションに1回のみ磁気分配場を更新することで効率的に統合し、精度を維持しながら速度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 磁化長さの変動を含むマイクロ磁気モンテカルロ法が、有限温度における自由エネルギーの正しいマクスウェル＝ボルツマン分布を再現できるか？
• RQ2 提案されたMMC法の計算時間と正確性について、動的sLLBアプローチと比較してどのように異なるか？
• RQ3 キラル磁性薄膜、グラニュラー磁性媒体、人工スピンアイスといった複雑な系を、有限温度で効率的にシミュレートできるか？
• RQ4 並列MMC設定において、1イテレーションに1回だけ磁気分配場を更新しても、熱力学的正確性に影響を及ぼさないか？
• RQ5 ゼロ温度エネルギー最小化ソルバーと比較して、有限温度状態の物理的正確性を維持したまま、計算時間の点で優れているか？

主な発見

• 楕円体問題において、GPU上でのsLLB動的アプローチ（14,342 s）と比較して、MMC法は最大19.2倍の高速化を達成（747 s）。
• CPU上での2次元楕円問題でも17.6倍の高速化を示し、異なるハードウェア環境でも一貫した性能向上を確認。
• 磁化長さのマクスウェル＝ボルツマン分布が正しく再現されており、統計力学的基礎の妥当性が裏付けられた。
• 1イテレーションに1回の磁気分配場更新でも、各受容移動後に更新する逐次法と区別がつかない結果が得られ、精度を損なわず計算効率が向上していることが確認された。
• 有限温度ヒステリシスループ、DMIを有するキラル磁性薄膜、グラニュラー媒体、人工スピンアイスのシミュレーションに成功し、広範な適用可能性を示した。
• アルゴリズムは完全に並列化可能でGPU実行と互換性があり、従来の動的手法では実現が困難だった大規模シミュレーションを可能にした。