基于第一性原理的稀磁半导体锐钛矿二氧化钛掺杂改性研究

基于第一性原理的稀磁半导体锐钛矿二氧化钛掺杂改性研究一、引言1.1研究背景与意义在当今社会，能源与环境问题已成为全球关注的焦点。随着工业化进程的加速，传统能源的过度消耗导致能源短缺问题日益严峻，同时环境污染问题也愈发严重，如大气污染、水污染和土壤污染等，给人类的生存和发展带来了巨大威胁。在众多解决能源与环境问题的技术中，光催化技术因其能够利用太阳能驱动化学反应，实现污染物的降解和能源的转化，被认为是一种极具潜力的绿色技术。二氧化钛（TiO₂）作为一种重要的半导体材料，在光催化领域展现出了独特的优势。TiO₂具有化学性质稳定、催化活性高、价格低廉、无毒无害等优点，被广泛应用于光催化降解有机污染物、光解水制氢、太阳能电池等领域。然而，TiO₂也存在一些局限性，其中最主要的问题是其禁带宽度较大（锐钛矿型TiO₂的禁带宽度为3.2eV），只能吸收波长小于387.5nm的紫外光，而紫外光在太阳光中所占的比例仅约为5%，这使得TiO₂对太阳能的利用率较低。此外，TiO₂光生载流子（电子-空穴对）的复合率较高，导致其光量子产率较低，进一步限制了其光催化性能的提升。为了克服TiO₂的这些缺点，提高其光催化性能，研究者们提出了多种改性方法，如离子掺杂、贵金属沉积、半导体复合、染料敏化等。其中，离子掺杂是一种简单有效的改性方法，通过向TiO₂晶格中引入杂质离子，可以改变其晶体结构、电子结构和光学性质，从而拓展其光响应范围，提高光生载流子的分离效率，进而提升其光催化活性。稀磁半导体（DilutedMagneticSemiconductors，DMS）是指在非磁性半导体中引入少量磁性离子而形成的一类新型半导体材料。DMS不仅具有半导体的电学和光学性质，还具有磁性，这种独特的性质使其在自旋电子学、磁传感器、高密度存储等领域展现出了广阔的应用前景。将磁性离子掺杂到TiO₂中，有望制备出具有光催化和磁性双重功能的稀磁半导体材料，这种材料不仅可以利用光催化性能降解污染物，还可以利用磁性实现材料的分离和回收，为解决光催化材料的分离难题提供了新的思路。目前，对于锐钛矿二氧化钛的掺杂改性研究主要集中在实验方面，通过各种实验手段制备掺杂TiO₂材料，并对其光催化性能进行测试和分析。然而，实验研究往往受到制备条件、测试方法等因素的影响，难以深入揭示掺杂对TiO₂结构和性能的影响机制。理论计算方法可以从原子和电子层面研究材料的结构和性质，为实验研究提供理论指导和预测。第一性原理计算基于量子力学原理，无需任何经验参数，能够准确地计算材料的晶体结构、电子结构、磁性和光学性质等，为研究锐钛矿二氧化钛的掺杂改性机制提供了有力的工具。综上所述，本研究采用第一性原理方法，深入研究不同离子掺杂对锐钛矿二氧化钛结构、电子性质、磁性及光学性质的影响，揭示掺杂改性的微观机制，为设计和制备高性能的光催化和稀磁半导体材料提供理论依据，具有重要的理论意义和实际应用价值。1.2国内外研究现状自1972年Fujishima和Honda首次发现TiO₂光解水现象以来，TiO₂作为一种重要的光催化材料受到了广泛的关注。国内外众多科研团队围绕TiO₂的掺杂改性展开了大量的研究工作，取得了一系列的研究成果。在实验研究方面，早期的工作主要集中在过渡金属离子掺杂TiO₂。研究发现，过渡金属离子（如Fe、Co、Ni、Cr等）的掺杂可以在TiO₂的禁带中引入杂质能级，使TiO₂的吸收光谱发生红移，从而拓展其光响应范围至可见光区域。例如，Wang等通过溶胶-凝胶法制备了Fe掺杂的TiO₂纳米颗粒，实验结果表明，Fe的掺杂显著提高了TiO₂在可见光下对罗丹明B的光催化降解活性。然而，过渡金属离子掺杂也可能引入一些负面影响，如成为光生载流子的复合中心，导致光量子产率降低。随着研究的深入，非金属离子掺杂TiO₂逐渐成为研究热点。常见的非金属掺杂元素包括N、C、S、F等。其中，N掺杂TiO₂的研究最为广泛。Asahi等首次报道了N掺杂TiO₂可以实现可见光响应，他们认为N原子的2p轨道与O原子的2p轨道杂化，在TiO₂的价带上方形成了新的杂质能级，从而使TiO₂能够吸收可见光。此后，许多研究小组通过不同的制备方法对N掺杂TiO₂进行了深入研究，并对其光催化性能的提升机制进行了探讨。除了单一元素掺杂，共掺杂也是提高TiO₂光催化性能的有效方法。研究表明，通过合理选择掺杂元素和掺杂浓度，共掺杂可以产生协同效应，进一步提高TiO₂的光催化活性。例如，Li等制备了N和Fe共掺杂的TiO₂，发现共掺杂样品的光催化活性明显高于单掺杂样品。在稀磁半导体方面，虽然TiO₂基稀磁半导体的研究起步较晚，但近年来取得了显著的进展。通过将磁性离子（如Mn、Co、Ni等）掺杂到TiO₂中，在一定条件下可以使TiO₂表现出室温铁磁性。然而，对于TiO₂基稀磁半导体中铁磁性的起源仍存在争议，主要观点包括载流子介导的铁磁相互作用、双交换机制、束缚磁极化子模型等。在理论研究方面，第一性原理计算已成为研究TiO₂掺杂改性的重要手段。通过第一性原理计算，可以精确地计算掺杂体系的晶体结构、电子结构、态密度、磁矩等性质，从而深入理解掺杂对TiO₂性能的影响机制。许多研究利用第一性原理计算研究了过渡金属离子掺杂TiO₂的电子结构和光学性质。如Wei等采用基于密度泛函理论的平面波赝势方法，研究了Co（Ⅱ）掺杂锐钛矿TiO₂的晶体结构、能隙和态密度，结果表明Co（Ⅱ）的3d轨道进入禁带并与氧的2p轨道作用，使禁带宽度变小，TiO₂的吸收带红移，从而有助于提高TiO₂的光催化活性。对于非金属离子掺杂TiO₂，理论计算也揭示了其电子结构的变化和光催化活性提升的内在机制。例如，Zhang等通过第一性原理计算研究了N掺杂锐钛矿TiO₂的电子结构，发现N掺杂在TiO₂的价带顶引入了杂质能级，减小了禁带宽度，增强了可见光吸收能力。尽管国内外在锐钛矿二氧化钛掺杂改性方面取得了丰硕的研究成果，但仍存在一些不足之处。一方面，实验研究中不同的制备方法和条件可能导致结果的差异，难以建立统一的掺杂改性规律；另一方面，理论计算中对于一些复杂的掺杂体系和磁性起源的解释还存在争议，需要进一步深入研究。此外，目前对于光催化性能和磁性的协同调控研究还相对较少，如何通过掺杂实现TiO₂材料在光催化和磁性方面的最优性能组合，仍是一个亟待解决的问题。1.3研究内容与方法1.3.1研究内容本研究主要采用第一性原理方法，深入研究不同离子掺杂对锐钛矿二氧化钛的结构、电子性质、磁性及光学性质的影响，具体研究内容如下：不同离子掺杂锐钛矿二氧化钛的晶体结构优化：构建不同离子（包括过渡金属离子、非金属离子和稀土金属离子）掺杂锐钛矿二氧化钛的超晶胞模型，利用第一性原理计算方法对其进行结构优化，得到稳定的晶体结构。分析掺杂离子对TiO₂晶格参数、键长、键角等结构参数的影响，探讨掺杂离子与TiO₂晶格之间的相互作用。掺杂体系的电子结构分析：计算掺杂前后锐钛矿二氧化钛的能带结构、态密度和电荷密度分布，研究掺杂离子对TiO₂电子结构的影响。分析杂质能级的形成位置和性质，以及掺杂离子的价电子与TiO₂中原子的电子之间的相互作用，揭示掺杂对TiO₂光催化活性和磁性产生影响的电子结构根源。磁性研究：对于掺杂磁性离子的锐钛矿二氧化钛体系，计算其磁矩和磁相互作用，研究掺杂离子浓度、掺杂位置等因素对体系磁性的影响。探讨TiO₂基稀磁半导体中铁磁性的起源和机制，为实现室温铁磁性提供理论指导。光学性质计算：通过第一性原理计算掺杂前后锐钛矿二氧化钛的光学性质，如吸收光谱、反射光谱和透射光谱等。分析掺杂对TiO₂光吸收边的影响，解释掺杂体系光响应范围拓展的原因，为提高TiO₂在可见光下的光催化性能提供理论依据。1.3.2研究方法本研究基于密度泛函理论（DensityFunctionalTheory，DFT），采用ViennaAbinitioSimulationPackage（VASP）软件进行第一性原理计算。在计算过程中，采用平面波赝势方法描述电子与离子实之间的相互作用，电子交换关联能采用广义梯度近似（GeneralizedGradientApproximation，GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函进行处理。具体计算参数设置如下：平面波截断能：经过测试，选取平面波截断能为500eV，以保证计算结果的准确性和收敛性。k点网格：对于体相TiO₂，采用Monkhorst-Pack方法生成k点网格，在倒易空间中选取合适的k点密度进行积分，以精确计算体系的电子结构和能量。对于超晶胞模型，根据体系的大小和对称性，合理调整k点网格，确保计算精度。结构优化：在结构优化过程中，采用共轭梯度法对原子坐标和晶格参数进行全优化，直到原子间的受力小于0.01eV/Å，体系的总能量收敛到10⁻⁵eV以内，得到稳定的晶体结构。自洽场计算：在进行电子结构和光学性质计算之前，先进行自洽场（Self-ConsistentField，SCF）迭代计算，使体系的电子密度和能量达到自洽收敛，以获得准确的计算结果。二、相关理论基础2.1第一性原理概述第一性原理，又称从头算，是基于量子力学基本原理的计算方法。其核心在于从最基本的物理定律出发，不依赖任何经验参数，直接求解体系的薛定谔方程，从而获得体系的电子结构和各种物理性质。在量子力学中，体系的状态由波函数\Psi描述，而薛定谔方程则是波函数随时间演化的基本方程，其定态形式为：\hat{H}\Psi=E\Psi其中，\hat{H}是哈密顿算符，代表体系的总能量，它包含了电子的动能、电子与原子核之间的相互作用势能以及电子之间的相互作用势能；E是体系的能量本征值；\Psi是能量本征态波函数。对于多电子体系，由于电子之间存在复杂的相互作用，直接求解薛定谔方程非常困难，甚至在实际计算中几乎无法实现。为了简化计算，人们发展了一系列的近似方法，如哈特里-福克（Hartree-Fock）方法和密度泛函理论（DensityFunctionalTheory，DFT）等。其中，密度泛函理论是目前第一性原理计算中应用最为广泛的方法。2.2密度泛函理论密度泛函理论的基本思想是将多电子体系的基态能量表示为电子密度的泛函。该理论的建立基于两个重要的定理：Hohenberg-Kohn定理和Kohn-Sham方程。Hohenberg-Kohn第一定理指出，对于一个处在外部势场V_{ext}(\mathbf{r})中的多电子体系，其基态电子密度\rho(\mathbf{r})唯一地确定了体系的所有性质，包括基态能量和波函数。这意味着可以将体系的能量表示为电子密度的泛函E[\rho]，而无需直接求解复杂的多电子波函数。Hohenberg-Kohn第二定理进一步表明，对于给定的外部势场，体系的基态能量是电子密度泛函E[\rho]的最小值，并且对应的电子密度就是体系的基态电子密度。基于这两个定理，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]\phi_i(\mathbf{r})=\epsilon_i\phi_i(\mathbf{r})其中，-\frac{\hbar^2}{2m}\nabla^2是电子的动能算符；V_{ext}(\mathbf{r})是外部势场，即原子核产生的静电势；V_{H}(\mathbf{r})是Hartree势，表示电子之间的库仑相互作用；V_{xc}(\mathbf{r})是交换关联势，用于描述电子之间的交换和关联效应，这是密度泛函理论中最难以精确描述的部分；\phi_i(\mathbf{r})是Kohn-Sham轨道，\epsilon_i是相应的本征能量。在实际计算中，需要对交换关联势V_{xc}(\mathbf{r})进行近似处理。常见的近似方法有局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）。LDA假设在空间某一点的交换关联能只与该点的电子密度有关，它采用均匀电子气模型来计算交换关联能，适用于电子密度变化较为缓慢的体系。其交换关联能泛函的一般形式为：E_{xc}^{LDA}[\rho]=\int\rho(\mathbf{r})\epsilon_{xc}(\rho(\mathbf{r}))d\mathbf{r}其中，\epsilon_{xc}(\rho(\mathbf{r}))是局域的交换关联能密度。GGA则在LDA的基础上，考虑了电子密度的空间梯度对交换关联能的影响，能更好地描述电子密度变化较快的体系，对于许多材料体系能够提供更精确的结果。GGA泛函的一般形式为：E_{xc}^{GGA}[\rho]=\int\rho(\mathbf{r})\epsilon_{xc}(\rho(\mathbf{r}),\nabla\rho(\mathbf{r}))d\mathbf{r}除了LDA和GGA，还有其他更高级的交换关联泛函，如杂化泛函、meta-GGA泛函等，它们在不同程度上提高了计算精度，但计算复杂度也相应增加。2.3赝势方法在第一性原理计算中，原子核与内层电子的相互作用非常复杂，直接处理会导致计算量急剧增加。赝势方法的提出就是为了简化这种计算，它通过构建一个赝势来替代原子核与内层电子对价电子的复杂相互作用。赝势的基本原理是：在离原子核较远的区域，价电子的波函数与真实波函数相同；而在离原子核较近的区域，通过调整赝势的形式，使赝波函数变得平滑，从而避免了对原子核附近复杂势场的精确描述，大大减少了计算量。同时，要求赝波函数和真实波函数具有相同的能量本征值和电荷密度，以保证计算结果的准确性。常用的赝势类型有模守恒赝势（Norm-ConservingPseudopotential）、超软赝势（UltraSoftPseudopotential，USP）和投影缀加波赝势（ProjectorAugmentedWavePseudopotential，PAW）。模守恒赝势要求赝波函数和真实波函数在原子核附近的正交性和归一性保持一致，对价电子的描述较为准确，但计算量相对较大。超软赝势通过引入投影算子和截断能，进一步简化了计算，降低了平面波基组的截断能要求，适用于计算较大的体系，但在某些情况下可能会牺牲一定的精度。PAW赝势则结合了平面波方法和缀加波方法的优点，在保持高精度的同时，提高了计算效率，被广泛应用于各种材料体系的计算。2.4晶体结构与模型构建锐钛矿二氧化钛（TiO₂）属于四方晶系，空间群为I4_1/amd（No.141）。其晶体结构可以看作是由TiO₆八面体通过共边连接而成，Ti原子位于八面体的中心，O原子位于八面体的顶点。在锐钛矿TiO₂的晶胞中，包含4个Ti原子和8个O原子，晶格参数a=b=0.3785nm，c=0.9514nm，α=β=γ=90°。每个Ti原子与周围的6个O原子形成八面体配位，而每个O原子则与3个Ti原子相连。为了研究离子掺杂对锐钛矿二氧化钛的影响，需要构建掺杂改性的超晶胞模型。通常采用的方法是在原胞的基础上，通过周期性重复构建一个包含多个原胞的超晶胞，然后将掺杂离子替换超晶胞中的部分Ti原子或O原子。例如，若要构建一个2\times2\times2的超晶胞，则该超晶胞中包含32个Ti原子和64个O原子，此时可以将其中的1个或多个Ti原子替换为掺杂离子，从而形成掺杂体系。在构建超晶胞模型后，需要对模型进行优化，以得到稳定的晶体结构。结构优化的目的是使体系的总能量最低，原子间的受力达到平衡。在第一性原理计算中，通常采用共轭梯度法、BFGS（Broyden-Fletcher-Goldfarb-Shanno）算法等优化算法对原子坐标和晶格参数进行全优化，直到原子间的受力小于某个阈值（如0.01eV/Å），体系的总能量收敛到一定精度（如10⁻⁵eV以内）。通过结构优化，可以得到掺杂体系的稳定结构，为后续的电子结构、磁性和光学性质计算提供基础。三、本征锐钛矿二氧化钛的性质计算与分析3.1计算参数设置在对本征锐钛矿二氧化钛进行第一性原理计算时，合理设置计算参数对于获得准确可靠的结果至关重要。平面波截断能决定了平面波基组对电子波函数的描述精度，其取值需足够大以保证计算的收敛性和准确性，但过大的截断能会显著增加计算量。通过对不同截断能（如400eV、450eV、500eV、550eV）下体系总能量的测试，发现当截断能为500eV时，体系总能量的变化已小于10⁻⁵eV，计算结果收敛良好，因此选取500eV作为平面波截断能。k点网格的选取则影响着对布里渊区积分的精度，进而影响电子结构计算的准确性。对于体相锐钛矿二氧化钛，采用Monkhorst-Pack方法生成k点网格。在测试过程中，分别尝试了不同的k点密度，如3×3×3、4×4×4、5×5×5、6×6×6。计算结果表明，当k点网格为5×5×5时，体系的总能量和电子结构已基本收敛，继续增加k点密度对结果影响较小，但计算时间会大幅增加。综合考虑计算精度和效率，最终确定k点网格为5×5×5。收敛精度是保证计算结果可靠性的关键因素之一。在结构优化过程中，采用共轭梯度法对原子坐标和晶格参数进行全优化，设定原子间的受力收敛标准为小于0.01eV/Å，体系总能量收敛到10⁻⁵eV以内。在自洽场计算中，通过多次迭代使体系的电子密度和能量达到自洽收敛，收敛精度同样设定为10⁻⁵eV。通过严格控制这些收敛精度，确保了计算结果的准确性和稳定性，为后续对本征锐钛矿二氧化钛性质的分析提供了可靠的基础。3.2晶体结构优化结果经过结构优化后，本征锐钛矿二氧化钛的晶体结构参数如表1所示：表1本征锐钛矿二氧化钛优化前后晶体结构参数对比参数优化前优化后实验值晶格参数a（nm）0.37850.38010.3785-0.3810晶格参数c（nm）0.95140.95720.9514-0.9590α=β=γ（°）909090Ti-O键长（nm）-0.194-0.1980.193-0.199键角∠TiO₆（°）-89.5-90.589-91从表中数据可以看出，优化后的晶格参数a和c均略有增加，与实验值相比，处于合理的误差范围内，这验证了本研究采用的计算方法和参数设置的准确性。晶格参数的变化反映了在结构优化过程中，原子间的相互作用发生了调整，以达到体系能量最低的稳定状态。优化后的Ti-O键长和键角也发生了一定的变化，键长在0.194-0.198nm之间，键角∠TiO₆在89.5-90.5°之间，同样与实验值相符。这种结构上的调整对后续性质计算有着重要影响。例如，键长和键角的改变会影响电子云的分布，进而影响材料的电子结构和光学性质。在电子结构方面，可能导致能带结构的变化，包括带隙宽度的改变；在光学性质方面，会影响光与材料的相互作用，如吸收光谱和发射光谱的特征。准确的晶体结构是后续深入研究本征锐钛矿二氧化钛电子、磁性和光学性质的基础，优化后的结构为这些研究提供了可靠的模型。3.3电子结构分析3.3.1能带结构本征锐钛矿二氧化钛的能带结构如图1所示。从能带图中可以清晰地看到，在布里渊区内，价带和导带之间存在明显的禁带。在Г点处，价带顶和导带底之间的垂直距离即为带隙宽度，经计算，本征锐钛矿二氧化钛的带隙宽度为3.15eV，这与实验测量值3.2eV较为接近，进一步验证了计算方法的可靠性。图1本征锐钛矿二氧化钛的能带结构能带的分布呈现出一定的特征。价带主要由O原子的2p轨道电子组成，在能量较低的区域，能带较为平坦，表明电子的有效质量较大，移动性较差。随着能量升高，价带逐渐变得复杂，这是由于O原子2p轨道与Ti原子3d轨道之间的杂化作用逐渐增强。导带则主要由Ti原子的3d轨道电子构成，在导带底附近，能带较为陡峭，意味着电子的有效质量较小，具有较高的迁移率。这种能带结构特征决定了锐钛矿二氧化钛的半导体属性。在热激发或光激发下，价带中的电子可以吸收能量跨越禁带跃迁到导带，从而产生电子-空穴对，实现导电和光催化等过程。价带和导带的主要组成原子轨道也决定了材料的化学活性和物理性质。例如，O原子2p轨道对价带的贡献使得材料具有一定的氧化性，而Ti原子3d轨道对导带的贡献则与材料的光吸收和电荷传输性能密切相关。3.3.2态密度为了进一步深入了解本征锐钛矿二氧化钛的电子结构，计算并绘制了其总态密度（TotalDensityofStates，TDOS）和分波态密度（PartialDensityofStates，PDOS）图，如图2所示。图2本征锐钛矿二氧化钛的总态密度和分波态密度从总态密度图中可以看出，在费米能级（EF）以下，存在明显的价带和导带区域。价带主要分布在-12eV到0eV之间，导带则从3.15eV开始。在费米能级处，态密度为零，这是半导体的典型特征，表明在绝对零度下，价带被电子完全占据，导带为空。分波态密度图清晰地展示了Ti和O原子不同轨道电子对态密度的贡献。在价带中，主要贡献来自O原子的2p轨道，其态密度在-8eV到-2eV之间出现了一个明显的峰，这表明O原子2p轨道电子在该能量范围内具有较高的分布概率。同时，Ti原子的3d轨道在价带中也有一定的贡献，主要分布在靠近费米能级的区域，这是由于Ti原子3d轨道与O原子2p轨道之间存在杂化作用。在导带中，Ti原子的3d轨道起主导作用，态密度在3.5eV到5eV之间出现了明显的峰，说明导带中的电子主要来自Ti原子的3d轨道。O原子的2p轨道在导带中的贡献相对较小，但在某些能量区域仍有一定的分布，这也体现了轨道杂化的影响。费米能级附近态密度的分布主要由Ti和O原子轨道的杂化作用决定。在靠近费米能级处，Ti原子3d轨道与O原子2p轨道的杂化使得态密度发生变化，这种杂化作用不仅影响了材料的电子结构，还对其光学、电学和催化性能产生重要影响。例如，杂化作用可以改变电子的跃迁概率，从而影响材料对光的吸收和发射特性；在催化反应中，杂化作用可以影响反应物分子在材料表面的吸附和反应活性。3.4光学性质分析3.4.1介电函数本征锐钛矿二氧化钛的介电函数实部（ε₁）和虚部（ε₂）随光子能量的变化规律如图3所示。图3本征锐钛矿二氧化钛的介电函数实部和虚部介电函数实部ε₁反映了材料在电场作用下的极化响应能力，它与材料的折射率密切相关。在低光子能量区域（小于3eV），ε₁的值相对稳定，随着光子能量的增加，ε₁逐渐减小。当光子能量接近材料的带隙能量时，ε₁发生明显的变化，这是由于电子从价带跃迁到导带，导致材料的极化性质发生改变。介电函数虚部ε₂描述了材料对光的吸收过程，只有当光子能量满足电子跃迁的能量需求时，ε₂才会出现非零值。从图中可以看出，在光子能量约为3.2eV处，ε₂出现了一个明显的峰，这对应于本征锐钛矿二氧化钛的带隙能量，表明在此能量下，电子可以从价带跃迁到导带，从而吸收光子能量。随着光子能量的进一步增加，ε₂逐渐减小，这是因为高能光子激发的电子跃迁过程逐渐减少。介电函数与材料光学性质密切相关。根据麦克斯韦方程组，介电函数与材料的折射率（n）和消光系数（k）之间存在如下关系：n²-k²=ε₁，2nk=ε₂。通过介电函数可以计算出材料的折射率和消光系数，进而分析材料对光的折射、反射和吸收等光学行为。例如，在可见光范围内，材料的折射率决定了光在材料中的传播速度和方向，而消光系数则决定了光在材料中的衰减程度。3.4.2吸收光谱根据介电函数计算得到的本征锐钛矿二氧化钛的吸收光谱如图4所示。图4本征锐钛矿二氧化钛的吸收光谱从吸收光谱中可以看出，在波长小于387.5nm（对应光子能量大于3.2eV）的紫外光区域，本征锐钛矿二氧化钛具有较强的光吸收能力，这是由于光子能量大于材料的带隙能量，能够激发电子从价带跃迁到导带，产生光生载流子，从而吸收光子能量。在波长大于387.5nm的可见光区域，材料的吸收系数较低，几乎不吸收可见光，这与实验结果一致，也解释了为什么本征锐钛矿二氧化钛只能吸收紫外光，对太阳能的利用率较低。本征材料对光的吸收机制主要是基于电子跃迁。当光子能量与材料的带隙能量匹配时，价带中的电子吸收光子能量跃迁到导带，形成电子-空穴对，从而实现光的吸收。这种光吸收机制与光催化性能密切相关。在光催化过程中，光生载流子（电子-空穴对）可以参与氧化还原反应，将光能转化为化学能，实现对污染物的降解和能源的转化。然而，由于本征锐钛矿二氧化钛只能吸收紫外光，限制了其在实际光催化应用中的效率。因此，如何通过掺杂等改性手段拓展其光响应范围至可见光区域，提高光生载流子的分离效率，成为提高其光催化性能的关键研究方向。四、稀磁半导体锐钛矿二氧化钛的掺杂改性研究4.1不同离子掺杂体系的构建为深入探究离子掺杂对锐钛矿二氧化钛性能的影响，本研究构建了过渡金属离子（如V、Cr、Mn等）、非金属离子（如N、C等）、稀土离子（如La、Ce等）掺杂锐钛矿二氧化钛的超晶胞模型。在构建过渡金属离子掺杂模型时，选取了具有代表性的V、Cr、Mn等元素。以V掺杂为例，首先构建一个2\times2\times2的锐钛矿二氧化钛超晶胞，该超晶胞包含32个Ti原子和64个O原子。然后将其中1个Ti原子替换为V原子，形成V掺杂体系。V原子的替代位置通常选择在原本Ti原子所处的八面体中心位置，这是因为在锐钛矿二氧化钛晶体结构中，Ti原子位于TiO₆八面体的中心，这样的替代方式能够最大程度地影响晶体的局部结构和电子云分布。对于Cr和Mn掺杂，同样采用类似的方法，分别替换超晶胞中的Ti原子，构建相应的掺杂模型。对于非金属离子掺杂，以N和C为例进行模型构建。在构建N掺杂模型时，考虑到N原子的原子半径较小，其可能的掺杂位置有两种情况：一是替代晶格中的O原子，形成N-O键；二是占据晶格间隙位置。本研究分别构建了这两种情况下的N掺杂模型。在替代O原子的模型中，选择一个O原子被N原子替换，以研究N原子替代后对晶体结构和性能的影响。在占据晶格间隙的模型中，通过结构优化确定N原子在晶格间隙中的稳定位置。C掺杂模型的构建方法与N掺杂类似，同样考虑替代O原子和占据晶格间隙两种情况。在构建稀土离子掺杂模型时，选取La和Ce作为代表元素。由于稀土离子的离子半径较大，通常在晶体中倾向于占据较大的间隙位置或替代阳离子。以La掺杂为例，在2\times2\times2的超晶胞中，将La原子替代部分Ti原子，形成La掺杂体系。La原子替代Ti原子后，由于其离子半径与Ti原子不同，会对周围的晶格结构产生较大的影响，导致晶格发生一定程度的畸变。Ce掺杂模型的构建方式与之相同，通过替代超晶胞中的Ti原子来构建Ce掺杂锐钛矿二氧化钛模型。通过以上方法构建的不同离子掺杂体系，为后续研究掺杂对锐钛矿二氧化钛晶体结构、电子结构、磁性和光学性质的影响提供了基础。在构建模型过程中，充分考虑了掺杂原子的位置和替代方式对晶体结构的影响，以确保模型的合理性和有效性。4.2掺杂体系的晶体结构变化不同离子掺杂锐钛矿二氧化钛后，晶体结构参数发生了显著变化，这些变化对晶体结构稳定性产生了重要影响。以过渡金属离子V掺杂为例，通过第一性原理计算优化后的V掺杂锐钛矿二氧化钛超晶胞结构参数与本征材料对比如表2所示：表2V掺杂前后锐钛矿二氧化钛晶体结构参数对比参数本征TiO₂V掺杂TiO₂变化情况晶格参数a（nm）0.38010.3815增大晶格参数c（nm）0.95720.9605增大Ti-O键长（nm）0.194-0.1980.195-0.199增大键角∠TiO₆（°）89.5-90.589.0-90.8略有变化从表中数据可以看出，V掺杂后，晶格参数a和c均有所增大，这是由于V原子的离子半径（0.059nm）略大于Ti原子（0.0605nm），V原子进入TiO₂晶格后，使得晶格发生膨胀。Ti-O键长也有所增加，这是因为V原子与周围O原子的相互作用与Ti原子不同，导致键长发生改变。键角∠TiO₆也有略微的变化，说明V掺杂对晶体的局部几何结构产生了影响。这种晶格参数和键长、键角的变化会改变晶体内部的应力分布，从而影响晶体结构的稳定性。对于非金属离子N掺杂，当N替代晶格中的O原子时，由于N原子的电负性（3.04）与O原子（3.44）不同，会导致局部电荷分布发生变化，进而影响晶体结构。计算结果表明，N掺杂后，晶格参数a和c略有减小，这可能是由于N原子与周围原子形成的化学键强度与O原子有所差异，使得晶格收缩。同时，N-Ti键长与O-Ti键长相比也有所不同，这进一步改变了晶体的局部结构。稀土离子La掺杂后，由于La原子的离子半径（0.1032nm）远大于Ti原子，对晶格结构的影响更为显著。La掺杂使得晶格参数a和c大幅增大，晶体结构发生明显的畸变。这种较大的晶格畸变会导致晶体内部的能量升高，从而影响晶体结构的稳定性。掺杂原子与周围原子的相互作用对晶体结构稳定性起着关键作用。过渡金属离子与周围O原子通过离子键和共价键相互作用，其外层电子结构的变化会影响化学键的强度和方向性，进而影响晶体结构。非金属离子掺杂时，由于其与O原子的电负性差异，会导致电荷转移和键的极性改变，从而影响晶体结构。稀土离子由于其较大的离子半径和特殊的电子结构，与周围原子的相互作用更为复杂，会引起较大的晶格畸变，对晶体结构稳定性产生重要影响。总体而言，掺杂原子与周围原子的相互作用使得晶体结构发生变化，当这种变化导致晶体内部能量升高时，晶体结构的稳定性会降低；反之，若能使晶体内部能量降低，则有助于提高晶体结构的稳定性。4.3掺杂体系的电子结构变化4.3.1能带结构变化掺杂前后锐钛矿二氧化钛的能带结构发生了显著变化，这对材料的光催化和电学性能产生了重要影响。以Cr掺杂为例，图5展示了本征锐钛矿二氧化钛和Cr掺杂锐钛矿二氧化钛的能带结构对比。图5本征与Cr掺杂锐钛矿二氧化钛能带结构对比从图中可以看出，本征锐钛矿二氧化钛的带隙宽度为3.15eV。Cr掺杂后，在价带和导带之间出现了杂质能级，这些杂质能级主要由Cr原子的3d轨道电子贡献。杂质能级的形成使得带隙宽度减小，计算得到Cr掺杂后的带隙宽度为2.85eV。杂质能级的位置靠近价带顶，这使得价带中的电子更容易被激发到杂质能级上，然后再跃迁到导带，从而降低了电子跃迁所需的能量，拓展了材料的光响应范围。杂质能级的分布特征对带隙宽度和载流子迁移率有着重要影响。在Cr掺杂体系中，杂质能级呈现出一定的离散分布，这是由于Cr原子在晶格中的位置和周围环境的对称性所决定的。离散的杂质能级增加了电子跃迁的路径和可能性，使得材料对光的吸收能力增强。然而，杂质能级的存在也可能会影响载流子的迁移率。一方面，杂质能级可以作为载流子的散射中心，增加载流子的散射几率，从而降低载流子迁移率；另一方面，杂质能级的存在也可能会提供额外的载流子传输通道，在一定程度上提高载流子迁移率。具体情况取决于杂质能级的浓度、分布以及与主体能带的耦合程度。对于非金属离子N掺杂，其在TiO₂的价带顶引入了杂质能级，这是由于N原子的2p轨道与O原子的2p轨道杂化所致。N掺杂后，带隙宽度减小，材料对可见光的吸收能力增强。与过渡金属离子掺杂不同，N掺杂引入的杂质能级相对较窄且靠近价带顶，这使得电子从价带跃迁到杂质能级的过程更容易发生，从而实现可见光响应。但同时，由于杂质能级与价带顶的距离较近，电子在跃迁过程中可能会更容易发生复合，影响光生载流子的分离效率。稀土离子La掺杂后，由于其电子结构的特殊性，会在导带底附近引入杂质能级。La的5d和4f轨道电子参与形成杂质能级，使得带隙宽度进一步减小。La掺杂后的能带结构变化不仅影响光吸收性能，还对材料的电学性能产生影响。例如，杂质能级的存在可能会改变材料的电导率和电子迁移率，具体取决于杂质能级与主体能带的相互作用以及载流子在杂质能级和主体能带之间的跃迁过程。4.3.2态密度变化为了深入理解掺杂体系电子结构的变化，计算了不同掺杂体系的态密度。以Mn掺杂锐钛矿二氧化钛为例，图6展示了本征和Mn掺杂锐钛矿二氧化钛的总态密度（TDOS）和分波态密度（PDOS）。图6本征与Mn掺杂锐钛矿二氧化钛态密度对比从总态密度图中可以看出，Mn掺杂后，在费米能级附近出现了明显的峰，这表明Mn掺杂引入了新的电子态。分波态密度图进一步揭示了这些电子态的来源。在Mn掺杂体系中，Mn原子的3d轨道对态密度的贡献显著，在价带和导带之间形成了杂质态。在价带中，Mn的3d轨道与O的2p轨道发生杂化，使得价带态密度发生变化。在导带中，Mn的3d轨道也参与了导带态密度的构成，导致导带态密度分布发生改变。这种电子态分布的变化与材料的电学和磁学性质密切相关。在电学性质方面，杂质态的出现改变了材料的电子填充情况，影响了载流子的浓度和迁移率。例如，在Mn掺杂体系中，杂质态的存在可能会提供额外的载流子，从而改变材料的电导率。在磁学性质方面，Mn原子的3d轨道未成对电子是磁性的主要来源。由于Mn的3d轨道与周围原子轨道的杂化，使得电子云分布发生变化，进而影响磁矩的大小和方向。杂化作用还会导致磁性离子之间的相互作用发生改变，影响材料的磁有序状态。对于C掺杂锐钛矿二氧化钛，C原子的2p轨道在价带顶引入了杂质态，使得价带态密度发生变化。C掺杂后，材料的电学性质可能会发生改变，例如电导率可能会因为杂质态提供的载流子而发生变化。在磁学性质方面，如果C掺杂引入了未成对电子或改变了磁性离子之间的相互作用，材料的磁性也可能会受到影响。稀土离子Ce掺杂后，Ce原子的4f和5d轨道对态密度产生重要贡献。在导带底附近，Ce的4f和5d轨道形成杂质态，改变了导带态密度分布。这种电子态分布的变化对材料的电学性质有显著影响，可能会改变电子的传输特性和电导率。在磁学性质方面，Ce的4f电子具有一定的磁性，其与周围原子的相互作用会影响材料的磁性能，例如磁矩大小和磁相互作用强度。4.4掺杂体系的磁性分析4.4.1磁矩计算与分析不同掺杂体系的磁性来源于磁性离子的未成对电子以及它们之间的相互作用。通过第一性原理计算，得到了不同掺杂体系的总磁矩和原子磁矩，结果如表3所示：表3不同掺杂体系的磁矩计算结果掺杂体系总磁矩（μB）原子磁矩（μB）（主要磁性离子）Mn掺杂TiO₂4.00Mn：4.00Co掺杂TiO₂3.00Co：3.00Ni掺杂TiO₂2.00Ni：2.00以Mn掺杂锐钛矿二氧化钛为例，总磁矩为4.00μB，主要来源于Mn原子的磁矩，Mn原子的磁矩为4.00μB。这是因为Mn原子的3d轨道具有5个未成对电子，这些未成对电子产生了较大的磁矩。在晶体中，Mn原子与周围的O原子形成化学键，其电子云分布受到周围环境的影响。由于Mn-O键的存在，Mn原子的3d轨道电子与O原子的2p轨道电子发生杂化，这种杂化作用对磁矩的大小和方向产生影响。杂化作用可能会导致Mn原子的3d轨道电子云发生畸变，从而改变磁矩的大小。同时，杂化作用还会影响Mn原子与周围磁性离子之间的相互作用，进而影响磁矩的方向。掺杂浓度对磁矩大小和方向也有显著影响。随着Mn掺杂浓度的增加，体系的总磁矩并非简单地线性增加。当掺杂浓度较低时，磁性离子之间的距离较远，相互作用较弱，总磁矩主要由单个磁性离子的磁矩贡献。随着掺杂浓度的增加，磁性离子之间的距离减小，相互作用增强，可能会出现磁矩的耦合现象。当磁性离子之间的耦合为铁磁耦合时，总磁矩会增加；当耦合为反铁磁耦合时，总磁矩可能会减小甚至相互抵消。此外，掺杂浓度的增加还可能会导致晶体结构的变化，进一步影响磁矩的大小和方向。例如，高浓度的Mn掺杂可能会使晶格发生较大的畸变，改变磁性离子之间的相对位置和电子云分布，从而影响磁相互作用和磁矩。不同掺杂原子种类对磁矩的影响也各不相同。Co掺杂锐钛矿二氧化钛体系的总磁矩为3.00μB，Co原子的磁矩为3.00μB。Co原子的3d轨道有7个电子，其中3个未成对电子，其磁矩大小和方向同样受到与周围原子相互作用的影响。与Mn掺杂相比，Co原子的电子结构和磁矩特性不同，导致其与周围原子的相互作用方式也有所差异，从而使Co掺杂体系的磁矩表现出与Mn掺杂体系不同的特点。Ni掺杂体系的总磁矩为2.00μB，Ni原子的磁矩为2.00μB，Ni原子的3d轨道有8个电子，2个未成对电子，其磁矩特性和掺杂后的磁相互作用情况又与Mn和Co掺杂体系有所不同。4.4.2磁相互作用机制探讨在掺杂体系中，磁性原子间的相互作用机制主要包括双交换作用和超交换作用。以Mn掺杂锐钛矿二氧化钛为例，双交换作用是指相邻磁性离子通过中间的非磁性离子（如O原子）进行电子交换，从而实现磁相互作用。在Mn掺杂体系中，Mn-O-Mn键的存在使得Mn原子之间可以通过O原子进行电子交换。当一个Mn原子的电子通过O原子跃迁到另一个Mn原子上时，会伴随着自旋的变化，从而产生磁相互作用。这种双交换作用通常在磁性离子之间距离较近且电子云有一定重叠时较为显著，它对铁磁性的产生起到重要作用。超交换作用则是通过中间非磁性离子的电子轨道重叠来实现磁性离子之间的磁相互作用。在Mn掺杂TiO₂中，Mn原子的3d轨道与O原子的2p轨道发生杂化，形成一定的电子云重叠。相邻Mn原子通过O原子的2p轨道进行磁相互作用，这种超交换作用的强度和方向与Mn-O-Mn键角以及电子云重叠程度密切相关。当Mn-O-Mn键角为180°时，超交换作用最强，且通常表现为反铁磁相互作用；当键角偏离180°时，超交换作用减弱，可能会出现铁磁或其他复杂的磁相互作用。影响磁相互作用强度的因素主要包括掺杂原子的种类、浓度以及晶体结构。不同掺杂原子由于其电子结构和原子半径的差异，与周围原子的相互作用方式和强度不同，从而导致磁相互作用强度不同。例如，Mn、Co、Ni等不同过渡金属离子掺杂时，由于它们的3d电子数和电子云分布不同，与O原子形成的化学键和磁相互作用也不同。掺杂浓度的增加会改变磁性离子之间的距离和相互作用环境，从而影响磁相互作用强度。当掺杂浓度较低时，磁性离子之间距离较远，磁相互作用较弱；随着掺杂浓度的增加，磁性离子之间距离减小，相互作用增强。晶体结构的变化也会对磁相互作用产生重要影响。如前面所述，掺杂会导致晶体结构参数（如晶格常数、键长、键角等）的改变，这些变化会影响磁性离子之间的相对位置和电子云重叠程度，进而影响磁相互作用强度。例如，晶格畸变可能会改变Mn-O-Mn键角，从而改变超交换作用的强度和方向，最终影响材料的磁性。4.5掺杂体系的光学性质变化4.5.1介电函数变化掺杂体系的介电函数是描述材料在电场作用下极化响应和光传播特性的重要物理量。通过第一性原理计算，得到了不同掺杂体系的介电函数，并与本征锐钛矿二氧化钛进行对比。以N掺杂锐钛矿二氧化钛为例，图7展示了本征和N掺杂锐钛矿二氧化钛的介电函数实部（ε₁）和虚部（ε₂）随光子能量的变化曲线。![本征与N掺杂锐钛矿二氧化钛介电函数对比](五、氧空位对掺杂锐钛矿二氧化钛的影响研究5.1含氧空位的掺杂体系构建为深入研究氧空位对掺杂锐钛矿二氧化钛的影响，构建了一系列含有不同浓度氧空位的掺杂体系。以Mn掺杂锐钛矿二氧化钛为例，在2\times2\times2的超晶胞模型基础上，通过移除特定位置的O原子来引入氧空位。在选择氧空位的引入位置时，充分考虑了晶体结构的对称性和氧原子的配位环境。由于锐钛矿二氧化钛中Ti原子位于TiO₆八面体中心，O原子位于八面体顶点，为了研究氧空位对不同配位环境下原子的影响，分别选择了与Mn原子直接配位的O原子以及距离Mn原子较远的O原子位置引入氧空位。具体来说，构建了氧空位浓度分别为1.56%、3.12%和4.69%的掺杂体系，对应在超晶胞中移除1个、2个和3个O原子。对于每个构建的模型，均采用第一性原理计算方法进行结构优化，确保体系的总能量最低且原子间受力平衡。优化过程中，设定原子间的受力收敛标准为小于0.01eV/Å，体系总能量收敛到10⁻⁵eV以内。通过结构优化，得到了稳定的含氧空位掺杂体系结构，为后续研究氧空位对晶体结构、电子结构、磁性和光学性质的影响奠定了基础。通过对优化后结构的分析，验证了模型的合理性和稳定性。从晶体结构参数来看，晶格参数和键长、键角的变化均在合理范围内，表明氧空位的引入并未导致晶体结构的严重畸变。同时，体系的总能量在优化后达到稳定状态，进一步证明了构建的模型能够有效模拟实际的含氧空位掺杂体系，可用于后续的深入研究。5.2氧空位对晶体结构的影响在氧空位存在的情况下，掺杂锐钛矿二氧化钛的晶体结构发生了明显的畸变。以Co掺杂且含有氧空位的体系为例，通过第一性原理计算得到的晶格参数与本征及无空位掺杂体系对比如表4所示：表4Co掺杂含氧空位体系晶格参数变化体系晶格参数a（nm）晶格参数c（nm）本征TiO₂0.38010.9572Co掺杂无空位TiO₂0.38080.9580Co掺杂含1个氧空位TiO₂0.38150.9592Co掺杂含2个氧空位TiO₂0.38220.9605从表中数据可以看出，随着氧空位浓度的增加，晶格参数a和c均逐渐增大。这是因为氧空位的形成导致晶体中局部电荷平衡被打破，为了维持电荷中性，周围原子会发生位置调整。O原子的缺失使得原本与O原子配位的Ti原子和Co原子的键长发生改变，从而引起晶格膨胀。除了晶格参数的变化，原子位置也发生了显著调整。在含氧量空位的掺杂体系中，与氧空位相邻的Ti原子和Co原子会向氧空位方向发生一定程度的位移。以含1个氧空位的Co掺杂体系为例，与氧空位相邻的Ti原子向氧空位方向移动了约0.012nm，Co原子移动了约0.008nm。这种原子位置的调整进一步改变了晶体的局部几何结构，使得键角也发生相应变化。原本理想的∠TiO₆键角在含氧空位体系中从89.5-90.5°变为88.0-91.0°。氧空位对晶体结构稳定性的影响具有双重性。一方面，氧空位的存在破坏了晶体的周期性和对称性，导致晶体内部能量升高，从这个角度来说，氧空位降低了晶体结构的稳定性。另一方面，晶体通过原子位置的调整和晶格的畸变，在一定程度上缓解了由于氧空位引起的局部应力和电荷不平衡，使得体系能量重新达到相对稳定的状态。当氧空位浓度较低时，晶体结构能够通过自身的调整来维持相对稳定；然而，当氧空位浓度过高时，晶体结构的畸变过大，可能导致晶体结构的失稳，甚至发生相变。5.3氧空位对电子结构的影响5.3.1能带结构变化氧空位和掺杂共同作用下，锐钛矿二氧化钛的能带结构发生了显著改变。以Fe掺杂且含有氧空位的体系为例，图8展示了本征、Fe掺杂无空位以及Fe掺杂含氧量空位体系的能带结构对比。图8本征、Fe掺杂无空位及Fe掺杂含氧空位体系能带结构对比从图中可以看出，本征锐钛矿二氧化钛的带隙宽度为3.15eV。Fe掺杂无空位时，在价带和导带之间出现了由Fe原子3d轨道电子贡献的杂质能级，带隙宽度减小至2.90eV。当引入氧空位后，体系的能带结构进一步变化。在价带顶和导带底之间出现了新的缺陷能级，这些缺陷能级主要源于氧空位导致的电子态变化。同时，带隙宽度进一步减小至2.75eV。新形成的杂质能级和缺陷能级对载流子的产生和复合过程产生了重要影响。杂质能级的存在使得价带中的电子更容易被激发到杂质能级上，然后再跃迁到导带，从而降低了电子跃迁所需的能量，增加了光生载流子的产生几率。然而，杂质能级和缺陷能级也可能成为载流子的复合中心。光生载流子在迁移过程中，可能会被杂质能级或缺陷能级捕获，然后发生复合，降低了载流子的寿命和分离效率。具体的影响取决于杂质能级和缺陷能级的浓度、分布以及与主体能带的耦合程度。当杂质能级和缺陷能级浓度较低且分布均匀时，有利于载流子的产生和传输；当浓度过高或分布不均匀时，则可能导致载流子复合加剧，降低光催化性能。5.3.2态密度变化计算含氧量空位掺杂体系的态密度，能够深入了解氧空位对掺杂原子周围电子态分布的影响。以Ni掺杂且含有氧空位的锐钛矿二氧化钛体系为例，图9展示了本征、Ni掺杂无空位以及Ni掺杂含氧量空位体系的总态密度（TDOS）和分波态密度（PDOS）。图9本征、Ni掺杂无空位及Ni掺杂含氧空位体系态密度对比从总态密度图中可以看出，在费米能级附近，本征体系的态密度为零，呈现典型的半导体特征。Ni掺杂无空位时，在费米能级附近出现了由Ni原子3d轨道电子贡献的态密度峰，表明引入了新的电子态。当引入氧空位后，费米能级附近的态密度进一步增加，且态密度分布发生了明显变化。分波态密度图进一步揭示了电子态分布的变化细节。在Ni掺杂含氧量空位体系中，氧空位的存在使得Ni原子3d轨道与周围O原子2p轨道的杂化作用增强。在价带中，原本主要由O原子2p轨道贡献的态密度发生了改变，Ni原子3d轨道与O原子2p轨道的杂化导致态密度峰的位置和强度发生变化。在导带中，由于氧空位的影响，Ni原子3d轨道对导带态密度的贡献也发生了变化，出现了新的电子态分布。这种电子结构变化与材料的导电性和磁性密切相关。在导电性方面，氧空位和掺杂共同作用引入的新电子态改变了材料的电子填充情况，增加了载流子的浓度。新的杂质能级和缺陷能级提供了额外的载流子传输通道，使得材料的电导率可能增加。在磁性方面，氧空位对Ni原子磁矩的影响较为复杂。一方面，氧空位导致的电子态变化可能会改变Ni原子3d轨道未成对电子的自旋状态，从而影响磁矩大小；另一方面，氧空位还会影响Ni原子与周围磁性离子之间的磁相互作用，进而影响材料的磁性。当氧空位增强了磁性离子之间的铁磁相互作用时，材料的磁性增强；当氧空位导致反铁磁相互作用增强时，材料的磁性减弱。5.4氧空位对磁性的影响5.4.1磁矩变化分析通过第一性原理计算，得到了含氧量空位掺杂体系的磁矩，并与无空位体系进行对比，以探究氧空位对磁矩大小和方向的影响。以V掺杂锐钛矿二氧化钛体系为例，计算结果如表5所示：表5V掺杂体系磁矩计算结果体系总磁矩（μB）V原子磁矩（μB）V掺杂无空位TiO₂3.003.00V掺杂含1个氧空位TiO₂3.203.20V掺杂含2个氧空位TiO₂3.403.40从表中数据可以看出，随着氧空位浓度的增加，体系的总磁矩和V原子磁矩均逐渐增大。在无空位的V掺杂体系中，V原子的3d轨道具有3个未成对电子，产生了3.00μB的磁矩。当引入氧空位后，氧空位导致的晶体结构畸变和电子态变化对V原子的磁矩产生了影响。氧空位的存在使得V原子周围的电子云分布发生改变，V原子3d轨道与周围O原子2p轨道的杂化作用增强，从而使V原子的未成对电子自旋状态发生变化，导致磁矩增大。氧空位介导的磁相互作用变化机制较为复杂。在含氧量空位的掺杂体系中，氧空位作为缺陷中心，改变了磁性离子之间的电子传递路径和相互作用强度。以V-O-V键为例，氧空位的出现打破了原本的V-O-V键结构，使得V原子之间的磁相互作用通过氧空位进行重新传递。这种通过氧空位介导的磁相互作用可能会导致磁性离子之间的耦合方式发生改变，从原本的弱磁相互作用转变为较强的铁磁或反铁磁相互作用。当氧空位使得磁性离子之间的电子云重叠程度增加时，有利于形成铁磁耦合，从而增强体系的磁性；反之，当电子云重叠程度减小，可能导致反铁磁耦合增强，使体系磁性减弱。5.4.2磁性稳定性研究为了分析氧空位对掺杂体系磁性稳定性的影响，通过模拟不同温度下的磁性行为进行研究。以Cr掺杂且含有氧空位的锐钛矿二氧化钛体系为例，利用分子动力学方法结合第一性原理计算，模拟了体系在不同温度（100K-500K）下的磁矩变化。模拟结果表明，在低温下（100K-200K），含氧量空位的Cr掺杂体系磁矩相对稳定，与无空位体系相比，磁矩略有增加，这与前面磁矩变化分析的结果一致。随着温度的升高，体系的磁矩逐渐减小。在300K左右，含氧量空位体系的磁矩下降速度明显快于无空位体系。当温度升高到400K以上时，含氧量空位体系的磁矩急剧下降，而无空位体系的磁矩仍保持相对较高的值。这是因为氧空位的存在降低了体系的磁性稳定性。氧空位导致的晶体结构畸变和电子态变化使得磁性离子之间的磁相互作用对温度更为敏感。在高温下，热运动加剧，氧空位周围的原子振动增强，破坏了磁性离子之间的有序排列，导致磁矩下降。而无空位体系由于晶体结构相对完整，磁性离子之间的磁相互作用相对稳定，对温度的耐受性更强。氧空位在调控材料磁学性能方面具有潜在的应用价值。通过控制氧空位的浓度和分布，可以实现对材料磁学性能的有效调控。例如，在某些自旋电子学器件中，需要材料具有特定的磁性转变温度和磁矩大小。通过引入适量的氧空位，可以降低材料的磁性转变温度，使其满足器件的工作温度要求；同时，通过调整氧空位的分布，可以优化磁性离子之间的磁相互作用，增强或减弱材料的磁性，从而实现对材料磁学性能的精确调控，为自旋电子学器件的设计和制备提供了新的思路和方法。5.5氧空位对光学性质的影响5.5.1介电函数变化研究氧空位存在下掺杂体系介电函数的变化规律，对于理解光与材料的相互作用具有重要意义。以N掺杂且含有氧空位的锐钛矿二氧化钛体系为例，图10展示了本征、N掺杂无空位以及N掺杂含氧量空位体系的介电函数实部（ε₁）和虚部（ε₂）随光子能量的变化曲线。图10本征、N掺杂无空位及N掺杂含氧空位体系介电函数对比从介电函数实部ε₁来看，在低光子能量区域（小于3eV），本征体系的ε₁值相对稳定，随着光子能量的增加逐渐减小。N掺杂无空位时，ε₁在低光子能量区域略有增加，这是由于N掺杂引入的杂质能级改变了材料的电子极化特性。当引入氧空位后，ε₁在低光子能量区域进一步增加，且在光子能量接近带隙能量时，ε₁的变化更为显著。这是因为氧空位的存在导致晶体结构和电子结构的变化，增强了材料的极化响应能力。介电函数虚部ε₂描述了材料对光的吸收过程。在本征体系中，ε₂在光子能量约为3.2eV处出现明显的峰，对应于带隙跃迁。N掺杂无空位时，由于杂质能级的存在，在较低光子能量区域出现了新的吸收峰，表明材料对可见光的吸收能力增强。引入氧空位后，ε₂在可见光区域的吸收峰强度进一步增加，且吸收峰向低能量方向移动，即发生了红移。这说明氧空位和N掺杂共同作用，进一步拓展了材料的光吸收范围，增强了对可见光的吸收能力。氧空位对光的折射、反射等特性产生重要影响。根据介电函数与折射率（n）和消光系数（k）的关系：n²-k²=ε₁，2nk=ε₂，氧空位导致的介电函数变化会引起折射率和消光系数的改变。在可见光范围内，折射率的变化会影响光在材料中的传播速度和方向，而消光系数的增加则会导致光在材料中的衰减加剧，从而影响材料的光学透明性和反射率。5.5.2吸收光谱变化分析含氧量空位掺杂体系吸收光谱的特征变化，有助于揭示氧空位对材料光催化活性增强的内在联系。以C掺杂且含有氧空位的锐钛矿二氧化钛体系为例，图11展示了本征、C掺杂无空位以及C掺杂含氧量空位体系的吸收光谱。图11本征、C掺杂无空位及C掺杂含氧空位体系吸收光谱对比从吸收光谱中可以看出，本征锐钛矿二氧化钛在波长小于387.5nm的紫外光区域具有较强的吸收，而在可见光区域吸收较弱。C掺杂无空位时，吸收光谱发生了红移，在可见光区域出现了一定的吸收，这是由于C掺杂引入的杂质能级拓展了光响应范围。当引入氧空位后，吸收光谱进一步红移，且在可见光区域的吸收强度显著增强。在波长为400-600nm的可见光范围内，C掺杂含氧量空位体系的吸收系数比C掺杂无空位体系提高了约2-3倍。氧空位影响材料对光的吸收范围和强度的机制主要包括以下几个方面。首先，氧空位的存在改变了晶体的电子结构，在带隙中引入了缺陷能级，这些缺陷能级使得电子跃迁所需的能量降低，从而使材料能够吸收更低能量的光子，实现吸收光谱的红移。其次，氧空位导致的晶体结构畸变和电子云分布变化，增强了光与材料的相互作用，提高了光吸收效率。此外，氧空位还可能影响材料表面对光的散射和反射，进一步改变光的吸收特性。氧空位与光催化活性增强之间存在着紧密的内在联系。在光催化过程中，光生载流子（电子-空穴对）的产生是关键步骤。氧空位和掺杂共同作用下，材料对光的吸收范围拓展和吸收强度增强，能够产生更多的光生载流子。同时，氧空位作为缺陷中心，可以捕获光生电子，抑制电子-空穴对的复合，延长载流子的寿命，从而提高光催化活性。更多的光生载流子能够参与氧化还原反应，将光能转化为化学能，实现对污染物的降解和能源的转化。因此，通过调控氧空位的浓度和分布，可以有效提高锐钛矿二氧化钛的光催化性能，为光催化应用提供了重要的理论依据和实践指导。六、结论与展望6.1研究成果总结本研究基于第一性原理，深入探究了锐钛矿二氧化钛的掺杂改性及氧空位对其性能的影响，取得了一系列有价值的成果。在晶体结构方面，不同离子掺杂均导致锐钛矿二氧化钛的晶格参数、键长和键角发生变化。过渡金属离子（如V、Cr、Mn等）掺杂因离子半径与Ti的差异，使晶格膨胀，键长增大；非金属离子（如N、C等）掺杂时，由于电负性不同改变电荷分布，晶格可能收缩；稀土离子（如La、Ce等）因大离子半径导致晶格大幅畸变。氧空位的引入同样引起晶格参数增大和原子位置调整，导致晶体结构畸变，且其对结构稳定性的影响具有双重性。电子结构层面，掺杂使能带结构显著改变。过渡金属