合成反铁磁斯格明子（SAF Skyrmion）的本征模式研究

《Scientific Reports》：Eigenmodes of synthetic antiferromagnetic skyrmions

　　

摘要专业翻译 研究人员利用微磁学本征值方法和衰减振荡（ringdown）模拟，系统研究了受限合成反铁磁（SAF）斯格明子（skyrmion）的激发模式。从铁磁（FM）单层中单个斯格明子的最低频模式（回旋模式（gyrotropic mode）和呼吸模式（breathing mode））出发，研究了反铁磁层间耦合如何改变SAF双层膜中的动力学行为。研究考虑了多种几何构型：方形和矩形受限结构中的单个SAF斯格明子、不等层厚结构以及含多个斯格明子的条带结构。反铁磁耦合显著改变了低频动力学特性。方形几何构型呈现出两个近简并的回旋模式，其中每一模式中两层均具有相同的旋转方向。而在矩形几何构型中，研究人员发现由两层中相反回旋方向产生的近线性SAF斯格明子平移模式（translational mode），这些平移模式成为SAF斯格明子链的特征低频激发。对于斯格明子链，研究人员识别出具有驻波状空间轮廓的集体平移模式和呼吸模式。除类铁磁呼吸模式外，SAF几何构型还支持两层反相振荡的呼吸模式。研究进一步展示了信号沿扩展SAF斯格明子链的传播，其传播速度比铁磁斯格明子链约快15%。这些结果为几何受限、层内耦合与层间耦合相互作用下SAF斯格明子的集体动力学提供了系统性的描述。

论文解读文章

一、研究背景与意义

磁性斯格明子（skyrmion）是具有拓扑非平凡特性的纳米尺度自旋纹理，在自旋电子学（spintronics）和磁子学（magnonics）领域具有广阔的应用前景，包括赛道存储器（racetrack memory）、信号处理用斯格明子晶体以及磁子逻辑门等。斯格明子或其集合体拥有丰富的动力学模式，包括回旋模式（gyrotropic mode，斯格明子结构区域围绕平衡位置轨道运动）和呼吸模式（breathing mode，斯格明子围绕中心周期性膨胀和收缩），这些模式及其集体激发已被广泛研究。

然而，铁磁斯格明子在实际器件应用中面临若干限制：非局域杂散场（stray field）使其在室温下的最小尺寸受限在数十纳米量级，且通常需要外加磁场来维持其稳定性。合成反铁磁体（synthetic antiferromagnet, SAF）作为一类能够承载斯格明子构型的新型材料体系，具有两个关键优势：其一，净磁化强度为零大幅降低了杂散场，有利于器件中磁性元件的紧密排列并减小斯格明子尺寸；其二，SAF斯格明子的净拓扑电荷为零，可消除斯格明子霍尔效应（skyrmion Hall effect），预测可实现沿电流方向的直线运动且速度优于铁磁对应体系。已有实验证明了SAF斯格明子在室温下的稳定化以及热激活扩散行为。

尽管单个SAF斯格明子在圆盘状几何中的回旋模式和呼吸模式已有研究，但包含多个斯格明子的SAF体系的集体动力学尚未得到系统探索。为此，研究人员将研究扩展至SAF材料中的多斯格明子体系，采用双层SAF模型，同时研究层内（斯格明子-斯格明子）和层间（反铁磁）耦合效应。

二、研究内容与主要结论

研究人员系统研究了受限SAF斯格明子及SAF斯格明子链的本征模式，涵盖多种几何构型：方形和矩形受限结构中的单个斯格明子、不等层厚结构、含两个、五个及二十五个斯格明子的条带结构。主要结论如下：

铁磁方形几何中的斯格明子。 在铁磁方形几何中，单个奈尔型（Néel）斯格明子呈现三个本征模式：顺时针回旋模式（0.84 GHz）、呼吸模式（15.68 GHz）和逆时针旋转模式（43.63 GHz）。回旋模式的旋转方向由斯格明子拓扑电荷的符号决定，改变斯格明子核心极性可反转旋转方向。

