基于第一性原理探究热电化合物输运及相关性质

基于第一性原理探究热电化合物输运及相关性质一、引言1.1研究背景与意义随着全球工业化进程的迅猛推进，能源问题愈发凸显，传统化石能源储量逐渐减少，且在使用过程中带来了严重的环境污染问题，如温室气体排放引发的全球气候变暖、酸雨等，开发可再生清洁能源迫在眉睫。热电材料作为一种能够实现热能与电能直接相互转换的功能材料，在温差发电和热电制冷等领域展现出巨大的应用潜力，受到了广泛关注。热电材料的应用可以有效地将废热转化为电能，提高能源利用效率，减少对环境的负面影响，符合可持续发展的理念。热电材料的工作原理基于塞贝克效应、珀尔帖效应和汤姆逊效应。其中，塞贝克效应是指当两种不同材料的导体构成闭合回路时，若在回路中存在温度差，就会产生电动势，这是温差发电的基础；珀尔帖效应则是在两种不同材料的导体之间施加电压，当电流通过时，两导体之间会产生温差，实现了电能向热能的转换，常用于热电制冷；汤姆逊效应是在温度梯度作用下，电流通过两种不同材料的导体时，会在两导体之间产生热量的吸收或释放。热电材料的性能通常用热电优值ZT来衡量，ZT=S²σT/κ，其中S为塞贝克系数，代表材料的热电转换效率；σ为电导率；T为绝对温度；κ为热导率。为了获得较高的热电优值ZT，材料需要具备高的塞贝克系数、高的电导率以及低的热导率。然而，这些性能参数之间往往存在相互制约的关系，这给高性能热电材料的研发带来了巨大挑战。在探索高性能热电材料的征程中，理论计算发挥着不可或缺的作用。第一性原理计算方法，作为一种基于量子力学基本原理的计算方法，能够从原子和电子层面深入探究材料的性质。它无需依赖实验数据或经验参数，仅凭借基本的物理常数，如电子质量、光速、质子中子质量等，就可以精确计算出材料体系的能量、电子结构等关键物理性质，为热电化合物的研究提供了坚实的理论支撑。通过第一性原理计算，科研人员可以深入了解热电化合物的电子结构、能带结构、声子谱等微观信息，进而揭示热电性能与材料微观结构之间的内在联系，为高性能热电材料的设计与开发提供重要的理论指导。这种从微观层面揭示材料性能本质的研究方法，有助于突破传统实验研究的局限性，加速新型热电材料的研发进程，为解决能源问题开辟新的道路。1.2研究目的与内容本研究旨在运用第一性原理计算方法，深入剖析热电化合物的输运及相关性质，揭示其内在物理机制，为高性能热电材料的设计与开发提供坚实的理论依据。具体研究内容如下：热电化合物的结构优化与稳定性分析：运用第一性原理计算方法，对多种热电化合物的晶体结构进行全面优化，精准确定其最稳定的结构形态。通过计算形成能、结合能等关键参数，深入分析热电化合物的稳定性，为后续研究奠定坚实基础。形成能是衡量化合物形成难易程度的重要指标，结合能则反映了原子间相互作用的强弱，对这些参数的深入研究有助于理解热电化合物的结构稳定性。电子结构与能带结构研究：深入探究热电化合物的电子结构和能带结构，详细分析其电子态密度、能带色散关系等重要性质。通过对这些性质的研究，揭示电子在材料中的运动规律，明确影响热电性能的关键电子因素。例如，电子态密度反映了电子在不同能量状态下的分布情况，能带色散关系则描述了电子能量与波矢之间的关系，这些信息对于理解热电材料的电学性能至关重要。输运性质的计算与分析：借助第一性原理计算和声子玻尔兹曼输运理论，精确计算热电化合物的电导率、热导率、塞贝克系数等输运性质。深入分析这些性质与材料结构、电子结构之间的内在联系，揭示输运过程中的物理机制。例如，电导率与电子的迁移率和载流子浓度密切相关，热导率则受到声子散射等因素的影响，通过对这些关系的研究，可以深入理解热电材料的输运性质。掺杂与缺陷对热电性能的影响：系统研究掺杂和缺陷对热电化合物性能的影响规律。通过在材料中引入不同类型的杂质原子或制造缺陷，观察其对电子结构、输运性质以及热电性能的调控作用。深入分析掺杂和缺陷对热电性能影响的微观机制，为通过材料改性提高热电性能提供理论指导。例如，掺杂可以改变材料的载流子浓度和能带结构，缺陷则可以影响声子散射，从而对热电性能产生显著影响。1.3国内外研究现状近年来，国内外学者在热电化合物的第一性原理研究方面取得了丰硕的成果。在电子结构与热电性能关系的研究上，众多学者通过第一性原理计算深入剖析了热电化合物的电子结构，发现电子态密度、能带结构等因素对热电性能有着关键影响。例如，对于某些化合物，其能带结构中的特定能级分布会影响电子的跃迁概率，进而影响电导率和塞贝克系数。在研究中还发现，一些材料具有特殊的电子结构，如多带结构或能带简并，这些特性能够增强电子的散射或促进电子的传输，从而对热电性能产生积极影响。在掺杂和缺陷调控热电性能的研究领域，研究者们通过第一性原理计算系统地研究了掺杂和缺陷对热电化合物性能的影响。研究表明，引入合适的杂质原子或制造特定类型的缺陷，可以有效地改变材料的电子结构和载流子浓度，进而优化热电性能。例如，在某些材料中，通过掺杂可以引入额外的载流子，提高电导率；而在另一些材料中，缺陷的存在可以增强声子散射，降低热导率，从而提高热电优值。在新型热电化合物的探索方面，研究者们利用第一性原理计算预测了多种具有潜在优异热电性能的新型化合物，并对其结构和性能进行了深入研究。这些新型化合物的发现为热电材料的发展提供了新的方向和选择。例如，通过计算模拟，发现了一些具有特殊晶体结构或原子排列的化合物，它们在理论上具有较高的热电优值，为实验合成和进一步研究提供了重要参考。然而，现有研究仍存在一些不足之处。一方面，第一性原理计算的精度和效率仍有待进一步提高。尽管当前的计算方法能够提供较为准确的结果，但在处理大规模体系或复杂物理过程时，计算资源的消耗仍然较大，计算时间较长，这限制了研究的范围和深度。另一方面，理论计算与实验研究的结合还不够紧密。虽然理论计算能够预测材料的性能，但实验验证的过程往往较为复杂，需要耗费大量的时间和资源。此外，实验条件的限制也可能导致理论与实验结果之间存在一定的偏差，如何更好地将理论计算与实验研究相结合，实现二者的相互验证和补充，仍然是一个亟待解决的问题。1.4研究方法与创新点本研究主要采用基于密度泛函理论（DFT）的第一性原理计算方法，借助VASP等计算软件开展研究工作。在计算过程中，运用投影缀加波（PAW）方法描述离子实与价电子之间的相互作用，采用广义梯度近似（GGA）处理交换关联能。为了获得准确的电子结构和能量信息，选取合适的平面波截断能和k点网格进行细致计算。同时，利用玻尔兹曼输运理论结合第一性原理计算结果，深入计算热电化合物的输运性质。在研究过程中，将理论计算与实验数据进行紧密对比和验证，确保研究结果的准确性和可靠性。本研究的创新点主要体现在以下几个方面：一是采用多尺度计算方法，从原子尺度的第一性原理计算出发，结合宏观尺度的输运理论，全面研究热电化合物的性能，这种多尺度的研究方法能够更深入地揭示热电性能的内在机制；二是在研究掺杂和缺陷对热电性能的影响时，不仅考虑了常见的杂质原子和简单缺陷，还创新性地引入了一些新型的掺杂策略和缺陷模型，为热电材料的性能调控提供了新的思路；三是将机器学习算法与第一性原理计算相结合，建立热电性能预测模型，通过对大量计算数据的学习和分析，快速筛选和预测具有潜在优异性能的热电化合物，提高研究效率，为新型热电材料的探索开辟了新的途径。二、热电化合物及第一性原理计算基础2.1热电化合物概述2.1.1热电化合物的定义与分类热电化合物，又称温差电材料，是一类基于塞贝克（Seebeck）效应及珀尔帖（Peltier）效应，能够实现热能与电能直接相互转换的功能材料。其独特的性质使得在能源领域展现出巨大的应用潜力，为解决能源问题提供了新的途径。热电化合物种类繁多，根据其组成元素和晶体结构的差异，常见的分类方式包括IV-VIA族、III-V族等化合物。IV-VIA族热电化合物中，碲化铅（PbTe）是典型代表，具有面心立方晶体结构，在中高温区（约500-800K）展现出优异的热电性能。其晶体结构中，铅原子和碲原子通过离子键和共价键相互作用，形成稳定的晶格结构，这种结构特性对其热电性能产生了重要影响。在一定温度范围内，PbTe的载流子迁移率较高，有利于提高电导率，同时通过合理的元素掺杂，如引入铊（Tl）、银（Ag）等杂质原子，可以有效地调控其载流子浓度，进一步优化热电性能。III-V族化合物中，砷化镓（GaAs）也是研究较多的热电材料。GaAs具有闪锌矿结构，由镓原子和砷原子交替排列组成。它在电子学领域应用广泛，同时在特定条件下也展现出一定的热电性能。GaAs的电子迁移率较高，这使得它在热电转换过程中能够有效地传输电子，但其热导率相对较高，限制了其热电优值的进一步提高。为了降低热导率，研究人员通过纳米结构化、合金化等方法，引入界面散射和声子散射中心，有效地降低了热导率，提高了热电性能。除了上述常见的分类，还有一些其他类型的热电化合物，如II-VI族化合物（如硫化镉CdS、硒化锌ZnSe等）、V-VI族化合物（如锑化铋Bi2Te3等）以及一些复杂的氧化物、硫族化合物等。这些化合物各自具有独特的晶体结构和电子特性，在不同的温度范围和应用场景中展现出不同的热电性能。2.1.2热电效应及性能参数热电效应主要包括Seebeck效应、Peltier效应和Thomson效应，这些效应是热电化合物实现热能与电能相互转换的基础。Seebeck效应，又称第一热电效应，是指当两种不同材料的导体构成闭合回路时，若在回路中存在温度差，就会产生电动势。这种效应的产生源于材料中电子的热运动。当存在温度差时，高温端的电子具有较高的能量，会向低温端扩散，从而在材料两端形成电荷积累，产生电势差。Seebeck效应是温差发电的基础，其产生的电动势大小与材料的Seebeck系数以及温度差成正比，Seebeck系数（S）的表达式为S=\frac{\DeltaV}{\DeltaT}，其中\DeltaV为产生的电势差，\DeltaT为温度差。Peltier效应，也称为第二热电效应，是在两种不同材料的导体之间施加电压，当电流通过时，两导体之间会产生温差。这是因为电流通过时，电子在两种材料的界面处会发生能量的转移，导致热量的吸收或释放。Peltier效应常用于热电制冷，其制冷量与电流大小和Peltier系数成正比，Peltier系数（\Pi）与Seebeck系数之间存在关系\Pi=ST，其中T为绝对温度。Thomson效应，即第三热电效应，是在温度梯度作用下，电流通过两种不同材料的导体时，会在两导体之间产生热量的吸收或释放。这种效应与材料的电子结构和温度梯度密切相关，热释放或吸收的大小随材料的导热性能和温度梯度的大小而变化。热电材料的性能通常用多个参数来衡量，其中塞贝克系数（S）、电导率（\sigma）、热导率（\kappa）和热电优值（ZT）是最为重要的参数。塞贝克系数反映了材料将热能转化为电能的能力，其大小与材料的电子结构、载流子浓度等因素有关。电导率表示材料传导电流的能力，它与载流子的迁移率和浓度密切相关，载流子迁移率越高、浓度越大，电导率越高。热导率衡量材料传导热量的能力，热导率越低，材料在保持温度差方面的能力越强，有利于提高热电性能。热电优值ZT是综合评价热电材料性能的关键参数，其计算公式为ZT=\frac{S^{2}\sigmaT}{\kappa}，其中T为绝对温度。ZT值越高，表明材料的热电转换效率越高，当ZT>1时，热电材料具有较好的应用前景；而当ZT>3时，热电材料有望在实际应用中与传统能源转换技术相竞争。在实际应用中，提高热电优值ZT是研发高性能热电材料的核心目标，需要通过优化材料的结构和成分，协同调控塞贝克系数、电导率和热导率等参数来实现。2.2第一性原理计算方法2.2.1理论基础第一性原理计算方法是基于量子力学基本原理的计算方法，其核心是通过求解多电子体系的薛定谔方程来获取体系的能量、电子结构等物理性质。薛定谔方程是量子力学的基本方程，它描述了微观粒子的运动状态，其一般形式为：H\psi=E\psi，其中H为哈密顿算符，它包含了体系中电子的动能、电子与原子核之间的相互作用能以及电子之间的相互作用能；\psi为波函数，它描述了电子在空间中的概率分布；E为体系的能量。对于多电子体系，由于电子之间存在复杂的相互作用，直接求解薛定谔方程非常困难，通常需要采用一些近似方法。在第一性原理计算中，常用的近似方法是密度泛函理论（DFT）。DFT的核心思想是将多电子体系的基态能量表示为电子密度的泛函，而不是波函数。通过引入交换关联泛函来描述电子之间的交换关联作用，将多电子问题转化为单电子问题，从而大大降低了计算复杂度。在DFT框架下，常用的交换关联泛函有局域密度近似（LDA）和广义梯度近似（GGA）。LDA假设电子气的交换关联能只与电子密度的局域值有关，计算简单，但对于一些体系的描述存在一定的误差；GGA则考虑了电子密度的梯度信息，对体系的描述更加准确，在许多情况下能够给出与实验值更为接近的结果。2.2.2计算软件与工具在第一性原理计算领域，有许多功能强大的计算软件，它们为科研人员深入研究材料的性质提供了有力的工具。其中，VASP（ViennaAbinitioSimulationPackage）和CASTEP（CambridgeSequentialTotalEnergyPackage）是两款应用广泛的软件，各自具有独特的特点和适用范围。VASP是由维也纳大学Hafner小组开发的一款功能强大的第一性原理计算软件，在材料科学、凝聚态物理、化学和纳米技术等众多领域得到了广泛应用。它采用平面波基组和赝势方法来处理电子-离子相互作用，能够精确求解多体薛定谔方程。在电子结构计算方面，VASP表现出色，能够精确预测材料的电子能带结构，通过分析能带结构，可以深入了解材料的导电性质，判断材料是导体、半导体还是绝缘体。同时，它还能进行态密度（DOS）分析，研究材料的电子态分布特征，为理解材料的电子性质提供重要信息。VASP支持自旋极化计算，这使得它在研究磁性材料性质时具有很大的优势，能够准确描述磁性材料中电子的自旋状态和磁相互作用。对于研究分子动力学，VASP具备强大的第一性原理分子动力学（AIMD）模拟能力，基于量子力学描述原子运动，能够研究材料在非零温度下的动力学行为，分析化学反应机制，模拟材料在不同条件下的相变过程以及研究点缺陷、位错等缺陷的动力学行为。然而，VASP也存在一些局限性，它在计算大体系时，对计算资源的需求较高，计算时间较长。CASTEP是MaterialsStudio软件中的一个重要模块，专门用于第一性原理计算。它基于密度泛函理论，采用平面波赝势方法，能够对晶体、表面、界面等各种材料体系进行精确的结构优化和性质计算。CASTEP的优势在于其具有友好的图形用户界面，操作相对简单，对于初学者来说更容易上手。在结构优化方面，CASTEP能够快速准确地确定材料的稳定结构，通过优化原子坐标和晶格参数，使体系的能量达到最低。它还可以计算材料的弹性常数、声子谱等性质，对于研究材料的力学性能和热学性能具有重要意义。此外，CASTEP支持多种交换关联泛函的选择，用户可以根据具体的研究体系和需求，灵活选择合适的泛函进行计算。不过，与VASP相比，CASTEP在某些复杂体系的计算精度上可能略逊一筹。除了VASP和CASTEP，还有一些其他的第一性原理计算软件，如QuantumESPRESSO、SIESTA等。QuantumESPRESSO是一款开源的计算软件，具有高效的计算性能和丰富的功能，能够进行电子结构计算、分子动力学模拟等多种计算任务。SIESTA则采用了线性组合原子轨道基组，在处理大分子和表面体系时具有独特的优势。这些软件各有特点，科研人员可以根据研究的具体需求和体系特点，选择合适的计算软件进行研究。2.2.3计算流程与关键参数设置第一性原理计算的流程通常包括模型构建、参数设置、计算任务提交和结果分析等几个关键步骤，每个步骤都对计算结果的准确性和可靠性有着重要影响。在模型构建阶段，需要根据研究对象的特点，选择合适的晶体结构或分子模型，并确定原子的坐标和晶格参数。对于晶体材料，需要明确其晶体结构类型，如面心立方、体心立方、六方密堆积等，并根据实验数据或理论知识确定晶格常数。对于分子体系，要准确构建分子的三维结构，考虑分子的几何构型和原子间的键长、键角等信息。在构建模型时，还需要考虑体系的周期性边界条件，通过设置超胞来模拟材料的宏观性质。例如，在研究晶体材料的表面性质时，通常会构建包含表面层和一定厚度体内层的超胞模型，以准确描述表面原子的电子结构和相互作用。参数设置是第一性原理计算中至关重要的环节，直接影响计算结果的精度和计算效率。其中，平面波截断能和k点网格是两个关键参数。平面波截断能决定了平面波基组的大小，它表示平面波展开中所包含的最大动能。截断能设置过低，会导致基组不完备，计算结果不准确；而截断能设置过高，则会增加计算量，浪费计算资源。因此，需要通过测试不同的截断能，观察体系能量、原子力等物理量的收敛情况，来确定合适的截断能。一般来说，对于大多数材料体系，截断能在300-500eV范围内能够得到较好的收敛结果。k点网格用于采样倒易空间，它决定了对布里渊区积分的精度。k点网格越密，对布里渊区的采样越精细，计算结果越准确，但计算量也会相应增加。通常采用Monkhorst-Pack方法来生成k点网格，对于晶体材料，需要根据晶体结构的对称性和体系大小，合理选择k点网格的密度。例如，对于简单立方结构的晶体，可采用较稀疏的k点网格；而对于复杂结构的晶体或研究体系对电子结构精度要求较高时，则需要采用更密的k点网格。在计算过程中，还需要选择合适的交换关联泛函，如前文所述的LDA和GGA等，不同的泛函对体系的描述能力不同，需要根据具体研究对象进行选择。完成模型构建和参数设置后，即可提交计算任务。计算任务通常在高性能计算集群上运行，利用多处理器并行计算来加速计算过程。在计算过程中，需要密切关注计算的运行状态，检查是否存在错误或异常情况。如果计算出现问题，需要及时分析原因，调整参数或模型，重新提交计算。计算任务完成后，需要对计算结果进行详细的分析。首先，检查体系的能量是否收敛，原子力是否满足收敛标准，以确保计算结果的可靠性。然后，分析电子结构相关的结果，如能带结构、态密度等，了解材料的电子特性和电学性质。通过能带结构可以判断材料的能带间隙大小，确定材料是半导体还是金属，以及分析电子在能带中的分布情况和跃迁特性。态密度分析则可以提供电子在不同能量状态下的分布信息，帮助理解材料的电子态密度与能量的关系。此外，还可以分析电荷密度分布，直观地了解电子在实空间中的分布情况，研究原子间的成键特性和电荷转移情况。对于涉及分子动力学模拟的计算，还需要分析原子的运动轨迹、速度分布等信息，研究材料的动力学行为和热力学性质。在第一性原理计算中，模型构建的准确性、参数设置的合理性以及结果分析的全面性和深入性，是获得可靠计算结果的关键。通过合理选择计算软件，严格按照计算流程进行操作，并根据研究对象的特点和需求优化关键参数设置，能够为热电化合物等材料的研究提供准确、深入的理论依据。三、热电化合物的电子结构与输运性质3.1能带结构与态密度分析3.1.1能带结构计算与分析以碲化铅（PbTe）这一典型的热电化合物为研究对象，运用基于密度泛函理论（DFT）的第一性原理计算方法，借助VASP软件对其能带结构展开精确计算。在计算过程中，选用广义梯度近似（GGA）来处理交换关联能，采用投影缀加波（PAW）方法描述离子实与价电子之间的相互作用，平面波截断能设定为500eV，k点网格采用Monkhorst-Pack方法生成，设置为8×8×8，以确保计算结果的准确性和可靠性。图1展示了PbTe的能带结构计算结果。从图中可以清晰地观察到，PbTe具有典型的半导体能带结构特征，存在明显的能带间隙。其价带顶位于布里渊区的L点，导带底位于Γ点，这种能带结构特性使得PbTe在一定条件下能够实现电子的激发和跃迁，从而参与电输运过程。在低温状态下，由于热激发能量较低，电子主要分布在价带中，难以跨越能带间隙进入导带，因此电导率较低。随着温度逐渐升高，热激发能量增大，更多的电子获得足够的能量从价带跃迁到导带，导带中的电子浓度增加，电导率随之提高。同时，温度的升高也会导致电子的散射增强，这在一定程度上会限制电导率的进一步提升。在高温区，当电子浓度增加带来的电导率提升效应与电子散射增强导致的电导率降低效应达到平衡时，电导率将趋于稳定。能带的色散关系对电子的迁移率有着重要影响。在PbTe的能带结构中，导带底附近的能带色散较为平缓，这意味着电子在该区域的有效质量较大，迁移率较低。而在价带顶附近，能带色散相对较陡，电子的有效质量较小，迁移率较高。这种能带色散关系的差异使得在不同的载流子浓度条件下，电子的输运性质会有所不同。当以电子作为载流子（n型掺杂）时，由于导带电子迁移率较低，电导率的提升相对较为困难；而当以空穴作为载流子（p型掺杂）时，价带空穴迁移率较高，更有利于提高电导率。通过合理的掺杂调控，可以改变载流子的类型和浓度，进而优化PbTe的电输运性质。3.1.2电子态密度计算与分析在完成PbTe能带结构计算的基础上，进一步运用相同的计算方法和参数设置，对其电子态密度进行深入计算。图2展示了PbTe的总态密度（TDOS）和分波态密度（PDOS）计算结果。从总态密度图中可以看出，在费米能级（EF）附近，态密度存在明显的峰值，这表明在该能量区域电子态较为丰富，电子参与输运的可能性较大。在价带区域，态密度主要由Te的p轨道贡献，而在导带区域，态密度则主要由Pb的s轨道和Te的p轨道共同贡献。这种电子态密度的分布特征与PbTe的化学键性质密切相关，Pb-Te键具有一定的离子性和共价性，使得电子在不同原子轨道之间的分布呈现出特定的规律。通过分析分波态密度，可以更清晰地了解各个原子轨道对电子态密度的贡献。在价带顶附近，Te的p轨道态密度占据主导地位，这是因为Te原子的电负性相对较大，其p轨道电子更容易受到外界因素的影响。当在PbTe中引入掺杂原子时，掺杂原子的电子态会与Pb和Te的电子态相互作用，从而改变电子态密度的分布。例如，当引入施主杂质（如I）时，杂质原子会向导带提供额外的电子，使得导带中的电子态密度增加，费米能级向导带移动；而当引入受主杂质（如Na）时，杂质原子会在价带中产生空穴，增加价带中的空穴态密度，费米能级向价带移动。这种费米能级的移动和电子态密度的变化会显著影响PbTe的电导率和塞贝克系数等输运性质。在热电材料中，态密度在费米能级处的陡峭变化（高dg(E)/dE）有助于提升塞贝克系数。从PbTe的态密度图中可以观察到，在费米能级附近，态密度的变化较为陡峭，这为PbTe具有一定的热电性能提供了电子结构基础。通过调整材料的成分、晶体结构或引入缺陷等方式，可以进一步优化态密度在费米能级处的分布，从而提高塞贝克系数，增强热电性能。3.2电导率与塞贝克系数计算3.2.1电导率计算方法与结果电导率是描述材料导电能力的重要物理量，对于热电化合物的性能评估具有关键意义。在本研究中，采用基于玻尔兹曼输运理论的方法来计算热电化合物的电导率。该方法基于第一性原理计算得到的电子结构信息，通过求解玻尔兹曼输运方程，能够准确地计算出材料在不同条件下的电导率。玻尔兹曼输运方程的一般形式为：\frac{\partialf}{\partialt}+\vec{v}\cdot\nabla_{\vec{r}}f+\vec{F}\cdot\nabla_{\vec{k}}f=\left(\frac{\partialf}{\partialt}\right)_{coll}，其中f为电子分布函数，\vec{v}为电子速度，\vec{r}为位置矢量，\vec{k}为波矢，\vec{F}为外力，\left(\frac{\partialf}{\partialt}\right)_{coll}为碰撞项。在计算电导率时，通常假设电子处于稳态，即\frac{\partialf}{\partialt}=0，并且忽略电子-电子相互作用，将碰撞项简化为弛豫时间近似，即\left(\frac{\partialf}{\partialt}\right)_{coll}=-\frac{f-f_0}{\tau}，其中f_0为平衡态电子分布函数，\tau为弛豫时间。通过这些近似处理，结合第一性原理计算得到的能带结构和态密度等信息，可以计算出电导率的表达式为：\sigma=\frac{e^2}{V}\sum_{\vec{k}}\tau(\vec{k})\frac{\partialf_0}{\partialE}\left|\frac{\partialE}{\partial\vec{k}}\right|^2，其中e为电子电荷，V为晶体体积，E为电子能量。以硫化铅（PbS）为例，展示电导率的计算结果。运用基于密度泛函理论的第一性原理计算方法，结合上述电导率计算方法，对不同温度和载流子浓度下的PbS电导率进行了详细计算。计算结果如图3所示。从图中可以看出，在低温范围内，电导率随温度的升高而逐渐增大。这是因为随着温度的升高，热激发产生的载流子数量增加，同时电子的迁移率也有所提高，从而导致电导率上升。在高温区域，电导率随温度的升高而逐渐降低。这是由于高温下晶格振动加剧，电子与声子的散射增强，使得电子的迁移率显著下降，尽管载流子浓度可能仍在增加，但迁移率的降低对电导率的影响更为显著，导致电导率最终下降。载流子浓度对电导率也有着重要影响。随着载流子浓度的增加，电导率呈现出明显的上升趋势。这是因为载流子数量的增多使得参与导电的粒子数增加，从而提高了材料的导电能力。在实际应用中，可以通过掺杂等手段来调控载流子浓度，进而优化热电化合物的电导率。当在PbS中掺入适量的施主杂质（如Cl）时，杂质原子会向导带提供额外的电子，增加载流子浓度，从而提高电导率；而当掺入受主杂质（如Li）时，会在价带中产生空穴，同样可以改变载流子浓度，对电导率产生影响。除了温度和载流子浓度外，晶体结构、缺陷等因素也会对电导率产生影响。不同的晶体结构会导致电子的散射机制和迁移率不同，从而影响电导率。例如，具有紧密堆积结构的晶体，电子的散射相对较少，迁移率较高，电导率也相对较大；而具有疏松结构的晶体，电子散射较多，迁移率较低，电导率则较小。缺陷的存在会引入额外的散射中心，降低电子的迁移率，从而降低电导率。点缺陷（如空位、间隙原子）和线缺陷（如位错）都会对电子的输运产生阻碍作用。在材料制备过程中，通过优化制备工艺，减少缺陷的产生，可以有效地提高电导率。3.2.2塞贝克系数计算方法与结果塞贝克系数是衡量热电材料将热能转化为电能能力的关键参数，其大小直接反映了材料的热电转换效率。本研究采用基于线性响应理论的方法来计算塞贝克系数。该方法基于第一性原理计算得到的电子结构信息，通过求解线性响应方程，能够准确地计算出材料的塞贝克系数。塞贝克系数的计算公式为：S=-\frac{\pi^2k_B^2T}{3e}\frac{\left.\frac{\partial\sigma(E)}{\partialE}\right|_{E=E_F}}{\sigma(E_F)}，其中k_B为玻尔兹曼常数，T为绝对温度，e为电子电荷，\sigma(E)为能量为E时的电导率，E_F为费米能级。从公式中可以看出，塞贝克系数与电导率随能量的变化率以及费米能级处的电导率密切相关。以碲化铋（Bi2Te3）为例，展示塞贝克系数的计算结果。运用基于密度泛函理论的第一性原理计算方法，结合上述塞贝克系数计算方法，对不同温度和载流子浓度下的Bi2Te3塞贝克系数进行了精确计算。计算结果如图4所示。从图中可以看出，塞贝克系数与温度和载流子浓度呈现出明显的依赖关系。在低温区域，塞贝克系数随温度的升高而逐渐增大。这是因为在低温下，电子的能量分布较为集中，随着温度的升高，电子的能量分布展宽，更多的电子能够参与热电转换过程，从而导致塞贝克系数增大。在高温区域，塞贝克系数随温度的升高而逐渐减小。这是由于高温下电子的散射增强，电导率随能量的变化率减小，根据塞贝克系数的计算公式，其值也会相应减小。载流子浓度对塞贝克系数的影响也十分显著。随着载流子浓度的增加，塞贝克系数呈现出先增大后减小的趋势。当载流子浓度较低时，增加载流子浓度可以提高电导率随能量的变化率，从而增大塞贝克系数；然而，当载流子浓度过高时，电导率随能量的变化率减小，同时费米能级处的电导率增大，根据塞贝克系数的计算公式，其值会逐渐减小。在实际应用中，需要通过精确调控载流子浓度，使塞贝克系数达到最佳值，以提高热电材料的性能。塞贝克系数与材料的电子结构密切相关。材料的能带结构、态密度分布等因素都会影响塞贝克系数的大小。在Bi2Te3中，其独特的层状晶体结构和电子能带结构使得在费米能级附近具有合适的态密度分布和能带色散关系，有利于提高塞贝克系数。通过调整材料的成分、晶体结构或引入缺陷等方式，可以改变电子结构，进而优化塞贝克系数。在Bi2Te3中掺入适量的Sb或Se等元素，可以改变其电子结构，调整载流子浓度和电导率随能量的变化率，从而提高塞贝克系数。3.3电子输运性质的影响因素3.3.1晶体结构的影响晶体结构是决定热电化合物电子输运性质的重要因素之一，不同的晶体结构会导致电子在材料中的运动方式和散射机制存在显著差异，进而对电导率、塞贝克系数等电子输运性质产生重要影响。以方钴矿结构和笼目结构的热电化合物为例进行分析。方钴矿结构的热电化合物（如CoSb3）具有较为规整的晶体结构，其原子排列呈现出一定的周期性。在这种结构中，电子的散射主要来源于声子散射和杂质散射。由于晶体结构的周期性较好，声子的传播较为有序，电子与声子的散射相对较弱。然而，当材料中存在杂质原子时，杂质原子会破坏晶体结构的周期性，引入额外的散射中心，从而增强电子的散射，降低电导率。方钴矿结构的能带结构相对较为简单，在费米能级附近的态密度分布和能带色散关系对塞贝克系数有着重要影响。通过合理的掺杂调控，可以改变费米能级的位置和态密度分布，进而优化塞贝克系数。笼目结构的热电化合物（如NaZn13）具有独特的三维开放框架结构，这种结构中存在较大的空隙和通道。电子在这种结构中的输运过程中，除了受到声子散射和杂质散射外，还会受到结构中缺陷和空隙的散射。由于结构的开放性，电子在空隙和通道中传播时，散射几率相对较大，这会导致电导率降低。然而，笼目结构的特殊电子结构使得在费米能级附近具有较高的态密度变化率，这有利于提高塞贝克系数。笼目结构中的原子振动模式较为复杂，声子的散射机制也与传统晶体结构不同，这对热导率也会产生影响。通过对笼目结构的优化和调控，如引入合适的客体原子填充空隙，可以改变电子的散射机制和晶体的振动特性，从而实现对电子输运性质的有效调控。晶体结构的对称性也会对电子输运性质产生影响。具有较高对称性的晶体结构，电子在各个方向上的输运性质较为一致，电导率和塞贝克系数的各向异性较小；而具有较低对称性的晶体结构，电子在不同方向上的散射机制和迁移率可能存在较大差异，导致电导率和塞贝克系数具有明显的各向异性。在一些层状晶体结构的热电化合物中，电子在层内和层间的输运性质往往不同，层内的电导率和塞贝克系数可能与层间存在显著差异。这种各向异性的电子输运性质在实际应用中需要加以考虑，通过合理的材料设计和制备工艺，可以充分利用各向异性特性，优化热电材料的性能。3.3.2掺杂的影响掺杂是调控热电化合物电子输运性质的重要手段之一，通过在材料中引入杂质原子，可以有效地改变材料的电子结构、载流子浓度和散射机制，从而实现对电导率、塞贝克系数等电子输运性质的精确调控。以碲化铅（PbTe）为例，研究掺杂对其电子输运性质的影响。当在PbTe中掺入施主杂质（如铊Tl）时，Tl原子会取代Pb原子的位置，由于Tl的价电子数比Pb多，会向导带提供额外的电子，从而增加载流子浓度。随着载流子浓度的增加，电导率显著提高。如图5所示，未掺杂的PbTe电导率较低，而掺入适量Tl后的PbTe电导率明显增大。同时，掺杂也会对塞贝克系数产生影响。由于载流子浓度的增加，费米能级向导带移动，根据塞贝克系数的计算公式，其值会发生变化。在一定的掺杂浓度范围内，塞贝克系数可能会先增大后减小。这是因为在掺杂初期，载流子浓度的增加使得电导率随能量的变化率增大，从而导致塞贝克系数增大；然而，当掺杂浓度过高时，电导率随能量的变化率减小，同时费米能级处的电导率增大，根据塞贝克系数的计算公式，其值会逐渐减小。当在PbTe中掺入受主杂质（如钠Na）时，Na原子会取代Pb原子的位置，由于Na的价电子数比Pb少，会在价带中产生空穴，从而改变载流子类型和浓度。与施主掺杂类似，受主掺杂也会对电导率和塞贝克系数产生影响。随着空穴浓度的增加，电导率会发生变化，同时塞贝克系数也会相应改变。在这种情况下，塞贝克系数的变化趋势与施主掺杂时有所不同，具体取决于空穴浓度的变化以及价带结构的改变。除了改变载流子浓度外，掺杂还会影响材料的散射机制。杂质原子的引入会破坏晶体结构的周期性，形成额外的散射中心，从而增强电子的散射。这种散射的增强会降低电子的迁移率，进而对电导率产生负面影响。然而，在某些情况下，适当的掺杂可以通过引入特定的散射机制，如共振散射，来优化电子输运性质。在一些热电材料中，通过掺杂引入与主体材料具有特定电子结构的杂质原子，可以实现共振散射，在一定程度上提高塞贝克系数，同时保持电导率在可接受的范围内。3.3.3温度的影响温度是影响热电化合物电子输运性质的关键因素之一，随着温度的变化，材料的电子结构、载流子浓度、散射机制等都会发生改变，从而导致电导率、塞贝克系数等电子输运性质呈现出不同的变化规律。在低温范围内，热电化合物的电导率主要受载流子浓度和电子迁移率的影响。由于低温下热激发能量较低，载流子主要来源于杂质电离。随着温度的升高，杂质电离程度逐渐增加，载流子浓度随之增大，电导率也逐渐提高。此时，电子的散射主要四、热电化合物的晶格热导率与热输运性质4.1声子色散与声子态密度4.1.1声子色散计算与分析以硫化铅（PbS）这一典型的热电化合物为研究对象，运用基于密度泛函理论（DFT）的第一性原理计算方法，借助VASP软件对其声子色散进行精确计算。在计算过程中，采用有限位移法来获取力常数矩阵，进而得到声子色散关系。为了确保计算结果的准确性，平面波截断能设置为500eV，k点网格采用Monkhorst-Pack方法生成，设置为8×8×8，在倒易空间中进行细致采样。图6展示了PbS的声子色散计算结果。从图中可以清晰地观察到，PbS的声子色散曲线呈现出丰富的特征。在低频区域，声子色散曲线主要由声学支主导，声学支声子代表了晶体整体的振动，其频率与波矢呈近似线性关系。随着波矢的增加，声学支声子的频率逐渐增大，这表明晶体的振动模式逐渐变得复杂。在高频区域，声子色散曲线主要由光学支主导，光学支声子对应原子相对于彼此的振动，其频率相对较高，且与波矢的关系较为复杂。声子色散曲线与晶格振动密切相关。晶格振动是晶体中原子在其平衡位置附近的微小热振动，而声子则是晶格振动的量子化表现。不同的声子模式对应着不同的晶格振动方式，声学支声子对应晶体的整体平移振动，而光学支声子对应原子间的相对振动。在PbS的晶格中，Pb原子和S原子通过离子键相互作用，形成稳定的晶格结构。当晶体受到外界扰动时，原子会偏离其平衡位置，产生晶格振动。声子色散曲线反映了晶格振动的频率与波矢之间的关系，通过分析声子色散曲线，可以深入了解晶格振动的特性。在PbS的声子色散曲线中，存在一些特殊的点和区域。在布里渊区的中心（Γ点），声学支声子的频率为零，这对应着晶体的静止状态。随着波矢向布里渊区边界移动，声学支声子的频率逐渐增大，在布里渊区边界处达到最大值。光学支声子在Γ点附近具有较高的频率，且随着波矢的变化，光学支声子的频率也会发生明显的变化。这些特殊的点和区域反映了PbS晶格振动的对称性和周期性。声子色散曲线还受到晶体结构、原子质量和原子间相互作用等因素的影响。PbS具有立方晶体结构，这种结构的对称性较高，使得声子色散曲线在不同方向上具有一定的对称性。原子质量的差异也会对声子色散曲线产生影响，Pb原子的质量较大，S原子的质量较小，这种质量差异导致了光学支声子和声学支声子的频率分布不同。原子间的相互作用强度也会影响声子色散曲线，较强的原子间相互作用会使声子的频率升高，而较弱的原子间相互作用会使声子的频率降低。4.1.2声子态密度计算与分析在完成PbS声子色散计算的基础上，进一步运用相同的计算方法和参数设置，对其声子态密度进行深入计算。声子态密度（PhononDensityofStates，PDOS）是指单位频率范围内的声子模式数，它反映了声子在不同频率下的分布情况。图7展示了PbS的声子态密度计算结果。从图中可以看出，PbS的声子态密度分布呈现出明显的特征。在低频区域，声子态密度主要由声学支声子贡献，声学支声子的态密度随着频率的增加而逐渐增大。在高频区域，声子态密度主要由光学支声子贡献，光学支声子的态密度在一定频率范围内呈现出峰值。声子态密度分布对晶格热导率有着重要影响。晶格热导率主要由声子的传输贡献，而声子的传输与声子态密度密切相关。在低频区域，声学支声子的态密度较大，声子的平均自由程较长，这使得声学支声子在热传输过程中起着重要作用。然而，在高频区域，光学支声子的态密度虽然也较大，但由于光学支声子的频率较高，声子之间的散射较强，导致光学支声子的平均自由程较短，对晶格热导率的贡献相对较小。在PbS中，声子态密度的分布还与原子的振动模式有关。通过分析分波声子态密度，可以更清晰地了解各个原子对声子态密度的贡献。在低频区域，Pb原子和S原子的振动模式较为相似，它们对声学支声子态密度的贡献较为接近。在高频区域，S原子的振动模式相对较为活跃，其对光学支声子态密度的贡献较大。当在PbS中引入杂质或缺陷时，声子态密度的分布会发生变化。杂质或缺陷的存在会破坏晶体的周期性结构，引入额外的散射中心，从而改变声子的散射机制和平均自由程。杂质原子的质量和原子间相互作用与主体原子不同，会导致声子态密度在某些频率范围内发生变化。这种变化会影响声子的传输，进而对晶格热导率产生影响。在PbS中掺入少量的Cl原子，Cl原子的质量和原子间相互作用与Pb和S原子不同，会使声子态密度在某些频率范围内发生变化，导致声子的散射增强，平均自由程减小，从而降低晶格热导率。4.2晶格热导率计算4.2.1计算方法与原理晶格热导率是描述材料中热量传导能力的重要物理量，它对于理解热电化合物的热输运性质至关重要。本研究采用基于Boltzmann输运方程的方法来计算热电化合物的晶格热导率。Boltzmann输运方程是描述非平衡态下粒子输运过程的基本方程，其一般形式为：\frac{\partialf}{\partialt}+\vec{v}\cdot\nabla_{\vec{r}}f+\vec{F}\cdot\nabla_{\vec{k}}f=\left(\frac{\partialf}{\partialt}\right)_{coll}，其中f为粒子分布函数，\vec{v}为粒子速度，\vec{r}为位置矢量，\vec{k}为波矢，\vec{F}为外力，\left(\frac{\partialf}{\partialt}\right)_{coll}为碰撞项。在计算晶格热导率时，将声子视为粒子，通过求解Boltzmann输运方程来得到声子的输运性质，进而计算出晶格热导率。在基于Boltzmann输运方程计算晶格热导率时，通常采用一些近似方法来简化计算。常用的近似方法包括弛豫时间近似和三声子散射近似。弛豫时间近似假设声子的散射过程可以用一个弛豫时间来描述，即声子在两次散射之间的平均自由时间。在弛豫时间近似下，Boltzmann输运方程可以简化为：f-f_0=-\tau(\vec{k})\vec{v}\cdot\nabla_{\vec{r}}f_0，其中f_0为平衡态声子分布函数，\tau(\vec{k})为弛豫时间。三声子散射近似假设声子之间的散射主要是三声子过程，即一个声子与另外两个声子相互作用，产生新的声子。在三声子散射近似下，可以通过计算三声子散射率来得到声子的弛豫时间。通过求解Boltzmann输运方程，结合上述近似方法，可以得到晶格热导率的表达式为：\kappa_{l}=\frac{1}{3}\int_{BZ}C_{v}(\omega)\vec{v}_{g}(\vec{k})\vec{v}_{g}(\vec{k})\tau(\vec{k})d\vec{k}，其中\kappa_{l}为晶格热导率，C_{v}(\omega)为声子比热，\vec{v}_{g}(\vec{k})为声子群速度，\tau(\vec{k})为弛豫时间，积分是在整个布里渊区进行。在实际计算中，需要首先通过第一性原理计算得到材料的力常数矩阵，进而计算出声子的频率、群速度和声子态密度等信息。利用这些信息，结合上述晶格热导率的计算公式，就可以计算出材料的晶格热导率。在VASP软件中，可以通过有限位移法计算力常数矩阵，然后利用Phonopy等软件计算声子的相关性质，并最终计算出晶格热导率。4.2.2计算结果与分析以碲化铋（Bi2Te3）为例，展示晶格热导率的计算结果。运用基于密度泛函理论的第一性原理计算方法，结合基于Boltzmann输运方程的晶格热导率计算方法，对不同温度下的Bi2Te3晶格热导率进行了详细计算。计算结果如图8所示。从图中可以看出，Bi2Te3的晶格热导率随温度的变化呈现出明显的规律。在低温区域，晶格热导率随着温度的升高而逐渐增大。这是因为在低温下，声子的平均自由程较大，声子之间的散射较弱，随着温度的升高，声子的能量增加，声子的平均自由程略有减小，但声子的比热增加更为显著，导致晶格热导率增大。在高温区域，晶格热导率随着温度的升高而逐渐减小。这是由于高温下声子的平均自由程减小，声子之间的散射增强，尽管声子的比热仍然随着温度的升高而增加，但散射增强对晶格热导率的降低作用更为明显，导致晶格热导率最终减小。除了温度，材料的晶体结构、原子质量和原子间相互作用等因素也会对晶格热导率产生重要影响。Bi2Te3具有层状晶体结构，这种结构使得声子在层内和层间的传输特性不同。层内原子间的相互作用较强，声子在层内的平均自由程较大，传输速度较快，因此层内的晶格热导率相对较大；而层间原子间的相互作用较弱，声子在层间的散射较强，平均自由程较小，传输速度较慢，导致层间的晶格热导率相对较小。原子质量的差异也会影响晶格热导率。Bi和Te原子的质量不同，这种质量差异导致声子的振动模式和散射机制不同。较重的原子振动频率较低，声子的平均自由程较短，散射较强，从而降低了晶格热导率。而较轻的原子振动频率较高，声子的平均自由程较长，散射较弱，有利于提高晶格热导率。在Bi2Te3中，Te原子的质量相对较大，这在一定程度上增加了声子的散射，降低了晶格热导率。原子间相互作用的强度也对晶格热导率有显著影响。较强的原子间相互作用会使声子的频率升高，平均自由程增大，从而提高晶格热导率；而较弱的原子间相互作用会使声子的频率降低，平均自由程减小，导致晶格热导率降低。在Bi2Te3中，层内原子间的相互作用较强，有利于声子的传输，而层间原子间的相互作用较弱，限制了声子的传输，这也是导致层内和层间晶格热导率差异的重要原因之一。4.3热输运性质的影响因素4.3.1原子质量与键强的影响原子质量和键强是影响热电化合物热输运性质的重要因素，它们通过改变声子的振动特性和散射机制，对晶格热导率产生显著影响。以硒化铅（PbSe）为例，分析原子质量对热输运性质的影响。PbSe由铅（Pb）和硒（Se）原子组成，Pb原子的质量相对较大，Se原子的质量相对较小。原子质量的差异导致声子的振动特性不同。较重的Pb原子振动频率较低，声子的能量较低，平均自由程较短。这是因为质量较大的原子在相同的力作用下，加速度较小，振动速度较慢，从而使得声子的频率降低。较低的声子频率导致声子之间的散射增强，平均自由程减小，进而降低了晶格热导率。较轻的Se原子振动频率较高，声子的能量较高，平均自由程相对较长。然而，由于PbSe中Pb原子的质量占比较大，整体上原子质量的差异使得晶格热导率受到抑制。键强对热输运性质也有着重要影响。在PbSe中，Pb-Se键具有一定的离子性和共价性，键强适中。较强的键强意味着原子间的相互作用较强，原子在平衡位置附近的振动较为困难，声子的频率较高。高频率的声子具有较高的能量，在传输过程中不容易被散射，平均自由程较大，有利于提高晶格热导率。相反，较弱的键强使得原子间的相互作用较弱，原子容易偏离平衡位置，声子的频率较低。低频率的声子能量较低，容易与其他声子或晶体中的缺陷发生散射，平均自由程减小，导致晶格热导率降低。在一些具有较弱化学键的热电化合物中，如某些有机热电材料，由于键强较弱，晶格热导率通常较低。原子质量和键强还会相互影响，共同作用于热输运性质。在一些化合物中，原子质量较大的同时，键强也可能较弱。在这种情况下，较重原子导致的声子频率降低和较弱键强导致的声子散射增强会相互叠加，进一步降低晶格热导率。而在另一些化合物中，原子质量较小且键强较强时，声子的频率较高且散射较弱，有利于提高晶格热导率。4.3.2晶格缺陷与杂质的影响晶格缺陷和杂质是影响热电化合物热输运性质的关键因素，它们通过改变晶体的结构和电子状态，对声子的散射和输运产生重要影响，进而改变晶格热导率。晶格缺陷包括点缺陷（如空位、间隙原子）、线缺陷（如位错）和面缺陷（如晶界）等。以碲化铋（Bi2Te3）为例，当晶体中存在空位缺陷时，空位周围的原子会发生弛豫，导致局部晶格结构发生畸变。这种畸变会破坏声子传播的周期性，使得声子在传播过程中与空位发生散射。声子与空位的散射概率与空位浓度和空位的大小有关。空位浓度越高，声子与空位相遇的概率越大，散射越强，平均自由程减小，晶格热导率降低。空位的大小也会影响散射强度，较大的空位会引起更大的晶格畸变，对声子的散射作用更强。位错是晶体中的线缺陷，它是由于晶体在生长或受力过程中原子排列的不规则而产生的。位错周围存在着应力场，声子在传播过程中会与位错发生相互作用。位错对声子的散射机制较为复杂，包括弹性散射和非弹性散射。弹性散射是指声子与位错相互作用后，声子的能量不变，只是传播方向发生改变；非弹性散射则是声子与位错相互作用后，声子的能量发生变化。位错的密度和类型会影响声子的散射强度。高密度的位错会提供更多的散射中心，增强声子的散射，降低晶格热导率。不同类型的位错（如刃型位错、螺型位错）对声子的散射作用也有所不同。晶界是晶体中的面缺陷，它是不同晶粒之间的界面。在Bi2Te3多晶材料中，晶界的存在会导致声子在晶界处发生散射。晶界处原子排列不规则，原子间的相互作用与晶粒内部不同，这使得声子在跨越晶界时会受到强烈的散射。晶界对声子的散射作用与晶界的结构、取向和晶界能等因素有关。具有高晶界能的晶界通常对声子的散射更强。通过控制晶界的结构和取向，可以调节晶界对声子的散射作用，从而实现对晶格热导率的调控。在纳米结构的Bi2Te3材料中，由于晶界数量增多，声子与晶界的散射增强，晶格热导率显著降低。杂质原子的引入也会对热电化合物的热输运性质产生重要影响。当在Bi2Te3中掺入杂质原子时，杂质原子的质量和电子结构与主体原子不同，会导致局部晶格结构和电子状态发生改变。杂质原子的质量差异会引起声子的散射，较重或较轻的杂质原子会使声子的振动频率发生变化，从而增加声子的散射概率。杂质原子的电子结构差异会影响电子-声子相互作用，进而影响声子的散射和输运。在Bi2Te3中掺入少量的Sb原子，Sb原子的电子结构与Bi和Te原子不同，会改变电子-声子相互作用，导致声子的散射增强，晶格热导率降低。4.3.3温度的影响温度是影响热电化合物热输运性质的关键因素之一，随着温度的变化，材料的晶格振动特性、声子散射机制等都会发生改变，从而导致晶格热导率呈现出不同的变化规律。以硅锗（SiGe）合金为例，在低温区域（通常低于德拜温度的1/3），晶格热导率随着温度的升高而逐渐增大。这是因为在低温下，声子的平均自由程主要受限于晶体中的杂质和缺陷散射。随着温度的升高，声子的能量增加，声子与杂质和缺陷的散射概率略有减小，平均自由程增大。同时，声子的比热随着温度的升高而逐渐增大。根据五、热电化合物输运性质的调控策略5.1结构调控5.1.1纳米结构化纳米结构化是调控热电化合物输运性质的一种重要策略，通过将材料制备成纳米结构，可以显著改变其物理性质，实现对热电性能的优化。当热电化合物被制成纳米结构时，材料的比表面积大幅增加，界面数量显著增多，这些界面能够有效地散射声子，从而降低晶格热导率。以碲化铋（Bi2Te3）基纳米材料为例，研究表明，纳米晶界对声子的散射作用比对电子的散射作用更强，这是因为声子的波长与纳米晶界的尺寸在同一数量级，声子更容易与晶界发生相互作用而被散射。而电子的平均自由程相对较短，在通过纳米晶界时，虽然也会受到一定程度的散射，但相比声子，其受到的影响较小。这种对声子和电子散射的差异性，使得在降低晶格热导率的同时，能够保持相对较高的电导率。纳米结构化还可以通过量子限域效应来影响电子的输运性质。在纳米尺度下，电子的运动受到限制，能级发生量子化，导致电子态密度分布发生变化。这种变化可以增强电子的散射，从而提高塞贝克系数。在纳米线结构的热电材料中，由于电子在一维方向上的运动受到限制，电子态密度在某些能量范围内出现尖锐的峰值，使得塞贝克系数显著提高。量子限域效应还可以改变电子的有效质量，进而影响电导率。当电子的有效质量增加时，电导率会降低；反之，电导率会增加。通过合理设计纳米结构的尺寸和形状，可以调控量子限域效应的强度，实现对电导率和塞贝克系数的协同优化。纳米结构化对热电化合物的热稳定性也有重要影响。由于纳米结构的比表面积大，表面原子的活性较高，在高温下，表面原子的扩散和反应速度加快，可能导致材料的结构和性能发生变化。因此，在制备纳米结构热电材料时，需要采取适当的措施来提高其热稳定性。可以通过表面修饰、包覆等方法，降低表面原子的活性，抑制表面原子的扩散和反应。也可以通过优化纳米结构的设计，提高材料的结构稳定性，减少高温下结构变化对热电性能的影响。5.1.2合金化合金化是调控热电化合物输运性质的另一种重要手段，通过将不同元素合金化，可以改变材料的晶体结构、电子结构和原子间相互作用，从而实现对热电性能的有效调控。合金化可以通过引入质量和尺寸差异来增强声子散射，降低晶格热导率。在硅锗（SiGe）合金中，Ge原子的质量比Si原子大，原子半径也略有不同。当Ge原子掺入Si晶格中时，会引起晶格的局部畸变，形成质量和尺寸波动。这种波动会散射声子，使得声子的平均自由程减小，从而降低晶格热导率。合金化还可以改变材料的晶体结构对称性，进一步影响声子的散射机制。在一些合金体系中，合金化会导致晶体结构从简单的立方结构转变为复杂的多晶结构，晶界和位错等缺陷增多，这些缺陷会提供更多的声子散射中心，从而更有效地降低晶格热导率。合金化还可以调节材料的电子结构，优化电导率和塞贝克系数。在某些热电化合物中，通过合金化引入特定的元素，可以改变材料的能带结构，调整费米能级的位置。在碲化铅（PbTe）中掺入铊（Tl），Tl原子的价电子数比Pb多，会向导带提供额外的电子，使费米能级向导带移动，增加载流子浓度，从而提高电导率。同时，合金化还可能改变能带的色散关系和态密度分布，影响电子的迁移率和塞贝克系数。当合金化导致能带结构发生变化，使电子的有效质量减小，迁移率增大时，电导率会进一步提高；而当合金化使费米能级附近的态密度分布发生改变，增加了电导率随能量的变化率时，塞贝克系数会增大。合金化过程中的原子相互作用也会对热电性能产生影响。不同元素之间的原子相互作用会改变化学键的性质和强度，从而影响电子和声子的输运。在一些合金体系中，原子之间可能形成共价键、离子键或金属键等不同类型的化学键，这些化学键的性质和强度会影响电子的离域程度和声子的振动特性。较强的化学键可能使电子更难移动，降低电导率，但同时也可能使声子的频率升高，增强声子的散射，降低晶格热导率。因此，在合金化过程中，需要综合考虑原子相互作用对电子和声子输运的影响，通过合理选择合金元素和控制合金成分，实现对热电性能的优化。5.2掺杂调控5.2.1施主与受主掺杂施主与受主掺杂是调控热电化合物电子结构和输运性质的重要手段，通过引入不同类型的杂质原子，可以有效地改变材料的载流子浓度和能带结构，从而对热电性能产生显著影响。当在热电化合物中引入施主杂质时，施主原子会向导带提供额外的电子，增加载流子浓度。以碲化铋（Bi2Te3）为例，当掺入铟（In）作为施主杂质时，In原子取代Bi原子的位置，由于In的价电子数比Bi多，会向导带释放一个额外的电子，使导带中的电子浓度增加。随着载流子浓度的增加，电导率显著提高。如图9所示，未掺杂的Bi2Te3电导率较低，而掺入适量In后的Bi2Te3电导率明显增大。施主掺杂还会影响塞贝克系数。由于载流子浓度的增加，费米能级向导带移动，根据塞贝克系数的计算公式，其值会发生变化。在一定的掺杂浓度范围内，塞贝克系数可能会先增大后减小。这是因为在掺杂初期，载流子浓度的增加使得电导率随能量的变化率增大，从而导致塞贝克系数增大；然而，当掺杂浓度过高时，电导率随能量的变化率减小，同时费米能级处的电导率增大，根据塞贝克系数的计算公式，其值会逐渐减小。受主掺杂则是在热电化合物中引入能够接受电子的杂质原子，从而在价带中产生空穴，改变载流子类型和浓度。在Bi2Te3中掺入钠（Na）作为受主杂质，Na原子取代Bi原子的位置，由于Na的价电子数比Bi少，会在价带中产生一个空穴，使价带中的空穴浓度增加。与施主掺杂类似，受主掺杂也会对电导率和塞贝克系数产生影响。随着空穴浓度的增加，电导率会发生变化，同时塞贝克系数也会相应改变。在这种情况下，塞贝克系数的变化趋势与施主掺杂时有所不同，具体取决于空穴浓度的变化以及价带结构的改变。施主与受主掺杂还会影响材料的散射机制。杂质原子的引入会破坏晶体结构的周期性，形成额外的散射中心，从而增强电子的散射。这种散射的增强会降低电子的迁移率，进而对电导率产生负面影响。然而，在某些情况下，适当的掺杂可以通过引入特定的散射机制，如共振散射，来优化电子输运性质。在一些热电材料中，通过掺杂引入与主体材料具有特定电子结构的杂质原子，可以实现共振散射，在一定程度上提高塞贝克系数，同时保持电导率在可接受的范围内。5.2.2共掺杂与协同效应共掺杂是指在热电化合物中同时引入两种或多种不同类型的杂质原子，通过它们之间的协同作用来调控材料的输运性质，实现对热电性能的进一步优化。共掺杂体系中存在着复杂的相互作用，这些相互作用不仅影响载流子的浓度和迁移率，还会改变材料的能带结构和声子散射机制，从而对热电性能产生独特的影响。以碲化铅（PbTe）为例，研究发现，当同时掺入铊（Tl）和钠（Na）进行共掺杂时，会产生显著的协同效应。Tl作为施主杂质，向导带提供额外的电子，增加电子浓度；而Na作为受主杂质，在价带中产生空穴。这两种杂质的引入，一方面通过改变载流子浓度，对电导率和塞贝克系数产生影响。在一定的共掺杂浓度范围内，由于电子和空穴浓度的协同变化，电导率和塞贝克系数可以得到同时优化。另一方面，共掺杂还会引入复杂的缺陷结构和应力场。Tl和Na原子的尺寸与Pb原子不同，它们的掺入导致晶格发生畸变，形成应力场。这种应力场会散射声子，降低晶格热导率。共掺杂引入的缺陷结构还可能改变电子的散射机制，进一步影响电导率和塞贝克系数。在这种共掺杂体系中，可能会出现一些特殊的散射机制，如缺陷共振散射，它可以在一定程度上提高塞贝克系数，同时对电导率的影响较小。共掺杂还可以通过调节杂质原子之间的相互作用，实现对能带结构的精细调控。在一些共掺杂体系中，不同杂质原子的电子态会相互耦合，形成新的能级结构。在PbTe中同时掺入铟（In）和锑（Sb），In和Sb的电子态会相互作用，在能带中引入新的能级。这些新能级的出现可以改变电子的跃迁路径和概率，从而影响电导率和塞贝克系数。通过合理选择共掺杂的杂质原子和控制它们的浓度比例，可以实现对能带结构的精确调控，优化热电性能。共掺杂体系中的协同效应还体现在对材料稳定性的影响上。一些共掺杂组合可以增强材料的结构稳定性，抑制杂质原子的偏析和扩散。在某些热电化合物中，单一掺杂时杂质原子容易在晶界或表面偏析，导致材料性能下降。而通过共掺杂，不同杂质原子之间的相互作用可以抑制这种偏析现象，提高材料的稳定性。共掺杂还可以改善材料的化学稳定性，增强材料对环境因素的耐受性。在一些易氧化的热电材料中，通过共掺杂引入具有抗氧化性的杂质原子，可以提高材料的抗氧化能力，延长其使用寿命。5.3外场调控5.3.1电场调控电场调控是一种有效的调控热电化合物输运性质的方法，通过施加外部电场，可以改变材料内部的电子分布和运动状态，进而影响热电性能。当在热电化合物上施加电场时，电场会对电子产生作用力，改变电子的运动轨迹和能量分布。以碲化铋（Bi2Te3）为例，在电场的作用下，电子会受到一个与电场方向相反的力，从而加速运动。这种加速运动会增加电子的迁移率，进而提高电导率。根据欧姆定律，电导率与电子迁移率和载流子浓度成正比，在载流子浓度不变的情况下，电子迁移率的增加会导致电导率增大。如图10所示，随着施加电场强度的增加，Bi2Te3的电导率逐渐增大。电场还会影响塞贝克系数。在热电材料中，塞贝克系数与电子的能量分布和散射机制密切相关。施加电场后，电子的能量分布会发生变化，导致电导率随能量的变化率改变，从而影响塞贝克系数。在一些情况下，电场可以使电子的能量分布更加均匀，减小电导率随能量的变化率，从而降低塞贝克系数。而在另一些情况下，电场可以增强电子的散射，增加电导率随能量的变化率，从而提高塞贝克系数。具体的影响取决于材料的电子结构和电场强度。电场调控还可以通过改变材料的能带结构来影响热电性能。在电场的作用下，材料的能带会发生弯曲，费米能级的位置也会发生变化。这种能带结构的改变会影响电子的跃迁概率和散射机制，进而影响电导率和塞贝克系数。在一些半导体热电材料中，施加电场可以使能带弯曲，形成内建电场，从而增强电子的散射，降低电导率。而在一些具有特殊能带结构的材料中，电场可以改变能带的对称性和色散关系，优化电子的输运性质，提高热电性能。5.3.2磁场调控磁场调控是研究热电化合物输运性质的另一个重要方向，磁场的施加会对材料中的电子产生洛伦兹力，从而影响电子的运动轨迹和能量分布，进而对热电性能产生影响。此外，磁场还会引发磁电耦合效应，进一步调控热电性能。当在热电化合物中施加磁场时，电子在磁场中运动时会受到洛伦兹力的作用，其运动轨迹会发生弯曲，形成螺旋状的运动路径。这种运动轨迹的改变会增加电子的散射概率，从而影响电导率。以碲化铅（PbTe）为例，在磁场的作用下，电子的散射增强，电导率降低。磁场对塞贝克系数也有影响。由于电子运动轨迹的改变，电子的能量分布和散射机制发生变化，导致塞贝克系数发生改变。在某些情况下，磁场可以使电子的能量分布更加不均匀，增加电导率随能量的变化率，从而提高塞贝克系数；而在另一些情况下，磁场会使电子的散射增强，降低电子的迁移率，从而降低塞贝克系数。具体的影响取决于材料的电子结构和磁场强度。磁场还会引发磁电耦合效应，这是指磁场和电场之间的相互作用对热电性能的影响。在一些具有特殊晶体结构的热电材料中，如某些氧化物热电材料，磁场的施加会改变材料的晶体结构对称性，导致电极化强度发生变化。这种电极化强度的变化会影响电子的输运性质，进而影响热电性能。在一些铁电-铁磁复合材料中，磁场和电场的共同作用可以产生协同效应，增强热电性能。通过调节磁场和电场的强度和方向，可以实现对热电性能的有效调控。磁电耦合效应还可以通过改变材料的能带结构来影响热电性能。在磁场和电场的作用下，材料的能带会发生复杂的变化，费米能级的位置和电子态密度分布也会改变。这些变化会影响电子的跃迁概率和散射机制，从而对电导率和塞贝克系数产生影响。在一些拓扑材料中，磁场和电场的作用可以激发拓扑表面态，这些表面态具有独特的电子输运性质，能够显著提高热电性能。通过研究磁电耦合效应，深入了解磁场对热电化合物输运性质的影响机制，为热电材料的性能优化提供了新的途径。六、案例分析与应用前景6.1典型热电化合物的输运性质研究6.1.1IV-VIA族化合物IV-VIA族化合物作为一类重要的热电材料，展现出独特的输运性质，受到了广泛的研究关注。以碲化铅（PbTe）为例，其具有面心立方晶体结构，在中高温区（约500-800K）展现出优异的热电性能。运用第一性原理计算方法，深入分析其输运性质及影响因素，能够为该类化合物的性能优化和应用提供重要的理论依据。从能带结构来看，PbTe具有典型的半导体能带结构，价带顶位于布里渊区的L点，导带底位于Γ点，存在明显的能带间隙。这种能带结构特性使得PbTe在一定条件下能够实现电子的激发和跃迁，从而参与电输运过程。在低温状态下，由于热激发能量较低，电子主要分布在价带中，难以跨越能带间隙进入导带，因此电导率较低。随着温度逐渐升高，热激发能量增大，更多的电子获得足够的能量从价带跃迁到导带，导带中的电子浓度增加，电导率随之提高。然而，温度的升高也会导致电子的散射增强，这在一定程度上会限制电导率的进一步提升。在高温区，当电子浓度增加带来的电导率提升效应与电子散射增强导致的电导率降低效应达到平衡时，电导率将趋于稳定。电子态密度分析进一步揭示了PbTe的电子特性。在费米能级（EF）附近，态密度存在明显的峰值，表明在该能量区域电子态较为丰富，电子参与输运的可能性较大。在价带区域，态密度主要由Te的p轨道贡献，而在导带区域，态密度则主要由Pb的s轨道和Te的p轨道共同贡献。这种电子态密度的分布特征与PbTe的化学键性质密切相关，Pb-Te键具有一定的离子性和共价性，使得电子在不同原子轨道之间的分布呈现出特定的规律。当在PbTe中引入掺杂原子时，掺杂原子的电子态会与Pb和Te的电子态相互作用，从而改变电子态密度的分布。引入施主杂质（如I）时，杂质原子会向导带提供额外的电子，使得导带中的电子态密度增加，费米能级向导带移动；而引入受主杂质（如Na）时，杂质原子会在价带中产生空穴，增加价带中的空穴态密度，费米能级向价带移动。这种费米能级的移动和电子态密度的变化会显著影响PbTe的电导率和塞贝克系数等输运性质。除了电子结构对输运性质的影响，晶体结构、掺杂和温度等因素也起着关键作用。PbTe的面心立方