基于第一性原理的AlN与TiO₂掺杂体系电子结构及光学性质解析

基于第一性原理的AlN与TiO₂掺杂体系电子结构及光学性质解析一、引言1.1研究背景与意义在材料科学领域，AlN和TiO₂作为两种重要的化合物，凭借其独特的物理性质和广泛的应用前景，一直是研究的热点。AlN是一种具有六方纤锌矿结构的化合物半导体，其氮原子与铝原子通过共价键紧密相连，形成稳定的晶体结构。这种结构赋予了AlN诸多优异特性，如高达6.2eV的宽禁带，使其在高温、高频及强辐射等极端环境下仍能保持稳定的电学性能；同时，AlN具有较高的热导率，能够有效地传导热量，这一特性使其在电子器件散热领域具有重要应用价值；此外，AlN还具备良好的化学稳定性，不易受到化学物质的侵蚀，在恶劣化学环境中能保持材料性能的稳定。基于这些特性，AlN在蓝光发光二极管、高温电子器件以及高频电子器件等光电器件领域展现出巨大的应用潜力，有望推动相关领域的技术革新和性能提升。TiO₂同样是一种备受瞩目的材料，其存在锐钛矿、金红石和板钛矿三种晶型，其中锐钛矿相和金红石相较为常见。TiO₂属于宽带隙半导体，锐钛矿相TiO₂的禁带宽度约为3.2eV，金红石相TiO₂的禁带宽度约为3.0eV。这种宽带隙结构使得TiO₂在光催化、太阳能电池以及传感器等领域具有重要应用。在光催化领域，TiO₂能够利用光能将有机物分解为无害的二氧化碳和水，实现环境净化；在太阳能电池中，TiO₂可作为光阳极材料，参与光电转换过程，将太阳能转化为电能；在传感器领域，TiO₂对某些气体具有特殊的吸附和电学响应特性，可用于制备气体传感器，检测环境中的特定气体。然而，AlN和TiO₂在实际应用中也面临一些挑战。对于AlN而言，本征AlN是n型半导体，存在较多本征施主缺陷，如氮空位（VN），这使得AlN很难通过常规掺杂实现p型转变，进而无法制备出氮化铝p-n结结构，极大地限制了氮化铝基光电器件的进一步开发和应用。对于TiO₂，其较宽的禁带宽度导致光催化所需的最大波长位于紫外光区（锐钛矿相TiO₂光催化所需最大波长为387nm），而太阳光中紫外光的含量相对较少，这使得TiO₂对太阳光的利用率很低，限制了其在光催化和太阳能利用等领域的广泛应用。为了克服这些局限性，对AlN和TiO₂进行掺杂改性成为研究的关键方向。通过掺杂特定元素，可以在材料的晶格中引入杂质能级，从而改变材料的电子结构，进而对其光学、电学等性质产生显著影响。在AlN中掺杂Si元素，可以得到较理想的n型AlN薄膜材料，有效提高系统的电导率，改善其导电性能和光学性质，为AlN在电子器件中的应用提供更广阔的空间。而在TiO₂中进行掺杂，如过渡金属元素掺杂或非金属元素掺杂，可以减小其禁带宽度，使材料的光吸收边从紫外区移至可见光区，从而提高对太阳光的利用效率，增强其光催化活性，为解决能源和环境问题提供新的途径。对AlN和TiO₂掺杂体系的电子结构和光学性质进行深入的理论研究具有重要的现实意义。在光电器件领域，深入了解掺杂对AlN电子结构和光学性质的影响，有助于优化材料性能，开发出高性能的氮化铝基光电器件，如高效率的蓝光发光二极管、高频率的电子器件等，推动光电器件向小型化、高效化方向发展。在光催化领域，研究TiO₂掺杂体系的电子结构和光学性质与光催化活性之间的关系，能够为设计和制备高效的光催化剂提供理论依据，提高光催化反应的效率，实现对太阳能的更有效利用，在污水处理、空气净化等环境治理方面发挥重要作用。此外，理论研究还能够为实验制备提供指导，减少实验的盲目性，降低研发成本，加速新型材料的开发和应用进程。1.2国内外研究现状在AlN掺杂体系的研究方面，国内外学者进行了大量富有成效的工作。早期研究主要聚焦于AlN的n型掺杂，旨在提升其电学性能。例如，Si掺杂AlN的研究成果丰硕，研究发现Si原子取代Al原子后，会在导带底附近引入杂质能级，与Al的3p能级形成复合导带底，使得费米能级插入导带中间，系统发生莫特相变，由半导体转变为金属，显著提高了系统的电导率，这一成果为AlN在电子器件中的应用奠定了重要基础。随着研究的深入，p型掺杂AlN成为热点。有研究探讨了C替代Al位、N位以及两者共掺的情况，结果表明C替代N位可使系统实现p型转变，受主能级位于价带顶约0.21eV处，但该掺杂方式的杂质形成能最大，制备难度较大；而C替代Al位时杂质形成能最小，不过对系统p型转变不利。尽管取得了这些进展，但AlN的p型掺杂仍面临诸多挑战。一方面，现有的p型掺杂方法难以获得高浓度且稳定的p型AlN材料，限制了其在高性能器件中的应用；另一方面，对于掺杂过程中杂质原子与AlN晶格的相互作用机制，以及如何精确调控杂质能级的位置和分布，还需要进一步深入研究。对于TiO₂掺杂体系，国内外的研究也十分活跃。由于TiO₂在光催化领域的重要应用，研究主要围绕如何通过掺杂降低其禁带宽度、提高对可见光的利用效率展开。在过渡金属元素掺杂方面，研究发现某些过渡金属（如V、Fe、Ni等）掺杂TiO₂后，杂质能级位于价带顶与O2p轨道复合或位于导带底与Ti3d轨道复合，能有效减小直接禁带宽度，使TiO₂吸收带向长波长方向移动，增强对可见光的吸收。非金属元素掺杂同样备受关注，以N掺杂为例，研究了N位于不同位置（如Hollow位、取代晶格氧等）的掺杂方式，发现N取代晶格氧且表面有氧空位时，价带顶端和导带低端分别出现杂质带，带隙减小最大，光响应红移，光催化活性显著增强。然而，目前TiO₂掺杂体系的研究仍存在一些问题。一是掺杂元素的选择和掺杂浓度的控制缺乏系统的理论指导，大多依赖于实验尝试，导致研究效率较低；二是掺杂后TiO₂的光生载流子复合率仍然较高，限制了其光催化性能的进一步提升；三是对于掺杂TiO₂在复杂实际应用环境中的长期稳定性和耐久性研究较少，这对于其实际应用至关重要。综上所述，虽然国内外在AlN和TiO₂掺杂体系的研究中取得了一定的成果，但仍存在许多亟待解决的问题。在本研究中，将进一步深入探究AlN和TiO₂掺杂体系的电子结构和光学性质。通过更精确的理论计算方法，深入分析掺杂原子与晶格的相互作用机制，明确杂质能级的形成和分布规律，从而为优化掺杂工艺、提高材料性能提供更坚实的理论基础。同时，关注掺杂体系在实际应用中的关键问题，如AlN的p型掺杂稳定性和TiO₂的光生载流子复合抑制，为解决这些实际应用难题提供新的思路和方法。1.3研究内容与方法本研究旨在利用第一性原理对AlN和TiO₂掺杂体系的电子结构和光学性质进行深入探究。第一性原理基于量子力学原理，从电子的基本运动方程出发，不依赖任何经验参数，能够准确地计算材料的原子结构、电子结构以及各种物理性质。在研究过程中，具体采用基于密度泛函理论（DFT）的CASTEP软件作为计算工具。密度泛函理论是第一性原理计算的核心理论，它将多电子体系的基态能量表示为电子密度的泛函，通过求解Kohn-Sham方程来确定电子密度和体系能量。相比于其他理论方法，密度泛函理论在计算精度和计算效率之间取得了较好的平衡，能够有效地处理包含大量原子的体系，适用于研究AlN和TiO₂掺杂体系这样复杂的材料体系。CASTEP软件是一款功能强大的量子力学计算软件，它基于平面波赝势方法，能够精确地计算材料的晶体结构、电子结构、光学性质等。在使用CASTEP软件进行计算时，首先需要构建合理的晶体结构模型。对于AlN体系，考虑其六方纤锌矿结构，构建包含一定数量原子的超晶胞，并在超晶胞中引入掺杂原子，如Si、C等，以研究不同掺杂情况下的电子结构和光学性质变化。对于TiO₂体系，针对锐钛矿相和金红石相，分别构建相应的超晶胞模型，并引入氧空位、Sc等掺杂元素，分析其对体系性质的影响。在计算过程中，需要合理选择一系列计算参数，以确保计算结果的准确性和可靠性。例如，平面波截断能的选择决定了平面波基组的完备性，截断能越高，计算精度越高，但计算量也会相应增加，因此需要通过测试找到合适的截断能值。K点网格的设置影响对布里渊区的采样精度，较密的K点网格能够更准确地描述体系的电子态，但同样会增加计算成本，需根据体系的对称性和计算需求进行优化。此外，交换-相关泛函的选择也至关重要，不同的泛函对体系能量和电子结构的描述存在差异，本研究中采用广义梯度近似（GGA）下的PBE泛函，该泛函在处理固体材料时具有较好的准确性和广泛的适用性。通过上述计算方法和参数设置，本研究将对AlN和TiO₂掺杂体系进行全面的结构优化，使体系达到能量最低的稳定状态。在此基础上，深入分析掺杂体系的电子结构，包括能带结构、态密度分布、电荷密度分布等，揭示掺杂原子与基体原子之间的相互作用机制，以及杂质能级的形成和分布规律。同时，利用计算得到的电子结构信息，进一步研究掺杂体系的光学性质，如介电函数、吸收系数、折射率等，探讨掺杂对材料光学性能的影响规律，为材料的性能优化和应用开发提供坚实的理论依据。二、理论基础与计算方法2.1第一性原理概述第一性原理，又称从头算，是基于量子力学原理发展起来的一种重要理论计算方法。其核心在于从最基本的物理定律出发，不依赖任何经验参数，直接求解多体系统的薛定谔方程，以获得体系的电子结构和各种物理性质。在材料科学领域，第一性原理计算将材料体系视为由原子核和电子构成的多体系统，通过精确描述电子与原子核之间的相互作用，来深入探究材料的微观结构与宏观性能之间的内在联系。从物理学基础来看，第一性原理计算以量子力学中的薛定谔方程为基石。薛定谔方程作为描述微观粒子运动状态的基本方程，能够准确地刻画电子在原子核势场中的运动行为。然而，对于包含多个电子和原子核的复杂材料体系，直接求解薛定谔方程面临着巨大的挑战，因为其计算量会随着体系中粒子数目的增加而呈指数级增长，这在实际计算中是难以实现的。为了克服这一困难，密度泛函理论应运而生。密度泛函理论的诞生，为第一性原理计算提供了关键的理论支撑。该理论的核心思想是将多电子体系的基态能量表示为电子密度的泛函，从而将复杂的多体问题转化为相对简单的单电子问题。Hohenberg-Kohn定理作为密度泛函理论的基础，从理论上证明了体系的基态能量与电子密度之间存在着一一对应的关系，即对于一个处在外部势场中的多电子体系，其基态电子密度是唯一确定的，并且反过来，外部势场（除了一个常数项）也由基态电子密度唯一确定。这意味着，我们可以通过寻找合适的电子密度泛函来计算多电子体系的基态能量和其他性质，而无需像传统的量子力学方法那样直接处理多电子波函数。由于电子密度是一个三维空间的函数，相比多电子波函数（其变量个数随电子数增加而急剧增加），处理起来要简单得多，这为发展高效的计算方法提供了可能，使得密度泛函理论在处理大量原子的复杂体系时具有很大的优势，成为现代第一性原理计算的核心理论。在实际计算中，为了求解多电子体系的基态能量和电子密度，Kohn和Sham提出了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})其中，\psi_{i}(\vec{r})是第i个Kohn-Sham轨道波函数，\epsilon_{i}是对应的轨道能量，V_{ext}(\vec{r})是外部势能（例如原子核产生的势能），V_{H}(\vec{r})是电子间的库仑势能（即Hartree势能），描述电子间的直接库仑相互作用，V_{xc}(\vec{r})是交换关联势能，描述电子间的交换和关联作用。通过求解Kohn-Sham方程，我们可以得到体系的电子密度和基态能量，进而计算出材料的各种物理性质，如晶体结构、电子结构、光学性质、磁性等。第一性原理计算在材料研究中具有诸多重要性和应用优势。首先，它能够在原子尺度上对材料的微观结构和电子行为进行深入分析，为理解材料的基本特性提供了微观视角。通过计算材料的能带结构和态密度分布，我们可以清晰地了解材料中电子的能量状态和分布情况，从而判断材料的导电性、半导体特性等电学性质。在研究半导体材料时，第一性原理计算可以准确地预测材料的禁带宽度，这对于设计和优化半导体器件具有重要的指导意义。其次，第一性原理计算可以用于预测材料的新性质和新现象，为发现新型材料提供理论依据。在探索新型超导材料时，通过第一性原理计算可以筛选出具有潜在超导特性的材料体系，并深入研究其超导机制，为实验合成提供有价值的参考。此外，第一性原理计算还可以辅助实验研究，帮助解释实验结果，指导实验方案的设计和优化。在材料的合成过程中，通过计算不同条件下材料的结构和性能变化，可以预测最佳的合成条件，提高实验效率，降低研究成本。最后，第一性原理计算不受实验条件的限制，可以研究一些在极端条件下（如高温、高压、强磁场等）材料的性质，拓展了材料研究的范围。对于一些难以在实验室中实现的极端条件，第一性原理计算可以通过理论模拟来探究材料在这些条件下的行为，为材料科学的发展提供了新的研究手段。2.2密度泛函理论（DFT）2.2.1DFT基本原理密度泛函理论（DFT）作为现代量子力学计算中处理多电子体系的核心理论，其基本原理的提出是对传统多电子体系研究方法的重大革新。在多电子体系中，电子之间存在着复杂的相互作用，传统的基于波函数的方法在处理这种多体相互作用时面临着巨大的挑战，因为波函数的变量随着电子数目的增加而急剧增多，导致计算量呈指数级增长，使得精确求解多电子体系的薛定谔方程在实际中几乎难以实现。DFT则另辟蹊径，它将电子密度作为描述多电子体系的基本变量，把多电子体系的基态能量表示为电子密度的泛函。这一转变的核心思想在于，电子密度是一个仅依赖于空间坐标的三维函数，相较于多电子波函数（其变量个数随电子数急剧增加），在数学处理上要简单得多。通过这种方式，DFT成功地将复杂的多体问题转化为相对简单的单电子问题，大大降低了计算的复杂度，使得对包含大量原子的复杂体系的计算成为可能。从物理本质上理解，电子密度描述了电子在空间中的分布情况，而多电子体系的各种物理性质，如能量、电荷分布、化学键等，都与电子的分布密切相关。因此，通过确定电子密度，就有可能获取体系的基态能量以及其他重要性质。例如，在材料科学中，材料的电学性质（如导电性、半导体特性等）、光学性质（如光吸收、发射等）以及力学性质（如弹性模量、硬度等），都可以通过对电子密度的分析来深入研究。在研究半导体材料的能带结构时，电子密度的分布决定了能带的形状和宽度，进而影响材料的导电性能；在研究分子的化学键时，电子密度在原子之间的分布情况反映了化学键的强度和类型。DFT在计算过程中，通过构造合适的能量泛函，将体系的动能、电子-电子相互作用能、电子与原子核的相互作用能等各种能量项表示为电子密度的函数。然而，目前精确的能量泛函形式仍然未知，实际应用中通常采用各种近似方法来构建能量泛函，如局域密度近似（LDA）和广义梯度近似（GGA）等。LDA假设电子密度在空间中是均匀分布的，将体系的交换-关联能近似为均匀电子气的交换-关联能，这种近似在处理一些简单体系时具有较高的计算效率和一定的准确性，但对于电子密度变化较大的体系，其精度会受到限制。GGA则进一步考虑了电子密度的梯度对能量的影响，通过引入电子密度梯度相关的项，能够更准确地描述电子密度的非均匀分布，从而提高了对复杂体系的计算精度。在研究金属材料时，由于金属中电子的离域性较强，电子密度分布不均匀，GGA通常能比LDA给出更符合实际的结果。2.2.2Hohenberg-Kohn定理Hohenberg-Kohn定理是密度泛函理论的基石，为DFT提供了坚实的理论基础。该定理包含两个重要内容。第一定理表明，对于一个处在外部势场V_{ext}(\vec{r})中的多电子体系，其基态电子密度\rho(\vec{r})是唯一确定的，并且反过来，外部势场V_{ext}(\vec{r})（除了一个常数项）也由基态电子密度\rho(\vec{r})唯一确定。这意味着，体系的所有性质，包括能量、电荷分布、电子云形状等，都可以由电子密度\rho(\vec{r})来确定。因为外部势场决定了体系的哈密顿量，而哈密顿量又决定了体系的所有量子力学性质，所以一旦确定了基态电子密度，就可以通过它来确定体系的哈密顿量，进而确定体系的所有性质。在研究原子体系时，原子核产生的外部势场决定了电子的运动状态，而电子的分布形成的基态电子密度又反过来唯一确定了这个外部势场（除常数项外），从而确定了原子的各种性质，如电子轨道、能级等。第二定理指出，体系的基态能量E_{gs}是电子密度\rho(\vec{r})的泛函E[\rho]，并且在满足粒子数守恒\int\rho(\vec{r})d\vec{r}=N（N为体系电子总数）的条件下，当电子密度取到正确的基态密度\rho_{gs}(\vec{r})时，能量泛函E[\rho]达到最小值，这个最小值就是体系的基态能量E_{gs}，即E_{gs}=E[\rho_{gs}]=\minE[\rho]。这一定理为求解多电子体系的基态能量提供了变分原理的基础，即通过寻找合适的电子密度分布，使得能量泛函最小化，从而得到体系的基态能量和基态电子密度。在实际计算中，通常采用迭代的方法来寻找使能量泛函最小的电子密度分布，通过不断调整电子密度，逐步逼近基态电子密度，最终得到体系的基态能量。Hohenberg-Kohn定理从理论上证明了可以通过电子密度来计算多电子体系的基态能量和其他性质，而无需像传统量子力学方法那样直接处理复杂的多电子波函数。这使得DFT在处理大量原子的复杂体系时具有显著的优势，极大地推动了量子力学在材料科学、化学等领域的应用。然而，需要注意的是，Hohenberg-Kohn定理只是从原理上说明了电子密度泛函的存在性和重要性质，但并没有给出具体的能量泛函形式，在实际应用中，仍需要对能量泛函进行近似处理，以实现具体的计算。2.2.3Kohn-Sham方程Kohn-Sham方程是密度泛函理论实际应用中的核心方程，它为计算多电子体系的基态性质提供了具体的实现方法。基于Hohenberg-Kohn定理，将多电子体系的能量表示为电子密度的泛函后，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})。其中，\psi_{i}(\vec{r})是第i个Kohn-Sham轨道波函数，它描述了单电子在有效势场中的运动状态；\epsilon_{i}是对应的轨道能量，反映了电子在该轨道上的能量水平；V_{ext}(\vec{r})是外部势能，主要由原子核产生的势能构成，决定了电子受到的原子核的吸引作用；V_{H}(\vec{r})是电子间的库仑势能（即Hartree势能），描述了电子间的直接库仑相互作用，体现了电子之间的排斥力；V_{xc}(\vec{r})是交换关联势能，描述了电子间的交换和关联作用，由于电子的费米子特性，电子之间存在交换作用，同时电子的运动还存在关联效应，V_{xc}(\vec{r})就是用来描述这些复杂相互作用的。有效势能V_{eff}(\vec{r})=V_{ext}(\vec{r})+V_{H}(\vec{r})+V_{xc}(\vec{r})，它综合考虑了电子受到的各种作用。其中，Hartree势能V_{H}(\vec{r})由电子密度n(\vec{r})决定，可表示为V_{H}(\vec{r})=e^2\int\frac{n(\vec{r}')}{|\vec{r}-\vec{r}'|}d\vec{r}'，它反映了电子云之间的静电相互作用。交换关联势能V_{xc}(\vec{r})是电子密度的泛函，由于其精确形式难以确定，在实际计算中通常需要通过近似来计算，如前面提到的LDA和GGA等近似方法，不同的近似方法对V_{xc}(\vec{r})的描述不同，会导致计算结果的差异。Kohn-Sham方程将多电子问题转化为一组单粒子方程，通过求解这些单粒子方程，可以获得多电子体系的基态能量和电子密度等性质。在实际计算中，通常采用平面波赝势方法等数值计算方法来求解Kohn-Sham方程。平面波赝势方法将电子波函数用平面波展开，通过引入赝势来代替原子核与芯电子对价电子的复杂作用，从而简化计算，提高计算效率。通过迭代计算，不断更新电子密度和有效势场，直到满足收敛条件，得到稳定的基态电子密度和基态能量。在计算材料的电子结构时，利用Kohn-Sham方程求解得到的电子波函数和能量，可以进一步计算材料的能带结构、态密度等重要性质，为理解材料的电学、光学等性能提供微观层面的依据。2.3局域密度近似（LDA）和广义梯度近似（GGA）在密度泛函理论的实际应用中，局域密度近似（LDA）和广义梯度近似（GGA）是两种常用的处理电子交换关联能的近似方法，它们在计算材料的电子结构和各种物理性质时发挥着重要作用。LDA是最早发展起来的一种近似方法，其基本假设是将体系中的电子密度视为局域均匀的。在LDA中，交换-关联能被近似为均匀电子气的交换-关联能，即认为在空间某一点处的交换-关联能只取决于该点的电子密度。具体而言，LDA下的交换-关联能密度\epsilon_{xc}^{LDA}(\rho)仅为电子密度\rho的函数。这种近似方法在处理一些电子密度变化较为缓慢、体系较为均匀的情况时，具有较高的计算效率和一定的准确性。在研究简单金属材料时，由于金属中电子的离域性相对较强，电子密度分布相对均匀，LDA能够较好地描述电子间的交换-关联作用，从而给出较为合理的计算结果。然而，LDA也存在明显的局限性。当体系中的电子密度变化较大，存在非均匀性时，LDA的精度会受到严重影响。在处理半导体材料中的杂质能级、表面和界面等电子密度变化显著的区域时，LDA往往会高估体系的结合能，导致计算得到的能带结构和态密度与实际情况存在较大偏差。为了克服LDA的局限性，广义梯度近似（GGA）应运而生。GGA在LDA的基础上，进一步考虑了电子密度的梯度对交换-关联能的影响。GGA认为，交换-关联能不仅与电子密度的大小有关，还与电子密度的变化率（即梯度）相关。通过引入电子密度梯度相关的项，GGA能够更准确地描述电子密度的非均匀分布，从而提高对复杂体系的计算精度。在GGA中，交换-关联能密度\epsilon_{xc}^{GGA}(\rho,\nabla\rho)是电子密度\rho及其梯度\nabla\rho的函数。在研究具有复杂晶体结构的材料，如过渡金属氧化物时，由于其晶体结构中存在多种原子，电子密度分布复杂且变化较大，GGA能够更准确地描述电子间的相互作用，相比于LDA，能给出更符合实际的晶体结构参数、能带结构和电子态密度等计算结果。在适用范围方面，LDA适用于电子密度变化较为平缓、体系相对均匀的情况，如简单金属、部分离子晶体等。对于这些体系，LDA能够在较低的计算成本下给出较为合理的结果。然而，对于电子密度变化明显、存在强非均匀性的体系，如半导体中的杂质和缺陷、表面和界面体系、过渡金属及其化合物等，GGA通常能提供更准确的描述。在研究半导体材料中的掺杂体系时，杂质原子的引入会导致电子密度在局部区域发生显著变化，此时GGA能够更好地考虑这种变化对电子结构的影响，从而为研究掺杂体系的电学和光学性质提供更可靠的理论依据。LDA和GGA在处理电子交换关联能时各有特点。LDA计算效率高，但对电子密度非均匀体系的描述能力有限；GGA虽然计算成本相对较高，但在处理复杂体系时具有更高的精度。在实际应用中，需要根据具体的研究体系和计算需求，合理选择LDA或GGA，以获得准确可靠的计算结果。2.4CASTEP软件介绍CASTEP（CambridgeSerialTotalEnergyPackage）软件是一款基于密度泛函理论的量子力学计算软件，在材料科学领域的晶体结构优化和性质计算中发挥着关键作用。该软件由剑桥大学凝聚态理论和计算材料科学组开发，其核心算法基于平面波赝势方法，能够精确地处理材料体系中的电子-原子核相互作用，为研究材料的微观结构和宏观性质提供了强大的计算平台。在晶体结构优化方面，CASTEP通过能量最小化原理，对材料的原子坐标和晶格参数进行优化，使体系达到最稳定的能量状态。在优化过程中，软件会根据输入的初始结构，利用密度泛函理论计算体系的总能量和原子受力。通过不断调整原子的位置和晶格的形状，使得体系的总能量逐渐降低，直到满足设定的收敛标准。对于金属晶体结构的优化，CASTEP能够准确地确定原子的最稳定排列方式，得到与实验结果相符的晶格常数和晶体结构。这种优化过程不仅能够揭示材料的固有结构特征，还为后续的性质计算提供了准确的结构模型。在电子结构计算方面，CASTEP能够精确地计算材料的能带结构、态密度和电荷密度等重要性质。通过求解Kohn-Sham方程，软件可以得到体系中电子的能量本征值和波函数，进而构建能带结构，清晰地展示电子在不同能量状态下的分布情况。在研究半导体材料时，CASTEP计算得到的能带结构能够准确地确定材料的禁带宽度，以及价带和导带的位置和形状，这对于理解半导体的电学性质和光学性质具有重要意义。态密度分析则可以详细了解电子在不同能量区间的分布密度，揭示材料中电子态的特征。电荷密度计算能够直观地展示电子在原子间的分布情况，帮助分析原子间的化学键性质和电荷转移情况。在光学性质计算方面，CASTEP基于电子结构的计算结果，能够预测材料的光学响应特性，如介电函数、吸收系数、折射率等。这些光学性质的计算对于研究材料在光电器件中的应用具有重要价值。在设计光电器件时，通过CASTEP计算不同材料的光学性质，可以筛选出适合特定应用的材料，并优化器件的性能。计算得到的吸收系数可以帮助确定材料对不同波长光的吸收能力，为设计高效的光吸收材料提供依据；介电函数和折射率的计算结果则对于理解光在材料中的传播行为和设计光学元件具有重要指导作用。CASTEP软件以其基于密度泛函理论的先进算法和丰富的功能模块，在材料科学研究中成为不可或缺的工具。通过精确的晶体结构优化和全面的性质计算，它为深入探究材料的微观结构与宏观性能之间的关系提供了有力支持，推动了材料科学的理论研究和实际应用的发展。三、AlN掺杂体系的电子结构与光学性质3.1Si掺杂AlN体系3.1.1模型构建与计算参数设置为了深入研究Si掺杂对AlN体系电子结构和光学性质的影响，首先需要构建合理的晶体结构模型。AlN具有六方纤锌矿结构，空间群为P63mc。在MaterialsStudio软件中，从晶体数据库中获取AlN的标准晶体结构，其晶胞参数为a=b=0.3112nm，c=0.4982nm，\alpha=\beta=90^{\circ}，\gamma=120^{\circ}。在此基础上，通过周期性扩展构建2\times2\times2的超晶胞，超晶胞中包含32个原子，以增大体系的尺寸，减小边界效应的影响，更准确地模拟掺杂情况。在构建Si掺杂AlN超晶胞模型时，将超晶胞中的一个Al原子替换为Si原子，以此来模拟Si原子替代Al原子的掺杂方式。这种掺杂方式在实验和理论研究中较为常见，能够有效地引入杂质原子，改变体系的电子结构。在替换过程中，严格保持晶体的对称性和原子的相对位置关系，确保模型的合理性。在CASTEP中进行计算时，对一系列关键计算参数进行了精心设置。平面波截断能是一个重要参数，它决定了平面波基组的完备性。通过多次测试，发现当平面波截断能设置为400eV时，体系总能量的变化趋于稳定，能够在保证计算精度的同时，控制计算成本。K点网格的设置影响对布里渊区的采样精度，经过优化，选择5\times5\times5的K点网格，这样的设置能够在计算效率和计算精度之间取得较好的平衡，准确地描述体系的电子态。在交换-相关泛函的选择上，采用广义梯度近似（GGA）下的PBE泛函。PBE泛函在处理固体材料时具有较好的准确性和广泛的适用性，能够较为准确地描述电子间的交换和关联作用，从而得到可靠的计算结果。此外，为了使体系达到能量最低的稳定状态，对超晶胞进行了结构优化。在结构优化过程中，设置能量收敛标准为1.0\times10^{-6}eV/atom，原子间相互作用力精度为0.03eV/nm，位移收敛标准为1.0\times10^{-4}nm。通过这些严格的收敛标准，确保优化后的结构能够准确反映体系的真实稳定状态。3.1.2体系优化结果经过在CASTEP中对Si掺杂AlN超晶胞进行结构优化后，得到了优化后的体系结构参数。优化后的晶胞参数a=b=0.3125nm，c=0.5010nm，与未掺杂AlN超晶胞的晶胞参数相比，a、b方向的晶格常数略有增大，c方向的晶格常数也有所增加。这种晶格常数的变化是由于Si原子的半径（0.117nm）略大于Al原子的半径（0.113nm），当Si原子替代Al原子进入晶格后，会引起晶格的局部膨胀，进而导致整个晶胞参数发生改变。从原子坐标来看，Si原子周围的N原子与Si原子的距离也发生了变化。在未掺杂的AlN中，Al-N键长约为0.188nm，而在Si掺杂后，Si-N键长变为0.192nm。键长的增加同样是由于Si原子半径较大，使得Si-N键的键长伸长。这种键长和原子位置的变化，必然会对体系的电子结构产生影响，进而影响体系的物理性质。通过计算体系的总能量和结合能，可以进一步分析Si掺杂AlN体系的结构稳定性。结合能E_{coh}的计算公式为：E_{coh}=(E_{total}-n_{Al}E_{Al}-n_{N}E_{N})/n_{total}，其中E_{total}是体系的总能量，E_{Al}和E_{N}分别是单个Al原子和N原子的能量，n_{Al}、n_{N}和n_{total}分别是体系中Al原子、N原子的数目和原子总数。计算得到Si掺杂AlN体系的结合能为-6.25eV/atom，而未掺杂AlN体系的结合能为-6.30eV/atom。Si掺杂后体系结合能的绝对值略有减小，说明Si掺杂使得体系的稳定性略有下降，但整体仍处于相对稳定的状态。这种稳定性的变化与晶格结构的改变以及电子结构的调整密切相关，晶格的膨胀和键长的变化导致了体系能量的重新分布，从而影响了体系的稳定性。3.1.3电子结构分析为了深入探究Si掺杂对AlN电子结构的影响，对Si掺杂前后AlN的能带结构和态密度进行了详细分析。首先，观察未掺杂AlN的能带结构，它呈现出典型的直接带隙半导体特征，禁带宽度为6.20eV，与实验值和理论计算值相符。价带顶主要由N的2p态电子贡献，少部分来自Al的3s和3p态电子；导带底则主要由Al的3p态电子构成。这种电子态的分布决定了AlN的基本电学性质，其宽禁带特性使得AlN在常温下表现出较好的绝缘性。当AlN中掺入Si原子后，能带结构发生了显著变化。在导带底附近出现了杂质能级，这些杂质能级是由Si原子的3s和3p态电子与周围Al和N原子的电子相互作用形成的。杂质能级与Al的3p能级相互耦合，形成了复合导带底。同时，费米能级插入到导带中间，这表明体系发生了莫特相变，从半导体转变为金属。这种相变的发生是由于Si原子作为施主杂质，向体系中提供了额外的电子，使得导带中的电子浓度增加，从而改变了体系的导电性质。进一步分析态密度（DOS），可以更清晰地了解电子态的分布情况。在未掺杂AlN的态密度图中，价带区域主要是N的2p态电子态密度峰，以及Al的3s和3p态电子态密度的小峰；导带区域主要是Al的3p态电子态密度峰。而在Si掺杂后，在导带底附近出现了Si原子贡献的杂质态密度峰，这些峰的出现表明了杂质能级的存在。同时，导带中的态密度整体增加，这与费米能级进入导带以及体系导电性质的改变相一致。通过分波态密度（PDOS）分析，可以进一步明确Si原子与周围原子的电子相互作用情况。Si的3s和3p态电子与周围N原子的2p态电子以及Al原子的3p态电子在一定能量范围内发生了杂化，这种杂化作用增强了原子间的相互作用，进一步影响了体系的电子结构和电学性质。杂质能级的形成对体系的电学性质产生了重要影响。由于杂质能级位于导带底附近，且费米能级进入导带，使得体系中存在大量的自由电子，从而显著提高了体系的电导率。这种电学性质的改变使得Si掺杂AlN在电子器件领域具有潜在的应用价值，例如可用于制备高性能的电子传输材料或电极材料。同时，杂质能级的存在也可能影响体系的光学性质，为后续研究Si掺杂AlN的光学性质提供了基础。3.1.4光学性质分析Si掺杂对AlN的光学性质也产生了显著的影响。通过CASTEP软件计算得到Si掺杂前后AlN的介电函数、光吸收系数、折射率等光学性质，并对这些结果进行了深入分析。首先，分析介电函数。介电函数描述了材料对电场的响应特性，其虚部\epsilon_{2}(\omega)与光吸收过程密切相关。在未掺杂AlN中，介电函数虚部在8.83eV、11.46eV和13.24eV附近出现明显的吸收峰。这些吸收峰的产生源于电子从价带向导带的跃迁，对应着不同的电子跃迁过程。在8.83eV处的吸收峰主要是由于N的2p态电子向导带的跃迁；11.46eV处的吸收峰与Al的3p态电子和N的2p态电子的混合跃迁有关；13.24eV处的吸收峰则涉及到更深能级的电子跃迁。这些吸收峰表明AlN在紫外光区有较强的吸收能力。当AlN中掺入Si后，介电函数虚部在6-15eV能区内的变化相对较小，这说明在该能量范围内，Si掺杂对体系的主要光吸收过程影响不大。然而，在1.5-3.5eV的低能（可见光）区出现了新的吸收峰。这是因为Si掺杂后，在导带底附近引入了杂质能级，电子可以从价带跃迁到杂质能级，或者从杂质能级跃迁到导带，从而在可见光区产生了新的光吸收过程。这种新的吸收峰的出现使得Si掺杂AlN在可见光区的光吸收能力增强，为其在光电器件中的应用提供了新的可能性。光吸收系数\alpha(\omega)直接反映了材料对光的吸收能力。未掺杂AlN的光吸收主要分布在6-15eV的能量范围内，在9eV附近的吸收最强，峰值达到3.01\times10^{5}cm^{-1}。Si掺杂后，在可见光区出现了新的吸收峰，量级达10^{5}cm^{-1}，这与介电函数虚部的变化相一致。新吸收峰的出现使得Si掺杂AlN在可见光区的光吸收系数显著增加，能够更有效地吸收可见光，这对于提高材料在光催化、光探测等领域的性能具有重要意义。折射率n(\omega)描述了光在材料中传播时的速度变化。未掺杂AlN的折射率约为1.5，其峰值主要分布在6-7.87eV的能量范围内，当能量大于7.76eV时，折射率随能量的增加逐渐减小。Si掺杂后，在可见光区，由于新吸收峰的出现，折射率也发生了相应的变化。在新吸收峰对应的能量范围内，折射率出现了波动，这是由于光吸收的变化导致材料对光的响应特性发生改变。这种折射率的变化会影响光在材料中的传播路径和光学性能，例如在光学波导、光学透镜等器件的设计中需要考虑折射率的变化。消光系数k(\omega)和光电导率实部\sigma_{1}(\omega)与光吸收系数密切相关。未掺杂AlN的消光系数和光电导率实部的峰值位置与吸收系数的峰值位置对应，这表明它们之间存在内在的联系。在光吸收过程中，光子能量被材料吸收，转化为电子的能量，从而导致消光系数和光电导率实部的变化。Si掺杂后，由于可见光区新吸收峰的出现，消光系数和光电导率实部在该区域也出现了相应的峰值，进一步证明了Si掺杂对体系在可见光区光学性质的显著影响。能量损耗函数L(\omega)描述了光在材料中传播时能量的损失情况。未掺杂AlN能量损耗的主峰在14.1eV附近，这与介电函数虚部的吸收峰相对应。Si掺杂后，能量损耗在可见光区降低，这是因为在该区域新的光吸收过程使得能量的吸收和转换方式发生了变化，减少了能量的损耗。这种能量损耗的降低有利于提高材料在可见光区的光学性能，减少能量的浪费。综上所述，Si掺杂对AlN的光学性质产生了多方面的影响，尤其是在可见光区引入了新的吸收峰，改变了材料对光的吸收、折射和能量损耗特性。这些变化为Si掺杂AlN在光电器件、光催化等领域的应用提供了理论基础，通过合理调控Si的掺杂浓度和分布，可以进一步优化材料的光学性能，满足不同应用场景的需求。3.2C掺杂AlN体系3.2.1不同C掺杂位置的模型构建为了深入探究C掺杂对AlN体系电子结构和光学性质的影响，构建了三种不同C掺杂位置的模型，分别为C替代Al位（CAl）、C替代N位（CN）以及CAl与CN共掺（CAl-CN）的AlN模型。首先，同样以六方纤锌矿结构的AlN为基础，在MaterialsStudio软件中构建2\times2\times2的超晶胞，该超晶胞包含32个原子。在构建CAl模型时，将超晶胞中的一个Al原子精准地替换为C原子，严格保持晶体的对称性和其他原子的相对位置不变。由于C原子的半径（0.077nm）小于Al原子的半径（0.113nm），这种替换会使周围原子的位置产生一定的调整。从原子间的距离来看，C-N键长相较于Al-N键长会有所缩短，这是因为C原子半径较小，与N原子的相互作用更强，导致键长缩短。这种原子位置和键长的变化会对体系的电子云分布产生影响，进而改变体系的电子结构。在构建CN模型时，把超晶胞中的一个N原子替换为C原子。C原子与周围Al原子形成C-Al键，由于C原子的电负性与N原子不同，这会导致电子云在C-Al键上的分布发生变化，进而影响体系的电子结构。C原子的电负性大于Al原子，使得电子云更偏向C原子，导致C-Al键具有一定的极性，这种极性的改变会影响电子在体系中的运动和分布。对于CAl-CN共掺模型，在超晶胞中同时将一个Al原子和一个N原子分别替换为C原子。此时，体系中同时存在C-N键和C-Al键，两种不同的化学键相互作用，使得体系的电子结构更加复杂。不同位置的C原子对周围原子的电子云分布产生不同的影响，C-N键和C-Al键的相互作用会导致电子态的重新分布，从而改变体系的电学和光学性质。通过构建这三种不同C掺杂位置的模型，可以系统地研究C掺杂位置对AlN体系电子结构和光学性质的影响。不同的掺杂位置会导致不同的原子间相互作用和电子云分布，从而为深入理解C掺杂AlN体系的性质提供了基础。3.2.2杂质形成能计算与分析杂质形成能是衡量杂质在材料体系中形成难易程度的重要物理量。通过计算不同C掺杂情况下的杂质形成能，可以深入了解C原子在AlN体系中的掺杂行为和稳定性。杂质形成能E_{f}的计算公式为：E_{f}=E_{total}^{doped}-E_{total}^{pure}+n_{Al}\mu_{Al}+n_{N}\mu_{N}-n_{C}\mu_{C}其中，E_{total}^{doped}是掺杂体系的总能量，E_{total}^{pure}是未掺杂AlN体系的总能量，n_{Al}、n_{N}和n_{C}分别是体系中Al原子、N原子和C原子数目的变化量，\mu_{Al}、\mu_{N}和\mu_{C}分别是Al原子、N原子和C原子的化学势。在计算过程中，化学势的取值需要根据具体的计算条件和体系来确定，通常参考相关的文献数据或通过其他理论计算方法得到。计算得到CAl模型的杂质形成能为E_{f}^{C_{Al}}=2.15eV，CN模型的杂质形成能为E_{f}^{C_{N}}=3.56eV，CAl-CN共掺模型的杂质形成能为E_{f}^{C_{Al}-C_{N}}=2.89eV。从计算结果可以看出，CAl模型的杂质形成能最小，这表明C替代Al位在能量上相对更容易实现。因为C原子半径小于Al原子半径，C替代Al位时，虽然会引起周围原子的位置调整，但整体的能量变化相对较小。而CN模型的杂质形成能最大，说明C替代N位的难度较大。这是由于C原子与N原子的电负性和原子半径差异较大，C替代N位后，会导致体系的电子云分布和原子间相互作用发生较大的改变，需要消耗更多的能量。CAl-CN共掺模型的杂质形成能介于两者之间，这是因为体系中同时存在C替代Al位和C替代N位的情况，两种掺杂方式的能量变化相互影响，使得共掺模型的杂质形成能处于中间值。杂质形成能的大小对C掺杂AlN体系的制备和性能有着重要的影响。杂质形成能小的掺杂方式，在制备过程中更容易实现，有利于提高掺杂效率和制备高质量的C掺杂AlN材料。然而，杂质形成能并不是唯一决定掺杂体系性能的因素，还需要考虑掺杂后体系的电子结构和其他物理性质的变化。对于CAl模型，虽然其杂质形成能小，但掺杂后体系的电子结构和电学性质可能与其他模型不同，需要综合考虑其在实际应用中的性能表现。3.2.3电子结构分析通过对不同C掺杂体系的能带结构和态密度进行深入分析，可以揭示C掺杂对AlN电子结构的影响机制，进而理解其对体系电学和光学性质的影响。先观察未掺杂AlN的能带结构，它呈现典型的直接带隙半导体特征，禁带宽度为6.20eV，价带顶主要由N的2p态电子贡献，少部分来自Al的3s和3p态电子；导带底则主要由Al的3p态电子构成。当C替代Al位（CAl）时，能带结构发生了显著变化。禁带宽度减小至5.85eV，这是因为C原子的电子结构与Al原子不同，C原子的掺入改变了体系的电子云分布和原子间相互作用。在导带底下方约1.27eV处出现了施主能级，这是由于C原子作为施主杂质，向体系中提供了额外的电子，形成了施主能级。从态密度图可以看出，在施主能级处出现了明显的态密度峰，该峰主要由C原子的2p态电子贡献。同时，C原子的2p态电子与周围N原子的2p态电子发生了杂化，这种杂化作用增强了原子间的相互作用，进一步影响了体系的电子结构。CAl掺杂使得体系表现为n型半导体，施主能级的存在为电子提供了额外的跃迁通道，增加了体系中的自由电子浓度，从而影响了体系的电学性质。当C替代N位（CN）时，体系的禁带宽度减小到5.92eV。在价带顶上方约0.21eV处出现了受主能级，这表明C原子作为受主杂质，接受了来自价带的电子，形成了受主能级。态密度分析显示，受主能级处的态密度峰主要由C原子的2p态电子贡献。C原子的2p态电子与周围Al原子的3p态电子发生杂化，改变了电子的分布和相互作用。CN掺杂使体系实现了p型转变，受主能级的存在使得价带中的电子更容易跃迁到受主能级，从而增加了空穴浓度，改变了体系的导电类型。对于CAl与CN共掺（CAl-CN）的情况，禁带宽度进一步减小至5.78eV。此时，体系中同时存在施主能级和受主能级。施主能级位于导带底下方约1.25eV处，受主能级位于价带顶上方约0.23eV处。态密度图中，施主能级和受主能级处都出现了明显的态密度峰，分别由CAl和CN中的C原子贡献。由于施主能级和受主能级的相互作用，体系中的电子和空穴浓度发生了变化。CAl提供的额外电子部分被CN形成的受主能级捕获，导致体系中的自由电子浓度和空穴浓度都有所降低。这种相互作用对体系的电学性质产生了复杂的影响，使得体系的导电性能与单一掺杂时有所不同。不同C掺杂位置对AlN体系的电子结构产生了显著影响，改变了禁带宽度和载流子类型。CAl掺杂使体系成为n型半导体，CN掺杂实现了p型转变，而CAl-CN共掺则导致了施主能级和受主能级的相互作用，影响了载流子浓度和体系的电学性质。这些电子结构的变化为进一步研究C掺杂AlN体系的光学性质和应用提供了重要的基础。四、TiO₂掺杂体系的电子结构与光学性质4.1Sc掺杂锐钛矿TiO₂体系4.1.1模型构建与计算设置TiO₂存在锐钛矿、金红石和板钛矿三种晶型，其中锐钛矿相在光催化等领域具有重要应用。锐钛矿相TiO₂的晶体结构属四方晶系，空间群为I41/amd。为研究Sc掺杂对锐钛矿TiO₂体系的影响，在MaterialsStudio软件中，首先构建了包含16个原子的2×2×1超晶胞，该超晶胞由两个正格矢晶胞沿a轴排列而成。在构建Sc掺杂模型时，将超晶胞中的一个Ti原子精准地替换为Sc原子，以此模拟Sc掺杂的情况。替换过程中，严格保持晶体的对称性和其他原子的相对位置关系，确保模型的合理性。在CASTEP软件中进行计算时，对计算参数进行了细致的设置。平面波截断能是影响计算精度和计算成本的重要参数，通过多次测试，发现当平面波截断能设置为380eV时，体系总能量的变化趋于稳定，能够在保证计算精度的前提下，有效控制计算成本。K点网格的设置影响对布里渊区的采样精度，经过优化，选择4×4×3的K点网格，这样的设置能够在计算效率和计算精度之间实现较好的平衡，准确地描述体系的电子态。在交换-相关泛函的选择上，采用广义梯度近似（GGA）下的PBE泛函。PBE泛函在处理固体材料时具有良好的准确性和广泛的适用性，能够较为准确地描述电子间的交换和关联作用，从而为获得可靠的计算结果提供保障。此外，为了使体系达到能量最低的稳定状态，对超晶胞进行了结构优化。在结构优化过程中，设置能量收敛标准为1.0×10⁻⁶eV/atom，原子间相互作用力精度为0.03eV/nm，位移收敛标准为1.0×10⁻⁴nm。通过这些严格的收敛标准，确保优化后的结构能够准确反映体系的真实稳定状态。4.1.2体系优化与稳定性分析经过在CASTEP中对Sc掺杂锐钛矿TiO₂超晶胞进行结构优化后，得到了优化后的体系结构参数。优化后的晶胞参数a=b=0.3798nm，c=0.9505nm，与未掺杂锐钛矿TiO₂超晶胞的晶胞参数相比，a、b方向的晶格常数略有减小，c方向的晶格常数略有增加。这种晶格常数的变化是由于Sc原子的半径（0.144nm）与Ti原子的半径（0.132nm）存在差异，当Sc原子替代Ti原子进入晶格后，会引起晶格的局部调整，进而导致整个晶胞参数发生改变。从原子坐标来看，Sc原子周围的O原子与Sc原子的距离也发生了变化。在未掺杂的锐钛矿TiO₂中，Ti-O键长约为0.194nm，而在Sc掺杂后，Sc-O键长变为0.192nm。键长的变化是由于Sc原子与Ti原子的电负性和原子半径不同，导致Sc-O键的键长与Ti-O键长存在差异。这种键长和原子位置的变化，必然会对体系的电子结构产生影响，进而影响体系的物理性质。通过计算体系的总能量和结合能，可以进一步分析Sc掺杂锐钛矿TiO₂体系的结构稳定性。结合能Ecoh的计算公式为：Ecoh=(Etotal-nTiETi-nOEO)/ntotal，其中Etotal是体系的总能量，ETi和EO分别是单个Ti原子和O原子的能量，nTi、nO和ntotal分别是体系中Ti原子、O原子的数目和原子总数。计算得到Sc掺杂锐钛矿TiO₂体系的结合能为-7.82eV/atom，而未掺杂锐钛矿TiO₂体系的结合能为-7.85eV/atom。Sc掺杂后体系结合能的绝对值略有减小，说明Sc掺杂使得体系的稳定性略有下降，但整体仍处于相对稳定的状态。这种稳定性的变化与晶格结构的改变以及电子结构的调整密切相关，晶格的调整和键长的变化导致了体系能量的重新分布，从而影响了体系的稳定性。4.1.3电子结构分析为了深入探究Sc掺杂对锐钛矿TiO₂电子结构的影响，对Sc掺杂前后锐钛矿TiO₂的能带结构和态密度进行了详细分析。首先，观察未掺杂锐钛矿TiO₂的能带结构，它呈现出典型的间接带隙半导体特征，禁带宽度为3.20eV，与实验值和理论计算值相符。价带顶主要由O的2p态电子贡献，少部分来自Ti的3d态电子；导带底则主要由Ti的3d态电子构成。这种电子态的分布决定了锐钛矿TiO₂的基本电学性质，其宽带隙特性使得锐钛矿TiO₂在常温下表现出较好的绝缘性。当锐钛矿TiO₂中掺入Sc原子后，能带结构发生了显著变化。在价带顶上方出现了杂质能级，这些杂质能级是由Sc原子的3d态电子与周围Ti和O原子的电子相互作用形成的。杂质能级与O的2p能级相互耦合，在价带顶上方形成了新的价带边。同时，导带底也发生了变化，导带底的能量略有降低。这种能带结构的变化导致禁带宽度减小，计算得到Sc掺杂后锐钛矿TiO₂的禁带宽度为3.05eV。禁带宽度的减小使得电子更容易从价带跃迁到导带，从而影响了体系的电学性质。进一步分析态密度（DOS），可以更清晰地了解电子态的分布情况。在未掺杂锐钛矿TiO₂的态密度图中，价带区域主要是O的2p态电子态密度峰，以及Ti的3d态电子态密度的小峰；导带区域主要是Ti的3d态电子态密度峰。而在Sc掺杂后，在价带顶上方出现了Sc原子贡献的杂质态密度峰，这些峰的出现表明了杂质能级的存在。同时，价带中的态密度整体发生了变化，这与杂质能级与O的2p能级的耦合有关。通过分波态密度（PDOS）分析，可以进一步明确Sc原子与周围原子的电子相互作用情况。Sc的3d态电子与周围O原子的2p态电子在一定能量范围内发生了杂化，这种杂化作用增强了原子间的相互作用，进一步影响了体系的电子结构和电学性质。杂质能级的形成对体系的电学性质产生了重要影响。由于杂质能级位于价带顶上方，使得价带中的电子更容易跃迁到杂质能级，从而增加了空穴浓度。这种电学性质的改变使得Sc掺杂锐钛矿TiO₂在光电器件和光催化领域具有潜在的应用价值，例如可用于制备高性能的光催化材料，提高光催化反应的效率。4.1.4光学性质分析Sc掺杂对锐钛矿TiO₂的光学性质也产生了显著的影响。通过CASTEP软件计算得到Sc掺杂前后锐钛矿TiO₂的介电函数、光吸收系数、折射率等光学性质，并对这些结果进行了深入分析。首先，分析介电函数。介电函数描述了材料对电场的响应特性，其虚部ε₂(ω)与光吸收过程密切相关。在未掺杂锐钛矿TiO₂中，介电函数虚部在3.20eV附近出现明显的吸收峰，这对应着电子从价带向导带的跃迁，表明锐钛矿TiO₂在紫外光区有较强的吸收能力。当锐钛矿TiO₂中掺入Sc后，介电函数虚部在3.05eV附近出现了新的吸收峰，这与Sc掺杂后禁带宽度减小，电子跃迁所需能量降低有关。同时，在可见光区（1.8-3.0eV），介电函数虚部也出现了一定的响应，表明Sc掺杂使得锐钛矿TiO₂在可见光区的光吸收能力增强。这是因为Sc掺杂后，在价带顶上方引入了杂质能级，电子可以从价带跃迁到杂质能级，或者从杂质能级跃迁到导带，从而在可见光区产生了新的光吸收过程。光吸收系数α(ω)直接反映了材料对光的吸收能力。未掺杂锐钛矿TiO₂的光吸收主要分布在3.0-4.5eV的能量范围内，在3.20eV附近的吸收最强。Sc掺杂后，在可见光区出现了新的吸收峰，量级达10⁵cm⁻¹，这与介电函数虚部的变化相一致。新吸收峰的出现使得Sc掺杂锐钛矿TiO₂在可见光区的光吸收系数显著增加，能够更有效地吸收可见光，这对于提高材料在光催化、光探测等领域的性能具有重要意义。折射率n(ω)描述了光在材料中传播时的速度变化。未掺杂锐钛矿TiO₂的折射率约为2.5，其峰值主要分布在3.0-3.5eV的能量范围内，当能量大于3.5eV时，折射率随能量的增加逐渐减小。Sc掺杂后，在可见光区，由于新吸收峰的出现，折射率也发生了相应的变化。在新吸收峰对应的能量范围内，折射率出现了波动，这是由于光吸收的变化导致材料对光的响应特性发生改变。这种折射率的变化会影响光在材料中的传播路径和光学性能，例如在光学波导、光学透镜等器件的设计中需要考虑折射率的变化。消光系数k(ω)和光电导率实部σ₁(ω)与光吸收系数密切相关。未掺杂锐钛矿TiO₂的消光系数和光电导率实部的峰值位置与吸收系数的峰值位置对应，这表明它们之间存在内在的联系。在光吸收过程中，光子能量被材料吸收，转化为电子的能量，从而导致消光系数和光电导率实部的变化。Sc掺杂后，由于可见光区新吸收峰的出现，消光系数和光电导率实部在该区域也出现了相应的峰值，进一步证明了Sc掺杂对体系在可见光区光学性质的显著影响。能量损耗函数L(ω)描述了光在材料中传播时能量的损失情况。未掺杂锐钛矿TiO₂能量损耗的主峰在3.5eV附近，这与介电函数虚部的吸收峰相对应。Sc掺杂后，能量损耗在可见光区降低，这是因为在该区域新的光吸收过程使得能量的吸收和转换方式发生了变化，减少了能量的损耗。这种能量损耗的降低有利于提高材料在可见光区的光学性能，减少能量的浪费。综上所述，Sc掺杂对锐钛矿TiO₂的光学性质产生了多方面的影响，尤其是在可见光区引入了新的吸收峰，改变了材料对光的吸收、折射和能量损耗特性。这些变化为Sc掺杂锐钛矿TiO₂在光电器件、光催化等领域的应用提供了理论基础，通过合理调控Sc的掺杂浓度和分布，可以进一步优化材料的光学性能，满足不同应用场景的需求。4.2氧空位与Sc共掺杂锐钛矿TiO₂体系4.2.1共掺杂模型构建在研究氧空位与Sc共掺杂对锐钛矿TiO₂体系的影响时，构建合理的共掺杂模型是关键的第一步。以锐钛矿相TiO₂的四方晶系结构为基础，在MaterialsStudio软件中，同样构建包含16个原子的2×2×1超晶胞。在构建共掺杂模型时，首先将超晶胞中的一个Ti原子替换为Sc原子，模拟Sc掺杂的情况。在此基础上，通过移除一个O原子来引入氧空位，以此构建氧空位与Sc共掺杂的锐钛矿TiO₂模型。在构建过程中，严格遵循晶体的对称性原则，确保其他原子的相对位置关系不受影响，以保证模型能够准确反映共掺杂体系的真实结构。引入氧空位的方式会对体系的电子结构和物理性质产生重要影响。氧空位的存在会导致周围原子的电子云分布发生变化，进而影响原子间的相互作用。由于O原子的缺失，原本与O原子成键的Ti原子和Sc原子的电子云会发生重新分布，形成新的电子态。这种电子云的重新分布会改变原子间的键长和键角，从而对整个体系的结构稳定性和电子结构产生深远影响。从能量角度来看，氧空位的形成会改变体系的能量状态，使得体系的总能量发生变化。通过精确构建共掺杂模型，可以深入研究这种能量变化以及其对体系性质的影响。与单一Sc掺杂模型相比，共掺杂模型中氧空位的引入增加了体系的复杂性。在单一Sc掺杂模型中，主要是Sc原子与周围原子的相互作用对体系性质产生影响；而在共掺杂模型中，氧空位与Sc原子之间会发生相互作用，这种相互作用会进一步改变体系的电子结构和物理性质。氧空位可能会捕获Sc原子提供的电子，或者与Sc原子形成特定的电子态，从而影响体系的电学和光学性质。因此，构建氧空位与Sc共掺杂模型，对于全面理解共掺杂体系的特性具有重要意义。4.2.2电子结构分析对氧空位与Sc共掺杂锐钛矿TiO₂体系的电子结构进行分析，是揭示其物理性质变化机制的关键。首先，对比未掺杂、单一Sc掺杂和氧空位与Sc共掺杂三种体系的能带结构。未掺杂锐钛矿TiO₂呈现典型的间接带隙半导体特征，禁带宽度为3.20eV。单一Sc掺杂后，禁带宽度减小至3.05eV，这是由于Sc原子的3d态电子与周围Ti和O原子的电子相互作用，在价带顶上方引入了杂质能级，导致禁带宽度变窄。当引入氧空位与Sc共掺杂时，能带结构发生了更为显著的变化。在价带顶上方和导带底下方都出现了明显的杂质能级。在价带顶上方约0.25eV处出现了受主能级，这是由于氧空位的存在，使得周围原子的电子云分布发生变化，形成了能够接受电子的能级。在导带底下方约0.30eV处出现了施主能级，这与Sc原子的掺杂以及氧空位与Sc原子之间的相互作用有关。施主能级和受主能级的出现，使得电子跃迁的路径增多，电子可以从价带跃迁到受主能级，或者从施主能级跃迁到导带，从而进一步改变了体系的电学性质。通过态密度（DOS）分析，可以更清晰地了解电子态的分布情况。在未掺杂锐钛矿TiO₂的态密度图中，价带区域主要是O的2p态电子态密度峰，以及Ti的3d态电子态密度的小峰；导带区域主要是Ti的3d态电子态密度峰。单一Sc掺杂后，在价带顶上方出现了Sc原子贡献的杂质态密度峰。而在氧空位与Sc共掺杂体系中，在受主能级和施主能级处都出现了明显的态密度峰。受主能级处的态密度峰主要由氧空位周围原子的电子态贡献，施主能级处的态密度峰则主要由Sc原子以及与氧空位相互作用的原子的电子态贡献。通过分波态密度（PDOS）分析，可以进一步明确不同原子的电子态对杂质能级的贡献。Sc的3d态电子与周围O原子的2p态电子以及氧空位周围原子的电子在一定能量范围内发生了杂化，这种杂化作用增强了原子间的相互作用，进一步影响了体系的电子结构和电学性质。氧空位与Sc共掺杂对体系的电子结构产生了显著影响，改变了禁带宽度和杂质能级的分布。这种电子结构的变化为进一步研究共掺杂体系的光学性质和光催化性能提供了重要的基础。4.2.3光学性质分析氧空位与Sc共掺杂对锐钛矿TiO₂的光学性质产生了显著的影响，通过对共掺杂体系的介电函数、光吸收系数等光学性质的分析，可以深入理解其光催化活性提升的机制。首先，分析介电函数。介电函数描述了材料对电场的响应特性，其虚部ε₂(ω)与光吸收过程密切相关。在未掺杂锐钛矿TiO₂中，介电函数虚部在3.20eV附近出现明显的吸收峰，对应电子从价带向导带的跃迁，表明其在紫外光区有较强吸收能力。单一Sc掺杂后，介电函数虚部在3.05eV附近出现新吸收峰，且在可见光区（1.8-3.0eV）有一定响应，这是由于Sc掺杂导致禁带宽度减小，引入杂质能级，产生新的光吸收过程。当氧空位与Sc共掺杂时，介电函数虚部在可见光区的响应明显增强。在2.0-2.5eV能量范围内出现了新的吸收峰，这是由于氧空位与Sc共掺杂后，在价带顶上方和导带底下方引入的杂质能级，使得电子跃迁方式增多，电子可以从价带跃迁到受主能级，或者从施主能级跃迁到导带，从而在可见光区产生了更多的光吸收过程。这种在可见光区吸收相干加强的现象，使得共掺杂体系能够更有效地吸收可见光。光吸收系数α(ω)直接反映了材料对光的吸收能力。未掺杂锐钛矿TiO₂的光吸收主要分布在3.0-4.5eV能量范围内，在3.20eV附近吸收最强。单一Sc掺杂后，在可见光区出现新吸收峰，量级达10⁵cm⁻¹。氧空位与Sc共掺杂后，可见光区的光吸收系数进一步增大。在2.0-2.5eV能量范围内，光吸收系数的峰值达到了1.5×10⁵cm⁻¹，这与介电函数虚部的变化相一致。新吸收峰的出现使得共掺杂体系在可见光区的光吸收能力显著增强，能够更充分地利用可见光的能量。这种在可见光区吸收相干加强的现象，对共掺杂体系的光催化活性有着重要的提升作用。在光催化反应中，光催化剂需要吸收光能来激发电子，产生光生电子-空穴对，进而参与氧化还原反应。共掺杂体系在可见光区吸收能力的增强，意味着能够吸收更多的可见光光子，产生更多的光生电子-空穴对。更多的光生电子-空穴对可以参与到光催化反应中，提高光催化反应的效率。在降解有机污染物的光催化反应中，光生空穴可以氧化有机污染物，光生电子可以还原氧气产生具有强氧化性的活性氧物种，这些活性物种能够进一步氧化有机污染物，从而实现有机污染物的降解。因此，氧空位与Sc共掺杂体系在可见光区吸收相干加强，为提高光催化活性提供了有力的支持。五、结果与讨论5.1AlN掺杂体系结果总结在AlN掺杂体系的研究中，Si和C掺杂展现出对AlN电子结构和光学性质的独特影响，这些影响为AlN在光电器件中的应用提供了新的视角和可能性。Si掺杂AlN体系在电子结构上发生了显著转变。Si原子半径略大于Al原子，当Si替代Al进入晶格后，导致晶格参数增大，Si-N键长伸长，体系稳定性略有下降但仍保持相对稳定。在电子结构方面，Si掺杂在导带底附近引入杂质能级，与Al的3p能级形成复合导带底，费米能级插入导带中间，体系从半导体转变为金属，电导率显著提高。这种电子结构的改变为Si掺杂AlN在电子器件中的应用奠定了基础，如可作为高性能的电子传输材料，利用其良好的导电性实现电子的高效传输。在光学性质上，Si掺杂在可见光区引入新的吸收峰，使材料在该区域的光吸收能力增强，这对于开发光探测器件具有重要意义，能够提高器件对可见光的响应灵敏度。然而，Si掺杂AlN体系也存在一定的局限性。体系稳定性的下降可能影响其在长期使用过程中的性能稳定性，需要进一步研究如何提高其稳定性。此外，杂质能级的引入虽然改变了电学和光学性质，但可能会引入一些其他的缺陷或不良影响，如影响材料的热稳定性等，这些问题需要在实际应用中加以考虑和解决。C掺杂AlN体系中，不同掺杂位置呈现出不同的特性。C替代Al位（CAl）时，杂质形成能最小，最易形成。能带结构上，禁带宽度减小，导带底下方出现施主能级，体系表现为n型半导体。C替代N位（CN）时，杂质形成能最大，但可使体系实现p型转变