基于微磁学新型模型的磁性斯格明子生成机制与应用探索

基于微磁学新型模型的磁性斯格明子生成机制与应用探索一、引言1.1研究背景与意义在现代凝聚态物理与自旋电子学的前沿探索中，磁性斯格明子（MagneticSkyrmion）作为一种具有独特拓扑性质的纳米磁畴结构，正逐渐成为研究的焦点。自20世纪60年代英国核物理学家TonySkyrme首次提出斯格明子的概念，用以描述非线性场理论中的非拓扑平庸类粒子解以来，这一概念在多个领域得到了广泛的应用和拓展。1989年，Bogdanov等人从理论上预测了磁性斯格明子的存在，推断在手性磁体中可能存在具有拓扑性质的自旋构型。直至2009年，Mühlbauer等人利用小角度中子散射实验首次证实了磁性斯格明子的存在，开启了对其深入研究的新篇章。磁性斯格明子具有一系列引人注目的特性，使其在自旋电子学等领域展现出巨大的应用潜力。它具有拓扑稳定性高的特点，这意味着其自旋结构在一定条件下能够保持稳定，不易受到外界干扰而发生改变，为信息的可靠存储提供了坚实的基础。其驱动电流密度低，相较于传统的磁存储单元，能够在较低的能耗下实现信息的写入、读取和传输，符合现代电子设备对低功耗的追求。磁性斯格明子还可用磁、电场和温度等多物理场进行调控，这种多场调控的特性为器件的设计和应用提供了更多的灵活性和可能性。在自旋电子学领域，磁性斯格明子有望引发一场新的技术革命。传统的金属氧化物半导体（CMOS）技术正逐渐逼近其物理极限，量子隧穿效应等问题限制了器件尺寸的进一步缩小和性能的提升。而磁性斯格明子的出现为解决这些问题提供了新的思路。例如，在信息存储方面，基于磁性斯格明子的赛道存储器具有高密度、高速度和低能耗的优势，有望成为下一代存储技术的有力竞争者。磁性斯格明子还可应用于逻辑运算、纳米振荡器、类脑运算和晶体管等领域，为实现高性能、低功耗的自旋电子学器件提供了可能。微磁学作为研究磁性材料微观磁性的重要理论方法，在磁性斯格明子的研究中发挥着不可或缺的作用。通过建立微磁学模型，能够深入理解磁性斯格明子的形成机制、稳定性条件以及其在外场作用下的动力学行为。传统的微磁学模型在解释一些复杂的磁性现象时存在一定的局限性，难以准确描述磁性斯格明子在强自旋-轨道耦合、多场相互作用等复杂条件下的行为。因此，发展新型的微磁学模型对于深入研究磁性斯格明子具有重要的理论和实际意义。本研究致力于通过建立微磁学新型模型来深入探究磁性斯格明子的产生机制及其特性。通过该模型，有望更准确地预测磁性斯格明子的形成条件和稳定性范围，为实验制备和实际应用提供理论指导。研究新型模型中不同参数对磁性斯格明子的影响，有助于揭示其内在物理规律，为进一步优化材料性能和器件设计提供科学依据。本研究对于推动磁性斯格明子在自旋电子学领域的应用，促进新型自旋电子学器件的发展具有重要的现实意义。1.2国内外研究现状磁性斯格明子的研究自其被实验证实以来，在国内外均取得了显著的进展，微磁学模型作为研究其特性和行为的重要工具，也得到了广泛的应用和深入的发展。在国外，早期的研究主要集中在利用传统微磁学模型来解释磁性斯格明子在简单体系中的形成和基本性质。Mühlbauer等人在首次证实磁性斯格明子存在的研究中，便借助了微磁学模拟来辅助理解其在MnSi材料中的形成机制，通过模拟自旋相互作用，揭示了Dzyaloshinskii-Moriya相互作用（DMI）在稳定斯格明子结构中的关键作用。此后，诸多研究基于Landau-Lifshitz-Gilbert（LLG）方程这一经典的微磁学方程，对磁性斯格明子在不同材料体系和外场条件下的行为进行了大量模拟。例如，在具有体DMI的材料体系中，研究人员通过调整微磁学模型中的参数，成功模拟出布洛赫型斯格明子的稳定存在，并分析了其在外加磁场和电流作用下的动力学行为，发现斯格明子能够在较低的电流密度下被驱动移动，这为其在自旋电子学器件中的应用提供了理论基础。随着研究的深入，国外学者开始关注如何改进微磁学模型以更准确地描述复杂的物理现象。针对磁性斯格明子在强自旋-轨道耦合体系中的行为，一些研究引入了更精确的自旋-轨道相互作用项到微磁学模型中，从而能够更好地解释实验中观察到的斯格明子的奇特动力学特性，如非对称的运动轨迹和与传统理论预测不同的响应行为。在多场耦合调控斯格明子的研究中，也发展了考虑电场、温度场等多物理场与磁场相互作用的新型微磁学模型，为实现对斯格明子的多维度精确操控提供了理论支持。在国内，磁性斯格明子的微磁学研究也呈现出蓬勃发展的态势。中国科学院物理研究所的科研团队在磁性斯格明子材料体系和微磁学模拟方面取得了一系列重要成果。他们通过微磁学模拟，在六角结构MnNiGa合金中研究双涡旋磁性斯格明子的形成条件和稳定性，发现磁各向异性和磁结构的不稳定特征在其中起到了关键作用。在二维磁性材料Fe3GeTe2中观察到布洛赫型磁性斯格明子自旋结构的研究中，微磁学模拟也被用于验证斯格明子的磁矩排列方式，进一步证实了理论分析的结果。清华大学的研究团队在利用微磁学模型研究磁性斯格明子方面也有突出贡献。赵永刚研究组利用微磁学仿真分析了受限圆盘的特殊边界条件对于斯格明子包产生和零场稳定性的机制，从能量角度解释了斯格明子包拓扑数对于圆盘尺寸的依赖性。江万军课题组借助微磁学模拟和蒂勒方程，研究了磁场梯度驱动下奈尔型斯格明子的准粒子动力学行为，揭示了斯格明子拓扑动力学特征与自旋的相似性。尽管国内外在利用微磁学模型研究磁性斯格明子方面取得了丰硕的成果，但仍存在一些不足和待解决的问题。一方面，现有的微磁学模型在处理一些极端条件下的磁性斯格明子行为时仍存在局限性，如在超高磁场或极低温度下，模型的准确性有待进一步提高。另一方面，对于复杂材料体系中多种相互作用的精确描述还不够完善，特别是在同时考虑多种微观机制的耦合效应时，模型的复杂性和计算量急剧增加，导致计算效率低下且结果的可靠性难以保证。在实验与理论模型的结合方面，虽然微磁学模拟能够为实验提供一定的指导，但在一些情况下，模拟结果与实验观测仍存在偏差，需要进一步优化模型参数和理论框架，以实现更准确的定量预测。1.3研究方法与创新点本研究综合运用理论分析、数值模拟等方法，旨在通过新型微磁学模型深入探究磁性斯格明子的产生机制与特性，在研究方法和模型构建上展现出显著的创新之处。在理论分析方面，深入剖析磁性斯格明子形成的物理机制。从微观层面出发，基于量子力学中的自旋-轨道相互作用理论，详细阐述Dzyaloshinskii-Moriya相互作用（DMI）在磁性斯格明子稳定存在中的关键作用。通过对海森伯交换相互作用、磁各向异性相互作用以及偶极相互作用等多种相互作用的综合考量，建立起一套完整的理论框架，用以解释在不同材料体系和外场条件下磁性斯格明子的形成条件和稳定性。例如，在分析具有体DMI的材料时，通过理论推导揭示了自旋-轨道耦合如何导致自旋结构的非共线排列，进而形成稳定的布洛赫型斯格明子；在研究具有界面DMI的薄膜体系时，从界面原子的电子结构出发，探讨了界面处自旋-轨道相互作用对奈尔型斯格明子形成的影响。数值模拟是本研究的重要手段之一。利用基于有限元方法的微磁学模拟软件，对磁性斯格明子在不同材料和外场下的行为进行精确模拟。在模拟过程中，将材料的微观参数，如交换刚度常数、DMI常数、磁各向异性常数等精确输入模型中，以确保模拟结果的准确性。通过改变这些参数，系统研究它们对磁性斯格明子的尺寸、稳定性、动力学行为等方面的影响。例如，在模拟外加磁场对磁性斯格明子的作用时，精确设置磁场的大小和方向，观察斯格明子在磁场下的自旋重排和运动轨迹；在研究电流驱动下斯格明子的动力学时，考虑自旋转移矩和自旋轨道矩的作用，模拟斯格明子在电流作用下的移动速度、方向以及可能出现的斯格明子霍尔效应等现象。本研究的创新点主要体现在新型微磁学模型的构建上。该模型在传统微磁学模型的基础上，引入了更精确的自旋-轨道相互作用项，以更准确地描述强自旋-轨道耦合体系中磁性斯格明子的行为。传统模型在处理这类体系时，往往忽略了自旋-轨道相互作用的高阶项以及其在不同晶体方向上的各向异性，导致对斯格明子奇特动力学特性的解释存在偏差。本研究通过引入基于第一性原理计算得到的自旋-轨道相互作用参数，能够更全面地考虑其对自旋结构的影响，从而准确解释实验中观察到的斯格明子的非对称运动轨迹、与传统理论预测不同的响应行为等现象。本模型还考虑了多物理场之间的耦合效应。在实际应用中，磁性斯格明子往往受到磁、电、温度等多种物理场的共同作用，而现有模型大多只考虑单一或两种场的作用。本研究建立的新型微磁学模型通过引入电场-磁场耦合项、温度-磁场耦合项等，能够更真实地模拟多场作用下磁性斯格明子的行为。例如，在模拟电场调控斯格明子的过程中，考虑电场通过改变材料的电子结构，进而影响自旋-轨道相互作用和磁各向异性，最终实现对斯格明子的有效调控；在研究温度对斯格明子的影响时，考虑温度引起的热涨落对自旋结构稳定性的影响，以及温度变化导致的材料磁性参数的改变对斯格明子行为的作用。本研究通过新型微磁学模型，在理论分析和数值模拟的基础上，有望为磁性斯格明子的研究提供更深入、准确的理解，为其在自旋电子学领域的应用提供更坚实的理论支持和技术指导。二、磁性斯格明子与微磁学理论基础2.1磁性斯格明子概述2.1.1定义与特性磁性斯格明子是一种具有独特拓扑性质的纳米磁畴结构，其定义源于量子场论中的拓扑孤子概念。在磁性材料中，磁性斯格明子表现为局域化的自旋涡旋结构，自旋围绕一个中心核心进行有序排列，形成类似于漩涡的构型。这种特殊的自旋排列方式赋予了磁性斯格明子一系列独特的特性，使其在自旋电子学等领域展现出巨大的应用潜力。拓扑保护是磁性斯格明子的核心特性之一。从拓扑学的角度来看，磁性斯格明子具有非零的拓扑电荷，这使得其自旋结构在一定条件下能够保持稳定，不易受到外界小的扰动而发生改变。具体而言，拓扑电荷是通过对斯格明子中自旋的空间分布进行积分计算得到的一个整数，它描述了自旋在空间中的缠绕情况。当外界扰动不足以克服斯格明子的拓扑能量壁垒时，其拓扑结构将保持不变，从而保证了信息存储的可靠性。这种拓扑保护特性使得磁性斯格明子在信息存储领域具有重要的应用价值，相较于传统的磁存储单元，能够有效降低数据丢失的风险。磁性斯格明子具有小尺寸的特性，其尺寸通常在纳米量级，一般为10-100nm。这一特性使其在实现高密度信息存储方面具有天然的优势。随着信息技术的飞速发展，对存储设备的密度要求越来越高，传统的存储技术由于物理尺寸的限制，逐渐难以满足不断增长的需求。而磁性斯格明子的小尺寸特性，使其能够在极小的空间内存储信息，为实现超高密度的存储提供了可能。例如，在未来的硬盘存储技术中，利用磁性斯格明子作为存储单元，有望大幅提高存储密度，实现存储容量的飞跃。低驱动电流是磁性斯格明子的又一显著优势。在传统的磁存储和自旋电子学器件中，畴壁的移动通常需要较大的电流密度来驱动，这不仅增加了器件的能耗，还可能导致器件发热等问题。而磁性斯格明子在较低的电流密度下就能被有效驱动，实验和理论研究表明，其驱动电流密度可比传统畴壁低几个数量级。这是因为磁性斯格明子与自旋极化电流之间存在特殊的相互作用，使得其在电流作用下能够稳定地移动。这种低驱动电流的特性使得基于磁性斯格明子的器件在运行过程中能耗大幅降低，符合现代电子设备对低功耗的追求，有助于延长设备的续航时间，减少能源消耗。除了上述特性外，磁性斯格明子还具有多场调控的能力。它可用磁、电场和温度等多物理场进行调控。在磁场作用下，磁性斯格明子的自旋结构会发生相应的变化，例如，通过改变外加磁场的大小和方向，可以实现斯格明子的产生、湮灭、移动和旋转等操作；电场则可以通过改变材料的电子结构，进而影响自旋-轨道相互作用和磁各向异性，实现对斯格明子的间接调控；温度的变化也会对磁性斯格明子的稳定性和动力学行为产生影响，当温度升高时，热涨落会增强，可能导致斯格明子的尺寸和稳定性发生改变。这种多场调控的特性为磁性斯格明子在复杂的应用场景中的灵活应用提供了可能，研究人员可以根据实际需求，选择合适的物理场对斯格明子进行精确调控，以实现特定的功能。2.1.2形成机制与应用前景磁性斯格明子的形成机制较为复杂，涉及多种磁性相互作用和物理条件。在自旋冰材料中，磁性斯格明子的形成与自旋之间的强交换相互作用以及磁各向异性密切相关。自旋冰材料具有特殊的晶体结构，其中的磁性离子形成了类似于冰的四面体网络。在这种结构中，由于自旋之间的竞争相互作用，使得磁矩倾向于形成一种特殊的“二元”排列，从而产生磁性斯格明子。这种特殊的排列方式满足了冰规则，即每个四面体顶点上的磁矩有两个指向四面体内部，两个指向外部，从而形成了稳定的斯格明子结构。在磁性多层膜体系中，Dzyaloshinskii-Moriya相互作用（DMI）在磁性斯格明子的形成中起着关键作用。DMI是一种存在于非中心对称材料中的不对称交换相互作用，它源于自旋-轨道耦合效应。在磁性多层膜中，由于界面的存在打破了空间反演对称性，使得DMI得以增强。DMI与铁磁交换相互作用、磁各向异性相互作用之间的竞争和平衡，促使自旋形成了具有手性的斯格明子结构。具体来说，DMI会导致相邻自旋之间产生一定的夹角，使得自旋排列不再是简单的共线排列，而是形成了类似于涡旋的结构，从而稳定了磁性斯格明子的存在。磁性斯格明子在自旋电子器件领域具有广阔的应用前景。在磁性随机存储器（MRAM）中，利用磁性斯格明子作为存储单元，有望实现高速、低功耗、高密度的信息存储。与传统的MRAM相比，基于斯格明子的MRAM具有更快的读写速度，这是因为斯格明子可以在低电流密度下快速移动，从而实现信息的快速写入和读取；其功耗更低，由于斯格明子的低驱动电流特性，减少了器件在运行过程中的能量消耗；还具有更高的存储密度，小尺寸的斯格明子能够在有限的空间内存储更多的信息。在自旋逻辑器件方面，磁性斯格明子可用于构建新型的逻辑门和电路。传统的逻辑器件主要基于电子的电荷属性，而自旋逻辑器件则利用电子的自旋属性来实现逻辑运算。磁性斯格明子的独特自旋结构和动力学特性使其能够实现一些传统逻辑器件难以实现的逻辑功能，例如，通过控制斯格明子的产生、湮灭和移动，可以实现布尔逻辑运算中的“与”“或”“非”等操作。基于斯格明子的自旋逻辑器件还具有低功耗、高速运行和抗干扰能力强等优点，有望为未来的计算机芯片技术带来新的突破。在量子计算领域，磁性斯格明子也展现出了潜在的应用价值。量子比特是量子计算的核心单元，需要具备良好的量子相干性和可操控性。磁性斯格明子的拓扑保护特性使其具有一定的量子特性，有望作为量子比特的候选者之一。研究人员可以通过对磁性斯格明子的自旋状态进行精确调控，实现量子比特的初始化、操作和测量等过程。例如，利用外部磁场和电场对斯格明子的自旋进行操控，实现量子比特的状态转换和逻辑运算。磁性斯格明子还可以与超导电路等其他量子系统相结合，构建混合量子计算平台，为量子计算的发展提供新的思路和方法。2.2微磁学理论基础2.2.1微磁学基本方程微磁学理论作为研究磁性材料微观磁性的重要工具，其基本方程是理解磁性现象的关键。微磁学基本方程主要包括静态微磁学方程和动态微磁学方程，它们从不同角度描述了磁性材料中磁矩的分布和演化规律。静态微磁学方程主要用于研究磁化强度矢量处于平衡状态时的空间形态分布。当温度远低于铁磁材料的居里温度时，磁化强度矢量的大小可近似等于材料的饱和磁化强度，此时磁矩形态的变化主要由局域磁矩的方向决定。在静态情况下，铁磁体的自由能主要包含多种能量项，通过对这些能量项进行分析和求解，可以得到磁矩的稳态分布。其核心思想是基于能量极小化原理，即系统总是倾向于处于能量最低的状态。在均匀外磁场H中，磁矩\vec{m}与外磁场的相互作用能（塞曼能）可表示为E_{Z}=-\mu_{0}M_{s}\vec{m}\cdot\vec{H}，其中\mu_{0}为真空磁导率，M_{s}为饱和磁化强度。当磁矩方向与外磁场方向一致时，塞曼能最低，系统处于最稳定状态。通过求解能量泛函的极小值，可以确定磁矩在平衡状态下的空间分布，从而得到磁性材料的磁畴结构等信息。动态微磁学方程则主要用于研究磁矩随时间变化的进动特性，其核心是朗道-栗弗席兹-吉尔伯特（LLG）方程：\frac{d\vec{m}}{dt}=-\gamma\vec{m}\times\vec{H}_{eff}+\frac{\alpha}{M_{s}}\vec{m}\times\frac{d\vec{m}}{dt}其中，t为时间，\gamma为旋磁比，\vec{H}_{eff}为有效磁场，它是由系统中各种能量对应的有效场叠加而成，M_{s}为饱和磁化强度，\alpha为材料阻尼系数，取值一般在0.01-0.05。方程右侧第一项描述的是磁矩绕有效磁场的进动行为，体现了磁矩在有效磁场作用下的旋转特性；第二项为阻尼项，描述了体系能量趋向最低能量时的速度，它使得磁矩的进动逐渐衰减，最终趋向于能量最低的稳定状态。在磁性材料受到外加脉冲磁场的作用时，磁矩会在有效磁场的作用下发生进动。当脉冲磁场施加的瞬间，磁矩会迅速偏离其初始平衡方向，开始绕有效磁场进动。随着时间的推移，由于阻尼项的存在，磁矩的进动幅度会逐渐减小，最终趋向于与有效磁场方向一致，达到新的平衡状态。在自旋电子学器件中，电流的注入会产生自旋转移矩效应，这可以通过在LLG方程中添加相应的项来描述。自旋转移矩效应会对磁矩的运动产生影响，从而实现对磁性器件的电学调控。LLG方程还可以写成朗道-栗弗席兹（LL）方程的形式，主要差别在于阻尼项的具体形式，两者在本质上是等价的，但在不同的应用场景中，根据计算的便利性和物理问题的特点，可以选择使用不同形式的方程。在处理一些需要精确描述磁矩进动过程的问题时，LLG方程由于其对阻尼项的更细致描述，可能更适合；而在一些对计算效率要求较高，且对阻尼项精度要求相对较低的情况下，LL方程可能更具优势。2.2.2磁性材料中的磁能磁性材料中的磁能是理解磁性斯格明子形成和行为的关键因素，不同形式的磁能相互作用决定了磁矩的分布和稳定性。交换能是磁性材料中最基本的能量形式之一，它源于电子间的交换相互作用，是一种量子力学效应。在铁磁材料中，相邻原子的电子自旋倾向于平行排列，以降低交换能。交换能的表达式为E_{ex}=-2A\sum_{i,j}\vec{S}_{i}\cdot\vec{S}_{j}，其中A为交换积分常数，\vec{S}_{i}和\vec{S}_{j}分别为相邻原子的自旋矢量。交换能倾向于使自旋保持平行排列，从而维持材料的铁磁特性。在磁性斯格明子的形成过程中，交换能与其他磁能相互竞争，影响着斯格明子的尺寸和稳定性。当交换能较强时，自旋更倾向于平行排列，不利于斯格明子的形成；反之，较弱的交换能则为斯格明子的形成提供了一定的条件。磁晶各向异性能是由于晶体结构的各向异性导致的磁性能差异所产生的能量。在晶体中，不同晶轴方向上的磁矩取向具有不同的能量，磁矩倾向于沿能量较低的易磁化轴方向排列。对于立方晶系的磁性材料，磁晶各向异性能可表示为E_{k}=K_{1}(\alpha_{1}^{2}\alpha_{2}^{2}+\alpha_{2}^{2}\alpha_{3}^{2}+\alpha_{3}^{2}\alpha_{1}^{2})+K_{2}\alpha_{1}^{2}\alpha_{2}^{2}\alpha_{3}^{2}，其中K_{1}和K_{2}为磁晶各向异性常数，\alpha_{1}、\alpha_{2}、\alpha_{3}为磁矩方向与晶轴的方向余弦。磁晶各向异性能对磁性斯格明子的影响主要体现在其对斯格明子磁矩取向的约束上，使得斯格明子的磁矩在满足其他条件的同时，还需考虑磁晶各向异性的作用，从而影响斯格明子的稳定性和拓扑结构。退磁能是由于磁性材料的形状和磁化不均匀性导致的磁矩间相互作用产生的能量。当磁性材料被磁化时，会在其表面和内部产生磁极，这些磁极会产生退磁场，退磁能就是磁矩在退磁场中的能量。退磁能的表达式为E_{d}=-\frac{1}{2}\mu_{0}\int\vec{M}\cdot\vec{H}_{d}dV，其中\vec{H}_{d}为退磁场，\vec{M}为磁化强度。退磁能倾向于使磁矩分布均匀，以降低能量。在磁性斯格明子的研究中，退磁能会影响斯格明子的边界条件和相互作用，对于多斯格明子体系，退磁能会导致斯格明子之间的相互排斥或吸引，从而影响它们的排列和稳定性。Dzyaloshinskii-Moriya相互作用能（DMI）是一种存在于非中心对称材料中的不对称交换相互作用能，它源于自旋-轨道耦合效应。DMI在相邻自旋之间产生不对称的交换相互作用，导致自旋形成非共线的手性结构，这对于磁性斯格明子的形成和稳定起着至关重要的作用。其表达式为E_{DMI}=D\sum_{i,j}\vec{S}_{i}\times\vec{S}_{j}\cdot\vec{n}_{ij}，其中D为DMI常数，\vec{n}_{ij}为相邻自旋连线方向的单位矢量。在具有DMI的材料中，DMI与交换能、磁晶各向异性能等相互竞争和平衡，促使自旋形成具有手性的斯格明子结构。在某些磁性多层膜中，由于界面的存在打破了空间反演对称性，使得DMI增强，从而有利于奈尔型斯格明子的形成。2.2.3微磁学模拟方法微磁学模拟是研究磁性材料微观磁性的重要手段，它能够通过数值计算的方法求解微磁学方程，从而深入了解磁性斯格明子等磁结构的特性和行为。有限差分法和有限元法是两种常用的微磁学模拟方法，它们在求解微磁学方程中各有优势和应用场景。有限差分法是将待模拟的整个磁性样品划分为一个个小的规则的微磁学单元，通常为立方体结构。在每个单元内，假设磁矩是均匀分布的，通过对微磁学方程进行离散化处理，将连续的空间和时间变量转化为离散的数值，从而进行数值求解。对于描述磁矩随时间变化的朗道-栗弗席兹-吉尔伯特（LLG）方程，在有限差分法中，会将时间划分为一系列小的时间步长\Deltat，空间划分为小的网格单元。在每个时间步长内，根据前一时刻各单元的磁矩和有效磁场，通过有限差分公式计算下一时刻各单元磁矩的变化，从而逐步模拟磁矩的动态演化过程。有限差分法的优点在于计算速度相对较快，因为其单元划分规则，计算过程相对简单，易于实现并行计算，能够在较短时间内完成大规模的模拟计算。在研究磁性斯格明子在均匀磁场下的动态行为时，利用有限差分法可以快速得到不同时刻斯格明子的位置、自旋结构等信息，为研究其动力学特性提供了高效的手段。由于其对曲面边界的处理能力较弱，在模拟具有复杂几何形状的磁性材料时，可能会引入较大的误差。有限元法则是将磁性材料划分为一系列形状更为灵活的单元，如三角形、四面体等，能够更好地拟合复杂的几何形状。在有限元法中，通过对每个单元内的磁矩进行插值逼近，将微磁学方程转化为矩阵形式的代数方程组进行求解。在处理具有曲面边界的磁性薄膜时，有限元法可以根据薄膜的实际形状精确划分单元，更准确地描述边界条件，从而提高模拟的精度。有限元法还能够更好地处理材料参数的空间变化，对于研究磁性材料中不同区域具有不同磁性能的情况具有优势。在模拟磁性多层膜时，不同层的材料参数如交换刚度常数、磁各向异性常数等可能不同，有限元法可以方便地考虑这些差异，更真实地模拟多层膜中的磁现象。其计算静磁相互作用（退磁场）时耗时较长，因为在计算退磁场时，需要对整个计算区域进行积分运算，计算量较大。在实际应用中，选择有限差分法还是有限元法，需要根据具体的研究对象和问题来决定。如果研究对象的几何形状简单，对计算速度要求较高，有限差分法可能是更好的选择；而对于几何形状复杂、对模拟精度要求较高的情况，有限元法更能满足需求。也可以结合两种方法的优点，采用混合方法进行模拟，以提高模拟的效率和精度。三、微磁学新型模型构建3.1新型模型的设计思路传统微磁学模型在描述磁性斯格明子的产生和特性时存在一定的局限性。在处理强自旋-轨道耦合体系时，传统模型往往难以准确描述自旋-轨道相互作用对斯格明子自旋结构的影响。由于其对自旋-轨道相互作用的描述过于简化，无法解释实验中观察到的斯格明子在强自旋-轨道耦合下的奇特动力学行为，如非对称的运动轨迹和与传统理论预测不同的响应行为。在多场相互作用的情况下，传统模型也难以全面考虑磁、电、温度等多种物理场对斯格明子的综合影响，导致对斯格明子在复杂外场条件下的行为预测出现偏差。为了克服这些局限性，本研究设计的新型微磁学模型引入了更精确的自旋-轨道相互作用项。基于第一性原理计算，能够得到更准确的自旋-轨道相互作用参数，这些参数能够全面考虑自旋-轨道相互作用的高阶项以及其在不同晶体方向上的各向异性。通过将这些参数引入模型中，能够更真实地描述自旋-轨道相互作用对自旋结构的影响，从而准确解释斯格明子在强自旋-轨道耦合体系中的奇特动力学特性。新型模型还考虑了多物理场之间的耦合效应。在实际应用中，磁性斯格明子通常会受到磁、电、温度等多种物理场的共同作用。为了更真实地模拟这种复杂的物理情况，新型模型引入了电场-磁场耦合项、温度-磁场耦合项等。电场-磁场耦合项考虑了电场通过改变材料的电子结构，进而影响自旋-轨道相互作用和磁各向异性，最终实现对斯格明子的有效调控。在一些具有电场调控特性的磁性材料中，电场的变化会导致材料内部电子云分布的改变，从而影响自旋-轨道相互作用的强度和方向，进一步影响斯格明子的自旋结构和稳定性。通过在模型中引入电场-磁场耦合项，可以准确描述这种电场对斯格明子的调控机制。温度-磁场耦合项则考虑了温度引起的热涨落对自旋结构稳定性的影响，以及温度变化导致的材料磁性参数的改变对斯格明子行为的作用。当温度升高时，热涨落会增强，可能导致斯格明子的尺寸和稳定性发生改变。温度的变化还会引起材料的磁性参数，如交换刚度常数、磁各向异性常数等发生变化，从而影响斯格明子的形成和动力学行为。新型模型通过引入温度-磁场耦合项，能够全面考虑这些温度相关的因素对斯格明子的影响，更准确地预测斯格明子在不同温度条件下的行为。新型模型在设计时还充分考虑了计算效率和可扩展性。采用了先进的数值计算方法和算法优化技术，以降低计算复杂度，提高计算速度，使得模型能够在合理的时间内完成大规模的模拟计算。在处理复杂材料体系和多物理场相互作用时，模型的可扩展性能够保证其能够方便地添加新的物理项和参数，以适应不同的研究需求。通过模块化的设计思路，将模型中的不同物理过程和相互作用项进行独立封装，使得在需要添加新的物理机制时，只需对相应的模块进行修改和扩展，而不会影响整个模型的稳定性和计算效率。3.2模型关键参数与特性分析新型微磁学模型中包含多个关键参数，这些参数对模型特性和磁性斯格明子的产生有着至关重要的影响。交换常数（A）是描述磁性材料中相邻原子自旋间交换相互作用强度的关键参数。在海森伯模型中，交换能E_{ex}=-2A\sum_{i,j}\vec{S}_{i}\cdot\vec{S}_{j}，其中\vec{S}_{i}和\vec{S}_{j}分别为相邻原子的自旋矢量。交换常数A的大小直接影响着自旋间的相互作用，A越大，自旋越倾向于平行排列，有利于维持铁磁态；A越小，自旋间的相互作用减弱，为非共线自旋结构的形成提供了可能，这对于磁性斯格明子的产生至关重要。在一些具有Dzyaloshinskii-Moriya相互作用（DMI）的材料中，当交换常数A与DMI常数达到一定的平衡时，自旋会形成具有手性的斯格明子结构。通过数值模拟研究发现，在特定的材料体系中，当交换常数A从1\times10^{-11}J/m增加到2\times10^{-11}J/m时，磁性斯格明子的尺寸会逐渐减小，这是因为较强的交换相互作用使得自旋更紧密地排列，从而压缩了斯格明子的空间范围。各向异性常数（K）体现了磁性材料在不同方向上磁性能的差异，磁晶各向异性能E_{k}=K_{1}(\alpha_{1}^{2}\alpha_{2}^{2}+\alpha_{2}^{2}\alpha_{3}^{2}+\alpha_{3}^{2}\alpha_{1}^{2})+K_{2}\alpha_{1}^{2}\alpha_{2}^{2}\alpha_{3}^{2}（以立方晶系为例），其中K_{1}和K_{2}为磁晶各向异性常数，\alpha_{1}、\alpha_{2}、\alpha_{3}为磁矩方向与晶轴的方向余弦。各向异性常数决定了磁矩的易磁化方向，对磁性斯格明子的稳定性和磁矩取向有着重要影响。在具有单轴各向异性的材料中，磁矩倾向于沿单轴方向排列，当存在其他相互作用（如DMI）时，会促使磁矩在满足各向异性的基础上形成斯格明子结构。若各向异性常数K增大，磁矩沿易磁化方向排列的趋势增强，这可能会提高磁性斯格明子的稳定性，但同时也可能限制斯格明子的形成和运动，因为磁矩的重新取向需要克服更大的能量壁垒。通过改变各向异性常数K的值，观察到在K较小时，磁性斯格明子更容易在外加磁场的作用下发生变形和移动；而当K增大到一定程度后，斯格明子变得更加稳定，对外加磁场的响应也变得更加迟缓。Dzyaloshinskii-Moriya相互作用常数（D）是新型微磁学模型中另一个关键参数，它描述了DMI的强度。DMI是一种存在于非中心对称材料中的不对称交换相互作用，其能量表达式为E_{DMI}=D\sum_{i,j}\vec{S}_{i}\times\vec{S}_{j}\cdot\vec{n}_{ij}，其中\vec{n}_{ij}为相邻自旋连线方向的单位矢量。DMI常数D的正负和大小决定了自旋的手性和非共线排列方式，对于磁性斯格明子的形成起着核心作用。在具有正DMI常数的材料中，自旋会形成特定手性的斯格明子结构，如布洛赫型斯格明子；而在负DMI常数的情况下，可能形成奈尔型斯格明子。当DMI常数D增大时，自旋的非共线排列趋势增强，更容易形成稳定的斯格明子结构，且斯格明子的尺寸会随着D的增大而减小。通过模拟不同DMI常数下磁性斯格明子的形成过程，发现当D从1\times10^{-3}J/m^{2}增加到2\times10^{-3}J/m^{2}时，在相同的模拟条件下，斯格明子能够更快地形成，且其边界更加清晰，表明其稳定性增强。饱和磁化强度（M_{s}）反映了磁性材料在饱和磁化状态下的磁化程度，在塞曼能E_{Z}=-\mu_{0}M_{s}\vec{m}\cdot\vec{H}中，M_{s}直接影响磁矩与外磁场的相互作用能量。饱和磁化强度M_{s}对磁性斯格明子的产生和特性也有显著影响，它决定了磁矩的大小，进而影响斯格明子的能量和稳定性。在相同的外场和其他参数条件下，M_{s}越大，磁矩与外磁场的相互作用越强，这可能会影响斯格明子在磁场中的行为，如运动速度和方向。对于一些基于磁性斯格明子的器件应用，M_{s}的大小还会影响器件的性能，如在磁性随机存储器中，M_{s}会影响信息的写入和读取速度。通过模拟发现，当饱和磁化强度M_{s}从1\times10^{5}A/m增加到2\times10^{5}A/m时，磁性斯格明子在相同外加磁场下的运动速度明显加快，这是因为更大的M_{s}使得磁矩与磁场的相互作用增强，从而获得更大的驱动力。3.3模型验证与对比分析为了验证新型微磁学模型的准确性和可靠性，将模拟结果与已有实验数据进行了详细对比。在磁性斯格明子的尺寸方面，实验上通过洛伦兹透射电子显微镜等技术，对具有Dzyaloshinskii-Moriya相互作用（DMI）的磁性薄膜中的斯格明子尺寸进行了精确测量。将新型模型的模拟结果与这些实验数据对比发现，在相同的材料参数和外场条件下，新型模型预测的斯格明子尺寸与实验测量值高度吻合。在特定的[Ta/CoFeB/MgO]多层膜体系中，实验测得在一定的DMI常数和外磁场下，磁性斯格明子的直径约为50nm，新型模型的模拟结果为48nm，误差在合理范围内，表明新型模型能够准确预测斯格明子的尺寸。在磁性斯格明子的稳定性研究中，实验通过测量斯格明子在不同温度和外磁场下的寿命来评估其稳定性。新型模型在考虑了温度-磁场耦合效应后，对斯格明子的稳定性进行了模拟。结果显示，模拟得到的斯格明子在不同条件下的稳定性变化趋势与实验结果一致。当温度升高时，实验中观察到斯格明子的稳定性逐渐下降，新型模型也准确地模拟出了这一趋势，通过计算斯格明子的能量变化，预测了其在高温下可能发生的湮灭现象，与实验中观察到的结果相符。还将新型模型的模拟结果与传统微磁学模型进行了对比分析，以凸显新型模型的优势。在处理强自旋-轨道耦合体系时，传统模型由于对自旋-轨道相互作用的描述不够精确，无法准确预测磁性斯格明子的动力学行为。在模拟电流驱动下斯格明子的运动时，传统模型预测斯格明子将沿着电流方向直线运动，而实际实验中观察到斯格明子会发生一定角度的偏转，即出现斯格明子霍尔效应。新型模型由于引入了更精确的自旋-轨道相互作用项，能够准确地模拟出斯格明子霍尔效应，预测的斯格明子运动轨迹与实验结果一致。在多场相互作用的情况下，传统模型同样存在局限性。当同时考虑磁场和电场对磁性斯格明子的作用时，传统模型难以全面考虑电场通过改变材料电子结构对斯格明子的间接影响。新型模型通过引入电场-磁场耦合项，能够准确描述电场对斯格明子的调控机制。在模拟电场调控斯格明子的实验中，新型模型能够准确预测斯格明子在电场作用下的自旋结构变化和运动行为，而传统模型的预测结果与实验相差较大。通过与已有实验数据和传统微磁学模型的对比分析，充分验证了新型微磁学模型在描述磁性斯格明子的产生、特性和动力学行为方面的准确性和优越性，为进一步深入研究磁性斯格明子提供了可靠的理论工具。四、基于新型模型的磁性斯格明子产生过程4.1模拟条件设定在利用新型微磁学模型模拟磁性斯格明子的产生过程时，首先确定模拟体系的几何结构为二维正方形薄膜，其边长设定为L=500nm，厚度为t=1nm。选择二维正方形薄膜结构是因为在实际的自旋电子学器件应用中，薄膜材料是常用的体系，且正方形结构便于进行数值计算和边界条件的设定，能够更直观地研究磁性斯格明子在二维平面内的形成和演化过程。模拟体系选用的材料为具有典型Dzyaloshinskii-Moriya相互作用（DMI）的[Ta/CoFeB/MgO]多层膜体系。这种材料体系在实验中已被广泛研究，且具有明确的磁性参数，有利于与实验结果进行对比验证。其中，交换常数A=1.5\times10^{-11}J/m，这一数值是基于该材料体系的大量实验测量和理论计算确定的，它反映了材料中相邻原子自旋间交换相互作用的强度。各向异性常数K=1\times10^{6}J/m^{3}，体现了材料在不同方向上磁性能的差异，决定了磁矩的易磁化方向。Dzyaloshinskii-Moriya相互作用常数D=3\times10^{-3}J/m^{2}，描述了DMI的强度，对磁性斯格明子的形成起着关键作用。饱和磁化强度M_{s}=1\times10^{6}A/m，反映了材料在饱和磁化状态下的磁化程度。外部条件的设定对于磁性斯格明子的产生至关重要。首先考虑外加磁场，初始外加磁场设置为沿薄膜法线方向（z方向），大小为H_{0}=100mT。选择这一磁场方向和大小是因为在前期的研究中发现，垂直于薄膜平面的磁场能够有效地诱导自旋的重新排列，促进磁性斯格明子的形成。在该磁场强度下，能够在一定程度上克服材料内部的能量壁垒，使得自旋结构发生变化，有利于观察到磁性斯格明子的产生过程。当磁场强度较小时，可能无法提供足够的能量来促使自旋结构的转变，难以形成稳定的斯格明子；而当磁场强度过大时，可能会导致自旋直接被磁化到与磁场方向一致，同样不利于斯格明子的形成。电流条件的设定为注入沿x方向的自旋极化电流，电流密度为j=1\times10^{11}A/m^{2}。自旋极化电流在磁性斯格明子的产生和动力学研究中起着重要作用，它可以通过自旋转移矩效应与自旋相互作用，从而驱动自旋的运动和重新排列。选择这一电流密度是因为在相关的实验和理论研究中发现，该电流密度既能够有效地驱动自旋的变化，又不会因为电流过大而导致材料的过热或其他不良效应。在较低的电流密度下，自旋转移矩效应较弱，可能无法有效地驱动自旋形成斯格明子结构；而过高的电流密度则可能会引起材料的焦耳热效应，影响材料的磁性和斯格明子的稳定性。在模拟过程中，还考虑了温度因素，将体系温度设定为T=300K，接近室温条件。温度对磁性材料的磁性有显著影响，在室温下进行模拟更符合实际应用场景。在较高温度下，热涨落会增强，可能导致自旋结构的不稳定，不利于磁性斯格明子的形成和稳定存在；而在低温下，虽然自旋结构相对稳定，但与实际应用条件相差较大。通过将温度设定为300K，能够在更接近实际应用的条件下研究磁性斯格明子的产生和特性。4.2斯格明子产生的动态过程模拟在模拟过程的初始阶段，体系处于均匀的初始磁矩状态，磁矩方向与外加磁场方向一致，均沿薄膜法线方向（z方向）。此时，体系的能量主要由塞曼能和交换能构成，由于磁矩方向一致，交换能处于相对较低的状态，而塞曼能也达到了一定的平衡。随着模拟时间的推进，注入沿x方向的自旋极化电流，电流密度为j=1\times10^{11}A/m^{2}。自旋极化电流通过自旋转移矩效应与自旋相互作用，使得部分区域的磁矩开始发生偏离。在自旋转移矩的作用下，磁矩受到一个与电流方向相关的力矩，这个力矩促使磁矩开始绕其平衡位置进动，从而导致部分区域的磁矩方向发生改变。当模拟时间达到t_1=1ns时，磁矩分布出现了明显的变化，在薄膜的局部区域开始出现自旋的非共线排列。由于Dzyaloshinskii-Moriya相互作用（DMI）的存在，这种非共线排列逐渐呈现出手性特征，形成了一些微小的自旋漩涡结构。DMI使得相邻自旋之间产生不对称的交换相互作用，导致自旋排列不再是简单的共线排列，而是形成了具有一定手性的结构。这些微小的自旋漩涡结构是磁性斯格明子的雏形，它们的出现标志着磁性斯格明子开始逐渐形成。随着模拟时间进一步增加到t_2=3ns，这些自旋漩涡结构逐渐演化和合并。在这个过程中，自旋之间的相互作用不断调整，能量也在不断地重新分布。一些相邻的自旋漩涡结构由于相互吸引而合并，使得斯格明子的尺寸逐渐增大。同时，由于体系中存在多种相互作用的竞争，如交换能、DMI和磁各向异性能等，斯格明子的结构也在不断地优化，以达到能量最低的稳定状态。在这个阶段，斯格明子的边界逐渐清晰，自旋的排列更加有序，形成了较为稳定的涡旋状结构。当模拟时间达到t_3=5ns时，磁性斯格明子基本形成稳定状态。此时，斯格明子的拓扑结构已经稳定，其自旋围绕中心核心呈规则的涡旋状排列，拓扑电荷保持稳定。在稳定状态下，斯格明子的能量达到了一个相对较低的平衡值，各种相互作用之间达到了平衡。交换能、DMI、磁各向异性能和塞曼能等相互制约，使得斯格明子的结构能够在一定条件下保持稳定。外加磁场和电流对斯格明子的稳定性也有一定的影响。当外加磁场发生变化时，斯格明子的磁矩会受到磁场力的作用，可能导致其自旋结构发生微小的调整，但由于其拓扑保护特性，在一定范围内，斯格明子的拓扑结构仍然能够保持稳定。电流的持续作用也会对斯格明子产生影响，在稳定状态下，电流驱动斯格明子以一定的速度沿着特定方向移动，其移动速度和方向与电流密度、材料参数等因素密切相关。4.3影响斯格明子产生的因素分析外加磁场作为影响磁性斯格明子产生的关键因素之一，其强度和方向对斯格明子的产生效率、稳定性和尺寸有着显著的影响。在强度方面，当外加磁场强度较小时，如低于一定的阈值，材料内部的自旋结构相对稳定，难以产生磁性斯格明子。这是因为较小的磁场无法提供足够的能量来克服自旋之间的相互作用能垒，使得自旋难以形成具有特殊拓扑结构的斯格明子。随着外加磁场强度的逐渐增加，达到一定强度后，磁场能够有效地破坏原有的自旋排列，促使自旋发生重排，从而为斯格明子的产生创造条件。当磁场强度继续增加时，斯格明子的稳定性会发生变化。较强的磁场会对斯格明子的自旋结构产生较大的作用力，可能导致斯格明子的尺寸减小，甚至在磁场强度超过一定值时，斯格明子会被破坏，其拓扑结构消失，自旋被磁化到与磁场方向一致。外加磁场的方向也对斯格明子的产生产生重要影响。在不同的磁场方向下，自旋受到的力矩方向和大小不同，从而影响斯格明子的形成和特性。当外加磁场方向与材料的易磁化轴方向一致时，自旋更容易在外加磁场的作用下发生重排，有利于斯格明子的产生；而当磁场方向与易磁化轴方向夹角较大时，自旋需要克服更大的各向异性能垒才能发生重排，这可能会降低斯格明子的产生效率，甚至导致无法形成稳定的斯格明子。在一些具有面内各向异性的磁性薄膜中，当外加磁场垂直于薄膜平面时，能够有效地诱导自旋在平面内发生重排，形成稳定的斯格明子结构；而当磁场平行于薄膜平面时，斯格明子的形成则相对困难，且形成的斯格明子结构可能不稳定。电流密度对磁性斯格明子的产生和运动也有着重要的影响。在自旋极化电流的作用下，电流通过自旋转移矩效应与自旋相互作用，从而影响斯格明子的产生和动力学行为。当电流密度较低时，自旋转移矩效应较弱，难以有效地驱动自旋形成斯格明子结构。随着电流密度的增加，自旋转移矩效应增强，能够提供足够的能量来驱动自旋的运动和重排，从而促进斯格明子的产生。在一定的电流密度范围内，斯格明子能够在外加电流的驱动下稳定地运动，其运动速度与电流密度成正比。当电流密度过高时，会产生焦耳热效应，导致材料温度升高，这可能会破坏斯格明子的稳定性，甚至使其湮灭。过高的电流密度还可能导致材料的损伤，影响其磁性和斯格明子的产生。材料参数是影响磁性斯格明子产生的内在因素，不同的材料参数会导致斯格明子的产生和特性存在差异。交换常数（A）反映了材料中相邻原子自旋间交换相互作用的强度，它对斯格明子的尺寸和稳定性有着重要影响。当交换常数较大时，自旋间的相互作用较强，自旋更倾向于平行排列，这会使得斯格明子的尺寸减小，稳定性增强；而交换常数较小时，自旋间的相互作用较弱，有利于斯格明子的形成，但斯格明子的稳定性可能会降低。各向异性常数（K）体现了材料在不同方向上磁性能的差异，决定了磁矩的易磁化方向。较大的各向异性常数会使磁矩更倾向于沿易磁化方向排列，这可能会影响斯格明子的形成和运动，因为磁矩的重新取向需要克服更大的能量壁垒；而较小的各向异性常数则相对有利于斯格明子的产生和运动。Dzyaloshinskii-Moriya相互作用常数（D）是决定磁性斯格明子形成的关键参数之一。DMI是一种存在于非中心对称材料中的不对称交换相互作用，其常数D的正负和大小决定了自旋的手性和非共线排列方式。在具有正DMI常数的材料中，自旋会形成特定手性的斯格明子结构，如布洛赫型斯格明子；而在负DMI常数的情况下，可能形成奈尔型斯格明子。当DMI常数D增大时，自旋的非共线排列趋势增强，更容易形成稳定的斯格明子结构，且斯格明子的尺寸会随着D的增大而减小。饱和磁化强度（M_{s}）反映了材料在饱和磁化状态下的磁化程度，它对斯格明子的能量和稳定性也有显著影响。较大的饱和磁化强度会使磁矩与外磁场的相互作用增强，从而影响斯格明子在磁场中的行为，如运动速度和方向。五、案例分析5.1案例一：[Ta/CoFeB/MgO]多层膜体系中斯格明子的产生[Ta/CoFeB/MgO]多层膜体系是一种在磁性斯格明子研究中被广泛关注的材料体系，其独特的结构和优异的磁性特性为斯格明子的产生提供了有利条件。该多层膜体系由Ta、CoFeB和MgO等材料交替堆叠而成，各层材料在体系中发挥着不同的作用。Ta层作为重金属层，具有较强的自旋-轨道耦合效应，能够在与CoFeB层的界面处诱导出Dzyaloshinskii-Moriya相互作用（DMI）。这种DMI对于磁性斯格明子的形成至关重要，它打破了自旋排列的对称性，促使自旋形成具有手性的结构，从而为斯格明子的产生创造了条件。CoFeB层是磁性层，其主要作用是提供磁性，其中的铁磁交换相互作用是维持自旋有序排列的重要因素。在没有其他干扰的情况下，CoFeB层中的自旋倾向于平行排列，形成铁磁态。MgO层则作为隔离层，其作用是分隔不同的磁性层，减少磁性层之间的磁耦合，同时也能够影响体系的磁各向异性。通过精确控制MgO层的厚度，可以调节体系的垂直磁各向异性，进而影响磁性斯格明子的稳定性和尺寸。利用新型微磁学模型对[Ta/CoFeB/MgO]多层膜体系中斯格明子的产生进行模拟，得到了一系列有价值的结果。在模拟过程中，设置体系的几何结构为二维正方形薄膜，边长为500nm，厚度为1nm。材料参数设置为交换常数A=1.5\times10^{-11}J/m，各向异性常数K=1\times10^{6}J/m^{3}，DMI常数D=3\times10^{-3}J/m^{2}，饱和磁化强度M_{s}=1\times10^{6}A/m。外部条件设定为初始外加磁场沿薄膜法线方向（z方向），大小为H_{0}=100mT，注入沿x方向的自旋极化电流，电流密度为j=1\times10^{11}A/m^{2}，体系温度为T=300K。模拟结果显示，在初始阶段，体系处于均匀的初始磁矩状态，磁矩方向与外加磁场方向一致。随着自旋极化电流的注入，在自旋转移矩效应和DMI的共同作用下，体系中的磁矩开始发生变化。在1ns时，部分区域的磁矩出现非共线排列，形成了一些微小的自旋漩涡结构，这些结构是斯格明子的雏形。随着时间的推进，到3ns时，自旋漩涡结构逐渐演化和合并，斯格明子的尺寸不断增大，边界逐渐清晰。当模拟时间达到5ns时，磁性斯格明子基本形成稳定状态，其自旋围绕中心核心呈规则的涡旋状排列，拓扑电荷保持稳定。为了验证模拟结果的准确性，将其与相关实验结果进行对比。实验上通过洛伦兹透射电子显微镜（LTEM）对[Ta/CoFeB/MgO]多层膜体系中的磁性斯格明子进行观察。实验结果显示，在类似的外加磁场和电流条件下，能够观察到稳定的磁性斯格明子，其尺寸和结构与模拟结果具有较高的契合度。实验测得的斯格明子直径约为40-50nm，而模拟得到的斯格明子直径在45nm左右，误差在合理范围内。在斯格明子的自旋结构方面，实验观察到的自旋排列方式与模拟结果中的涡旋状结构一致，进一步验证了新型微磁学模型在预测[Ta/CoFeB/MgO]多层膜体系中斯格明子产生和特性方面的准确性。5.2案例二：二维Janus磁体MnSTe中斯格明子的调控二维Janus磁体MnSTe由于其独特的原子结构和电子特性，为磁性斯格明子的研究提供了一个新颖的平台。在这种材料中，原子的非对称排列导致了较强的自旋-轨道耦合效应，进而产生了明显的Dzyaloshinskii-Moriya相互作用（DMI），这为斯格明子的形成和调控创造了有利条件。通过新型微磁学模型的模拟，研究发现改变外部磁场是调控二维Janus磁体MnSTe中斯格明子的有效方法之一。当施加沿垂直于薄膜平面方向的外加磁场时，随着磁场强度的逐渐增加，斯格明子的直径呈现出规律性的变化。模拟结果表明，在一定的磁场强度范围内，斯格明子的直径会随着磁场强度的增大而逐渐减小。这是因为外加磁场的增强会对斯格明子的自旋结构产生更大的作用力，使得自旋更紧密地排列，从而压缩了斯格明子的空间范围。当磁场强度达到一定阈值时，斯格明子的拓扑荷甚至会发生改变。在特定的模拟条件下，当磁场强度从0逐渐增加到0.5T时，观察到斯格明子的直径从初始的约20nm减小到10nm左右，并且在磁场强度达到0.4T时，斯格明子的拓扑荷从+1转变为-1，这表明通过改变外加磁场的强度，可以实现对斯格明子拓扑性质的有效调控。改变材料的参数也是调控斯格明子的重要手段。在二维Janus磁体MnSTe中，调整交换常数（A）会对斯格明子的稳定性和尺寸产生显著影响。当交换常数增大时，自旋间的相互作用增强，自旋更倾向于平行排列，这会使得斯格明子的尺寸减小，稳定性增强。通过模拟将交换常数A从1\times10^{-11}J/m增加到2\times10^{-11}J/m，发现斯格明子的尺寸明显减小，同时在受到外界干扰时，斯格明子的结构更加稳定，不易发生变形和湮灭。调整各向异性常数（K）也会对斯格明子的特性产生影响。各向异性常数决定了磁矩的易磁化方向，当各向异性常数增大时，磁矩沿易磁化方向排列的趋势增强，这可能会提高斯格明子的稳定性，但同时也可能限制斯格明子的形成和运动。在模拟中，将各向异性常数K从5\times10^{5}J/m^{3}增加到1\times10^{6}J/m^{3}，发现斯格明子在受到外加磁场或电流驱动时，其运动速度明显降低，且在相同的外界条件下，斯格明子更难发生变形和移动，表明其稳定性得到了提高。通过新型微磁学模型对二维Janus磁体MnSTe中斯格明子的调控研究，深入了解了外部条件和材料参数对斯格明子的影响机制，为进一步实现对斯格明子的精确调控和在自旋电子学器件中的应用提供了理论基础。5.3案例对比与经验总结在[Ta/CoFeB/MgO]多层膜体系和二维Janus磁体MnSTe这两个案例中，磁性斯格明子的产生和调控存在着诸多异同点。从相同点来看，Dzyaloshinskii-Moriya相互作用（DMI）在两个案例中都对磁性斯格明子的产生起到了关键作用。在[Ta/CoFeB/MgO]多层膜体系中，Ta层与CoFeB层的界面处诱导出的DMI促使自旋形成手性结构，从而产生斯格明子；在二维Janus磁体MnSTe中，原子的非对称排列导致较强的自旋-轨道耦合效应，进而产生明显的DMI，为斯格明子的形成创造了条件。外加磁场在两个案例中都是调控斯格明子的重要手段。通过改变外加磁场的强度和方向，可以改变斯格明子的尺寸、稳定性和拓扑荷等特性。在[Ta/CoFeB/MgO]多层膜体系中，合适强度和方向的外加磁场能够诱导自旋重排，促进斯格明子的产生和稳定；在二维Janus磁体MnSTe中，随着外加磁场强度的变化，斯格明子的直径会发生改变，甚至拓扑荷也会发生转变。两个案例也存在一些不同之处。在材料体系方面，[Ta/CoFeB/MgO]多层膜体系是一种多层膜结构，各层材料的特性和相互作用共同影响着斯格明子的产生和特性；而二维Janus磁体MnSTe是一种二维材料，其原子结构和电子特性具有独特性，与多层膜体系有着明显的区别。在斯格明子的调控方式上，除了外加磁场外，[Ta/CoFeB/MgO]多层膜体系中还可以通过注入自旋极化电流来调控斯格明子，电流通过自旋转移矩效应与自旋相互作用，影响斯格明子的产生和运动；而在二维Janus磁体MnSTe中，主要通过改变材料参数（如交换常数、各向异性常数等）来调控斯格明子，不同的材料参数会对斯格明子的稳定性和尺寸产生显著影响。综合两个案例，利用微磁学新型模型产生和调控斯格明子的普适性规律和经验如下：DMI是产生斯格明子的关键因素之一，在材料设计中，应寻找或创造能够增强DMI的条件，如利用材料的界面效应或原子结构的非对称性来诱导DMI。外加磁场是调控斯格明子的有效手段，通过精确控制外加磁场的强度和方向，可以实现对斯格明子多种特性的调控。材料参数对斯格明子的产生和特性有着重要影响，在研究和应用中，需要深入了解不同材料参数的作用，通过调整材料参数来优化斯格明子的性能。在实际应用中，还可以结合多种调控手段，如同时利用外加磁场和电流，或者同时改变材料参数和外加磁场等，以实现对斯格明子的更精确调控。六、结论与展望6.1研究成果总结本研究成功构建了一种新型微磁学模型，该模型在传统模型的基础上进行了创新和优化，为磁性斯格明子的研究提供了更精确、更全面的理论工具。通过引入基于第一性原理计算得到的更精确的自旋-轨道相互作用项，全面考虑了自旋-轨道相互作用的高阶项以及其在不同晶体方向上的各向异性，能够准确描述强自旋-轨道耦合体系中磁性斯格明子的行为，有效解释了实验中观察到的斯格明子的非对称运动轨迹、与传统理论预测不同的响应行为等奇特动力学特性。考虑了多物理场之间的耦合效应，引入电场-磁场耦合项、温度-磁场耦合项等，能够更真实地模拟磁、电、温度等多物理场共同作用下磁性斯格明子的行为，为多场调控斯格明子提供了理论支持。利用新型微磁学模型对磁性斯格明子的产生过程进行了深入模拟研究。设定了二维正方形薄膜的模拟体系，选用[Ta/CoFeB/MgO]多层膜体系作为模拟材料，详细设置了材料参数和外部条件，包括