基于第一性原理的二维材料磁性特性及调控机制研究

基于第一性原理的二维材料磁性特性及调控机制研究一、引言1.1研究背景与意义二维材料作为材料科学领域的前沿研究对象，近年来受到了广泛的关注。这类材料通常由单个或少数几个原子层组成，其独特的原子排列和电子结构赋予了它们许多新奇的物理特性，这些特性与传统的三维材料截然不同，为材料科学的发展开辟了新的道路。二维材料具有高载流子迁移率、高比表面积、优异的力学性能和独特的光学性质等特点，使其在电子学、能源、传感器、催化等众多领域展现出巨大的应用潜力。例如，在电子学领域，二维材料可用于制造高性能的晶体管、逻辑电路和传感器，有望推动集成电路向更小尺寸、更高性能的方向发展。在能源领域，二维材料可应用于太阳能电池、锂离子电池和超级电容器等，提高能量转换和存储效率。在传感器领域，二维材料对某些气体分子具有特殊的吸附和电学响应特性，可用于制备高灵敏度的气体传感器。在催化领域，二维材料的高比表面积和特殊的电子结构使其成为潜在的高效催化剂。在二维材料的诸多特性中，磁性是一个备受关注的研究方向。磁性在现代科技中扮演着至关重要的角色，从日常使用的电子设备到高端的信息技术领域，都离不开磁性材料的应用。对于二维材料而言，研究其磁性不仅有助于深入理解低维体系中的磁相互作用和量子物理现象，还为自旋电子学、量子计算等新兴领域的发展提供了关键的材料基础。在自旋电子学中，利用电子的自旋属性来存储和处理信息，需要具有特定磁性的材料作为核心元件。二维材料由于其原子级的厚度和独特的电子结构，能够提供长的自旋弛豫时间和自旋扩散长度，这对于实现高效的自旋输运和调控至关重要。如果能够在二维材料中实现室温铁磁性，将极大地推动自旋电子学器件的发展，使其能够在常温下稳定工作，从而广泛应用于计算机存储、通信等领域。在量子计算领域，二维磁性材料中的量子比特有望提供更高的计算速度和更低的能耗，为解决复杂的科学问题提供强大的计算能力。此外，二维材料的磁性还与其他物理性质，如电学、光学性质之间存在着密切的耦合关系，研究这些耦合效应有助于开发新型的多功能材料和器件，实现单一材料在多个领域的应用，从而提高材料的利用效率和器件的性能。然而，研究二维材料的磁性面临着诸多挑战。由于二维材料的原子层数少，表面原子的比例相对较高，表面效应和量子尺寸效应显著，这使得二维材料的磁性表现出与三维材料不同的特点，传统的磁性理论和研究方法在二维体系中往往不再适用。此外，二维材料的制备过程中容易引入缺陷和杂质，这些因素会对磁性产生复杂的影响，增加了研究的难度。实验上精确测量二维材料的磁性也面临着技术难题，由于样品的尺寸小、信号弱，需要高灵敏度的测量技术和设备。在这样的背景下，第一性原理计算方法成为研究二维材料磁性的重要手段。第一性原理计算基于量子力学原理，从电子和原子核的相互作用出发，无需引入过多的经验参数，能够直接计算材料的电子结构、能量和各种物理性质。通过第一性原理计算，可以深入了解二维材料中原子的排列方式、电子的分布状态以及磁相互作用的本质，为实验研究提供理论指导和预测。例如，计算不同原子掺杂或缺陷情况下二维材料的磁性变化，有助于设计具有特定磁性的二维材料；研究二维材料在外部电场、磁场作用下的磁性响应，为开发新型的磁电、磁光器件提供理论依据。第一性原理计算还可以模拟二维材料在不同制备条件下的结构和磁性，帮助优化制备工艺，提高材料的质量和性能。因此，开展二维材料磁性相关问题的第一性原理研究，对于揭示二维材料磁性的微观机制、探索新型二维磁性材料以及推动其在自旋电子学等领域的应用具有重要的科学意义和实际应用价值。1.2二维材料概述二维材料是指电子仅可在两个维度的非纳米尺度（1-100nm）上自由运动（平面运动）的材料，其原子排列方式呈现出与传统三维材料截然不同的特征。传统三维材料的原子在三维空间中周期性排列，原子间的相互作用在各个方向上都较为显著；而二维材料仅由单个或少数几个原子层组成，原子主要在二维平面内相互作用，层间相互作用相对较弱。这种独特的原子结构赋予了二维材料许多新颖的电子特性。从电子特性方面来看，二维材料中的电子态密度分布与三维材料有很大差异。在三维材料中，电子态密度是连续变化的，而在二维材料中，由于电子在一个方向上的运动受到限制，电子态密度呈现出阶梯状分布。以石墨烯为例，其电子具有线性色散关系，形成狄拉克锥结构，电子表现出类似无质量狄拉克费米子的行为，这使得石墨烯具有超高的载流子迁移率，在室温下可达10^5cm^2/(V·s)以上，远远超过传统半导体材料。二维过渡金属二硫化物（如MoS_2、WS_2等）具有半导体特性，其能带结构与块体材料相比发生了显著变化，且带隙随层数的减少而增大。这种独特的电子特性为二维材料在电子学领域的应用奠定了基础。二维材料还具有许多独特的物理化学性质。高载流子迁移率使得二维材料在高速电子器件中具有巨大的应用潜力，有望实现更快的运算速度和更低的功耗。二维材料通常具有大比表面积，这一特性使其在能源存储与转换、催化、传感器等领域表现出色。例如，在锂离子电池中，二维材料作为电极材料能够提供更多的活性位点，有利于锂离子的快速嵌入和脱出，从而提高电池的充放电性能和循环稳定性；在催化领域，大比表面积可以增加催化剂与反应物的接触面积，提高催化反应效率。二维材料还展现出优异的力学性能，如石墨烯具有极高的强度和柔韧性，能够承受较大的拉伸应力而不发生破裂，这使得二维材料在柔性电子器件中具有广阔的应用前景，可用于制备可弯曲、可折叠的电子设备。根据组成元素和结构的不同，二维材料可以分为多种类型。石墨烯是最早被发现和研究的二维材料，它由碳原子以sp^2杂化形成六角型蜂窝状晶格，具有优异的电学、力学和热学性能。由于其零带隙的特性，在电子学应用中需要通过一些手段（如掺杂、施加电场等）来打开带隙。过渡金属二硫化物（TMDCs）是另一类重要的二维材料，其化学式为MX_2（M为过渡金属，如Mo、W等；X为硫族元素，如S、Se、Te等），具有层状结构，层内原子通过共价键结合，层间通过范德华力相互作用。MoS_2在光电器件领域具有潜在应用价值，它可以吸收特定波长的光，产生光生载流子，可用于制备光电探测器、发光二极管等。六方氮化硼（h-BN）与石墨烯具有相似的结构，由硼原子和氮原子组成六角形蜂窝状晶格，但其具有高达5.2eV的宽带隙，是一种优良的绝缘体，具有高的热导率和化学稳定性，可作为二维材料异质结构中的绝缘层或衬底材料。黑磷是一种新型二维材料，具有类似褶皱的层状结构，其独特之处在于具有直接带隙，且带隙随层数变化可在0.3-2eV范围内调节，在光电子学和半导体器件领域具有重要的应用前景，可用于制造高性能的晶体管、光电探测器等。二维材料家族还包括过渡金属碳化物（MXenes）、二维金属有机框架（2D-MOFs）等，它们各自具有独特的结构和性能，在不同领域展现出潜在的应用价值。1.3第一性原理计算方法第一性原理计算方法是以量子力学为基础，从微观层面研究材料性质的重要手段。其核心在于通过求解薛定谔方程来描述材料中电子和原子核的相互作用，进而揭示材料的各种物理性质。薛定谔方程是量子力学的基本方程，它描述了微观粒子的波函数随时间和空间的演化。对于一个由N个电子和M个原子核组成的多粒子体系，其薛定谔方程可表示为：H\Psi(\mathbf{r}_1,\mathbf{r}_2,\cdots,\mathbf{r}_N,\mathbf{R}_1,\mathbf{R}_2,\cdots,\mathbf{R}_M,t)=i\hbar\frac{\partial\Psi(\mathbf{r}_1,\mathbf{r}_2,\cdots,\mathbf{r}_N,\mathbf{R}_1,\mathbf{R}_2,\cdots,\mathbf{R}_M,t)}{\partialt}其中，H是体系的哈密顿算符，包含电子的动能、电子与原子核之间的库仑相互作用能、电子之间的库仑相互作用能以及原子核之间的库仑相互作用能等项；\Psi是体系的波函数，它包含了体系中所有粒子的位置和自旋信息；\mathbf{r}_i和\mathbf{R}_j分别表示第i个电子和第j个原子核的位置矢量；\hbar是约化普朗克常数。然而，直接求解多粒子体系的薛定谔方程是极其困难的，因为其波函数是一个关于所有粒子坐标的高维函数。为了简化计算，人们引入了一系列近似方法，其中密度泛函理论（DFT）是目前第一性原理计算中应用最为广泛的理论框架。DFT的核心思想是将多电子体系的基态能量表示为电子密度的泛函，即体系的基态能量E可以写成电子密度\rho(\mathbf{r})的函数：E=E[\rho(\mathbf{r})]。根据Hohenberg-Kohn定理，对于一个给定的外部势场（由原子核产生），电子密度\rho(\mathbf{r})与体系的基态能量之间存在一一对应的关系，这意味着只要知道了电子密度，就可以确定体系的基态能量和其他物理性质。在DFT中，通常使用Kohn-Sham方程来求解电子密度。Kohn-Sham方程将多电子问题转化为一组单电子方程，通过引入一个等效的单电子势场（包含了电子之间的交换关联作用），使得每个电子在这个等效势场中独立运动。Kohn-Sham方程的形式为：\left[-\frac{\hbar^2}{2m}\nabla^2+V_{ext}(\mathbf{r})+V_{H}(\mathbf{r})+V_{xc}(\mathbf{r})\right]\psi_i(\mathbf{r})=\epsilon_i\psi_i(\mathbf{r})其中，-\frac{\hbar^2}{2m}\nabla^2是电子的动能算符；V_{ext}(\mathbf{r})是外部势场，即原子核产生的静电势；V_{H}(\mathbf{r})是Hartree势，表示电子之间的库仑相互作用的平均场；V_{xc}(\mathbf{r})是交换关联势，描述了电子之间的交换作用和关联作用；\psi_i(\mathbf{r})是第i个单电子波函数；\epsilon_i是对应的本征能量。电子密度\rho(\mathbf{r})可以通过单电子波函数表示为：\rho(\mathbf{r})=\sum_{i=1}^{N}\vert\psi_i(\mathbf{r})\vert^2在实际计算中，交换关联势V_{xc}(\mathbf{r})的精确形式是未知的，需要采用近似方法来处理。常用的近似方法有局域密度近似（LDA）和广义梯度近似（GGA）等。LDA假设交换关联能密度只与电子密度的局域值有关，即E_{xc}^{LDA}[\rho]=\int\epsilon_{xc}^{LDA}(\rho(\mathbf{r}))\rho(\mathbf{r})d\mathbf{r}，其中\epsilon_{xc}^{LDA}(\rho)是局域交换关联能密度，它是通过对均匀电子气的计算得到的。GGA则在LDA的基础上考虑了电子密度的梯度信息，认为交换关联能密度不仅与电子密度的局域值有关，还与电子密度的变化率有关，其交换关联能泛函的形式更为复杂。尽管这些近似方法在一定程度上能够准确描述材料的性质，但对于一些强关联体系或具有弱相互作用的体系，它们的计算结果可能存在较大误差，需要进一步发展更精确的交换关联泛函。第一性原理计算在研究材料的电子结构和磁性等性质方面具有诸多优势。由于其基于量子力学原理，从最基本的物理规律出发，无需引入过多的经验参数，能够直接计算材料的电子结构和各种物理性质，从而避免了因经验参数的不确定性而带来的误差。通过第一性原理计算，可以得到材料中电子的波函数、能量本征值、电子态密度等信息，进而深入了解材料的电子结构特征。在研究二维材料的磁性时，可以计算磁性原子的磁矩、磁相互作用能以及自旋极化等信息，揭示磁性的微观起源和磁相互作用机制。第一性原理计算能够深入到原子层面分析材料的性质。它可以精确描述原子的位置、原子间的相互作用以及电子在原子间的分布情况，从而对材料的结构和性质进行原子级别的精准控制。在设计新型二维磁性材料时，可以通过改变原子的种类、掺杂原子的位置和浓度等方式，利用第一性原理计算预测材料的磁性变化，为实验制备提供理论指导。第一性原理计算还可以模拟材料在不同外部条件下（如温度、压力、电场、磁场等）的性质变化，研究材料的性能与外部条件之间的关系，为材料的应用提供理论依据。二、二维材料磁性的理论基础2.1磁性的基本概念磁性起源于物质内部电子的运动，主要包括电子的自旋和轨道运动。电子的自旋是其固有属性，可看作是电子绕自身轴的旋转，产生自旋磁矩；电子绕原子核的轨道运动则产生轨道磁矩。原子的磁矩是电子轨道磁矩和自旋磁矩的矢量和，而物质的磁性则由大量原子磁矩的相互作用和排列方式决定。常见的磁性类型包括铁磁性、反铁磁性和亚铁磁性，它们各自具有独特的特点。铁磁性材料在居里温度以下，原子磁矩会自发地沿同一方向排列，形成磁畴。在没有外加磁场时，各磁畴的磁化方向随机分布，宏观上不显示磁性；当施加外加磁场时，磁畴会逐渐转向与外磁场方向一致，材料表现出很强的磁性，且在外磁场撤去后仍能保持一定的磁性，即具有剩磁。铁、钴、镍等金属及其合金是典型的铁磁性材料，广泛应用于变压器铁芯、电机、磁存储等领域。反铁磁性材料中，相邻原子的磁矩大小相等、方向相反，相互抵消，宏观上不显示磁性。例如，MnO是一种典型的反铁磁性材料，其磁矩在晶格中呈反平行排列。反铁磁性材料的磁性相对较弱，通常在特定温度（奈尔温度）以下才表现出反铁磁有序，当温度高于奈尔温度时，反铁磁性消失，材料转变为顺磁性。亚铁磁性材料中，不同原子的磁矩大小不相等，且呈反平行排列，虽然磁矩不能完全抵消，但仍具有一定的净磁矩，宏观上表现出磁性。常见的亚铁磁性材料有铁氧体，如Fe_3O_4，其化学式可写成Fe^{2+}Fe_2^{3+}O_4，其中Fe^{2+}和Fe^{3+}的磁矩反平行排列，但由于磁矩大小不同，存在净磁矩。亚铁磁性材料的磁性介于铁磁性和反铁磁性之间，具有较高的电阻率，常用于高频磁性器件。磁矩是描述磁性强弱和方向的物理量，它等于电流强度与电流回路所围面积的乘积再乘以单位法向量。对于单个原子，磁矩是电子轨道磁矩和自旋磁矩的矢量和，单位为玻尔磁子（\mu_B），其数值约为9.274\times10^{-24}J/T。在材料中，磁矩的分布和相互作用决定了材料的磁性行为。磁化强度是单位体积内的磁矩矢量和，用于描述材料被磁化的程度和方向，单位是安培每米（A/m）。磁化强度不仅反映了材料内部磁矩的总体取向，还与材料的磁性状态密切相关。当材料处于未磁化状态时，磁化强度为零；在磁化过程中，随着外加磁场的增加，磁化强度逐渐增大，当达到饱和磁化状态时，磁化强度达到最大值。居里温度是磁性材料的一个重要特征温度，对于铁磁性和亚铁磁性材料，当温度低于居里温度时，材料表现出铁磁性或亚铁磁性；当温度高于居里温度时，热运动使原子磁矩的取向变得无序，材料的磁性消失，转变为顺磁性。不同磁性材料的居里温度各不相同，例如，铁的居里温度约为1043K，钴的居里温度约为1394K，镍的居里温度约为627K。居里温度的高低反映了材料中磁相互作用的强弱，磁相互作用越强，居里温度越高。2.2二维材料中的磁性机制二维材料中磁性的产生涉及多种复杂机制，自旋-轨道耦合和超交换相互作用是其中较为关键的两种。自旋-轨道耦合是指电子的自旋与其轨道运动之间的相互作用。在二维材料中，由于原子层数少，电子的轨道和自旋自由度受到的限制与三维材料不同，使得自旋-轨道耦合效应在二维材料的磁性中扮演着重要角色。这种耦合作用可以导致电子的自旋方向与轨道运动方向相互关联，进而影响原子的磁矩和材料的磁性。以二维过渡金属二硫化物（TMDCs）为例，如MoS_2，其中的过渡金属原子（如Mo）具有较强的自旋-轨道耦合。在MoS_2的晶体结构中，Mo原子的d电子轨道与S原子的p电子轨道相互作用，这种相互作用使得电子的自旋-轨道耦合增强。自旋-轨道耦合导致了MoS_2能带结构的自旋劈裂，即在相同的能量下，电子的自旋向上和自旋向下状态具有不同的能量，这种自旋劈裂对MoS_2的磁性和电学性质产生了显著影响，为其在自旋电子学器件中的应用提供了基础。超交换相互作用是通过中间非磁性原子介导的磁性原子之间的间接磁相互作用。在二维材料中，当磁性原子通过非磁性原子（如氧、硫等）相互连接时，超交换相互作用起着重要作用。其原理基于磁性原子的电子与中间非磁性原子的电子之间的量子力学交换作用。以二维铁磁材料CrI_3为例，Cr原子是磁性原子，I原子为非磁性原子。Cr原子之间通过I原子发生超交换相互作用，使得相邻Cr原子的磁矩呈现出特定的排列方式，从而形成铁磁性。具体来说，Cr原子的3d电子与I原子的p电子发生杂化，通过这种杂化，Cr原子之间产生了间接的磁相互作用，这种相互作用使得Cr原子的磁矩能够在一定温度范围内保持平行排列，使CrI_3表现出铁磁性。超交换相互作用的强度和性质与磁性原子和中间非磁性原子的种类、原子间的距离以及晶体结构等因素密切相关。在不同的二维材料中，由于这些因素的差异，超交换相互作用的效果也会有所不同，从而导致材料呈现出不同的磁性特征。二维材料的磁性与电子结构之间存在着紧密的联系。电子结构决定了原子磁矩的大小和方向。原子的磁矩主要由未成对电子的自旋磁矩贡献，而电子在原子轨道中的分布情况决定了未成对电子的数量和自旋取向。在二维材料中，由于原子的排列方式和原子间的相互作用与三维材料不同，电子的轨道分布和自旋状态也会发生变化，进而影响原子磁矩。以石墨烯为例，本征石墨烯的碳原子以sp^2杂化形成六角型蜂窝状晶格，电子在这种晶格结构中形成了离域的π电子云，由于电子的自旋-轨道耦合作用较弱，且电子的自旋取向相互抵消，本征石墨烯表现为抗磁性。然而，当石墨烯中引入杂质原子（如氮、硼等）或存在缺陷时，会打破电子的自旋对称性，导致出现未成对电子，从而产生磁矩。电子的巡游性对二维材料的磁性也有重要影响。在一些二维材料中，电子具有一定的巡游性，即电子可以在晶格中自由移动。电子的巡游性会影响磁性原子之间的相互作用，使得磁相互作用不再局限于近邻原子之间，而是可以通过电子的传导在更大范围内发生。这种长程的磁相互作用对于材料的磁性有序状态的形成和稳定性具有重要意义。在二维铁磁材料中，电子的巡游性可以促进磁矩的平行排列，增强材料的铁磁性；而在反铁磁材料中，电子的巡游性可能会导致磁矩的反平行排列更加稳定。2.3影响二维材料磁性的因素原子结构是影响二维材料磁性的关键因素之一，其中原子种类和排列方式起着重要作用。不同种类的原子具有不同的电子构型和磁矩，这直接影响着二维材料的磁性。例如，过渡金属原子通常具有未填满的d电子壳层，这些未成对电子使得过渡金属原子具有较大的磁矩。在二维过渡金属硫族化合物（如MoS_2、WS_2等）中，过渡金属原子（如Mo、W）的磁矩对材料的磁性产生重要贡献。通过改变原子种类，如在MoS_2中用W原子部分替代Mo原子，会导致材料电子结构的变化，进而影响其磁性。由于W和Mo的电子构型和电负性存在差异，替代后会改变原子间的电荷分布和电子云重叠程度，使得自旋-轨道耦合和超交换相互作用发生变化，从而改变材料的磁性。原子的排列方式也显著影响二维材料的磁性。在二维材料中，原子的不同排列方式会导致不同的晶体结构和原子间距离，进而影响电子的运动和相互作用。以石墨烯为例，本征石墨烯中碳原子以规则的六角型蜂窝状晶格排列，电子的自旋-轨道耦合作用较弱，表现为抗磁性。然而，当石墨烯的晶格结构发生缺陷或变形时，如出现空位、边缘缺陷等，原子的排列方式被破坏，会导致电子的自旋对称性被打破，从而产生磁矩。在锯齿型边缘的石墨烯纳米带中，边缘原子的排列方式使得电子的自旋在边缘处发生特殊的耦合作用，导致纳米带出现磁性，且磁矩的大小和方向与纳米带的宽度和边缘结构密切相关。晶体结构对二维材料的磁性同样有着重要影响，其中晶格常数和对称性是两个关键方面。晶格常数决定了原子间的距离，而原子间距离的变化会对磁相互作用产生显著影响。在二维铁磁材料CrI_3中，晶格常数的改变会影响Cr原子之间的超交换相互作用。当晶格常数增大时，Cr原子之间的距离增大，超交换相互作用减弱，使得磁矩之间的耦合作用变弱，材料的磁性可能会发生变化，如居里温度降低；反之，当晶格常数减小时，超交换相互作用增强，磁矩之间的耦合作用变强，材料的磁性可能增强，居里温度升高。晶体的对称性也对二维材料的磁性有着深刻影响。不同的晶体对称性会导致电子在晶体中的运动状态和相互作用方式不同，从而影响材料的磁性。具有高对称性的晶体结构，如六方晶系的CrI_3，其原子排列的对称性使得电子的轨道和自旋分布相对较为规则，磁相互作用具有一定的规律性。在这种情况下，材料的磁性表现出较好的各向同性，即磁性在各个方向上的性质较为相似。而当晶体对称性降低时，如通过引入缺陷或施加外部应力等方式破坏晶体的对称性，会导致电子的轨道和自旋分布发生变化，磁相互作用的各向异性增强。在一些二维材料中，当晶体对称性被破坏后，会出现磁各向异性，使得材料在不同方向上的磁性性质不同，这对于材料在自旋电子学等领域的应用具有重要意义，可以通过控制磁各向异性来实现对自旋的有效调控。外部条件对二维材料的磁性也有着不可忽视的影响，温度、压力和磁场是常见的外部因素。温度对二维材料磁性的影响主要体现在热运动对原子磁矩排列的干扰上。随着温度的升高，原子的热运动加剧，原子磁矩的取向变得更加无序，这会削弱材料的磁性。对于铁磁性二维材料，当温度升高到居里温度时，热运动足以破坏原子磁矩的自发有序排列，材料的铁磁性消失，转变为顺磁性。以二维铁磁材料Fe_3GeTe_2为例，其居里温度约为220K，当温度低于居里温度时，材料表现出铁磁性，原子磁矩呈有序排列；当温度高于居里温度时，原子磁矩的有序排列被破坏，材料转变为顺磁性。压力作为一种外部条件，能够改变二维材料的晶格结构和原子间距离，从而对磁性产生影响。当对二维材料施加压力时，晶格会发生压缩或畸变，原子间的相互作用增强或改变。在一些二维过渡金属氧化物中，施加压力会使原子间距离减小，电子云的重叠程度增加，导致自旋-轨道耦合和超交换相互作用发生变化，进而改变材料的磁性。研究发现，在一定压力范围内，某些二维过渡金属氧化物的磁矩会随着压力的增加而增大，这是由于压力导致原子间的磁相互作用增强，使得磁矩更容易有序排列。磁场是调控二维材料磁性的重要外部手段。当二维材料处于外加磁场中时，原子磁矩会受到磁场的作用而发生取向变化。对于顺磁性二维材料，在外加磁场的作用下，原子磁矩会顺着磁场方向排列，材料表现出磁性，且磁化强度随着磁场强度的增加而增大。对于铁磁性二维材料，外加磁场可以使磁畴的取向发生改变，从而影响材料的磁性。在弱磁场下，磁畴会逐渐转向与外磁场方向一致，材料的磁化强度逐渐增加；当磁场强度达到一定程度时，所有磁畴都沿外磁场方向排列，材料达到饱和磁化状态。磁场还可以诱导二维材料产生一些特殊的磁性现象，如磁致伸缩效应，即材料在磁场作用下发生尺寸变化，这种效应在一些二维磁性材料中被观察到，对于开发新型磁传感器和磁驱动器件具有潜在的应用价值。三、基于第一性原理的二维材料磁性研究方法3.1第一性原理计算软件及工具在二维材料磁性的第一性原理研究中，VASP和CASTEP是两款广泛应用的计算软件，它们各自具有独特的功能特点、适用范围和优势。VASP（ViennaAb-initioSimulationPackage）是一款基于密度泛函理论的计算软件，在材料科学领域应用极为广泛。它能够精确地计算材料的电子结构、能量和各种物理性质。在计算效率方面，VASP采用了平面波赝势方法（PWPM）和投影缀加波（PAW）方法，这些方法能够有效地处理电子与原子核之间的相互作用，从而提高计算效率。在研究二维过渡金属硫族化合物（TMDCs）的磁性时，VASP可以快速准确地计算出材料的电子结构和磁矩。VASP还具备处理复杂体系的能力，能够对包含缺陷、杂质或施加外场的二维材料体系进行计算。在研究掺杂对二维材料磁性的影响时，VASP可以精确地模拟掺杂原子与基体原子之间的相互作用，从而深入分析掺杂对磁性的影响机制。CASTEP（CambridgeSequentialTotalEnergyPackage）是另一款重要的第一性原理计算软件，它基于平面波赝势方法，在材料模拟方面具有强大的功能。CASTEP的一个显著优势是其图形化界面，这使得操作更加直观和便捷，对于初学者来说更容易上手。在构建二维材料的晶体结构模型时，用户可以通过CASTEP的图形界面直接进行原子的添加、删除和位置调整等操作。CASTEP在计算材料的晶体结构和电子性质方面表现出色。它可以精确地优化二维材料的晶格参数和原子坐标，得到稳定的晶体结构。在研究二维材料的电子结构时，CASTEP能够准确地计算出能带结构、态密度等信息，为分析材料的电学和磁性性质提供重要依据。在使用VASP进行二维材料磁性计算时，有一些经验和注意事项。合理设置计算参数至关重要。例如，在计算磁性体系时，需要正确设置自旋极化相关的参数，如ISPIN（自旋极化标识）、MAGMOM（磁矩设置）等。对于铁磁性二维材料，通常将ISPIN设置为2，表示开启自旋极化计算，并根据具体情况合理设置MAGMOM的值，以准确描述原子的磁矩。选择合适的赝势也非常关键。赝势的质量会直接影响计算结果的准确性。对于过渡金属元素，需要选择能够准确描述其电子结构的赝势。在研究二维过渡金属氧化物时，应选择考虑了过渡金属d电子特性的赝势，以确保计算结果能够反映材料的真实磁性。使用CASTEP进行二维材料磁性计算时，同样需要注意一些问题。在构建模型时，要确保模型的合理性和准确性。对于二维材料，需要正确设置周期性边界条件，以模拟材料的无限周期性。在研究二维材料的表面吸附对磁性的影响时，要合理构建吸附模型，包括选择合适的吸附原子、吸附位置和吸附构型等。在分析计算结果时，要充分利用CASTEP提供的各种可视化工具，如能带结构、态密度的可视化等。通过这些可视化工具，可以更直观地理解材料的电子结构和磁性性质，从而深入分析计算结果。3.2计算模型的构建在基于第一性原理研究二维材料的磁性时，构建合理的计算模型是获得准确结果的关键步骤。这一过程需要依据二维材料的晶体结构进行细致的设计，充分考虑多种因素对计算结果的影响。超胞模型是研究二维材料磁性常用的计算模型之一。超胞是在原胞的基础上，通过一定的扩展方式构建而成的更大的晶胞。选择合适的超胞模型对于准确描述二维材料的磁性至关重要。在构建超胞模型时，需要考虑材料的对称性和周期性。以二维过渡金属硫族化合物（TMDCs）为例，其晶体结构通常具有六方对称性。在构建超胞模型时，应保持这种对称性，以确保计算结果的准确性。对于MoS_2，其原胞包含一个Mo原子和两个S原子，在构建超胞时，可以通过在二维平面内重复原胞的方式来扩大晶胞尺寸。常见的超胞扩展方式有2×2×1、3×3×1等，其中前面两个数字表示在二维平面内的扩展倍数，最后一个数字表示在垂直于二维平面方向的扩展倍数。不同的超胞扩展倍数会对计算结果产生影响。较小的超胞扩展倍数计算量相对较小，但可能无法充分考虑原子间的长程相互作用；较大的超胞扩展倍数虽然能更准确地描述原子间的相互作用，但计算量会大幅增加。因此，需要根据具体研究目的和计算资源，选择合适的超胞扩展倍数。在研究MoS_2中单个缺陷对磁性的影响时，2×2×1的超胞模型可能就足够了；而在研究缺陷间的相互作用时，则可能需要更大的超胞模型，如3×3×1或4×4×1。原子坐标的设置也是计算模型构建中的重要环节。准确的原子坐标能够保证计算模型的正确性，进而影响计算结果的可靠性。在设置原子坐标时，需要参考实验数据或已有的理论研究成果。对于一些常见的二维材料，如石墨烯、MoS_2等，其原子坐标在文献中已有详细报道，可以直接引用。对于新合成的二维材料或结构较为复杂的二维材料，可能需要通过实验手段（如X射线衍射、透射电子显微镜等）来确定原子坐标。在设置原子坐标时，还需要注意坐标的精度。一般来说，原子坐标的精度应达到一定的量级，以确保计算结果的准确性。在VASP计算中，原子坐标的精度通常设置为小数点后6位或更高。如果原子坐标的精度过低，可能会导致计算结果出现偏差，例如，在计算二维材料的晶格常数时，原子坐标的误差可能会导致晶格常数的计算值与实验值不符。周期性边界条件是计算模型构建中必须考虑的因素之一。二维材料在二维平面内具有无限的周期性，为了模拟这种周期性，需要在计算模型中施加周期性边界条件。周期性边界条件的原理是将超胞在各个方向上进行无限重复，使得超胞中的原子与相邻超胞中的原子相互作用。在VASP和CASTEP等计算软件中，都提供了设置周期性边界条件的功能。通过设置周期性边界条件，可以有效地减少计算量，同时保证计算结果能够反映二维材料的真实性质。在研究二维材料的电子结构时，周期性边界条件能够使电子在整个二维平面内自由运动，从而得到准确的电子态密度和能带结构。如果不施加周期性边界条件，电子将被限制在超胞内，无法体现二维材料的周期性特征，导致计算结果与实际情况不符。真空层的设置对于二维材料的计算也非常重要。由于二维材料只有原子级的厚度，为了避免相邻超胞之间的相互作用，需要在垂直于二维平面的方向上设置一定厚度的真空层。真空层的厚度需要根据具体情况进行合理选择。如果真空层厚度过小，相邻超胞之间可能会产生相互作用，影响计算结果的准确性；如果真空层厚度过大，会增加计算量，延长计算时间。一般来说，真空层的厚度应设置在10-20Å之间。对于一些具有较大层间相互作用的二维材料，如双层石墨烯，可能需要适当增加真空层的厚度，以确保相邻层之间的相互作用可以忽略不计。在研究二维材料的吸附性质时，真空层的厚度也会影响吸附原子与二维材料之间的相互作用。如果真空层过薄，吸附原子可能会与相邻超胞中的原子发生相互作用，导致计算得到的吸附能不准确。3.3计算参数的选择与优化在基于第一性原理的二维材料磁性研究中，计算参数的选择与优化至关重要，它直接影响到计算结果的准确性和可靠性。交换关联泛函、平面波截断能、k点网格密度等参数在计算过程中起着关键作用，下面将对这些参数进行详细分析。交换关联泛函是密度泛函理论中的核心部分，它描述了电子之间的交换和关联相互作用。目前常用的交换关联泛函有局域密度近似（LDA）和广义梯度近似（GGA）等。LDA假设交换关联能只与电子密度的局域值有关，它在处理一些简单体系时能够取得较好的结果。对于金属体系，LDA能够较为准确地描述电子的行为，因为金属中的电子云分布相对较为均匀。然而，对于具有复杂电子结构的体系，如过渡金属化合物，LDA的计算结果可能会出现偏差。这是因为LDA忽略了电子密度的梯度信息，而在过渡金属化合物中，电子密度的变化较为剧烈，这种忽略会导致对电子间相互作用的描述不够准确。GGA在LDA的基础上考虑了电子密度的梯度信息，认为交换关联能不仅与电子密度的局域值有关，还与电子密度的变化率有关。因此，GGA在处理具有复杂电子结构的体系时表现更为出色。在研究二维过渡金属硫族化合物（TMDCs）的磁性时，GGA能够更准确地描述过渡金属原子与硫族原子之间的电子相互作用，从而得到更符合实际情况的磁性性质。对于MoS_2，GGA计算得到的能带结构和磁矩与实验结果更为接近。不同的GGA泛函，如PBE（Perdew-Burke-Ernzerhof）、PW91（Perdew-Wang1991）等，在计算结果上也会存在一定差异。PBE泛函在计算二维材料的晶格常数和电子结构时，通常能够给出与实验值较为吻合的结果，因此在二维材料的研究中应用较为广泛。在选择交换关联泛函时，需要根据具体的研究体系和目的进行权衡，选择最适合的泛函以获得准确的计算结果。平面波截断能是平面波赝势方法中的一个重要参数，它决定了平面波基组的大小。在计算过程中，电子的波函数用平面波展开，平面波截断能越高，参与展开的平面波数量就越多，对电子波函数的描述就越精确。然而，过高的平面波截断能会导致计算量急剧增加，计算时间大幅延长。因此，需要在计算精度和计算效率之间找到平衡。在研究二维材料时，通常需要通过测试来确定合适的平面波截断能。以石墨烯为例，一般情况下，平面波截断能在400-500eV之间能够得到较为准确的结果。当平面波截断能低于400eV时，对电子波函数的描述不够精确，会导致计算得到的晶格常数、电子结构等结果与实际情况存在较大偏差。而当平面波截断能高于500eV时，虽然计算精度会有所提高，但计算量的增加会使计算成本大幅上升，且精度提升的幅度并不明显。对于不同的二维材料，由于其原子结构和电子结构的差异，合适的平面波截断能也会有所不同。在研究二维过渡金属氧化物时，由于过渡金属原子的电子结构较为复杂，可能需要更高的平面波截断能才能准确描述其电子行为，一般建议在500-600eV之间进行测试和选择。k点网格密度用于描述在倒易空间中对布里渊区的采样密度。k点网格越密，对布里渊区的采样就越精细，计算结果就越准确。但同样地，过密的k点网格会增加计算量。在二维材料的计算中，通常采用Monkhorst-Pack方法来生成k点网格。对于二维材料，由于其在二维平面内具有周期性，而在垂直方向上的周期性较弱（通常通过设置真空层来模拟），因此在k点网格的设置上，二维平面内的k点密度需要相对较高，而垂直方向上的k点密度可以适当降低。在研究单层MoS_2时，在二维平面内采用3×3×1的k点网格能够较好地平衡计算精度和计算效率。当k点网格过疏，如采用2×2×1的k点网格时，会导致对布里渊区的采样不足，计算得到的能带结构和态密度等结果会出现明显的误差，无法准确反映材料的电子性质。而当k点网格过密，如采用5×5×1的k点网格时，虽然计算精度会提高，但计算量会显著增加，计算时间也会大幅延长。对于不同尺寸的超胞模型，k点网格密度也需要进行相应的调整。随着超胞尺寸的增大，为了保证计算精度，k点网格需要适当加密。对于一个较大尺寸的MoS_2超胞模型（如4×4×1），可能需要采用4×4×1的k点网格来获得准确的计算结果。为了更直观地说明参数优化的重要性，下面通过对比不同参数下的计算结果进行分析。以二维铁磁材料CrI_3为例，分别采用LDA和GGA（PBE）泛函进行计算。在LDA计算中，得到的CrI_3的晶格常数与实验值相比存在一定偏差，且计算得到的磁矩也与实验值有较大差异。而采用GGA（PBE）泛函计算时，晶格常数和磁矩与实验值更为接近，能够更准确地描述CrI_3的结构和磁性。在平面波截断能的对比中，当平面波截断能为350eV时，计算得到的CrI_3的能带结构出现明显的畸变，无法准确反映其电子结构；当平面波截断能提高到450eV时，能带结构得到明显改善，能够更准确地描述电子的能量状态。在k点网格密度的对比中，采用2×2×1的k点网格计算时，CrI_3的态密度曲线出现较多噪声，无法准确反映电子态的分布情况；而采用3×3×1的k点网格计算时，态密度曲线更加平滑，能够准确地反映电子态的分布。这些对比结果表明，合理选择和优化计算参数能够显著提高计算结果的准确性，从而为二维材料磁性的研究提供可靠的理论依据。在进行二维材料磁性的第一性原理计算时，必须充分重视计算参数的选择与优化，以获得准确、可靠的计算结果。3.4计算结果的分析与验证在完成二维材料磁性的第一性原理计算后，对计算结果进行深入分析与验证是确保研究可靠性和有效性的关键环节。通过对电子结构、磁矩分布和能带结构等结果的分析，能够揭示二维材料磁性的微观机制和物理本质。对于电子结构的分析，主要关注电子的分布状态和能级结构。通过计算得到的电子密度分布，可以直观地了解电子在原子周围的聚集情况。在二维过渡金属硫族化合物（TMDCs）中，电子在过渡金属原子和硫族原子周围的分布存在明显差异。在MoS_2中，Mo原子周围的电子密度相对较高，这是由于Mo原子的d电子轨道与S原子的p电子轨道相互作用，使得电子在Mo原子周围聚集。这种电子分布的差异对材料的磁性产生重要影响，因为电子的自旋和轨道运动与电子密度密切相关。分析能级结构可以了解电子的能量状态。通过计算得到的能带结构和态密度图，可以确定材料的价带、导带以及能隙大小。在二维铁磁材料中，能隙的存在与否以及大小对材料的磁性和电学性质有着重要影响。如果能隙较小，电子容易在价带和导带之间跃迁，可能导致材料的磁性不稳定；而较大的能隙则有助于保持材料的磁性。磁矩分布是分析二维材料磁性的重要指标。通过计算得到的磁矩分布，可以了解磁性原子的磁矩大小和方向。在二维铁磁材料中，磁性原子的磁矩通常呈现出平行排列的状态，从而形成宏观的铁磁性。在CrI_3中，Cr原子的磁矩沿同一方向排列，使得材料表现出铁磁性。分析磁矩分布还可以研究磁相互作用的范围和强度。通过计算不同原子之间的磁矩耦合，可以确定磁相互作用是短程的还是长程的，以及相互作用的强度大小。在一些二维材料中，磁相互作用可能只在近邻原子之间发生，而在另一些材料中，磁相互作用可能通过电子的巡游性在较大范围内传播。能带结构是反映材料电子能量与波矢关系的重要物理量，对分析二维材料的磁性和电学性质具有重要意义。通过分析能带结构，可以了解电子在晶体中的运动状态和能量分布。在二维材料中，能带结构的特点与材料的晶体结构和原子间相互作用密切相关。对于具有六角形蜂窝状晶格结构的石墨烯，其能带结构在K点附近形成狄拉克锥，电子具有线性色散关系，这使得石墨烯具有独特的电学性质。而在二维过渡金属氧化物中，由于过渡金属原子的d电子与氧原子的p电子之间的强相互作用，能带结构会出现明显的分裂和展宽，这对材料的磁性和电学性质产生重要影响。分析能带结构还可以确定材料的磁性转变温度。通过计算不同温度下的能带结构，可以观察到随着温度的升高，能带结构的变化情况。当温度接近居里温度时，能带结构会发生明显的变化，这是由于热运动导致原子磁矩的无序排列，从而影响了电子的能量状态。为了验证第一性原理计算结果的准确性，通常需要将计算结果与实验结果或其他理论计算结果进行对比。与实验结果对比是最直接的验证方法。在研究二维材料的磁性时，可以将计算得到的磁矩、居里温度、磁化强度等物理量与实验测量值进行比较。在研究二维铁磁材料Fe_3GeTe_2时，通过第一性原理计算得到的居里温度与实验测量值进行对比。如果计算值与实验值相符，说明计算方法和模型是可靠的；如果存在偏差，则需要进一步分析原因。计算值与实验值存在差异可能是由于计算模型的简化、交换关联泛函的近似性、实验测量误差等多种因素导致的。在计算模型方面，可能忽略了一些实际存在的因素，如杂质、缺陷等，这些因素可能会对材料的磁性产生重要影响。在交换关联泛函方面，目前常用的近似方法可能无法准确描述电子之间的相互作用，从而导致计算结果与实际情况存在偏差。为了解决这些问题，可以采用更精确的计算模型和交换关联泛函，或者结合实验测量结果对计算模型进行修正。与其他理论计算结果对比也是验证计算结果准确性的重要方法。不同的理论计算方法可能采用不同的近似和假设，通过对比不同方法的计算结果，可以更全面地了解材料的性质。在研究二维材料的磁性时，可以将基于密度泛函理论的第一性原理计算结果与基于紧束缚模型、海森堡模型等其他理论方法的计算结果进行对比。如果不同方法的计算结果相互吻合，说明计算结果具有较高的可靠性；如果存在差异，则需要分析差异产生的原因。不同方法的计算结果存在差异可能是由于方法本身的局限性、模型参数的选择等因素导致的。在紧束缚模型中，对电子相互作用的描述相对简单，可能无法准确反映复杂的电子结构；而在海森堡模型中，主要关注磁矩之间的相互作用，对电子的具体行为描述不够详细。在进行对比时，需要充分考虑不同方法的特点和适用范围，以便更准确地评估计算结果的准确性。四、二维材料磁性的第一性原理研究实例4.1石墨烯基二维材料的磁性研究4.1.1石墨烯的本征磁性及调控本征石墨烯由碳原子以sp^2杂化形成六角型蜂窝状晶格，电子在这种高度对称且均匀的结构中形成离域的π电子云。由于电子的自旋-轨道耦合作用极其微弱，且电子的自旋取向在整体上相互抵消，使得本征石墨烯不具备磁性，表现为抗磁性。这种抗磁性源于电子在外加磁场中的拉莫尔进动，产生与外加磁场方向相反的感应磁矩，但这种抗磁性非常微弱，在实际应用中难以被利用。为了赋予石墨烯磁性，科研人员开展了大量研究，其中掺杂和引入缺陷是两种主要的调控方式。在掺杂方面，氮、硼等原子的掺杂备受关注。以氮掺杂为例，氮原子的外层电子结构与碳原子不同，当氮原子替代石墨烯中的碳原子时，会引入额外的电子。第一性原理计算表明，氮原子的孤对电子会打破石墨烯原有的电子自旋对称性，从而产生局域磁矩。研究发现，随着氮掺杂浓度的增加，石墨烯体系的总磁矩呈现出先增大后减小的趋势。当氮掺杂浓度较低时，每个氮原子引入的局域磁矩相互独立，体系总磁矩随着氮原子数量的增加而增大；然而，当氮掺杂浓度过高时，氮原子之间的相互作用增强，导致部分磁矩相互抵消，体系总磁矩反而减小。不同的氮掺杂位置也会对磁性产生显著影响。在石墨烯的六角晶格中，氮原子处于不同的位置时，与周围碳原子的电子云重叠程度不同，进而影响磁矩的大小和方向。理论计算预测，当氮原子处于特定的晶格位置时，能够形成稳定的铁磁耦合，使得石墨烯体系表现出宏观的铁磁性。硼掺杂同样能够改变石墨烯的磁性。硼原子的价电子数比碳原子少一个，掺杂后会在石墨烯中产生空穴。这些空穴会改变电子的分布状态，导致电子自旋的重新排列，从而产生磁矩。与氮掺杂不同的是，硼掺杂诱导的磁矩与氮掺杂诱导的磁矩在大小和方向上可能存在差异。通过第一性原理计算可以深入分析硼掺杂石墨烯的电子结构变化，发现硼原子周围的电子云密度降低，空穴的存在使得电子的自旋极化增强，从而产生磁矩。在研究硼掺杂对石墨烯磁性的影响时，还需要考虑硼原子与其他杂质或缺陷的协同作用。在实际制备过程中，很难避免引入其他杂质或缺陷，这些因素可能会与硼掺杂相互作用，进一步改变石墨烯的磁性。引入缺陷也是调控石墨烯磁性的有效手段。空位是一种常见的缺陷类型，当石墨烯晶格中的碳原子缺失时，会形成空位缺陷。第一性原理计算显示，空位周围的碳原子会发生重构，以降低体系的能量。这种重构导致电子云分布的改变，使得空位周围的碳原子产生未成对电子，从而形成局域磁矩。在具有单个空位的石墨烯体系中，空位周围的几个碳原子会表现出明显的自旋极化，形成稳定的磁矩。研究还发现，空位缺陷的浓度和分布对石墨烯的磁性有重要影响。当空位浓度较低时，各个空位产生的磁矩相互独立，体系的磁性较弱；随着空位浓度的增加，空位之间的相互作用增强，磁矩之间可能发生耦合，使得体系的磁性增强。如果空位分布不均匀，会导致磁性的各向异性，即在不同方向上石墨烯的磁性表现不同。边缘缺陷也是影响石墨烯磁性的重要因素。在石墨烯纳米带中，锯齿型边缘和扶手椅型边缘具有不同的电子结构和磁性特征。锯齿型边缘的石墨烯纳米带由于边缘碳原子的悬键，容易产生局域磁矩。第一性原理计算表明，锯齿型边缘的碳原子形成了特殊的自旋排列，使得纳米带在边缘处呈现出铁磁性。而且，纳米带的宽度和边缘的化学修饰对磁性也有显著影响。随着纳米带宽度的增加，边缘磁矩之间的相互作用增强，可能导致磁性的增强或改变。对锯齿型边缘进行化学修饰，如加氢或加氟等，可以改变边缘碳原子的电子结构，从而调控磁性。加氢修饰可以饱和边缘碳原子的悬键，减少未成对电子的数量，降低磁矩；而加氟修饰可能会引入额外的电子或空穴，改变电子的自旋状态，从而影响磁性。扶手椅型边缘的石墨烯纳米带通常磁性较弱，但通过引入特定的缺陷或掺杂原子，也可以诱导出磁性。在扶手椅型边缘引入空位缺陷或氮掺杂原子，能够打破边缘的电子对称性，产生局域磁矩。通过第一性原理计算可以精确分析不同修饰和缺陷情况下扶手椅型边缘石墨烯纳米带的电子结构和磁性变化，为设计具有特定磁性的石墨烯基材料提供理论指导。4.1.2石墨烯与其他材料复合体系的磁性石墨烯与过渡金属复合形成的二维体系展现出独特的磁性特征，这主要源于复合材料界面处强烈的电子相互作用。以石墨烯与铁（Fe）复合体系为例，第一性原理计算表明，在界面处，石墨烯的π电子与Fe原子的d电子发生杂化。这种杂化作用使得电子云的分布发生改变，电子的自旋状态也相应调整。由于Fe原子具有较大的磁矩，其磁矩通过界面处的电子杂化作用传递给石墨烯，从而诱导石墨烯产生磁性。计算结果显示，在石墨烯/Fe复合体系中，石墨烯靠近Fe的一侧电子自旋极化增强，形成了一定的磁矩。这种磁矩的大小和方向与Fe原子的磁矩以及界面处的电子杂化强度密切相关。当Fe原子在石墨烯表面的吸附方式不同时，界面处的电子杂化程度也会有所差异，进而导致石墨烯磁矩的变化。如果Fe原子以特定的吸附构型在石墨烯表面形成有序排列，能够增强电子杂化作用，使得石墨烯的磁矩增大。在石墨烯与磁性氧化物复合体系中，如石墨烯与氧化铁（Fe_2O_3）复合，界面处的电子相互作用同样对磁性产生重要影响。Fe_2O_3具有一定的磁性，其晶体结构中的Fe原子磁矩在一定条件下会呈现出有序排列。当石墨烯与Fe_2O_3复合时，界面处的电子会发生重新分布。由于石墨烯的高导电性和特殊的电子结构，它能够与Fe_2O_3中的电子发生相互作用，影响Fe_2O_3的电子态和磁矩。通过第一性原理计算发现，在石墨烯/Fe_2O_3复合体系中，界面处的电子转移使得Fe_2O_3的磁矩发生改变，同时石墨烯也被诱导产生一定的磁矩。这种相互作用还会导致复合体系的电子能带结构发生变化，出现新的电子态。这些新的电子态对复合体系的磁性和电学性质都有重要影响，例如，可能会改变复合体系的居里温度和磁导率。这些石墨烯基复合体系在多个领域展现出潜在的应用价值。在磁性传感器领域，利用石墨烯与过渡金属或磁性氧化物复合体系的磁性对外部磁场的敏感响应，可以制备高灵敏度的磁性传感器。由于石墨烯具有高载流子迁移率和大比表面积，能够快速传递电子信号，与磁性材料复合后，能够增强对磁场的感知能力。当外部磁场发生变化时，复合体系的磁矩会相应改变，从而引起电阻或电容等电学性质的变化，通过检测这些电学信号的变化，就可以实现对磁场的精确测量。在自旋电子学器件方面，石墨烯基复合体系的独特磁性和电学性质使其有望成为构建新型自旋电子学器件的关键材料。在自旋过滤器中，利用复合体系中电子的自旋极化特性，可以实现对特定自旋方向电子的选择性传输，从而提高自旋电子学器件的性能。在自旋存储器件中，通过调控复合体系的磁性状态，可以实现信息的存储和读取，有望提高存储密度和读写速度。4.2过渡金属二硫化物（TMDs）的磁性研究4.2.1TMDs的晶体结构与磁性关系过渡金属二硫化物（TMDs）具有典型的MX₂型结构，其中M代表过渡金属，如钼（Mo）、钨（W）等；X为硫族元素，常见的有硫（S）、硒（Se）等。这种结构呈现出层状特征，每一层由一个过渡金属原子层夹在两个硫族原子层之间，通过共价键紧密结合。层与层之间则通过较弱的范德华力相互作用，使得TMDs可以容易地被剥离成单层或少数层结构。在TMDs中，常见的晶体结构包括2H相和1T相。2H相是一种六方晶系结构，其原子排列具有高度的对称性。在2H-MoS₂中，钼原子位于六方晶格的中心，硫原子则分布在钼原子的上下两侧，形成三棱柱配位结构。这种结构使得电子在晶体中的分布较为均匀，自旋-轨道耦合作用相对较弱。从磁性角度来看，2H相的TMDs通常表现出较弱的磁性或非磁性。这是因为在这种结构中，原子磁矩之间的相互作用较弱，难以形成有序的磁排列。通过第一性原理计算发现，2H-MoS₂的总磁矩接近于零，其电子结构中的自旋向上和自旋向下状态几乎简并，磁相互作用能较小。1T相则属于四方晶系结构，与2H相相比，其原子排列发生了明显的变化。在1T-MoS₂中，钼原子的配位环境从三棱柱变为八面体，这种结构变化导致电子云的分布发生改变，自旋-轨道耦合作用增强。由于结构的改变，1T相的TMDs中原子磁矩之间的相互作用增强，更容易形成磁有序状态。第一性原理计算表明，1T-MoS₂具有一定的固有磁矩，相邻钼原子之间存在反铁磁相互作用。这种反铁磁相互作用是通过硫原子作为中间媒介，发生超交换相互作用实现的。在1T相中，由于钼原子的d电子与硫原子的p电子之间的杂化作用更强，使得超交换相互作用得以增强，从而稳定了反铁磁结构。不同晶体结构对TMDs磁性的影响机制主要源于原子间的距离、电子云的重叠程度以及自旋-轨道耦合强度的变化。在2H相中，原子间的距离和电子云重叠程度相对较小，自旋-轨道耦合作用较弱，导致磁相互作用较弱。而在1T相中，原子间的距离和电子云重叠程度发生改变，自旋-轨道耦合作用增强，使得磁相互作用增强，从而产生不同的磁性表现。研究还发现，晶体结构的变化会影响TMDs的电子能带结构，进而影响磁性。在2H相中，能带结构相对简单，电子的巡游性较强；而在1T相中，由于结构的变化，能带结构变得更加复杂，出现了一些局域化的电子态，这些局域化电子态对磁性的形成和稳定起到了重要作用。4.2.2缺陷与掺杂对TMDs磁性的影响在过渡金属二硫化物（TMDs）中，缺陷和掺杂是调控其磁性的重要手段，它们通过改变材料的电子结构来影响磁性。硫空位是TMDs中常见的缺陷类型之一。当TMDs晶体中出现硫空位时，会导致局部电荷分布的改变。以MoS_2为例，硫空位的存在使得周围的钼原子的配位环境发生变化，原本与硫原子成键的电子云分布被打破。第一性原理计算表明，硫空位周围的钼原子会出现未成对电子，从而产生局域磁矩。这些局域磁矩之间通过超交换相互作用产生耦合，可能形成铁磁或反铁磁排列。研究发现，在含有硫空位的MoS_2体系中，当硫空位浓度较低时，局域磁矩之间的相互作用较弱，体系可能表现出顺磁性；随着硫空位浓度的增加，局域磁矩之间的耦合作用增强，可能会出现铁磁或反铁磁有序。这是因为随着硫空位浓度的增加，相邻硫空位之间的距离减小，超交换相互作用增强，使得磁矩之间的排列更加有序。金属空位同样会对TMDs的磁性产生显著影响。在MoS_2中，钼空位的出现会导致周围硫原子的电子云重新分布。由于钼原子的缺失，周围硫原子的电子失去了与之成键的对象，电子云发生畸变。第一性原理计算显示，钼空位周围的硫原子会出现自旋极化，形成局域磁矩。与硫空位不同的是，钼空位诱导的局域磁矩的大小和方向与周围硫原子的配位情况密切相关。在一些情况下，钼空位周围的硫原子形成的局域磁矩可能会相互耦合，形成稳定的磁有序结构。研究还发现，金属空位的存在会改变TMDs的能带结构。钼空位的出现会在能带中引入新的缺陷能级，这些缺陷能级与价带和导带相互作用，影响电子的占据情况，进而影响磁性。如果缺陷能级上的电子具有未成对自旋，就会对材料的磁性产生贡献。掺杂也是调控TMDs磁性的有效方法。过渡金属掺杂是常见的掺杂方式之一。当在MoS_2中掺入其他过渡金属原子（如Cr、Mn等）时，掺杂原子的电子结构与钼原子不同，会引入额外的电子或空穴。第一性原理计算表明，掺杂原子的未成对电子会与周围原子的电子发生相互作用，形成新的磁耦合。在Cr掺杂的MoS_2中，Cr原子的3d电子与周围硫原子的电子发生杂化，产生了强的磁相互作用。这种磁相互作用使得Cr掺杂的MoS_2体系表现出铁磁性。掺杂原子的浓度和分布对磁性也有重要影响。当掺杂浓度较低时，掺杂原子之间的相互作用较弱，体系的磁性主要由单个掺杂原子的磁矩决定；随着掺杂浓度的增加，掺杂原子之间的相互作用增强，可能会形成磁畴，导致体系的磁性增强。稀土元素掺杂同样能够改变TMDs的磁性。以MoS_2中掺入稀土元素Gd为例，Gd原子具有多个未成对的4f电子。第一性原理计算发现，Gd原子的4f电子与MoS_2中的电子发生弱的相互作用，这种相互作用虽然较弱，但由于4f电子的局域性和高自旋特性，使得Gd掺杂的MoS_2体系产生了一定的磁矩。稀土元素的掺杂还可能影响TMDs的晶体结构和电子能带结构。由于稀土元素的离子半径较大，掺入后可能会引起晶格畸变，从而改变原子间的距离和电子云的重叠程度，进一步影响磁性。Gd掺杂可能会导致MoS_2的晶格常数发生变化，这种变化会影响原子间的磁相互作用，进而改变材料的磁性。通过缺陷和掺杂来设计具有特定磁性的TMDs材料具有很大的潜力。可以根据实际需求，精确控制缺陷的类型、浓度和分布，以及掺杂原子的种类、浓度和位置。如果需要制备具有铁磁性的TMDs材料，可以通过引入适量的硫空位或过渡金属掺杂来实现；如果需要制备具有特殊磁各向异性的材料，可以通过控制缺陷或掺杂的分布来实现。通过第一性原理计算可以预先模拟不同缺陷和掺杂情况下TMDs的磁性变化，为实验制备提供理论指导，从而提高材料设计的效率和成功率。4.3其他二维材料的磁性研究4.3.1二维金属有机框架（MOFs）的磁性二维金属有机框架（MOFs）是一类由金属离子或金属簇与有机配体通过配位键自组装形成的二维材料。其结构特点在于具有周期性的多孔网络结构，这种结构赋予了二维MOFs独特的物理化学性质。在二维MOFs中，金属离子或金属簇作为节点，有机配体则作为连接臂，二者通过配位键相互连接，形成了具有规则排列的二维平面结构。在Co_3(HITP)_2中，钴离子（Co^{2+}）作为节点，与2,3,6,7,10,11-六氨基三苯（HITP）有机配体通过配位键连接，形成了具有Kagome晶格形式的二维结构。这种结构中，有机配体的长度和几何形状决定了二维MOFs的孔道大小和形状，而金属离子的种类和价态则影响着材料的电子结构和磁性。二维MOFs的磁性主要源于金属离子的未成对电子以及有机配体与金属离子之间的相互作用。金属离子的磁矩是二维MOFs磁性的重要来源。过渡金属离子通常具有未填满的d电子壳层，这些未成对电子使得金属离子具有固有磁矩。在二维MOFs中，金属离子的磁矩通过有机配体的桥连作用相互耦合，从而产生宏观的磁性。在基于铁离子（Fe^{3+}）的二维MOFs中，Fe^{3+}的3d电子中有5个未成对电子，具有较大的磁矩。这些磁矩通过有机配体的桥连作用，在一定条件下可以形成铁磁或反铁磁排列，使得二维MOFs表现出相应的磁性。有机配体与金属离子之间的相互作用对二维MOFs的磁性有着至关重要的影响。通过第一性原理计算可以深入研究这种相互作用。计算结果表明，有机配体的电子云与金属离子的电子云发生重叠，形成了配位键。这种配位键的形成不仅影响了金属离子的电子结构，还改变了金属离子磁矩之间的耦合方式。在一些二维MOFs中，有机配体的π电子与金属离子的d电子发生杂化，使得金属离子之间的磁相互作用增强，有利于形成稳定的磁有序结构。不同的有机配体具有不同的电子结构和空间构型，它们与金属离子的相互作用也会有所不同。含有氮、氧等杂原子的有机配体，由于杂原子的电负性较大，能够与金属离子形成更强的配位键，从而对磁性产生更大的影响。研究还发现，有机配体的长度和柔韧性也会影响金属离子之间的磁相互作用。较长或较柔韧的有机配体可能会导致金属离子之间的距离增大或磁相互作用的方向发生变化，从而改变二维MOFs的磁性。二维MOFs在磁性存储和催化等领域展现出了巨大的应用潜力。在磁性存储方面，二维MOFs的磁性可以通过外部磁场进行调控，且其具有较高的磁各向异性。这种特性使得二维MOFs有望成为下一代高密度磁性存储材料的候选者。通过控制二维MOFs的结构和组成，可以调节其磁各向异性，从而实现对磁矩取向的精确控制。在磁性存储器件中，利用二维MOFs的磁各向异性，可以提高存储单元的稳定性和读写速度。在催化领域，二维MOFs的大比表面积和丰富的活性位点使其具有良好的催化性能。磁性二维MOFs还可以通过外加磁场实现催化剂的分离和回收。在一些催化反应中，将磁性二维MOFs作为催化剂，反应结束后，可以通过施加外部磁场将催化剂从反应体系中分离出来，实现催化剂的重复利用，降低生产成本。4.3.2二维黑磷及其衍生物的磁性二维黑磷具有类似褶皱的层状晶体结构，由磷原子通过共价键相互连接形成。在黑磷的晶体结构中，磷原子呈锯齿状排列，形成了类似于蜂窝状的结构。这种结构赋予了二维黑磷独特的电子特性，使其具有直接带隙，且带隙随层数的变化可在0.3-2eV范围内调节。与石墨烯等二维材料不同，二维黑磷的电子结构具有明显的各向异性，这是由于其晶体结构在不同方向上的原子排列和键长存在差异。在平行于层平面的方向上，电子的运动受到一定的限制，导致电子的有效质量和迁移率在不同方向上有所不同。这种电子结构的各向异性对二维黑磷的电学和光学性质产生了重要影响。通过化学修饰和与其他材料复合等方式，可以诱导二维黑磷及其衍生物产生磁性。化学修饰是调控二维黑磷磁性的一种有效手段。氧化是常见的化学修饰方法之一。当二维黑磷被氧化时，表面的磷原子与氧原子发生化学反应，形成磷-氧键。第一性原理计算表明，氧化过程会导致黑磷表面的电子结构发生变化，产生未成对电子，从而形成局域磁矩。在轻度氧化的二维黑磷中，表面的磷原子被部分氧化，形成了磷氧化物基团。这些基团周围的电子云分布发生改变，使得电子自旋极化，产生了局域磁矩。研究还发现，氧化程度对二维黑磷的磁性有重要影响。随着氧化程度的增加，局域磁矩的数量和强度可能会发生变化。如果氧化程度过高，可能会导致黑磷的晶体结构被破坏，磁性反而减弱。卤化也是一种有效的化学修饰方法。当二维黑磷与卤素原子（如氟、氯、溴等）发生反应时，卤素原子会吸附在黑磷表面，与磷原子形成化学键。第一性原理计算显示，卤化过程会改变黑磷的电子结构，使得电子的自旋状态发生变化，从而产生磁性。在氟化的二维黑磷中，氟原子的电负性较大，与磷原子形成的化学键使得电子云向氟原子偏移，导致黑磷表面的电子自旋极化，形成局域磁矩。不同的卤素原子对二维黑磷磁性的影响可能会有所不同。由于卤素原子的电负性和原子半径不同，它们与黑磷的相互作用强度和方式也会存在差异，从而导致产生的磁性特征不同。与其他材料复合是诱导二维黑磷磁性的另一种重要方式。当二维黑磷与磁性材料复合时，通过界面处的电子相互作用，可以诱导黑磷产生磁性。二维黑磷与铁磁金属（如铁、钴、镍等）复合时，在界面处，黑磷的电子与金属原子的电子发生杂化。这种杂化作用使得黑磷的电子结构发生改变，电子自旋极化，从而产生磁矩。第一性原理计算表明，在二维黑磷/铁复合体系中，界面处的电子转移使得黑磷靠近铁的一侧出现自旋极化，形成了一定的磁矩。这种磁矩的大小和方向与铁原子的磁矩以及界面处的电子杂化强度密切相关。二维黑磷与磁性氧化物（如氧化铁、氧化钴等）复合时，同样可以通过界面处的电子相互作用诱导黑磷产生磁性。在二维黑磷/氧化铁复合体系中，氧化铁的磁性通过界面传递给黑磷，使得黑磷被诱导产生磁矩。这种复合体系的磁性还受到复合比例、界面结构等因素的影响。这些磁性黑磷材料在电子学和传感器等领域具有广阔的应用前景。在电子学领域，磁性黑磷材料的独特磁性和电学性质使其有望用于制造自旋电子学器件。在自旋晶体管中，利用磁性黑磷材料的自旋极化特性，可以实现对电子自旋的有效调控，从而提高晶体管的性能。磁性黑磷材料还可以用于制造磁性存储器，通过控制磁矩的取向来存储信息，有望提高存储密度和读写速度。在传感器领域，磁性黑磷材料对某些气体分子具有特殊的吸附和电学响应特性。当磁性黑磷材料吸附特定气体分子时，其磁性和电学性质会发生变化。利用这种特性，可以制备高灵敏度的气体传感器。在检测氨气时，磁性黑磷材料吸附氨气分子后，磁矩和电阻会发生明显变化，通过检测这些变化就可以实现对氨气的快速、准确检测。五、二维材料磁性的调控与应用前景5.1二维材料磁性的调控方法5.1.1电场调控电场调控二维材料磁性的原理基于电场诱导的电荷转移和能带结构变化。当在二维材料上施加电场时，会导致材料内部的电荷分布发生改变。在二维过渡金属硫族化合物（TMDs）中，如MoS_2，施加电场会使电子在材料内部重新分布。第一性原理计算表明，电场会改变MoS_2中过渡金属原子（如Mo）与硫原子之间的电子云分布，导致原子磁矩的变化。这是因为电场的作用使得电子的能量状态发生改变，电子在不同原子轨道之间的分布发生调整，从而影响了原子的磁矩。电场还会引起二维材料能带结构的变化。对于具有半导体特性的二维材料，如MoS_2，电场可以改变其能带的宽度和位置，进而影响材料的磁性。当施加电场时，MoS_2的导带和价带会发生移动，带隙大小也会改变。这种能带结构的变化会影响电子的占据情况，从而改变材料的磁性。如果电场导致带隙减小，电子更容易激发到导带，可能会改变材料的磁矩。在相关的第一性原理计算和实验研究中，取得了许多有价值的成果。理论计算预测，在双层CrI_3中，通过施加电场可以实现层间磁耦合的调控。当施加正电场时，双层CrI_3的层间磁耦合由反铁磁转变为铁磁；当施加负电场时，层间磁耦合保持反铁磁状态。这一预测得到了实验的验证。实验上，通过在双层CrI_3上施加电场，利用磁光克尔效应等技术测量其磁性变化，结果与理论计算相符。在二维材料与衬底组成的异质结构中，电场调控也展现出重要作用。以石墨烯与SiO_2衬底组成的异质结构为例，通过在衬底上施加电场，可以调控石墨烯的磁性。实验发现，随着电场强度的增加，石墨烯的磁矩逐渐增大，这是由于电场诱导的电荷转移改变了石墨烯的电子结构，从而产生了磁性。电场调控在实现二维材料磁性动态调控方面具有显著优势。电场调控具有快速响应的特点。相比于其他调控方法，如温度调控，电场调控可以在短时间内实现对磁性的改变。这使得电场调控在高速自旋电子学器件中具有潜在应用价值。在自旋逻辑器件中，利用电场对二维材料磁性的快速调控，可以实现信息的快速处理和传输。电场调控具有可精确控制的优点。通过调节电场的强度和方向，可以精确地控制二维材料的磁性状态。这为制备高性能的磁性器件提供了可能。在磁存储器件中，利用电场的精确调控，可以实现对存储单元磁性状态的精确写入和读取，提高存储密度和读写速度。然而，电场调控也面临一些挑战。电场调控的效果受到材料与衬底之间界面性质的影响。如果界面存在缺陷或杂质，会影响电荷的转移和电场的分布，从而降低电场调控的效果。在制备二维材料与衬底的异质结构时，需要精确控制界面的质量，以提高电场调控的效率。电场调控还受到材料本身性质的限制。对于一些二维材料，由于其电子结构的稳定性，电场对其磁性的调控作用可能较弱。在一些具有强共价键的