ZnO基稀磁半导体电子结构与磁特性：基于第一性原理与蒙特卡洛模拟的探究

ZnO基稀磁半导体电子结构与磁特性：基于第一性原理与蒙特卡洛模拟的探究一、绪论1.1半导体自旋电子学概述半导体自旋电子学是一门研究半导体材料中电子自旋属性及其应用的新兴交叉学科，其发展历程充满了探索与突破。20世纪80年代，随着纳米科学技术的兴起，科学家们逐渐意识到电子自旋在电子器件中的潜在应用价值，自旋电子学的概念应运而生。1988年，巴西学者Baibich在[Fe/Cr]周期性多层膜中发现了巨磁电阻（GMR）效应，当施加外磁场时，其电阻变化率高达50%，这一发现犹如一颗启明星，为自旋电子学的发展照亮了道路，成为该领域发展的重要里程碑，引发了科研人员对电子自旋相关特性的深入研究，半导体自旋电子学也由此逐渐崭露头角。半导体自旋电子学主要包含两个重要分支领域。其一为半导体磁电子学，它致力于将磁性功能融入半导体之中，比如构建磁性半导体，或者制备半导体与磁性材料的复合体。这一领域的研究成果有望直接推动半导体器件的革新，如光绝缘体、磁传感器、非挥发性存储等器件的实现。倘若磁性或自旋能够被光或电场精准控制，那么创造出前所未有的全新器件也并非遥不可及。另一个领域是半导体量子自旋电子学，主要聚焦于自旋的量子力学特性在半导体中的实际应用。在应用领域方面，半导体自旋电子学展现出了巨大的潜力。在数据存储领域，基于半导体自旋电子学原理的磁性随机存储器（MRAM）成为研究热点。与传统的动态随机存取存储器（DRAM）和闪存相比，MRAM具有非易失性、高速读写、低功耗以及无限次擦写等显著优势。例如，EverspinTechnologies公司已经推出了基于自旋转移力矩（STT）技术的MRAM产品，在工业控制、汽车电子等领域得到了应用，能够有效提高数据存储的安全性和读写速度。在量子计算领域，半导体自旋电子学也为实现量子比特提供了新的途径。电子的自旋可以作为量子比特的候选者之一，利用半导体材料中的量子点或量子阱结构，可以精确地操控电子自旋的状态，实现量子比特的初始化、单比特门操作和多比特门操作等。荷兰代尔夫特理工大学的科研团队在基于半导体自旋量子比特的量子计算研究方面取得了重要进展，成功实现了多个量子比特的纠缠，为未来量子计算机的发展奠定了基础。在半导体激光器领域，自旋极化的注入可以显著提高激光器的性能。通过控制注入载流子的自旋方向，可以增强激光器的发光效率、降低阈值电流，并实现对激光偏振特性的精确控制。例如，在垂直腔面发射激光器（VCSEL）中，利用自旋极化的电子注入，可以实现单模、低阈值的激光发射，在光通信和光存储等领域具有广阔的应用前景。半导体自旋电子学凭借其独特的物理特性和广泛的应用前景，在现代电子技术中占据着愈发重要的地位，为未来电子器件的发展和信息技术的革新提供了强大的动力和支撑。1.2稀磁半导体的研究进展1.2.1稀磁半导体的概念与特点稀磁半导体（DilutedMagneticSemiconductors，DMS），又被称作半磁半导体，是指在III-V族、II-VI族、II-V族或IV-VI族化合物半导体中，由磁性过渡族金属离子（如Mn、Fe、Co等）或稀土金属离子部分替代非磁性阳离子所形成的一类新型半导体材料。其化学式通常可表示为A_{1-x}M_{x}B，其中A为非磁性阳离子，M为磁性离子，B为阴离子，x代表磁性离子的掺杂浓度。与传统半导体相比，稀磁半导体不仅保留了半导体的电学特性，还引入了磁性，实现了在同一材料中对电子电荷和自旋两种自由度的利用，这使其具备了许多独特的物理性质和潜在应用价值。稀磁半导体具有一些显著的特点。理想状态下，其居里温度（T_C）应高于500K，这样在常温环境中便能保持稳定的铁磁性，满足自旋电子器件在室温下工作的需求。其铁磁性与形成载流子的杂质能带的自旋分裂紧密相关联，这种关联为实现对电子自旋的精确操控提供了物理基础。再者，稀磁半导体能够进行n型和p型掺杂，如同传统半导体一样，可以通过掺杂来灵活调节其电学性能，从而满足不同器件应用的需求。此外，它还应具备高的迁移率和自旋散射长度，这有助于提高电子在材料中的输运效率，减少能量损耗，进而提升器件的性能。稀磁半导体还具有磁光效应和反常霍尔效应等独特的物理效应。磁光效应使得其在光通信、光存储等领域展现出潜在的应用价值，例如利用巨法拉第旋转效应可制作光隔离器、磁光调制器等光电器件；反常霍尔效应则为研究材料的磁性和电子输运特性提供了重要的手段，通过对反常霍尔效应的测量和分析，可以深入了解材料中载流子的自旋极化和散射机制。1.2.2国内外研究现状自稀磁半导体的概念提出以来，在国内外都引发了广泛且深入的研究。在室温铁磁性研究方面，众多科研团队投入大量精力，旨在寻找或制备出居里温度高于室温的稀磁半导体材料。(Ga,Mn)As作为III-V体系稀磁半导体的典型代表，是研究最为广泛的材料之一。然而，在这类材料中，由于(Ga^{3+},Mn^{2+})的异价掺杂同时引入了自旋和电荷，使得Mn的含量难以有效提高，进而严重制约了材料居里温度的提升。为了突破这一限制，中国科学院物理研究所/北京凝聚态物理国家研究中心靳常青领导的研究团队在2011年提出了自旋电荷解耦掺杂实现稀磁半导体机制，成功实现了电荷与自旋掺杂机制的分别调控。沿着这一创新的材料设计思路，他们发现了居里温度高达180K的(Ba,K)(Zn,Mn)_2As_2（简称BZA）新型稀磁半导体，并通过进一步的研究，将其可控居里温度刷新至230K。近期，该团队运用高压热处理技术，成功将BZA的居里温度提高到260K，接近室温。高压热处理能够突破材料固溶度的限制，稳定晶格畸变，从而提高了BZA的Mn掺杂浓度，达到了常压下Mn浓度的2倍，同时保持了[MnAs]四面体基元的规则构型。在磁性来源研究方面，目前存在多种理论模型用于解释稀磁半导体的磁性起源。其中，描述绝缘体中磁性的直接交换作用和超交换作用、载流子媒介交换作用以及描述部分氧化物中掺杂磁性的束缚磁极化子模型较为常见。传统铁磁金属之间的铁磁耦合通常用直接交换作用机制来解释，而金属氧化物、硫化物、氟族化合物以及铁氧体中的反铁磁性或亚铁磁性则采用超交换作用机制来描述。超交换作用又被称为间接交换作用，其特点是磁性原子间的相互作用通过中间的阴离子间接完成，交换积分符号取决于金属氧化物的键角和过渡金属的d电子组态。载流子媒介交换作用认为，稀磁半导体中的磁性离子与载流子之间存在强烈的自旋-自旋交换相互作用，这种相互作用在外加电场或磁场的影响下，会显著改变载流子的行为，从而产生异于半导体基质的磁性特性。束缚磁极化子模型则主要用于解释部分氧化物稀磁半导体中，由于磁性离子的掺杂形成束缚磁极化子，进而导致材料呈现出磁性的现象。在制备方法研究上，国际上不断探索和改进各种制备技术，以获得高质量的稀磁半导体材料。分子束外延（MBE）技术凭借其在原子尺度上精确控制外延膜厚、掺杂和界面平整度的优势，成为制备高质量DMS薄膜的重要手段。例如，邓加军等人利用MBE技术在GaAs上成功生长出厚度为500nm的GaAs:Mn薄膜。当Mn的质量分数为7%时，晶体质量良好，呈现闪锌矿结构，且随着温度从低到高变化，铁磁转变温度从65K上升到115K。金属有机化学气相沉积（MOCVD）法也被广泛应用于制备II-VI及III-V族的稀磁半导体，该方法能够合成组分按任意比例组成的人工合成材料，并精确控制薄膜厚度到原子级。国内在稀磁半导体研究领域也取得了一系列重要成果。在制备方法方面，中国科学院物理研究所的研究团队通过不断优化制备工艺，在新型稀磁半导体材料的制备上取得了突破性进展。在性能研究方面，众多科研机构对稀磁半导体的电学、磁学和光学性能进行了深入研究，为其实际应用提供了坚实的理论基础。1.2.3研究中存在的问题及解决方法当前稀磁半导体研究中仍面临诸多关键问题，其中居里温度偏低是制约其实际应用的核心难题之一。大部分已研究的稀磁半导体材料居里温度远低于室温，导致在常温环境下无法保持稳定的铁磁性，严重限制了其在自旋电子器件中的应用。以(Ga,Mn)As为例，由于其异价掺杂导致的Mn含量限制，使得居里温度难以提升，无法满足室温下自旋电子器件的工作要求。磁性来源的不确定性也是研究中面临的重要挑战。尽管存在多种理论模型来解释稀磁半导体的磁性起源，但不同材料体系中的磁性来源复杂多样，尚未形成统一且明确的理论解释。在一些材料中，磁性可能是由多种交换作用共同作用的结果，这使得准确确定磁性来源变得极为困难，进而影响了对材料磁性的有效调控和优化。为解决居里温度偏低的问题，科研人员从多个角度展开探索。在材料设计方面，提出了自旋电荷解耦掺杂等创新机制，通过分别调控电荷和自旋掺杂，为提高居里温度开辟了新途径。如前文所述的(Ba,K)(Zn,Mn)_2As_2新型稀磁半导体，就是基于这一机制设计并取得了居里温度提升的显著成果。在制备工艺上，采用高压热处理等技术，突破材料固溶度限制，稳定晶格畸变，从而提高磁性离子掺杂浓度，提升居里温度。针对磁性来源不确定的问题，一方面通过高精度的实验技术，如同步辐射X射线衍射、穆斯堡尔谱等，深入研究材料的微观结构和电子态，以获取更准确的磁性信息；另一方面，结合第一性原理计算等理论方法，对不同材料体系中的磁性相互作用进行模拟和分析，从理论层面揭示磁性起源，为实验研究提供指导。通过实验与理论的紧密结合，逐步明确不同材料体系中磁性的来源和作用机制，为实现对稀磁半导体磁性的精准调控奠定基础。二、理论基础与研究方法2.1第一性原理计算方法2.1.1理论基础第一性原理计算方法的理论根基是量子力学，其核心在于对体系薛定谔方程的求解。量子力学作为描述微观粒子运动规律的理论，从根本上揭示了原子之间相互作用进而构成固体的本质。在固体物理学领域，量子力学成功阐释了物质为何会呈现出半导体、导体、绝缘体等不同特性，为材料科学的发展奠定了坚实的理论基础。然而，固体是包含10^{23}数量级粒子的多粒子系统，直接应用量子理论会使物理方程变得极为复杂，难以求解。因此，在将量子理论应用于固体系统时，需要采用一系列近似和简化方法。绝热近似（Born-Oppenheimer近似）是其中关键的一步，它基于原子核质量远大于电子质量这一事实，将电子的运动和原子核的运动分离开来。在固体中，原子核的运动速度相对较慢，可近似认为电子在瞬间固定的原子核势场中运动。以金属晶体为例，原子核构成晶格骨架，电子在晶格的周期性势场中运动。通过绝热近似，多粒子系统得以简化为多电子系统，大大降低了计算的复杂度。Hartree-Fock近似则进一步将多电子问题简化为仅以单电子波函数（分子轨道）为基本变量的单粒子问题。它假设每个电子都在其他电子的平均场中独立运动，通过自洽场迭代的方式求解单电子波函数。在计算原子的电子结构时，可将每个电子的波函数看作是在原子核和其他电子的平均势场中运动的结果。不过，该近似中波函数的行列式表示导致求解需要巨大的计算量，对于分子体系尚可作为研究的出发点，但在研究固态体系时存在局限性。1964年，Hohenberg和Kohn提出了严格的密度泛函理论（DensityFunctionalTheory，DFT）。DFT建立在非均匀电子气理论基础之上，其核心思想是将系统中的电子密度作为基本变量，而非传统的电子波函数。该理论指出，体系的基态能量是电子密度的唯一泛函，通过对电子密度的变分求解，可得到体系的基态能量和电子结构。这一理论的提出，为多电子体系的研究提供了全新的视角和方法。1965年，Kohn和Sham提出了Kohn-Sham方程，将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题与单电子Kohn-Sham方程。Kohn-Sham方程的形式类似于单电子的薛定谔方程，其中包含了有效势，该有效势涵盖了电子与原子核的相互作用、电子之间的库仑相互作用以及交换关联相互作用。在实际应用DFT时，需要对电子间的交换关联作用进行近似。局域密度近似（LocalDensityApproximation，LDA）假设在空间某点的交换关联能只与该点的电子密度有关，且等于具有相同电子密度的均匀电子气的交换关联能。这种近似在处理电子密度变化缓慢的体系时表现良好，例如金属和简单半导体。广义梯度近似（GeneralizedGradientApproximation，GGA）则考虑了电子密度的梯度对交换关联能的影响，能够更准确地描述电子密度变化较为剧烈的体系，如分子和一些过渡金属氧化物。GGA通过引入与电子密度梯度相关的项，对交换关联能进行修正，提高了计算的精度。2.1.2实现过程利用第一性原理计算材料性质，首先需确定体系模型。对于晶体材料，由于其具有周期对称性，需明确晶胞和晶胞中原子的坐标。在确定模型前，通常要先确定平衡晶格常数。以半导体材料为例，可通过改变三个基矢的大小，使单胞体积在一定范围内变化（如81%-119%）。运用第一性原理计算不同体积模型的能量，再通过拟合Murnaghan方程：E(V)=E_0+\frac{B_0V_0}{B_0^\prime}\left[\left(\frac{V_0}{V}\right)^{B_0^\prime}-1\right]+B_0V\left(1-\frac{V_0}{V}\right)，其中E(V)是体系能量，E_0为基态能量，V_0是基态平衡体积，B_0是体模量，B_0^\prime是体模量对压强的导数，V是原胞体积。通过拟合该方程，可得到晶体的晶格常数以及单胞的能量，确定体系的基态能量和平衡体积，从而得到模型三个基矢的大小，确定晶体的平衡晶格常数。在求得晶格常数的基础上，根据研究的物理问题确定体系中包含原子数目的多少，建立第一性原理计算模型。模型通常选取一个或几个单胞（超单胞），且选取的模型需具有三个基矢方向，保证沿着三个基矢方向平移能构成无限大的晶体。输入的原子坐标有笛卡尔坐标（Cartesiancoordinates）和分数坐标（fractionalcoordinates）两种形式。笛卡尔坐标直接表示原子在三维空间中的位置，而分数坐标则是相对于晶胞基矢的比例坐标。对于具有面心立方结构的晶体，若取一个单胞作为研究模型，三个基矢分别为\vec{a}=(a,0,0)，\vec{b}=(0,a,0)，\vec{c}=(0,0,a)（a为晶格常数），单胞中原子的坐标既可用笛卡尔坐标表示，也可用分数坐标表示，二者可通过公式\vec{r}=x\vec{a}+y\vec{b}+z\vec{c}（其中\vec{r}为笛卡尔坐标，(x,y,z)为分数坐标）相互转换。完成模型构建后，设置计算参数，包括平面波截断能量、K点网格密度等。平面波截断能量决定了平面波基组的大小，截断能量越高，计算精度越高，但计算量也越大。K点网格密度则影响对布里渊区的采样精度，密度越高，对电子态的描述越准确。在计算ZnO基稀磁半导体时，合理设置平面波截断能量为500eV，K点网格采用6\times6\times6的密度，能在保证计算精度的同时，控制计算成本。随后进行自洽场迭代计算，从初始猜测的电子密度出发，求解Kohn-Sham方程得到新的电子密度，再将新的电子密度代入方程，反复迭代直至电子密度和能量收敛。当相邻两次迭代的能量差小于设定的收敛阈值（如10^{-6}eV/atom）时，认为计算收敛，得到体系的基态电子结构和能量。2.1.3应用软件在第一性原理计算领域，VASP（ViennaAb-initioSimulationPackage）是一款常用且功能强大的计算程序包。它基于密度泛函理论，采用平面波赝势方法，能够高效地计算材料的电子结构、几何结构和各种物理性质。VASP的使用流程相对规范。首先，准备输入文件，包括晶体结构信息（如原子坐标、晶格常数等）、赝势文件以及计算参数设置文件。可通过文本编辑器创建这些文件，将晶体结构信息按照特定格式写入坐标文件。接着，将准备好的输入文件提交到计算集群或服务器上运行。在运行过程中，VASP会读取输入文件，根据设置的参数进行计算，并将计算过程中的信息和结果输出到日志文件和其他结果文件中。计算完成后，通过分析输出文件，获取体系的能量、电子结构、电荷密度等信息。可利用专门的可视化软件（如VESTA）读取输出文件，直观地查看晶体结构和电荷密度分布等。VASP具有诸多优势。在计算效率方面，它采用了高效的算法和并行计算技术，能够充分利用计算资源，大大缩短计算时间。在处理大规模体系时，通过并行计算，可将计算任务分配到多个处理器核心上同时进行，显著提高计算速度。在精度方面，VASP提供了多种交换关联泛函和赝势选择，用户可根据研究体系的特点选择合适的参数，以获得高精度的计算结果。对于过渡金属掺杂的体系，选择合适的泛函和赝势，能准确描述体系的磁性和电子结构。VASP还具备广泛的适用性，可用于计算各种晶体材料、表面、界面以及分子体系的性质，在材料科学、化学、物理等多个领域都有广泛的应用。2.2蒙特卡洛模拟方法2.2.1理论基础蒙特卡洛模拟方法源于20世纪40年代美国在研制原子弹的“曼哈顿计划”，数学家冯・诺伊曼和乌拉姆等人为解决裂变物质的中子随机扩散问题而提出。其基本思想是基于概率统计理论，通过大量的随机抽样和统计计算来模拟复杂系统的行为或求解数学物理问题。该方法以摩纳哥的蒙特卡洛赌场命名，寓意其依赖于随机性和概率，如同赌博中的随机过程。在研究材料磁性时，蒙特卡洛模拟主要基于伊辛模型（Isingmodel）。伊辛模型将材料中的原子磁矩简化为只取向上（+1）或向下（-1）两种状态的自旋变量。以铁磁材料为例，原子磁矩间存在相互作用，当相邻原子磁矩同向时，体系能量降低，表现出铁磁性。体系的能量可表示为：E=-J\sum_{i,j}S_iS_j-\muH\sum_{i}S_i，其中J为交换相互作用常数，S_i和S_j分别表示第i和第j个原子的自旋，\mu为磁矩，H为外磁场。通过蒙特卡洛模拟，从大量随机的自旋状态中统计体系的能量、磁化强度等物理量，进而研究材料的磁性。如在模拟铁磁材料的居里温度时，随着温度升高，原子热运动加剧，自旋的无序度增加，当温度达到居里温度时，材料的铁磁性消失，通过蒙特卡洛模拟不同温度下体系的磁化强度变化，可确定居里温度。2.2.2实现过程蒙特卡洛模拟的实现步骤较为清晰。首先是初始条件设定，需确定模拟体系的大小，比如在模拟二维材料的磁性时，确定晶格的行数和列数。同时，随机初始化每个原子的自旋状态，使自旋向上和向下的原子随机分布。设置模拟的温度、外磁场强度等参数，温度决定了原子热运动的剧烈程度，外磁场则会影响原子磁矩的取向。随机抽样环节，采用梅森旋转算法等伪随机数生成器产生随机数。该算法生成的伪随机数具有良好的统计特性，接近理想的均匀分布。利用这些随机数来决定自旋的翻转。在给定的温度下，根据玻尔兹曼分布计算自旋翻转的概率P=\exp(-\DeltaE/kT)，其中\DeltaE为自旋翻转前后体系能量的变化，k为玻尔兹曼常数，T为温度。若生成的随机数小于P，则自旋发生翻转。状态更新是基于随机抽样的结果，对体系中原子的自旋状态进行更新。每次更新后，重新计算体系的能量和磁化强度等物理量。通过多次重复随机抽样和状态更新过程，使体系达到稳定状态。当体系的能量和磁化强度等物理量在多次更新后变化很小时，认为体系达到稳定。此时，对体系的物理量进行统计平均，得到模拟结果。如在模拟铁磁材料的磁滞回线时，逐渐改变外磁场强度，在每个外磁场下进行蒙特卡洛模拟，统计不同外磁场下体系的磁化强度，从而得到磁滞回线。2.2.3在本研究中的应用在研究ZnO基稀磁半导体时，蒙特卡洛模拟从构建合适的模型入手。考虑到ZnO的晶体结构为纤锌矿结构，构建包含磁性离子（如Mn、Co等）掺杂的ZnO超晶胞模型。确定超晶胞的大小和磁性离子的掺杂位置与浓度。设置模拟参数，如温度范围从低温到高于居里温度，步长为一定值（如5K），以研究材料磁性随温度的变化。外磁场强度也设置多个值，从0逐渐增加到饱和磁场，再逐渐减小，以获取磁滞回线。模拟的重点在于准确模拟磁性离子与ZnO基质中原子之间的交换相互作用，以及磁性离子之间的相互作用。通过调整交换相互作用常数J的值，使其与第一性原理计算结果或实验数据相匹配。在模拟过程中，重点关注体系的磁化强度、磁矩分布等物理量的变化。模拟的目标是深入了解ZnO基稀磁半导体的磁性起源和磁性随温度、外磁场的变化规律。通过模拟不同掺杂浓度下材料的磁性，分析掺杂浓度对居里温度、磁化强度等磁性参数的影响。还可以通过模拟探究材料中缺陷对磁性的影响，为实验研究提供理论指导和预测。2.3相关物理参量2.3.1交换相互作用交换相互作用是一种量子力学效应，在ZnO基稀磁半导体的磁性形成中扮演着关键角色。从本质上讲，它源于电子的全同性原理和泡利不相容原理。在ZnO晶体中，当磁性离子（如Mn、Co等）替代部分Zn离子时，磁性离子的d电子与ZnO基质中O的2p电子之间会产生强烈的交换相互作用。这种相互作用使得电子的自旋状态发生关联，进而影响材料的磁性。以Mn掺杂的ZnO基稀磁半导体为例，Mn离子的3d电子与O的2p电子通过p-d交换相互作用，形成了特定的自旋排列。在这种交换相互作用下，电子的自旋会倾向于保持特定的相对取向，从而产生净磁矩。这种交换相互作用的强度和性质决定了材料是呈现铁磁性还是反铁磁性。当交换相互作用使得相邻磁性离子的自旋平行排列时，材料表现出铁磁性；反之，当自旋反平行排列时，材料呈现反铁磁性。交换相互作用还与材料的晶体结构密切相关。ZnO的纤锌矿结构中，原子的排列方式影响了电子云的重叠程度，进而影响交换相互作用的强度。在不同的晶体方向上，交换相互作用的大小可能存在差异，这也导致材料的磁性在不同方向上表现出各向异性。2.3.2铁磁性和反铁磁性铁磁性是指材料在无外磁场时能够自发磁化，且磁化方向保持一致的特性。在ZnO基稀磁半导体中，若材料表现出铁磁性，其内部的磁性离子自旋会在一定温度范围内自发地平行排列，形成磁畴。当施加外磁场时，磁畴会逐渐转向外磁场方向，使得材料的磁化强度迅速增加。在一定的外磁场强度下，材料会达到饱和磁化状态，此时磁化强度不再随外磁场的增加而显著变化。反铁磁性则是指材料中相邻磁性离子的自旋呈反平行排列，整体磁矩为零的特性。在ZnO基稀磁半导体中，若存在反铁磁性相互作用，磁性离子的自旋会交替排列，形成反铁磁结构。这种结构在宏观上不表现出磁性，但在特定的实验条件下，如通过中子散射等技术，可以探测到其内部的反铁磁有序。判断ZnO基稀磁半导体是表现出铁磁性还是反铁磁性，主要通过磁滞回线测量、磁性温度依赖性测量等实验手段。磁滞回线中，若存在明显的剩磁和矫顽力，则表明材料具有铁磁性；而反铁磁性材料的磁滞回线通常较为平坦，无明显的剩磁和矫顽力。在磁性温度依赖性测量中，铁磁性材料的磁化强度会随温度升高而逐渐降低，当温度达到居里温度时，铁磁性消失；反铁磁性材料则在奈尔温度以下表现出反铁磁有序，超过奈尔温度后，反铁磁性消失，材料转变为顺磁性。三、非磁性金属Cu掺杂ZnO的研究3.1电子结构分析利用基于密度泛函理论的第一性原理计算方法，对非磁性金属Cu掺杂ZnO体系的电子结构展开深入研究。在计算过程中，构建了包含特定原子数目的ZnO超晶胞模型，并精确确定Cu原子在超晶胞中的替代位置。通过精心设置平面波截断能量、K点网格密度等关键计算参数，确保计算结果的准确性和可靠性。在平面波截断能量的设置上，经过多次测试和验证，选取了400eV这一数值，使得在保证计算精度的同时，有效控制了计算资源的消耗。对于K点网格密度，采用了较为精细的5\times5\times5的设置，以实现对布里渊区的充分采样，从而更准确地描述体系的电子态。在成功构建模型并设置好参数后，进行自洽场迭代计算。经过多轮迭代，最终获得了体系的电子结构信息。从计算得到的能带结构来看，本征ZnO呈现出典型的半导体能带特征，价带顶主要由O的2p电子构成，导带底则主要由Zn的4s电子组成，其禁带宽度为3.37eV。当Cu掺杂后，能带结构发生了显著变化。在费米能级附近出现了新的杂质能级，这一杂质能级主要是由Cu的3d电子与O的2p电子之间强烈的杂化作用所产生。通过对态密度的进一步分析，可以更清晰地了解这种杂化作用。从总态密度图中可以看出，在费米能级附近，出现了明显的态密度峰，这表明该区域存在电子态的富集。分波态密度图则显示，Cu的3d态密度与O的2p态密度在费米能级附近有显著的重叠，这直接证明了Cu-3d电子和O-2p电子之间存在着强烈的杂化作用。这种杂化作用使得电子云在Cu和O原子周围发生重新分布，从而对体系的磁性和电学性质产生了重要影响。从电子云分布的角度来看，在本征ZnO中，Zn原子和O原子之间通过离子键和共价键相互作用，形成了稳定的晶体结构。电子云在Zn和O原子周围呈现出较为规则的分布。当Cu掺杂后，由于Cu-3d电子与O-2p电子的杂化，电子云在Cu和O原子周围发生了明显的变形和重叠。在Cu原子周围，电子云密度增加，且与O原子的电子云相互交织，形成了一种更为复杂的电子云分布模式。这种电子云分布的变化不仅影响了原子间的化学键性质，还对体系的电子输运和磁性产生了深远影响。通过对电荷密度差的分析，可以进一步揭示电子的转移情况。电荷密度差图显示，在Cu掺杂后，部分电子从Zn原子转移到了Cu原子和O原子周围。这表明Cu的掺杂导致了体系内电荷的重新分布。具体而言，由于Cu的电负性与Zn存在差异，Cu原子吸引了周围的电子，使得电子云向Cu原子偏移。这种电子转移进一步增强了Cu-3d电子与O-2p电子之间的杂化作用，从而对体系的磁性产生了重要影响。对电子结构的分析表明，非磁性金属Cu掺杂ZnO体系中，Cu-3d电子和O-2p电子之间的强烈杂化作用是体系具有独特物理性质的关键因素。这种杂化作用导致了费米能级附近杂质能级的出现，电子云分布的改变以及电荷的重新分布，为进一步理解该体系的磁性和电学性质奠定了坚实的基础。3.2磁性来源探究3.2.1Cu-3d电子和O-2p电子的杂化作用在深入剖析非磁性金属Cu掺杂ZnO体系的磁性来源时，Cu-3d电子和O-2p电子之间强烈的杂化作用是关键因素。通过对电子结构的精细分析，从多个角度证实了这种杂化作用对体系磁性的重要贡献。从轨道重叠的层面来看，在ZnO晶体中，O原子的2p轨道与Zn原子的4s、4p轨道存在一定程度的重叠，形成了稳定的化学键。当Cu原子替代Zn原子后，Cu的3d轨道与O的2p轨道发生了更为显著的重叠。利用电子轨道分析软件对体系的轨道重叠情况进行模拟，结果显示，在Cu掺杂的区域，Cu-3d轨道与O-2p轨道的重叠积分明显增大，表明两者之间的相互作用增强。这种轨道重叠导致了电子云在Cu和O原子之间的重新分布，使得电子的运动范围扩展到了Cu和O原子周围，从而形成了一种更为复杂的电子云分布模式。从电子态密度的角度进一步分析，在费米能级附近，Cu的3d态密度与O的2p态密度出现了明显的重叠峰。通过对不同掺杂浓度下体系的态密度进行对比研究发现，随着Cu掺杂浓度的增加，费米能级附近的态密度峰强度逐渐增强，这意味着Cu-3d电子和O-2p电子之间的杂化作用逐渐增强。这种杂化作用使得电子在不同原子轨道之间的跃迁变得更加容易，从而对体系的磁性产生了重要影响。为了更直观地展示这种杂化作用，通过可视化软件绘制了体系的电荷密度分布图。在图中可以清晰地看到，在Cu原子和O原子之间，电荷密度明显增加，形成了一个电荷富集区域。这表明Cu-3d电子和O-2p电子之间存在着强烈的相互作用，电子云在两者之间发生了明显的重叠。从化学键的角度分析，Cu-3d电子和O-2p电子的杂化作用改变了原子间的化学键性质。在本征ZnO中，Zn-O键主要表现为离子键和共价键的混合。当Cu掺杂后，由于Cu-3d电子与O-2p电子的杂化，Cu-O键的共价成分增加，离子成分相对减少。通过键长和键角的分析也证实了这一点，在Cu掺杂后，Cu-O键的键长和键角与本征ZnO中的Zn-O键相比发生了明显的变化，这进一步表明了化学键性质的改变。这种化学键性质的改变对体系的电子结构和磁性产生了深远影响，使得体系的磁性发生了显著变化。3.2.2局域磁矩与巡游电子在非磁性金属Cu掺杂ZnO体系中，局域磁矩和巡游电子的存在及其相互作用对磁性的产生起着至关重要的作用。局域磁矩主要来源于磁性离子的未成对电子。在Cu掺杂ZnO体系中，Cu原子的3d轨道上存在未成对电子，这些未成对电子形成了局域磁矩。通过对体系中Cu原子的电子结构分析，发现其3d轨道上的电子自旋向上和自旋向下的电子数存在差异，这种差异导致了净磁矩的产生，从而形成了局域磁矩。不同掺杂浓度下，局域磁矩的大小和方向会发生变化。随着Cu掺杂浓度的增加，局域磁矩的大小逐渐增大，这是因为更多的Cu原子引入了更多的未成对电子，从而增强了局域磁矩。局域磁矩之间的相互作用也会影响体系的磁性。当局域磁矩之间的相互作用为铁磁耦合时，它们的磁矩方向趋于一致，从而增强了体系的磁性；当相互作用为反铁磁耦合时，局域磁矩的方向相反，体系的磁性会受到削弱。巡游电子则是在整个晶体中自由移动的电子。在Cu掺杂ZnO体系中，由于Cu-3d电子和O-2p电子的杂化作用，使得电子的离域性增强，形成了巡游电子。通过电子态密度分析和电子云分布的研究，发现费米能级附近的电子具有较高的迁移率，能够在晶体中自由移动，这些电子即为巡游电子。巡游电子与局域磁矩之间存在着相互作用。巡游电子可以通过与局域磁矩的交换相互作用，影响局域磁矩的取向。当巡游电子与局域磁矩发生交换相互作用时，巡游电子的自旋会与局域磁矩的自旋发生耦合，从而改变局域磁矩的取向。这种相互作用在低温下更为显著，因为低温下巡游电子的能量较低，更容易与局域磁矩发生耦合。巡游电子还可以通过传导电流的方式，对体系的磁性产生影响。当有电流通过体系时，巡游电子的运动方向会受到磁场的影响，从而产生一个附加的磁场，这个附加磁场会与局域磁矩相互作用，进一步影响体系的磁性。在Cu掺杂ZnO体系中，局域磁矩和巡游电子之间的相互作用是一个动态的过程。在不同的温度和磁场条件下，它们之间的相互作用会发生变化，从而导致体系磁性的变化。在低温下，局域磁矩之间的铁磁耦合作用较强，巡游电子与局域磁矩的耦合也较为稳定，体系表现出较强的铁磁性。随着温度的升高，热运动加剧，局域磁矩的取向逐渐变得无序，巡游电子与局域磁矩的耦合也受到影响，体系的磁性逐渐减弱。当温度达到居里温度时，局域磁矩的取向完全无序，体系的铁磁性消失，转变为顺磁性。3.3居里温度的计算与分析3.3.1第一性原理结合蒙特卡洛计算为了深入探究非磁性金属Cu掺杂ZnO体系的居里温度，采用第一性原理结合蒙特卡洛模拟的方法进行计算。首先，通过第一性原理计算获取体系中磁性离子之间的磁性耦合强度。在计算过程中，利用基于密度泛函理论的VASP软件，精确计算不同原子间的交换相互作用能。对于Cu掺杂ZnO体系，重点计算Cu原子与相邻O原子以及Cu原子之间的交换相互作用能。通过对不同原子构型下体系能量的计算，得到磁性离子间的交换相互作用常数J。在获得磁性耦合强度后，将其作为重要参数输入到蒙特卡洛模拟中。基于伊辛模型，构建二维或三维的晶格模型来模拟ZnO基稀磁半导体体系。在模拟过程中，设定每个晶格位点上的自旋变量，其取值为+1或-1，分别代表自旋向上和自旋向下。利用蒙特卡洛方法中的Metropolis算法，随机选择晶格位点并尝试翻转其自旋，根据体系能量的变化和玻尔兹曼分布来决定自旋是否真正翻转。在每个温度下，进行大量的蒙特卡洛步模拟，使体系达到平衡状态。通过统计体系在不同温度下的磁化强度，得到磁化强度随温度的变化曲线。当磁化强度随温度变化曲线出现明显的下降转折点时，该转折点所对应的温度即为居里温度。对于三种不同掺杂浓度的Cu掺杂ZnO体系，分别进行上述计算过程。在掺杂浓度为2.5%的体系中，经过第一性原理计算得到磁性离子间的交换相互作用常数J_1。将J_1代入蒙特卡洛模拟中，经过多次模拟和统计分析，得到该体系的居里温度T_{C1}。同样地，对于掺杂浓度为5%和7.5%的体系，分别得到交换相互作用常数J_2、J_3以及对应的居里温度T_{C2}、T_{C3}。通过这种方法，能够较为准确地计算出不同掺杂浓度下Cu掺杂ZnO体系的居里温度，为进一步研究掺杂浓度对居里温度的影响提供了数据支持。3.3.2居里温度与掺杂浓度的关系通过对三种不同掺杂浓度下Cu掺杂ZnO体系居里温度的计算，深入分析居里温度与掺杂浓度之间的关系。研究发现，随着Cu掺杂浓度的增加，体系的居里温度呈现出先升高后降低的变化趋势。在掺杂浓度较低时，如2.5%的掺杂浓度下，体系的居里温度相对较低。这是因为此时体系中磁性离子的数量较少，磁性离子之间的相互作用较弱，热运动相对容易破坏自旋的有序排列，导致居里温度较低。当掺杂浓度增加到5%时，体系的居里温度显著升高。这是由于随着磁性离子数量的增加，磁性离子之间的交换相互作用增强，自旋之间更容易形成有序排列，从而提高了体系的铁磁稳定性，使得居里温度升高。在这个过程中，Cu-3d电子与O-2p电子之间的杂化作用也随着掺杂浓度的增加而增强，进一步促进了自旋的有序排列，对居里温度的提升起到了积极作用。当掺杂浓度继续增加到7.5%时，体系的居里温度反而下降。这可能是由于过高的掺杂浓度导致磁性离子之间的距离过近，产生了较强的反铁磁相互作用，破坏了自旋的有序排列，从而降低了体系的铁磁稳定性，使得居里温度降低。过高的掺杂浓度还可能引入更多的缺陷和杂质，这些缺陷和杂质会干扰电子的传输和自旋的相互作用，对体系的磁性产生负面影响，进而降低居里温度。居里温度与掺杂浓度之间的这种关系对材料的实际应用具有重要影响。在自旋电子器件的应用中，需要材料具有较高的居里温度，以保证器件在室温下能够稳定工作。因此，在制备Cu掺杂ZnO基稀磁半导体材料时，需要精确控制掺杂浓度，选择合适的掺杂浓度范围，以获得较高的居里温度，满足实际应用的需求。对于需要在高温环境下工作的自旋电子器件，应选择居里温度较高的掺杂浓度；而对于对磁性均匀性要求较高的应用场景，则需要避免过高的掺杂浓度，以防止反铁磁相互作用的产生。3.4磁性耦合机制为了深入探究非磁性金属Cu掺杂ZnO体系的磁性耦合机制，通过计算电子转移情况来分析电子之间的相互作用特性。利用第一性原理计算方法，对体系中电子在不同原子间的转移进行了精确计算。计算结果表明，在Cu掺杂ZnO体系中，电子存在远程转移现象。具体而言，电子能够在Cu原子与相邻的O原子之间，以及不同的Cu原子之间进行远程转移。这种远程转移特性使得电子的运动范围不再局限于单个原子周围，而是扩展到了更大的空间尺度。从电子云分布的角度来看，电子的远程转移导致电子云在Cu和O原子周围发生了明显的重叠和扩展。在Cu原子周围，电子云密度较高，且与相邻O原子的电子云相互交织。随着距离的增加，电子云密度逐渐降低，但仍然存在一定的电子云分布，这表明电子能够在较大范围内进行转移。通过对不同原子间电子转移路径的分析，发现电子可以通过Cu-O键以及Cu-Cu键进行远程转移。在Cu-O键中，由于Cu-3d电子和O-2p电子的杂化作用，电子云在Cu和O原子之间形成了一个连续的分布区域，使得电子能够顺利地在两者之间转移。在Cu-Cu键中，虽然距离相对较远，但由于电子的离域性，仍然存在一定的电子转移概率。结合局域磁矩和巡游电子的作用进一步分析，发现该体系的磁性耦合是一种RKKY（Ruderman-Kittel-Kasuya-Yosida）耦合。如前文所述，局域磁矩主要来源于Cu原子的3d轨道上的未成对电子，这些未成对电子形成了具有一定方向和大小的局域磁矩。巡游电子则是在整个晶体中自由移动的电子，它们可以通过与局域磁矩的交换相互作用，影响局域磁矩的取向。在RKKY耦合机制中，巡游电子的自旋极化会在空间中产生一个振荡的自旋密度分布，这个振荡的自旋密度分布会与局域磁矩相互作用，从而导致局域磁矩之间产生间接的耦合作用。在Cu掺杂ZnO体系中，由于电子的远程转移，巡游电子能够在不同的局域磁矩之间传递信息，使得局域磁矩之间的耦合作用得以实现。当巡游电子经过一个局域磁矩时，其自旋会受到局域磁矩的影响而发生极化。这种极化的巡游电子再经过另一个局域磁矩时，会将第一个局域磁矩的信息传递给第二个局域磁矩，从而使两个局域磁矩之间产生耦合作用。这种耦合作用的大小和方向与巡游电子的自旋极化状态以及局域磁矩的相对位置有关。当巡游电子的自旋极化与局域磁矩的方向一致时，会增强局域磁矩之间的耦合作用，使它们的磁矩更倾向于平行排列，从而表现出铁磁性；当巡游电子的自旋极化与局域磁矩的方向相反时，会削弱局域磁矩之间的耦合作用，甚至导致它们的磁矩反平行排列，表现出反铁磁性。通过对电子转移的计算和分析，明确了非磁性金属Cu掺杂ZnO体系中电子的远程转移特性。结合局域磁矩和巡游电子的作用，证实了该体系的磁性耦合是一种RKKY耦合。这种耦合机制使得局域磁矩之间通过巡游电子的介导产生相互作用，从而决定了体系的磁性性质。对磁性耦合机制的深入理解，为进一步研究Cu掺杂ZnO体系的磁性行为以及开发基于该体系的自旋电子器件提供了重要的理论基础。四、非金属元素C掺杂ZnO的研究4.1电子结构与磁性基态运用基于密度泛函理论的第一性原理计算方法，对非金属元素C掺杂ZnO体系展开深入研究。构建包含特定原子数目的ZnO超晶胞模型，确保模型能够准确反映体系的真实结构。在模型中，精心设定C原子替代ZnO中O原子的位置，以模拟实际的掺杂情况。对平面波截断能量和K点网格密度等关键计算参数进行细致优化，通过多次测试和对比，最终确定平面波截断能量为450eV，K点网格采用6\times6\times6的密度，以保证计算结果的准确性和可靠性。通过自洽场迭代计算，成功获得体系的电子结构。从能带结构来看，本征ZnO呈现典型的半导体能带特征，其价带顶主要由O的2p电子构成，导带底主要由Zn的4s电子组成，禁带宽度约为3.37eV。当C掺杂后，能带结构发生显著变化。在费米能级附近出现了新的杂质能级，这一杂质能级主要源于C的2p电子与ZnO中Zn的4s、4p电子以及O的2p电子之间的杂化作用。通过对态密度的深入分析，进一步揭示了这种杂化作用的本质。在总态密度图中，费米能级附近出现明显的态密度峰，表明该区域存在电子态的富集。分波态密度图则清晰显示，C的2p态密度与Zn的4s、4p态密度以及O的2p态密度在费米能级附近存在显著重叠，充分证实了C与Zn、O原子之间存在强烈的杂化作用。这种杂化作用使得电子云在C、Zn和O原子周围发生重新分布，对体系的磁性和电学性质产生重要影响。为了进一步探究体系的磁性基态，对体系的总磁矩和原子磁矩进行精确计算。计算结果表明，C掺杂ZnO体系具有稳定的铁磁基态，总磁矩不为零。通过对各原子磁矩的分析发现，C原子周围的磁矩对体系总磁矩的贡献较大。这是因为C的2p电子与周围原子的电子发生杂化后，电子的自旋状态发生改变，形成了具有一定方向和大小的磁矩。在C掺杂的ZnO体系中，C原子的2p电子与相邻O原子的2p电子以及Zn原子的4s、4p电子相互作用，使得C原子周围的电子云分布发生变化，从而产生了明显的磁矩。这种磁矩的产生与C原子在ZnO晶格中的位置以及与周围原子的相互作用密切相关。不同位置的C原子，其与周围原子的距离和电子云重叠程度不同，导致产生的磁矩大小和方向也有所差异。通过对不同掺杂位置的C原子磁矩进行计算和分析，发现当C原子替代靠近Zn原子的O原子时，产生的磁矩较大，对体系总磁矩的贡献更为显著。这一结果表明，C原子在ZnO晶格中的位置对体系的磁性基态具有重要影响，通过合理控制C原子的掺杂位置，可以有效调控体系的磁性。4.2居里温度与平面磁矩分布4.2.1居里温度计算为了准确计算非金属元素C掺杂ZnO体系的居里温度，采用第一性原理与蒙特卡洛模拟相结合的方法。通过第一性原理计算，精确获取体系中原子间的交换耦合常数，以此作为蒙特卡洛模拟的关键输入参数。在蒙特卡洛模拟中，基于伊辛模型构建体系模型，将体系划分为多个晶格位点，每个位点代表一个原子，位点上的自旋变量取值为+1或-1，分别对应自旋向上和自旋向下两种状态。利用梅森旋转算法生成伪随机数，以此决定自旋的翻转。在给定温度下，依据玻尔兹曼分布计算自旋翻转的概率，若生成的随机数小于该概率，则自旋发生翻转。通过多次重复这一过程，使体系达到稳定状态。经过模拟计算，得到了不同掺杂浓度下C掺杂ZnO体系的居里温度。当C掺杂浓度为2.5%时，体系的居里温度约为350K；当掺杂浓度提高到5%时，居里温度升高至约400K；而当掺杂浓度达到7.5%时，居里温度进一步提升至约450K。与其他体系的居里温度进行对比分析，如Mn掺杂ZnO体系，其居里温度通常在200-300K之间，C掺杂ZnO体系展现出了相对较高的居里温度。这种差异源于不同掺杂原子与ZnO基质之间的相互作用不同。C原子的2p电子与ZnO中Zn和O原子的电子形成了独特的杂化作用，增强了体系中磁性离子间的交换相互作用，从而提高了居里温度。在一些传统的稀磁半导体体系中，由于磁性离子间的相互作用较弱，热运动容易破坏自旋的有序排列，导致居里温度较低。而C掺杂ZnO体系中，通过C原子的掺杂，有效增强了磁性离子间的相互作用，使得体系在较高温度下仍能保持自旋的有序排列，进而提高了居里温度。4.2.2平面磁矩分布研究通过深入分析C掺杂ZnO体系的平面磁矩分布，能够进一步揭示局域磁矩的存在和分布特点。利用第一性原理计算，获取体系中每个原子的磁矩信息。将体系划分为多个平面，统计每个平面内原子磁矩的矢量和，得到平面磁矩分布。在C掺杂ZnO体系中，发现存在明显的局域磁矩。在C原子周围，磁矩呈现出较强的局域化分布。这是由于C原子的2p电子与周围原子的电子发生杂化，导致电子云分布发生变化，从而形成了具有一定方向和大小的局域磁矩。通过对不同平面的磁矩分布进行分析，发现局域磁矩的方向并非完全一致，而是存在一定的角度分布。在某些平面内，局域磁矩呈现出有序排列的趋势，使得该平面内的净磁矩较大；而在其他平面内，局域磁矩的方向较为杂乱，净磁矩相对较小。不同掺杂浓度下，平面磁矩分布也存在差异。随着C掺杂浓度的增加，局域磁矩的数量增多，且分布范围扩大。在低掺杂浓度下，局域磁矩主要集中在C原子周围，相互之间的作用较弱；而在高掺杂浓度下，局域磁矩之间的相互作用增强，形成了更大范围的磁有序区域。这种平面磁矩分布的变化与体系的磁性密切相关。局域磁矩的有序排列是体系具有铁磁性的重要基础，随着掺杂浓度的增加，局域磁矩之间的相互作用增强，促进了磁有序区域的形成和扩大，从而增强了体系的铁磁性。4.3电子转移与磁性耦合倾向为了深入探究非金属元素C掺杂ZnO体系的磁性耦合机制，对体系中电子的转移情况进行了细致计算。利用第一性原理计算方法，精确分析了C原子与相邻Zn、O原子之间的电子转移。计算结果表明，C原子与周围原子之间存在明显的电子转移现象。具体而言，C原子的2p电子与相邻O原子的2p电子以及Zn原子的4s、4p电子之间发生了电子转移。这种电子转移导致电子云在C、Zn和O原子周围发生了重新分布。通过对电荷密度差的分析，可以清晰地看到在C原子与相邻原子之间，电荷密度发生了显著变化。在C-O键之间，电荷密度增加，表明电子云在这一区域出现了聚集；而在C-Zn键之间，电荷密度也有一定程度的改变，显示出电子的转移对化学键性质产生了影响。为了进一步分析磁性耦合倾向，计算了S、P、d三个轨道上自旋向上和自旋向下电子数的差值。在C掺杂ZnO体系中，C原子的2p轨道以及与C原子相邻的O原子和Zn原子的相关轨道上，自旋向上和自旋向下电子数存在明显差异。这种差异导致了净磁矩的产生，进而影响了磁性耦合。从计算结果来看，自旋向上和自旋向下电子数的差值在空间上呈现出一定的分布规律。在C原子周围，这种差值较为显著，表明局域磁矩较大；随着距离C原子的增加，差值逐渐减小，局域磁矩也相应减弱。这种变化趋势与RKKY耦合机制中局域磁矩通过巡游电子相互作用的特点相符合。在RKKY耦合中，巡游电子的自旋极化会在空间中产生一个振荡的自旋密度分布，这个分布与局域磁矩相互作用，导致局域磁矩之间产生间接的耦合作用。在C掺杂ZnO体系中，由于电子的转移，形成了一定数量的巡游电子。这些巡游电子的自旋极化与局域磁矩相互作用，使得局域磁矩之间的耦合作用呈现出远程的特性。当巡游电子经过一个局域磁矩时，其自旋会受到局域磁矩的影响而发生极化。这种极化的巡游电子再经过另一个局域磁矩时，会将第一个局域磁矩的信息传递给第二个局域磁矩，从而使两个局域磁矩之间产生耦合作用。这种耦合作用的大小和方向与巡游电子的自旋极化状态以及局域磁矩的相对位置有关。当巡游电子的自旋极化与局域磁矩的方向一致时，会增强局域磁矩之间的耦合作用，使它们的磁矩更倾向于平行排列，从而表现出铁磁性；当巡游电子的自旋极化与局域磁矩的方向相反时，会削弱局域磁矩之间的耦合作用，甚至导致它们的磁矩反平行排列，表现出反铁磁性。综合电子转移情况和自旋轨道电子数差值的分析结果，可以判断该稀磁半导体的磁性耦合倾向于RKKY耦合。五、Al和S掺杂ZnO的研究5.1非磁性结果分析通过第一性原理计算，对Al和S掺杂ZnO体系的电子结构进行深入剖析，发现其呈现非磁性的特性。从电子结构角度来看，在本征ZnO中，Zn原子的电子构型为[Ar]3d^{10}4s^{2}，O原子的电子构型为[He]2s^{2}2p^{4}。在形成ZnO晶体时，Zn的4s电子与O的2p电子相互作用，形成稳定的化学键，使体系的电子结构处于相对稳定的状态。当Al掺杂ZnO时，Al原子替代部分Zn原子。Al的电子构型为[Ne]3s^{2}3p^{1}，其外层电子数比Zn少。掺杂后，Al的3s和3p电子与周围O原子的2p电子相互作用。通过对态密度的分析，发现Al掺杂后，在费米能级附近并没有出现明显的自旋极化现象，即自旋向上和自旋向下的电子态密度几乎相同。这表明Al的掺杂没有引入局域磁矩，体系整体呈现非磁性。从电子云分布的角度来看，Al-O键的形成并没有导致电子云的自旋极化分布，电子云在空间中的分布相对均匀，没有形成具有磁矩的局域区域。当S掺杂ZnO时，S原子替代部分O原子。S的电子构型为[Ne]3s^{2}3p^{4}，与O的电子构型有一定相似性。S掺杂后，S的3p电子与周围Zn原子的4s和4p电子相互作用。同样，通过态密度分析，在费米能级附近也未出现明显的自旋极化，自旋向上和自旋向下的电子态密度基本一致，说明S的掺杂也没有产生局域磁矩，体系保持非磁性。从电荷密度差的角度分析，S-Zn键的形成过程中，电荷的转移并没有导致自旋极化的产生，电子的分布没有呈现出与磁性相关的特征。对比本征ZnO以及Cu、C掺杂体系，在本征ZnO中，由于原子间的相互作用，电子结构处于稳定的非磁性状态。而Cu掺杂体系中，Cu-3d电子和O-2p电子之间强烈的杂化作用导致了自旋极化，产生了局域磁矩，使体系具有磁性。C掺杂体系中，C的2p电子与Zn、O原子的电子杂化，形成了具有磁矩的局域区域，体系呈现铁磁性。Al和S掺杂ZnO体系中，无论是Al与O的相互作用，还是S与Zn的相互作用，都没有引发电子结构的自旋极化，无法形成有效的局域磁矩，这是导致体系呈现非磁性的根本原因。5.2电子转移特性对Al和S掺杂ZnO体系中电子转移情况的研究，有助于深入理解其物理特性。利用第一性原理计算，精确分析了电子在不同原子间的转移路径和数量。计算结果显示，在该体系中电子存在远程转移现象。以Al掺杂为例，电子能够在Al原子与相邻的O原子之间进行远程转移。这是由于Al原子的外层电子与O原子的电子之间存在相互作用，使得电子云在两者之间发生了重叠和扩展。在S掺杂的情况下，电子同样可以在S原子与周围的Zn原子之间实现远程转移。通过对电荷密度差的分析，可以清晰地观察到电子云在原子间的重新分布。在Al-O键和S-Zn键周围，电荷密度发生了明显的变化，表明电子在这些区域存在转移现象。虽然存在电子的远程转移，但体系中不存在局域磁矩。在Cu掺杂ZnO体系中，由于Cu-3d电子和O-2p电子之间的杂化作用，产生了局域磁矩，使得体系具有磁性。而在C掺杂ZnO体系中，C原子的2p电子与Zn、O原子的电子杂化，也形成了具有磁矩的局域区域，导致体系呈现铁磁性。在Al和S掺杂ZnO体系中，尽管电子存在远程转移，但Al与O、S与Zn之间的相互作用并没有引发电子结构的自旋极化，无法形成有效的局域磁矩。从电子轨道的角度分析，Al和S原子的电子轨道与周围原子的电子轨道相互作用时，没有产生明显的自旋极化现象，电子的自旋向上和自旋向下的分布较为均匀，这使得体系整体不具备磁性。这种电子转移特性和局域磁矩的差异，导致Al和S掺杂ZnO体系与其他磁性掺杂体系在磁性表现上存在显著不同。5.3对稀磁半导体磁性来源的启示通过对Al和S掺杂ZnO体系的研究，结合前文对Cu、C掺杂体系的分析，为理解稀磁半导体磁性来源和磁耦合机制提供了多方面的启示。在磁性来源方面，金属原子与非金属原子之间的杂化作用至关重要。在Cu掺杂ZnO体系中，Cu-3d电子和O-2p电子之间强烈的杂化作用，使得电子云在Cu和O原子周围发生重新分布，产生了局域磁矩，从而使体系具有磁性。C掺杂ZnO体系中，C的2p电子与Zn、O原子的电子杂化，形成了具有磁矩的局域区域，导致体系呈现铁磁性。而Al和S掺杂ZnO体系中，由于Al与O、S与Zn之间的相互作用未引发电子结构的自旋极化，无法形成有效的局域磁矩，体系呈现非磁性，这从反面证明了杂化作用对磁性产生的关键影响。从磁耦合机制角度来看，RKKY耦合在稀磁半导体中具有重要作用。在Cu掺杂ZnO体系中，通过计算电子转移发现存在巡游电子，巡游电子与局域磁矩磁性耦合导致了铁磁性，证实了磁性耦合是一种RKKY耦合。C掺杂ZnO体系中，通过计算S、P、d三个轨道上自旋向上和自旋向下电子数的差值，判断其磁性耦合倾向于RKKY耦合。Al和S掺杂ZnO体系虽无磁性，但电子存在远程转移现象，这与RKKY耦合中巡游电子的作用相关，进一步说明巡游电子在磁耦合机制中的重要性。这表明在稀磁半导体中，磁耦合机制与电子的转移特性以及局域磁矩和巡游电子的相互作用密切相关。通过对不同掺杂体系的研究，可以更深入地理解稀磁半导体的磁性来源和磁耦合机制，为进一步开发和优化稀磁半导体材料提供理论基础。六、ZnO基稀磁半导体磁性来源与磁耦合机制总结6.1磁性来源综合分析通过对非磁性金属Cu掺杂ZnO体系和非金属元素C掺杂ZnO体系的深入研究，发现金属原子和非金属原子之间的杂化作用在稀磁半导体磁性来源中起着关键作用。在Cu掺杂ZnO体系中，Cu-3d电子和O-2p电子之间的杂化作用表现得极为显著。从电子云分布角度来看，这种杂化使得电子云在Cu和O原子周围发生明显的变形和重叠，形成了更为复杂的电子云分布模式。在Cu原子周围，电子云密度增加，且与O原子的电子云相互交织。从态密度分析结果可知，在费米能级附近，Cu的3d态密度与O的2p态密度出现了明显的重叠峰，随着Cu掺杂浓度的增加，该重叠峰强度逐渐增强，进一步证明了杂化作用的增强。这种杂化作用导致了电子在不同原子轨道之间的跃迁变得更加容易，使得体系中产生了局域磁矩。这些局域磁矩的存在是体系具有磁性的重要基础。在C掺杂ZnO体系中，C的2p电子与Zn、O原子的电子之间同样存在强烈的杂化作用。通过对态密度的分析，清晰地看到C的2p态密度与Zn的4s、4p态密度以及O的2p态密度在费米能级附近存在显著重叠。从电子云分布情况来看，在C原子与相邻Zn、O原子之间，电子云发生了明显的重新分布。这种杂化作用使得C原子周围的电子云分布发生变化，从而产生了明显的磁矩。不同位置的C原子，由于其与周围原子的距离和电子云重叠程度不同，导致产生的磁矩大小和方向也有所差异。对比这两个体系，发现金属原子（如Cu）和非金属原子（如C）与ZnO基质中的原子形成杂化作用时，虽然具体的原子轨道参与杂化有所不同，但都能导致电子云分布的改变和局域磁矩的产生。这种局域磁矩的产生是稀磁半导体具有磁性的重要前提。而且，随着掺杂原子浓度的增加，杂化作用增强，局域磁矩也相应增强。在Cu掺杂体系中，随着Cu掺杂浓度从2.5%增加到7.5%，体系的磁性逐渐增强，居里温度也随之发生变化；在C掺杂体系中，随着C掺杂浓度的增加，局域磁矩数量增多，分布范围扩大，体系的铁磁性也得到增强。6.2磁耦合机制的统一认识通过对非磁性金属Cu掺杂ZnO体系和非金属元素C掺杂ZnO体系的深入研究，发现RKKY耦合在这两个体系中均发挥着重要作用。在Cu掺杂ZnO体系中，通过对电子转移的计算和分析，明确了存在巡游电子，且电子能够在Cu原子与相邻的O原子之间，以及不同的Cu原子之间进行远程转移。这种远程转移特性使得巡游电子能够在不同的局域磁矩之间传递信息，从而实现了局域磁矩之间的耦合作用。结合局域磁矩和巡游电子的作用，证实了该体系的磁性耦合是一种RKKY耦合。在这种耦合机制中，巡游电子的自旋极化会在空间中产生一个振荡的自旋密度分布，这个振荡的自旋密度分布会与局域磁矩相互作用，从而导致局域磁矩之间产生间接的耦合作用。当巡游电子经过一个局域磁矩时，其自旋会受到局域磁矩的影响而发生极化。这种极化的巡游电子再经过另一个局域磁矩时，会将第一个局域磁矩的信息传递给第二个局域磁矩，从而使两个局域磁矩之间产生耦合作用。这种耦合作用的大小和方向与巡游电子的自旋极化状态以及局域磁矩的相对位置有关。当巡游电子的自旋极化与局域磁矩的方向一致时，会增强局域磁矩之间的耦合作用，使它们的磁矩更倾向于平行排列，从而表现出铁磁性；当巡游电子的自旋极化与局域磁矩的方向相反时，会削弱局域磁矩之间的耦合作用，甚至导致它们的磁矩反平行排列，表现出反铁磁性。在C掺杂ZnO体系中，通过计算S、P、d三个轨道上自旋向上和自旋向下电子数的差值，判断其磁性耦合倾向于RKKY耦合。从电子转移情况来看，C原子与周围的Zn、O原子之间存在明显的电子转移现象，导致电子云在C、Zn和O原子周围发生了重新分布。这种电子转移形成了一定数量的巡游电子，巡游电子与局域磁矩之间的相互作用使得局域磁矩之间的耦合作用呈现出远程的特性。这与RKKY耦合中巡游电子介导局域磁矩之间相互作用的特点相符合。对比这两个体系，虽然掺杂原子不同，原子间的杂化作用细节有所差异，但在磁耦合机制方面，都体现出RKKY耦合的特征。这表明RKKY耦合在ZnO基稀磁半导体中具有一定的普遍性。无论是金属原子掺杂还是非金属原子掺杂，只要体系中存在巡游电子和局域磁矩，且巡游电子能够在不同的局域磁矩之间传递信息，就有可能通过RKKY耦合实现局域磁矩之间的相互作用，从而决定体系的磁性性质。这