SrTiO₃空位缺陷及替位掺杂的第一性原理研究：电子结构与性能调控机制

SrTiO₃空位缺陷及替位掺杂的第一性原理研究：电子结构与性能调控机制一、引言1.1SrTiO₃材料概述1.1.1SrTiO₃的基本结构与特性SrTiO₃，即钛酸锶，是一种重要的钙钛矿型氧化物，具有独特的结构和优异的性能，在材料科学领域备受关注。其晶体结构属于立方晶系，空间群为Pm-3m，具有高度的对称性。在这种结构中，Sr²⁺离子位于立方晶胞的顶点，Ti⁴⁺离子占据晶胞的体心位置，而O²⁻离子则处于晶胞棱边的中点，形成了一个规则的八面体配位结构，其中Ti⁴⁺离子被六个O²⁻离子所包围。这种结构赋予了SrTiO₃许多特殊的物理性质。从晶体学参数来看，其晶格常数a=b=c，在室温下约为3.905Å，这一精确的晶格常数对于维持其内部原子间的相互作用和晶体的稳定性起着关键作用。在电子结构方面，SrTiO₃表现出典型的半导体特性，其禁带宽度约为3.2eV，这一数值使得它在光电器件应用中具有重要的价值。价带顶主要由O2p态电子贡献，而导带底则主要由Ti3d态电子主导。这种电子结构的分布决定了SrTiO₃的光学和电学响应特性。在光学性质上，由于其宽禁带特性，SrTiO₃对可见光具有较高的透过率，使其在透明导电电极和光探测器等光电器件中展现出潜在的应用前景。在电学性质方面，本征的SrTiO₃是一种绝缘体，但其电学性能可以通过掺杂等手段进行精确调控，从而实现从绝缘态到半导体甚至金属态的转变。1.1.2SrTiO₃的应用领域及研究现状SrTiO₃凭借其优异的性能，在多个领域展现出广泛的应用价值和研究潜力。在晶界电容器领域，SrTiO₃基陶瓷由于其高介电常数和低介电损耗的特性，成为制备高性能晶界层电容器的关键材料。通过精细控制制备工艺和优化微观结构，可以显著提高电容器的储能密度和稳定性，满足现代电子设备对小型化、高性能电容器的需求。在氧敏传感器方面，SrTiO₃对氧气具有良好的吸附和反应活性，其电阻值会随着周围氧气浓度的变化而发生显著改变。基于这一特性，SrTiO₃被广泛应用于制备高灵敏度、快速响应的氧敏传感器，用于环境监测、工业过程控制等领域。作为光开关材料，SrTiO₃的光学性能可通过外部电场或光场进行有效调控。在光通信和光信息处理领域，这种可调控的光学特性使其成为实现高速、低能耗光开关的理想选择。在高温超导薄膜衬底方面，SrTiO₃与高温超导材料之间具有良好的晶格匹配和化学兼容性。这使得它成为生长高质量高温超导薄膜的首选衬底材料，对于提高超导薄膜的性能和稳定性，推动高温超导技术的实际应用具有重要意义。随着科技的不断进步，SrTiO₃在超晶格和下一代超大规模集成器件中也展现出巨大的潜在应用价值。其独特的电子结构和物理性质为构建新型量子器件和高性能集成电路提供了新的可能性。目前，对于SrTiO₃的研究涵盖了多个方面。在电子结构研究中，科学家们通过理论计算和实验测量相结合的方法，深入探究其能带结构、态密度以及电子-声子相互作用等微观特性，以揭示其物理本质和内在规律。在半导化掺杂研究中，致力于寻找合适的掺杂元素和优化掺杂工艺，实现对SrTiO₃电学性能的精确调控，以满足不同应用场景的需求。在导电模型研究方面，深入探讨SrTiO₃的导电机制，建立准确的导电模型，为材料的电学性能优化提供理论指导。对于缺陷态和表面/界面态的研究也取得了重要进展。通过先进的表征技术，如扫描隧道显微镜（STM）、X射线光电子能谱（XPS）等，深入研究缺陷和表面/界面态对材料性能的影响机制，为材料的性能改进和应用拓展提供了关键的理论支持。尽管在SrTiO₃的研究方面已经取得了显著成果，但仍存在许多亟待解决的问题和挑战。例如，如何进一步提高SrTiO₃基材料的性能稳定性和可靠性，如何实现对其微观结构和性能的精确调控，以及如何拓展其在新兴领域的应用等。这些问题的解决将为SrTiO₃材料的发展和应用带来新的机遇和突破。1.2空位缺陷与替位掺杂对材料性能的影响1.2.1空位缺陷的形成与作用机制在SrTiO₃晶体中，空位缺陷的形成是一个复杂的物理过程，受到多种因素的影响。从热力学角度来看，温度的升高会增加原子的热振动能量，当原子获得足够的能量时，就有可能脱离其原本的晶格位置，从而形成空位。在晶体生长过程中，由于原子排列的不完美或生长条件的不均匀，也可能导致空位的产生。化学计量比的偏离同样是产生空位缺陷的重要原因。当SrTiO₃中的Sr、Ti和O原子的比例偏离理想的化学计量比时，为了维持晶体的电中性，就会在晶格中形成相应的空位。O空位的形成会对SrTiO₃的电子结构和性能产生显著影响。当O原子离开晶格位置后，会在其周围留下两个未配对的电子，这些电子会在导带中形成一个缺陷能级。这个缺陷能级的出现使得体系的费米能级进入导带中，从而使材料表现出n型掺杂的特性。随着O空位浓度的增加，进入导带中的费米能级位置会逐渐加深，材料的导电性也会随之增强。在一些研究中发现，适量的O空位可以提高SrTiO₃的光催化活性，因为O空位可以作为电子-空穴对的捕获中心，延长载流子的寿命，从而促进光催化反应的进行。Sr空位的引入同样会改变SrTiO₃的电子结构。由于Sr原子的缺失，会导致周围电荷分布的重新调整，使得体系的费米能级进入价带中，表现出p型掺杂的特性。随着Sr空位浓度的增大，费米能级进入价带的位置会更深，材料的电学性能也会发生相应的变化。Ti空位的形成对SrTiO₃的影响更为显著。Ti原子在晶体结构中起着关键的作用，其缺失会导致晶体结构的局部畸变和电子结构的重大改变。Ti空位的引入会使体系的费米能级进入价带，并且在价带顶附近引入缺陷能级。与Sr空位和O空位相比，Ti空位对体系导电性的改善效果更为明显，因为Ti空位所产生的缺陷能级更容易提供载流子，从而增强材料的导电能力。1.2.2替位掺杂的原理及效果替位掺杂是一种通过将杂质原子引入到晶体晶格中，替代原有原子位置的方法，从而改变材料的几何结构和电子结构，实现对材料性能的调控。在SrTiO₃中，不同元素的替位掺杂会产生不同的效果。当Y替位Sr时，由于Y³⁺离子的价态与Sr²⁺离子不同，会在晶体中引入额外的电子。这些电子会进入导带，使得体系Sr₁₋ₓYₓTiO₃（X=0.125，0.25，0.33）的费米能级进入导带中。在导带底出现大量由掺杂原子贡献的自由载流子-电子，这些自由电子的存在显著提高了材料的电导率，从而有效改善了SrTiO₃的导电性能。这种n型掺杂机制在许多电子器件应用中具有重要意义，例如在制备半导体器件时，可以通过Y替位Sr掺杂来调节材料的导电性能，满足不同器件对电学性能的要求。当In或Al替位Ti时，会导致体系的电子结构发生变化，表现出p型掺杂的特性。以SrTi₁₋ₓInₓO₃（SrTi₁₋ₓAlₓO₃）（X=0.125，0.25，0.33）为例，In³⁺或Al³⁺离子的引入会在价带顶附近引入缺陷能级。这些缺陷能级可以捕获价带中的电子，使得费米能级进入价带中，从而增加了价带中的空穴浓度。费米能级处电子态密度不为零，表明材料中存在着一定数量的载流子，这使得体系的导电能力增强。这种p型掺杂对于一些需要空穴传导的器件，如p型半导体器件和某些光电器件，具有重要的应用价值。通过精确控制替位掺杂的元素种类、掺杂浓度和掺杂位置，可以实现对SrTiO₃电学性能的精细调控。不同的掺杂元素和浓度会对材料的能带结构、载流子浓度和迁移率等产生不同程度的影响，从而为材料的性能优化和应用拓展提供了广阔的空间。在实际应用中，可以根据具体的需求选择合适的替位掺杂方案，以实现材料性能的最大化。1.3第一性原理计算方法介绍1.3.1第一性原理计算的基本原理第一性原理计算是一种基于量子力学原理的计算方法，它直接从基本的物理定律出发，无需借助任何经验参数，通过求解多体薛定谔方程来计算体系的电子结构和相关性质。在固体材料中，体系由原子核和电子组成，它们之间存在着复杂的相互作用。第一性原理计算的核心就是精确描述这些相互作用，从而揭示材料的微观结构和物理性质。对于一个包含N个电子和M个原子核的体系，其哈密顿量可以表示为：H=\sum_{i=1}^{N}\left(-\frac{\hbar^2}{2m_e}\nabla_{i}^{2}-\sum_{A=1}^{M}\frac{Z_Ae^2}{r_{iA}}\right)+\sum_{i\ltj}^{N}\frac{e^2}{r_{ij}}+\sum_{A\ltB}^{M}\frac{Z_AZ_Be^2}{R_{AB}}其中，第一项表示电子的动能，第二项表示电子与原子核之间的库仑吸引能，第三项表示电子之间的库仑排斥能，第四项表示原子核之间的库仑排斥能。\hbar是约化普朗克常数，m_e是电子质量，e是电子电荷，Z_A是原子核A的电荷数，r_{iA}是电子i与原子核A之间的距离，r_{ij}是电子i和电子j之间的距离，R_{AB}是原子核A和原子核B之间的距离。求解薛定谔方程H\Psi=E\Psi，其中\Psi是体系的波函数，E是体系的能量。然而，由于多体相互作用的复杂性，直接求解多电子体系的薛定谔方程是非常困难的，甚至是不可能的。为了简化计算，通常采用一些近似方法。其中，最常用的近似方法是密度泛函理论（DFT）。DFT的基本思想是将多电子体系的能量表示为电子密度的泛函，而不是波函数的泛函。通过引入交换-关联泛函来描述电子之间的交换和关联相互作用，从而将多电子问题转化为单电子问题。在DFT框架下，体系的能量可以表示为：E[\rho]=T[\rho]+V_{ion-ele}[\rho]+V_{ele-ele}[\rho]+E_{xc}[\rho]其中，T[\rho]是电子的动能泛函，V_{ion-ele}[\rho]是电子与离子之间的相互作用能泛函，V_{ele-ele}[\rho]是电子之间的库仑相互作用能泛函，E_{xc}[\rho]是交换-关联能泛函。通过求解Kohn-Sham方程，可以得到体系的电子密度和能量。1.3.2常用的第一性原理计算软件与方法在SrTiO₃的研究中，基于密度泛函理论的VASP（ViennaAb-initioSimulationPackage）软件是一款广泛应用的第一性原理计算软件。VASP采用平面波赝势方法，将电子的波函数用平面波基组展开，同时使用赝势来描述离子实与电子之间的相互作用。这种方法在保证计算精度的同时，大大提高了计算效率，使得对复杂体系的计算成为可能。在使用VASP进行计算时，首先需要构建SrTiO₃的晶体结构模型。根据SrTiO₃的晶体结构特点，确定晶胞参数和原子坐标，并将这些信息输入到计算软件中。然后，选择合适的交换-关联泛函。常用的交换-关联泛函包括局域密度近似（LDA）和广义梯度近似（GGA）等。不同的交换-关联泛函对计算结果会产生一定的影响，因此需要根据具体的研究目的和体系特点选择合适的泛函。在计算过程中，还需要设置一系列的计算参数，如截断能、k点网格等。截断能决定了平面波基组的大小，影响计算的精度和效率。k点网格的选择则影响对布里渊区的采样精度，进而影响计算结果的准确性。通过VASP软件，可以进行几何结构优化，以获得体系的最低能量结构。在结构优化过程中，软件会自动调整原子的位置和晶胞的形状，使得体系的总能量最小化。还可以计算体系的电子结构，如能带结构、态密度等。能带结构可以直观地展示电子在晶体中的能量分布情况，而态密度则反映了不同能量状态下电子的分布密度。这些电子结构信息对于理解SrTiO₃的物理性质和电子行为具有重要意义。1.4研究目的与意义1.4.1目的本研究旨在运用第一性原理计算方法，深入且系统地探究SrTiO₃的空位缺陷及替位掺杂对其电子结构和性能的影响规律。通过精确构建不同空位缺陷类型（如O空位、Sr空位、Ti空位）和不同替位掺杂体系（如Y替位Sr、In或Al替位Ti等）的原子模型，利用基于密度泛函理论的计算软件，细致计算和分析体系的几何结构、电子结构（包括能带结构、态密度等）以及相关物理性能。通过对这些计算结果的深入研究，揭示空位缺陷和替位掺杂对SrTiO₃电子结构和性能影响的内在机制，明确不同缺陷和掺杂类型与材料性能之间的定量关系。在此基础上，为通过合理引入空位缺陷和替位掺杂来优化SrTiO₃材料的性能，提供坚实的理论依据和精确的指导，以满足不同应用领域对SrTiO₃材料性能的多样化需求。1.4.2意义对SrTiO₃空位缺陷及替位掺杂的研究具有重要的理论和实际意义。从理论层面来看，深入研究SrTiO₃的空位缺陷和替位掺杂对其电子结构和性能的影响，有助于我们更全面、深入地理解SrTiO₃材料的本征特性和物理本质。通过精确揭示缺陷和掺杂与材料性能之间的内在联系，可以进一步完善和丰富钙钛矿型氧化物材料的理论体系，为后续相关材料的研究提供重要的理论参考和研究思路。这不仅有助于深入探讨材料中的电子-电子相互作用、电子-晶格相互作用等微观物理过程，还能为新型功能材料的设计和开发提供理论指导。在实际应用方面，SrTiO₃作为一种具有广泛应用前景的功能材料，其性能的优化对于拓展其应用领域至关重要。通过本研究，可以为SrTiO₃材料在晶界电容器、氧敏传感器、光开关、高温超导薄膜衬底等领域的应用提供性能优化的理论依据。通过合理调控空位缺陷和替位掺杂，可以显著提高SrTiO₃材料的电学性能、光学性能和催化性能等，从而满足不同应用场景对材料性能的严格要求。这有助于推动SrTiO₃材料在电子、能源、环境等领域的实际应用，促进相关技术的发展和创新，为解决实际工程问题提供有效的材料解决方案。对SrTiO₃的研究还可以为其他钙钛矿型氧化物材料的研究和应用提供借鉴，推动整个功能材料领域的发展。二、理论基础与计算方法2.1密度泛函理论2.1.1密度泛函理论的基本概念密度泛函理论（DensityFunctionalTheory，DFT）是一种在量子力学框架下用于研究多电子体系电子结构的重要理论方法。在多电子体系中，电子之间存在着复杂的相互作用，这使得直接求解多电子体系的薛定谔方程变得极其困难。DFT的核心思想是将多电子体系的基态能量表示为电子密度的泛函，从而将复杂的多电子问题转化为相对简单的单电子问题。Hohenberg-Kohn定理是DFT的基石。第一定理表明，对于一个处在外部势场V_{ext}(r)中的多电子体系，其基态电子密度\rho(r)是唯一确定的，并且反过来，外部势场V_{ext}(r)（除了一个常数项）也由基态电子密度\rho(r)唯一确定。这意味着体系的所有性质都可以由电子密度\rho(r)来确定，因为外部势场决定了体系的哈密顿量，进而决定了体系的所有量子力学性质。第二定理指出，体系的基态能量E_{gs}是电子密度\rho(r)的泛函E[\rho]，并且在满足粒子数守恒\int\rho(r)dr=N（N为体系电子总数）的条件下，当电子密度取到正确的基态密度\rho_{gs}(r)时，能量泛函E[\rho]达到最小值，这个最小值就是体系的基态能量E_{gs}，即E_{gs}=E[\rho_{gs}]=\minE[\rho]。基于Hohenberg-Kohn定理，DFT将多电子体系的能量表示为电子密度的泛函，其一般形式为：E[\rho]=T[\rho]+V_{ion-ele}[\rho]+V_{ele-ele}[\rho]+E_{xc}[\rho]其中，T[\rho]是电子的动能泛函，描述了电子的运动动能；V_{ion-ele}[\rho]是电子与离子之间的相互作用能泛函，体现了电子与原子核之间的库仑吸引作用；V_{ele-ele}[\rho]是电子之间的库仑相互作用能泛函，反映了电子之间的排斥作用；E_{xc}[\rho]是交换-关联能泛函，用于描述电子之间的交换和关联效应。交换能是由于电子的费米子特性，相同自旋的电子不能占据同一量子态而产生的能量；关联能则是由于电子之间的库仑相互作用，使得电子的运动相互关联而产生的能量。在实际计算中，T[\rho]、V_{ion-ele}[\rho]和V_{ele-ele}[\rho]都有较为明确的表达式，而E_{xc}[\rho]的精确形式目前还无法准确给出，通常需要采用近似方法来处理。2.1.2Kohn-Sham方程的推导与应用Kohn-Sham方程是DFT在实际计算中的核心方程，它为计算多电子体系的基态性质提供了一种有效的途径。为了求解多电子体系的能量，Kohn和Sham引入了一个无相互作用的参考体系，其电子密度与真实多电子体系的电子密度相同。通过变分原理，对能量泛函求极值，得到Kohn-Sham方程。首先，将多电子体系的总能量泛函E[\rho]进行分解。假设体系中有N个电子，电子密度为\rho(r)，则总能量可以表示为：E[\rho]=T_{s}[\rho]+\intV_{ext}(r)\rho(r)dr+\frac{1}{2}\int\frac{\rho(r)\rho(r')}{|r-r'|}drdr'+E_{xc}[\rho]其中，T_{s}[\rho]是无相互作用参考体系中电子的动能泛函；\intV_{ext}(r)\rho(r)dr是电子与外部势场的相互作用能；\frac{1}{2}\int\frac{\rho(r)\rho(r')}{|r-r'|}drdr'是电子之间的库仑相互作用能，也称为Hartree能；E_{xc}[\rho]是交换-关联能泛函。对于无相互作用参考体系，其电子的运动可以用一组单粒子波函数\psi_{i}(r)来描述，电子密度\rho(r)可以表示为：\rho(r)=\sum_{i=1}^{N}|\psi_{i}(r)|^{2}根据变分原理，对能量泛函E[\rho]关于单粒子波函数\psi_{i}(r)求变分，并满足波函数的正交归一化条件\int\psi_{i}^{*}(r)\psi_{j}(r)dr=\delta_{ij}，可以得到Kohn-Sham方程：\left[-\frac{\hbar^{2}}{2m}\nabla^{2}+V_{ext}(r)+V_{H}(r)+V_{xc}(r)\right]\psi_{i}(r)=\epsilon_{i}\psi_{i}(r)其中，-\frac{\hbar^{2}}{2m}\nabla^{2}是电子的动能算符；V_{ext}(r)是外部势场；V_{H}(r)是Hartree势，定义为V_{H}(r)=\int\frac{\rho(r')}{|r-r'|}dr'，它描述了电子之间的平均库仑相互作用；V_{xc}(r)是交换-关联势，定义为V_{xc}(r)=\frac{\deltaE_{xc}[\rho]}{\delta\rho(r)}，它包含了电子之间的交换和关联效应；\epsilon_{i}是第i个单粒子态的能量。在SrTiO₃的计算中，应用Kohn-Sham方程来确定体系的电子态和能量。首先，根据SrTiO₃的晶体结构，确定原子的位置和外部势场V_{ext}(r)。然后，通过迭代求解Kohn-Sham方程，得到体系的电子密度\rho(r)和单粒子波函数\psi_{i}(r)。在迭代过程中，不断更新电子密度和势场，直到能量和电子密度收敛。当收敛后，根据得到的电子密度和单粒子波函数，可以计算体系的总能量、能带结构、态密度等性质。通过分析这些性质，可以深入了解SrTiO₃的电子结构和物理特性，为研究空位缺陷和替位掺杂对其性能的影响提供理论基础。2.2平面波赝势方法2.2.1平面波基组的选择与特点在第一性原理计算中，选择合适的基组来展开电子的波函数是非常关键的步骤，平面波基组由于其独特的性质而被广泛应用。平面波基组具有完备性，这意味着任何一个满足一定条件的函数都可以用平面波的线性组合来精确表示。对于电子的波函数，平面波基组能够提供一种通用且精确的描述方式。根据傅里叶变换的原理，在倒易空间中，平面波是一系列具有不同波矢的正弦和余弦函数的组合，这些波矢的取值可以覆盖整个倒易空间。通过调整平面波的系数，就可以构建出与真实电子波函数非常接近的近似波函数。这种完备性使得平面波基组在处理各种不同类型的体系时都具有很强的适应性，无论是简单的原子体系还是复杂的晶体结构，都能够用平面波基组进行有效的描述。平面波基组还具有正交性，不同波矢的平面波之间相互正交。具体来说，对于两个平面波\exp(ik_{1}\cdotr)和\exp(ik_{2}\cdotr)，当k_{1}\neqk_{2}时，它们在整个空间上的积分满足\int\exp(i(k_{1}-k_{2})\cdotr)dr=0。这种正交性为计算带来了很大的便利，在求解Kohn-Sham方程等计算过程中，可以利用平面波的正交性简化数学运算，减少计算的复杂性。在计算电子密度和能量等物理量时，正交性使得不同平面波之间的交叉项为零，从而可以将复杂的积分运算转化为对单个平面波的简单运算，大大提高了计算效率。平面波基组的另一个重要特点是它的平移对称性。由于晶体具有周期性结构，电子的波函数也具有相应的周期性。平面波\exp(ik\cdotr)在晶体的周期性边界条件下，满足\exp(ik\cdot(r+R))=\exp(ik\cdotr)\exp(ik\cdotR)，其中R是晶体的晶格矢量。这意味着平面波在晶体中的平移不变性与晶体的周期性结构相匹配，能够很好地描述晶体中电子的行为。这种平移对称性使得在计算晶体的电子结构时，可以利用晶体的周期性来减少计算量。通过只计算一个晶胞内的电子状态，然后利用平移对称性将结果扩展到整个晶体，从而大大降低了计算的复杂度。在实际计算中，平面波基组的截断能是一个重要的参数。截断能决定了参与计算的平面波的数量和能量范围。当截断能较低时，参与计算的平面波数量较少，计算速度较快，但可能无法准确描述电子波函数的细节，导致计算结果的精度较低。当截断能较高时，参与计算的平面波数量增多，可以更精确地描述电子波函数，但计算量也会显著增加。因此，需要根据具体的计算需求和体系特点，选择合适的截断能。对于SrTiO₃体系，通常需要通过一系列的测试计算，来确定能够保证计算精度且计算效率较高的截断能值。一般来说，随着截断能的增加，计算结果会逐渐收敛到一个稳定的值，当计算结果的变化小于一定的阈值时，就可以认为选择的截断能是合适的。2.2.2赝势的作用与构建赝势是平面波赝势方法中的一个关键概念，它在简化计算过程和提高计算效率方面起着重要作用。在晶体中，原子核和内层电子（芯电子）构成离子实，价电子在离子实产生的势场中运动。由于芯电子的波函数在离子实附近变化非常剧烈，存在许多节点，这使得在计算中直接处理芯电子会导致计算量极大且计算过程复杂。赝势的作用就是通过构建一个假想的势场来代替真实的离子实势场，从而简化计算。赝势的构建基于以下原理：在离子实外部区域，价电子的波函数变化相对平缓，与自由电子的平面波相近。而在离子实内部，由于芯电子的存在，真实的势场使得价电子波函数产生剧烈振荡。赝势的构建目标是在不改变价电子在离子实外部区域的波函数和能量本征值的前提下，使离子实内部的势场变得平缓，从而避免对芯电子波函数的复杂处理。具体来说，赝势通常由两部分组成：一部分是原子核的库仑吸引势，另一部分是对芯电子排斥作用的等效描述。这两部分相互作用使得总的赝势在离子实内部区域变得较为平坦，同时保证在离子实外部区域与真实势场的一致性。在构建赝势时，需要满足一些条件。首先，在芯区之外（通常以一个特定的芯半径r_c为界），赝势应与全电子势完全符合，即等于原子核与芯层电子的总电荷产生的库仑势。这确保了在离子实外部区域，价电子感受到的势场与真实情况一致，从而保证了计算结果的准确性。其次，用赝势方法计算的价电子的赝波函数与对应的全电子方法计算的波函数应具有同样的本征能量，并且在芯区之外两种波函数的形状和幅度都相同。与全电子波函数不同，赝波函数是无径向节点的平缓变化的函数，这使得在计算中可以用相对较少的平面波来展开赝波函数，大大减少了计算量。为了更好地体现芯层电子的作用，各种赝势还需满足许多其他要求，例如赝波函数和全电子波函数的对数微商在芯半径r_c处相等。目前，有多种构建赝势的方法，应用较为广泛的两种从头算赝势是模守恒赝势（norm-conservingpseudopotential）和超软赝势（ultrasoftpseudopotential）。模守恒赝势要求赝波函数和全电子波函数在离子实外部区域具有相同的电荷密度积分，从而保证了赝势在不同化学环境中的可移植性。超软赝势则通过引入投影算子等技术，进一步降低了对平面波基组的要求，使得在较低的截断能下也能获得较高的计算精度。在SrTiO₃的计算中，选择合适的赝势对于计算结果的准确性和计算效率至关重要。不同的赝势构建方法和参数设置会对计算结果产生一定的影响。通过对比不同赝势下的计算结果，选择能够准确描述SrTiO₃体系电子结构且计算效率较高的赝势。在构建SrTiO₃的赝势时，需要考虑Sr、Ti和O原子的电子结构特点，合理确定芯半径、价电子态等参数，以确保赝势能够准确反映离子实与价电子之间的相互作用。2.3计算模型与参数设置2.3.1SrTiO₃晶体模型的构建构建准确的SrTiO₃晶体模型是进行第一性原理计算的基础。SrTiO₃具有立方晶系的钙钛矿结构，空间群为Pm-3m。在构建晶体模型时，首先需要确定晶胞的选取。通常选取一个最小的重复单元作为晶胞，对于SrTiO₃，其晶胞参数在室温下约为a=b=c=3.905Å。在晶胞中，Sr²⁺离子位于立方晶胞的顶点，每个顶点的Sr²⁺离子被8个晶胞共享，因此每个晶胞中Sr²⁺离子的实际数量为8\times\frac{1}{8}=1个。Ti⁴⁺离子占据晶胞的体心位置，每个晶胞中只有1个Ti⁴⁺离子。O²⁻离子处于晶胞棱边的中点，每条棱边的O²⁻离子被4个晶胞共享，因此每个晶胞中O²⁻离子的数量为12\times\frac{1}{4}=3个。这样，一个SrTiO₃晶胞中包含1个Sr²⁺离子、1个Ti⁴⁺离子和3个O²⁻离子，符合其化学计量比。为了准确描述原子在晶胞中的位置，需要确定原子坐标。以晶胞的一个顶点为原点，建立直角坐标系。对于Sr²⁺离子，其坐标为(0,0,0)；Ti⁴⁺离子的坐标为(\frac{1}{2},\frac{1}{2},\frac{1}{2})；O²⁻离子有3个，分别位于棱边中点，其坐标分别为(\frac{1}{2},0,0)、(0,\frac{1}{2},0)和(0,0,\frac{1}{2})。在构建模型时，还需要考虑晶体的对称性和周期性边界条件。由于SrTiO₃晶体具有高度的对称性，利用其对称性可以减少计算量。在计算中，通过设置周期性边界条件，使得晶胞在三维空间中无限重复，从而模拟出宏观晶体的性质。这意味着在计算一个晶胞内的电子结构时，需要考虑晶胞边界上原子与相邻晶胞中原子的相互作用。通过周期性边界条件，将相邻晶胞的原子信息引入到当前晶胞的计算中，保证了计算结果的准确性。为了研究SrTiO₃的空位缺陷和替位掺杂，需要在上述理想晶体模型的基础上进行相应的修改。对于空位缺陷模型，例如构建O空位模型时，只需将晶胞中某个O²⁻离子的坐标删除，即可模拟O原子缺失的情况。对于Sr空位和Ti空位模型，也采用类似的方法。在构建替位掺杂模型时，以Y替位Sr为例，将晶胞中位于顶点的Sr²⁺离子的坐标替换为Y³⁺离子的坐标。同样，对于In或Al替位Ti的模型，将晶胞体心位置的Ti⁴⁺离子坐标替换为In³⁺或Al³⁺离子的坐标。在构建这些缺陷和掺杂模型时，需要注意保持体系的电中性。当引入的杂质原子价态与被替代原子价态不同时，会导致体系电荷分布的变化。为了维持电中性，可以通过调整体系中的电子数三、SrTiO₃的本征结构与电子性质3.1SrTiO₃的晶体结构优化3.1.1优化方法与过程在对SrTiO₃晶体结构进行优化时，本研究采用了共轭梯度法（ConjugateGradient,CG）。共轭梯度法是一种广泛应用于求解线性方程组和优化问题的迭代算法，在晶体结构优化中具有重要的应用价值。其核心思想是通过构建一系列共轭方向，使得在每次迭代中沿着这些方向进行搜索，从而逐步逼近体系的最低能量状态。在SrTiO₃的结构优化中，共轭梯度法通过不断调整原子的位置和晶胞的参数，使得体系的总能量逐渐降低，直至达到最小值。具体优化过程如下：首先，基于SrTiO₃的理想晶体结构模型，确定原子的初始位置和晶胞参数。在SrTiO₃的立方晶系结构中，晶胞参数设定为a=b=c=3.905Å，这是基于其在室温下的典型晶格常数。原子坐标的初始设定为：Sr²⁺离子位于晶胞顶点，坐标为(0,0,0)；Ti⁴⁺离子占据晶胞体心，坐标为(\frac{1}{2},\frac{1}{2},\frac{1}{2})；O²⁻离子处于晶胞棱边中点，坐标分别为(\frac{1}{2},0,0)、(0,\frac{1}{2},0)和(0,0,\frac{1}{2})。这些初始参数为结构优化提供了基础。在优化过程中，以体系的总能量作为优化目标函数。根据共轭梯度法的原理，首先计算体系在当前原子位置和晶胞参数下的梯度。梯度反映了体系能量随原子位置和晶胞参数变化的方向和速率。通过计算梯度，确定一个搜索方向，使得沿着这个方向移动原子和调整晶胞参数能够最大程度地降低体系的能量。在确定搜索方向后，需要确定步长。步长决定了在搜索方向上移动的距离。通过精确计算，确定一个合适的步长，使得在当前搜索方向上，体系能量能够得到最大程度的降低。在每次迭代中，根据计算得到的搜索方向和步长，更新原子的位置和晶胞参数。然后，重新计算体系的能量和梯度，进入下一次迭代。这个过程不断重复，直到体系能量收敛到一个稳定的值，即达到最小值。在收敛过程中，设置了严格的收敛标准，如能量变化小于10^{-5}eV/atom，原子间力小于0.01eV/Å等。只有当体系满足这些收敛标准时，才认为结构优化完成。3.1.2优化结果与分析经过共轭梯度法的优化，得到了稳定的SrTiO₃晶体结构。优化后的晶格常数为a=b=c=3.902Å。与实验值（约3.905Å）相比，相对误差仅为(3.905-3.902)/3.905×100\%\approx0.077\%。这一微小的误差表明，基于共轭梯度法的结构优化计算结果与实验值高度吻合，验证了计算方法的准确性。从键长来看，优化后的Sr-O键长为2.745Å，Ti-O键长为1.958Å。在SrTiO₃的晶体结构中，Sr-O键和Ti-O键的长度对于维持晶体的稳定性和电子结构起着重要作用。这些键长的优化结果与理论预期相符，进一步证明了优化结果的可靠性。对于键角，O-Ti-O键角优化后为180^{\circ}，这与SrTiO₃的理想晶体结构中八面体配位的几何特征一致。在理想的SrTiO₃结构中，Ti⁴⁺离子被六个O²⁻离子以八面体形式包围，O-Ti-O键角应为180^{\circ}。优化后的键角结果表明，计算得到的晶体结构保持了良好的几何对称性。将优化结果与其他理论计算结果进行对比，也显示出较好的一致性。在一些相关研究中，采用不同计算方法得到的SrTiO₃晶格常数和键长等参数与本研究结果相近。这种一致性进一步验证了本研究采用的共轭梯度法优化结构的可靠性。通过对优化结果的分析，可以得出结论：本研究采用的基于共轭梯度法的第一性原理计算方法，能够准确地描述SrTiO₃的晶体结构。优化后的晶体结构参数与实验值和其他理论计算结果高度吻合，为后续对SrTiO₃电子结构和性能的研究提供了坚实的基础。这些准确的结构参数将有助于深入理解SrTiO₃的物理性质和电子行为，为进一步研究空位缺陷和替位掺杂对其性能的影响提供了可靠的前提。3.2SrTiO₃的电子结构计算3.2.1能带结构分析通过第一性原理计算，得到了SrTiO₃的能带结构，如图[X]所示。从能带图中可以清晰地看到，SrTiO₃的能带分布呈现出明显的特征。价带顶位于布里渊区的R点，导带底位于Γ点。这表明电子从价带激发到导带需要跨越一定的能量间隙，即带隙。计算得到的带隙大小为3.10eV。与实验值（约3.2eV）相比，计算结果略小，但相对误差在可接受范围内。这种误差主要是由于计算过程中采用的近似方法，如密度泛函理论中的交换-关联泛函近似等。尽管存在一定误差，但计算得到的带隙大小仍然能够准确反映SrTiO₃的半导体特性。根据能带结构的特征，可以判断SrTiO₃为间接带隙半导体。在间接带隙半导体中，电子从价带顶跃迁到导带底时，需要同时满足能量和动量守恒。这意味着电子的跃迁过程需要声子的参与，因为声子可以提供动量。在SrTiO₃中，由于价带顶和导带底位于不同的k点，电子跃迁时需要吸收或发射声子来满足动量守恒条件。这种间接带隙特性对SrTiO₃的光学和电学性能产生了重要影响。在光学性质方面，间接带隙使得SrTiO₃在吸收光子时，光子能量不仅要满足带隙能量，还需要考虑声子的能量和动量。这导致SrTiO₃的光吸收效率相对较低，光发射过程也相对复杂。在电学性质方面，间接带隙影响了电子和空穴的复合过程，使得载流子的寿命相对较长。这在一些应用中，如光电器件的稳定性和响应速度等方面，具有重要的意义。对导带底和价带顶的电子态进行分析，发现导带底主要由Ti3d态电子贡献。Ti原子的3d轨道在导带底附近形成了一系列的能级，这些能级对电子的激发和传输起着关键作用。价带顶主要由O2p态电子贡献。O原子的2p轨道在价带顶附近形成了密集的能级分布，决定了价带顶的电子性质。这种电子态的分布与SrTiO₃的晶体结构和原子间的相互作用密切相关。在SrTiO₃的晶体结构中，Ti原子与周围的O原子形成了强的化学键，这种化学键的相互作用导致了电子态的分布特征。了解这些电子态的分布对于深入理解SrTiO₃的电子结构和物理性质具有重要意义。3.2.2态密度分析计算得到的SrTiO₃的总态密度（TDOS）和各原子的分波态密度（PDOS）如图[X]所示。通过对态密度的分析，可以深入了解SrTiO₃中电子的分布和能量状态。从总态密度图中可以看出，在费米能级E_F以下，存在着明显的价带，而在E_F以上，为导带。价带和导带之间存在着明显的能隙，这与能带结构分析中得到的带隙结果一致。在价带顶附近，态密度主要由O2p态电子贡献。O原子的2p轨道在价带顶形成了一个尖锐的峰，表明在这个能量范围内，O2p态电子的密度较高。这是由于O原子的电负性较大，吸引电子的能力较强，使得O2p态电子在价带顶占据了主导地位。在导带底附近，态密度主要由Ti3d态电子贡献。Ti原子的3d轨道在导带底形成了一系列的峰，表明Ti3d态电子在导带底的能量状态较为丰富。这是由于Ti原子的3d轨道与周围O原子的2p轨道相互作用，形成了导带底的能级结构。对分波态密度的进一步分析，可以更清晰地了解不同原子轨道对电子态的贡献。在-7eV到-4eV的能量范围内，Sr4d态电子对总态密度有一定的贡献。Sr原子的4d轨道与周围原子的轨道相互作用，在这个能量区间形成了一定的电子态分布。在-2eV到0eV的能量范围内，除了O2p态电子的贡献外，Ti3s和3p态电子也对总态密度有一定的贡献。这些态电子的贡献使得价带顶附近的电子结构更加复杂。在2eV到5eV的能量范围内，除了Ti3d态电子的贡献外，O2s态电子也对总态密度有一定的贡献。这些态电子的相互作用决定了导带底附近的电子结构。通过态密度分析，可以得出结论：SrTiO₃的电子结构主要由O2p态电子和Ti3d态电子主导。O2p态电子对价带顶的贡献和Ti3d态电子对导带底的贡献，决定了SrTiO₃的半导体特性和电子行为。了解这些电子态的分布和贡献，对于理解SrTiO₃的物理性质和性能优化具有重要意义。3.2.3Mulliken布居分析为了深入研究SrTiO₃中原子间的电荷转移和化学键的性质，进行了Mulliken布居分析。Mulliken布居分析是一种常用的方法，用于计算原子的电荷分布和原子间的成键特征。通过该分析，可以得到原子的净电荷、原子间的重叠布居数等重要信息。分析结果表明，在SrTiO₃中，Sr原子带有正电荷，其净电荷约为+1.25。这是因为Sr原子的电负性相对较小，在与周围原子形成化学键时，电子云偏向于电负性较大的原子，使得Sr原子失去部分电子而带有正电荷。Ti原子的净电荷约为+0.78。Ti原子与周围的O原子形成了强的化学键，在这个过程中，Ti原子也失去了部分电子，但由于其与O原子之间存在一定程度的共价键成分，使得Ti原子的净电荷相对较小。O原子带有负电荷，其净电荷约为-0.68。O原子的电负性较大，吸引电子的能力较强，在与Sr和Ti原子形成化学键时，获得了部分电子而带有负电荷。从原子间的重叠布居数来看，Ti-O键的重叠布居数较大，约为0.45。这表明Ti-O键具有较强的共价键成分。在SrTiO₃的晶体结构中，Ti原子与六个O原子形成了八面体配位结构，Ti-O键的共价键性质对晶体的稳定性和电子结构起着关键作用。Sr-O键的重叠布居数相对较小，约为0.18，表明Sr-O键主要表现为离子键性质。Sr原子与O原子之间的离子键作用，使得SrTiO₃晶体具有一定的离子晶体特征。通过Mulliken布居分析，可以清晰地了解SrTiO₃中原子间的电荷转移和化学键的性质。这些信息进一步揭示了SrTiO₃的电子结构和价键特征，为理解其物理性质和性能提供了重要的微观基础。与能带结构和态密度分析结果相结合，可以更全面地认识SrTiO₃的电子结构和物理性质。四、SrTiO₃的空位缺陷研究4.1不同空位缺陷模型的构建4.1.1O空位缺陷模型在构建含O空位的SrTiO₃模型时，基于优化后的SrTiO₃晶体结构，利用超胞技术来确定O空位的位置和浓度。首先，选择合适的超胞尺寸，通常采用2×2×2的超胞，该超胞包含32个原子（8个Sr原子、8个Ti原子和16个O原子）。这种尺寸的超胞既能有效降低空位之间的相互作用，又能保证计算的准确性和可行性。在确定超胞后，通过删除超胞中特定位置的O原子来引入O空位。考虑到SrTiO₃晶体结构的对称性，选择在Ti-O八面体中的不同位置引入O空位。由于Ti-O八面体在晶体中呈规则排列，O原子在八面体中的位置具有一定的对称性。分别选择八面体顶点、棱边中点等不同位置的O原子进行删除，以研究不同位置O空位对体系的影响。在引入O空位时，还需考虑不同的空位浓度。设置空位浓度为0.125、0.25和0.33等不同情况。当空位浓度为0.125时，在超胞中删除2个O原子；当空位浓度为0.25时，删除4个O原子；当空位浓度为0.33时，删除5个O原子。通过这种方式，构建了一系列不同空位浓度和位置的O空位缺陷模型。在构建模型过程中，还需考虑体系的电中性。由于O原子的缺失会导致体系电荷分布的变化，为了维持电中性，需要对体系中的电子数进行相应调整。通过增加体系中的电子数来补偿O原子缺失所带来的电荷变化，确保模型的合理性。4.1.2Sr空位缺陷模型构建含Sr空位的SrTiO₃模型时，同样基于优化后的晶体结构和超胞技术。选用2×2×2的超胞作为基础，在该超胞中，通过删除特定位置的Sr原子来引入Sr空位。由于Sr原子位于晶胞的顶点，在引入Sr空位时，考虑到晶体的对称性，选择不同顶点位置的Sr原子进行删除。分别删除一个、两个或三个顶点位置的Sr原子，以实现不同浓度的Sr空位引入。当删除一个Sr原子时，空位浓度为0.125；删除两个Sr原子时，空位浓度为0.25；删除三个Sr原子时，空位浓度为0.33。在引入Sr空位后，体系的电荷分布会发生改变。为了保持体系的电中性，根据Sr原子的价态（+2价），相应地减少体系中的电子数。通过这种方式，确保模型在引入Sr空位后仍能保持电中性，从而保证模型的合理性和准确性。与O空位缺陷模型类似，Sr空位缺陷模型的构建也需要考虑空位之间的相互作用以及计算的可行性。合适的超胞尺寸和空位浓度设置，能够有效减少空位之间的相互干扰，同时保证计算结果能够准确反映Sr空位对体系的影响。4.1.3Ti空位缺陷模型构建含Ti空位的SrTiO₃模型时，以优化后的晶体结构为基础，采用2×2×2的超胞。在超胞中，通过删除体心位置的Ti原子来引入Ti空位。由于Ti原子在超胞中位于体心位置，且数量相对较少，删除体心位置的Ti原子能够有效引入Ti空位。通过控制删除Ti原子的数量来调整空位浓度。当删除一个Ti原子时，空位浓度为0.125；删除两个Ti原子时，空位浓度为0.25；删除三个Ti原子时，空位浓度为0.33。在引入Ti空位后，由于Ti原子的价态为+4价，为了维持体系的电中性，需要相应地减少体系中的电子数。根据Ti原子的价态和删除的Ti原子数量，精确计算并调整体系的电子数，确保模型的电中性。在构建Ti空位缺陷模型时，同样需要考虑晶体结构的对称性以及空位之间的相互作用。合理选择超胞尺寸和空位浓度，能够避免空位之间的过度相互作用，使得计算结果能够准确反映Ti空位对体系的影响。通过精确构建不同空位缺陷类型和浓度的SrTiO₃模型，为后续研究空位缺陷对材料几何结构、电子结构以及形成能等方面的影响提供了坚实的基础。4.2空位缺陷对几何结构的影响4.2.1晶格常数的变化通过第一性原理计算，得到了不同空位缺陷模型的晶格常数，并与本征SrTiO₃的晶格常数进行对比，以分析空位引入对晶格常数的影响。计算结果表明，O空位的引入会使SrTiO₃的晶格常数发生变化。当O空位浓度为0.125时，晶格常数从本征的3.902Å增大到3.910Å；当O空位浓度增加到0.25时，晶格常数进一步增大至3.918Å；当O空位浓度达到0.33时，晶格常数为3.925Å。这是因为O原子的缺失导致Ti-O八面体结构发生畸变。在SrTiO₃晶体中，Ti-O八面体通过共用顶点的O原子相互连接，形成稳定的晶体结构。当O空位出现时，原本与O原子相连的Ti原子失去了部分配位，导致Ti-O键长发生变化。为了保持结构的稳定性，晶体结构会进行调整，使得晶格常数增大。随着O空位浓度的增加，更多的Ti-O八面体结构受到影响，晶格常数的增大趋势更加明显。Sr空位的引入对晶格常数也有显著影响。当Sr空位浓度为0.125时，晶格常数减小至3.895Å；当Sr空位浓度为0.25时，晶格常数进一步减小到3.888Å；当Sr空位浓度为0.33时，晶格常数为3.880Å。Sr原子位于晶胞的顶点，其主要作用是维持晶体结构的稳定性和电荷平衡。当Sr空位出现时，晶体结构的对称性受到破坏，为了补偿Sr原子的缺失，周围原子会向空位处发生一定程度的偏移。这种原子的偏移导致晶体结构的收缩，从而使得晶格常数减小。随着Sr空位浓度的增加，更多的原子参与到结构调整中，晶格常数的减小趋势更加显著。Ti空位的引入同样会改变晶格常数。当Ti空位浓度为0.125时，晶格常数减小至3.890Å；当Ti空位浓度为0.25时，晶格常数为3.882Å；当Ti空位浓度为0.33时，晶格常数减小到3.875Å。Ti原子在SrTiO₃晶体中处于体心位置，是Ti-O八面体的中心原子，对晶体结构的稳定性起着关键作用。当Ti空位出现时，周围的O原子失去了中心原子的约束，会向空位处靠拢。这种原子的移动导致晶体结构的收缩，使得晶格常数减小。随着Ti空位浓度的增加，更多的O原子受到影响，晶格常数的减小幅度也随之增大。通过对不同空位缺陷模型晶格常数变化的分析，可以看出空位缺陷的引入会显著改变SrTiO₃的晶体结构。这种晶格常数的变化不仅反映了空位对晶体结构的影响，还会进一步影响材料的电子结构和物理性能。4.2.2原子位置的调整空位缺陷的存在不仅会导致晶格常数的变化，还会引起原子位置的调整。通过分析原子间距离和键角的变化，可以深入了解空位周围原子的弛豫情况以及这种弛豫对晶体结构稳定性的影响。在含O空位的SrTiO₃模型中，当O空位出现时，周围的Ti原子和O原子会发生明显的弛豫。以O空位浓度为0.125的模型为例，与O空位相邻的Ti-O键长从本征的1.958Å增大到2.010Å。这是因为O原子的缺失使得Ti原子周围的电子云分布发生改变，Ti原子对剩余O原子的吸引力减弱，导致Ti-O键长增大。与该O空位相邻的O-Ti-O键角也从本征的180°减小到175°。这种键角的变化是由于Ti-O键长的改变以及周围原子的弛豫所导致的。随着O空位浓度的增加，更多的Ti-O键和O-Ti-O键受到影响，键长和键角的变化更加明显。这种原子位置的调整会影响晶体结构的稳定性。键长和键角的变化会导致晶体内部的应力分布发生改变。当应力积累到一定程度时，可能会引发晶体结构的相变或产生其他缺陷，从而影响晶体结构的稳定性。在含Sr空位的SrTiO₃模型中，Sr空位的出现同样会导致周围原子的弛豫。以Sr空位浓度为0.125的模型为例，与Sr空位相邻的Sr-O键长从本征的2.745Å减小到2.700Å。这是因为Sr原子的缺失使得周围O原子对Sr空位的吸引力增强，导致Sr-O键长减小。与该Sr空位相邻的O-Sr-O键角也发生了变化，从本征的90°增大到92°。这种键角的变化是由于Sr-O键长的改变以及周围原子的重新分布所引起的。随着Sr空位浓度的增加，更多的Sr-O键和O-Sr-O键受到影响，键长和键角的变化更加显著。这种原子位置的调整同样会对晶体结构的稳定性产生影响。键长和键角的变化会改变晶体内部的相互作用力，使得晶体结构的稳定性发生变化。如果原子位置的调整导致晶体内部的相互作用力失衡，可能会导致晶体结构的不稳定。在含Ti空位的SrTiO₃模型中，Ti空位的出现会使周围的O原子发生明显的弛豫。以Ti空位浓度为0.125的模型为例，与Ti空位相邻的O-O键长从本征的2.745Å减小到2.680Å。这是因为Ti原子的缺失使得周围O原子之间的排斥力相对增大，导致O-O键长减小。与该Ti空位相邻的O-O-O键角也从本征的180°减小到172°。这种键角的变化是由于O-O键长的改变以及周围O原子的重新排列所导致的。随着Ti空位浓度的增加，更多的O-O键和O-O-O键受到影响，键长和键角的变化更加明显。这种原子位置的调整对晶体结构的稳定性具有重要影响。Ti原子的缺失和周围O原子的弛豫会破坏Ti-O八面体的结构完整性，从而影响晶体结构的稳定性。如果晶体结构的稳定性受到严重破坏，可能会导致材料性能的劣化。通过对不同空位缺陷模型中原子位置调整的分析，可以看出空位缺陷会引起晶体结构的局部畸变，这种畸变会影响晶体结构的稳定性。深入研究原子位置的调整对于理解空位缺陷对SrTiO₃材料性能的影响具有重要意义。4.3空位缺陷对电子结构的影响4.3.1能带结构的改变通过第一性原理计算，得到了含空位缺陷的SrTiO₃的能带结构，并与本征SrTiO₃的能带结构进行对比，以分析空位缺陷对能带结构的影响。计算结果表明，O空位的引入会使SrTiO₃的能带结构发生显著变化。在本征SrTiO₃中，价带顶位于布里渊区的R点，导带底位于Γ点，带隙大小为3.10eV。当引入O空位后，体系的费米能级进入导带中。以O空位浓度为0.125的模型为例，费米能级进入导带的深度约为0.15eV。随着O空位浓度的增加，费米能级进入导带的位置逐渐加深。当O空位浓度为0.25时，费米能级进入导带的深度约为0.25eV；当O空位浓度为0.33时，费米能级进入导带的深度约为0.35eV。这表明O空位的引入使得体系中出现了额外的电子，这些电子进入导带，从而改变了体系的电学性质，使材料表现出n型掺杂的特性。在禁带中也产生了缺陷能级。这些缺陷能级位于导带底附近，随着O空位浓度的增加，缺陷能级的数量也逐渐增多。这些缺陷能级的存在会影响电子的跃迁过程，从而对材料的光学和电学性能产生影响。在光吸收过程中，电子可以从缺陷能级跃迁到导带，从而改变材料的光吸收特性。Sr空位的引入同样会改变SrTiO₃的能带结构。当引入Sr空位后，体系的费米能级进入价带中。以Sr空位浓度为0.125的模型为例，费米能级进入价带的深度约为0.10eV。随着Sr空位浓度的增加，费米能级进入价带的位置逐渐加深。当Sr空位浓度为0.25时，费米能级进入价带的深度约为0.20eV；当Sr空位浓度为0.33时，费米能级进入价带的深度约为0.30eV。这表明Sr空位的引入使得体系中出现了空穴，这些空穴进入价带，从而使材料表现出p型掺杂的特性。在价带顶附近引入了缺陷能级。这些缺陷能级随着Sr空位浓度的增加而增多，它们的存在会影响价带中电子的状态，从而对材料的电学性能产生影响。在电导率方面，这些缺陷能级可以提供额外的载流子，从而改变材料的导电性能。Ti空位的引入对SrTiO₃的能带结构影响更为显著。当引入Ti空位后，体系的费米能级进入价带中，且在价带顶附近引入了较多的缺陷能级。以Ti空位浓度为0.125的模型为例，费米能级进入价带的深度约为0.18eV。随着Ti空位浓度的增加，费米能级进入价带的位置逐渐加深，缺陷能级的数量也明显增多。当Ti空位浓度为0.25时，费米能级进入价带的深度约为0.30eV；当Ti空位浓度为0.33时，费米能级进入价带的深度约为0.40eV。与Sr空位和O空位相比，Ti空位对体系导电性的改善效果更为明显。这是因为Ti空位所产生的缺陷能级更容易提供载流子，从而增强了材料的导电能力。通过对不同空位缺陷模型能带结构变化的分析，可以看出空位缺陷会显著改变SrTiO₃的电子结构，从而对材料的电学和光学性能产生重要影响。4.3.2态密度的变化分析空位缺陷引起的态密度变化，有助于深入了解空位缺陷对SrTiO₃电子结构的影响。通过计算含空位缺陷的SrTiO₃的总态密度（TDOS）和各原子的分波态密度（PDOS），并与本征SrTiO₃的态密度进行对比，观察在缺陷能级处态密度的分布情况以及不同原子轨道对缺陷态密度的贡献。在含O空位的SrTiO₃中，O空位的引入导致在导带底附近出现了新的态密度峰。以O空位浓度为0.125的模型为例，在导带底约0.1-0.3eV的能量范围内，出现了明显的态密度峰。这些新的态密度峰主要由Ti3d态电子贡献。由于O空位的存在，使得Ti原子周围的电子云分布发生改变，Ti3d态电子的能量状态发生变化，从而在导带底附近产生了新的态密度峰。随着O空位浓度的增加，这些新的态密度峰的强度逐渐增强，表明更多的电子占据了这些缺陷能级。在价带顶附近，态密度也发生了一些变化。O2p态电子的态密度略有下降，这是因为O原子的缺失导致O2p态电子的数量减少。在含Sr空位的SrTiO₃中，Sr空位的引入使得在价带顶附近出现了新的态密度峰。以Sr空位浓度为0.125的模型为例，在价带顶约-0.1--0.3eV的能量范围内，出现了明显的态密度峰。这些新的态密度峰主要由O2p态电子贡献。由于Sr空位的存在，使得周围O原子的电子云分布发生改变，O2p态电子的能量状态发生变化，从而在价带顶附近产生了新的态密度峰。随着Sr空位浓度的增加，这些新的态密度峰的强度逐渐增强，表明更多的空穴占据了这些缺陷能级。在导带底附近，态密度变化相对较小。Sr4d态电子的态密度略有下降，这是因为Sr原子的缺失导致Sr4d态电子的数量减少。在含Ti空位的SrTiO₃中，Ti空位的引入导致在价带顶附近出现了大量的新态密度峰。以Ti空位浓度为0.125的模型为例，在价带顶约-0.2--0.5eV的能量范围内，出现了五、SrTiO₃的替位掺杂研究5.1替位掺杂模型的构建5.1.1Y替位Sr掺杂模型构建Y替位Sr的SrTiO₃掺杂模型时，首先基于优化后的SrTiO₃立方晶系晶体结构。选用2×2×2的超胞，此超胞中包含8个Sr原子、8个Ti原子和16个O原子。为引入Y替位Sr的掺杂，在超胞中选取特定顶点位置的Sr原子，将其替换为Y原子。在确定掺杂浓度时，通过控制替换的Sr原子数量来实现。当掺杂浓度为0.125时，在超胞的8个Sr原子中，替换1个Sr原子为Y原子；当掺杂浓度为0.25时，替换2个Sr原子；当掺杂浓度为0.33时，替换3个Sr原子。在构建模型过程中，考虑到Y³⁺离子的价态为+3，而被替换的Sr²⁺离子价态为+2。为保持体系的电中性，在引入Y原子后，体系会额外引入一个电子。通过在计算中调整电子数，确保体系满足电中性条件。这样构建的Y替位Sr掺杂模型，能够准确模拟Y掺杂对SrTiO₃体系的影响。不同掺杂浓度的模型可以用于研究掺杂浓度与体系性能之间的关系，为后续分析Y替位Sr掺杂对SrTiO₃几何结构、电子结构以及物理性能的影响提供了基础。5.1.2In替位Ti掺杂模型在构建In替位Ti的SrTiO₃掺杂模型时，同样以优化后的SrTiO₃晶体结构为基础，采用2×2×2的超胞。在超胞中，Ti原子位于体心位置。为实现In替位Ti的掺杂，选择体心位置的Ti原子进行替换。在确定掺杂浓度方面，当掺杂浓度为0.125时，在超胞的8个Ti原子中，替换1个Ti原子为In原子；当掺杂浓度为0.25时，替换2个Ti原子；当掺杂浓度为0.33时，替换3个Ti原子。由于In³⁺离子的价态为+3，被替换的Ti⁴⁺离子价态为+4。为维持体系的电中性，在引入In原子后，体系会产生一个空穴。在计算过程中，通过调整电子数来体现空穴的产生，保证体系的电中性。在构建模型时，还需考虑In原子与周围原子的相互作用。In原子的半径和电子云分布与Ti原子不同，这会导致In-O键的键长和键角与Ti-O键有所差异。在优化结构时，充分考虑这种差异，使模型能够准确反映In替位Ti掺杂后的体系结构。通过精确构建In替位Ti掺杂模型，可以深入研究In掺杂对SrTiO₃体系的影响，为探究In掺杂对SrTiO₃电子结构和物理性能的改变提供可靠的模型基础。5.1.3Al替位Ti掺杂模型构建Al替位Ti的SrTiO₃掺杂模型，以优化后的SrTiO₃晶体结构为起始，采用2×2×2的超胞。在超胞中，确定Ti原子的位置后，进行Al原子对Ti原子的替换。对于掺杂浓度的控制，当掺杂浓度为0.125时，将超胞中8个Ti原子中的1个替换为Al原子；当掺杂浓度为0.25时，替换2个Ti原子；当掺杂浓度为0.33时，替换3个Ti原子。由于Al³⁺离子价态为+3，与被替换的Ti⁴⁺离子价态不同。为保证体系的电中性，在引入Al原子后，体系会产生空穴，通过调整电子数来实现电中性。Al原子的半径和电子结构与Ti原子存在差异，这会对体系的结构和电子环境产生影响。在构建模型时，考虑到Al原子与周围O原子之间的相互作用。Al-O键的键长和键角会因Al原子的特性而发生变化。在结构优化过程中，精确计算和调整这些参数，使模型能够准确反映Al替位Ti掺杂后的结构特征。通过构建不同掺杂浓度的Al替位Ti掺杂模型，可以系统地研究Al掺杂对SrTiO₃体系的影响，为分析Al掺杂对SrTiO₃电子结构和物理性能的作用提供有效的模型支持。5.2替位掺杂对几何结构的影响5.2.1晶格常数的变化通过第一性原理计算，得到不同替位掺杂模型的晶格常数，并与本征SrTiO₃的晶格常数进行对比，以分析掺杂原子的引入对晶格常数的影响。计算结果显示，Y替位Sr掺杂会使SrTiO₃的晶格常数发生变化。当Y掺杂浓度为0.125时，晶格常数从本征的3.902Å减小到3.898Å；当掺杂浓度增加到0.25时，晶格常数进一步减小至3.895Å；当掺杂浓度达到0.33时，晶格常数为3.892Å。这是因为Y³⁺离子的半径（0.90Å）略小于Sr²⁺离子的半径（1.18Å）。当Y替位Sr后，较小的Y³⁺离子进入晶格，使得周围原子间的距离减小，为了维持晶体结构的稳定性，晶格常数会相应减小。随着Y掺杂浓度的增加，更多的Y³⁺离子替换Sr²⁺离子，晶格常数的减小趋势更加明显。In替位Ti掺杂对晶格常数也有显著影响。当In掺杂浓度为0.125时，晶格常数增大到3.908Å；当掺杂浓度为0.25时，晶格常数进一步增大至3.915Å；当掺杂浓度为0.33时，晶格常数为3.922Å。In³⁺离子的半径（0.80Å）大于Ti⁴⁺离子的半径（0.605Å）。当In替位Ti后，较大的In³⁺离子进入晶格，导致周围原子间的距离增大，从而使晶格常数增大。随着In掺杂浓度的增加，更多的In³⁺离子替换Ti⁴⁺离子，晶格常数的增大趋势更加显著。Al替位Ti掺杂同样会改变晶格常数。当Al掺杂浓度为0.125时，晶格常数增大到3.906Å；当掺杂浓度为0.25时，晶格常数增大至3.913Å；当掺杂浓度为0.33时，晶格常数为3.920Å。Al³⁺离子的半径（0.535Å）略大于Ti⁴⁺离子的半径。当Al替位Ti后，会使周围原子间的距离增大，导致晶格常数增大。随着Al掺杂浓度的增加，更多的Al³⁺离子替换Ti⁴⁺离子，晶格常数的增大趋势更加明显。通过对不同替位掺杂模型晶格常数变化的分析，可以看出掺杂原子的半径和价态对晶格常数有重要影响。这种晶格常数的变化不仅反映了掺杂对晶体结构的影响，还会进一步影响材料的电子结构和物理性能。5.2.2键长与键角的变化替位掺杂不仅会导致晶格常数的改变，还会引起键长和键角的变化。在Y替位Sr的掺杂体系中，以掺杂浓度为0.125的模型为例，Y-O键长为2.680Å，相比本征Sr-O键长2.745Å有所减小。这是因为Y³⁺离子半径小于Sr²⁺离子半径，Y替位Sr后，Y与O之间的距离缩短。与Y原子相邻的O-Y-O键角为178°，相比本征O-Sr-O键角180°略有减小。随着Y掺杂浓度的增加，更多的Y-O键形成，键长和键角的变化趋势更加明显。这种键长和键角的变化会影响晶体内部的应力分布。键长的减小和键角的变化会导致晶体内部局部应力增加。当应力积累到一定程度时，可能会影响晶体的稳定性，甚至引发晶体结构的相变。在In替位Ti的掺杂体系中，以掺杂浓度为0.125的模型为例，In-O键长为2.010Å，大于本征Ti-O键长1.958Å。这是由于In³⁺离子半径大于Ti⁴⁺离子半径，In替位Ti后，In与O之间的距离增大。与In原子相邻的O-In-O键角为176°，相比本征O-Ti-O键角180°有所减小。随着In掺杂浓度的增加，In-O键长和O-In-O键角的变化更加显著。这种键长和键角的变化对晶体结构的局部和整体都有影响。在局部，会改变In-O八面体的结构；在整体上，会影响晶体的堆积方式和对称性，从而影响材料的力学性能和电学性能。在Al替位Ti的掺杂体系中，以掺杂浓度为0.125的模型为例，Al-O键长为1.980Å，大于本征Ti-O键长。这是因为Al³⁺离子半径略大于Ti⁴⁺离子半径，Al替位Ti后，Al与O之间的距离增大。与Al原子相邻的O-Al-O键角为177°，相比本征O-Ti-O键角180°略有减小。随着Al掺杂浓度的增加，Al-O键长和O-Al-O键角的变化更加明显。这种键长和键角的变化会对材料的电学性能产生潜在影响。键长和键角的改变会影响电子云的分布，从而影响电子的传输和材料的导电性能。通过对不同替位掺杂体系键长和键角变化的分析，可以深入了解掺杂对晶体结构的影响机制，以及这种影响对材料性能的潜在作用。5.3替位掺杂对电子结构的影响5.3.1能带结构的变化通过第一性原理计算，得到替位掺杂后的SrTiO₃的能带结构，并与本征SrTiO₃的能带结构进行对比，以分析掺杂原子对能带结构的影响。计算结果表明，Y替位Sr掺杂会使SrTiO₃的能带结构发生显著变化。在本征SrTiO₃中，价带顶位于布里渊区的R点，导带底位于Γ点，带隙大小为3.10eV。当Y替位Sr掺杂后，体系的费米能级进入导带中。以Y掺杂浓度为0.12