Z型结纳米轨道中垂直磁各向异性工程调控斯格明子热稳定性与脉冲电流驱动动力学

《Nanomaterials》：Perpendicular Magnetic Anisotropy Engineering for Skyrmion Thermal Stability and Pulsed-Current-Driven Dynamics in a Z-Junction Nanotrack

　　

磁性斯格明子（magnetic skyrmions）是纳米尺度自旋电子器件的潜在信息载体，但其在受限结构中的应用需要同时控制拓扑存活与电流驱动运动。研究人员采用微磁学模拟研究了Z型结纳米轨道中的垂直磁各向异性（perpendicular magnetic anisotropy, PMA）工程。在300 K下，针对结深50 nm、结宽5 nm、脉冲电流密度1.8 × 10¹² A/m²的条件，将各向异性常数K_u从0.70变化至1.00 MJ/m³。每个有限温度条件均考虑五个独立的热实现。在63 ns观测区间内，K_u = 0.70–0.90 MJ/m³未观察到湮灭，而0.95和1.00 MJ/m³的全部五个实现均在约9–10 ns内发生湮灭。0.90 MJ/m³情形尽管受到强烈压缩，仍保持拓扑完整。额外的面内网格细化重现了0.90 MJ/m³的存活与1.00 MJ/m³的湮灭。随后研究人员考察了所选紧凑参考条件K_u = 0.85 MJ/m³在0至500 K范围内的行为，并以0 K的确定性参考为基准。轨迹显示在可用观测范围内斯格明子类拓扑态得以保持，同时伴随温度依赖的位置、半径和速度变化。这些结果将有限时间拓扑存活与动力学可重复性区分开来。温度扫描描述了固定参数模型内的热噪声效应；它并未建立热激活寿命、普适最优各向异性或实验器件工作温度范围。

**研究背景与问题**

磁性斯格明子（magnetic skyrmions）是一类具有连续旋转局域磁化和非零拓扑电荷的纳米尺度非共线自旋织构，其粒子性行为、紧凑尺寸以及可被产生、输运和湮灭的能力使其成为高密度存储、逻辑及非常规计算器件的潜在信息载体。斯格明子基器件的可行性依赖于受限磁性纳米结构中孤立斯格明子的可靠成核与受控运动。在超薄磁性纳米轨道中，斯格明子态由交换相互作用、垂直磁各向异性（PMA）、退磁能与界面Dzyaloshinskii–Moriya相互作用（DMI）之间的平衡决定。然而，受限纳米轨道中的可靠斯格明子输运仍面临挑战：电流驱动的斯格明子会产生纵向与横向运动，由此产生的斯格明子霍尔效应（skyrmion Hall effect）使其向轨道一侧偏转，增加边缘相互作用与湮灭概率，在窄轨道、弯曲段、结区等几何非均匀区域尤为显著。几何工程因此成为控制斯格明子运动的重要手段，但几何不连续并非被动障碍，而是产生与斯格明子尺寸、轨迹和内部结构相互作用的非均匀磁能景观。温度进一步引入限制，因为热涨落持续扰动局域磁化、斯格明子边界与电流驱动轨迹。垂直磁各向异性是控制面外磁化能量偏好并改变反转斯格明子核尺寸的关键设计参数：适度增加K_u可产生更紧凑的斯格明子，但过度各向异性可能减小斯格明子尺寸及其容纳电流诱导和热诱导变形的能力。此前研究已确立磁各向异性、饱和磁化强度、DMI、温度和纳米线尺寸对斯格明子成核与热稳定性的影响，并单独研究了Z型结纳米线中横向畴壁的手性依赖钉扎与脱钉，但这些研究未能确定斯格明子如何响应各向异性控制的压缩、脉冲电流激发与Z型结非对称受限的联合效应。本研究即针对这一耦合展开。

**研究内容与结论**

