基于第一性原理探究单层二碲化钼畴界缺陷特性及影响

基于第一性原理探究单层二碲化钼畴界缺陷特性及影响一、引言1.1研究背景与意义2004年，曼彻斯特大学Geim小组成功分离出单原子层的石墨材料——石墨烯，这一突破性成果开启了二维材料的研究热潮。二维材料，作为电子仅可在两个维度的非纳米尺度（1-100nm）上自由运动的材料，展现出诸多与传统三维材料截然不同的特性。其特殊的原子结构和电子态分布，赋予了材料在电学、磁学、光学和力学等多方面的优异性能，在半导体器件、能源存储与转换、传感器、催化等众多领域呈现出巨大的应用潜力。二维过渡金属二硫化物（TMDs）是二维材料家族中的重要成员，具有独特的晶体结构和电子性质。这类材料由过渡金属原子（如Mo、W等）与硫族原子（如S、Se、Te等）通过共价键结合而成，呈现出层状结构。在TMDs中，单层二碲化钼（MoTe₂）凭借其独特的物理性质，成为了研究的焦点之一。MoTe₂体材料是间接带隙半导体，然而，当厚度减薄至单层时，量子限域效应使其带隙转变为约1.10eV的直接带隙。这种直接带隙特性，不仅拓宽了其光谱响应范围至近红外波段，还使得单层MoTe₂在光电器件应用中展现出巨大优势。此外，室温下单层MoTe₂的载流子迁移率比MoS₂更高，这为其在高速电子学器件中的应用提供了可能。这些优异的性质，使得单层MoTe₂在晶体管、超导体、谷光电子器件以及光电探测器等领域具有广阔的应用前景。例如，在晶体管应用中，较高的载流子迁移率可提高器件的开关速度和工作频率；在谷光电子器件中，其独特的能谷特性有望实现信息的高效编码与传输。高质量的二维材料是实现高性能器件应用的基础。在实际制备过程中，无论是化学气相沉积（CVD）还是化学气相传输（CVT）等方法，都难以避免地会引入各种本征缺陷。点缺陷、畴界缺陷等在二维材料中普遍存在。这些缺陷的存在，对二维材料的电学、光学、磁学等性质产生着深远的影响。值得注意的是，缺陷对材料性质的影响并非完全负面。例如，在本征过渡金属硫化物中，由于其无磁性，在自旋电子学领域的应用受到极大限制。而二维TMDs部分缺陷处会产生较强的局域磁矩，这为自旋电子器件的应用开辟了新的途径。又如，某些缺陷可以作为活性位点，增强材料在催化反应中的活性；在光学性质方面，缺陷可能导致发光特性的改变，为新型发光材料的开发提供了可能。因此，深入研究二维材料中缺陷的性质，对于理解材料的本征特性、优化材料性能以及拓展材料应用具有至关重要的意义。单层MoTe₂中本征畴界缺陷的磁性和电子性质的研究相对较少。畴界缺陷作为一种特殊的线缺陷，其原子排列的不规则性会导致局部电子结构和能量状态的改变，进而影响材料的宏观性质。通过对单层MoTe₂畴界缺陷的深入研究，可以揭示缺陷与材料性能之间的内在联系，为调控材料性质提供理论依据。从理论层面来看，研究畴界缺陷有助于深化对二维材料中电子-晶格相互作用、电子关联效应等基础物理问题的理解，丰富二维材料的物理内涵；从应用角度出发，掌握畴界缺陷对材料性能的影响规律，能够为高质量单层MoTe₂材料的制备工艺优化提供指导，提高材料的质量和一致性，降低生产成本，推动其在各类器件中的实际应用。同时，通过对畴界缺陷的有效调控，还可能开发出具有独特性能的新型材料，满足未来科技发展对高性能材料的需求。1.2国内外研究现状自2004年石墨烯被首次成功制备以来，二维材料凭借其独特的原子结构和物理性质，在全球范围内引发了广泛而深入的研究热潮。二维过渡金属二硫化物（TMDs）作为二维材料家族的重要成员，由于其在电学、光学、力学等方面展现出的优异性能，在下一代纳米电子设备、能源存储与转换、传感器等众多领域具有重要的应用价值，成为了研究的焦点之一。在众多TMDs材料中，二碲化钼（MoTe₂）因其特殊的晶体结构和电子性质，受到了国内外科研人员的高度关注。在国外，美国麻省理工学院的研究团队在二维材料的基础研究方面取得了多项重要成果。他们通过先进的实验技术，深入研究了二维材料中电子与晶格的相互作用，揭示了一些独特的物理现象。例如，在对MoTe₂的研究中，利用高分辨率电子显微镜，观察到了MoTe₂晶体中原子的精确排列方式，以及在外部电场作用下电子态的变化，为理解MoTe₂的电学性质提供了直观的实验依据。哈佛大学的科研人员则致力于探索MoTe₂在量子器件中的应用潜力，通过巧妙的设计和制备工艺，成功构建了基于MoTe₂的原型量子器件，并对其性能进行了系统研究，发现了一些与传统器件不同的量子特性，为未来量子计算和通信技术的发展提供了新的思路。在国内，清华大学的研究小组在二维材料的生长和制备技术方面取得了突破性进展。他们通过优化化学气相沉积（CVD）工艺，成功制备出高质量、大面积的单层MoTe₂薄膜，并且在薄膜的均匀性和结晶质量方面达到了国际领先水平。这一成果为MoTe₂在实际器件中的应用奠定了坚实的材料基础。中国科学院物理研究所的科研团队则在MoTe₂的物理性质研究方面做出了重要贡献。他们运用多种先进的实验手段，如角分辨光电子能谱（ARPES）、扫描隧道显微镜（STM）等，对MoTe₂的电子结构、能带特性等进行了深入研究，精确测量了MoTe₂的能带结构和电子态密度，揭示了其能带反转和拓扑性质，为MoTe₂在拓扑量子计算领域的应用提供了理论支持。在二维材料的缺陷研究方面，国外的研究起步较早，取得了一系列具有影响力的成果。英国剑桥大学的研究团队通过理论计算和实验验证相结合的方法，系统研究了二维材料中点缺陷对电子输运性质的影响。他们建立了精确的理论模型，能够准确预测点缺陷的形成能和迁移率，以及对材料电学性能的影响机制。美国斯坦福大学的科研人员则专注于研究二维材料中畴界缺陷的结构和性质，利用先进的透射电子显微镜技术，成功观察到了畴界缺陷的原子结构和分布情况，为理解畴界缺陷的形成机制提供了关键的实验证据。国内的科研团队在二维材料缺陷研究领域也取得了显著的成绩。北京大学的研究小组通过第一性原理计算，深入研究了二维材料中缺陷与杂质之间的相互作用，揭示了杂质对缺陷形成和演化的影响规律，为通过杂质工程调控二维材料性能提供了理论指导。中国科学技术大学的科研团队则在实验上通过对二维材料进行精确的缺陷控制，制备出具有特定性能的二维材料器件，并对其性能进行了详细研究，展示了缺陷调控在优化二维材料器件性能方面的巨大潜力。然而，目前对于单层MoTe₂中本征畴界缺陷的磁性和电子性质的研究相对较少。虽然国内外学者在二维材料的整体研究中取得了众多成果，但针对单层MoTe₂这一特定材料的畴界缺陷研究仍存在许多空白。畴界缺陷作为一种特殊的线缺陷，其原子排列的不规则性会导致局部电子结构和能量状态的改变，进而影响材料的宏观性质。现有的研究主要集中在点缺陷对单层MoTe₂性质的影响，而对于畴界缺陷的研究还不够系统和深入。在磁性方面，虽然已经知道二维TMDs部分缺陷处会产生局域磁矩，但对于单层MoTe₂畴界缺陷处磁矩的产生机制、大小和方向的调控等问题，尚未有明确的结论。在电子性质方面，畴界缺陷对电子的散射机制、载流子迁移率以及能带结构的影响等，都有待进一步深入研究。深入研究单层MoTe₂中本征畴界缺陷的磁性和电子性质具有重要的理论和实际意义。从理论角度来看，这有助于深化对二维材料中电子-晶格相互作用、电子关联效应等基础物理问题的理解，丰富二维材料的物理内涵。从实际应用角度出发，掌握畴界缺陷对材料性能的影响规律，能够为高质量单层MoTe₂材料的制备工艺优化提供指导，提高材料的质量和一致性，降低生产成本，推动其在各类器件中的实际应用。1.3研究内容与方法本研究聚焦于单层二碲化钼（MoTe₂）畴界缺陷，运用第一性原理计算方法，深入剖析其电子性质和磁性，具体研究内容如下：单层MoTe₂的结构优化与稳定性分析：构建单层MoTe₂的初始原子结构模型，采用基于密度泛函理论（DFT）的第一性原理计算方法，对其进行结构优化。通过能量最小化算法，使体系达到稳定状态，得到晶格常数、原子坐标等结构参数，并与实验值和其他理论计算结果进行对比，验证模型的准确性。计算单层MoTe₂的形成能，评估其在不同条件下的稳定性，为后续研究奠定基础。畴界缺陷模型的构建与分类：依据实验观测和理论分析，构建多种类型的单层MoTe₂畴界缺陷模型，如晶界错配位错型畴界、堆垛层错型畴界等。考虑不同的畴界取向、原子排列方式以及缺陷间距，全面系统地研究畴界缺陷的结构特征对材料性质的影响。对所构建的畴界缺陷模型进行结构优化，确保其处于能量最低的稳定状态，分析缺陷区域的原子位移、键长和键角变化，揭示畴界缺陷的原子结构细节。畴界缺陷对电子性质的影响：利用第一性原理计算方法，研究畴界缺陷对单层MoTe₂电子结构的影响，包括能带结构、态密度、电荷密度分布等。分析畴界缺陷导致的能带畸变、能隙变化以及杂质能级的出现，探讨这些变化对材料电学性能的影响机制。通过计算电子局域函数（ELF）和差分电荷密度，研究畴界缺陷附近电子的局域化程度和电荷转移情况，揭示电子与缺陷的相互作用规律。结合载流子迁移率的计算，评估畴界缺陷对单层MoTe₂载流子输运性质的影响，为其在电子器件中的应用提供理论依据。畴界缺陷对磁性的影响：探究单层MoTe₂畴界缺陷处的磁性起源和性质，计算缺陷体系的总磁矩和原子磁矩，分析磁矩的分布和大小与缺陷结构的关系。通过自旋极化密度泛函理论（SP-DFT）计算，研究畴界缺陷对自旋极化状态的影响，揭示磁性与电子结构之间的内在联系。探讨不同类型畴界缺陷的磁性耦合机制，分析铁磁、反铁磁或其他磁性状态的稳定性，为基于单层MoTe₂的自旋电子学器件应用提供理论指导。外部电场对畴界缺陷性质的调控：在施加外部电场的条件下，研究单层MoTe₂畴界缺陷的电子性质和磁性的变化规律。通过改变电场的强度和方向，观察能带结构、态密度、电荷密度分布以及磁矩等物理量的响应，分析外部电场对畴界缺陷性质的调控机制。探索利用外部电场实现对单层MoTe₂畴界缺陷电学和磁学性能有效调控的可能性，为新型器件的设计提供新思路。本研究主要采用基于密度泛函理论的第一性原理计算方法，具体使用的计算软件为ViennaAbinitioSimulationPackage（VASP）。该软件基于平面波赝势方法，能够精确描述原子和电子之间的相互作用，广泛应用于材料科学领域的理论研究。在计算过程中，采用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）交换关联泛函，以准确描述电子的交换和关联能。平面波截断能设置为500eV，以保证计算精度。采用Monkhorst-Pack方法对布里渊区进行k点采样，对于超晶胞结构，k点网格设置为3×3×1，以平衡计算精度和计算效率。结构优化过程中，原子受力收敛标准设置为0.01eV/Å，能量收敛标准设置为10⁻⁵eV，以确保体系达到稳定状态。通过这些计算参数的合理设置，能够准确地研究单层MoTe₂畴界缺陷的电子性质和磁性。二、第一性原理与相关理论基础2.1第一性原理概述第一性原理（FirstPrinciples），最初源于量子力学领域，是一种基于基本物理规律和原理，从最基本的条件和规则出发，不依赖于经验参数或实验数据，通过求解薛定谔方程等基本方程来描述和模拟物质与系统行为的计算方法。其基本思想可追溯到古希腊哲学家亚里士多德的哲学理念，他认为每个系统中都存在一个最基本的、不可简化的命题，这一命题构成了系统的基础，是决定事物本质的不变法则。在现代科学语境下，第一性原理被广泛应用于物理学、化学、材料科学等众多领域，成为研究微观世界和揭示物质基本性质的重要工具。在材料科学研究中，第一性原理计算方法具有独特的优势。传统的材料研究往往依赖于实验试错法，这种方法不仅耗时费力，而且成本高昂。而第一性原理计算则为材料研究提供了一种全新的思路和方法。它能够从原子和电子层面深入理解材料的结构与性能之间的内在联系，预测材料的各种性质，如电子结构、晶体结构、力学性质、磁性、光学性质等。通过第一性原理计算，研究人员可以在理论层面上对材料进行设计和优化，筛选出具有潜在应用价值的材料，为实验研究提供指导，从而大大缩短材料研发周期，降低研发成本。第一性原理计算的核心在于对薛定谔方程的求解。对于由多个原子和电子组成的体系，其哈密顿量（Hamiltonian）可以表示为电子的动能、电子与原子核之间的相互作用势能以及电子之间的相互作用势能之和。薛定谔方程的一般形式为：H\psi=E\psi其中，H为哈密顿算符，\psi为波函数，E为体系的能量本征值。波函数\psi描述了体系中电子的状态，通过求解薛定谔方程，可以得到体系的能量本征值和波函数，进而计算出体系的各种物理性质。然而，对于实际的多电子体系，直接求解薛定谔方程是极其困难的，甚至是不可能的，因为电子之间存在着复杂的相互作用，使得方程的求解变得极为复杂。为了简化计算，人们引入了一系列的近似方法，其中最常用的是密度泛函理论（DensityFunctionalTheory，DFT）。DFT的核心思想是将多电子体系的能量表示为电子密度的泛函，而不是波函数。Hohenberg-Kohn定理为DFT提供了坚实的理论基础，该定理表明，基态系统的所有物理性质都由电子密度唯一决定，能量与电子密度为一一映射关系。基于这一定理，Kohn-Sham提出了Kohn-Sham方程，将多电子问题转化为单电子问题，从而大大简化了计算过程。在实际计算中，通过选择合适的交换关联泛函来近似描述电子之间的交换和关联能，进而求解Kohn-Sham方程，得到体系的电子结构和能量等物理性质。在具体计算过程中，第一性原理计算还涉及到诸多细节和技术。例如，为了描述原子与电子之间的相互作用，常采用赝势方法。由于原子核周围的电子分布较为复杂，尤其是内层电子与原子核之间的强相互作用，给计算带来了很大困难。赝势方法通过构建一个等效的势场，来代替真实的原子核-电子相互作用势，从而简化计算。在赝势中，将原子核及其内层电子视为一个整体，即离子实，只考虑外层价电子与离子实之间的相互作用，这样可以有效地减少计算量，同时保证计算结果的准确性。此外，在处理晶体结构时，常采用平面波基组展开电子波函数。平面波具有标准正交化和能量单一性的特点，对任何原子都适用且等同看待空间中的任何区域，不需要修正重叠误差。通过将电子波函数用平面波基组展开，可以将复杂的薛定谔方程转化为矩阵方程进行求解。但在实际计算中，由于平面波基组的数量是无限的，为了控制计算量，需要引入截断能的概念，只保留能量小于截断能的平面波进行计算。2.2密度泛函理论（DFT）密度泛函理论（DensityFunctionalTheory，DFT）是第一性原理计算中最为关键的理论基础之一，它的出现极大地推动了材料科学、化学等领域的理论研究进展。其核心思想在于将多电子体系的能量表示为电子密度的泛函，而非传统的波函数。这一理论的提出，为解决多电子体系的复杂问题提供了一种全新的思路和方法。Hohenberg-Kohn定理是密度泛函理论的基石，该定理包含两个重要内容。其一，基态系统的所有物理性质都由电子密度唯一决定，这意味着只要确定了电子密度，就可以通过相应的泛函计算出体系的各种物理性质，如能量、电荷密度分布、能带结构等。这种一一对应的关系，使得研究人员可以将复杂的多电子体系问题转化为对电子密度的研究，大大简化了计算过程。其二，对应于电子密度的变分原理表明，任意近似电子密度所对应的能量值都大于等于基态对应的真正密度所决定的能量值。这一原理为寻找体系的基态能量提供了一种有效的方法，即通过不断调整电子密度，使得体系能量达到最小值，从而得到基态能量和基态电子密度。在密度泛函理论中，体系的总能量可以表示为多个部分的和，即E[\rho]=T[\rho]+E_{ne}[\rho]+J[\rho]+K[\rho]。其中，T[\rho]代表电子的动能，它描述了电子在体系中的运动状态；E_{ne}[\rho]是核子-电子势能，体现了原子核与电子之间的相互作用；J[\rho]为库仑积分，反映了电子之间的库仑排斥作用；K[\rho]是交换积分，考虑了电子的交换关联效应。在这些能量项中，E_{ne}[\rho]和J[\rho]的形式是已知的，可以通过经典的静电学理论进行计算。然而，T[\rho]和K[\rho]的具体形式较为复杂，难以精确确定，这也是密度泛函理论在实际应用中需要解决的关键问题。为了求解多电子体系的薛定谔方程，Kohn-Sham提出了Kohn-Sham方程，将多电子问题巧妙地转化为单电子问题。该方程的核心思想是引入一组等效的单电子波函数，通过求解这些单电子波函数来得到体系的电子结构和能量。具体来说，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})是第i个单电子波函数，\epsilon_{i}是对应的本征能量。通过求解Kohn-Sham方程，可以得到单电子波函数和本征能量，进而计算出体系的电子密度、总能量等物理量。交换关联泛函在密度泛函理论中起着至关重要的作用，它直接影响着计算结果的准确性。由于交换关联效应的复杂性，目前尚无法精确求解交换关联泛函，因此在实际计算中通常采用各种近似方法。常见的交换关联泛函包括局域密度近似（LocalDensityApproximation，LDA）、广义梯度近似（GeneralizedGradientApproximation，GGA）等。LDA假设局域电子密度可以被认为是均匀电子气，或等效地说，电子密度是随空间缓慢变化的函数。在这种近似下，交换关联能只与电子密度的局域值有关。LDA在处理一些简单体系时表现出较好的性能，能够给出与实验结果较为接近的计算结果。然而，对于一些复杂体系，如过渡金属氧化物、强关联体系等，LDA的计算结果往往与实验值存在较大偏差。GGA则对LDA进行了改进，它考虑了电子密度的梯度信息，认为交换关联能不仅与电子密度的局域值有关，还与电子密度的变化率有关。GGA中的Perdew-Wang（PW91）泛函通过修改Vosko，Wilk，andNusair（VWN）的泛函形式，加入了高阶项，能够更准确地描述电子的交换关联效应。在处理一些具有复杂电子结构的体系时，GGA通常能够给出比LDA更精确的计算结果。例如，在研究过渡金属化合物的电子结构和磁性时，GGA能够更好地描述过渡金属离子的d电子之间的相互作用，从而得到与实验更为相符的结果。除了LDA和GGA，还有一些其他的交换关联泛函，如混合泛函（HybridFunctionals）等。混合泛函将Hartree-Fock理论中的精确交换项与密度泛函理论中的交换关联项相结合，试图在提高计算精度的同时，保持计算效率。在一些对精度要求较高的计算中，混合泛函能够发挥重要作用，但由于其计算量相对较大，在实际应用中受到一定的限制。2.3平面波赝势方法平面波赝势方法（PlaneWavePseudopotentialMethod）是基于密度泛函理论的第一性原理计算中广泛应用的一种方法，它为解决多电子体系的复杂计算问题提供了一种高效且精确的途径。在固体材料中，原子核周围存在内层电子和外层价电子。内层电子紧密束缚于原子核，其波函数在原子核附近变化极为剧烈，这使得在计算中精确描述其行为需要大量的平面波基函数，从而导致计算量呈指数级增长。为了克服这一困难，平面波赝势方法引入了赝势的概念。赝势通过构建一个等效势场，将原子核及其内层电子视为一个整体，即离子实，仅考虑外层价电子与离子实之间的相互作用。在赝势的作用下，价电子感受到的是一个经过平滑处理的势场，其波函数在离子实区域的变化变得相对平缓，不再需要大量的平面波来描述。这不仅大大减少了计算所需的平面波数量，降低了计算量，还能够保证计算结果的准确性。例如，在对过渡金属材料的计算中，内层电子的复杂波函数会给计算带来极大挑战，而引入赝势后，能够有效地简化计算过程，使得对这类材料的研究成为可能。平面波赝势方法中，采用平面波基组展开电子波函数是其另一个关键步骤。根据晶体的空间平移对称性，布洛赫（Bloch）定理表明，能带电子的波函数总可以写成\psi_{k}(r)=\varphi_{k}(r)e^{ik\cdotr}的形式，其中k是电子波矢，\varphi_{k}(r)是具有晶体平移周期性的周期函数。在周期性边界条件下，单粒子轨道波函数可以用平面波基展开为\psi(r)=\frac{1}{\sqrt{N_{cell}}}\sum_{G}\alpha(G)e^{i(K+G)\cdotr}，其中\frac{1}{\sqrt{N_{cell}}}是归一化因子，N_{cell}是原胞体积；G是原胞的倒格矢，K是第一布里渊区的波矢，\alpha(G)是展开系数。平面波具有标准正交化和能量单一性的特点，对任何原子都适用，且等同看待空间中的任何区域，不需要修正重叠误差。这使得平面波函数基组适用于许多体系，为求解Kohn-Sham方程提供了一种高效的方案。同时，平面波基能很方便地采用快速傅里叶变换（FFT）技术，使能量、力等的计算在实空间和倒空间迅速转换，进一步提高了计算效率。在实际计算中，由于平面波基组的数量是无限的，为了控制计算量，需要引入截断能（EnergyCutoff）的概念。给定一个截断能E_{cut}=\frac{\hbar^{2}(G+K)^{2}}{2m}，对G的求和可以限制在(G+K)^{2}/2\leqE_{cut}的范围内，即要求用于展开的波函数的能量小于E_{cut}。截断能的选择对计算结果的精度有着重要影响。当截断能较低时，参与计算的平面波数量较少，计算量小，但可能会导致计算结果的精度不足；随着截断能的增加，参与计算的平面波数量增多，计算精度会提高，但计算量也会相应增大。因此，在实际计算中，需要通过测试不同的截断能，找到一个既能保证计算精度，又能控制计算量的合适值。例如，在对硅晶体的计算中，通过逐步增加截断能，观察体系总能量、电子密度等物理量的变化，当这些物理量在不同截断能下的变化小于一定阈值时，即可确定该截断能为合适的值。平面波赝势方法通过引入赝势简化了原子核与电子之间的相互作用，采用平面波基组展开电子波函数，并通过合理选择截断能平衡计算精度和计算量，为研究材料的电子结构和性质提供了一种强大的工具。在本研究中，采用平面波赝势方法能够准确地描述单层二碲化钼及其畴界缺陷体系中电子的行为，为深入研究其电子性质和磁性奠定了坚实的计算基础。2.4计算软件与参数设置本研究采用维也纳从头算模拟软件包（ViennaAbinitioSimulationPackage，VASP）进行第一性原理计算。VASP是一款基于密度泛函理论的计算软件，它在材料科学、凝聚态物理等领域有着广泛的应用。该软件采用平面波赝势方法，能够精确描述原子和电子之间的相互作用，通过求解Kohn-Sham方程得到体系的电子结构和能量等物理性质。其优势在于计算精度高、可处理的体系类型广泛，能够有效处理包括分子、晶体、表面等多种体系，为研究单层二碲化钼及其畴界缺陷体系提供了有力的工具。在计算过程中，为了准确描述电子的交换和关联能，采用广义梯度近似（GeneralizedGradientApproximation，GGA）下的Perdew-Burke-Ernzerhof（PBE）交换关联泛函。PBE泛函考虑了电子密度的梯度信息，相比于局域密度近似（LocalDensityApproximation，LDA），能够更准确地描述电子的交换关联效应，在处理具有复杂电子结构的体系时表现出更好的性能。例如，在研究过渡金属化合物的电子结构和磁性时，PBE泛函能够更准确地描述过渡金属离子的d电子之间的相互作用，从而得到与实验更为相符的结果。平面波截断能的选择对计算结果的精度和计算效率有着重要影响。经过一系列的测试计算，确定平面波截断能设置为500eV。在测试过程中，逐步增加截断能，观察体系总能量、电子密度等物理量的变化。当截断能较低时，参与计算的平面波数量较少，计算量小，但体系总能量和电子密度等物理量的计算结果与较高截断能下的结果存在较大偏差。随着截断能的增加，计算结果逐渐收敛，当截断能达到500eV时，体系总能量和电子密度等物理量在不同截断能下的变化小于一定阈值，表明此时的计算结果已经收敛，能够满足计算精度的要求。同时，500eV的截断能在保证计算精度的前提下，也能控制计算量在可接受的范围内，确保计算能够在合理的时间内完成。对于布里渊区的k点采样，采用Monkhorst-Pack方法。在处理超晶胞结构时，k点网格设置为3×3×1。k点采样的密度对计算结果的精度也有影响，较密的k点网格能够提高计算精度，但同时也会增加计算量。通过测试不同的k点网格密度，发现3×3×1的k点网格能够在平衡计算精度和计算效率方面表现出较好的性能。在该k点网格设置下，计算得到的能带结构、态密度等物理量与更密k点网格下的计算结果相比，差异在可接受范围内，同时计算时间明显缩短，能够满足本研究对计算精度和效率的需求。在结构优化过程中，设置原子受力收敛标准为0.01eV/Å，能量收敛标准为10⁻⁵eV。当原子受力小于0.01eV/Å且能量收敛达到10⁻⁵eV时，认为体系达到了稳定状态。这样的收敛标准能够确保优化后的结构处于能量最低的稳定状态，为后续的电子性质和磁性研究提供准确的结构基础。通过对单层二碲化钼及其畴界缺陷体系的结构优化，得到了稳定的原子结构，为深入研究其电子性质和磁性奠定了坚实的基础。三、单层二碲化钼的结构与性质3.1单层二碲化钼的晶体结构二碲化钼（MoTe₂）作为一种典型的过渡金属二硫族化合物，具有独特的晶体结构。在体相结构中，MoTe₂存在两种常见的相，分别为2H相和1T'相。2H相属于六方晶系，空间群为P6₃/mmc，其晶体结构由Mo原子层和Te原子层交替堆叠而成，通过较弱的范德华力相互作用结合在一起。在这种结构中，每个Mo原子被六个Te原子以三棱柱的方式配位，形成了一个稳定的结构单元。1T'相则属于单斜晶系，空间群为C2/m，它与2H相的区别在于Mo原子的配位环境发生了变化，从三棱柱配位转变为八面体配位。在1T'相中，Mo原子的电子结构和化学键性质也相应发生改变，导致其具有与2H相不同的物理性质，如电学、光学和磁性等。当MoTe₂的厚度减薄至单层时，其原子排列方式和结构特征展现出更为独特的性质。单层MoTe₂具有类似于蜂巢状的结构，每个Mo原子位于六边形的中心，周围被六个Te原子环绕，形成了一个高度对称的平面结构。这种结构中，Mo-Te键主要表现为共价键特征，Mo原子通过与周围Te原子共享电子，形成了稳定的化学键，从而维持了单层MoTe₂的平面结构。通过第一性原理计算，我们可以得到单层MoTe₂的晶格常数等结构参数。计算结果表明，其晶格常数a和b的值相近，约为0.345nm，这与实验测量结果以及其他理论计算结果相符，进一步验证了计算方法和模型的准确性。在单层MoTe₂中，Mo-Te键的键长和键角对其结构稳定性和物理性质有着重要影响。通过结构优化计算，得到Mo-Te键的平均键长约为0.265nm。这一键长反映了Mo和Te原子之间的相互作用强度，适中的键长使得Mo和Te原子能够有效地共享电子，形成稳定的共价键。同时，MoTe₆八面体的键角也具有特定的值，例如Te-Mo-Te键角约为90°和180°，这种键角的分布决定了MoTe₂的晶体结构对称性，进而影响其电子结构和物理性质。在研究电子性质时，键角的大小会影响电子云的分布和电子的跃迁概率，从而对材料的电学和光学性质产生影响。层内原子间的相互作用除了共价键外，还存在着一定程度的离子键成分。虽然MoTe₂主要以共价键结合，但由于Mo和Te原子的电负性差异，使得电子云在Mo和Te原子之间存在一定程度的偏移，从而产生了部分离子键特性。这种离子键成分对MoTe₂的物理性质有着重要影响，在光学性质方面，离子键成分会影响材料的极化率和光吸收特性；在电学性质方面，它会影响载流子的迁移率和散射机制。通过计算差分电荷密度，可以直观地观察到电子云在Mo和Te原子之间的分布情况，进一步分析离子键成分对材料性质的影响。从差分电荷密度图中可以看出，在Mo-Te键附近，电子云出现了明显的偏移，靠近Te原子一侧的电子云密度相对较高，这表明Te原子对电子的吸引能力较强，从而体现了离子键的特性。3.2单层二碲化钼的电子结构通过第一性原理计算，我们对单层二碲化钼的电子结构进行了深入研究，包括能带结构、态密度和电荷分布，以揭示其内在的电子特性。图1展示了单层二碲化钼的能带结构。从图中可以清晰地看到，在高对称点处，导带和价带之间存在明显的带隙，这表明单层二碲化钼是一种半导体材料。计算得到的带隙值约为1.10eV，这与实验测量结果以及其他理论计算结果相符。在能带结构中，导带底位于Γ点，价带顶则位于K点，这种导带底和价带顶在动量空间中的位置关系决定了其电子跃迁特性。在Γ点附近，导带和价带的色散关系较为平缓，这意味着电子在这些区域的有效质量较大，运动相对较为缓慢；而在K点附近，能带的色散关系较为陡峭，电子的有效质量较小，具有较高的迁移率。这种能带结构特征对单层二碲化钼的电学性质有着重要影响，在电子输运过程中，电子在不同区域的迁移率差异会导致电流的不均匀分布，进而影响器件的性能。为了进一步分析电子结构，我们计算了单层二碲化钼的总态密度（TDOS）和分波态密度（PDOS），结果如图2所示。总态密度反映了体系中所有电子态在能量空间的分布情况。从图中可以看出，在费米能级附近，态密度存在明显的峰值，这表明在该能量区域内存在较多的电子态，电子的活动较为活跃。分波态密度则展示了不同原子轨道对总态密度的贡献。Mo原子的4d轨道和Te原子的5p轨道在费米能级附近的态密度贡献较大，这说明在单层二碲化钼中，Mo-Te键主要是由Mo原子的4d轨道和Te原子的5p轨道相互作用形成的。在价带区域，Te原子的5p轨道贡献占据主导地位，这表明价带中的电子主要来源于Te原子的5p电子；而在导带区域，Mo原子的4d轨道贡献相对较大，说明导带中的电子主要与Mo原子的4d电子相关。这种原子轨道的贡献分布，不仅决定了Mo-Te键的性质，还对材料的光学、电学等性质产生重要影响，在光学吸收过程中，电子从价带跃迁到导带，涉及到这些原子轨道之间的电子转移，因此原子轨道的贡献分布会影响光吸收的波长和强度。通过计算差分电荷密度，我们研究了单层二碲化钼中电荷的分布情况，结果如图3所示。差分电荷密度图能够直观地展示电子在原子间的转移和分布情况。从图中可以看出，在Mo-Te键附近，电子云出现了明显的偏移，靠近Te原子一侧的电子云密度相对较高。这表明在Mo-Te键中，电子更倾向于分布在Te原子周围，这是由于Te原子的电负性比Mo原子大，对电子的吸引能力更强，从而导致电子云向Te原子偏移。这种电荷分布的不均匀性，使得Mo-Te键具有一定的极性，进而影响材料的电学性质，如载流子的迁移率和散射机制。由于电荷分布的不均匀，载流子在运动过程中会受到额外的散射作用，从而降低载流子的迁移率。同时，这种电荷分布的不均匀性也会影响材料的光学性质，在非线性光学过程中，电荷分布的不均匀会导致材料的极化率发生变化，从而影响非线性光学响应。3.3单层二碲化钼的电学性质单层二碲化钼的电学性质与其电子结构密切相关，载流子迁移率和电导率是衡量其电学性能的重要指标。载流子迁移率反映了载流子在材料中运动的难易程度，它与材料的能带结构、杂质和缺陷等因素密切相关。在单层二碲化钼中，导带和价带的色散关系对载流子迁移率有着重要影响。在导带底和价带顶附近，能带的色散关系较为陡峭，这意味着电子和空穴在这些区域具有较小的有效质量，能够相对自由地运动，从而具有较高的迁移率。理论计算表明，室温下单层二碲化钼的电子迁移率可达100-200cm²/(V・s)，这一数值相对较高，使其在高速电子学器件应用中具有潜在优势。例如，在晶体管应用中，较高的载流子迁移率可提高器件的开关速度和工作频率，有助于实现高性能的集成电路。然而，实际制备的单层二碲化钼中往往存在各种本征缺陷，如点缺陷和畴界缺陷等，这些缺陷会对载流子迁移率产生显著影响。点缺陷会破坏晶体的周期性结构，导致电子散射增强，从而降低载流子迁移率。畴界缺陷作为一种线缺陷，其原子排列的不规则性会在局部区域产生额外的散射中心，进一步阻碍载流子的运动。通过第一性原理计算研究畴界缺陷对载流子迁移率的影响，发现畴界缺陷附近的原子位移和键长变化会导致电子波函数的畸变，增加电子与缺陷的相互作用，使得载流子迁移率明显下降。这一结果表明，在制备高质量单层二碲化钼材料时，需要严格控制畴界缺陷的密度和分布，以提高材料的电学性能。电导率是描述材料导电能力的物理量，它与载流子浓度和迁移率密切相关，其计算公式为\sigma=nq\mu，其中\sigma为电导率，n为载流子浓度，q为载流子电荷量，\mu为载流子迁移率。在本征单层二碲化钼中，载流子浓度主要由热激发产生的电子-空穴对决定。由于其带隙相对较小，在一定温度下会有一定数量的电子从价带激发到导带，形成导电载流子。然而，与金属相比，本征半导体的载流子浓度较低，导致其电导率相对较小。通过对单层二碲化钼电导率的计算，得到其本征电导率在10⁻⁵-10⁻³S/cm的量级，这表明本征单层二碲化钼是一种中等导电性的材料。在实际应用中，常常通过掺杂等手段来调控单层二碲化钼的电学性质。掺杂可以引入额外的载流子，从而显著提高载流子浓度，进而提高电导率。当在单层二碲化钼中掺入施主杂质时，杂质原子会向导带提供额外的电子，使导带中的电子浓度增加，从而形成n型半导体，电导率相应增大；当掺入受主杂质时，杂质原子会在价带中产生空穴，增加空穴浓度，形成p型半导体，同样会提高电导率。通过精确控制掺杂类型和浓度，可以实现对单层二碲化钼电学性质的有效调控，满足不同器件应用的需求。在制备场效应晶体管时，可以通过适当的掺杂来优化器件的电学性能，提高其开关比和导通电流。3.4单层二碲化钼的光学性质本研究采用基于密度泛函理论的第一性原理方法，深入研究了单层二碲化钼的光学性质。通过计算其介电函数，进而导出吸收系数、反射率等光学参数，全面分析了单层二碲化钼的光学特性。介电函数是描述材料光学性质的重要物理量，它反映了材料对电磁波的响应能力，与材料的电子结构密切相关。在固体材料中，介电函数的虚部\varepsilon_{2}(\omega)可通过计算材料的带间跃迁概率得到，其表达式为：\varepsilon_{2}(\omega)=\frac{4\pie^{2}}{\Omega\omega^{2}}\sum_{n,n',k}\vert\langle\psi_{n,k}\vert\vec{r}\vert\psi_{n',k}\rangle\vert^{2}\delta(E_{n',k}-E_{n,k}-\hbar\omega)其中，e为电子电荷，\Omega为晶体体积，\omega为光频率，\psi_{n,k}和\psi_{n',k}分别为初始态和末态的电子波函数，E_{n,k}和E_{n',k}为对应的能量，\langle\psi_{n,k}\vert\vec{r}\vert\psi_{n',k}\rangle为偶极矩阵元，\delta函数用于保证能量守恒，即只有当末态能量与初始态能量之差等于光子能量\hbar\omega时，带间跃迁才会发生。介电函数的实部\varepsilon_{1}(\omega)则可通过Kramers-Kronig关系从虚部计算得到：\varepsilon_{1}(\omega)=1+\frac{2}{\pi}P\int_{0}^{\infty}\frac{\omega'\varepsilon_{2}(\omega')}{\omega'^{2}-\omega^{2}}d\omega'其中P表示积分的主值。通过上述公式，我们可以从理论上计算出单层二碲化钼的介电函数。图4展示了单层二碲化钼的介电函数虚部\varepsilon_{2}(\omega)随光子能量的变化曲线。从图中可以看出，在低光子能量区域，介电函数虚部存在多个明显的峰值。这些峰值对应着不同的带间跃迁过程，反映了单层二碲化钼的电子结构特征。通过对能带结构和态密度的分析可知，较低能量的峰值主要源于价带顶附近的电子跃迁到导带底附近的能级，这些跃迁过程涉及到Mo原子的4d轨道和Te原子的5p轨道之间的电子转移。随着光子能量的增加，其他高能级之间的跃迁也逐渐出现，导致介电函数虚部的变化。根据介电函数，我们进一步计算了单层二碲化钼的吸收系数\alpha(\omega)，其计算公式为：\alpha(\omega)=\frac{\omega}{c}\sqrt{\frac{\varepsilon_{1}^{2}(\omega)+\varepsilon_{2}^{2}(\omega)-\varepsilon_{1}(\omega)}{2}}其中c为真空中的光速。图5给出了吸收系数随光子能量的变化曲线。从图中可以看出，单层二碲化钼在可见光和近红外波段具有较高的吸收系数。在1.0-2.5eV的光子能量范围内，吸收系数可达10^{5}cm^{-1}量级，这表明单层二碲化钼对该波段的光具有较强的吸收能力。这种良好的光吸收特性使得单层二碲化钼在光电探测器、光催化等领域具有潜在的应用价值。在光电探测器中，高吸收系数意味着能够更有效地吸收光子，产生更多的光生载流子，从而提高探测器的响应灵敏度。在发光特性方面，单层二碲化钼的直接带隙特性使其在光致发光过程中具有较高的发光效率。当受到能量大于其带隙的光激发时，价带中的电子被激发到导带，形成电子-空穴对。这些电子-空穴对在复合过程中会以光子的形式释放能量，产生光致发光现象。理论计算和实验研究表明，单层二碲化钼的光致发光峰位于近红外波段，与吸收光谱的特征相呼应。这种近红外发光特性在光通信、生物成像等领域具有潜在的应用前景。在光通信中，近红外波段的发光可以用于光信号的传输和调制，实现高速、低损耗的光通信；在生物成像中，近红外光具有较好的穿透性，能够对生物组织进行深层成像，有助于疾病的早期诊断和治疗。单层二碲化钼的光学性质使其在光电器件领域展现出巨大的应用潜力。在光电探测器方面，其高吸收系数和直接带隙特性使其能够高效地吸收光子并产生光生载流子，有望实现对可见光和近红外光的高灵敏度探测。通过与其他材料结合，构建异质结结构，可以进一步优化光电探测器的性能，提高其响应速度和探测精度。在发光二极管（LED）应用中，单层二碲化钼的近红外发光特性可用于开发新型的近红外LED，满足光通信、生物医学等领域对近红外光源的需求。通过对材料进行适当的掺杂和结构设计，可以调控其发光波长和发光效率，实现高性能的近红外LED器件。四、单层二碲化钼畴界缺陷模型构建4.1畴界缺陷的形成机制在二维材料的制备过程中，畴界缺陷的形成是一个复杂的过程，受到多种因素的影响。在单层二碲化钼的制备中，常见的化学气相沉积（CVD）和化学气相传输（CVT）等方法，由于生长过程中的原子扩散、吸附和反应等微观过程的随机性，不可避免地会引入畴界缺陷。以化学气相沉积法为例，在生长初期，衬底表面会形成多个二维晶核。这些晶核的取向往往是随机的，随着生长的进行，不同取向的晶核逐渐长大并相互融合。当两个具有不同晶格取向的晶粒相遇时，由于晶格的不匹配，在它们的交界处就会形成畴界缺陷。在化学气相沉积生长单层二碲化钼时，衬底表面的原子排列方式与二碲化钼的晶格结构存在差异，这会导致在成核过程中，不同区域的二碲化钼晶核在生长方向上出现偏差。当这些晶核长大后相互连接，就会在连接区域形成畴界，畴界处的原子需要通过调整自身位置和键长来适应不同晶格取向的过渡，从而产生了畴界缺陷。化学气相传输法制备单层二碲化钼时，反应气体在高温下分解产生的钼原子和碲原子在衬底表面沉积并反应生成二碲化钼。在这个过程中，反应气体的浓度分布、衬底表面的温度均匀性以及原子的扩散速率等因素都会影响二碲化钼的生长。如果衬底表面温度不均匀，在温度较高的区域，原子的扩散速率较快，生长速率也较快；而在温度较低的区域，生长速率则较慢。这种生长速率的差异会导致不同区域生长的二碲化钼晶体在取向和晶格参数上出现差异，当它们相互融合时，就容易形成畴界缺陷。衬底的性质对畴界缺陷的形成也有着重要影响。衬底的表面粗糙度、晶格匹配度以及化学活性等因素都会影响二碲化钼的成核和生长过程，进而影响畴界缺陷的形成。当衬底表面粗糙度较大时，会增加二碲化钼成核的随机性，导致不同取向的晶核数量增多，从而增加了畴界缺陷形成的概率。如果衬底与二碲化钼的晶格匹配度较差，在生长过程中会产生较大的晶格应力，为了释放这种应力，晶体在生长过程中会出现晶格畸变和取向调整，这也容易导致畴界缺陷的形成。生长温度和气体流量等生长条件对畴界缺陷的形成同样具有显著影响。生长温度决定了原子的扩散速率和化学反应速率。在较低的生长温度下，原子的扩散速率较慢，可能会导致原子在沉积过程中无法及时找到合适的晶格位置，从而形成较多的缺陷；而在过高的生长温度下，原子的扩散过于剧烈，会增加晶核生长的随机性，也容易导致畴界缺陷的产生。气体流量则影响着反应气体在衬底表面的浓度分布和反应速率。如果气体流量不均匀，会导致衬底表面不同区域的反应速率不一致，进而影响二碲化钼的生长均匀性，增加畴界缺陷形成的可能性。4.2畴界缺陷模型的建立为了深入研究单层二碲化钼中畴界缺陷对其电子性质和磁性的影响，我们采用基于密度泛函理论的第一性原理方法，构建了一系列畴界缺陷模型。首先，考虑到单层二碲化钼的六边形晶格结构，我们构建了具有不同取向的畴界缺陷模型。在实际生长过程中，由于原子的随机排列和生长条件的不均匀性，畴界的取向会呈现出多样性。通过构建不同取向的畴界模型，能够更全面地研究畴界取向对材料性质的影响。对于晶界错配位错型畴界，我们通过调整相邻晶粒的相对位置和取向，使得两个晶粒在交界处出现原子的错配和位错。具体来说，选取两个具有不同晶格取向的单层二碲化钼原胞，将它们沿着特定的方向拼接在一起，形成畴界。在拼接过程中，保持原子总数不变，通过结构优化使得体系能量达到最低。图6展示了一种典型的晶界错配位错型畴界模型的原子结构，从图中可以清晰地看到畴界处原子的排列不规则性，Mo原子和Te原子的位置发生了明显的偏移，键长和键角也发生了变化。对于堆垛层错型畴界，我们模拟了不同原子层之间的堆垛顺序错误。在理想的单层二碲化钼中，原子层具有特定的堆垛顺序，而堆垛层错会破坏这种顺序，导致畴界的形成。通过人为地改变原子层的堆垛顺序，构建堆垛层错型畴界模型。在结构优化过程中，采用共轭梯度法对原子位置进行弛豫，使体系达到能量最低的稳定状态。在优化过程中，不断调整原子的位置，使得原子间的相互作用力达到平衡，从而得到稳定的畴界结构。图7给出了堆垛层错型畴界模型的原子结构示意图，从图中可以观察到堆垛层错处原子的排列混乱，这会对电子的传输和磁性产生重要影响。为了研究畴界缺陷浓度对材料性质的影响，我们构建了不同缺陷间距的畴界缺陷模型。通过改变畴界之间的距离，来模拟不同浓度的畴界缺陷。在构建过程中，保持畴界的类型和取向不变，逐步减小畴界之间的距离，增加缺陷浓度。图8展示了不同缺陷间距的畴界缺陷模型，随着缺陷间距的减小，畴界缺陷的浓度逐渐增加，畴界之间的相互作用也逐渐增强。这种相互作用会导致电子态的进一步变化，影响材料的电学和磁学性质。当缺陷间距较小时，畴界之间的电子云会发生重叠，使得电子的散射增强，载流子迁移率降低；同时，畴界之间的磁性相互作用也会增强，可能导致材料的磁性发生变化。通过构建这些不同类型和浓度的畴界缺陷模型，为后续研究畴界缺陷对单层二碲化钼电子性质和磁性的影响提供了基础。在建立模型的过程中，我们充分考虑了实际生长过程中可能出现的各种情况，使得模型更接近真实的材料体系。这些模型的建立，有助于深入理解畴界缺陷在单层二碲化钼中的形成机制和作用规律，为优化材料性能和开发新型器件提供理论指导。4.3模型的优化与验证在构建了单层二碲化钼畴界缺陷模型后，对这些模型进行结构优化是确保模型准确性和可靠性的关键步骤。结构优化的目的是使体系达到能量最低的稳定状态，从而得到最合理的原子位置和结构参数。在优化过程中，采用共轭梯度法对原子位置进行弛豫。共轭梯度法是一种常用的优化算法，它通过迭代计算，不断调整原子的位置，使得体系的能量逐渐降低，直到满足收敛条件。在VASP软件中，设置原子受力收敛标准为0.01eV/Å，能量收敛标准为10⁻⁵eV。当原子所受的力小于0.01eV/Å且体系能量收敛达到10⁻⁵eV时，认为体系达到了稳定状态。通过这种严格的收敛标准，能够确保优化后的结构处于能量最低的稳定状态，为后续的电子性质和磁性研究提供准确的结构基础。以晶界错配位错型畴界模型为例，在优化前，畴界处原子的位置由于晶格的不匹配而处于较高能量状态，原子间的相互作用力不平衡。经过结构优化后，畴界处的原子通过调整位置和键长，使得原子间的相互作用力达到平衡，体系能量降低。从优化前后的原子结构对比图（图9）中可以清晰地看到，优化后畴界处的原子排列更加有序，键长和键角也发生了合理的变化，更符合实际的物理情况。为了验证所构建模型的准确性，将计算得到的结构参数与实验结果或其他理论计算结果进行对比。在晶格常数方面，实验测得的单层二碲化钼的晶格常数a和b约为0.345nm，通过第一性原理计算得到的晶格常数与实验值相符，误差在可接受范围内，这表明所采用的计算方法和模型能够准确地描述单层二碲化钼的晶体结构。对于畴界缺陷处的原子结构和键长等参数，虽然目前实验上直接测量存在一定困难，但与已有的理论研究结果进行对比，发现计算结果具有一致性。一些理论研究通过高分辨率透射电子显微镜模拟等方法，得到了畴界缺陷处的原子结构信息，本研究中优化后的模型与这些理论模拟结果在原子排列和键长变化趋势上相吻合，进一步验证了模型的可靠性。通过对单层二碲化钼畴界缺陷模型进行结构优化，并与实验和其他理论计算结果进行对比验证，确保了模型的准确性和可靠性，为后续深入研究畴界缺陷对单层二碲化钼电子性质和磁性的影响奠定了坚实的基础。这些优化后的模型能够更真实地反映实际材料中畴界缺陷的结构特征，有助于揭示畴界缺陷与材料性能之间的内在联系，为二维材料的性能优化和器件应用提供有力的理论支持。五、单层二碲化钼畴界缺陷的电子性质5.1缺陷对能带结构的影响为了深入探究畴界缺陷对单层二碲化钼电子性质的影响，我们首先对完美的单层二碲化钼以及含有畴界缺陷的单层二碲化钼进行了能带结构计算，计算结果如图10所示。图10(a)展示了完美单层二碲化钼的能带结构。在完美结构中，能带结构呈现出典型的半导体特征，导带和价带之间存在明显的带隙，计算得到的带隙值约为1.10eV，这与之前的研究结果一致。在高对称点处，导带底位于Γ点，价带顶位于K点，这种导带底和价带顶在动量空间中的位置关系决定了其电子跃迁特性。在Γ点附近，导带和价带的色散关系较为平缓，这意味着电子在这些区域的有效质量较大，运动相对较为缓慢；而在K点附近，能带的色散关系较为陡峭，电子的有效质量较小，具有较高的迁移率。当引入畴界缺陷后，能带结构发生了显著变化，如图10(b)所示。畴界缺陷的存在破坏了晶体的周期性结构，导致能带出现明显的畸变。在畴界缺陷附近，导带和价带的能量发生了偏移，带隙也发生了改变。具体来说，缺陷处的原子排列不规则，使得电子的波函数受到干扰，电子态发生变化。在某些情况下，缺陷会导致杂质能级的出现，这些杂质能级位于带隙中，对电子的输运和光学性质产生重要影响。从图中可以看出，在价带顶和导带底之间出现了一些离散的能级，这些能级就是由于畴界缺陷产生的杂质能级。这些杂质能级的存在，为电子的跃迁提供了额外的通道，可能会改变材料的电学和光学性质。为了进一步分析畴界缺陷对能带结构的影响机制，我们对缺陷处的原子结构和电子云分布进行了研究。通过差分电荷密度分析发现，畴界缺陷处的原子发生了明显的位移，键长和键角也发生了变化。这些结构变化导致了电子云的重新分布，使得电子与缺陷之间的相互作用增强。在缺陷处，电子云出现了局域化现象，电子的运动受到限制，从而影响了能带结构。由于缺陷处的原子位移，使得Mo-Te键的键长发生变化，导致电子云在Mo和Te原子之间的分布发生改变，进而影响了电子的能量状态，使得能带结构发生畸变。畴界缺陷的类型和浓度对能带结构的影响也存在差异。对于不同类型的畴界缺陷，如晶界错配位错型畴界和堆垛层错型畴界，它们对能带结构的影响方式和程度有所不同。晶界错配位错型畴界主要通过破坏原子的周期性排列，导致电子散射增强，从而使能带发生畸变；而堆垛层错型畴界则主要通过改变原子层之间的堆垛顺序，影响电子的层间耦合，进而改变能带结构。随着畴界缺陷浓度的增加，能带的畸变程度也会加剧，杂质能级的数量和分布也会发生变化。当缺陷浓度较低时，杂质能级相对较少，且分布较为离散；当缺陷浓度增加时，杂质能级的数量增多，它们之间的相互作用也会增强，可能会形成杂质能带，进一步影响材料的电学性质。5.2缺陷对态密度的影响态密度（DensityofStates，DOS）能够直观地反映体系中电子态在能量空间的分布情况，对于研究材料的电子性质具有重要意义。通过计算完美单层二碲化钼以及含有畴界缺陷的单层二碲化钼的态密度，我们深入分析了畴界缺陷对电子态分布的影响。图11展示了完美单层二碲化钼和含有畴界缺陷的单层二碲化钼的总态密度（TDOS）和分波态密度（PDOS）。在完美单层二碲化钼的态密度图中，总态密度在费米能级附近呈现出明显的特征。在价带区域，态密度逐渐增大，表明在该能量范围内电子态的分布较为密集；而在导带区域，态密度相对较小，这与能带结构中价带和导带的特征相对应。从分波态密度来看，Mo原子的4d轨道和Te原子的5p轨道在费米能级附近的态密度贡献较大，这说明Mo-Te键主要是由Mo原子的4d轨道和Te原子的5p轨道相互作用形成的。在价带中，Te原子的5p轨道贡献占据主导地位，这表明价带中的电子主要来源于Te原子的5p电子；而在导带中，Mo原子的4d轨道贡献相对较大，说明导带中的电子主要与Mo原子的4d电子相关。当引入畴界缺陷后，态密度发生了显著变化。在总态密度中，费米能级附近的态密度出现了明显的峰和谷，这是由于畴界缺陷导致了电子态的重新分布。在缺陷处，原子排列的不规则性使得电子的波函数受到干扰，电子态发生改变，从而在态密度中表现出明显的变化。从分波态密度来看，Mo原子和Te原子的态密度分布也发生了变化。在缺陷附近，Mo原子的4d轨道和Te原子的5p轨道的态密度出现了新的峰和谷，这表明缺陷对Mo-Te键的电子结构产生了显著影响。Mo原子4d轨道的态密度在某些能量区域出现了增强，这可能是由于缺陷导致了Mo原子周围电子云的重新分布，使得4d轨道的电子态发生了变化。为了进一步分析缺陷态在态密度中的表现，我们关注了缺陷处产生的杂质能级。在含有畴界缺陷的单层二碲化钼的态密度图中，可以明显观察到在带隙中出现了一些离散的态密度峰，这些峰对应着杂质能级。这些杂质能级的存在，为电子的占据和跃迁提供了额外的可能性。在电子占据态方面，杂质能级可能被电子占据，从而改变了材料的电子填充情况。当杂质能级被电子占据时，会导致材料的电学性质发生变化，可能会增加材料的导电性，也可能会改变材料的载流子类型。在未占据态方面，杂质能级的存在为电子的激发提供了新的途径。当材料受到光激发或其他外界作用时，电子可以从价带跃迁到杂质能级，或者从杂质能级跃迁到导带，这会影响材料的光学性质，如光吸收和发光特性。畴界缺陷对电子占据态和未占据态的影响还与缺陷的类型和浓度密切相关。对于不同类型的畴界缺陷，如晶界错配位错型畴界和堆垛层错型畴界，它们在态密度中产生的杂质能级的位置和强度可能不同。晶界错配位错型畴界由于其原子排列的不规则性，可能会在带隙中产生较多且分布较为离散的杂质能级；而堆垛层错型畴界则可能由于层间堆垛顺序的改变，导致杂质能级的分布相对集中在某些特定的能量区域。随着畴界缺陷浓度的增加，杂质能级的数量也会相应增加，它们之间的相互作用也会增强。当杂质能级的数量较多且相互作用较强时，可能会形成杂质能带，进一步影响电子的占据和跃迁行为，从而对材料的电学和光学性质产生更为显著的影响。5.3缺陷处的电荷分布与局域态通过计算差分电荷密度，我们深入研究了单层二碲化钼畴界缺陷处的电荷重新分布情况。差分电荷密度能够直观地展示在形成畴界缺陷后，电子在原子间的转移和分布变化。图12展示了含有畴界缺陷的单层二碲化钼的差分电荷密度图。在畴界缺陷处，可以明显观察到电荷的重新分布。与完美的单层二碲化钼相比，缺陷处的电子云发生了显著变化。在畴界附近，电子云出现了局域化现象，电子更倾向于聚集在某些特定的原子周围。这是由于畴界处的原子排列不规则，导致原子间的化学键发生变化，电子的运动受到限制，从而使得电子云在这些区域聚集。Mo原子和Te原子之间的键长和键角在畴界处发生了改变，这影响了电子云在它们之间的分布。一些Te原子周围的电子云密度明显增加，这表明电子在这些位置出现了局域化，形成了相对稳定的电子态。这种电荷重新分布对材料的电学性质产生了重要影响。畴界缺陷处的电荷局域化会导致载流子的散射增强。由于缺陷处电子云的局域化，载流子在通过畴界区域时，会与局域化的电子云相互作用，从而增加了散射概率，降低了载流子迁移率。畴界缺陷处的电荷分布不均匀还会导致局部电场的产生。这种局部电场会对载流子的运动产生额外的作用力，进一步影响载流子的输运过程，使得材料的电学性能变得更加复杂。局域态的形成与电荷重新分布密切相关。在畴界缺陷处，由于原子排列的不规则性和电荷的局域化，形成了一些局域态。这些局域态的电子被限制在特定的原子周围，具有相对稳定的能量状态。通过分析电子局域函数（ELF），我们可以更直观地了解局域态的分布情况。电子局域函数能够描述电子在空间中的局域化程度，其值在0-1之间，值越接近1，表示电子的局域化程度越高。在畴界缺陷处，电子局域函数的值明显增大，表明电子在这些区域的局域化程度较高，形成了明显的局域态。局域态的性质对材料的物理性质有着重要影响。在光学性质方面，局域态可以作为电子跃迁的中间态，影响材料的光吸收和发光特性。当材料受到光激发时，电子可以从价带跃迁到局域态，或者从局域态跃迁到导带，这会导致光吸收谱和发光谱中出现新的峰。在磁性方面，局域态中的电子具有一定的自旋特性，可能会对材料的磁性产生贡献。如果局域态中的电子自旋取向一致，就可能会产生局部的磁矩，从而影响材料的整体磁性。畴界缺陷处的电荷重新分布导致了局域态的形成，这些局域态对材料的电学、光学和磁性等性质产生了重要影响。深入研究缺陷处的电荷分布和局域态，有助于揭示畴界缺陷对单层二碲化钼物理性质的影响机制，为二维材料的性能优化和器件应用提供理论支持。5.4电子性质与材料性能的关联畴界缺陷引起的电子性质变化对单层二碲化钼的电学、光学和催化性能产生了深远影响，这些影响与材料在实际应用中的性能密切相关。在电学性能方面，畴界缺陷导致的能带畸变和杂质能级的出现，显著影响了材料的电导率和载流子迁移率。畴界缺陷处的原子排列不规则，破坏了晶体的周期性结构，使得电子在传输过程中受到散射增强，从而降低了载流子迁移率。杂质能级的存在也会改变材料的载流子浓度和类型，进而影响电导率。当杂质能级被电子占据时，会增加材料的电子浓度，可能使材料从本征半导体转变为n型半导体；反之，若杂质能级为空穴占据，则可能使材料转变为p型半导体。这种电学性能的改变，对基于单层二碲化钼的电子器件性能有着重要影响。在晶体管应用中，载流子迁移率的降低会导致器件的开关速度变慢，工作频率降低，从而限制了器件的性能提升；而电导率的变化则会影响器件的导通电流和功耗，进而影响整个电路的性能。在光学性能方面，畴界缺陷对材料的光吸收和发光特性产生了显著影响。由于畴界缺陷导致的能带结构变化，使得材料的光吸收谱和发光谱发生改变。畴界缺陷处的杂质能级为电子的跃迁提供了额外的通道，导致光吸收峰的位置和强度发生变化。在光吸收过程中，电子可以从价带跃迁到杂质能级，或者从杂质能级跃迁到导带，从而增加了光吸收的途径，改变了材料对不同波长光的吸收能力。在发光特性方面，缺陷处的局域态可以作为电子跃迁的中间态，影响材料的发光效率和发光波长。当电子从导带跃迁到局域态，再从局域态跃迁到价带时，会发出特定波长的光，从而改变材料的发光颜色和强度。这些光学性能的变化，使得单层二碲化钼在光电器件应用中具有不同的表现。在光电探测器中，光吸收特性的改变可能会影响探测器的响应灵敏度和探测波长范围；在发光二极管中，发光特性的变化则会影响器件的发光效率和颜色纯度。在催化性能方面，畴界缺陷的存在为材料提供了更多的活性位点，从而增强了材料的催化活性。畴界缺陷处的原子排列不规则，导致原子的配位不饱和，使得这些原子具有较高的化学活性。在催化反应中，反应物分子更容易吸附在这些活性位点上，发生化学反应。在一些氧化还原反应中，畴界缺陷处的Mo原子或Te原子可以作为活性中心，促进反应物分子的吸附和活化，从而加快反应速率。通过计算吸附能和反应能垒，可以进一步了解畴界缺陷对催化反应的影响机制。吸附能的大小反映了反应物分子与材料表面的结合强度，而反应能垒则决定了反应的难易程度。研究发现，在含有畴界缺陷的单层二碲化钼表面，反应物分子的吸附能明显降低，反应能垒也相应减小，这表明畴界缺陷能够显著增强材料的催化活性。六、单层二碲化钼畴界缺陷的磁性6.1磁性的起源与理论分析磁性是材料的一种重要物理性质，其起源与材料的电子结构密切相关。在单层二碲化钼中，本征的MoTe₂材料通常表现为非磁性，但当引入畴界缺陷后，部分缺陷处会产生磁性，这为其在自旋电子学领域的应用提供了新的可能性。从理论上来说，磁性的产生主要源于电子的自旋和轨道运动。在原子中，电子的自旋和轨道角动量会产生磁矩，当这些磁矩在材料中有序排列时，就会表现出宏观的磁性。在单层二碲化钼中，畴界缺陷的存在破坏了晶体的周期性结构，导致电子云的分布发生变化，从而使得部分电子的自旋和轨道角动量发生改变，产生了局域磁矩。畴界处的原子位移和键长变化会导致电子与原子核之间的相互作用发生改变，进而影响电子的自旋和轨道运动，使得电子的磁矩不再相互抵消，产生了净磁矩。在单层二碲化钼中，Mo原子的4d轨道和Te原子的5p轨道对磁性的贡献较为显著。畴界缺陷处，Mo原子和Te原子的配位环境发生变化，使得这些轨道上的电子态发生改变。由于原子位移，Mo原子的4d轨道与周围Te原子的5p轨道之间的相互作用增强，导致电子云的分布发生变化，部分电子的自旋取向发生改变，从而产生了局域磁矩。通过计算原子磁矩，我们发现畴界缺陷处的Mo原子和Te原子的磁矩明显不同于完美结构中的原子磁矩，这进一步证实了缺陷对原子磁矩的影响。为了更深入地理解磁性的起源，我们考虑了电子的交换关联效应。在密度泛函理论中，交换关联泛函描述了电子之间的交换和关联作用。在畴界缺陷处，由于电子云的局域化，电子之间的交换关联效应增强，这对磁性的产生和稳定起到了重要作用。通过选择合适的交换关联泛函，如广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函，能够更准确地描述电子之间的相互作用，从而更好地解释磁性的起源。从晶体场理论的角度来看，畴界缺陷处的原子结构变化会导致晶体场的畸变。晶体场的变化会影响电子在不同轨道上的能量分布，从而改变电子的自旋状态。在单层二碲化钼中，畴界缺陷处的Mo原子周围的晶体场发生畸变，使得Mo原子的4d轨道分裂，电子在这些轨道上的填充方式发生改变，进而导致磁矩的产生。这种基于晶体场理论的分析，为理解磁性的起源提供了另一个重要的视角，与基于电子结构和交换关联效应的分析相互补充，共同揭示了单层二碲化钼畴界缺陷处磁性的产生机制。6.2不同畴界缺陷的磁性特征为了深入了解单层二碲化钼中畴界缺陷的磁性特征，我们对不同类型的畴界缺陷进行了研究，包括晶界错配位错型畴界和堆垛层错型畴界，并分析了不同缺陷浓度下的磁性变化。对于晶界错配位错型畴界，计算结果表明，畴界缺陷处会产生明显的局域磁矩。图13展示了晶界错配位错型畴界处的原子磁矩分布情况。从图中可以看出，在畴界附近，Mo原子和Te原子的磁矩发生了显著变化。部分Mo原子的磁矩可达0.5μB以上，而Te原子的磁矩也有一定程度的增加。这是由于畴界处原子排列的不规则性，导致电子云分布发生改变，使得Mo原子和Te原子的电子自旋和轨道角动量发生变化，从而产生了局域磁矩。堆垛层错型畴界的磁性特征与晶界错配位错型畴界有所不同。图14给出了堆垛层错型畴界处的原子磁矩分布。在堆垛层错处，虽然也有部分原子产生了磁矩，但磁矩的大小和分布与晶界错配位错型畴界存在差异。堆垛层错型畴界处的磁矩主要集中在层错区域的边缘，且磁矩的大小相对较小，一般在0.2-0.3μB之间。这是因为堆垛层错主要影响了原子层之间的堆垛顺序，对原子的配位环境和电子云分布的影响方式与晶界错配位错型畴界不同，从而导致磁性特征的差异。随着畴界缺陷浓度的增加，材料的总磁矩也会发生变化。通过构建不同缺陷间距的畴界缺陷模型，研究了缺陷浓度对总磁矩的影响。结果发现，当缺陷间距较大时，畴界缺陷之间的相互作用较弱，总磁矩随着缺陷浓度的增加近似呈线性增长。随着缺陷间距的减小，畴界缺陷之间的相互作用逐渐增强，总磁矩的增长趋势逐渐变缓。当缺陷间距减小到一定程度时，畴界缺陷之间的相互作用变得非常强，可能会导致磁矩的相互抵消，使得总磁矩不再增加，甚至出现减小的趋势。不同畴界缺陷的磁性耦合机制也存在差异。在晶界错配位错型畴界中，由于畴界处原子的排列较为复杂，磁矩之间的耦合方式较为多样化，既有铁磁耦合，也有反铁磁耦合。而在堆垛层错型畴界中，磁矩之间主要表现为较弱的铁磁耦合。这种磁性耦合机制的差异，进一步影响了材料的宏观磁性。在晶界错配位错型畴界中，由于磁矩耦合方式的多样性，材料的宏观磁性可能会表现出较为复杂的行为；而在堆垛层错型畴界中，由于主要是铁磁耦合，材料的宏观磁性相对较为简单，可能更容易表现出铁磁性。6.3磁性与电子结构的关系单层二碲化钼畴界缺陷处的磁性与电子结构之间存在着紧密的内在联系，这种联系深刻影响着材料的物理性质。从电子云分布的角度来看，畴界缺陷的存在导致了电子云的重新分布。畴界处原子排列的不规则性使得电子云在某些区域出现局域化现象。这种局域化的电子云分布对磁性产生了重要影响。当电子云在某个原子周围局域化时，该原子的电子自旋和轨道角动