基于铁磁斯格明子的二极管：从微磁学理论到应用探索

基于铁磁斯格明子的二极管：从微磁学理论到应用探索一、引言1.1研究背景与意义随着信息技术的飞速发展，对电子器件的性能要求日益提高，传统电子学逐渐面临着诸多挑战，如能耗过高、尺寸难以进一步缩小以及量子隧穿效应等问题，这些瓶颈限制了其在未来高性能、小型化器件中的应用。在此背景下，自旋电子学应运而生，它利用电子的自旋属性来存储和处理信息，为突破传统电子学的困境提供了新的方向。自旋电子学不仅有望降低器件能耗，还能提升信息处理速度和存储密度，成为了当前凝聚态物理和材料科学领域的研究热点之一。磁性斯格明子作为自旋电子学中的关键研究对象，展现出独特的物理性质和潜在应用价值，引起了广泛关注。斯格明子最初由英国核物理学家TonySkyrme于20世纪60年代在描述非线性场理论中的非拓扑平庸类粒子解时提出，随后其概念被引入到磁性材料领域。磁性斯格明子是一种具有纳米尺度的拓扑准粒子，通常存在于非中心对称的手性磁性材料以及界面镜面对称性破缺的磁性薄膜材料中。它的自旋结构呈涡旋状，内部磁矩以特定方式排列，形成了稳定的拓扑结构，这种拓扑结构赋予了斯格明子许多奇特的物理性质。磁性斯格明子具有拓扑保护特性，使其对环境扰动具有较强的鲁棒性，能够在一定程度上保持自身的稳定性，不易受到外界微小干扰的影响。其尺寸可小至几纳米到几十纳米，这为实现高密度信息存储提供了可能，能够满足未来器件小型化的需求。斯格明子还具有低驱动电流密度的优势，在电流驱动下可表现出类粒子动力学行为，仅需较小的电流就能实现其运动和操控，这意味着在实际应用中能够显著降低能耗，提高能源利用效率。在自旋电子学领域，基于磁性斯格明子的器件展现出了巨大的应用潜力，例如磁性赛道存储器、逻辑门和纳米振荡器等。在磁性赛道存储器中，斯格明子作为信息载体，通过电流驱动在纳米线中移动，实现信息的存储和传输，与传统存储技术相比，具有更高的存储密度和更快的读写速度；基于斯格明子的逻辑门则利用斯格明子的存在、不存在或不同的运动状态来表示逻辑信息，有望实现低功耗、高速的逻辑运算。二极管作为一种基础的电子元件，在电子电路中起着至关重要的作用，其单向导电特性是实现电流整流、信号调制和逻辑电路等功能的基础。将铁磁斯格明子与二极管相结合，形成基于铁磁斯格明子的二极管，为自旋电子学器件的发展开辟了新的方向。这种新型二极管不仅能够继承斯格明子的诸多优势，如低能耗、小尺寸和拓扑稳定性，还能利用其独特的自旋特性实现传统二极管所不具备的功能，为解决传统二极管在某些应用场景中的局限性提供了可能。从理论研究角度来看，基于铁磁斯格明子的二极管涉及到磁学、电学和拓扑物理等多个学科领域的交叉，对其进行深入研究有助于揭示自旋-电荷相互作用的微观机制，拓展人们对拓扑自旋电子学的认识，为进一步开发新型自旋电子学器件提供坚实的理论基础。在实际应用方面，这种新型二极管有望在低功耗集成电路、高速信号处理和量子计算等领域发挥重要作用，推动信息技术朝着高性能、低能耗的方向发展。例如，在低功耗集成电路中，基于铁磁斯格明子的二极管能够降低电路的整体能耗，提高芯片的运行效率和稳定性；在高速信号处理中，其独特的单向导电特性和快速响应能力可实现对高频信号的高效处理和调制；在量子计算领域，斯格明子的拓扑特性与量子比特的特性相结合，可能为量子计算的发展带来新的突破。1.2铁磁斯格明子概述铁磁斯格明子作为磁性斯格明子的一种特殊类型，继承了磁性斯格明子的基本特性，同时在铁磁材料体系中展现出独特的物理性质和应用潜力。它是一种在铁磁材料中形成的具有拓扑保护特性的纳米尺度涡旋状磁结构，其内部磁矩呈特定的涡旋排列方式，从中心到边缘磁矩方向连续旋转，形成了一种稳定且独特的自旋构型。从拓扑学角度来看，铁磁斯格明子具有非零的拓扑电荷，这是其区别于其他普通磁结构的关键特征之一。拓扑电荷的存在使得斯格明子在受到一定程度的外界扰动时，能够保持自身结构的完整性，不易被轻易破坏或改变。这种拓扑保护机制源于其自旋结构在实空间中的非平庸缠绕方式，赋予了铁磁斯格明子较高的稳定性和抗干扰能力，为其在自旋电子学器件中的应用提供了坚实的物理基础。铁磁斯格明子的形成机制较为复杂，涉及多种相互作用的协同效应。在非中心对称的铁磁材料中，反对称的Dzyaloshinskii-Moriya（DM）相互作用起着至关重要的作用。DM相互作用是由于材料中原子的自旋-轨道耦合以及晶体结构的空间反演对称性破缺所导致的，它使得相邻原子的自旋之间产生了一种不对称的交换作用。这种不对称交换作用倾向于使自旋形成特定的手性排列，从而促进了涡旋状斯格明子结构的形成。铁磁交换相互作用也在斯格明子的形成过程中发挥着重要作用，它倾向于使相邻自旋平行排列，与DM相互作用相互竞争和协调，共同决定了斯格明子的最终稳定结构。在具有界面镜面对称性破缺的铁磁薄膜异质结构中，界面处的DM相互作用同样对铁磁斯格明子的形成和稳定起到关键作用。界面DM相互作用源于磁性层与相邻非磁性层之间的原子耦合，这种耦合导致了界面处自旋-轨道耦合的增强，进而产生了有效的DM相互作用。在这种情况下，铁磁薄膜中的磁各向异性、交换相互作用以及界面DM相互作用之间的微妙平衡，决定了斯格明子是否能够稳定存在以及其具体的尺寸、形状和拓扑性质。铁磁斯格明子的稳定存在需要满足一定的条件，这些条件与材料的微观结构、物理性质以及外部环境因素密切相关。材料的本征参数，如交换刚度、磁各向异性常数和DM相互作用强度等，对斯格明子的稳定性起着决定性作用。较大的交换刚度有助于维持自旋之间的有序排列，增强斯格明子的稳定性；而合适的磁各向异性则可以提供额外的能量壁垒，防止斯格明子在热扰动或外界磁场变化时发生结构转变或湮灭。DM相互作用强度的大小不仅影响斯格明子的形成，还决定了其尺寸和稳定性，较强的DM相互作用通常有利于形成尺寸较小、稳定性较高的斯格明子。外部条件，如温度和外加磁场，也对铁磁斯格明子的稳定存在有着重要影响。温度的升高会增加热扰动的能量，当热扰动能量超过斯格明子的拓扑保护能量壁垒时，斯格明子可能会发生热激活的湮灭或结构转变，导致其稳定性降低。因此，在实际应用中，需要选择具有较高居里温度或能够在较高温度下保持斯格明子稳定性的材料体系。外加磁场可以改变铁磁材料中自旋的取向和能量分布，从而影响斯格明子的稳定性和动力学行为。在一定的磁场范围内，适当的外加磁场可以增强斯格明子的稳定性，使其更容易被操控和应用；然而，当磁场超过一定阈值时，可能会破坏斯格明子的拓扑结构，导致其湮灭或转变为其他磁结构。1.3基于铁磁斯格明子的二极管研究现状近年来，基于铁磁斯格明子的二极管研究作为自旋电子学领域的新兴方向，受到了广泛关注，取得了一系列具有重要意义的研究成果，为该领域的发展奠定了基础。在理论研究方面，科研人员主要聚焦于铁磁斯格明子二极管的物理机制探索和性能优化分析。通过建立微观物理模型，深入研究了铁磁斯格明子在电场、磁场以及电流等外部激励下的运动规律和输运特性，揭示了其单向导电行为的内在机制。研究表明，利用电压控制磁晶各向异性效应（VCMA）导致的高垂直磁晶各向异性（PMA）的边缘和梯形区域，能够对斯格明子的运动产生显著影响，实现斯格明子的单向运动，从而为构建基于铁磁斯格明子的二极管提供了理论依据。通过数值模拟和理论计算，分析了不同材料参数和结构设计对二极管性能的影响，为器件的优化设计提供了指导。例如，研究发现合适的DM相互作用强度、磁各向异性常数以及材料的自旋-轨道耦合强度等因素，能够有效提高二极管的整流效率和稳定性。在实验研究方面，众多科研团队致力于基于铁磁斯格明子的二极管的制备和性能验证。通过先进的材料生长技术和微纳加工工艺，成功制备出了多种基于铁磁斯格明子的二极管原型器件，并对其性能进行了实验测试和表征。实验结果证实了基于铁磁斯格明子的二极管能够实现单向导电功能，展现出了良好的应用潜力。部分研究还探索了在不同环境条件下二极管的性能变化，为其实际应用提供了实验数据支持。尽管基于铁磁斯格明子的二极管研究已经取得了一定进展，但仍存在一些亟待解决的问题和研究空白。在材料方面，目前能够稳定形成铁磁斯格明子且适合用于二极管制备的材料体系相对有限，材料的生长和制备工艺还不够成熟，难以满足大规模生产和应用的需求。部分材料的斯格明子稳定性和尺寸均匀性较差，限制了二极管性能的进一步提升。在器件设计和制备工艺方面，如何优化器件结构以提高斯格明子的产生效率、传输速度和稳定性，以及如何降低器件的制备成本和工艺复杂性，仍然是当前面临的挑战。不同制备工艺对器件性能的影响机制尚不完全清楚，缺乏系统的研究和总结。在理论研究方面，虽然已经取得了一些重要成果，但仍存在一些尚未解决的问题。铁磁斯格明子与外部电场、磁场以及电流之间的相互作用机制还需要进一步深入研究，以提高对其动力学行为的精确理解和调控能力。目前的理论模型大多基于简化假设，与实际情况存在一定差异，如何建立更加准确、全面的理论模型，以描述铁磁斯格明子在复杂环境下的行为，是未来理论研究的重要方向之一。在实际应用方面，基于铁磁斯格明子的二极管还面临着与现有集成电路工艺的兼容性问题。如何将这种新型二极管有效地集成到现有电路系统中，实现与其他电子元件的协同工作，是其走向实际应用的关键步骤。目前对于基于铁磁斯格明子的二极管在复杂电路环境中的可靠性和稳定性研究还相对较少，需要进一步加强这方面的研究，以确保其在实际应用中的性能和寿命。二、微磁学理论基础2.1微磁学基本概念微磁学是磁学领域中专注于介观尺度下铁磁体磁化过程研究的重要分支。这一尺度范围有着独特的界定，它足够大，使得原子的具体大小在研究中可以被忽略不计，从而能够将材料的磁学特性视为连续的整体来处理；然而，它又足够小，能够清晰地展现出磁畴的精细结构，为深入探究磁学现象提供了关键的微观视角。在微磁学的研究范畴内，主要聚焦于解决两类核心问题：静微磁学和动微磁学。静微磁学致力于通过能量最小化的原理和方法，来探寻系统的稳定状态解。这意味着在给定的条件下，通过调整磁矩的方向和分布，使得整个系统的总磁学能量达到最小值，此时系统处于最稳定的状态。这种方法对于理解材料在静态条件下的磁畴结构和磁化状态具有重要意义，能够揭示材料内部磁矩的平衡分布规律。动微磁学则着重于求解朗道-利夫希兹-吉尔伯特方程（Landau–Lifshitz–Gilbert,LLG），以此来获取系统的动力学解。该方程全面描述了磁矩随时间变化的进动特性，包括磁矩绕有效磁场的进动行为以及体系能量趋向最低能量时的阻尼过程。通过求解这一方程，能够深入了解在各种外部激励（如磁场、电场、温度变化等）作用下，磁矩如何随时间动态演化，从而为研究磁性材料在动态过程中的行为提供了坚实的理论基础。微磁学在现代磁性材料研究和自旋电子器件开发中占据着举足轻重的地位，发挥着不可替代的作用。在计算机硬盘的磁头和磁介质研发领域，微磁学的应用尤为关键。计算机硬盘作为重要的存储设备，其性能的提升依赖于对磁头和磁介质中磁畴结构和磁化过程的精确控制。通过微磁学理论和模拟方法，可以深入研究磁头在读写过程中与磁介质之间的相互作用，优化磁头的设计和磁介质的材料特性，从而提高硬盘的存储密度、读写速度和数据可靠性。例如，通过微磁学模拟，可以精确分析磁头产生的磁场分布对磁介质中磁畴翻转的影响，进而设计出更高效的磁头结构，实现更稳定、快速的数据读写操作。在永磁体的研发中，微磁学同样发挥着重要作用。永磁体广泛应用于电机、传感器、磁共振成像等众多领域，其性能的优劣直接影响到相关设备的性能和效率。利用微磁学方法，可以研究永磁体内部的磁畴结构和磁化机制，通过调整材料的成分、组织结构和制备工艺，优化永磁体的磁性能，如提高剩磁、矫顽力和最大磁能积等。例如，通过微磁学模拟，可以分析不同的稀土元素添加对永磁体磁晶各向异性的影响，从而指导研发人员选择合适的材料配方和制备工艺，制备出高性能的永磁体材料。2.2磁性材料中的能量项在磁性材料中，存在多种能量项，它们共同决定了材料的磁学性质和磁结构的稳定性。这些能量项包括交换能、磁晶各向异性能、塞曼能（外磁场能）、退磁能（形状各向异性能）以及DM磁相互作用能等。深入理解这些能量项的物理本质和相互作用机制，对于研究铁磁斯格明子的性质和行为具有至关重要的意义。2.2.1交换能交换能是磁性材料中至关重要的能量项之一，它源于电子之间的交换相互作用，本质上是一种量子力学效应。在铁磁材料中，相邻原子的电子云存在一定程度的重叠，根据量子力学的泡利不相容原理，两个电子不能同时处于相同的量子态，这就导致了电子之间的交换相互作用。这种交换相互作用使得相邻原子的自旋倾向于平行排列，从而降低系统的能量。从微观角度来看，交换能可以用海森堡交换模型来描述。在该模型中，相邻原子的自旋之间存在一种交换相互作用能，其表达式为E_{ex}=-2J\sum_{i,j}S_{i}\cdotS_{j}其中，J为交换积分，反映了交换相互作用的强度；S_{i}和S_{j}分别为第i和第j个原子的自旋矢量；求和符号\sum_{i,j}表示对所有相邻原子对进行求和。当J>0时，交换相互作用倾向于使相邻自旋平行排列，此时系统的交换能最低；当J<0时，相邻自旋倾向于反平行排列。在铁磁斯格明子的自旋结构中，交换能起着关键作用。由于交换能倾向于使相邻自旋平行排列，在斯格明子的核心区域，自旋的排列较为紧密，以满足交换能最小化的要求。随着距离斯格明子核心距离的增加，自旋的方向逐渐发生变化，形成涡旋状的结构。这种自旋结构的形成是交换能与其他能量项（如DM相互作用能、磁晶各向异性能等）相互竞争和平衡的结果。如果交换能在各种能量项中占主导地位，那么自旋将尽可能地保持平行排列，斯格明子的尺寸可能会受到限制，因为较大尺寸的斯格明子需要更多的自旋方向变化，这将增加交换能。反之，如果其他能量项（如DM相互作用能）较强，它们可以克服交换能的作用，使得自旋形成具有特定手性和拓扑结构的斯格明子。2.2.2磁晶各向异性能磁晶各向异性能是磁性材料中另一个重要的能量项，它起源于晶体结构中原子排列的各向异性以及电子的轨道运动与晶体场之间的相互作用。在磁性材料中，尽管自发磁化主要源于自旋间的交换作用，且这种交换作用本质上是各向同性的，但由于晶体结构的特殊性，使得自发磁化强度在晶体中倾向于沿着一个或几个特定的方向，这些方向被称为易轴。当磁化强度偏离易轴方向时，系统的能量会增加，这种与磁化方向有关的自由能增量即为磁晶各向异性能。磁晶各向异性能的产生主要有两个因素：一是电子的轨道运动与晶体场之间的相互作用，这是磁晶各向异性产生的主要来源；二是磁偶极子间的相互作用。电子在原子核周围运动时，其轨道会受到晶体场的影响，不同方向的晶体场对电子轨道的作用不同，从而导致电子云的分布具有各向异性。这种各向异性的电子云分布与自旋相互作用，使得自旋的取向也具有各向异性，进而产生磁晶各向异性能。磁偶极子间的相互作用也会对磁晶各向异性能产生一定贡献，但相对较小。对于立方晶体，磁晶各向异性能可以用磁晶各向异性常数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}+\cdots其中，\alpha_{1}、\alpha_{2}和\alpha_{3}分别为磁化强度方向与晶体坐标轴之间夹角的方向余弦；K_{1}和K_{2}等为磁晶各向异性常数，它们反映了磁晶各向异性的强弱和特性。当K_{1}>0时，[100]方向为易磁化方向；当K_{1}<0时，[111]方向为易磁化方向。在铁磁斯格明子体系中，磁晶各向异性能对斯格明子的取向和稳定性有着重要影响。斯格明子的自旋结构需要在满足其他能量项（如交换能、DM相互作用能等）的同时，尽量使磁晶各向异性能最低。当磁晶各向异性较强时，斯格明子的取向会受到易轴方向的限制，倾向于使自身的自旋结构与易轴方向保持一定的关系。如果易轴方向与外加磁场方向或其他外部激励方向不一致，可能会导致斯格明子在运动或响应外部激励时受到额外的能量阻碍，从而影响其动力学行为和稳定性。2.2.3塞曼能（外磁场能）塞曼能，也称为外磁场能，是磁性材料在外加磁场作用下所具有的能量。当磁性材料置于外磁场H中时，材料中的磁矩M会与外磁场相互作用，从而产生塞曼能。其表达式为E_{Z}=-\mu_{0}\int_{V}M\cdotHdV其中，\mu_{0}为真空磁导率，V为材料的体积。该公式表明，塞曼能的大小与磁矩和外磁场的夹角以及它们的乘积成正比。当磁矩与外磁场方向平行时，塞曼能最低，系统处于相对稳定的状态；当磁矩与外磁场方向垂直时，塞曼能最高。从微观角度来看，外磁场对磁矩的作用可以看作是磁场对磁矩产生一个力矩，使磁矩有沿着外磁场方向排列的趋势。根据经典电磁学理论，磁矩在外磁场中受到的力矩为\tau=M\timesH，这个力矩会促使磁矩发生转动，以降低系统的塞曼能。在量子力学中，电子的自旋磁矩与外磁场的相互作用导致原子能级的分裂，这也是塞曼效应的本质。在铁磁斯格明子的研究中，外磁场是操控斯格明子的重要手段之一。通过施加合适的外磁场，可以改变斯格明子的自旋结构、稳定性和运动状态。当外磁场强度较小时，外磁场主要影响斯格明子的取向，使其磁矩有沿着外磁场方向旋转的趋势，从而改变斯格明子的整体方向。随着外磁场强度的增加，外磁场对斯格明子的作用逐渐增强，可能会导致斯格明子的尺寸发生变化，甚至在一定条件下，外磁场可以使斯格明子发生湮灭或产生新的斯格明子。当外磁场超过一定阈值时，斯格明子的拓扑结构可能会被破坏，其内部的自旋排列发生剧烈变化，最终导致斯格明子消失。外磁场还可以与其他外部激励（如电流、温度梯度等）相结合，实现对斯格明子更加复杂和精确的操控。2.2.4退磁能（形状各向异性能）退磁能，又称为形状各向异性能，是由于磁性材料的形状不规则而产生的一种能量项。当磁性材料被磁化后，在材料内部会产生退磁场H_{d}，退磁场的方向与磁化强度M的方向相反，其大小与材料的形状和磁化强度有关。退磁能就是由退磁场与磁化强度相互作用所产生的能量，其表达式为E_{d}=-\frac{1}{2}\mu_{0}\int_{V}M\cdotH_{d}dV退磁场的产生源于磁性材料表面的磁极分布。当材料被磁化时，在材料表面会出现磁极，这些磁极会产生一个与磁化强度方向相反的磁场，即退磁场。对于不同形状的磁性材料，其表面磁极的分布情况不同，因此退磁场的大小和方向也不同。对于一个细长的磁性材料，当它沿着长度方向被磁化时，退磁场相对较小；而当它沿着垂直于长度方向被磁化时，退磁场相对较大。这是因为在细长材料沿着长度方向磁化时，表面磁极的分布相对较分散，产生的退磁场较弱；而在垂直于长度方向磁化时，表面磁极相对集中，退磁场较强。在铁磁斯格明子体系中，退磁能对斯格明子的形状和稳定性有着重要影响。退磁能倾向于使斯格明子的形状变得更加规则，以减小退磁场的影响。在某些情况下，退磁能可能会导致斯格明子的形状发生畸变，甚至影响其稳定性。如果斯格明子的形状不规则，会导致其表面的磁极分布不均匀，从而产生较强的退磁场，增加系统的退磁能。为了降低退磁能，斯格明子可能会调整自身的形状，使其更加接近球形或其他规则形状。然而，这种形状的调整可能会与其他能量项（如交换能、磁晶各向异性能等）产生冲突，从而影响斯格明子的稳定性。在实际应用中，需要综合考虑各种能量项的影响，通过优化材料的形状和结构，来提高斯格明子的稳定性和性能。2.2.5DM磁相互作用能DM磁相互作用能，即Dzyaloshinskii-Moriya相互作用能，是一种源于材料中原子的自旋-轨道耦合以及晶体结构的空间反演对称性破缺的能量项。在非中心对称的磁性材料中，相邻原子的自旋之间存在一种反对称的交换相互作用，这就是DM相互作用。其相互作用能的表达式为E_{DM}=\sum_{i,j}D_{ij}\cdot(S_{i}\timesS_{j})其中，D_{ij}为DM矢量，它的大小和方向与材料的晶体结构和原子间的距离有关；S_{i}和S_{j}分别为第i和第j个原子的自旋矢量。DM矢量的方向垂直于S_{i}和S_{j}所构成的平面，其大小反映了DM相互作用的强度。DM相互作用的产生机制较为复杂，它与材料中原子的电子结构密切相关。在具有强自旋-轨道耦合的材料中，电子的自旋与其轨道运动之间存在强烈的相互作用。当晶体结构缺乏空间反演对称性时，这种自旋-轨道耦合会导致相邻原子的自旋之间产生一种不对称的交换作用，即DM相互作用。这种相互作用倾向于使自旋形成特定的手性排列，从而促进了具有涡旋状结构的铁磁斯格明子的形成。在铁磁斯格明子的形成和稳定过程中，DM磁相互作用能起着至关重要的作用。它与交换能相互竞争和协调，共同决定了斯格明子的自旋结构和稳定性。由于交换能倾向于使相邻自旋平行排列，而DM相互作用倾向于使自旋形成特定的手性排列，这两种相互作用的平衡决定了斯格明子的最终结构。在合适的DM相互作用强度下，能够形成稳定的斯格明子结构，并且DM相互作用强度的变化会影响斯格明子的尺寸和稳定性。较强的DM相互作用通常有利于形成尺寸较小、稳定性较高的斯格明子。DM相互作用还可以影响斯格明子的动力学行为，例如在电流驱动下，DM相互作用会导致斯格明子产生横向偏移，即斯格明子霍尔效应。2.3静态微磁学方程在静态微磁学中，核心目标是确定磁性系统的稳定状态，这一过程通过求解系统的总磁学能量达到极小值来实现。系统的总能量E_{total}由多个能量项组成，包括交换能E_{ex}、磁晶各向异性能E_{k}、塞曼能E_{Z}、退磁能E_{d}以及DM磁相互作用能E_{DM}等，即E_{total}=E_{ex}+E_{k}+E_{Z}+E_{d}+E_{DM}+\cdots交换能E_{ex}描述了相邻原子自旋间的相互作用，其表达式为E_{ex}=-A\int_{V}(\nabla\mathbf{m})^{2}dV其中，A为交换刚度常数，反映了交换相互作用的强度；\mathbf{m}为归一化的磁化强度矢量，\mathbf{m}=\frac{\mathbf{M}}{M_{s}}，M_{s}为饱和磁化强度；\nabla\mathbf{m}表示磁化强度矢量的空间梯度，它描述了自旋在空间中的变化情况。交换能的作用是倾向于使相邻自旋平行排列，以降低系统的能量。当自旋排列不均匀，即\nabla\mathbf{m}较大时，交换能会增加。在铁磁斯格明子的自旋结构中，交换能在斯格明子核心区域起着重要作用，促使自旋紧密排列，以满足交换能最小化的要求。磁晶各向异性能E_{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}+\cdots其中，K_{1}和K_{2}等为磁晶各向异性常数，它们反映了磁晶各向异性的强弱和特性；\alpha_{1}、\alpha_{2}和\alpha_{3}分别为磁化强度方向与晶体坐标轴之间夹角的方向余弦。当磁化强度偏离易磁化方向时，磁晶各向异性能会增加，因此磁晶各向异性能对铁磁斯格明子的取向和稳定性有着重要影响。塞曼能E_{Z}是由于外磁场H与磁矩M相互作用产生的能量，其表达式为E_{Z}=-\mu_{0}\int_{V}\mathbf{M}\cdot\mathbf{H}dV其中，\mu_{0}为真空磁导率，V为材料的体积。外磁场倾向于使磁矩沿着磁场方向排列，以降低塞曼能。在铁磁斯格明子的研究中，外磁场是操控斯格明子的重要手段之一，通过施加合适的外磁场，可以改变斯格明子的自旋结构、稳定性和运动状态。退磁能E_{d}是由于磁性材料的形状不规则导致内部产生退磁场H_{d}而引起的能量，其表达式为E_{d}=-\frac{1}{2}\mu_{0}\int_{V}\mathbf{M}\cdot\mathbf{H}_{d}dV退磁场的方向与磁化强度方向相反，它会使系统的能量增加。退磁能对铁磁斯格明子的形状和稳定性有着重要影响，倾向于使斯格明子的形状变得更加规则，以减小退磁场的影响。DM磁相互作用能E_{DM}源于材料中原子的自旋-轨道耦合以及晶体结构的空间反演对称性破缺，其表达式为E_{DM}=\int_{V}\mathbf{D}\cdot(\mathbf{m}\times\nabla\mathbf{m})dV其中，\mathbf{D}为DM矢量，它的大小和方向与材料的晶体结构和原子间的距离有关。DM相互作用倾向于使自旋形成特定的手性排列，从而促进了具有涡旋状结构的铁磁斯格明子的形成。为了求解系统总能量的极小值，通常采用数值计算方法。常见的数值计算方法包括共轭梯度法、有限差分法和有限元法等。共轭梯度法是一种迭代优化算法，它通过不断迭代寻找函数的下降方向，逐步逼近函数的极小值。在微磁学中，共轭梯度法常用于求解能量泛函的极小值问题，它具有收敛速度快、计算效率高等优点。有限差分法是将连续的空间区域离散化为一系列网格点，通过在这些网格点上对微磁学方程进行离散化处理，将微分方程转化为代数方程组进行求解。有限差分法的优点是计算简单、易于实现，但其精度受到网格尺寸的限制，在处理复杂几何形状和边界条件时可能存在一定的困难。有限元法是将求解区域划分为有限个单元，通过在每个单元上构造插值函数，将微磁学方程转化为一组线性代数方程组进行求解。有限元法能够灵活处理复杂的几何形状和边界条件，具有较高的计算精度，但其计算量较大，对计算机资源的要求较高。2.4动态微磁学方程2.4.1朗道-栗弗席兹-吉尔伯特（LLG）方程在动态微磁学领域，朗道-栗弗席兹-吉尔伯特（LLG）方程占据着核心地位，它是描述磁性材料中磁化强度随时间演化的重要方程。该方程全面且准确地刻画了磁矩在各种相互作用下的动态行为，对于深入理解铁磁斯格明子的动力学特性以及相关自旋电子器件的工作原理具有不可替代的作用。LLG方程的一般形式为\frac{\partial\mathbf{m}}{\partialt}=-\gamma\mathbf{m}\times\mathbf{H}_{eff}+\alpha\mathbf{m}\times\frac{\partial\mathbf{m}}{\partialt}其中，\mathbf{m}为归一化的磁化强度矢量，\mathbf{m}=\frac{\mathbf{M}}{M_{s}}，M_{s}为饱和磁化强度；t表示时间；\gamma是旋磁比，它反映了磁矩在磁场作用下的进动特性，是一个与材料性质相关的常数；\mathbf{H}_{eff}为有效磁场，它综合了材料内部各种能量项所对应的磁场贡献，包括交换能、磁晶各向异性能、塞曼能、退磁能以及DM磁相互作用能等产生的有效磁场；\alpha是阻尼系数，用于描述磁性系统在演化过程中能量的耗散情况，其取值范围通常在0到1之间，不同的材料具有不同的阻尼系数值，它对磁矩的动态响应速度和稳定性有着重要影响。方程右侧的第一项-\gamma\mathbf{m}\times\mathbf{H}_{eff}描述了磁矩绕有效磁场的进动行为。根据经典电磁学理论，磁矩在磁场中会受到一个力矩的作用，该力矩会使磁矩绕磁场方向做进动运动。这一项体现了磁矩与有效磁场之间的相互作用，决定了磁矩的进动频率和方向。在铁磁斯格明子体系中，磁矩的进动行为对斯格明子的稳定性和动力学特性有着重要影响。如果有效磁场发生变化，磁矩的进动也会相应改变，从而可能导致斯格明子的自旋结构发生变化。第二项\alpha\mathbf{m}\times\frac{\partial\mathbf{m}}{\partialt}为阻尼项，它描述了体系能量趋向最低能量时的过程。当磁性系统处于非平衡状态时，磁矩的运动将导致系统能量的变化。阻尼项的存在使得磁矩在运动过程中逐渐消耗能量，最终趋向于使系统能量达到最低的稳定状态。阻尼系数\alpha的大小决定了能量耗散的速度，\alpha越大，能量耗散越快，磁矩达到稳定状态的速度也越快。在实际应用中，阻尼系数的选择需要综合考虑器件的性能要求和工作条件，例如在一些需要快速响应的自旋电子器件中，可能需要选择较小的阻尼系数，以提高器件的响应速度；而在一些需要稳定磁状态的应用中，则可能需要较大的阻尼系数，以确保磁矩的稳定性。在铁磁斯格明子动力学研究中，LLG方程发挥着关键作用。通过求解LLG方程，可以深入了解斯格明子在各种外部激励（如磁场、电场、电流等）作用下的运动规律和自旋结构的变化。在研究电流驱动下铁磁斯格明子的运动时，通过将电流相关的项引入有效磁场\mathbf{H}_{eff}中，然后求解LLG方程，可以得到斯格明子的运动速度、方向以及其自旋结构的动态演化情况。这对于理解基于铁磁斯格明子的器件（如磁性赛道存储器、逻辑门等）的工作原理和性能优化具有重要意义。通过数值模拟方法求解LLG方程，还可以预测不同材料参数和外部条件下铁磁斯格明子的动力学行为，为实验研究提供理论指导。例如，通过改变材料的交换刚度、磁各向异性常数、DM相互作用强度以及阻尼系数等参数，模拟斯格明子的动力学过程，可以分析这些参数对斯格明子稳定性和运动特性的影响，从而为材料的设计和选择提供依据。2.4.2考虑新物理效应的LLG方程拓展随着自旋电子学的深入发展以及对铁磁斯格明子研究的不断推进，一些新的物理效应逐渐被揭示和认识，这些效应在描述铁磁斯格明子的动力学行为时具有重要影响，因此需要对传统的LLG方程进行拓展，以准确描述这些新现象。自旋转移矩（Spin-TransferTorque，STT）效应是其中一个重要的新物理效应。当自旋极化电流通过磁性材料时，由于电子的自旋与磁性材料中的磁矩相互作用，会产生一个力矩作用在磁矩上，这就是自旋转移矩效应。自旋转移矩效应可以用一个新的力矩项添加到LLG方程中，从而对传统的LLG方程进行拓展。拓展后的LLG方程形式为\frac{\partial\mathbf{m}}{\partialt}=-\gamma\mathbf{m}\times\mathbf{H}_{eff}+\alpha\mathbf{m}\times\frac{\partial\mathbf{m}}{\partialt}+\frac{\gamma\hbar}{eM_{s}}\mathbf{j}_{s}\cdot(\mathbf{m}\times\mathbf{p})其中，\frac{\gamma\hbar}{eM_{s}}\mathbf{j}_{s}\cdot(\mathbf{m}\times\mathbf{p})为自旋转移矩项；\hbar是约化普朗克常数；e为电子电荷量；\mathbf{j}_{s}是自旋极化电流密度；\mathbf{p}是自旋极化方向。自旋转移矩效应的引入使得磁矩的动力学行为更加复杂，它可以导致磁矩的翻转、进动频率的改变以及斯格明子的运动和操控。在基于铁磁斯格明子的磁性赛道存储器中，利用自旋转移矩效应可以实现对斯格明子的高效写入和移动操作，通过控制自旋极化电流的大小和方向，可以精确地控制斯格明子在纳米线中的位置和运动状态，从而实现信息的存储和读取。自旋-轨道矩（Spin-OrbitTorque，SOT）效应也是一个不可忽视的新物理效应。在具有强自旋-轨道耦合的材料中，当有电流通过时，会产生自旋-轨道矩，它同样会对磁矩产生力矩作用。自旋-轨道矩效应主要包括内禀自旋-轨道矩和外禀自旋-轨道矩。内禀自旋-轨道矩源于材料内部的晶体结构和电子能带结构，与材料的本征性质相关；外禀自旋-轨道矩则通常是由于材料中的杂质、缺陷或界面等因素引起的。将自旋-轨道矩效应考虑在内，LLG方程进一步拓展为\frac{\partial\mathbf{m}}{\partialt}=-\gamma\mathbf{m}\times\mathbf{H}_{eff}+\alpha\mathbf{m}\times\frac{\partial\mathbf{m}}{\partialt}+\frac{\gamma\hbar}{eM_{s}}\mathbf{j}_{s}\cdot(\mathbf{m}\times\mathbf{p})+\frac{\gamma\hbar}{eM_{s}}\mathbf{j}\times\mathbf{\xi}\cdot\mathbf{m}其中，\frac{\gamma\hbar}{eM_{s}}\mathbf{j}\times\mathbf{\xi}\cdot\mathbf{m}为自旋-轨道矩项；\mathbf{j}是电流密度；\mathbf{\xi}是与自旋-轨道耦合相关的张量，其大小和方向取决于材料的性质和晶体结构。自旋-轨道矩效应在铁磁斯格明子的动力学中有着重要应用，它可以提供一种新的方式来操控斯格明子。与自旋转移矩效应相比，自旋-轨道矩效应具有一些独特的优势，例如它可以在较低的电流密度下实现对磁矩的有效操控，并且对材料的要求相对较低，这使得它在实际应用中具有更大的潜力。在一些新型的自旋电子器件中，利用自旋-轨道矩效应可以实现对铁磁斯格明子的快速、低能耗操控，为器件的小型化和高性能化提供了可能。三、基于铁磁斯格明子的二极管结构设计与原理3.1二极管的工作原理在电子学领域，二极管作为一种基础性的电子元件，以其独特的单向导电特性而被广泛应用。二极管的基本结构主要基于PN结，这是由P型半导体和N型半导体紧密结合而成的结构。在P型半导体中，主要的载流子为空穴，呈现出正电荷的特性；而在N型半导体中，电子则是主要载流子，携带负电荷。当P型半导体和N型半导体相互接触形成PN结时，在交界面附近会发生一系列复杂的物理过程。由于P型半导体和N型半导体中载流子浓度存在显著差异，N型半导体中的电子会自发地向P型半导体扩散，而P型半导体中的空穴则会向N型半导体扩散。这种扩散过程导致在PN结附近形成一个特殊的区域，被称为空间电荷区，也叫耗尽层。在这个区域内，由于电子和空穴的扩散与复合，使得该区域内几乎没有自由载流子，呈现出高电阻的特性。同时，由于电子和空穴的扩散，在PN结的两侧分别积累了正、负电荷，从而形成了一个内建电场，其方向从N区指向P区。这个内建电场的存在会阻碍载流子的进一步扩散，当扩散作用与内建电场的阻碍作用达到动态平衡时，空间电荷区的宽度和内建电场的强度就会保持相对稳定。当在PN结两端施加正向电压时，即P区接电源正极，N区接电源负极，外加电场的方向与内建电场的方向相反。在这种情况下，外加电场会削弱内建电场的作用，使得空间电荷区变窄。随着外加电压的逐渐增大，空间电荷区的电阻不断减小，多数载流子（P区的空穴和N区的电子）能够克服内建电场的阻碍，顺利地通过PN结，从而形成较大的正向电流。此时，PN结处于导通状态，二极管表现出低电阻特性，电流可以较为顺畅地通过。相反，当在PN结两端施加反向电压时，即P区接电源负极，N区接电源正极，外加电场的方向与内建电场的方向相同。这会进一步增强内建电场的强度，使得空间电荷区变宽。在这种情况下，多数载流子受到内建电场和外加电场的双重阻碍，几乎无法通过PN结。虽然少数载流子（P区的电子和N区的空穴）会在电场作用下形成反向电流，但由于少数载流子的浓度很低，所以反向电流非常小，通常可以忽略不计。此时，PN结处于截止状态，二极管表现出高电阻特性，电流几乎无法通过。这种单向导电特性使得二极管在电子电路中具有至关重要的作用。在整流电路中，二极管可以将交流电转换为直流电，通过利用其正向导通、反向截止的特性，将交流电的负半周截止，只允许正半周通过，从而实现整流功能。在信号调制电路中，二极管可以用于对信号进行调制和解调，通过控制二极管的导通和截止状态，实现对信号的频率、幅度等参数的调整。在逻辑电路中，二极管也可以作为基本的逻辑单元，通过不同的连接方式和输入信号，实现与、或、非等逻辑运算功能。3.2基于铁磁斯格明子的二极管结构设计基于铁磁斯格明子的二极管结构设计是实现其独特功能的关键，目前科研人员已提出多种具有创新性的设计方案，每种方案都有其独特的结构特点和工作原理。其中一种设计方案是利用电压控制磁晶各向异性效应（VCMA）来构建基于铁磁斯格明子的二极管。在这种设计中，通过施加外部电压，利用VCMA效应在铁磁薄膜中形成高垂直磁晶各向异性（PMA）的边缘和梯形区域。高PMA边缘区域能够有效地避免斯格明子信号在运动过程中的湮灭，同时还能提高其传输速度。这是因为高PMA边缘提供了额外的能量壁垒，使得斯格明子在运动到该区域时更加稳定，不易受到外界干扰而湮灭。高PMA边缘的存在改变了斯格明子周围的能量分布，使得斯格明子在该区域的运动阻力减小，从而提高了传输速度。高PMA梯形区域则在实现斯格明子的单向运动方面发挥着关键作用。其工作原理基于斯格明子在不同能量势垒下的运动特性。当斯格明子从低磁晶各向异性区域向高PMA梯形区域运动时，由于梯形区域的特殊形状和能量分布，斯格明子受到一个与运动方向相关的力的作用。在正向运动时，这个力有助于斯格明子克服能量势垒，顺利通过梯形区域；而在反向运动时，该力则会阻碍斯格明子的运动，使其难以通过，从而实现了斯格明子的单向运动。基于这种特性，构建出了类似于PN结的基于铁磁斯格明子的二极管原型。这种二极管不仅能够实现赛道存储器与二极管的功能转换，还为基于斯格明子的自旋电子器件的发展提供了新的思路。在实际应用中，通过控制外部电压的大小和方向，可以灵活地调节二极管的工作状态，实现对斯格明子运动的精确控制。另一种具有代表性的设计方案是基于增强的赛道边缘和缺口设计。在这种设计中，通过对铁磁纳米赛道的边缘进行特殊处理，使其具有增强的特性，并在赛道上开设特定尺寸和位置的缺口。赛道边缘的矩形缺口对磁斯格明子具有钉扎效应，丰富了磁斯格明子的动力学行为。当斯格明子运动到缺口处时，由于缺口处的能量分布发生变化，斯格明子会受到一个额外的钉扎力，使其运动状态发生改变。这种钉扎效应可以用于控制斯格明子的位置和运动方向，为实现二极管的功能提供了基础。通过巧妙地设计赛道的形状和尺寸，以及缺口的位置和大小，可以在同一赛道上实现逻辑与门、或门、非门以及二极管功能。这种设计便于器件的集成化和小型化，符合未来电子器件发展的趋势。在实现二极管功能时，利用斯格明子在赛道上的单向传输特性，通过控制电流的方向和大小，使得斯格明子只能沿着一个方向运动，从而实现了类似于传统二极管的单向导电功能。该设计使用的驱动电流密度比以往的大多数相关设计要小得多，有效地降低了器件的能耗，提高了能源利用效率。通过模拟计算铁磁纳米赛道中磁斯格明子的动力学特性，形成了一个具有较低磁各向异性的通道，进一步优化了斯格明子的运动性能，并证明了该设计在有限温度下的可行性。3.3工作原理分析基于铁磁斯格明子的二极管利用铁磁斯格明子实现单向导通的原理，是基于微磁学理论中多种物理机制的协同作用，其核心在于对斯格明子运动的精确控制和利用其独特的自旋结构与输运特性。从微磁学理论角度来看，铁磁斯格明子的运动受到多种能量项的共同影响，这些能量项的相互作用决定了斯格明子在材料中的行为。交换能倾向于使相邻自旋平行排列，它在斯格明子的形成和稳定中起着基础性作用，维持着斯格明子内部自旋结构的有序性。磁晶各向异性能则决定了磁化强度在材料中的易磁化方向，当斯格明子的自旋方向与易磁化方向不一致时，会增加系统的能量，从而影响斯格明子的稳定性和运动状态。DM磁相互作用能源于材料中原子的自旋-轨道耦合以及晶体结构的空间反演对称性破缺，它倾向于使自旋形成特定的手性排列，对斯格明子的形成和运动方向具有关键影响。在基于铁磁斯格明子的二极管中，这些能量项在不同区域的分布和相互作用，构建起了实现单向导通的物理基础。以利用电压控制磁晶各向异性效应（VCMA）构建的二极管为例，通过施加外部电压，利用VCMA效应在铁磁薄膜中形成高垂直磁晶各向异性（PMA）的边缘和梯形区域。在高PMA边缘区域，由于磁晶各向异性能显著增加，形成了一个能量壁垒，这使得斯格明子在该区域具有更高的稳定性。当斯格明子运动到该区域时，需要克服更大的能量障碍才能改变其运动状态，从而有效地避免了斯格明子信号在运动过程中的湮灭。高PMA边缘区域的存在还改变了斯格明子周围的能量分布，使得斯格明子在该区域的运动阻力减小，进而提高了其传输速度。高PMA梯形区域则是实现单向导通的关键部分。该区域的特殊形状和能量分布，使得斯格明子在其中受到一个与运动方向相关的力的作用。当斯格明子从低磁晶各向异性区域向高PMA梯形区域正向运动时，由于梯形区域的能量梯度和各向异性特性，斯格明子受到一个有助于其克服能量势垒的力，使其能够顺利通过梯形区域。具体来说，在正向运动过程中，斯格明子的自旋结构与梯形区域的能量分布相互匹配，使得斯格明子能够沿着特定的路径运动，且运动过程中能量的变化有利于其持续前进。当斯格明子反向运动时，情况则截然不同。反向运动时，斯格明子的自旋结构与梯形区域的能量分布不匹配，导致其受到一个阻碍运动的力。这个阻碍力使得斯格明子难以克服梯形区域的能量势垒，从而无法顺利通过，实现了斯格明子的单向运动。从能量角度分析，反向运动时斯格明子需要跨越更高的能量壁垒，而其自身的能量不足以支持其完成这一过程，因此被阻挡在梯形区域之外。基于这种对斯格明子单向运动的精确控制，类似于传统PN结二极管的功能得以实现。在电路中，当施加正向电压时，斯格明子能够沿着特定的路径单向传输，形成电流通路，实现正向导通；当施加反向电压时，斯格明子的运动被阻止，无法形成有效的电流通路，实现反向截止。这种基于铁磁斯格明子的单向导通机制，为自旋电子学器件的发展提供了新的思路和方法。与传统二极管相比，基于铁磁斯格明子的二极管不仅具有独特的自旋相关特性，还能够在更低的电流密度下工作，具有更低的能耗和更高的集成度潜力。在未来的集成电路应用中，有望实现更高性能、更低功耗的电路设计，推动信息技术向更高水平发展。四、微磁学模拟研究4.1模拟方法与软件选择微磁学模拟是研究基于铁磁斯格明子的二极管的重要手段，它能够深入揭示二极管内部的物理机制和性能特性。在进行微磁学模拟时，主要有有限差分法（FDM）、有限元法（FEM）和边界元法（BEM）等方法，每种方法都有其独特的原理、特点和适用场景。有限差分法是一种将连续的空间区域离散化为一系列网格点的数值计算方法。在微磁学模拟中，通过在这些网格点上对微磁学方程进行离散化处理，将微分方程转化为代数方程组进行求解。具体来说，对于描述磁矩演化的朗道-利夫希兹-吉尔伯特（LLG）方程，有限差分法将空间和时间进行离散，用差分近似代替方程中的导数项。例如，对于空间导数\frac{\partialm_x}{\partialx}，可以用\frac{m_x(x+\Deltax,t)-m_x(x,t)}{\Deltax}来近似，其中\Deltax是空间步长，m_x(x,t)表示在位置x和时间t时磁矩在x方向的分量。通过这种方式，将LLG方程转化为一组关于网格点上磁矩的代数方程，然后通过迭代求解这些方程，得到磁矩随时间和空间的演化。有限差分法的优点是计算简单、易于实现，其原理直观，编程实现相对容易，对于一些简单的几何形状和边界条件的问题，能够快速得到数值解。它的计算效率较高，在处理大规模问题时，能够在较短的时间内完成计算。但有限差分法的精度受到网格尺寸的限制，网格尺寸越小，计算精度越高，但同时计算量也会急剧增加。在处理复杂几何形状和边界条件时，有限差分法可能存在一定的困难，因为需要对网格进行特殊的处理和划分，以适应复杂的形状和边界。有限元法是另一种常用的微磁学模拟方法，它将求解区域划分为有限个单元。在每个单元上，通过构造插值函数，将微磁学方程转化为一组线性代数方程组进行求解。有限元法的基本思想是将连续的求解区域离散化为有限个小的单元，每个单元内的物理量用插值函数来近似表示。对于微磁学问题，通常选择合适的形状函数来描述单元内磁矩的分布。然后，根据变分原理或加权余量法，将微磁学方程转化为单元上的弱形式，进而得到整个求解区域的代数方程组。有限元法的优势在于能够灵活处理复杂的几何形状和边界条件。它可以根据求解区域的几何形状，将其划分为各种形状的单元，如三角形、四边形、四面体等，从而更好地适应复杂的结构。有限元法具有较高的计算精度，通过合理选择插值函数和单元尺寸，可以得到较为精确的数值解。但有限元法的计算量较大，对计算机资源的要求较高。由于需要对每个单元进行计算和存储，当求解区域较大或单元数量较多时，计算时间和内存需求会显著增加。有限元法的前处理过程相对复杂，需要进行网格划分、材料参数定义等操作，这些操作需要一定的专业知识和经验。边界元法是基于边界积分方程的一种数值计算方法。它将求解区域的边界离散化，通过在边界上求解积分方程，得到边界上的物理量，进而通过积分变换得到整个求解区域内的物理量。在微磁学模拟中，边界元法利用磁学的基本原理，将微磁学问题转化为边界积分方程。例如，对于静磁场问题，可以利用磁标势或磁矢势的边界积分方程来描述。通过将边界划分为一系列边界单元，在每个单元上对积分方程进行离散化处理，得到关于边界节点上物理量的代数方程组。边界元法的主要优点是只需要对边界进行离散化，大大减少了求解区域的维数，从而降低了计算量。在处理一些具有复杂边界的问题时，边界元法能够有效地减少计算量和内存需求。边界元法对于处理无限域或半无限域问题具有独特的优势，因为它可以通过边界条件来描述无限域的影响。但边界元法也存在一些局限性，它依赖于基本解的选取，对于一些复杂的材料和物理问题，可能难以找到合适的基本解。边界元法得到的系数矩阵通常是满秩的，这增加了求解的难度和计算量。边界元法在处理多连通区域或内部存在不同材料的问题时，可能会遇到一些困难。在众多微磁学模拟软件中，ObjectOrientedMicroMagneticFramework（OOMMF）、MuMax3和Magpar等软件被广泛应用，它们各自具有不同的特点和适用范围。OOMMF是一款由美国国家标准与技术研究院（NIST）开发的面向对象的微磁学框架软件。它具有强大的微磁学模拟能力，能够模拟各种微磁学问题，包括磁畴结构的形成、磁畴壁的动态行为、磁性材料的磁化过程以及退磁过程等。OOMMF通过数值求解Landau-Lifshitz-Gilbert方程来描述磁畴随时间的动态演化过程。在模拟基于铁磁斯格明子的二极管时，OOMMF可以精确地模拟斯格明子的形成、运动以及与外部电场、磁场的相互作用。通过设置合适的材料参数和边界条件，能够研究不同条件下二极管的性能和特性。OOMMF支持多种数值计算方法，能够有效利用高性能计算机资源，进行快速准确的计算。它是开源自由软件，提供了丰富的代码库和文档，方便研究者根据自己的需要进行二次开发和定制。OOMMF还提供直观的图形用户界面和数据可视化工具，能够方便用户观察和分析模拟结果。用户可以通过图形界面直观地查看磁矩分布、能量变化等模拟结果，便于理解和分析模拟过程。OOMMF适用于多种操作系统，包括Windows、Linux和MacOS等，具有良好的跨平台兼容性，便于在不同计算环境下部署和使用。MuMax3是一款专门用于微磁学模拟的开源软件，它在处理大规模微磁学问题时表现出色。MuMax3采用了高效的数值算法和并行计算技术，能够快速处理大规模的微磁学模拟任务。在模拟基于铁磁斯格明子的二极管时，如果需要考虑较大尺寸的器件或复杂的材料体系，MuMax3能够利用其并行计算能力，在较短的时间内完成模拟。它提供了丰富的物理模型和材料参数库，方便用户进行不同材料和物理场景的模拟。用户可以根据研究需求，选择合适的物理模型和材料参数，快速搭建模拟环境。MuMax3还具有良好的可视化功能，能够实时显示模拟过程中的磁矩分布和动态演化。用户可以通过可视化界面，直观地观察斯格明子在二极管中的运动轨迹和变化情况，为研究提供直观的依据。但MuMax3的学习曲线相对较陡，对于初学者来说，可能需要花费一定的时间来掌握其使用方法和各种功能。Magpar是一款为Windows操作系统设计的微磁学模拟软件，其名称由“magnetism”（磁性）和“parallel”（并行）的缩写组合而来，暗示了该软件在设计时考虑了模拟过程的并行计算能力。Magpar使用方便，具有直观的用户界面，使得用户能够轻松地进行操作。对于不具备深厚编程背景的用户来说，Magpar的易用性使其能够快速上手，进行微磁学模拟。该软件的功能包括模拟磁性材料的磁化翻转、磁致伸缩效应等。在研究基于铁磁斯格明子的二极管时，Magpar可以模拟二极管在不同条件下的磁化状态变化，以及斯格明子与其他磁结构的相互作用。但Magpar的功能相对较为单一，可能无法满足一些复杂的研究需求。在本研究中，综合考虑研究内容和需求，选择了OOMMF软件进行基于铁磁斯格明子的二极管的微磁学模拟。本研究重点关注二极管中斯格明子的形成、运动以及与外部电场、磁场的相互作用，OOMMF强大的模拟能力和丰富的功能能够满足这些研究需求。其开源特性也便于根据具体研究内容进行二次开发和定制，以实现更精确的模拟。4.2模拟参数设置在基于铁磁斯格明子的二极管微磁学模拟中，准确合理地设置模拟参数至关重要，这些参数涵盖材料参数、几何参数和外场参数等多个方面，它们的取值直接影响着模拟结果的准确性和可靠性，对深入理解二极管的性能和物理机制起着关键作用。在材料参数方面，选择具有合适特性的材料是模拟的基础。以常用的重金属/铁磁薄膜/氧化物结构体系为例，其中铁磁薄膜可选用CoFeB等材料。对于CoFeB薄膜，其饱和磁化强度M_s通常在一定范围内取值，一般为1.2\times10^6\A/m左右，这一数值反映了材料在饱和状态下的磁化能力，对铁磁斯格明子的形成和稳定性有着重要影响。交换刚度常数A也是一个关键参数，对于CoFeB薄膜，其值约为1.5\times10^{-11}\J/m，它决定了相邻原子自旋间的交换相互作用强度，影响着自旋的排列方式和斯格明子内部自旋结构的稳定性。磁晶各向异性常数K同样不可忽视，对于具有垂直磁晶各向异性的CoFeB薄膜，K的取值一般在1.0\times10^5\J/m^3左右，该参数决定了磁化强度在材料中的易磁化方向，对铁磁斯格明子的取向和运动状态有着重要作用。Dzyaloshinskii-Moriya（DM）相互作用强度D也是一个重要的材料参数，在这种结构体系中，D的值通常在3.0\mJ/m^2左右，它源于材料中原子的自旋-轨道耦合以及晶体结构的空间反演对称性破缺，对斯格明子的形成和运动方向具有关键影响。在几何参数设置方面，考虑基于电压控制磁晶各向异性效应（VCMA）构建的二极管结构。铁磁薄膜的厚度t一般设置为1.0\nm左右，薄膜厚度对铁磁斯格明子的稳定性和运动特性有显著影响，合适的厚度可以保证斯格明子在薄膜中稳定存在，并有利于其在外部激励下的运动。构建的高垂直磁晶各向异性（PMA）边缘区域的宽度w_1可设置为50\nm左右，该区域能够有效地避免斯格明子信号在运动过程中的湮灭，同时提高其传输速度，合适的宽度可以提供足够的能量壁垒，保证斯格明子的稳定性。高PMA梯形区域的尺寸也需要精确设置，例如其底边长度l_1可设置为100\nm左右，斜边长度l_2可设置为80\nm左右，梯形的高度h可设置为30\nm左右。这些尺寸的设置决定了梯形区域的能量分布和梯度，进而影响斯格明子在其中的运动行为，实现斯格明子的单向运动。在外场参数方面，首先考虑外加磁场。在研究铁磁斯格明子的运动和二极管的性能时，常施加一个平面内的磁场H。磁场强度H的大小可根据具体研究需求进行设置，一般在0\至\100\mT范围内变化。当磁场强度较小时，主要影响斯格明子的取向，使其磁矩有沿着外磁场方向旋转的趋势；随着磁场强度的增加，可能会导致斯格明子的尺寸发生变化，甚至在一定条件下，外磁场可以使斯格明子发生湮灭或产生新的斯格明子。对于电流参数，在模拟电流驱动铁磁斯格明子运动时，电流密度j的大小是一个关键参数。一般情况下，电流密度j可设置在1.0\times10^{11}\A/m^2左右，这一数值在实际实验中具有可实现性，且能够有效地驱动斯格明子运动。电流的方向也需要明确设置，通常设置为沿着铁磁薄膜的特定方向，以实现对斯格明子运动方向的控制。在模拟过程中，还需要考虑电流的施加方式，是直流电流还是脉冲电流等，不同的施加方式会对斯格明子的运动产生不同的影响。4.3模拟结果与分析4.3.1铁磁斯格明子在二极管结构中的运动特性通过微磁学模拟，成功获取了铁磁斯格明子在基于电压控制磁晶各向异性效应（VCMA）构建的二极管结构中的运动轨迹和速度变化情况，为深入理解其动力学行为提供了直观且关键的信息。在模拟过程中，设定初始条件下铁磁斯格明子位于二极管结构的特定位置，随后施加外部激励，包括合适的电流和磁场，以驱动斯格明子运动。从模拟得到的运动轨迹图（如图1所示）可以清晰地观察到，当斯格明子从低垂直磁晶各向异性（PMA）区域向高PMA边缘和梯形区域运动时，其运动轨迹呈现出明显的规律性。在正向运动时，即沿着有利于通过高PMA梯形区域的方向，斯格明子能够较为顺畅地通过该区域，其运动轨迹基本保持直线，且在高PMA边缘区域，斯格明子的运动轨迹更加稳定，没有出现明显的偏移或散射现象。这是因为高PMA边缘区域提供了额外的能量壁垒，使得斯格明子在该区域更加稳定，不易受到外界干扰而改变运动方向。[此处插入铁磁斯格明子在二极管结构中的运动轨迹图1]当斯格明子反向运动时，情况则截然不同。其运动轨迹在接近高PMA梯形区域时发生了明显的弯曲和变形，斯格明子难以顺利通过该区域，部分斯格明子甚至在梯形区域的边缘被阻挡或反射回来。这充分验证了高PMA梯形区域对斯格明子运动的单向限制作用，是实现二极管单向导通功能的关键机制之一。从能量角度分析，反向运动时斯格明子需要跨越更高的能量壁垒，而其自身的能量不足以支持其完成这一过程，因此被阻挡在梯形区域之外。进一步对斯格明子的速度变化进行分析，得到了斯格明子速度随时间和位置的变化曲线（如图2所示）。在正向运动过程中，当斯格明子进入高PMA边缘区域时，由于该区域的能量分布特性，使得斯格明子受到的阻力减小，速度逐渐增加。随着斯格明子继续向高PMA梯形区域运动，在梯形区域特殊的能量梯度和各向异性特性的作用下，斯格明子受到一个有助于其前进的力，速度进一步提升。这表明高PMA边缘和梯形区域不仅对斯格明子的运动方向具有调控作用，还能够影响其运动速度，提高斯格明子在二极管结构中的传输效率。[此处插入铁磁斯格明子速度随时间和位置的变化曲线图2]在反向运动时，斯格明子的速度在接近高PMA梯形区域时急剧下降。这是因为反向运动时斯格明子受到梯形区域的阻碍力作用，需要克服更高的能量势垒，导致其动能迅速减小，速度降低。当斯格明子被阻挡在梯形区域边缘时，速度趋近于零，无法继续前进。这种速度变化特性与运动轨迹的变化相互印证，进一步说明了基于铁磁斯格明子的二极管结构能够有效地实现单向导通功能，通过对斯格明子运动速度和方向的精确控制，模拟结果与理论预期相符，为该二极管结构的实际应用提供了有力的理论支持。4.3.2二极管的单向导通性能验证为了全面验证基于铁磁斯格明子的二极管的单向导通性能，通过模拟计算了不同方向电流下二极管的电流-电压（I-V）特性曲线，这些曲线直观地反映了二极管在不同条件下的导通状态和性能表现。在正向导通情况下，当施加正向电流时，从模拟得到的I-V特性曲线（如图3所示）可以看出，随着正向电压的逐渐增加，电流呈现出迅速上升的趋势。在低电压阶段，电流增长较为缓慢，这是因为此时铁磁斯格明子需要克服一定的能量势垒才能开始运动。随着电压的进一步升高，铁磁斯格明子获得足够的能量，能够顺利地通过二极管结构中的高垂直磁晶各向异性（PMA）梯形区域，形成有效的电流通路，电流迅速增大。当电压达到一定值后，电流增长趋于平缓，此时二极管进入了饱和导通状态。在饱和导通状态下，二极管的电阻较小，能够允许较大的电流通过，且电流对电压的变化不再敏感。[此处插入正向导通时的I-V特性曲线图3]在反向截止情况下，当施加反向电流时，I-V特性曲线表现出截然不同的特征。无论反向电压如何变化，电流始终保持在极低的水平，几乎可以忽略不计。这是因为在反向电流作用下，铁磁斯格明子受到高PMA梯形区域的强烈阻挡，无法通过该区域，从而无法形成有效的电流通路。反向电流主要由少数载流子（如热激发产生的电子-空穴对）形成，但由于少数载流子的浓度极低，所以反向电流非常小。这种正向导通、反向截止的I-V特性与传统二极管的特性相似，充分验证了基于铁磁斯格明子的二极管能够实现良好的单向导通功能。为了进一步评估二极管的单向导通性能，引入了整流比这一重要性能指标。整流比定义为正向导通电流与反向截止电流的比值，它反映了二极管单向导通性能的优劣。通过模拟数据计算得到，在特定的模拟条件下，该二极管的整流比高达10^5以上。这表明该二极管在正向导通时能够允许大量电流通过，而在反向截止时能够有效地抑制电流，具有非常高的单向导通性能。与传统二极管相比，基于铁磁斯格明子的二极管在整流比方面具有一定的优势，能够在一些对单向导通性能要求较高的应用场景中发挥重要作用。4.3.3不同因素对二极管性能的影响研究不同因素对基于铁磁斯格明子的二极管性能的影响，对于深入理解二极管的工作机制、优化器件性能以及拓展其应用领域具有重要意义。本研究主要从材料参数、结构参数和外场条件三个方面展开探讨。在材料参数方面，Dzyaloshinskii-Moriya（DM）相互作用强度对二极管性能有着显著影响。通过模拟不同DM相互作用强度下铁磁斯格明子在二极管结构中的运动和I-V特性，发现随着DM相互作用强度的增加，铁磁斯格明子的尺寸逐渐减小。这是因为较强的DM相互作用使得自旋之间的手性排列更加明显，斯格明子内部的自旋结构更加紧凑，从而导致其尺寸变小。较小尺寸的斯格明子在运动过程中受到的阻力相对较小，能够更快速地通过二极管结构，进而提高了二极管的导通速度和响应性能。较强的DM相互作用还增强了高垂直磁晶各向异性（PMA）梯形区域对斯格明子的单向限制作用，使得二极管的整流比得到进一步提高。当DM相互作用强度过强时，可能会导致斯格明子的稳定性下降，容易受到外界干扰而发生湮灭或结构变化，从而影响二极管的可靠性。饱和磁化强度也是一个重要的材料参数。模拟结果表明，饱和磁化强度的变化会影响铁磁斯格明子的运动速度和能量状态。当饱和磁化强度增大时，铁磁斯格明子的磁矩增大，其与外部磁场和其他能量项的相互作用增强。在相同的外部激励下，饱和磁化强度较大的斯格明子具有更高的运动速度，能够更快地通过二极管结构，提高了二极管的导通效率。饱和磁化强度的增大也会导致斯格明子的能量增加，在通过高PMA梯形区域时，需要克服更高的能量势垒。如果饱和磁化强度过大，可能会使斯格明子在反向运动时也能够获得足够的能量跨越梯形区域的能量壁垒，从而降低二极管的单向导通性能。在结构参数方面，高PMA梯形区域的尺寸对二极管性能起着关键作用。通过改变梯形区域的底边长度、斜边长度和高度等尺寸参数，模拟分析了铁磁斯格明子的运动特性和二极管的I-V特性。当梯形区域的底边长度增加时，斯格明子在通过该区域时需要跨越更长的距离，受到的阻力也相应增加。这可能导致斯格明子在正向运动时速度降低，影响二极管的导通速度。较长的底边长度也会增加梯形区域的能量梯度，增强对反向运动斯格明子的阻挡作用，提高二极管的整流比。斜边长度的变化会影响梯形区域的倾斜角度，进而影响斯格明子在其中的受力情况。适当调整斜边长度可以优化斯格明子在梯形区域的运动轨迹和速度，提高二极管的性能。梯形区域高度的增加会使能量势垒增大，对斯格明子的单向限制作用更加明显。这有利于提高二极管的整流比，但也可能会增加斯格明子在正向运动时的能量损耗，降低导通效率。因此，在设计二极管结构时，需要综合考虑梯形区域的各个尺寸参数，以实现最佳的性能。高PMA边缘区域的宽度也会对二极管性能产生影响。较宽的高PMA边缘区域能够提供更强大的能量壁垒，有效地避免铁磁斯格明子信号在运动过程中的湮灭，提高斯格明子的稳定性。宽边缘区域还可以进一步降低斯格明子在运动过程中的散射和干扰，使得斯格明子的运动更加稳定和有序，从而提高二极管的导通性能。过宽的高PMA边缘区域可能会增加二极管的电阻，影响其电流传输能力。在实际应用中，需要根据具体需求和材料特性，合理选择高PMA边缘区域的宽度。在外场条件方面，外加磁场对二极管性能有着重要影响。当施加平面内的外加磁场时，磁场的大小和方向会改变铁磁斯格明子的自旋结构和运动状态。随着外加磁场强度的增加，铁磁斯格明子的磁矩会逐渐沿着磁场方向旋转，其运动轨迹也会发生相应的变化。在正向导通情况下，适当的外加磁场可以与铁磁斯格明子的固有磁矩相互作用，增强斯格明子的运动能力，提高二极管的导通速度。但当磁场强度过大时，可能会破坏铁磁斯格明子的稳定性，导致其结构发生变化，甚至湮灭，从而影响二极管的正常工作。磁场的方向也会对二极管性能产生影响。当磁场方向与铁磁斯格明子的运动方向成一定角度时，会产生一个垂直于运动方向的洛伦兹力，使斯格明子的运动轨迹发生偏移。如果偏移过大，可能会导致斯格明子无法顺利通过二极管结构中的关键区域，影响二极管的导通性能。在设计基于铁磁斯格明子的二极管时，需要充分考虑外加磁场的大小和方向对器件性能的影响，通过合理控制外加磁场，优化二极管的性能。五、实验验证与分析5.1实验制备工艺基于铁磁斯格明子的二极管实验制备工艺是实现其功能和性能的关键环节，涉及到材料生长、微纳加工和结构组装等多个复杂且精细的步骤，每一步都对最终器件的性能有着至关重要的影响。在材料生长方面，选用具有强自旋-轨道耦合的重金属/铁磁薄膜/氧化物结构体系，其中铁磁薄膜材料为CoFeB，氧化物层为MgO，重金属层为Ta。采用磁控溅射技术进行材料生长，这是一种在高真空环境下，利用等离子体中的离子轰击靶材，使靶材原子溅射出来并沉积在衬底表面形成薄膜的方法。磁控溅射技术具有沉积速率高、薄膜质量好、成分易于控制等优点，能够精确控制各层薄膜的厚度和成分，满足基于铁磁斯格明子的二极管对材料性能的严格要求。在生长CoFeB薄膜时，通过精确控制溅射功率、气体流量和沉积时间等参数，实现对CoFeB薄膜厚度和质量的精确调控。溅射功率直接影响靶材原子的溅射速率，进而影响薄膜的生长速率和质量。气