研究人员通过微磁学模拟考察了Z型结纳米轨道中垂直磁各向异性对斯格明子存活与脉冲电流驱动动力学的影响，采用拓扑电荷、存活与湮灭时间、核尺寸、轨迹和磁构型进行分析，特别区分了在观测区间内存活的强压缩斯格明子与失去拓扑的构型，并通过温度依赖分析进一步考察拓扑保持如何与受限运动和内部变形的变化共存。结果表明，在模拟的300 K条件下，K_u = 0.70–0.90 MJ/m³的每个取值全部五个实现均在63 ns观测区间内存活，而0.95和1.00 MJ/m³的全部五个实现均在约9–10 ns内湮灭；0.90 MJ/m³处的强压缩与拓扑存活共存。当面内单元宽度从2.0 nm减至0.5 nm时，额外网格比较重现了0.90 MJ/m³的存活与1.00 MJ/m³的湮灭。斯格明子向结区运动，在驱动期间经历受限位置偏移，并在脉冲序列结束后向后弛豫。所选K_u = 0.85 MJ/m³条件为温度依赖研究提供了紧凑的存活参考，而非已证明的联合最优。在所考察温度范围内，拓扑保持与位置、半径和速度行为的变化共存，非单调速度曲线未确立热共振。这些结论涉及指定模型中的有限时间行为：模拟区间内的存活并非实测热激活寿命，近似最小半径尺度并非普适坍缩阈值；由于磁学系数固定且焦耳热未自洽计算，规定温度扫描不能确定实验器件的工作范围或长期保持能力。定量结果适用于所研究的d = 50 nm、λ = 5 nm的Z型结，不同结尺寸可改变受限、变形及与边界的相互作用，并可能移动各向异性依赖的存活与湮灭结果，需系统几何研究后方可移植；结果亦特定于J = 1.8 × 10¹² A/m²，未进行电流密度扫描，故不能推断定量电流依赖阈值。超薄磁性多层膜中已证实室温Néel斯格明子及电流驱动运动，这些平台提供了实验背景，但其自旋轨道驱动响应并非对本文Zhang–Li模型的直接标定；特定Z型结的制备与成像（包括边缘粗糙度和加热效应）仍是实验验证所必需的。DMI、饱和磁化强度、交换刚度和阻尼同样影响核结构与驱动响应，由于这些参数固定，所选各向异性不应视为跨材料的普适值，需系统改变这些参数、几何与驱动协议以建立更广泛的设计准则。

**关键技术方法**

研究采用GPU加速的MuMax3微磁学模拟软件包，基于有限差分离散化求解时间和空间依赖的磁化动力学，通过包含绝热与非绝热自旋转移矩项的Landau–Lifshitz–Gilbert方程描述磁化时间演化，有效场由总磁能泛函变分获得，总能量包含交换能、各向异性能、界面DMI能、退磁能与Zeeman能。模拟结构为总尺寸600 × 140 × 4 nm³的垂直磁各向异性纳米轨道，含Z型结（结深d = 50 nm、结宽λ = 5 nm），计算区域划分为300 × 70 × 2个有限差分单元（单元尺寸2 × 2 × 2 nm³），并采用1.0 × 1.0 × 2.0 nm³与0.5 × 0.5 × 2.0 nm³网格进行面内离散化敏感性比较；几何以图像定义形状导入，所有方向施加开放边界条件。材料参数固定为M_s = 8.75 × 10⁵ A/m、A = 1.0 × 10^?11 J/m、D = 1.5 × 10^?3 J/m²、α = 0.05，易轴沿(0,0,1)，沿+z方向施加0.1 T均匀稳定磁场，K_u在0.70–1.00 MJ/m³间系统变化。初始Néel型斯格明子经0 K、零电流、K_u = 0.60 MJ/m³弛豫后逐步增加各向异性制备，并在0 K、零电流下演化1 ns确认稳定性。室温效应通过MuMax3中基于涨落–耗散定理的随机热场引入，工作温度300 K，热平衡1 ns后另加2 ns零电流区间确认稳定；每个各向异性及每个非零温度均以不同随机种子进行五次独立模拟，0 K为确定性参考。驱动采用沿+x方向的单极方波脉冲电流，J_max = 1.8 × 10¹² A/m²，自旋极化率P = 0.40、非绝热参数ξ = 0.05，脉冲频率1.3 GHz、占空比50%，脉冲序列施加50 ns，随后零电流继续10 ns考察向后弛豫、钉扎、存活或湮灭，分析室温区间共63 ns。斯格明子性质通过拓扑电荷Q（Q ≈ ?1为完整斯格明子，Q → 0为湮灭后拓扑平庸态）、最大连通反转核区质心位置、等效半径（与核面积相同的圆半径）及其观测区间最小值、湮灭时间（|Q| < 0.1首次出现并持续失去斯格明子类平台）评估，未达湮灭判据记为右删失观测τ_ann > 63 ns。温度依赖模拟在K_u = 0.85 MJ/m³下于0、100、200、300、400、500 K进行，拓扑电荷轨迹经101点中心移动平均滤波，坐标轨迹经七点三阶Savitzky–Golay滤波后数值微分计算速度分量与速率，受限运动区间为10 ≤ t ≤ 53 ns，并统计均方根速率、95百分位速率及速度分量时间标准差。

**研究结果**

3.1 室温各向异性筛选与两阶段策略：第一阶段评估300 K下各向异性依赖的存活与压缩，第二阶段在所选K_u = 0.85 MJ/m³条件下改变规定温度；共同几何与驱动协议使这些响应可在同一模型内比较，而观测窗口和有限种子数限制了统计与寿命推断。

