各向异性铁磁体中的多模磁孤子动力学及其不稳定性研究——面向自旋电子学应用

《Journal of Magnetism and Magnetic Materials》：Multi-mode magnetic soliton dynamics and instability in anisotropic ferromagnets for spintronic applications

　　R. Ravichandran | Raja Mohan O.S. | E. Parasuraman | Deepthi KVBL
非线性系统中心，金奈技术学院，金奈，泰米尔纳德邦 600 069，印度

**摘要**
本研究探讨了各向异性海森堡铁磁体中的非线性自旋激发的特性，旨在对磁孤子进行分类并评估其在自旋电子设备中的潜力。朗道-利夫希茨（L-L）方程控制着磁孤子的动力学行为，通过使用立体投影技术，该方程可简化为一种导数形式的非线性薛定谔（DNLS）型方程。通过使用Hirota的双线性方法对各种扭结-反扭结孤子碰撞进行解析求解，对铁磁介质进行了建模。此外，我们还提出了一种控制方案，用于调整各向异性能量密度、交换刚度、各向异性常数和饱和磁化参数，这些参数控制着孤子的速度并影响碰撞位置。调制不稳定性（MI）分析确定了特定扰动波数的稳定性和不稳定性区域，从而支持了理论预测。我们的发现揭示了各向异性铁磁体中具有自旋激发的微妙扭结-反扭结碰撞结构。所开发的激发表现出强烈的参数依赖性，为未来的自旋电子技术带来了 optimistic 前景。

**引言**
铁磁纳米结构中的磁孤子是现代自旋电子设备中的关键信息载体。尽管有许多可用的操控技术，但其实际应用仍然具有挑战性。随着量子信息系统越来越多地使用铁磁自旋网络，磁孤子因能够长距离保持其形状和能量而受到重视 [1]、[2]、[3]。因此，由海森堡自旋链控制的磁孤子动力学近年来受到了广泛关注。海森堡型模型现在在凝聚态物理学中起着重要作用，包括磁化切换和磁孤子的动力学行为，特别是在受内在磁相互作用（如交换作用和各向异性）支配的系统中。然而，在本研究中，我们将分析范围限制在内在磁相互作用上。L-L方程描述了一个具有丰富物理特性的经典非线性动力系统 [4]、[5]、[6]。除了各向同性情况外，各向异性扩展也通过黎曼-希尔伯特方法和逆散射方法得到了支持。尽管存在多种描述非线性自旋动力学的数学模型，但离散的海森堡自旋链有效地捕捉了交换耦合和磁各向异性 [7]。

**近期进展**
近期文献报道了在磁性介质中研究非线性孤子的最新进展。许多研究使用具有高阶非线性的非线性薛定谔方程来探讨磁孤子的存在和行为 [8]、[9]。这些研究展示了弹性碰撞、 breather 形成和不稳定性驱动的孤子动力学方面的有趣结果。分析和数值结果都强调了高阶非线性对磁孤子形状和稳定性特性的影响。然而，现有文献通常考虑简化的非线性，或假设交换作用、各向异性和磁化的影响是独立的 [10]、[11]。随着关键参数（尤其是耦合常数和饱和磁化）对垂直磁各向异性的显著影响，磁设备的发展前景十分广阔。在铁磁自旋链中，局域激发可以形成长寿命的内在模式，这需要单离子各向异性和双二次交换作用，突显了高阶各向异性的重要性 [12]、[13]。此外，Satoru Hayami 还表明，中心对称磁体可以支持多个量子态。因此，推进自旋电子学的发展需要精确且高效地控制磁孤子，这一控制受到内在材料属性和外部磁振子诱导效应的驱动 [14]、[15]。这些贡献通常通过铁磁共振线宽来区分，线宽的变化反映了各向异性和磁振子相关效应。这种状态会产生具有量子能量的电磁辐射，其能量可通过磁场通过玻尔磁振子和普朗克常数进行调节 [16]、[17]。磁孤子研究已经从纯粹的理论研究发展为重要的实验和数值发展。从铁磁薄膜、纳米线和多层磁系统中获得了磁波孤子和局域磁激发的实验证据，表明了它们在自旋电子应用中的物理存在性和重要性 [18]、[19]。同时，基于L-L模型的微磁模拟提供了关于非线性自旋动力学的详细信息。这使得可以与实验进行定量比较，并在现实材料条件下验证磁孤子的稳定性和传播特性。尽管取得了这些进展，仍需要通过分析全面理解材料的不同内在参数如何影响孤子的非线性动力学 [20]、[21]。现有的分析模型要么单独考虑不同参数，要么假设交换作用、各向异性和磁化的独立影响 [10]、[11]。因此，需要一个统一的分析框架来理解这些内在材料参数如何影响孤子。在这方面，本研究通过将L-L方程简化为广义DNLS模型，发展了一种分析方法。然后利用该模型系统地研究了多参数依赖下的孤子动力学，并通过MI分析进行了研究。利用DNLS模型所做的分析预测与这些解在磁性系统中的物理相关性之间建立了更清晰的联系 [22]、[23]。

