基于第一性原理探究掺杂立方BaTiO₃电子结构与光学性能的内在关联

基于第一性原理探究掺杂立方BaTiO₃电子结构与光学性能的内在关联一、引言1.1研究背景与意义在材料科学的广阔领域中，探索具有独特物理性质和潜在应用价值的材料始终是研究的核心主题。钛酸钡（BaTiO₃）作为一种典型的钙钛矿结构化合物，凭借其优异的铁电、压电、介电以及光学性能，在电子元器件、传感器、光电器件等众多领域展现出不可或缺的地位，被誉为“电子陶瓷工业的支柱”。纯净的BaTiO₃在室温下通常呈现四方相结构，然而，通过特定的制备工艺和条件，能够获得立方相的BaTiO₃。立方相BaTiO₃在某些应用场景中具有独特的优势，例如在一些对晶体对称性要求较高的光电器件中，立方相结构能够提供更加均匀的光学响应。同时，由于其结构特点，立方相BaTiO₃在理论研究中也成为了重要的模型体系，有助于深入理解钙钛矿结构材料的本征物理性质。随着现代科技的飞速发展，对材料性能的要求日益苛刻。为了进一步拓展BaTiO₃的应用范围并提升其性能，掺杂改性成为了一种行之有效的手段。通过向立方BaTiO₃晶格中引入特定的杂质原子（即掺杂），可以精确调控其电子结构。这种调控能够改变材料内部的电荷分布、能带结构以及电子的跃迁方式，进而对材料的光学性能产生显著影响。在光催化领域，合适的掺杂可以拓宽BaTiO₃对光的吸收范围，提高光生载流子的分离效率，从而增强其光催化活性，为解决环境污染和能源短缺等问题提供新的材料选择；在光电器件应用中，掺杂能够优化材料的发光特性、光电转换效率等关键性能指标，推动发光二极管、光电探测器等器件向高性能、小型化方向发展。例如，在光通信领域，具备特定光学性能的掺杂立方BaTiO₃可用于制造光调制器、光开关等核心器件，提高光信号的处理速度和传输效率，满足日益增长的高速通信需求。研究掺杂立方BaTiO₃的电子结构与光学性能，从基础研究层面来看，有助于深入揭示掺杂原子与BaTiO₃晶格之间的相互作用机制，丰富和完善材料电子结构与光学性能之间的内在联系理论。这种深入的理解不仅对于BaTiO₃基材料体系具有重要意义，还能够为其他钙钛矿结构材料以及功能材料的设计与开发提供理论指导和研究范式。从应用研究角度出发，精确掌握掺杂对立方BaTiO₃性能的影响规律，能够为新型光电器件的设计与制备提供坚实的材料基础。通过有针对性地选择掺杂元素和优化掺杂浓度，可以定制出具有特定功能和优异性能的材料，满足不同领域对高性能材料的迫切需求，推动相关产业的技术升级和创新发展。1.2研究目的本研究旨在通过深入探究掺杂对立方BaTiO₃电子结构与光学性能的影响，揭示其内在规律与作用机制，为该材料在光电器件、光催化等领域的应用提供坚实的理论支撑与数据基础。具体研究目的如下：揭示掺杂对电子结构的影响机制：借助先进的理论计算方法，如基于密度泛函理论（DFT）的第一性原理计算，精确分析不同掺杂元素（如过渡金属元素、稀土元素等）以及不同掺杂浓度下，立方BaTiO₃晶格内部电子云分布、能带结构（包括导带、价带位置及禁带宽度的变化）、电子态密度等电子结构特征的演变规律。例如，通过计算分析掺杂原子的电子轨道与BaTiO₃原有原子轨道的杂化情况，解释掺杂导致电子结构改变的微观机制，明确掺杂原子在晶格中所起的电子供体或受体作用，以及其对载流子浓度和迁移率的影响。阐明电子结构与光学性能的关联：建立立方BaTiO₃电子结构与光学性能之间的定量关系，深入理解光吸收、光发射、光折射等光学过程与电子跃迁、能带结构之间的内在联系。从理论层面分析不同掺杂体系中，由于电子结构变化所引起的光学吸收边的移动、光致发光峰的位移和强度变化等现象的物理本质。以光吸收为例，根据电子结构计算结果，解释掺杂如何改变材料对不同波长光的吸收能力，进而为调控材料的光吸收特性提供理论依据。优化材料性能并指导应用：基于对掺杂立方BaTiO₃电子结构与光学性能的研究成果，有针对性地提出材料性能优化策略。通过合理选择掺杂元素和精确控制掺杂浓度，实现对立方BaTiO₃光学性能的精准调控，满足不同应用场景对材料光学性能的特殊要求。在光催化领域，根据研究结果设计并制备具有高效光催化活性的掺杂立方BaTiO₃材料，提高其对太阳能的利用效率和对有机污染物的降解能力；在光电器件应用中，研发出具有特定发光波长、高发光效率和快速响应速度的掺杂立方BaTiO₃发光二极管或光电探测器等器件，推动相关领域的技术进步与创新发展。1.3国内外研究现状在材料科学领域，立方BaTiO₃由于其独特的晶体结构和物理性质，一直是国内外研究的重点对象。对立方BaTiO₃的研究最初主要聚焦于其晶体结构与基本物理性能的探索。早期研究借助X射线衍射（XRD）、透射电子显微镜（TEM）等实验手段，精确确定了立方BaTiO₃的晶体结构参数，为后续研究奠定了坚实基础。随着研究的深入，人们逐渐认识到立方BaTiO₃在光电器件、传感器等领域的潜在应用价值，从而引发了对其性能优化与改性的广泛研究。近年来，掺杂作为一种有效调控立方BaTiO₃性能的手段，受到了国内外科研人员的高度关注。在国外，诸多研究团队利用先进的理论计算方法和实验技术，深入探究了掺杂对立方BaTiO₃电子结构与光学性能的影响。美国某科研团队运用基于密度泛函理论（DFT）的第一性原理计算，系统研究了稀土元素掺杂对立方BaTiO₃电子结构的影响机制。结果表明，稀土元素的掺杂能够在立方BaTiO₃的禁带中引入新的能级，这些能级的位置和性质与掺杂元素的种类和浓度密切相关。这种能级的引入改变了材料内部的电子跃迁方式，进而对其光学性能产生显著影响。在实验方面，日本的研究人员通过溶胶-凝胶法制备了过渡金属掺杂的立方BaTiO₃薄膜，并利用光谱分析技术详细研究了其光学吸收和发光特性。研究发现，过渡金属掺杂能够有效拓宽立方BaTiO₃薄膜的光吸收范围，增强其在可见光区域的吸收强度，同时还能显著提高材料的光致发光效率，为其在光电器件中的应用提供了有力的实验依据。在国内，相关研究也取得了丰硕成果。众多科研机构和高校围绕掺杂立方BaTiO₃展开了深入研究，涵盖了理论计算、材料制备与性能表征等多个方面。一些研究团队采用第一性原理计算与实验相结合的方法，全面研究了不同掺杂体系下立方BaTiO₃的电子结构与光学性能之间的内在联系。通过精确控制掺杂元素的种类和浓度，成功实现了对立方BaTiO₃光学带隙的有效调控，使其能够满足不同光电器件的应用需求。在材料制备工艺方面，国内研究人员不断创新，开发出多种制备高质量掺杂立方BaTiO₃材料的方法，如熔盐法、水热法等。这些方法不仅能够精确控制材料的微观结构和掺杂均匀性，还能有效降低材料的制备成本，为掺杂立方BaTiO₃的大规模应用奠定了基础。例如，国内某高校利用熔盐法制备出了具有立方体形貌的BaTiO₃纳米粒子，该粒子结晶良好、纯度高，为后续的掺杂研究提供了优质的基体材料。尽管国内外在掺杂立方BaTiO₃的研究方面已经取得了显著进展，但仍存在一些亟待解决的问题。一方面，对于某些复杂掺杂体系下立方BaTiO₃的电子结构与光学性能之间的微观作用机制，尚未完全明晰，需要进一步深入研究。另一方面，如何在保证材料性能的前提下，实现掺杂立方BaTiO₃的大规模、低成本制备，也是目前面临的一个重要挑战。此外，随着科技的不断发展，对掺杂立方BaTiO₃在新型光电器件、光催化等领域的应用研究提出了更高的要求，需要进一步拓展其应用范围，挖掘其潜在应用价值。二、理论基础与研究方法2.1密度泛函理论密度泛函理论（DensityFunctionalTheory，DFT）是一种研究多电子体系电子结构的量子力学方法，在现代材料计算中占据着举足轻重的地位。其核心思想是将体系的能量表述为电子密度的泛函，通过对电子密度的求解来获取体系的各种性质，这一理论的诞生为多电子体系的研究开辟了新的道路。在传统的量子力学中，描述多电子体系需要求解复杂的多电子波函数，其变量数目随着电子数的增加而急剧增多，计算量呈指数级增长，这使得精确求解多电子体系的薛定谔方程变得极为困难。而DFT的出现巧妙地避开了这一难题，它以电子密度作为基本变量，取代了复杂的多电子波函数。电子密度是指在空间中不同位置上电子出现的概率密度分布，它仅仅是三个空间变量的函数，与多电子波函数相比，在概念和实际处理上都更加简便。DFT的理论基础建立在Hohenberg-Kohn定理之上。Hohenberg-Kohn第一定理指出，对于一个处在外部静电势中的多电子体系，其基态能量是电子密度的唯一泛函。这意味着，只要确定了体系的电子密度，就能够唯一确定体系的基态能量，从而为通过电子密度来研究体系性质提供了理论依据。Hohenberg-Kohn第二定理进一步证明，通过将体系能量对基态密度进行变分求极小值，就可以得到体系的基态能量。这两个定理的提出，为DFT奠定了坚实的理论基础，使得基于电子密度的计算成为可能。在实际应用中，DFT最普遍的实现方式是通过Kohn-Sham方法。Kohn-Sham方法将复杂的多体问题简化为一个没有相互作用的电子在有效势场中运动的问题。这个有效势场不仅包含了外部势场，还考虑了电子间库仑相互作用的影响，其中交换和相关作用是处理的关键难点。交换作用描述了电子由于全同性而产生的相互回避效应，相关作用则体现了电子之间的瞬时相互关联。目前，虽然没有精确求解交换相关能E_{XC}的方法，但科学家们提出了多种近似方法来逼近精确的泛函，其中最常用的是局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）。LDA近似假设体系中某点的交换相关能只与该点的电子密度有关，通过均匀电子气模型来计算体系的交换能，相关能部分则采用对自由电子气进行拟合的方法处理。GGA近似则在LDA的基础上进一步考虑了电子密度梯度的影响，能够更好地描述非均匀电子体系，在许多情况下提高了计算精度。在材料计算领域，DFT展现出了诸多显著优势。从计算效率上看，与传统的基于多电子波函数的方法相比，DFT在处理多电子体系时计算量大幅降低，能够在可接受的时间内对较大规模的体系进行计算。这使得研究人员能够对复杂的材料体系，如包含大量原子的晶体结构、具有复杂表面和界面的材料等进行深入研究。从适用范围来讲，DFT可广泛应用于原子、分子、固体、表面等不同尺度的体系，无论是研究简单的小分子的电子结构，还是探索复杂固体材料的物理性质，DFT都能提供有效的理论计算手段。在预测能力方面，DFT能够预测材料的多种性质，包括电子结构（如能带结构、电子态密度等）、光学性质（如光吸收、光发射等）、力学性质（如弹性常数、硬度等）以及化学反应活性等。这些预测结果为材料的设计、性能优化以及新功能材料的开发提供了重要的理论指导，帮助研究人员在实验之前对材料的性质和行为进行深入理解和预测，从而有针对性地开展实验研究，大大提高了材料研究的效率和成功率。2.2绝热近似与Hartree近似在多粒子体系的理论研究中，绝热近似（Born-Oppenheimer）和Hartree近似是两种重要的近似方法，它们在简化多粒子体系的复杂性、使得理论计算得以实现方面发挥了关键作用，为后续密度泛函理论等方法的应用奠定了基础。2.2.1绝热近似（Born-Oppenheimer）绝热近似，又称Born-Oppenheimer近似，其核心依据源于原子核与电子在质量和运动速度上的巨大差异。由于原子核的质量远远大于电子质量（通常原子核质量约为电子质量的数千倍），这使得电子的运动速度远快于原子核的运动速度。在实际的体系中，当电子快速运动时，原子核的位置几乎可以看作是瞬间固定的；反之，在考虑原子核的运动时，电子的具体空间分布对其影响也可以忽略不计。基于此，绝热近似将固体系统的总哈密顿量（无外场时）分解为原子核相关的部分和电子相关的部分。在求解电子态时，假定原子核处于瞬时位置固定不动，此时电子的薛定谔方程中仅涉及电子的坐标和相互作用；而在计算核的运动时，不考虑电子的具体空间分布，原子核的薛定谔方程仅与原子核的坐标和相互作用有关。通过这种近似处理，原本复杂的多粒子体系的运动问题被拆分成了相对独立的电子运动问题和原子核运动问题，大大降低了计算的难度。在材料计算领域，绝热近似有着广泛的应用。在研究晶体的电子结构时，借助绝热近似可以将注意力集中在电子在固定原子核势场中的运动状态，从而方便地求解电子的波函数和能量本征值，进而得到材料的能带结构、电子态密度等重要的电子结构信息。在研究材料的光学性质时，由于光与材料的相互作用主要涉及电子的跃迁，绝热近似下得到的电子结构信息能够为分析光吸收、光发射等光学过程提供坚实的基础。通过计算不同能级之间的电子跃迁概率和能量变化，能够准确预测材料对不同波长光的吸收和发射特性。在研究分子晶体的振动和转动光谱时，绝热近似使得我们可以分别考虑分子中电子的运动和原子核的振动、转动，从而清晰地解析出光谱的特征峰，为分子结构的确定和分子间相互作用的研究提供重要依据。2.2.2Hartree近似Hartree近似是一种将多电子问题简化为单电子问题的重要方法。在多电子体系中，电子之间存在着复杂的相互作用，使得直接求解多电子体系的薛定谔方程变得极为困难。Hartree近似的基本思路是将每个电子看作是在其他电子所产生的平均势场中独立运动。具体而言，假设体系中有N个电子，对于其中的第i个电子，Hartree近似认为它受到的其他电子的作用可以用一个平均的、自洽的势场V_{i}(r_{i})来描述，这个势场是由除第i个电子之外的其他N-1个电子的电荷分布所产生的。这样，多电子体系的薛定谔方程就可以被分解为N个单电子的薛定谔方程，每个单电子方程只涉及一个电子的坐标和这个平均势场。通过迭代求解这些单电子方程，不断调整平均势场，直到前后两次迭代得到的结果自洽，即达到收敛条件，从而得到每个单电子的波函数和能量。在材料科学的研究中，Hartree近似有着诸多应用。在研究金属材料的电子结构时，通过Hartree近似可以计算出金属中电子的分布和能量状态，进而解释金属的导电性、热导率等物理性质。在研究半导体材料时，Hartree近似能够帮助我们理解半导体的能带结构以及杂质对能带结构的影响，为半导体器件的设计和优化提供理论指导。例如，在研究掺杂半导体时，可以通过Hartree近似计算出掺杂原子引入的杂质能级以及这些能级对电子分布和输运性质的影响，从而深入了解掺杂半导体的电学性能变化规律。在研究分子材料时，Hartree近似可以用于计算分子的电子结构和化学键性质，解释分子的化学反应活性和稳定性，为药物分子设计、催化剂开发等领域提供重要的理论支持。2.2.3二者在材料计算中的作用绝热近似和Hartree近似在材料计算中相互配合，共同为复杂材料体系的研究提供了可行的途径。绝热近似将多粒子体系中的原子核和电子运动分离，使得我们能够先专注于电子在固定原子核势场中的行为，为后续的电子结构计算奠定基础。而Hartree近似进一步将多电子问题简化为单电子问题，使得我们可以通过求解相对简单的单电子方程来获取多电子体系的相关信息。这两种近似方法的结合，大大降低了材料计算的复杂性，使得我们能够对各种材料的电子结构、光学性能、力学性能等进行有效的理论研究和预测。在研究掺杂立方BaTiO₃时，绝热近似帮助我们将体系中的原子核和电子分开考虑，专注于电子结构的变化；Hartree近似则将多电子体系简化为单电子问题，方便我们计算掺杂原子对电子态的影响，从而深入理解掺杂对立方BaTiO₃电子结构与光学性能的作用机制。2.3计算软件与工具在本研究中，我们选用MaterialsStudio软件作为主要的计算工具。MaterialsStudio是一款功能强大且应用广泛的材料模拟软件，由美国Accelrys公司开发，专为材料科学领域的研究而设计，在材料模拟计算领域占据着重要地位。MaterialsStudio具备丰富的功能模块，这些模块相互配合，为材料研究提供了全方位的模拟手段。其中，MaterialsVisualizer是整个软件的图形化界面和核心基础环境，它集成了构建分子、晶体、界面、表面以及高分子材料结构所需的各类工具。研究人员可以通过该模块轻松地搭建各种复杂的材料结构模型，并对模型进行直观的操作、细致的观察以及深入的分析。它不仅能够处理文本、图形和表格等多种格式的数据，还支持与其他常用的字处理、电子表格和演示软件进行数据交换，极大地提高了数据处理和共享的便利性。例如，在构建立方BaTiO₃的晶体结构模型时，借助MaterialsVisualizer的直观操作界面，研究人员可以精确地定义原子的种类、位置和排列方式，快速搭建出符合研究需求的初始模型。在电子结构计算方面，MaterialsStudio中的Dmol3模块发挥着关键作用。Dmol3是一个基于密度泛函理论（DFT）的量子力学模拟模块，它能够提供精确的电子结构计算。该模块采用数值原子轨道基组，能够高效地处理各种分子和材料体系。在研究掺杂立方BaTiO₃的电子结构时，Dmol3模块可以准确计算体系的电子态密度、能带结构等重要电子结构信息。通过这些计算结果，我们能够深入了解掺杂原子对立方BaTiO₃电子结构的影响机制，如掺杂原子如何改变材料的电子云分布、如何影响能带的位置和宽度等。在计算过程中，Dmol3模块还可以灵活地选择不同的交换关联泛函，如局域密度近似（LDA）、广义梯度近似（GGA）及其各种衍生泛函，以满足不同精度要求的计算任务。对于光学性能的模拟，MaterialsStudio同样提供了相应的功能支持。结合电子结构计算结果，利用软件中的相关模块可以进一步计算材料的光学性质，如光吸收系数、光发射谱等。通过这些计算，我们能够从理论层面解释掺杂立方BaTiO₃在光学性能上的变化，如为什么某些掺杂体系会增强材料对特定波长光的吸收能力，或者改变材料的发光特性等。这种理论与实验相结合的研究方法，能够为实验研究提供有力的理论指导，帮助实验人员更好地理解实验现象，优化实验方案。MaterialsStudio还具有用户友好的界面设计，其操作符合微软标准，研究人员可以通过简单的对话框进行计算任务的设置和提交，即使是非计算机专业的科研人员也能快速上手。同时，该软件采用灵活的Client-Server结构，其核心模块Visualizer运行于客户端PC，而计算模块则运行于服务器端。这种结构使得用户可以将计算作业提交到网络上的任何一台服务器上进行计算，充分利用网络资源，提高计算效率。在处理大规模计算任务时，用户可以根据服务器的负载情况和自身需求，灵活选择合适的服务器进行计算，大大提高了研究工作的灵活性和效率。2.4光学性质理论基础材料的光学性质是其与光相互作用时所表现出的特性，这些性质不仅与材料的微观结构密切相关，还在众多领域有着广泛的应用，如光电器件、光通信、光催化等。在研究掺杂立方BaTiO₃的光学性质时，介电函数、吸收光谱、反射率等是重要的研究对象，它们各自有着独特的理论基础和计算方法。2.4.1介电函数介电函数（DielectricFunction）是描述材料在电场作用下电学响应特性的重要物理量，它在材料光学性质的研究中起着核心作用。从微观层面来看，当光照射到材料上时，光所携带的电场会与材料中的电子发生相互作用。这种相互作用会导致电子云发生畸变，电子会在电场的作用下产生位移，从而形成感应电偶极矩。介电函数正是用来量化这种材料对光场的响应程度，它反映了材料内部电子云在光场作用下的极化行为。介电函数通常表示为一个复数，即\varepsilon(\omega)=\varepsilon_{1}(\omega)+i\varepsilon_{2}(\omega)。其中，实部\varepsilon_{1}(\omega)被称为电容率，它主要描述材料在光场作用下储存电能的能力。当光照射到材料上时，一部分能量会被材料以电能的形式储存起来，\varepsilon_{1}(\omega)的值越大，说明材料储存电能的能力越强。在光学应用中，电容率与材料的折射率密切相关，通过折射率，我们可以进一步了解光在材料中的传播速度和方向的变化。虚部\varepsilon_{2}(\omega)则表示材料对光的吸收特性。当光与材料相互作用时，部分光的能量会被材料吸收，转化为材料内部电子的能量，从而导致光的强度减弱。\varepsilon_{2}(\omega)的值越大，说明材料对光的吸收能力越强。这种吸收过程与材料内部电子的跃迁过程紧密相关，当光的能量与材料中电子的能级差相匹配时，电子会吸收光子的能量从低能级跃迁到高能级，从而导致光的吸收。在理论计算中，介电函数可以通过电子结构计算结果结合一定的理论模型来获取。基于密度泛函理论（DFT）的计算方法能够得到材料的电子结构信息，如电子态密度、能带结构等。利用这些信息，我们可以通过Kramers-Kronig关系来计算介电函数。Kramers-Kronig关系是一种基于因果律和解析函数性质的数学关系，它将介电函数的实部和虚部联系起来。通过已知的介电函数虚部（通常可以通过计算电子跃迁概率得到），利用Kramers-Kronig关系就可以计算出介电函数的实部。这种计算方法在材料光学性质的理论研究中得到了广泛应用，为深入理解材料的光学响应机制提供了有力的工具。2.4.2吸收光谱吸收光谱（AbsorptionSpectrum）是材料光学性质的重要表征之一，它直观地展示了材料对不同波长光的吸收能力。当光照射到材料上时，光子的能量会与材料中的电子相互作用。如果光子的能量恰好等于材料中电子的能级差，电子就会吸收光子的能量，从低能级跃迁到高能级，这个过程就导致了光的吸收。吸收光谱就是通过测量材料对不同波长光的吸收程度，以波长或能量为横坐标，以吸收强度为纵坐标绘制而成的曲线。吸收光谱的形状和特征与材料的电子结构密切相关。在掺杂立方BaTiO₃中，由于掺杂原子的引入，会在材料的能带结构中引入新的能级。这些新能级的位置和性质会影响电子的跃迁过程，从而改变材料的吸收光谱。如果掺杂原子在禁带中引入了浅能级，那么在较低能量的光子照射下，电子就可以从价带跃迁到这些浅能级，从而导致在较低能量区域出现新的吸收峰。材料的晶体结构、晶格缺陷等因素也会对吸收光谱产生影响。晶体结构的对称性会影响电子跃迁的选择定则，从而影响吸收峰的强度和位置；晶格缺陷可能会导致局域电子态的出现，这些局域态也会参与光的吸收过程，使吸收光谱发生变化。在理论计算方面，吸收光谱可以通过计算材料的光吸收系数来得到。光吸收系数\alpha(\omega)与介电函数的虚部\varepsilon_{2}(\omega)密切相关，它们之间的关系可以表示为\alpha(\omega)=\frac{\omega}{c}\sqrt{\frac{\varepsilon_{2}(\omega)}{2}}，其中\omega是光的角频率，c是真空中的光速。通过前面介绍的介电函数计算方法，我们可以得到\varepsilon_{2}(\omega)，进而计算出光吸收系数\alpha(\omega)。将光吸收系数随波长或能量的变化关系绘制出来，就得到了材料的吸收光谱。这种理论计算方法能够从微观层面解释吸收光谱的形成机制，为实验研究提供理论指导。2.4.3反射率反射率（Reflectivity）是指光照射到材料表面时，被反射的光强度与入射光强度的比值，它反映了材料表面对光的反射能力。反射率的大小与材料的光学性质、表面状态以及光的入射角等因素密切相关。在光学应用中，反射率是一个关键参数，它直接影响着光在材料表面的传输和利用效率。在太阳能电池中，为了提高对太阳光的吸收效率，需要尽量降低电池表面的反射率，以减少光的损失；而在一些光学反射镜中，则需要提高反射率，以增强光的反射效果。从理论上讲，反射率可以通过菲涅尔公式（FresnelEquations）来计算。菲涅尔公式描述了光在两种不同介质界面上反射和折射时的电场强度和光强的变化关系。对于垂直入射的光，反射率R与材料的复折射率n=n_{r}+in_{i}（其中n_{r}是实折射率，n_{i}是虚折射率）之间的关系可以表示为R=\left|\frac{n-1}{n+1}\right|^{2}。而复折射率又与介电函数相关，n=\sqrt{\varepsilon}。通过这些关系，我们可以从材料的介电函数计算出反射率。在实际应用中，材料的表面状态对反射率有着显著影响。如果材料表面粗糙，光在表面会发生漫反射，这会使反射光的方向变得杂乱无章，从而降低了反射光的强度和方向性。表面的杂质、氧化层等也会改变材料的光学性质，进而影响反射率。在研究掺杂立方BaTiO₃的反射率时，不仅要考虑材料内部的电子结构对介电函数的影响，还要考虑材料表面的微观结构和化学组成对反射率的影响。通过对这些因素的综合分析，可以深入理解掺杂立方BaTiO₃的反射特性，为其在光学器件中的应用提供理论支持。三、立方相BaTiO₃的本征性质3.1理论模型构建在对立方相BaTiO₃的电子结构与光学性能进行深入研究之前，构建准确合理的理论模型是至关重要的基础工作。本研究采用基于密度泛函理论（DFT）的第一性原理计算方法，借助MaterialsStudio软件中的Dmol3模块来完成理论模型的构建。首先，明确立方相BaTiO₃的晶体结构特点。立方相BaTiO₃属于钙钛矿结构，其空间群为Pm-3m，具有高度的对称性。在这种结构中，Ba原子位于晶胞的八个顶点，Ti原子占据晶胞的体心位置，而O原子则处于晶胞六个面的面心。这种原子排列方式决定了立方相BaTiO₃的基本物理性质，也是构建理论模型的重要依据。利用MaterialsStudio软件中的MaterialsVisualizer模块，按照立方相BaTiO₃的晶体结构参数，精确搭建其初始晶胞模型。在搭建过程中，严格定义Ba、Ti、O三种原子的种类、坐标位置以及相互之间的键长、键角等几何参数。为了确保模型的准确性，参考大量的实验数据和理论研究成果，对所设置的结构参数进行反复验证和优化。根据文献报道，立方相BaTiO₃的晶格常数约为4.0Å，在构建模型时，将晶格常数设定在这一数值附近，并通过后续的结构优化过程使其达到能量最低的稳定状态。在构建模型时，还需要考虑电子相互作用的描述方法。如前文所述，密度泛函理论通过将体系能量表述为电子密度的泛函来处理多电子问题。在Dmol3模块中，采用广义梯度近似（GGA）来描述电子的交换关联能。GGA考虑了电子密度梯度的影响，相较于局域密度近似（LDA），能够更准确地描述电子的行为，从而提高计算结果的精度。同时，选择合适的基组对于准确描述原子的电子云分布和相互作用至关重要。本研究采用双数值加极化基组（DNP），该基组能够较好地描述原子的价电子和内层电子，在保证计算精度的前提下，有效控制计算量，提高计算效率。为了进一步验证模型的可靠性，对构建好的立方相BaTiO₃晶胞模型进行结构优化。在优化过程中，以体系总能量最低和原子受力最小为目标，通过不断调整原子的位置和晶胞的形状，使模型达到稳定状态。当原子间的相互作用力小于设定的收敛标准（通常为0.002Ha/Å），且体系总能量的变化小于一定阈值（如1×10⁻⁵Ha）时，认为结构优化收敛，此时得到的模型即为稳定的立方相BaTiO₃理论模型。通过与实验测得的晶格常数、原子间距离等结构参数进行对比，验证模型的准确性。如果计算结果与实验值存在偏差，分析可能的原因，如交换关联泛函的选择、基组的完备性等，并对模型进行相应的调整和优化，直至计算结果与实验数据相符或接近。3.2BaTiO₃的晶体结构与电子结构立方相BaTiO₃属于典型的钙钛矿结构，其晶体结构具有高度的对称性，空间群为Pm-3m。在这种结构中，Ba原子位于晶胞的八个顶点位置，每个Ba原子被八个晶胞所共有，因此在一个晶胞中Ba原子的有效数目为8×1/8=1。Ti原子占据晶胞的体心，其在晶胞中的位置独一无二，每个晶胞仅包含一个Ti原子。O原子处于晶胞六个面的面心，由于每个面心的O原子被两个晶胞所共有，所以一个晶胞中O原子的有效数目为6×1/2=3。这种原子排列方式使得立方相BaTiO₃的化学式为BaTiO₃，符合钙钛矿结构的通式ABO₃（其中A代表Ba，B代表Ti，O为氧原子）。从化学键的角度来看，在立方相BaTiO₃中，Ba²⁺离子半径较大，它与周围的O²⁻离子之间形成了相对较弱的离子键。这种离子键的形成主要源于Ba²⁺离子与O²⁻离子之间的静电吸引作用，其键长相对较长，键能相对较低。而Ti⁴⁺离子半径较小，它与O²⁻离子之间除了存在离子键成分外，还具有一定程度的共价键特征。这是因为Ti⁴⁺离子的外层电子结构与O²⁻离子的外层电子结构存在一定的轨道重叠，使得它们之间的化学键具有部分共价性。这种Ti-O键的共价性对立方相BaTiO₃的电子结构和物理性质产生了重要影响，它使得Ti-O键的键长相对较短，键能相对较高，增强了晶体结构的稳定性。利用MaterialsStudio软件中的Dmol3模块对立方相BaTiO₃的电子结构进行深入分析，得到了其能带结构和态密度分布。从能带结构来看，立方相BaTiO₃具有明显的价带和导带，二者之间存在一定宽度的禁带。价带主要由O2p轨道电子构成，这是因为O原子的电负性较大，其2p轨道电子具有相对较低的能量，在晶体中形成了价带。导带则主要由Ti3d轨道电子贡献，Ti原子的3d轨道电子在晶体中具有相对较高的能量，成为导带的主要组成部分。禁带宽度是材料电子结构的重要参数之一，它决定了材料的电学和光学性质。对于立方相BaTiO₃，计算得到的禁带宽度约为3.2eV（具体数值可能会因计算方法和参数设置的不同而略有差异），这表明立方相BaTiO₃在常温下是一种半导体材料。进一步分析态密度（DOS），可以更清晰地了解立方相BaTiO₃中电子的能量分布情况。在费米能级（EF）以下，存在着明显的价带态密度峰，主要集中在O2p轨道的能量范围内，这与前面关于价带主要由O2p轨道电子构成的结论一致。在费米能级以上，导带态密度逐渐上升，主要源于Ti3d轨道电子的贡献。通过对态密度的积分，可以得到不同能量区间内的电子占据数，从而深入了解电子在各个能级上的分布情况。在某些能量区间内，态密度可能会出现峰值，这意味着在这些能量处电子出现的概率较高，相应的能级对材料的物理性质可能具有重要影响。在靠近导带底的能量区域，态密度的变化情况会影响电子的激发和跃迁过程，进而影响材料的电学和光学性能。3.3BaTiO₃的光学性质立方相BaTiO₃的光学性质是其重要的物理特性之一，与材料的电子结构密切相关。通过基于密度泛函理论的第一性原理计算，结合MaterialsStudio软件中的相关模块，对立方相BaTiO₃的介电函数、吸收光谱、反射率等光学性质进行深入研究，能够揭示其光学行为的微观机制。3.3.1介电函数介电函数是描述材料在电场作用下电学响应特性的关键物理量，对于理解立方相BaTiO₃的光学性质具有重要意义。计算得到的立方相BaTiO₃的介电函数同样是一个复数，即\varepsilon(\omega)=\varepsilon_{1}(\omega)+i\varepsilon_{2}(\omega)。在低能量区域，介电函数的实部\varepsilon_{1}(\omega)相对较为稳定，随着光子能量的增加，\varepsilon_{1}(\omega)逐渐发生变化。这种变化与材料内部电子的极化响应密切相关。在低能量下，电子的极化主要是由原子核周围电子云的轻微畸变引起的，这种极化响应相对较弱，使得\varepsilon_{1}(\omega)变化较小。而当光子能量增加时，电子的极化程度增强，电子云的畸变更加明显，从而导致\varepsilon_{1}(\omega)发生显著变化。介电函数的虚部\varepsilon_{2}(\omega)则反映了材料对光的吸收特性。在某些特定的光子能量处，\varepsilon_{2}(\omega)会出现明显的峰值。这些峰值对应着材料内部电子的特定跃迁过程。通过与电子结构计算结果相结合分析发现，这些峰值主要是由于电子从价带顶附近的能级跃迁到导带底附近的能级所导致的。由于立方相BaTiO₃的价带主要由O2p轨道电子构成，导带主要由Ti3d轨道电子贡献，因此这些电子跃迁过程本质上是O2p电子向Ti3d电子的跃迁。当光子能量与这些能级之间的能量差相匹配时，电子就会吸收光子的能量发生跃迁，从而在\varepsilon_{2}(\omega)上表现出峰值。这种电子跃迁过程不仅决定了材料对光的吸收特性，还与材料的其他光学性质，如吸收光谱、反射率等密切相关。3.3.2吸收光谱立方相BaTiO₃的吸收光谱直观地展示了其对不同波长光的吸收能力，是研究其光学性质的重要方面。计算得到的吸收光谱表明，立方相BaTiO₃在紫外光区域具有较强的吸收能力。这主要是因为在紫外光区域，光子能量较高，能够激发电子从价带跃迁到导带，从而产生强烈的光吸收。随着波长的增加，进入可见光区域，吸收强度逐渐减弱。这是由于可见光的光子能量相对较低，不足以激发电子进行有效的跃迁，因此光吸收相对较弱。吸收光谱中的吸收峰位置与电子结构中的能级跃迁紧密相关。在前面分析介电函数时提到的电子从价带顶到导带底的跃迁，在吸收光谱中对应着主要的吸收峰。由于立方相BaTiO₃的禁带宽度约为3.2eV，根据光子能量与波长的关系E=hc/\lambda（其中E为光子能量，h为普朗克常数，c为光速，\lambda为波长），可以计算出对应于禁带宽度的吸收波长约为388nm，处于紫外光区域。这与吸收光谱中在紫外光区域出现的强吸收峰位置相符合。材料中的杂质、缺陷等因素也可能会在吸收光谱中引入额外的吸收峰。这些杂质或缺陷可能会在禁带中引入新的能级，使得电子能够在这些新能级与价带或导带之间进行跃迁，从而产生新的吸收峰。在实际应用中，通过控制材料中的杂质和缺陷含量，可以有效地调控立方相BaTiO₃的吸收光谱，以满足不同的光学应用需求。3.3.3反射率反射率是描述光在材料表面反射特性的重要参数，它反映了材料表面对光的反射能力。对于立方相BaTiO₃，其反射率与材料的介电函数以及光的入射角等因素密切相关。在垂直入射的情况下，根据菲涅尔公式，反射率R与材料的复折射率n=n_{r}+in_{i}（其中n_{r}是实折射率，n_{i}是虚折射率）之间的关系为R=\left|\frac{n-1}{n+1}\right|^{2}，而复折射率又与介电函数相关，n=\sqrt{\varepsilon}。通过计算得到的介电函数，可以进一步计算出立方相BaTiO₃的反射率。在不同波长的光照射下，立方相BaTiO₃的反射率呈现出一定的变化规律。在紫外光区域，由于材料对光的吸收较强，反射率相对较低。这是因为在紫外光区域，光子能量高，光与材料相互作用时，大部分能量被材料吸收，用于激发电子跃迁，从而导致反射光的强度较弱，反射率较低。在可见光区域，反射率相对较高。这是因为可见光的光子能量相对较低，材料对光的吸收较弱，光在材料表面反射的比例相对较大，因此反射率较高。材料的表面状态，如表面粗糙度、表面杂质等，对反射率也有显著影响。如果材料表面粗糙，光在表面会发生漫反射，使得反射光的方向变得杂乱无章，从而降低了反射光的强度和方向性，导致反射率降低。表面存在杂质或氧化层时，会改变材料表面的光学性质，进而影响反射率。在研究立方相BaTiO₃的反射率时，需要综合考虑材料的电子结构、介电函数以及表面状态等因素，以全面理解其反射特性。3.3.4光学性质与电子结构的关系立方相BaTiO₃的光学性质与电子结构之间存在着紧密的内在联系。从电子跃迁的角度来看，材料的光吸收和发射过程本质上是电子在不同能级之间的跃迁。在立方相BaTiO₃中，电子从价带跃迁到导带时，会吸收光子的能量，从而产生光吸收现象；反之，当电子从导带跃迁回价带时，会发射出光子，产生光发射现象。这种电子跃迁过程与材料的能带结构密切相关。能带结构中的价带和导带的位置、宽度以及能级分布情况，直接决定了电子跃迁的可能性和能量变化。禁带宽度的大小决定了光吸收的起始波长，只有当光子能量大于禁带宽度时，电子才能从价带跃迁到导带，从而产生光吸收。介电函数作为描述材料电学响应特性的物理量，与电子结构也有着深刻的关联。介电函数的实部和虚部分别反映了材料对光场的极化响应和吸收特性，而这些特性归根结底是由材料内部电子的行为所决定的。在光场的作用下，电子云会发生畸变，产生感应电偶极矩，这一过程与电子的轨道分布和能级结构密切相关。电子的跃迁过程会导致光的吸收，从而在介电函数的虚部表现出相应的峰值。通过对立方相BaTiO₃电子结构的深入分析，可以从微观层面解释其光学性质的本质，为进一步优化材料的光学性能提供理论依据。3.4本征性质小结立方相BaTiO₃具有独特的晶体结构，属于钙钛矿结构，空间群为Pm-3m。其晶体结构中，Ba原子位于晶胞顶点，Ti原子处于体心，O原子分布在面心，这种原子排列方式赋予了立方相BaTiO₃一定的物理性质。在化学键方面，Ba²⁺与O²⁻形成离子键，Ti⁴⁺与O²⁻之间除离子键外还具有一定共价键特征，这对材料的稳定性和电子结构产生重要影响。从电子结构角度来看，立方相BaTiO₃的价带主要由O2p轨道电子构成，导带主要由Ti3d轨道电子贡献，禁带宽度约为3.2eV，表现出半导体特性。通过对态密度的分析，进一步明确了电子在不同能级上的分布情况，为理解材料的电学和光学性质提供了微观基础。在光学性质方面，立方相BaTiO₃的介电函数、吸收光谱和反射率等特性与电子结构密切相关。介电函数的实部和虚部分别反映了材料对光场的极化响应和吸收特性，其虚部的峰值对应着电子从价带顶到导带底的跃迁过程。吸收光谱表明材料在紫外光区域有较强吸收，主要是由于光子能量足以激发电子跃迁，而在可见光区域吸收较弱。反射率则与材料的介电函数和光的入射角有关，在紫外光区域因吸收强而反射率低，在可见光区域则相对较高。立方相BaTiO₃的本征性质为后续研究掺杂对其电子结构与光学性能的影响提供了重要的基础和对比。通过深入理解其本征性质，能够更好地分析掺杂原子如何改变材料的电子结构和光学性质，从而为材料的性能优化和应用开发提供有力的理论支持。四、Ti位替代掺杂的影响4.1掺杂体系模型构建为深入探究Ti位替代掺杂对立方BaTiO₃电子结构与光学性能的影响，基于前文构建的立方相BaTiO₃本征结构模型，进一步构建Ti位替代掺杂体系模型。本研究选取过渡金属元素（如Cr、Mn、Fe、Co、Ni等）作为掺杂原子，这些元素因其独特的电子结构，在掺杂过程中能够与立方BaTiO₃晶格产生复杂的相互作用，从而对材料的电子结构和光学性能产生显著影响。在MaterialsStudio软件的MaterialsVisualizer模块中，对已优化的立方BaTiO₃晶胞模型进行操作。以Cr元素掺杂为例，首先在晶胞中选定位于体心位置的Ti原子，将其替换为Cr原子。为了模拟不同的掺杂浓度，构建不同比例的掺杂体系模型。当设定掺杂浓度为5%时，由于晶胞中Ti原子的总数为1，在一个晶胞中替换1个Ti原子为Cr原子，即可实现5%的掺杂浓度。若晶胞中Ti原子总数为20（可通过超胞构建获得更大的晶胞体系，以更精确地模拟掺杂情况），则替换1个Ti原子为Cr原子时，掺杂浓度为5%。通过这种方式，依次构建不同掺杂元素（Cr、Mn、Fe、Co、Ni等）在不同掺杂浓度（如2.5%、5%、7.5%等）下的立方BaTiO₃体系模型。在构建模型过程中，保持其他原子（Ba、O）的种类和位置不变，仅对Ti原子进行替换。同时，对构建好的掺杂体系模型进行结构优化，以确保模型的稳定性和合理性。在优化过程中，采用与本征结构模型优化相同的收敛标准，即原子间的相互作用力小于0.002Ha/Å，体系总能量的变化小于1×10⁻⁵Ha。通过结构优化，使掺杂原子在晶格中找到最稳定的位置，同时调整周围原子的位置，以最小化体系的能量。优化后的模型能够准确反映掺杂原子与立方BaTiO₃晶格之间的相互作用和几何结构关系，为后续的电子结构与光学性能计算提供可靠的基础。4.2晶格结构变化过渡金属元素在Ti位的掺杂会对立方BaTiO₃的晶格结构产生显著影响，这种影响主要体现在晶格常数和空间群的变化上，进而改变材料的物理性质。从晶格常数的变化来看，当Cr、Mn、Fe、Co、Ni等过渡金属元素替代Ti位时，由于这些掺杂原子与Ti原子的离子半径存在差异，会导致晶格发生畸变，从而引起晶格常数的改变。以Cr掺杂为例，Cr³⁺的离子半径（约0.615Å）略小于Ti⁴⁺的离子半径（约0.605Å）。当Cr原子取代Ti原子进入晶格后，由于其较小的离子半径，会使得周围原子向其靠近，从而导致晶格常数略微减小。根据第一性原理计算结果，在5%的Cr掺杂浓度下，立方BaTiO₃的晶格常数相较于本征结构减小了约0.01Å。随着Cr掺杂浓度的增加，晶格常数减小的趋势更为明显。在10%的Cr掺杂浓度下，晶格常数减小了约0.02Å。这是因为更多的Cr原子进入晶格，进一步加剧了晶格的畸变，使得晶格常数进一步减小。对于Mn掺杂，Mn²⁺的离子半径（约0.83Å）大于Ti⁴⁺的离子半径。当Mn原子取代Ti原子后，会使周围原子向外扩张，导致晶格常数增大。在5%的Mn掺杂浓度下，计算得到立方BaTiO₃的晶格常数相较于本征结构增大了约0.015Å。同样，随着Mn掺杂浓度的增加，晶格常数增大的幅度也会增加。在10%的Mn掺杂浓度下，晶格常数增大了约0.03Å。这种晶格常数随掺杂浓度的变化规律与掺杂原子和Ti原子的离子半径差异密切相关，离子半径差异越大，晶格常数的变化幅度就越大。从空间群的角度分析，在低掺杂浓度下，如2.5%的掺杂浓度时，立方BaTiO₃的空间群仍保持为Pm-3m，晶体结构的对称性未发生明显改变。这是因为少量的掺杂原子对整个晶格结构的影响较小，晶格能够维持其原有的高度对称性。然而，当掺杂浓度增加到一定程度时，情况会发生变化。当Cr、Mn等掺杂原子的浓度达到10%时，部分体系的空间群会发生改变。对于高浓度的Cr掺杂体系，由于Cr原子与周围原子之间的相互作用较强，可能会导致晶体结构的局部畸变加剧，从而使空间群从Pm-3m转变为其他对称性较低的空间群，如P4mm。这种空间群的转变意味着晶体结构的对称性降低，原子的排列方式发生了较大的变化，进而会对材料的物理性质产生深远影响。空间群的改变会影响材料的电学性质，使得材料的电导率、介电常数等发生变化；在光学性质方面，可能会导致材料的光吸收、光发射等特性发生改变。4.3形成能与结构稳定性形成能是衡量掺杂体系稳定性的关键物理量，它对于理解掺杂原子在立方BaTiO₃晶格中的存在状态以及掺杂体系的结构稳定性具有重要意义。通过基于密度泛函理论（DFT）的第一性原理计算，能够准确获取不同掺杂体系的形成能，进而深入分析掺杂原子占据Ti位的难易程度以及掺杂体系的稳定性。形成能的计算公式如下：E_{form}=E_{doped}-E_{pristine}+\sum_{i}n_{i}\mu_{i}其中，E_{form}表示掺杂体系的形成能，E_{doped}是掺杂后的体系总能量，E_{pristine}为未掺杂的立方BaTiO₃体系总能量，n_{i}代表第i种原子的数目变化，\mu_{i}则是第i种原子的化学势。化学势的取值与原子所处的化学环境密切相关，在实际计算中，通常参考相应元素的稳定相来确定化学势。对于过渡金属元素（如Cr、Mn、Fe、Co、Ni等），其化学势一般以金属单质的能量为参考；对于Ba、Ti、O等原子，其化学势则根据它们在立方BaTiO₃中的化学环境进行合理设定。以Cr掺杂立方BaTiO₃体系为例，在5%的掺杂浓度下，计算得到的形成能为-0.25eV。这表明Cr原子在该掺杂浓度下，进入立方BaTiO₃晶格占据Ti位的过程是一个能量降低的过程，即Cr掺杂体系具有一定的稳定性。形成能为负值，意味着掺杂过程在能量上是有利的，体系更倾向于形成掺杂结构。当Cr掺杂浓度增加到10%时，形成能变为-0.18eV。随着掺杂浓度的增加，形成能的绝对值减小，说明掺杂体系的稳定性有所降低。这是因为高掺杂浓度下，掺杂原子之间的相互作用增强，可能导致晶格畸变加剧，从而使体系的能量升高，稳定性下降。对于Mn掺杂体系，在5%的掺杂浓度时，形成能为-0.32eV，相较于相同掺杂浓度下的Cr掺杂体系，其形成能的绝对值更大。这说明在该掺杂浓度下，Mn原子进入立方BaTiO₃晶格占据Ti位比Cr原子更容易，Mn掺杂体系具有更高的稳定性。这一现象与Mn和Cr原子的电子结构以及它们与立方BaTiO₃晶格的相互作用有关。Mn原子的电子结构使得它与周围原子之间能够形成更稳定的化学键，从而降低了体系的能量，提高了掺杂体系的稳定性。当Mn掺杂浓度增加到10%时，形成能变为-0.25eV，同样随着掺杂浓度的增加，稳定性有所下降。通过对不同过渡金属元素掺杂立方BaTiO₃体系形成能的计算和分析，可以得出：在低掺杂浓度下，Cr、Mn、Fe、Co、Ni等过渡金属元素都能够以一定的稳定性占据Ti位，形成掺杂体系。然而，随着掺杂浓度的升高，掺杂原子之间的相互作用增强，晶格畸变加剧，导致形成能升高，掺杂体系的稳定性逐渐降低。不同掺杂元素的形成能存在差异，这反映了它们占据Ti位的难易程度不同。Mn元素在低掺杂浓度下具有较低的形成能，说明其占据Ti位相对容易，形成的掺杂体系稳定性较高；而Cr、Fe、Co、Ni等元素的形成能相对较高，占据Ti位的难度相对较大，掺杂体系的稳定性也相对较低。这些结果为进一步理解掺杂对立方BaTiO₃电子结构与光学性能的影响提供了重要的基础，因为掺杂体系的稳定性直接关系到掺杂原子在晶格中的存在状态和电子结构的变化，进而影响材料的光学性能。4.4电子结构改变过渡金属元素在Ti位的掺杂对立方BaTiO₃的电子结构产生了显著且复杂的影响，这种影响主要体现在杂质能级的出现以及禁带宽度的变化等方面，这些变化从微观层面深刻地改变了材料的电学和光学性质。当Cr、Mn、Fe、Co、Ni等过渡金属元素掺杂到立方BaTiO₃的Ti位时，在其能带结构中引入了新的杂质能级。以Cr掺杂体系为例，计算结果表明，在立方BaTiO₃的禁带中出现了由Cr原子的3d电子轨道贡献的杂质能级。这些杂质能级的位置与Cr原子的电子结构以及它与周围原子的相互作用密切相关。Cr原子的3d电子具有多个不同的能级，在掺杂体系中，由于与BaTiO₃晶格中的O原子形成化学键，其3d电子的能级发生了分裂和位移，从而在禁带中形成了特定位置的杂质能级。这些杂质能级的出现为电子的跃迁提供了新的通道。在光激发或电激发的作用下，电子可以从价带跃迁到这些杂质能级，或者从杂质能级跃迁到导带。这种新的跃迁路径改变了材料内部的电荷传输和能量转换过程，进而对材料的电学和光学性质产生重要影响。对于Mn掺杂体系，同样在禁带中引入了杂质能级。Mn原子的3d电子具有半充满的稳定结构，在掺杂后，其与周围原子的相互作用使得3d电子的能级发生变化，从而在禁带中形成杂质能级。与Cr掺杂体系不同的是，Mn掺杂体系中杂质能级的位置和性质有所差异。由于Mn原子的电子结构特点，其形成的杂质能级更靠近导带底。这使得电子从杂质能级跃迁到导带所需的能量相对较小，在一定程度上增加了电子跃迁的概率。这种差异导致Mn掺杂体系在电学和光学性质上表现出与Cr掺杂体系不同的特征。在电学性质方面，Mn掺杂体系可能具有更高的电导率，因为电子更容易被激发到导带参与导电；在光学性质方面，由于杂质能级靠近导带底，Mn掺杂体系可能对特定波长的光具有更强的吸收能力，从而改变材料的光学吸收光谱。掺杂对立方BaTiO₃禁带宽度的影响也是研究的重点之一。随着Cr、Mn、Fe、Co、Ni等过渡金属元素的掺杂，禁带宽度发生了明显的变化。以Cr掺杂为例，随着Cr掺杂浓度的增加，禁带宽度逐渐减小。在5%的Cr掺杂浓度下，禁带宽度从本征立方BaTiO₃的约3.2eV减小到约3.0eV；当Cr掺杂浓度增加到10%时，禁带宽度进一步减小至约2.8eV。这种禁带宽度的减小主要是由于杂质能级的引入以及Cr原子与周围原子之间的相互作用导致的。杂质能级的出现使得价带和导带之间的能量间隔减小，同时Cr原子与O原子之间的化学键变化也影响了电子的能量分布，进一步导致禁带宽度的收缩。对于Mn掺杂体系，禁带宽度同样随着掺杂浓度的增加而减小。在5%的Mn掺杂浓度下，禁带宽度减小到约3.05eV；当Mn掺杂浓度增加到10%时，禁带宽度减小至约2.9eV。与Cr掺杂体系相比，Mn掺杂体系在相同掺杂浓度下禁带宽度的减小幅度相对较小。这是因为Mn原子与Cr原子的电子结构和离子半径存在差异，导致它们与BaTiO₃晶格的相互作用方式不同。Mn原子的电子结构使得它与周围原子形成的化学键相对较弱，对禁带宽度的影响相对较小。禁带宽度的变化对立方BaTiO₃的电学和光学性质具有重要意义。在电学性质方面，禁带宽度的减小使得电子更容易从价带激发到导带，从而增加了材料的电导率。在光学性质方面，禁带宽度的减小导致材料对光的吸收边向长波方向移动，即红移现象。这意味着材料对更长波长的光具有吸收能力，从而改变了材料的光学吸收范围。在光催化应用中，禁带宽度的减小可以使材料吸收更多的可见光，提高光催化效率；在光电器件中，禁带宽度的变化会影响器件的发光波长和发光效率等性能。4.5光学性能变化Ti位过渡金属元素掺杂对立方BaTiO₃的光学性能产生了显著影响，这种影响在介电性能、吸收边位置、吸收系数和反射率等方面均有体现，这些变化与掺杂导致的电子结构改变密切相关。从介电性能方面来看，掺杂使得立方BaTiO₃的介电函数发生明显变化。以Cr掺杂体系为例，在低频区域，掺杂后的介电函数实部相较于本征立方BaTiO₃有所增大。这是因为Cr原子的掺杂引入了新的电子态，这些电子态在低频电场作用下能够更有效地参与极化过程，从而增强了材料的极化响应，导致介电函数实部增大。在高频区域，由于电子跃迁等过程的变化，介电函数的虚部也出现了明显的峰值位移和强度变化。Cr掺杂后，在某些特定的高频光子能量处，介电函数虚部的峰值强度增强，这意味着材料对这些频率的光吸收能力增强。这种增强与Cr掺杂引入的杂质能级有关，电子在这些杂质能级与价带或导带之间的跃迁概率增加，导致光吸收增强。掺杂对立方BaTiO₃吸收边位置的影响十分显著，随着Cr、Mn、Fe、Co、Ni等过渡金属元素的掺杂，吸收边发生了明显的移动。以Mn掺杂体系为例，计算结果表明，Mn掺杂使得立方BaTiO₃的吸收边向长波方向移动，即发生了红移现象。这是由于Mn掺杂导致禁带宽度减小，电子从价带跃迁到导带所需的能量降低。根据光子能量与波长的关系E=hc/\lambda（其中E为光子能量，h为普朗克常数，c为光速，\lambda为波长），能量降低意味着能够吸收波长更长的光，从而使吸收边红移。这种吸收边的红移现象在光催化应用中具有重要意义，它使得材料能够吸收更多的可见光，提高光催化效率。在降解有机污染物的光催化反应中，吸收边红移后的掺杂立方BaTiO₃能够利用更多的可见光能量来激发电子跃迁，产生更多的光生载流子，从而增强对有机污染物的降解能力。掺杂还导致立方BaTiO₃吸收系数和反射率发生变化。在吸收系数方面，不同掺杂元素和掺杂浓度下，吸收系数在不同波长范围内呈现出不同的变化趋势。对于Fe掺杂体系，在可见光区域，随着Fe掺杂浓度的增加，吸收系数逐渐增大。这是因为Fe掺杂引入的杂质能级在可见光区域提供了更多的电子跃迁通道，使得电子更容易吸收光子的能量发生跃迁，从而增强了对可见光的吸收。在反射率方面，掺杂后材料的反射率与本征立方BaTiO₃相比也有所改变。在某些波长范围内，掺杂体系的反射率降低，这意味着更多的光能够进入材料内部，参与光与材料的相互作用。在紫外光区域，Cr掺杂体系的反射率相较于本征结构有所降低，这使得更多的紫外光能够被材料吸收，用于激发电子跃迁等过程。材料的表面状态对反射率也有重要影响，掺杂过程可能会改变材料的表面微观结构，从而进一步影响反射率。如果掺杂导致材料表面粗糙度增加，光在表面会发生漫反射，使得反射光的方向变得杂乱无章，从而降低反射率。4.6Ti位掺杂小结Ti位过渡金属元素掺杂对立方BaTiO₃的结构、电子结构和光学性能产生了多方面的显著影响。在结构方面，掺杂原子与Ti原子的离子半径差异导致晶格常数发生变化，低掺杂浓度下空间群保持稳定，高浓度时部分体系空间群改变，结构对称性降低。形成能分析表明，低掺杂浓度时掺杂体系具有一定稳定性，高浓度下稳定性下降，不同掺杂元素占据Ti位的难易程度和稳定性存在差异。电子结构上，掺杂引入杂质能级，为电子跃迁提供新通道，改变了电子传输和能量转换过程。同时，禁带宽度随掺杂浓度增加而减小，影响材料的电学和光学性质，如使电导率增加，光吸收边红移。光学性能上，掺杂导致介电函数在低频和高频区域均发生变化，吸收边红移，在光催化等应用中有利于吸收更多可见光。吸收系数和反射率也随掺杂元素和浓度改变，在不同波长范围呈现不同变化趋势，材料表面状态也会影响反射率。这些研究结果为深入理解掺杂立方BaTiO₃的性质以及其在光电器件、光催化等领域的应用提供了重要的理论依据。五、Ba位替代掺杂的影响5.1掺杂体系模型构建为深入研究Ba位替代掺杂对立方BaTiO₃电子结构与光学性能的影响，本研究基于立方相BaTiO₃的本征结构模型，借助MaterialsStudio软件构建Ba位替代掺杂体系模型。在选择掺杂元素时，综合考虑了元素的电子结构、离子半径以及化学性质等因素，选取了La、Sr等元素作为Ba位的掺杂原子。La³⁺离子半径为1.36Å，与Ba²⁺离子半径1.61Å较为接近，在掺杂过程中能够较好地融入晶格，同时其独特的电子结构（5d¹6s²）可能会对立方BaTiO₃的电子结构产生显著影响。Sr²⁺离子半径为1.44Å，与Ba²⁺离子半径也具有一定的相似性，且Sr元素在一些材料体系中表现出良好的掺杂效果，能够有效调控材料的性能。在MaterialsStudio软件的MaterialsVisualizer模块中，对已优化的立方BaTiO₃晶胞模型进行操作。以La掺杂为例，在晶胞中选定位于顶点位置的Ba原子，将其替换为La原子。为模拟不同的掺杂浓度，构建不同比例的掺杂体系模型。当设定掺杂浓度为5%时，由于晶胞中Ba原子的总数为1（考虑晶胞的等效性，一个晶胞中Ba原子的有效数目为1），在一个晶胞中替换1个Ba原子为La原子，即可实现5%的掺杂浓度。若构建更大的超胞模型，例如包含20个Ba原子的超胞（通过合理的晶胞拓展获得），则替换1个Ba原子为La原子时，掺杂浓度为5%。通过这种方式，依次构建不同掺杂元素（La、Sr等）在不同掺杂浓度（如2.5%、5%、7.5%等）下的立方BaTiO₃体系模型。在构建模型过程中，严格保持其他原子（Ti、O）的种类和位置不变，仅对Ba原子进行替换。构建完成后，对掺杂体系模型进行结构优化，采用与本征结构模型优化相同的收敛标准，即原子间的相互作用力小于0.002Ha/Å，体系总能量的变化小于1×10⁻⁵Ha。通过结构优化，使掺杂原子在晶格中找到最稳定的位置，同时调整周围原子的位置，以最小化体系的能量。优化后的模型能够准确反映掺杂原子与立方BaTiO₃晶格之间的相互作用和几何结构关系，为后续的电子结构与光学性能计算提供可靠的基础。5.2晶格结构变化当La、Sr等元素在Ba位掺杂立方BaTiO₃时，会对其晶格结构产生显著影响，这种影响主要体现在晶格畸变和对称性变化等方面，进而改变材料的物理性质。从晶格畸变的角度来看，由于La³⁺离子半径（1.36Å）和Sr²⁺离子半径（1.44Å）与Ba²⁺离子半径（1.61Å）存在差异，当它们取代Ba位时，会导致晶格发生不同程度的畸变。以La掺杂为例，由于La³⁺离子半径小于Ba²⁺离子半径，当La原子进入晶格占据Ba位后，周围的O原子会向La原子靠近，从而使得晶格常数减小。根据第一性原理计算结果，在5%的La掺杂浓度下，立方BaTiO₃的晶格常数相较于本征结构减小了约0.02Å。随着La掺杂浓度的增加，晶格常数减小的趋势更为明显。在10%的La掺杂浓度下，晶格常数减小了约0.035Å。这种晶格常数的减小是由于更多的La原子进入晶格，进一步加剧了晶格的畸变，使得晶格整体收缩。对于Sr掺杂体系，由于Sr²⁺离子半径也小于Ba²⁺离子半径，同样会导致晶格常数减小。在5%的Sr掺杂浓度下，晶格常数减小了约0.015Å；当Sr掺杂浓度增加到10%时，晶格常数减小了约0.025Å。与La掺杂体系相比，Sr掺杂导致的晶格常数减小幅度相对较小。这是因为Sr²⁺离子半径与Ba²⁺离子半径的差异相对较小，所以对晶格结构的影响程度也相对较弱。从对称性变化方面分析，在低掺杂浓度下，如2.5%的掺杂浓度时，立方BaTiO₃的空间群仍保持为Pm-3m，晶体结构的对称性未发生明显改变。这是因为少量的掺杂原子对整个晶格结构的影响较小，晶格能够维持其原有的高度对称性。然而，当La、Sr等掺杂原子的浓度增加到一定程度时，情况会发生变化。当La掺杂浓度达到10%时，部分体系的空间群会发生改变。由于La原子与周围原子之间的相互作用较强，可能会导致晶体结构的局部畸变加剧，从而使空间群从Pm-3m转变为其他对称性较低的空间群，如P4mm。这种空间群的转变意味着晶体结构的对称性降低，原子的排列方式发生了较大的变化。空间群的改变会对材料的物理性质产生深远影响，它可能会改变材料的电学性质，使得材料的电导率、介电常数等发生变化；在光学性质方面，可能会导致材料的光吸收、光发射等特性发生改变。晶格结构的变化还会影响材料的力学性质。晶格畸变和对称性变化会改变原子间的键长和键角，从而影响材料的弹性常数和硬度。在高掺杂浓度下，由于晶格畸变加剧，原子间的相互作用力发生变化，材料的弹性常数可能会减小，硬度可能会降低。这种力学性质的变化在实际应用中需要加以考虑，特别是在一些对材料力学性能要求较高的场合，如电子器件的封装材料等。5.3形成能与结构稳定性形成能是评估Ba位掺杂体系稳定性的关键物理量，它对于理解掺杂原子在立方BaTiO₃晶格中的存在状态以及掺杂体系的稳定性具有重要意义。通过基于密度泛函理论（DFT）的第一性原理计算，能够准确获取不同掺杂体系的形成能，进而深入分析掺杂原子占据Ba位的难易程度以及掺杂体系的稳定性。形成能的计算公式为：E_{form}=E_{doped}-E_{pristine}+\sum_{i}n_{i}\mu_{i}其中，E_{form}表示掺杂体系的形成能，E_{doped}是掺杂后的体系总能量，E_{pristine}为未掺杂的立方BaTiO₃体系总能量，n_{i}代表第i种原子的数目变化，\mu_{i}则是第i种原子的化学势。化学势的取值与原子所处的化学环境密切相关，在实际计算中，通常参考相应元素的稳定相来确定化学势。对于La、Sr等掺杂元素，其化学势一般以金属单质的能量为参考；对于Ba、Ti、O等原子，其化学势则根据它们在立方BaTiO₃中的化学环境进行合理设定。以La掺杂立方BaTiO₃体系为例，在5%的掺杂浓度下，计算得到的形成能为-0.30eV。这表明La原子在该掺杂浓度下，进入立方BaTiO₃晶格占据Ba位的过程是一个能量降低的过程，即La掺杂体系具有一定的稳定性。形成能为负值，意味着掺杂过程在能量上是有利的，体系更倾向于形成掺杂结构。当La掺杂浓度增加到10%时，形成能变为-0.22eV。随着掺杂浓度的增加，形成能的绝对值减小，说明掺杂体系的稳定性有所降低。这是因为高掺杂浓度下，掺杂原子之间的相互作用增强，可能导致晶格畸变加剧，从而使体系的能量升高，稳定性下降。对于Sr掺杂体系，在5%的掺杂浓度时，形成能为-0.35eV，相较于相同掺杂浓度下的La掺杂体系，其形成能的绝对值更大。这说明在该掺杂浓度下，Sr原子进入立方BaTiO₃晶格占据Ba位比La原子更容易，Sr掺杂体系具有更高的稳定性。这一现象与Sr和La原子的电子结构以及它们与立方BaTiO₃晶格的相互作用有关。Sr原子的电子结构使得它与周围原子之间能够形成更稳定的化学键，从而降低了体系的能量，提高了掺杂体系的稳定性。当Sr掺杂浓度增加到10%时，形成能变为-0.28eV，同样随着掺杂浓度的增加，稳定性有所下降。通过对不同元素掺杂立方BaTiO₃体系形成能的计算和分析，可以得出：在低掺杂浓度下，La、Sr等元素都能够以一定的稳定性占据Ba位，形成掺杂体系。然而，随着掺杂浓度的升高，掺杂原子之间的相互作用增强，晶格畸变加剧，导致形成能升高，掺杂体系的稳定性逐渐降低。不同掺杂元素的形成能存在差异，这反映了它们占据Ba位的难易程度不同。Sr元素在低掺杂浓度下具有较低的形成能，说明其占据Ba位相对容易，形成的掺杂体系稳定性较高；而La元素的形成能相对较高，占据Ba位的难度相对较大，掺杂体系的稳定性也相对较低。这些结果为进一步理解掺杂对立方BaTiO₃电子结构与光学性能的影响提供了重要的基础，因为掺杂体系的稳定性直接关系到掺杂原子在晶格中的存在状态和电子结构的变化，进而影响材料的光学性能。5.4电子结构改变当La、Sr等元素在Ba位掺杂立方BaTiO₃时，会对其电子结构产生显著影响，这种影响主要体现在电子云分布、能带结构以及态密度等方面，进而改变材料的电学和光学性质。从电子云分布的角度来看，由于La³⁺和Sr²⁺的电子结构与Ba²⁺存在差异，掺杂后会导致晶格中电子云分布发生变化。以La掺杂为例，La³⁺具有5d¹6s²的电子结构，当它取代Ba²⁺进入晶格后，其5d电子会参与到与周围O原子的化学键形成中。La³⁺的5d电子与O2p电子