3.2 垂直磁各向异性对斯格明子存活的影响：在T = 300 K、d = 50 nm、λ = 5 nm、J = 1.8 × 10¹² A/m²下，K_u = 0.70–0.90 MJ/m³的每个采样各向异性全部五个实现均存活至63 ns；0.90 MJ/m³时核强烈压缩但观测存活分数仍为5/5。0.95和1.00 MJ/m³时全部五个实现均湮灭。因此增加垂直各向异性并不产生有限时间存活的单调改善；0.90 MJ/m³存活与两个更高值湮灭的区分适用于采样条件，该离散扫描既未确定连续边界的精确位置，也未确定平衡相图。

3.3 垂直磁各向异性对存活与湮灭时间的影响：K_u = 0.70–0.90 MJ/m³的上箭头表示63 ns观测区间内未观察到湮灭，构成右删失观测，仅能确立τ > 63 ns，不指定63 ns事件时间，也不确立存活情形间的寿命排序。0.95和1.00 MJ/m³的五个实现均在约9–10 ns内湮灭，其均值与标准差由各向异性的五个单独湮灭时间计算；两个高各向异性组间的微小间隔未用于确立统计分辨的湮灭时间差异。这些观测涉及施加脉冲协议下的存活，不提供激活势垒、尝试频率或外推至模拟区间之外的热激活寿命。

3.4 拓扑电荷演化：K_u = 0.70、0.85、0.90 MJ/m³时轨迹保持在斯格明子类平台Q ≈ ?1附近；0.95和1.00 MJ/m³时平台丢失，Q趋近于零，随后持续接近零将湮灭与否则存活的斯格明子的瞬态变形区分开来。图中9.54和9.21 ns的垂直参考线分别为0.95和1.00 MJ/m³的阈值穿越参考时间，并非系综平均事件时间；显示轨迹间的变化不改变存活与湮灭的定性区分，补充图S1的额外网格比较重现了0.90和1.00 MJ/m³的这一区分。

3.5 垂直磁各向异性对斯格明子轨迹的影响：电流约3 ns施加后斯格明子向Z型结运动，在10 ≤ t ≤ 53 ns期间于几何不连续处反复前后运动，而非持续穿过结区；这些偏移的幅度和时序随各向异性变化。0.90 MJ/m³处相对较大的偏移与拓扑电荷保持共存，因此仅位置变化不构成湮灭证据；结诱导受限与自旋转移矩、热涨落和内部变形共同决定斯格明子轨迹。脉冲序列约53 ns结束后，斯格明子中心位置向较小纵向位置移动，表明远离结区的向后弛豫；0.80 MJ/m³位置轨迹约在61 ns结束，未外推至63 ns，该较短位置记录本身不决定单独记录的拓扑电荷信号的端点。

3.6 垂直磁各向异性对最小等效半径的影响：K_u从0.70增至0.85 MJ/m³时，所绘半径从约16 nm减至约9 nm而斯格明子在观测区间内存活；0.90 MJ/m³时达到约5.2 nm而未失去斯格明子类拓扑平台。0.95和1.00 MJ/m³的构型在湮灭前也达到小的可检测半径，但其半径并非存活斯格明子的平衡尺寸；存活与湮灭情形相似的约5 nm尺度表明半径本身不是唯一的坍缩判据，其解释依赖于磁构型、几何、离散化和可检测反转核的定义。增加K_u有利于面外磁化并减小反转核空间范围，交换、DMI、静磁相互作用和外加场也贡献于变形构型，因此观测到的压缩作为该平衡的一部分讨论，而非作为单独能垒的测量或DMI在某一普适半径处变得不足的证明；Q(t)网格比较未独立确定数值半径值的收敛性。

3.7 各向异性依赖存活的磁构型证据：初始核随各向异性增加变得更紧凑；后续构型显示与结区相互作用期间的变形，包括较高各向异性下结内边界附近的强非对称压缩。图7中面板(f)和(h)保留局域反转核，代表接近湮灭过程的构型而非完全湮灭的最终态，拓扑丢失由Q(t)演化单独确立。构型与电荷轨迹共同与边界影响的湮灭过程一致，但未分离出唯一的微观路径或测量其能垒。核心区别在于斯格明子可经历显著压缩而保留拓扑，而较高各向异性下电流驱动的边界相互作用反而可导致斯格明子类态的不可逆丢失；该解释不赋予普适临界半径，也不将内禀霍尔偏转作为结相互作用的唯一来源。

