[论文解读] 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.

研究动机与目标

• 开发一种用于计算扩展磁性系统中传播自旋波本征模的频域微磁学方法。
• 解决时域微磁模拟在检测低频和简并模式方面的局限性。
• 通过引入波矢依赖的偶极-偶极相互作用，将动态矩阵方法扩展至处理传播平面波。
• 实现对具有空间非均匀性的磁介质中自旋波模式的研究，例如具有非共线自旋结构或变化材料参数的体系。
• 通过与薄膜的解析模型及条带的极限情况对比，对方法进行验证。

提出的方法

• 该方法基于线性化Landau-Lifshitz方程，并采用有限差分法将磁性介质离散化为若干单元。
• 构建一个动态矩阵，编码作用于每个单元磁化强度上的有效磁场，包括交换相互作用和偶极-偶极相互作用。
• 偶极-偶极相互作用通过波矢依赖的平面波动态退磁张量建模，该张量由磁静势的积分表达式推导而来。
• 动态退磁张量依赖于自旋波的几何形状、相对位置和波矢，将经典退磁因子概念推广至振荡磁化的情形。
• 采用复矩阵对角化方法（通过Eigen库实现）提取自旋波系统的本征频率和本征模。
• 该方法使用C++实现，借助GNU科学库进行数值积分，支持具有任意面内非均匀性的扩展薄膜和薄条带。

实验结果

研究问题

• RQ1 如何将动态矩阵方法扩展至计算扩展磁性薄膜和条带中的传播自旋波模式，而非仅限于受限纳米结构中的驻波？
• RQ2 对于给定波矢的自旋波，偶极-偶极相互作用的正确表述是什么？如何将其整合入有限差分微磁框架？
• RQ3 动态偶极相互作用在多大程度上诱导自旋波模式的非互易性，即反向传播波具有不同频率？
• RQ4 该方法如何准确描述磁性非均匀体系中的自旋波模式，例如双异质结构或畴壁波导？
• RQ5 数值结果能否针对特定几何结构（如具有已知自旋波解的薄膜）与解析模型进行验证？

主要发现

• 该方法成功计算了扩展磁性薄膜和条带中的传播自旋波模式，数值结果与薄膜中后向体积波和表面波的解析解高度一致。
• 对于厚度为40 nm、交换刚度为零且面内磁场为0.1 T的薄膜，当|k| = 60 rad/μm时，数值模式分布与Damon-Eshbach模型的解析预测高度吻合。
• 在k → 0极限下，动态退磁张量正确地退化为静态退磁因子，验证了与已有理论的一致性。
• 当单元高度b趋于无穷大时，动态退磁因子收敛于无限厚平板的值，验证了该方法在准二维极限下的行为正确性。
• 该方法揭示了非均匀体系中自旋波色散的非互易性，例如双异质结构中的表面波以及畴壁沿线的波导模式，其中反向传播的模式表现出不同的频率。
• 采用C++和Eigen库的实现支持高效的复矩阵对角化，使该方法可扩展用于真实磁子结构的模拟。