基于第一性原理探究二维MoS₂与MnO磁学性质的微观机制与应用潜力

基于第一性原理探究二维MoS₂与MnO磁学性质的微观机制与应用潜力一、引言1.1研究背景与意义近年来，二维材料由于其独特的原子结构和电子特性，在凝聚态物理和材料科学领域引起了广泛关注。自2004年石墨烯被首次成功剥离以来，二维材料的研究取得了飞速发展，众多具有优异电学、光学、力学和热学性能的二维材料被相继发现和研究。二维材料的原子级厚度赋予其与体相材料截然不同的物理性质，例如，石墨烯具有超高的载流子迁移率和优异的力学性能，使其在高速电子器件和柔性电子学领域展现出巨大的应用潜力；过渡金属二硫属化合物（TMDs），如MoS₂、WS₂等，具有直接带隙，在光电器件，如光电探测器、发光二极管等方面具有广阔的应用前景。磁性材料作为一类重要的功能材料，在现代科技中发挥着不可或缺的作用。从日常生活中的电子设备到高端的信息技术、能源存储与转换、生物医学等领域，磁性材料都有着广泛的应用。例如，在数据存储领域，硬盘、磁带等存储介质依赖磁性材料来记录和读取信息；在电力传输中，变压器中的铁芯利用磁性材料的高磁导率来实现高效的电磁能量转换；在生物医学领域，磁性纳米颗粒被用于磁共振成像（MRI）对比增强、药物靶向递送和磁热疗等。随着科技的不断进步，对磁性材料性能的要求也日益提高，探索具有新颖磁学性质和高性能的磁性材料成为材料科学领域的研究热点之一。二维MoS₂作为一种典型的过渡金属二硫属化合物，具有独特的晶体结构和电学性质。其晶体结构由一层钼原子（Mo）夹在两层硫原子（S）之间形成三明治结构，层间通过范德华力相互作用。这种结构赋予MoS₂一些特殊的性质，如在单层状态下具有直接带隙（约1.8eV），而体相MoS₂则为间接带隙。由于其原子级厚度和特殊的电子结构，二维MoS₂在电学、光学和力学等方面表现出与传统材料不同的特性，使其在纳米电子学、光电器件、传感器和催化等领域具有潜在的应用价值。然而，本征的二维MoS₂通常不具有磁性，通过对其进行磁性调控，引入磁性，有望拓展其在自旋电子学等领域的应用，为实现新型的自旋相关器件提供可能。MnO是一种典型的过渡金属氧化物磁性材料，具有反铁磁性质。MnO的晶体结构为岩盐结构，其中Mn²⁺离子和O²⁻离子交替排列。在低温下，MnO表现出反铁磁有序，其奈尔温度（TN）约为118K。MnO的磁学性质使其在磁存储、磁传感器和自旋电子学等领域具有潜在的应用价值。对MnO的磁学性质进行深入研究，不仅有助于理解其基本的磁相互作用机制，而且对于开发基于MnO的高性能磁性材料和器件具有重要的指导意义。对二维MoS₂和MnO磁学性质的研究具有重要的理论和实际意义。从理论角度来看，二维MoS₂和MnO的磁学性质研究涉及到电子结构、磁相互作用、晶体场效应等多个重要的物理概念和理论。通过研究二维MoS₂的磁性调控机制，可以深入理解电子关联、自旋轨道耦合等量子效应在低维体系中的作用，丰富和完善低维磁性物理理论。对于MnO，研究其磁学性质与晶体结构、电子结构之间的关系，有助于揭示过渡金属氧化物磁性材料中复杂的磁相互作用机制，为磁性材料的理论研究提供重要的实验和理论依据。在实际应用方面，二维MoS₂和MnO的磁学性质研究成果具有广泛的应用前景。对于二维MoS₂，引入磁性后，有望开发出基于二维MoS₂的新型自旋电子学器件，如自旋场效应晶体管、磁随机存取存储器等，这些器件具有低功耗、高速读写和高集成度等优点，将为下一代信息技术的发展提供关键支撑。在传感器领域，利用二维MoS₂的磁性和电学性质的耦合效应，可以制备高灵敏度的磁传感器，用于生物分子检测、磁场探测等。对于MnO，深入了解其磁学性质有助于优化其在磁存储和磁传感器等领域的性能，提高存储密度和传感灵敏度。此外，MnO在生物医学领域的潜在应用，如磁热疗和磁共振成像对比增强等，也需要对其磁学性质有深入的认识和精确的调控。因此，开展二维MoS₂和MnO磁学性质的研究，对于推动材料科学、凝聚态物理和相关应用技术的发展具有重要的意义。1.2国内外研究现状在二维MoS₂磁学性质的研究方面，国内外研究取得了诸多进展。理论研究上，通过第一性原理计算，许多学者对其本征磁性以及掺杂、缺陷等对磁性的影响展开探讨。例如，一些研究表明，本征二维MoS₂由于其原子结构和电子分布特点，通常呈现非磁性，但引入特定杂质原子进行掺杂时，可能诱导出磁性。国内的一些研究团队通过第一性原理计算发现，过渡金属原子如Mn、Fe等掺杂到MoS₂中，会产生局域磁矩，且磁矩大小和方向与掺杂原子的种类、位置以及周围原子环境密切相关。在实验研究中，科研人员尝试多种方法制备具有磁性的二维MoS₂材料。化学气相沉积（CVD）法是常用的制备手段之一，通过精确控制生长条件，可实现对MoS₂晶体结构和缺陷的调控，进而影响其磁学性质。国外有研究利用CVD法在特定衬底上生长MoS₂薄膜，并通过后续的离子注入等工艺引入磁性离子，成功观测到MoS₂薄膜的磁响应。在MnO磁学性质的研究领域，国内外学者也进行了大量工作。理论研究深入分析了MnO的晶体结构与磁相互作用之间的关系。MnO的岩盐结构中，Mn²⁺离子的3d电子对其磁学性质起关键作用，通过晶体场理论和自旋相关理论，可解释其反铁磁有序的形成机制。国内研究团队利用先进的理论计算方法，对MnO在不同压力、温度条件下的磁学性质变化进行模拟，揭示了外部条件对其磁相互作用的影响规律。在实验方面，研究手段不断丰富。中子散射技术是研究MnO磁结构的重要工具，通过中子与MnO样品中原子核和磁矩的相互作用，可精确测定其磁有序结构和磁激发谱。国外研究人员利用中子散射实验，详细研究了MnO在低温下的反铁磁相变过程，确定了其奈尔温度以及磁有序转变的微观机制。此外，磁性测量技术如超导量子干涉仪（SQUID）磁强计等，用于精确测量MnO的磁化强度、磁滞回线等宏观磁学性质，为深入理解其磁行为提供了重要数据。尽管国内外在二维MoS₂和MnO磁学性质研究方面取得了显著成果，但仍存在一些不足和待解决问题。对于二维MoS₂，目前在磁性调控的精确性和稳定性方面仍有待提高。虽然通过掺杂等手段可诱导磁性，但如何实现对磁矩大小、方向以及磁性均匀性的精确控制，以满足实际应用中对磁性材料性能一致性的要求，仍是研究的难点。在制备工艺方面，现有的制备方法大多存在制备过程复杂、成本高、产量低等问题，限制了具有磁性的二维MoS₂材料的大规模应用。此外，对于二维MoS₂与衬底或其他材料复合时的界面磁学性质研究还相对较少，而界面磁学性质对其在自旋电子学器件中的应用性能有着重要影响。对于MnO，虽然对其基本磁学性质有了较为深入的了解，但在高温、高压等极端条件下的磁学性质以及与其他材料复合后的磁学性能变化研究还不够充分。随着应用需求的不断提高，需要进一步探索如何通过材料设计和制备工艺优化，提高MnO的磁性能，如提高其磁导率、降低磁损耗等。在MnO基复合材料的研究中，如何实现不同组分之间的协同作用，以获得具有优异综合性能的磁性材料，也是亟待解决的问题。此外，对于MnO在生物医学等新兴应用领域的磁学性质与生物兼容性之间的关系研究还处于起步阶段，需要进一步深入探索。1.3研究内容与方法本研究聚焦于二维MoS₂和MnO的磁学性质，运用基于密度泛函理论的第一性原理计算方法，深入探究其微观结构与磁学性能之间的内在联系，具体研究内容如下：二维MoS₂和MnO的晶体结构优化与稳定性分析：构建二维MoS₂和MnO的初始晶体结构模型，借助第一性原理计算对其进行结构优化，获取体系能量最低的稳定结构。通过分析优化后的晶格参数、原子坐标以及键长、键角等结构信息，明确二维MoS₂和MnO的晶体结构特征。进一步计算结构的形成能、结合能等物理量，评估其结构稳定性，深入探讨结构稳定性与原子间相互作用、电子云分布之间的关联。二维MoS₂和MnO的电子结构计算与分析：基于优化后的稳定结构，计算二维MoS₂和MnO的电子结构，包括能带结构、态密度（总态密度和分波态密度）以及电荷密度分布等。通过对能带结构的分析，确定材料的能带宽度、带隙大小和类型（直接带隙或间接带隙），深入理解电子在材料中的能量分布和运动状态。结合态密度和电荷密度分布，研究不同原子轨道对电子态的贡献，揭示电子在原子间的分布和转移情况，从电子层面阐释材料的电学和磁学性质起源。二维MoS₂和MnO的磁学性质计算与分析：运用第一性原理计算方法，详细研究二维MoS₂和MnO的磁学性质，包括磁矩大小、磁矩方向、磁相互作用类型（铁磁、反铁磁或亚铁磁）以及磁各向异性等。针对二维MoS₂，重点研究通过掺杂、缺陷等手段引入磁性的机制，分析掺杂原子种类、浓度、位置以及缺陷类型对磁学性质的影响规律。对于MnO，深入探讨其反铁磁性质与晶体结构、电子结构之间的内在关系，研究外部条件（如压力、温度）对其磁学性质的调控作用。本研究采用基于密度泛函理论（DFT）的第一性原理计算方法，该方法以量子力学为基础，从电子和原子核的相互作用出发，通过求解多电子体系的薛定谔方程，精确计算材料的电子结构和各种物理性质。在计算过程中，无需借助任何经验参数，仅需输入材料的原子种类和坐标等基本信息，便可从理论上预测材料的结构和性质。本研究选用VASP（ViennaAbinitioSimulationPackage）软件进行计算，VASP是一款广泛应用于材料科学领域的第一性原理计算软件，具有计算精度高、计算速度快、功能强大等优点。它采用平面波赝势方法，能够有效处理周期性边界条件下的晶体结构，在计算材料的电子结构、磁学性质、光学性质等方面表现出色。在计算中，选用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）交换关联泛函，以准确描述电子之间的交换关联相互作用。同时，合理设置平面波截断能和k点网格，确保计算结果的准确性和收敛性。二、理论基础与计算方法2.1第一性原理基本理论2.1.1密度泛函理论（DFT）密度泛函理论（DensityFunctionalTheory，DFT）是一种研究多电子体系电子结构的量子力学方法，在材料科学、化学等领域有着广泛的应用。其核心思想是将多电子体系的基态能量表示为电子密度的泛函，从而将复杂的多体问题简化为相对简单的电子密度函数问题。在传统的量子力学方法中，如Hartree-Fock方法，是基于复杂的多电子波函数来描述体系的状态。然而，多电子波函数依赖于3N个变量（N为电子数，每个电子包含三个空间变量），这使得计算量随着电子数的增加呈指数级增长，在实际应用中面临巨大的挑战。DFT的主要目标是用电子密度取代波函数作为研究的基本量，电子密度仅是三个空间变量的函数，在概念和计算上都更为方便处理。Hohenberg-Kohn定理是DFT的基石。Hohenberg-Kohn第一定理指出，对于任何外部势场作用下的多电子体系，其基态能量是该体系电子密度的唯一泛函。这意味着，只要确定了电子密度分布，就可以确定体系的所有基态物理量，如基态能量、电子结构等，从而避免了直接计算复杂的多体波函数。Hohenberg-Kohn第二定理进一步阐述了如何通过变分原理确定体系的基态电子密度。该定理表明，对于任意一个体系的基态能量泛函，其对应的真实电子密度是使该泛函取最小值的电子密度，即在所有可能的电子密度分布中，真实基态电子密度使得体系的能量达到最低。通过这个变分原理，我们可以将原本需要解决的多体波函数问题转化为寻找电子密度泛函极小值的问题，极大地简化了计算过程。在实际应用中，Kohn-Sham方程是实现DFT计算的关键。Kohn和Sham在1965年提出了一种近似的处理方式，引入一组虚拟的非相互作用电子来代替实际的相互作用电子，并构建一个与真实体系等效的非相互作用体系。Kohn-Sham方程的表达式为：\left[-\frac{\hbar^{2}}{2m}\nabla^{2}+V_{ext}(\vec{r})+V_{H}(\vec{r})+V_{xc}(\vec{r})\right]\psi_{i}(\vec{r})=\epsilon_{i}\psi_{i}(\vec{r})其中，-\frac{\hbar^{2}}{2m}\nabla^{2}是动能算符，V_{ext}(\vec{r})是外部势场，V_{H}(\vec{r})是Hartree势，描述电子之间的库仑相互作用，V_{xc}(\vec{r})是交换-关联势，用于描述电子之间的交换和关联效应，\psi_{i}(\vec{r})是Kohn-Sham轨道，\epsilon_{i}是对应的本征能量。通过求解Kohn-Sham方程，可以得到Kohn-Sham轨道和本征能量，进而计算出体系的电子密度n(\vec{r})=\sum_{i=1}^{N_{occ}}|\psi_{i}(\vec{r})|^{2}，其中N_{occ}是占据轨道的数目。体系的总能量可以表示为：E=\sum_{i=1}^{N_{occ}}\epsilon_{i}-\intd\vec{r}n(\vec{r})V_{H}(\vec{r})-\intd\vec{r}n(\vec{r})V_{xc}(\vec{r})+E_{xc}[n]其中，E_{xc}[n]是交换-关联能泛函。处理交换-关联能E_{xc}是KSDFT中的难点，目前并没有精确求解的方法。常见的近似方法有局域密度近似（LDA）和广义梯度近似（GGA）。LDA近似使用均匀电子气来计算体系的交换能（均匀电子气的交换能是可以精确求解的），而相关能部分则采用对自由电子气进行拟合的方法来处理。GGA则在LDA的基础上，考虑了电子密度的梯度信息，能更准确地描述电子间的交换和关联作用。在本研究中，选用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）交换关联泛函，以准确描述电子之间的交换关联相互作用。2.1.2赝势平面波方法在基于密度泛函理论的第一性原理计算中，赝势平面波方法是一种常用且有效的计算方法。实际求解Kohn-Sham方程时，由于原子核产生的势场项在原子中心是发散的，波函数在靠近原子核的区域变化剧烈，需要采用大量的平面波展开才能准确描述电子的行为，这使得计算成本变得非常大。为了降低计算量，赝势平面波方法引入了赝势的概念，并采用平面波基组来展开晶体波函数。平面波是自由电子气的本征函数，具有标准正交化和能量单一性的特点，对任何原子都适用且等同对待空间中的任何区域，不需要修正重叠误差。根据晶体的空间平移对称性，布洛赫（Bloch）定理证明，能带电子的波函数\psi(\vec{r},\vec{k})总是能够写成\psi(\vec{r},\vec{k})=u(\vec{r},\vec{k})e^{i\vec{k}\cdot\vec{r}}，其中\vec{k}是电子波矢，u(\vec{r},\vec{k})是具有晶体平移周期性的周期函数。采用周期性边界条件后，单粒子轨道波函数可以用平面波基展开为\psi(\vec{r})=\frac{1}{\sqrt{\Omega}}\sum_{\vec{G}}\varphi(\vec{G})e^{i(\vec{K}+\vec{G})\cdot\vec{r}}，式中\frac{1}{\sqrt{\Omega}}是归一化因子，\Omega是原胞体积；\vec{G}是原胞的倒格矢，\vec{K}是第一Brillouin区的波矢，\varphi(\vec{G})是展开系数。然而，由于电子数目的无限性，从原则上讲\vec{K}矢量的个数是无限的，每个\vec{K}矢量处的电子波函数展开形式包含的平面波数量也是无限多的。基于计算成本的考虑，实际计算中只能取有限个平面波数。通常的做法是确定一个截断能量E_{cut}，给定截断能E_{cut}=\frac{\hbar^{2}(\vec{G}+\vec{K})^{2}}{2m}，对\vec{G}的求和可以限制在(\vec{G}+\vec{K})^{2}/2\leqE_{cut}的范围内，即要求用于展开的波函数的能量小于E_{cut}。通过增加截断能量可以减小因函数基组不完备而产生的总能量计算误差。赝势的引入是赝势平面波方法的关键。由于全电子的计算工作量大，非全电子的计算通常具有优势。当不考虑原子的芯电子时，原子成为“赝原子”（原子芯+价电子），这时价电子运动所处的势场就相当于来自原子芯的“赝势”（原本是原子芯内所有电子提供的势场）。可以证明，将薛定谔方程中的势能换成赝势，则存在相应的赝波函数，使体系的本征值不变。赝势的作用是在保持计算精度的前提下，将原子核附近快速振荡的真实波函数替换为相对平滑的赝波函数，从而减少平面波基组的数量，降低计算成本。同时，赝势能够保证在某个半径r_{c}之外，赝原子的能量本征值以及赝波函数都和“整个原子”时的解一致，在r_{c}半径之内，赝势应当尽量平滑，使波函数的振荡很小，且赝波函数没有节点。赝势平面波方法具有较高的计算效率，使之日益发展成为有效的计算方法。该方法结合了平面波基组的优点和赝势的近似处理，既能够利用平面波基组方便地采用快速傅里叶变换（FFT）技术，使能量、力等的计算在实空间和倒空间快速转换，减少计算工作量，又能够通过赝势的引入降低计算成本，适用于计算较大体系的基态属性。在本研究中，采用赝势平面波方法，能够有效地对二维MoS₂和MnO体系进行计算，准确地获取其电子结构和磁学性质。2.2计算软件与模型构建2.2.1VASP软件简介VASP（ViennaAbinitioSimulationPackage）是一款在材料计算领域应用广泛且功能强大的软件，由维也纳大学Hafner小组开发。它基于量子力学原理，能够精确地模拟材料的电子结构和各种物理性质。VASP的核心优势在于其采用密度泛函理论（DFT）来描述电子之间的相互作用，通过求解Kohn-Sham方程，可获得材料中电子的分布和能量状态，进而计算出材料的多种性质。在功能方面，VASP具备处理多种体系的能力，无论是原子、分子、团簇，还是纳米线、薄膜、晶体等材料体系，它都能通过周期性边界条件（或超原胞模型）进行有效的模拟。在计算材料的结构参数时，VASP能够精确地优化材料的晶体结构，确定键长、键角、晶格常数以及原子的精确位置等信息。通过计算材料的状态方程和力学性质，VASP可以得到体弹性模量、弹性常数等关键力学参数，为材料的力学性能评估提供重要依据。在电子结构计算方面，VASP能够给出材料的能级分布、电荷密度分布、能带结构以及电子态密度等信息，这些信息对于理解材料的电学、光学和磁学性质起着至关重要的作用。在研究材料的光学性质时，VASP可计算材料的介电函数、吸收光谱等光学参数，为光电器件的设计和优化提供理论支持。对于磁学性质的计算，VASP能够确定材料的磁矩大小、磁矩方向、磁相互作用类型以及磁各向异性等，为磁性材料的研究和应用提供关键数据。此外，VASP还支持从头分子动力学模拟，能够研究材料在原子尺度下的动态行为，如原子的扩散、化学反应过程等。在材料计算领域，VASP软件的应用极为广泛。在凝聚态物理研究中，科研人员利用VASP研究各种新型材料的电子结构和物理性质，探索材料的新奇量子现象，如拓扑绝缘体的拓扑性质、高温超导材料的超导机制等。在材料科学研究中，VASP被用于材料的设计和优化，通过模拟不同成分和结构的材料性能，筛选出具有优异性能的材料，为新型材料的研发提供理论指导。例如，在开发新型超导材料时，研究人员使用VASP计算不同元素掺杂对材料电子结构和超导性能的影响，从而指导实验合成具有更高超导转变温度的材料。在化学领域，VASP可用于研究化学反应过程和催化剂的作用机制，通过模拟分子在催化剂表面的吸附和反应过程，深入了解催化反应的微观机理，为新型催化剂的设计提供依据。总之，VASP软件凭借其强大的功能和高精度的计算结果，成为材料计算领域不可或缺的工具，为众多科研和应用领域提供了重要的理论支持和技术手段。2.2.2二维MoS₂模型构建二维MoS₂具有独特的原子结构，其晶体结构属于六方晶系。在构建二维MoS₂模型时，首先考虑其基本的原子排列方式。二维MoS₂的基本结构单元由一层钼原子（Mo）夹在两层硫原子（S）之间组成，形成类似于三明治的结构。Mo原子位于中间层，与上下两层的S原子通过共价键相互连接。在这种结构中，Mo-S键长以及S-Mo-S键角等结构参数对于描述其原子结构和物理性质至关重要。通过查阅相关文献以及实验数据，确定构建二维MoS₂模型的晶格参数。二维MoS₂的晶格常数a和b相等，通常取值约为3.16Å，这一数值与实验测量结果相吻合。晶格常数c表示层间距离，对于二维MoS₂，由于层间主要通过范德华力相互作用，c值相对较大，一般取值约为6.15Å。在建立模型时，选用合适的超胞大小以满足计算需求和精度要求。通常采用2×2×1的超胞，这样的超胞尺寸既能保证计算的准确性，又能在合理的计算资源范围内完成计算任务。在该超胞中，包含了一定数量的Mo和S原子，具体原子数目根据超胞的构建方式和原子排列确定。在这个2×2×1的超胞中，Mo原子的数量为4个，S原子的数量为8个。在模型中，明确各个原子的坐标位置，以准确描述原子在超胞中的相对位置关系。根据六方晶系的对称性和晶格参数，确定Mo原子和S原子在超胞中的坐标，使得模型能够准确反映二维MoS₂的原子结构特征。为了避免二维MoS₂模型在计算过程中与相邻层发生不必要的相互作用，在垂直于二维平面的方向上添加足够厚的真空层。一般来说，添加15Å左右的真空层能够有效地隔离相邻层之间的相互影响，确保计算结果仅反映二维MoS₂本身的性质。通过以上步骤，成功构建了用于第一性原理计算的二维MoS₂模型，该模型准确地描述了二维MoS₂的原子结构和晶格参数，为后续研究其电子结构和磁学性质奠定了坚实的基础。2.2.3MnO模型构建MnO的晶体结构为典型的岩盐结构（NaCl型结构），属于立方晶系。在岩盐结构中，Mn²⁺离子和O²⁻离子交替排列，形成面心立方晶格。每个Mn²⁺离子周围被六个O²⁻离子以八面体配位的方式包围，同样，每个O²⁻离子也被六个Mn²⁺离子以八面体配位的方式包围。这种离子排列方式决定了MnO的晶体结构特征和物理性质。确定MnO模型的晶格参数是构建模型的关键步骤之一。根据大量的实验研究和理论计算结果，MnO的晶格常数a通常取值约为4.44Å。在构建模型时，采用1×1×1的原胞作为基本单元，该原胞中包含4个Mn²⁺离子和4个O²⁻离子。在原胞中，Mn²⁺离子和O²⁻离子的坐标位置根据立方晶系的对称性确定。例如，在面心立方晶格中，Mn²⁺离子可以位于(0,0,0)、(0.5,0.5,0)、(0.5,0,0.5)和(0,0.5,0.5)等位置，而O²⁻离子则位于(0.25,0.25,0.25)、(0.75,0.75,0.25)、(0.75,0.25,0.75)和(0.25,0.75,0.75)等位置，这样的坐标设置符合MnO的岩盐结构特征。在进行第一性原理计算时，为了满足周期性边界条件的要求，确保计算结果的准确性和可靠性，对构建好的MnO原胞模型进行适当的扩展。可以根据实际计算需求，将原胞扩展为超胞，如2×2×2的超胞。在扩展后的超胞中，原子的坐标位置相应地进行调整，以保持晶体结构的对称性和完整性。通过以上方法，成功构建了MnO的晶体结构模型，该模型准确地反映了MnO的晶体结构特征和原子坐标信息，为后续深入研究MnO的电子结构和磁学性质提供了可靠的模型基础。三、二维MoS₂磁学性质的第一性原理研究3.1本征二维MoS₂的结构与电子性质3.1.1晶体结构分析本征二维MoS₂具有独特的晶体结构，属于六方晶系。其基本结构单元呈现出典型的“三明治”构型，由一层钼原子（Mo）夹在两层硫原子（S）之间构成。在这种结构中，Mo原子与上下两层的S原子通过共价键紧密相连，形成稳定的原子排列。从晶体对称性角度来看，二维MoS₂具有高度的对称性。其晶格点阵为二维的简单六角结构，晶格常数a和b相等，通常取值约为3.16Å，这一数值与众多实验测量结果以及理论计算值高度吻合。晶格常数c表示层间距离，由于层间主要通过范德华力相互作用，c值相对较大，一般约为6.15Å。这种层状结构使得二维MoS₂在平面内具有较好的稳定性和各向同性，而在垂直于平面方向上则表现出较弱的相互作用。在原子间的键长和键角方面，Mo-S键长约为0.2411nm，这一长度反映了Mo和S原子之间的共价键强度。S-Mo-S键角最小约为46.21度，这种特定的键角决定了MoS₂的原子空间排列方式，对其电子结构和物理性质产生重要影响。在原胞结构中，虽然晶体为二维结构，但原胞并非二维。以常见的2×2×1超胞为例，其中包含4个Mo原子和8个S原子。Mo原子位于平行四边形的四个角（较小内角为60度），原胞内部的两个S原子分别处于三个Mo原子（正三角形）的正上方和正下方。这种原子排列方式形成了稳定的晶体结构，为二维MoS₂的各种物理性质奠定了基础。3.1.2电子结构计算通过基于密度泛函理论的第一性原理计算，得到了本征二维MoS₂的能带结构和态密度图。能带结构是描述材料中电子能量与波矢关系的重要物理量，它反映了电子在晶体中的运动状态和能量分布。在二维MoS₂的能带结构中，价带顶和导带底均位于布里渊区的K点，这表明二维MoS₂是直接带隙半导体。计算得到的带隙宽度约为1.8eV，这一数值与实验测量结果以及其他理论计算结果相符。直接带隙半导体具有较高的光吸收和发射效率，使得二维MoS₂在光电器件应用中具有独特的优势。从态密度图可以更直观地了解电子在不同能量状态下的分布情况。总态密度反映了材料中所有电子态在能量空间的分布，而分波态密度则能够进一步分析不同原子轨道对电子态的贡献。在二维MoS₂中，价带主要由S原子的3p轨道和Mo原子的4d轨道电子贡献，导带则主要由Mo原子的4d轨道电子贡献。在费米能级附近，S原子的3p轨道和Mo原子的4d轨道电子态密度存在明显的重叠，这表明在该能量区域，电子在S和Mo原子之间存在较强的相互作用。这种电子相互作用对二维MoS₂的电学和光学性质有着重要影响。在光吸收过程中，电子可以从价带跃迁到导带，由于价带和导带在K点的直接带隙特性以及电子态密度的分布情况，使得二维MoS₂对特定波长的光具有较强的吸收能力。3.1.3磁学性质分析理论计算和实验研究均表明，本征二维MoS₂通常不具有磁性。这一磁性特征与其晶体结构和电子分布密切相关。在二维MoS₂的晶体结构中，原子的电子云分布呈现出高度的对称性。Mo原子和S原子通过共价键结合，电子在原子之间的分布相对均匀，没有明显的未成对电子。根据磁性的产生原理，磁性通常源于材料中原子的未成对电子产生的磁矩。在本征二维MoS₂中，由于电子的配对情况良好，不存在未成对电子，因此无法形成有效的磁矩，导致材料整体表现为非磁性。从电子结构角度进一步分析，二维MoS₂的能带结构和态密度分布也支持其无磁性的结论。在能带结构中，价带和导带的电子填充情况使得电子的总自旋为零，没有净磁矩产生。态密度图中，不同原子轨道的电子态分布没有出现导致磁性的特殊情况。与具有磁性的材料相比，例如过渡金属铁（Fe），其原子具有未成对的3d电子，这些未成对电子的自旋产生磁矩，使得Fe具有铁磁性。而二维MoS₂中不存在类似的未成对电子结构，因此不具备磁性。3.2掺杂对二维MoS₂磁学性质的影响3.2.1常见掺杂元素选择为了调控二维MoS₂的磁学性质，选择合适的掺杂元素至关重要。常见的掺杂元素包括过渡金属元素（如Mn、Fe、Co、Ni等）和非金属元素（如N、P等）。这些元素的选择基于多方面的考虑。过渡金属元素具有未填满的d电子壳层，其d电子的自旋和轨道角动量对材料的磁性产生重要影响。例如，Mn原子具有5个未成对的3d电子，其磁矩较大，当Mn原子掺杂到二维MoS₂中时，有可能通过与MoS₂中的电子相互作用，诱导出磁性。Fe原子同样具有多个未成对的3d电子，其掺杂有望在MoS₂中引入局域磁矩，并可能通过磁相互作用形成长程磁有序。非金属元素如N、P等的掺杂也备受关注。N原子的电负性与S原子相近，但价电子数不同，掺杂后可能改变MoS₂的电子结构和电荷分布，从而影响其磁性。N原子掺杂可能在MoS₂中引入额外的电子或空穴，这些载流子与MoS₂的本征电子相互作用，产生自旋极化，进而诱导出磁性。P原子的外层电子结构与S原子不同，P掺杂可能导致MoS₂晶格中的电子云分布发生变化，影响电子的自旋状态，为磁性的调控提供新的途径。从理论和实验研究来看，不同元素掺杂对MoS₂磁性的影响呈现出多样化的特征。过渡金属掺杂往往能引入较大的局域磁矩，且磁矩的大小和方向与掺杂原子的种类、浓度以及在MoS₂晶格中的位置密切相关。一些研究表明，随着过渡金属掺杂浓度的增加，MoS₂的磁矩可能逐渐增大，材料可能从非磁性转变为磁性。然而，过高的掺杂浓度可能导致晶格畸变加剧，影响材料的稳定性和磁性的均匀性。非金属元素掺杂对MoS₂磁性的影响相对较为复杂，不仅涉及电子结构的改变，还可能影响MoS₂的化学活性和晶体结构的稳定性。在某些情况下，非金属掺杂可能通过与MoS₂形成特定的化学键或电子态，间接影响其磁性。3.2.2掺杂体系的结构优化对掺杂后的二维MoS₂体系进行结构优化是研究其磁学性质的重要基础。以Mn掺杂为例，在构建的2×2×1的MoS₂超胞中，将一个Mo原子替换为Mn原子。通过基于密度泛函理论的第一性原理计算，利用VASP软件对掺杂体系进行结构优化，使体系达到能量最低的稳定状态。优化后的结构参数表明，掺杂原子对晶格结构产生了显著影响。与本征二维MoS₂相比，Mn掺杂后，Mo-S键长发生了变化。原本Mo-S键长约为0.2411nm，掺杂后，与Mn原子相邻的S原子的Mo-S键长变为约0.2450nm。这种键长的变化是由于Mn原子的原子半径与Mo原子不同，以及Mn原子与周围S原子之间的电子相互作用与Mo-S之间的相互作用存在差异。Mn原子的原子半径略大于Mo原子，当Mn原子替代Mo原子后，为了适应新的原子尺寸和电子云分布，周围的S原子会发生位置调整，导致Mo-S键长改变。键角也发生了明显的变化。S-Mo-S键角在本征MoS₂中最小约为46.21度，Mn掺杂后，与Mn原子相关的S-Mn-S键角变为约47.50度。这种键角的改变进一步影响了MoS₂的晶体结构对称性。原本二维MoS₂具有高度的六角对称性，Mn掺杂后，由于键长和键角的变化，晶格的局部对称性被打破，导致晶体结构的整体对称性降低。这种结构变化对MoS₂的电子结构和磁学性质产生了深远影响。结构的变化会导致电子云分布的改变，进而影响电子的能量状态和自旋状态，最终影响材料的磁性。3.2.3电子结构与磁学性质变化通过第一性原理计算，对比掺杂前后二维MoS₂的电子结构和磁学性质，发现掺杂导致了显著的变化。从能带结构来看，本征二维MoS₂是直接带隙半导体，带隙宽度约为1.8eV。当Mn掺杂后，能带结构发生了明显的变化。在费米能级附近出现了新的杂质能级，这些杂质能级是由Mn原子的3d电子与MoS₂中的电子相互作用产生的。杂质能级的出现改变了电子的跃迁路径和能量状态，使得材料的电学和磁学性质发生改变。由于杂质能级的存在，电子在跃迁过程中可能会发生自旋翻转，从而产生磁矩。从态密度图可以更清晰地看到电子结构的变化。本征MoS₂的价带主要由S原子的3p轨道和Mo原子的4d轨道电子贡献，导带主要由Mo原子的4d轨道电子贡献。Mn掺杂后，在费米能级附近，Mn原子的3d轨道电子态密度出现了明显的峰值。这表明Mn原子的3d电子对费米能级附近的电子态产生了重要贡献，且与MoS₂中的电子发生了强烈的相互作用。这种相互作用导致了电子自旋的重新分布，使得体系产生了净磁矩。Mn原子的3d电子与周围S原子的电子形成了特定的自旋耦合状态，从而使整个体系表现出磁性。从磁学性质角度分析，本征二维MoS₂不具有磁性，而Mn掺杂后，体系产生了明显的磁矩。计算得到的Mn掺杂MoS₂体系的磁矩约为3.0μB，磁矩主要来源于Mn原子的3d电子。Mn原子的5个未成对3d电子在MoS₂晶格中形成了局域磁矩，且这些局域磁矩之间通过RKKY相互作用等机制产生了长程磁有序。RKKY相互作用是通过传导电子作为媒介，使得不同磁性原子之间产生磁相互作用。在Mn掺杂MoS₂体系中，Mn原子的局域磁矩通过与周围S原子的电子相互作用，再通过传导电子的间接作用，与其他Mn原子的局域磁矩相互耦合，形成了长程磁有序，使材料表现出宏观磁性。3.3不同掺杂浓度下的磁学性质演变3.3.1浓度梯度设置为深入探究掺杂浓度对二维MoS₂磁学性质的影响，精心设计并构建了一系列具有不同掺杂浓度的二维MoS₂模型。在模型构建过程中，以2×2×1的MoS₂超胞为基础，通过精确控制替换Mo原子的数量来实现不同的掺杂浓度。具体而言，设置了5%、10%、15%和20%这四个具有代表性的掺杂浓度梯度。在5%掺杂浓度的模型中，每20个原子中有1个Mo原子被掺杂原子（如Mn）替换；在10%掺杂浓度下，每10个原子中有1个Mo原子被替换；15%掺杂浓度对应每20个原子中有3个Mo原子被替换；20%掺杂浓度则是每5个原子中有1个Mo原子被替换。这种浓度梯度的设置既涵盖了低浓度掺杂的情况，能够研究少量掺杂原子对体系的影响，又包含了较高浓度掺杂的情形，以探索在较高掺杂水平下体系磁学性质的变化规律。同时，这种梯度设置也与以往相关研究中常用的浓度范围相呼应，便于与其他研究结果进行对比和分析。3.3.2磁矩与磁性变化规律通过第一性原理计算，深入分析了不同掺杂浓度下二维MoS₂体系的磁矩变化情况。计算结果表明，随着掺杂浓度的逐渐增加，体系的磁矩呈现出明显的上升趋势。在5%的低掺杂浓度下，体系的磁矩相对较小，约为0.5μB。这是因为在低掺杂浓度时，掺杂原子的数量较少，它们之间的磁相互作用相对较弱，且与周围MoS₂晶格中的电子相互作用范围也有限，导致体系整体磁矩较小。随着掺杂浓度提升至10%，磁矩增大到约1.2μB。此时，掺杂原子数量的增加使得它们之间的距离相对减小，磁相互作用增强，更多的电子参与到自旋极化过程中，从而导致体系磁矩显著增大。当掺杂浓度达到15%时，磁矩进一步增大至约2.0μB。在这个浓度下，掺杂原子在体系中的分布更加密集，它们之间的磁相互作用进一步增强，形成了更广泛的自旋极化区域，使得体系磁矩进一步增大。在20%的高掺杂浓度下，磁矩达到约3.0μB。此时，掺杂原子的相互作用占据主导地位，体系中形成了较为稳定的磁有序结构，磁矩达到相对较大的值。从磁性变化规律来看，随着掺杂浓度的增加，体系的磁性逐渐增强。在低掺杂浓度阶段，体系的磁性较弱，表现出一定的顺磁性特征。随着掺杂浓度的提高，体系逐渐向铁磁性转变。当掺杂浓度达到一定程度时，体系呈现出明显的铁磁性，具有较高的磁化强度和磁滞回线。这种磁性演变规律对二维MoS₂在自旋电子学领域的应用具有重要影响。在低掺杂浓度下，虽然体系磁性较弱，但可能适用于一些对磁性要求不高，但需要利用其电学和光学性质与微弱磁性耦合的应用场景，如某些传感器的设计。而在高掺杂浓度下，具有明显铁磁性的二维MoS₂则有望应用于自旋存储器件中，利用其稳定的磁有序结构来实现信息的存储和读取。四、MnO磁学性质的第一性原理研究4.1MnO的晶体结构与电子性质4.1.1晶体结构特征MnO具有典型的岩盐结构，属于立方晶系，空间群为Fm-3m。其晶体结构中，Mn²⁺离子和O²⁻离子交替排列，形成面心立方晶格。在这种结构中，每个Mn²⁺离子被六个O²⁻离子以八面体配位的方式紧密包围，同理，每个O²⁻离子也被六个Mn²⁺离子以八面体配位的方式环绕。这种离子排列方式赋予MnO晶体高度的对称性和稳定性。通过实验测量和理论计算，确定MnO的晶格常数a约为4.44Å。这一晶格常数反映了MnO晶体中离子间的平均距离，对其物理性质有着重要影响。在面心立方晶格中，Mn²⁺离子的坐标位置可以表示为(0,0,0)、(0.5,0.5,0)、(0.5,0,0.5)和(0,0.5,0.5)等，而O²⁻离子的坐标则为(0.25,0.25,0.25)、(0.75,0.75,0.25)、(0.75,0.25,0.75)和(0.25,0.75,0.75)等。这种精确的坐标设置与MnO的岩盐结构特征高度契合，确保了晶体结构的准确性。为了满足第一性原理计算中周期性边界条件的要求，本研究对构建好的MnO原胞模型进行了适当扩展。将原胞扩展为2×2×2的超胞，在扩展后的超胞中，原子的坐标位置相应地进行了调整。例如，Mn²⁺离子在扩展超胞中的坐标会根据原胞的重复和对称性进行相应的变化，以保持晶体结构的完整性和对称性。通过这种方式，成功构建了适用于计算的MnO晶体结构模型，为后续深入研究MnO的电子结构和磁学性质奠定了坚实的基础。4.1.2电子结构计算结果运用基于密度泛函理论的第一性原理计算方法，对MnO的电子结构进行了深入研究，得到了其能带结构和态密度图。MnO的能带结构较为复杂，反映了其内部电子的能量分布和运动状态。价带顶和导带底之间存在一定的能量间隔，表明MnO是一种具有一定带隙的半导体材料。计算得到的带隙宽度约为3.6eV，这一数值与实验测量结果以及其他理论计算结果基本相符。在能带结构中，价带主要由O原子的2p轨道电子和Mn原子的3d轨道电子贡献，导带则主要由Mn原子的3d轨道电子贡献。这种电子轨道的贡献分布表明，MnO的电子结构中，O原子和Mn原子之间存在着强烈的相互作用。态密度图能够更直观地展示电子在不同能量状态下的分布情况。从总态密度图可以看出，在费米能级附近，电子态密度存在明显的峰值，这与能带结构中价带顶和导带底的位置相对应。分波态密度图进一步揭示了不同原子轨道对电子态的贡献。在价带区域，O原子的2p轨道电子态密度占据主导地位，同时Mn原子的3d轨道电子态密度也有一定的贡献。这表明在价带中，O原子的2p电子与Mn原子的3d电子之间存在着较强的杂化作用。在导带区域，Mn原子的3d轨道电子态密度起主要作用，说明导带中的电子主要来源于Mn原子的3d轨道。这种电子态密度的分布特征与能带结构的分析结果相互印证，进一步说明了MnO的电子结构特点。4.1.3磁学性质分析MnO本征具有反铁磁性质，这一特性与其晶体结构和电子结构密切相关。在MnO的晶体结构中，Mn²⁺离子的3d电子对其磁学性质起着关键作用。Mn²⁺离子的3d轨道上有5个电子，这些电子的自旋和轨道角动量会产生磁矩。在反铁磁状态下，相邻Mn²⁺离子的磁矩方向相反，从而导致整体磁矩相互抵消，宏观上表现为零磁矩。从电子结构角度分析，MnO的反铁磁性质源于其电子间的交换相互作用。在MnO晶体中，Mn²⁺离子之间通过O²⁻离子产生超交换相互作用。这种超交换相互作用使得相邻Mn²⁺离子的磁矩呈现反平行排列，形成反铁磁有序结构。具体来说，O²⁺离子的2p电子与相邻Mn²⁺离子的3d电子之间存在着特定的轨道重叠和电子云分布，导致了磁矩的反平行排列。这种超交换相互作用是MnO反铁磁性质的重要来源。MnO的奈尔温度（TN）约为118K。在奈尔温度以下，MnO保持反铁磁有序状态；当温度升高超过奈尔温度时，热运动的能量足以破坏磁矩的反平行排列，MnO会逐渐转变为顺磁状态。这种温度对磁学性质的影响与MnO的晶体结构和电子结构的热稳定性密切相关。随着温度的升高，原子的热振动加剧，电子云的分布和相互作用也会发生变化，从而导致磁有序状态的改变。四、MnO磁学性质的第一性原理研究4.2外界因素对MnO磁学性质的影响4.2.1压力作用下的磁学性质变化为深入探究压力对MnO磁学性质的影响，运用第一性原理计算方法，模拟了不同压力条件下MnO的晶体结构、电子结构以及磁学性质。在计算过程中，系统地设置了一系列压力值，包括0GPa（常压）、10GPa、20GPa、30GPa和40GPa，以全面研究压力从低到高变化时MnO性质的演变规律。随着压力的逐渐增加，MnO的晶体结构发生了显著变化。晶格常数a呈现出单调减小的趋势。在常压下，MnO的晶格常数a约为4.44Å，当压力增加到10GPa时，晶格常数a减小至约4.38Å。这是因为压力的作用使得Mn²⁺离子和O²⁻离子之间的距离缩短，离子间的相互作用增强，从而导致晶格收缩。当压力进一步增加到40GPa时，晶格常数a减小至约4.20Å。这种晶格常数的变化对MnO的电子结构和磁学性质产生了深远的影响。从电子结构角度来看，压力的改变导致了能带结构和态密度的变化。随着压力的增加，MnO的能带宽度逐渐变宽。在常压下，价带顶和导带底之间的能量间隔（带隙）约为3.6eV。当压力达到10GPa时，带隙宽度减小至约3.4eV。这是由于压力使得离子间的距离减小，电子云的重叠程度增加，电子的能量状态发生改变，导致带隙减小。态密度图也显示出明显的变化。在费米能级附近，随着压力的增加，电子态密度的峰值发生了移动和变化。在常压下，价带主要由O原子的2p轨道和Mn原子的3d轨道电子贡献，导带主要由Mn原子的3d轨道电子贡献。随着压力的增加，O原子的2p轨道和Mn原子的3d轨道电子态密度的重叠程度增加，表明电子在O和Mn原子之间的相互作用增强。压力对MnO的磁学性质也产生了重要影响。计算结果表明，随着压力的增加，MnO的磁矩逐渐减小。在常压下，MnO的磁矩主要来源于Mn²⁺离子的3d电子，每个Mn²⁺离子的磁矩约为5μB。当压力增加到10GPa时，Mn²⁺离子的磁矩减小至约4.8μB。这是因为压力改变了Mn²⁺离子的电子云分布和自旋状态，使得磁矩减小。随着压力的进一步增加，磁矩继续减小。当压力达到40GPa时，Mn²⁺离子的磁矩减小至约4.2μB。压力对MnO的磁各向异性也有影响。在常压下，MnO具有一定的磁各向异性，随着压力的增加，磁各向异性逐渐减小。这是由于压力使得晶体结构的对称性发生变化，从而影响了磁各向异性。4.2.2温度对磁学性质的影响采用准谐近似方法，深入研究了温度对MnO磁学性质的影响。在研究过程中，系统地考虑了温度从0K到300K的变化范围，以全面揭示温度对MnO磁学性质的影响规律。随着温度的升高，MnO的晶体结构逐渐发生变化。晶格常数a呈现出逐渐增大的趋势。在0K时，MnO的晶格常数a约为4.43Å，当温度升高到300K时，晶格常数a增大至约4.46Å。这是因为温度升高，原子的热振动加剧，原子间的平均距离增大，导致晶格膨胀。这种晶格常数的变化对MnO的电子结构和磁学性质产生了重要影响。从电子结构角度来看，温度的变化导致了能带结构和态密度的改变。随着温度的升高，MnO的能带宽度逐渐变窄。在0K时，价带顶和导带底之间的能量间隔（带隙）约为3.65eV。当温度升高到300K时，带隙宽度减小至约3.55eV。这是由于温度升高，原子的热振动增强，电子云的分布变得更加弥散，电子的能量状态发生变化，导致带隙减小。态密度图也显示出明显的变化。在费米能级附近，随着温度的升高，电子态密度的峰值发生了移动和变化。在0K时，价带主要由O原子的2p轨道和Mn原子的3d轨道电子贡献，导带主要由Mn原子的3d轨道电子贡献。随着温度的升高，O原子的2p轨道和Mn原子的3d轨道电子态密度的重叠程度略有增加，表明电子在O和Mn原子之间的相互作用略有增强。温度对MnO的磁学性质影响显著。随着温度的升高，MnO的磁矩逐渐减小。在0K时，MnO处于反铁磁有序状态，Mn²⁺离子的磁矩为5μB。当温度升高到100K时，由于热运动的影响，部分Mn²⁺离子的磁矩方向发生了偏离，导致整体磁矩略有减小，每个Mn²⁺离子的磁矩约为4.9μB。当温度接近奈尔温度（TN约为118K）时，热运动的能量足以破坏磁矩的反平行排列，磁矩急剧减小。当温度超过奈尔温度时，MnO转变为顺磁状态，磁矩趋近于零。通过对MnO居里温度的计算，发现其居里温度约为118K，与实验值相符。在居里温度以下，MnO保持反铁磁有序，磁矩呈现出规则的反平行排列。当温度逐渐接近居里温度时，热涨落的影响逐渐增强，磁矩的反平行排列逐渐被破坏，磁有序度逐渐降低。当温度超过居里温度时，磁矩的反平行排列完全被破坏，MnO进入顺磁状态，磁矩不再具有长程有序性。4.3MnO与其他材料复合后的磁学性能4.3.1复合材料模型构建为研究MnO与其他材料复合后的磁学性能，选择与石墨烯构建复合材料模型。在构建模型时，首先考虑两者的结构特点和相互作用方式。MnO具有岩盐结构，属于立方晶系，而石墨烯是由碳原子组成的二维蜂窝状晶格结构。在VASP软件中，采用超胞方法构建MnO/石墨烯复合材料模型。对于MnO，选取2×2×2的超胞，以充分体现其晶体结构特征和原子间相互作用。对于石墨烯，构建一个包含16个碳原子的2×2超胞，以保证模型的完整性和计算的准确性。将MnO超胞和石墨烯超胞在垂直方向上进行组合，使MnO与石墨烯之间形成异质结构。在构建过程中，仔细调整两者的相对位置和取向，以模拟实际的复合情况。通过多次优化和尝试，确定MnO的（111）面与石墨烯的平面平行，且MnO中的Mn原子与石墨烯中的C原子之间的距离为3.5Å，这一距离既保证了两者之间存在一定的相互作用，又避免了原子间的过度重叠。为了准确描述MnO与石墨烯之间的相互作用，在计算中采用了范德华力修正。由于MnO与石墨烯之间主要通过范德华力相互作用，这种修正能够更真实地反映两者之间的结合情况。在VASP计算中，选用DFT-D3方法来考虑范德华力的影响。通过这种方法，对MnO/石墨烯复合材料模型进行结构优化，使体系达到能量最低的稳定状态。优化后的结构参数表明，MnO与石墨烯之间的结合较为稳定，界面处的原子位置和键长、键角等结构参数发生了一定的变化。这些变化反映了MnO与石墨烯之间的相互作用对结构的影响，为后续研究复合材料的磁学性质奠定了基础。4.3.2复合体系的磁学性质分析通过第一性原理计算，深入分析MnO/石墨烯复合体系的磁学性质。计算结果表明，复合体系的磁学性质与单一的MnO相比发生了显著变化。从磁矩角度来看，在单一的MnO中，Mn²⁺离子的磁矩主要来源于其3d电子，每个Mn²⁺离子的磁矩约为5μB。在MnO/石墨烯复合体系中，由于MnO与石墨烯之间的界面相互作用，Mn²⁺离子的磁矩发生了改变。部分Mn²⁺离子的磁矩减小至约4.8μB。这是因为石墨烯中的π电子与MnO中的Mn²⁺离子的3d电子发生了一定程度的相互作用，导致电子云分布发生变化，从而影响了Mn²⁺离子的磁矩。从磁相互作用类型来看，单一的MnO具有反铁磁性质，相邻Mn²⁺离子的磁矩方向相反。在复合体系中，磁相互作用变得更为复杂。虽然整体上仍然以反铁磁相互作用为主，但在MnO与石墨烯的界面附近，由于界面处原子的电子结构和自旋状态的改变，出现了一定程度的铁磁相互作用。这种铁磁相互作用的出现与界面处的电荷转移和自旋极化有关。通过电荷密度差分图可以清晰地看到，在MnO与石墨烯的界面处，存在明显的电荷转移现象，电子从MnO向石墨烯转移，导致界面处的电子自旋极化，从而产生了铁磁相互作用。这种界面相互作用对磁性的影响具有重要的应用潜力。在自旋电子学领域，MnO/石墨烯复合体系的这种独特磁学性质可用于设计新型的自旋电子学器件。利用界面处的铁磁相互作用和反铁磁相互作用的共存，可以实现对自旋极化电流的有效调控，为开发高性能的自旋场效应晶体管和磁随机存取存储器等器件提供了新的思路。在传感器领域，复合体系的磁学性质变化使其对外部磁场的响应更加敏感，有望用于制备高灵敏度的磁传感器，用于检测微弱的磁场变化和生物分子的磁信号。五、二维MoS₂与MnO磁学性质的对比与应用前景5.1两者磁学性质的对比分析5.1.1磁性来源差异二维MoS₂的磁性来源主要依赖于外部引入的杂质原子或缺陷。本征的二维MoS₂由于其原子结构和电子分布的高度对称性，电子配对良好，不存在未成对电子，因此不具有磁性。当通过掺杂过渡金属元素（如Mn、Fe等）或引入缺陷时，会改变MoS₂的电子结构。以Mn掺杂为例，Mn原子具有未填满的3d电子壳层，其3d电子的自旋和轨道角动量会产生磁矩。Mn原子掺杂到MoS₂晶格中后，与周围的Mo和S原子的电子发生相互作用，导致电子云分布发生变化，产生局域磁矩。这种局域磁矩之间通过RKKY相互作用等机制，在一定条件下可以形成长程磁有序，从而使二维MoS₂表现出磁性。MnO的磁性则源于其晶体结构中Mn²⁺离子的3d电子。MnO具有岩盐结构，Mn²⁺离子在晶体中呈周期性排列。Mn²⁺离子的3d轨道上有5个电子，这些电子的自旋和轨道角动量产生磁矩。在MnO中，相邻Mn²⁺离子之间通过O²⁻离子产生超交换相互作用。这种超交换相互作用使得相邻Mn²⁺离子的磁矩呈现反平行排列，形成反铁磁有序结构。在这种结构中，虽然每个Mn²⁺离子都具有磁矩，但由于相邻磁矩的反平行排列，宏观上整体磁矩相互抵消，表现为零磁矩。与二维MoS₂通过外部引入磁性的方式不同，MnO的磁性是其本征晶体结构和电子结构所决定的固有属性。5.1.2磁学性能特点比较在饱和磁化强度方面，两者表现出明显差异。对于二维MoS₂，其饱和磁化强度与掺杂元素的种类、浓度以及缺陷的类型和密度密切相关。在低掺杂浓度下，由于磁性原子数量较少，局域磁矩之间的相互作用较弱，饱和磁化强度相对较低。随着掺杂浓度的增加，磁性原子之间的距离减小，相互作用增强，饱和磁化强度逐渐增大。在某些高掺杂浓度的情况下，二维MoS₂的饱和磁化强度可以达到一定的值，但总体来说，由于其磁性主要依赖于外部引入，饱和磁化强度相对有限。MnO作为反铁磁材料，其宏观饱和磁化强度为零。在反铁磁状态下，相邻Mn²⁺离子的磁矩方向相反，相互抵消，使得在宏观上没有净磁矩。只有在施加足够强的外磁场，克服反铁磁相互作用，使磁矩方向发生改变时，才会表现出一定的磁化强度，但这与一般意义上的饱和磁化强度概念不同。从居里温度来看，二维MoS₂由于其低维结构和磁性来源的特殊性，其居里温度通常较低。在一些掺杂的二维MoS₂体系中，居里温度可能在几十K到几百K之间，这限制了其在高温环境下的磁性应用。而MnO的奈尔温度（类似于反铁磁材料的居里温度）约为118K。虽然MnO在奈尔温度以下保持反铁磁有序，但与一些具有较高居里温度的铁磁材料相比，其磁有序状态对温度的稳定性相对较低。当温度升高接近奈尔温度时，热运动的能量足以破坏磁矩的反平行排列，MnO会逐渐转变为顺磁状态。五、二维MoS₂与MnO磁学性质的对比与应用前景5.2在自旋电子学领域的应用潜力5.2.1二维MoS₂在自旋电子器件中的应用设想二维MoS₂由于其独特的电学和磁学性质，在自旋电子器件领域展现出广阔的应用前景。在自旋晶体管的设计中，二维MoS₂可作为沟道材料。传统的晶体管主要基于电荷输运原理，而自旋晶体管则利用电子的自旋自由度来实现信息的处理和传输。二维MoS₂的原子级厚度使其具有较强的量子限域效应和自旋轨道耦合作用。通过外部电场的调控，可以有效地控制电子的自旋极化方向和输运特性。当施加一定的电场时，二维MoS₂中的电子自旋可以发生取向变化，从而实现自旋电流的产生和调控。这种基于自旋的信息传输方式具有低功耗、高速和高集成度的优势，有望为未来的高速计算和低功耗电子设备提供关键技术支持。在磁存储领域，二维MoS₂也具有潜在的应用价值。利用二维MoS₂的磁性和电学性质的耦合效应，可以设计新型的磁随机存取存储器（MRAM）。在MRAM中，信息的存储和读取是通过磁性材料的磁化方向来实现的。二维MoS₂的磁性可以通过掺杂等手段进行调控，使其具备稳定的磁矩和磁各向异性。通过控制外部磁场或电场，可以改变二维MoS₂的磁化方向，从而实现信息的写入和读取。与传统的存储技术相比，基于二维MoS₂的MRAM具有非易失性、高速读写和低功耗等优点，能够满足现代信息技术对存储设备高性能、小型化和低功耗的需求。5.2.2MnO在自旋电子学中的潜在应用MnO在自旋电子学领域也展现出潜在的应用可能性。在自旋阀的设计中，MnO可作为反铁磁层。自旋阀是一种基于巨磁电阻效应的自旋电子学器件，由两层铁磁层和一层非磁性层组成。MnO的反铁磁性质使其能够有效地钉扎相邻铁磁层的磁矩方向，从而提高自旋阀的磁电阻效应和稳定性。在自旋阀中，当外部磁场发生变化时，两层铁磁层的磁矩方向会发生相对变化，导致电阻发生显著变化。MnO的反铁磁层可以通过与铁磁层之间的交换耦合作用，稳定铁磁层的磁矩方向，使得自旋阀在不同的磁场条件下能够保持稳定的电阻状态，提高器件的可靠性和性能。在磁性隧道结方面，MnO也具有应用潜力。磁性隧道结是由两层铁磁层和一层绝缘层组成的三明治结构。当两层铁磁层的磁矩方向平行时，自旋电子可以透过绝缘层进行隧穿，形成低阻状态；当磁矩方向反平行时，隧穿会被阻挡，形成高阻状态。MnO可以作为绝缘层或磁性层应用于磁性隧道结中。作为绝缘层时，MnO的晶体结构和电子性质可以影响隧道结的隧穿特性和磁电阻效应。其带隙结构和电子云分布会影响电子的隧穿概率，从而影响隧道结的电阻变化。作为磁性层时，MnO的反铁磁性质可以与其他铁磁层形成特定的磁相互作用，实现对自旋电子的有效调控。然而，MnO在自旋电子学应用中也面临一些挑战。其反铁磁性质使得磁矩的调控相对困难，需要更强的外场或特殊的调控手段来实现磁矩的改变。MnO与其他材料的界面兼容性也是需要解决的问题，界面处的晶格失配和电子相互作用可能会影响器件的性能和稳定性。五、二维MoS₂与MnO磁学性质的对比与应用前景5.3在其他领域的潜在应用探讨5.3.1传感器应用二维MoS₂在传感器领域展现出巨大的应用潜力。在磁传感器方面，由于其原子级厚度和独特的电学与磁学性质，二维MoS₂对微弱磁场具有高度敏感性。当外部磁场发生变化时，二维MoS₂的电子结构会相应改变，导致其电学性质如电阻、电容等发生变化。这种磁电耦合效应使得二维MoS₂可用于制备高灵敏度的磁传感器，能够检测极其微弱的磁场信号。在生物医学领域，可用于检测生物分子的磁信号，实现对生物分子的高灵敏度检测和分析。在生物传感器应用中，二维MoS₂的大比表面积和良好的生物兼容性使其成为理想的生物传感材料。通过表面功能化修饰，可将特定的生物识别分子（如抗体、核酸等）固定在二维MoS₂表面。当目标生物分子与修饰后的二维MoS₂表面的识别分子发生特异性结合时，会引起二维MoS₂电学或光学性质的变化，从而实现对目标生物分子的检测。利用二维MoS₂修饰的场效应晶体管构建生物传感器，可实现对蛋白质、DNA等生物分子的高灵敏度检测，检测限可达到皮摩尔级甚至更低。MnO在传感器领域也具有独特的应用优势。基于MnO的磁电阻效应，可制备磁电阻传感器。MnO的反铁磁性质使其在磁场变化时，电阻会发生显著变化，这种磁电阻效应可用于检测磁场的强度和方向变化。在一些精密测量领域，如地磁测量、生物磁信号检测等，MnO磁电阻传感器能够提供高精度的磁场测量结果。MnO还可与其他材料复合，制备多功能传感器。将MnO与纳米线材料复合，可制备出具有高灵敏度的压力和磁场双响应传感器。在压力作用下，纳米线的形变会导致MnO的磁学性质发生变化，同时MnO对磁场的响应也会改变，通过检测这种双重变化，可实现对压力和磁场的同时检测，拓展了传感器的应用范围。5.3.2能源存储与转换应用在能源存储领域，二维MoS₂展现出潜在的应用价值。二维MoS₂具有较高的理论比容量，可作为锂离子电池电极材料。其独特的层状结构能够为锂离子的嵌入和脱出提供丰富的通道，有利于提高电池的充放电性能。在二维MoS₂中引入磁性后，磁学性质与电学性质的耦合可能进一步优化电池的性能。磁性的存在可能影响锂离子在材料中的扩散速率和电子传输效率，从而提高电池的充放电倍率和循环稳定性。一些研究表明，通过掺杂磁性原子制备的磁性二维MoS₂作为锂离子电池电极材料，在高电流密度下仍能保持较高的比容量，循环稳定性也得到显著提升。MnO在电磁感应领域具有重要的应用潜力。MnO的反铁磁性质使其在磁场变化时能够产生感应电流，可用于制备电磁感应器件。在变压器中，MnO可作为磁性材料用于增强电磁感应效率，提高能量转换效率。将MnO与其他磁性材料复合，可进一步优化电磁感应器件的性能。将MnO与铁氧体材料复合，制备的复合磁性材料在电磁感应过程中，能够实现更好的磁导率匹配和更低的磁损耗，提高变压器等电磁感应器件的性能和稳定性。六、结论与展望6.1研究成果总结本研究运用基于密度泛函理论的第一性原理计算方法，对二维MoS₂和MnO的磁学性质展开深入探究，成功揭示了二者磁学性质的内在机制及影响因素，取得了一系列具有重要理论与实践意义的研究成果。在二维MoS₂的研究中，明确了本征二维MoS₂由于其原子结构和电子分布的高度对称性，电子配对良好，不存在未成对电子，因此呈现非磁性。通过引入掺杂和缺陷等外部调控手段，成功实现了对其磁学性质的有效调控。当过渡金属元素如Mn掺杂到二维MoS₂中时，由于Mn原子具有未填满的3d电子壳层，其3d电子的自旋和轨道角动量产生磁矩，与周围的Mo和S原子的电子发生相互作用，改变了电子云分布，从而诱导出磁性。随着掺杂浓度的逐渐增加，体系的磁矩呈现出明显的上升趋势。在低掺杂浓度下，由于磁性原子数量较少，局域磁矩之间的相互作用较弱，体系磁矩较小。随着掺杂浓度的提升，磁性原子之间的距离减小，相互作用增强，更多的电子参与到自旋极化过程中，使得体系磁矩显著增大。通过对不同掺杂浓度下二维MoS₂体系的磁学性质演变进行系统研究，揭示了掺杂浓度与磁矩、磁性变化之间的定量关系，为二维MoS₂在自旋电子学等领域的应用提供了关键的理论依据。对于MnO，深入剖析了其晶体结构、电子结构与磁学性质之间的内在联系。MnO具有典型的岩盐结构，属于立方晶系，在这种结构中，Mn²⁺离子和O²⁻离子交替排列，形成面心立方晶格。每个Mn²⁺离子被六个O²⁻离子以八面体配位的方式包围，这种离子排列方式决定了MnO的晶体结构特征和物理性质。MnO的磁性源于其晶体结构中Mn²⁺离子的3d电子。Mn²⁺离子的3d轨道上有5个电子，这些电子的自旋和轨道角动量产生磁矩。相邻Mn²⁺离子之间通过O²⁻离子产生超交换相互作用，使得相邻Mn²⁺离子的磁矩呈现反平行排列，形成反铁磁有序结构。在反铁磁状态下，虽然每个Mn²⁺离子都具有磁矩，但由于相邻磁矩的反平行排列，宏观上整体磁矩相互抵消，表现为零磁矩。研究还发现，外界因素如压力和温度对MnO的磁学性质具有显著影响。随着压力的增加，MnO的晶格常数减小，能带宽度变宽，磁矩逐渐减小，磁各向异性也逐渐减小。这是因为压力改变了离子间的距离和电子云分布，影响了电子的能量状态和自旋状态。随着温度的升高，MnO的晶格常数增大，能带宽度变窄，磁矩逐渐减小。当温度接近奈尔温度（TN约为118K）时，热运动的能量足以破坏磁矩的反平行排列，磁矩急剧减小。当温度超过奈尔温度时，MnO转变为顺磁状态，磁矩趋近于零。此外，通过构建MnO与石墨烯的复合材料模型，研究发现复合体系的磁学性质与单一的MnO相比发生了显著变化。由于MnO与石墨烯之间的界面相互作用，Mn²⁺离子的磁矩发生了改变