**尽管许多研究使用DNLS型方程来描述铁磁材料中发生的非线性自旋激发，但现有的大多数DNLS型方程采用传统的立方非线性和第二阶色散。本研究引入了新的术语，更准确地考虑了来自立体投影的高阶非线性色散耦合系数。本研究中生成的模型包括了 |?|2?xx 和 ?(|?|2)xx 等项，这些项在典型的DNLS方程中并未考虑。这些额外项的加入允许考虑交换刚度、各向异性能量和饱和磁化对整体非线性结构的综合影响。这与之前的研究不同，后者在为磁性系统进行DNLS简化时分别处理了交换刚度、各向异性能量和饱和磁化。因此，目前的公式使得可以研究多个物理参数函数下的孤子动力学。因此，本文的新颖之处在于提出了一个高阶DNLS方程。此外，我们证明了额外参数的加入会引发标准DNLS方程中未观察到的非线性行为。**

**研究结构**
第2节采用适当的假设和立体投影将L-L方程简化为DNLS方程，该方程控制着各向异性铁磁体中磁孤子的演化。第3节使用Hirota双线性方法构建DNLS方程的一孤子和双孤子解，并研究单轴各向异性铁磁体中的扭结-反扭结相互作用。第4节分析了铁磁自旋系统中载波波的MI，第5节总结了研究结果。虽然海森堡铁磁体中的孤子传播已有充分记录，但很少有研究从理论上分析交换刚度、各向异性常数、各向异性能量密度和饱和磁化如何影响孤子的形成和碰撞行为。现有的基于DNLS的模型通常独立研究这些参数。相比之下，本研究通过MI分析发展了一个统一的分析框架，证明这些材料常数如何共同决定了孤子的振幅、宽度、碰撞角度和相互作用弹性。

**主要贡献**
- 推导了在立体投影方法下由L-L方程产生的高阶非线性色散项的DNLS方程。
- 提出了一种连贯的方法，用于分析交换刚度、各向异性能量和饱和磁化对非线性自旋波的影响。
- 通过应用Hirota双线性方法，解析获得了DNLS方程的精确单孤子和双孤子解，从而研究了扭结-反扭结、X型和Y型孤子碰撞。
- 分析了孤子动力学随参数变化的情况，包括相互作用和MI方面。

**章节片段**
**DNLS模型中的磁孤子工程**
我们分析了一个具有交换常数和单轴各向异性常数的1D各向异性铁磁体，其中磁化动力学受到交换作用和各向异性的内在相互作用的控制，经典磁化向量P控制自旋动力学。具体来说，我们分析了一个一维铁磁系统，其中经典磁化向量P控制自旋动力学。系统的总能量来自交换作用和磁各向异性，表达为……

**磁孤子碰撞动力学的分析框架**
使用Hirota双线性方法构建的海森堡铁磁模型在描述扭结-反扭结孤子对方面的有效性得到了评估 [27]。这种形式产生了方程（5）的不同单孤子和双孤子解。对孤子碰撞的详细研究提供了对非线性波动力学的洞察，近期研究揭示了在磁性、凝聚态和等离子体介质中的多种相互作用行为 [28]、[29]。为了研究这些非线性现象……

**调制不稳定性对磁孤子演化的影响**
色散和非线性之间的平衡控制着许多物理系统，特别是在铁磁介质中，MI导致边带的指数放大 [35]、[36]。在此过程中，连续波上的小扰动会导致其分裂成一系列周期性脉冲，从而形成孤子、breather和MI模式。MI描述了连续波对小扰动的响应，并标志着局部模式形成的开始。

**结论**
我们研究了一维海森堡铁磁体中的孤子驱动磁化动力学，以阐明对自旋电子设备相关的机制。通过使用立体投影，L-L方程被映射到DNLS型模型上，从而通过Hirota的双线性方法得到了精确的分析双孤子解。我们的分析表明，包括可控的扭结-反扭结对在内的局域磁孤子表现出由内在因素决定的不同X型和Y型碰撞模式。

**作者署名声明**
R. Ravichandran：撰写——原始草稿，形式分析，概念化。
Raja Mohan O.S.：可视化，方法学。
E. Parasuraman：验证，监督。
Deepthi KVBL：软件，资源。

**利益冲突声明**
本文不存在已知的利益冲突。

**致谢**
R. Ravichandran 深表感谢印度泰米尔纳德邦金奈技术学院非线性系统中心的宝贵帮助，包括提供实验室资源、技术讨论和手稿支持——资助编号CIT/CNS/2025/RP-015。