[論文レビュー] Propagating spin-wave normal modes: A dynamic matrix approach using plane-wave demagnetizating tensors

Y. Henry, O. Gladii+1 more|arXiv (Cornell University)|Nov 18, 2016

Magnetic properties of thin films被引用数 8

ひとこと要約

本稿では、動的行列法と平面波磁化率テンソルを用いた有限差分磁気モーメント法を提案し、広がった磁性膜およびストリップにおける励起するスピン波正規モードの計算を可能にする。交換相互作用および長距離に及ぶ磁気双極子相互作用を高精度に扱えるため、二層膜における表面波やドメインウォールに沿った波ガイドモードなど、磁性的に不均一な系における非可逆的スピン波モードの研究が可能となり、解析モデルとの完全な数値的検証が行われた。

ABSTRACT

We present a finite-difference micromagnetic approach for determining the normal modes of spin-waves propagating in extended magnetic films and strips, which is based on the linearized Landau-Lifshitz equation and uses the dynamic matrix method. The model takes into account both short range exchange interactions and long range dipole-dipole interactions. The latter are accounted for through plane-wave dynamic demagnetization factors, which depend not only on the geometry and relative positions of the magnetic cells, as usual demagnetization factors do, but also on the wave vector of the propagating waves. Such a numerical model is most relevant when the spin-wave medium is spatially inhomogeneous perpendicular to the direction of propagation, either in its magnetic properties or in its equilibrium magnetic configuration. We illustrate this point by studying surface spin-waves in magnetic bilayer films and spin-waves channelized along magnetic domain walls in perpendicularly magnetized strips. In both cases, dynamic dipolar interactions produce non-reciprocity effects, where counter-propagative spin-waves have different frequencies.

研究の動機と目的

• 広がった磁性系における励起スピン波正規モードを計算する周波数領域磁気モーメント法の開発を目的とする。
• 時間領域磁気モーメントシミュレーションの限界を克服し、低周波数および簡約状態のモードを検出可能とする。
• 波数依存の磁気双極子相互作用を組み込むことで、動的行列法を平面波の伝搬に拡張する。
• 非同一磁気的構造や物性が空間的に変化する系におけるスピン波モードの研究を可能とする。
• 薄膜およびストリップの極限ケースに対する解析モデルと比較することで、手法の妥当性を検証する。

提案手法

• 本手法は線形化されたランダウ＝リフシッツ方程式に基づき、磁性体をセルに有限差分離散化する。
• 各セルの磁化に作用する有効磁場（交換相互作用および磁気双極子相互作用を含む）を記述する動的行列を構築する。
• 磁気双極子相互作用は、磁気ポテンシャルの積分表現から導出された波数依存の平面波動的磁化率テンソルを用いてモデル化する。
• 動的磁化率テンソルは、系の幾何形状、相対的位置およびスピン波の波数に依存し、振動する磁化状態に一般化された古典的磁化率係数の概念を拡張する。
• 複素行列の固有値分解（Eigenライブラリを用いて）により、スピン波系の固有振動数および固有モードを抽出する。
• C++を用いて実装され、数値積分にはGNU科学計算ライブラリを用い、任意の面内不均一性を持つ広がった膜および薄いストリップの両方をサポートする。

実験結果

リサーチクエスチョン

• RQ1 動的行列法を、封止されたナノ構造における定常波にとどまらず、広がった磁性膜およびストリップにおける励起スピン波モードの計算にどのように拡張できるか？
• RQ2 与えられた波数を持つスピン波に対して、磁気双極子相互作用を正しく定式化し、有限差分磁気モーメントフレームワークにどのように統合できるか？
• RQ3 動的双極子相互作用が、逆方向に伝搬する波が異なる周波数を持つという非可逆的効果をどれほど引き起こすか？
• RQ4 磁性的に不均一な系、例えば二層膜やドメインウォール波ガイドなどにおけるスピン波モードを、本手法はどのように高精度に記述できるか？
• RQ5 特定の幾何形状（既知のスピン波解が得られる薄膜など）に対して、数値結果が解析モデルと一致するか検証可能か？

主な発見

• 本手法は、広がった磁性膜およびストリップにおける励起スピン波モードを正確に計算でき、バックワードボリューム波および表面波の両方において、数値結果と解析解の間で極めて良好な一致が得られた。
• 交換剛性がゼロで、面内磁場が0.1 Tの40 nm厚の膜において、|k| = 60 rad/μmのモードプロファイルが、Damon-Eshbachモデルによる解析的予測と高い精度で一致した。
• k → 0の極限において、動的磁化率テンソルは静的磁化率係数に正確に還元され、従来の理論と整合性があることが確認された。
• セルの高さbが非常に大きい極限では、動的磁化率係数が無限大スラブの値に収束し、準2次元領域における手法の振る舞いが妥当であることが検証された。
• 本手法は、二層膜における表面波やドメインウォールに沿った波ガイドモードなど、磁性的に不均一な系における非可逆的スピン波分散を明らかにした。逆方向に伝搬するモードが異なる周波数を持つことが観測された。
• C++およびEigenライブラリを用いた実装により、複素行列の固有値分解が効率的に行えるようになり、実用的なマグノニクス構造へのスケーラビリティが確保された。