SAF方形几何中的斯格明子。 在SAF方形几何中，两层中的斯格明子通过层间RKKY耦合相互作用，呈现两个近简并的回旋模式（4.68 GHz和4.75 GHz），每个模式中两层斯格明子具有相同的旋转方向。由于顶层斯格明子（核心向下）偏好顺时针旋转而底层斯格明子（核心向上）偏好逆时针旋转，系统存在"挫败"（frustration）效应：非偏好旋转方向的层具有减小的轨迹振幅。当层厚不对称时（顶层2 nm、底层1 nm），近简并性被解除，模式分裂为2.1 GHz和8.28 GHz两个回旋模式。此外，SAF体系支持同相呼吸模式（20.37 GHz）和反相呼吸模式（34.12 GHz），后者因RKKY耦合在大界面面积上的高能量代价而具有更高频率。

SAF矩形几何中的斯格明子：平移模式。 当方形对称性被打破（L_x=32 nm，L_y=40 nm），两个回旋模式被平移模式取代：每个层中的斯格明子沿其偏好方向回旋（顶层顺时针、底层逆时针），两层椭圆轨迹的x分量相互抵消，使得SAF斯格明子核心位置呈现近线性运动。沿样本长轴（y方向）的平移模式频率较低，沿短轴（x方向）的平移模式频率较高。随层间耦合强度增大，椭圆轨迹变得更加细长，最终呈现线性运动。在L_x=L_y附近观察到典型的避免交叉（avoided crossing）行为，两种平移模式在对称点附近杂化并恢复为回旋模式。

双斯格明子系统。 在64 nm×40 nm条带中容纳两个斯格明子对，前四个模式（A-D）为平移模式（3-6 GHz），后四个模式（E-H）为呼吸模式（18-35 GHz）。平移模式中相邻斯格明子可呈现同相或反相运动；呼吸模式中，层间同相与反相、左右斯格明子间同相与反相的组合决定了模式的能量高低。随层间耦合强度减小，模式频率整体降低，且部分模式发生重排交叉。

五斯格明子系统。 在五斯格明子条带中，不同类型模式开始形成能带。水平平移模式遵循固定端驻波模型（振幅分布为sin函数形式），而垂直平移模式因条带端部对垂直方向无约束而遵循自由端驻波模型（振幅分布为cos函数形式），且垂直平移模式的模式排序与水平平移模式相反，表明相邻斯格明子间的有效耦合在两种方向上的符号不同。

信号传播。 在25个斯格明子的长链中，约20 GHz处形成顶层-底层同相呼吸模式带。通过高斯脉冲激发最左侧斯格明子，信号沿条带传播并在边缘反射，测得传播速度约300 m/s，最大群速度约330 m/s，比相似几何的铁磁斯格明子链（约255-280 m/s）快约15%，与SAF体系约15%更大的带宽一致。

三、关键技术与方法

研究人员主要采用了以下关键技术方法：微磁学本征值方法（eigenvalue method），通过求解无阻尼线性化朗道-利夫希茨-吉尔伯特（Landau-Lifshitz-Gilbert, LLG）方程获得本征值和本征向量；衰减振荡（ringdown）方法，通过施加外部激发场模拟完整动力学并计算功率谱密度（power spectral density, PSD）；相位-振幅映射可视化技术，用于直观展示各本征模式的空间分布特征；Thiele模型用于描述低频回旋模式；驻波模型用于拟合多斯格明子体系中各斯格明子的集体激发振幅分布。模拟采用双层SAF模型，层间通过RKKY反铁磁耦合连接。

四、讨论与结论

研究人员系统研究了受限合成反铁磁斯格明子及其链的本征模式。在铁磁方形几何中，恢复了此前报道的特征回旋模式和呼吸模式。在SAF方形几何中，发现两个近简并的回旋模式，层内旋转方向一致，并存在同相和反相两种呼吸模式。打破方形对称性后，回旋模式转变为平移模式，两层中斯格明子以相反方向回旋，产生有效的线性运动。对于斯格明子链，识别出具有驻波状空间轮廓的集体平移模式和呼吸模式，其中水平平移遵循固定端边界条件，垂直平移遵循自由端边界条件。SAF斯格明子链的信号传播速度比铁磁对应体系快约15%，归因于SAF体系更大的带宽。这些结果为几何受限、层内耦合与层间耦合相互作用下SAF斯格明子的集体动力学提供了系统描述，对SAF斯格明子在自旋电子学和磁子学器件中的潜在应用具有重要意义。