3.8 K_u = 0.85 MJ/m³下温度依赖的斯格明子存活与动力学：在0至500 K的可用观测范围内斯格明子类电荷平台始终存在，无持续向湮灭态的转变；平均电荷随温度升高变得更负，在有限离散样品中这些非整数值不应解释为额外斯格明子或整数拓扑数的增加。纵向中心位置显示向结区运动后受限位置偏移，脉冲序列约53 ns结束后斯格明子向后弛豫；高温下偏移变得更不规则，但其详细序列并非随温度单调，400 K位置轨迹约在61.9 ns结束，未用于推断63 ns时的最终位置。等效半径显示接近结区时初始压缩，随后反复膨胀与收缩；修正后的数值汇总给出受限区间平均半径在0 K约10.08 nm、500 K约13.24 nm，这些值指10 ≤ t ≤ 53 ns而非最小半径或最终弛豫构型；高温下较大的尺寸变化与观测区间内拓扑保持共存。有限温度下曲线及其离散度描述五实现分析，0 K模拟为确定性参考故不显示实现间不确定度，T > 0时曲线为五个独立热实现的系综均值、离散度由其标准差表示；温度比较限于规定噪声的固定参数模型，并非高达500 K实验热运行的证明。

3.9 斯格明子速率与速度涨落的温度依赖：由于受限斯格明子反复改变方向，小的有符号平均速度可与显著速率和涨落共存，因此较大的均方根速率代表更强的局域运动，而不一定是更有效的穿过结区输运。显示的中心值总体随温度增加但呈非单调行为，最大值出现在400 K；纵向涨落在整个温度范围内始终大于横向涨落。这些趋势表征受限斯格明子的温度依赖动力学响应，400 K处观测到的最大值不被解释为统计确立的最优或热共振的证据。速度分量由斯格明子中心坐标的时间导数获得，直接以m/s报告（1 nm/ns = 1 m/s），无额外速度归一化；σ(v_x)和σ(v_y)分别表征纵向和横向速度分量的时间涨落。图8与图9共同表明斯格明子类态的保持可与位置、尺寸和局域运动的显著变化共存，确立特定耦合或共振机制需要单独的时间分辨或频域分析。

**讨论与结论**

讨论部分强调，有限时间拓扑存活与动力学可重复性应加以区分：幸存情形τ > 63 ns表示模拟区间内未观察到湮灭，湮灭情形的时间由|Q| < 0.1后持续失去斯格明子平台定义，半径描述湮灭前的可检测核；“过渡”表示强压缩下的存活，而非中间存活分数。最小等效斯格明子半径随垂直各向异性增加而系统减小，湮灭情形的半径对应拓扑坍缩前的最小可检测核半径而非平衡半径。结论部分翻译如下：研究人员研究了Z型结纳米轨道中垂直磁各向异性对斯格明子存活与脉冲电流驱动动力学的影响。对于模拟的300 K条件，K_u = 0.70–0.90 MJ/m³的每个取值全部五个实现均在63 ns观测区间内存活，而0.95和1.00 MJ/m³的全部五个实现均在约9–10 ns内湮灭；0.90 MJ/m³处的强压缩与拓扑存活共存。当面内单元宽度从2.0 nm减至0.5 nm时，额外网格比较重现了0.90 MJ/m³的存活与1.00 MJ/m³的湮灭。斯格明子向结区运动，在驱动期间经历受限位置偏移，并在脉冲序列结束后向后弛豫。所选K_u = 0.85 MJ/m³条件为温度依赖研究提供了紧凑的存活参考，而非已证明的联合最优。在所考察温度范围内，拓扑保持与位置、半径和速度行为的变化共存；非单调速度曲线未确立热共振。这些结论涉及指定模型中的有限时间行为：模拟区间内的存活并非实测热激活寿命，近似最小半径尺度并非普适坍缩阈值。由于磁学系数固定且焦耳热未自洽计算，规定温度扫描不能确定实验器件的工作范围或长期保持能力。定量结果适用于所研究的d = 50 nm、λ = 5 nm的Z型结；不同结尺寸可改变受限、变形及与边界的相互作用，并可能移动各向异性依赖的存活与湮灭结果，需系统几何研究后方可将定量结果移植至其他设计。结果亦特定于J = 1.8 × 10¹² A/m²：增加电流增强自旋转移矩，可能增强边界附近的位移与变形，减小电流则可能减弱位移或增强钉扎；未进行电流密度扫描，故不推断定量电流依赖阈值。超薄磁性多层膜中已证实室温Néel斯格明子及电流驱动运动，这些平台提供了实验背景，但其自旋轨道驱动响应并非对本文Zhang–Li模型的直接标定；特别是名义5 nm结特征要求对边缘几何、厚度和电流分布的严苛控制，特定Z型结的制备与成像（包括边缘粗糙度和加热效应）仍是实验验证所必需的。DMI、饱和磁化强度、交换刚度和阻尼同样影响核结构与驱动响应；由于这些参数固定，所选各向异性不应视为跨材料的普适值，需系统改变这些参数、几何与驱动协议以建立更广泛的设计准则。