基于第一性原理的CuMO₂体系电子结构与热电性质深度剖析

基于第一性原理的CuMO₂体系电子结构与热电性质深度剖析一、绪论1.1研究背景与意义随着全球经济的快速发展，能源消耗与日俱增，传统化石能源的日益枯竭以及其使用带来的环境污染问题，如温室气体排放导致的全球气候变暖、酸雨等，已成为人类社会可持续发展面临的严峻挑战。开发高效、清洁的新能源技术和提高能源利用效率迫在眉睫。热电材料作为一种能够实现热能和电能直接相互转换的功能材料，在能源领域展现出巨大的应用潜力，引起了广泛关注。热电材料的工作原理基于塞贝克效应（Seebeckeffect）、帕尔贴效应（Peltiereffect）和汤姆逊效应（Thomsoneffect）。其中，塞贝克效应是指当两种不同的导体或半导体组成的回路中存在温度梯度时，会产生热电势，从而实现热能到电能的转换，这是热电发电的基础；帕尔贴效应则是当电流通过两种不同导体或半导体的界面时，会在界面处产生吸热或放热现象，可用于实现固态制冷；汤姆逊效应是指在存在温度梯度的单一导体或半导体中，当电流通过时，导体会吸收或放出热量。这些效应使得热电材料在能源转换领域具有独特的优势，例如可利用工业废热、汽车尾气余热等低品位热源进行发电，提高能源利用效率，降低环境污染；也可通过帕尔贴效应实现固态制冷，具有无噪音、无振动、易于维护等优点，适用于精密仪器、电子设备等场景的制冷需求。目前，虽然已有一些热电材料得到了应用，如碲化铋基材料在低温制冷领域、铅基材料在中温废热回收领域、硅基材料在高温工业废热回收和太阳能热电转换等领域有一定应用，但现有热电材料的热电转换效率普遍较低，限制了其大规模应用。因此，开发新型高性能热电材料成为该领域的研究重点。在众多新型热电材料的研究中，CuMO₂（M为过渡金属元素，如Cr、Fe、Mn等）体系由于其独特的晶体结构和电子结构，展现出潜在的优良热电性能，成为研究热点之一。CuMO₂体系属于铜铁矿结构氧化物，其晶体结构具有准二维层状超晶格结构，在这种结构中，八面体[MO₆]层和六方Cu原子层沿c轴交替排列，层间距离较大形成开放式的层状结构。这种结构特点使得该体系在电学和热学性能方面具有独特的优势，例如较大的层间距离可能会影响电子的传输和散射，从而对电导率和Seebeck系数产生影响；同时，开放式的层状结构也可能会影响声子的传播，进而影响热导率。通过对CuMO₂体系电子结构和热电性质的深入研究，有望揭示其内在的物理机制，为进一步优化其热电性能提供理论依据，从而推动高性能热电材料的发展，在能源转换领域发挥更大的作用，缓解当前的能源和环境危机。1.2热电材料研究现状热电材料的研究历史可以追溯到19世纪初，1821年德国物理学家ThomasSeebeck发现了塞贝克效应，1834年JeanPeltier发现了帕尔贴效应，1851年WilliamThomson发现了汤姆逊效应，这些效应的发现为热电材料的研究奠定了理论基础。此后，热电材料的研究经历了多个阶段的发展。早期主要围绕金属材料进行研究，但由于金属材料的热电转换效率很低，研究进展缓慢。直到20世纪中叶，随着半导体材料的开发和应用，研究人员发现了Seebeck系数大于100μV/K的半导体材料，热电材料的研究进入了一个新的高潮。经过长期的研究与发展，目前常见的热电材料主要包括以下几类：碲化铋基热电材料：主要由碲化铋（Bi₂Te₃）及其固溶体组成。在低温范围（通常低于200℃），碲化铋基材料具有优异的热电性能，其热电优值ZT可达到1左右。这得益于其独特的层状晶体结构，层间的弱相互作用使得电子传输相对容易，同时对声子有较强的散射作用，降低了热导率。由于其良好的热电性能，碲化铋基材料被广泛应用于温差发电和电子冷却领域，如在小型冰箱、饮水机、电子设备散热等方面有实际应用。例如，在一些便携式冷藏设备中，碲化铋基热电制冷模块可以实现高效的制冷效果，且具有无噪音、无振动的优点。铅基热电材料：以铅（Pb）为主要元素，通过添加其他元素形成固溶体或化合物。在中温范围（200-600℃），铅基热电材料展现出较高的热电优值。例如，PbTe是一种典型的铅基热电材料，通过合适的掺杂和微观结构调控，其ZT值在中温区能达到较高水平。铅基材料成本相对较低，制备工艺也较为成熟。在汽车尾气余热利用领域，铅基热电材料可将尾气中的废热转化为电能，为汽车电子设备供电，提高能源利用效率；在工业废热回收方面，也有一定的应用潜力。硅基热电材料：以硅（Si）为主要元素，通过掺杂其他元素或形成化合物来改善其热电性能。硅基热电材料在高温范围（高于600℃）具有较好的热电性能，同时具备机械强度高、化学稳定性好的优点。在工业废热回收中，硅基热电材料可用于高温热源的能量转换；在太阳能热电转换领域，可与太阳能电池结合，进一步提高太阳能的利用效率。例如，在一些高温工业炉的余热回收系统中，硅基热电材料制成的发电模块能够稳定运行，实现废热的有效利用。氧化物热电材料：包括铜铁矿结构氧化物（如CuMO₂）、钙钛矿结构氧化物、钴酸盐类氧化物等。氧化物热电材料具有使用温度高、不怕氧化、无污染、使用寿命长、制备方便等优点。以NaCo₂O₄为代表的钴酸盐类氧化物在中温区有良好的热电性能；而铜铁矿结构氧化物CuMO₂由于其独特的晶体结构和电子结构，在热电性能方面展现出潜在的优势，成为研究热点之一。氧化物热电材料在中高温区的热电发电领域具有重要的应用潜力，尤其是在对环境友好和稳定性要求较高的场合。方钴矿型热电材料：通式为AB₃，其中A是金属元素（如Ir、Co、Rh、Fe等），B是V族元素（如P、As、Sb等）。方钴矿化合物具有立方晶系晶体结构，每个晶胞内有两个较大的笼状孔隙。这种特殊的结构使得方钴矿型热电材料可以通过在笼状孔隙中填充其他原子来调控其热电性能。填充原子可以散射声子，降低热导率，同时对电子传输影响较小。通过合理的填充和掺杂，方钴矿型热电材料在中高温区能获得较高的ZT值，在废热回收和空间电源等领域具有潜在的应用价值。虽然热电材料的研究取得了一定进展，但目前热电材料的热电转换效率仍有待提高，成本也相对较高，限制了其大规模应用。为了突破这些瓶颈，科研人员不断探索新的材料体系和制备工艺，如纳米结构热电材料、复合材料、低维热电材料等，通过纳米效应、界面散射等机制来优化热电性能；同时，借助第一性原理计算等理论方法，深入理解热电材料的电子结构和输运机制，为材料的设计和优化提供理论指导。1.3CuMO₂体系热电材料研究状况CuMO₂体系热电材料属于铜铁矿结构氧化物，其晶体结构具有独特的准二维层状超晶格结构。在这种结构中，八面体[MO₆]层和六方Cu原子层沿c轴交替排列。层间距离较大，形成开放式的层状结构，这一结构特点对其电学和热学性能有着重要影响。较大的层间距离使得电子在传输过程中可能会受到不同程度的散射，从而影响电导率；同时，这种开放式结构对声子的散射机制也与其他结构不同，可能导致热导率发生变化，这些因素都为其热电性能的研究提供了丰富的物理内涵。在热电性能研究方面，研究人员通过实验和理论计算等多种方法对CuMO₂体系进行了深入探索。在实验方面，研究发现通过对CuMO₂体系进行元素掺杂可以有效地调控其热电性能。例如，对CuFeO₂进行Zn掺杂，研究发现适量的Zn掺杂可以在不显著降低电导率的情况下，提高Seebeck系数，从而提高功率因子。这是因为Zn的掺入改变了体系的电子结构，调整了载流子浓度和迁移率，进而影响了热电性能。在热导率方面，研究发现通过引入纳米结构或控制微观结构，可以增加声子散射，降低热导率。制备纳米晶CuCrO₂材料，其晶格热导率明显低于常规块体材料，这是由于纳米晶界对声子有强烈的散射作用。在理论计算方面，第一性原理计算在研究CuMO₂体系热电性质中发挥了重要作用。通过第一性原理计算可以深入了解体系的电子结构，如能带结构、态密度等，从而揭示热电性能的内在物理机制。计算CuFeO₂的电子结构发现，其能带结构具有一定的特点，费米能级附近的电子态分布对电导率和Seebeck系数有重要影响。通过理论计算还可以预测不同掺杂或结构变化对热电性能的影响，为实验研究提供理论指导。理论计算预测在CuCrO₂中进行特定元素的掺杂可以提高其热电性能，后续的实验研究在一定程度上验证了这一预测。在应用方面，CuMO₂体系热电材料由于其潜在的优良热电性能，在多个领域展现出应用潜力。在温差发电领域，利用其在温度梯度下实现热能到电能的转换特性，可用于回收工业废热、汽车尾气余热等低品位热源的能量，提高能源利用效率。在一些工业生产过程中，会产生大量的废热，若能利用CuMO₂体系热电材料将这些废热转化为电能，不仅可以减少能源浪费，还能降低生产成本。在固态制冷领域，基于帕尔贴效应，CuMO₂体系热电材料也有潜在的应用前景，可用于小型制冷设备、电子设备的散热等场景。由于其具有无噪音、无振动、易于维护等优点，在对制冷设备体积和稳定性要求较高的场合具有独特优势。尽管CuMO₂体系热电材料的研究取得了一定进展，但仍面临一些挑战。目前该体系热电材料的热电转换效率与实际应用需求相比仍有差距，需要进一步优化材料的性能。在制备工艺方面，现有的制备方法可能存在成本较高、制备过程复杂等问题，限制了其大规模生产和应用。因此，未来的研究需要在提高热电性能和改进制备工艺等方面深入开展，以推动CuMO₂体系热电材料从实验室研究走向实际应用。1.4铜铁矿氧化物结构铜铁矿氧化物具有独特的晶体结构，其化学式通常表示为ABO₂，其中A为一价金属离子（如Cu⁺、Ag⁺等），B为三价金属离子（如Cr³⁺、Fe³⁺、Mn³⁺等）。其晶体结构属于六方晶系，空间群为P6₃/mmc，具有准二维层状超晶格结构。在这种结构中，八面体[BO₆]层和六方A原子层沿c轴交替排列。具体来说，B原子位于八面体的中心，被六个O原子包围形成[BO₆]八面体，这些八面体通过共边或共角的方式连接成层状结构；而A原子则位于六方密堆积的位置，与[BO₆]层交替排列，形成沿c轴方向的层状结构。层间距离较大，形成开放式的层状结构，这种结构特点使得铜铁矿氧化物在电子和热电性质方面表现出独特的性能。这种晶体结构对电子性质有着显著的影响。较大的层间距离使得电子在层间的传输受到一定阻碍，导致电子的散射增加，从而影响电导率。由于A原子和B原子的电子结构不同，以及[BO₆]八面体的存在，使得铜铁矿氧化物的能带结构较为复杂。在费米能级附近，电子态密度的分布与晶体结构密切相关。对于CuFeO₂，[FeO₆]八面体中的Fe原子的3d电子与O原子的2p电子之间存在较强的相互作用，这种相互作用会影响电子的跃迁和散射，进而影响电导率和Seebeck系数。同时，晶体结构中的缺陷和杂质也会对电子性质产生影响，如氧空位的存在会改变电子的分布和传输路径，从而改变材料的电学性能。在热电性质方面，铜铁矿氧化物的晶体结构同样起着关键作用。热导率是衡量热电材料性能的重要参数之一，它与声子的传播密切相关。铜铁矿氧化物的开放式层状结构对声子的散射机制与其他结构不同。较大的层间距离使得声子在层间传播时容易受到散射，从而降低了热导率。[BO₆]八面体的振动模式也会影响声子的散射，这些八面体的振动可以与声子相互作用，增加声子的散射几率，进一步降低热导率。在Seebeck系数方面，晶体结构决定了电子的分布和传输特性，进而影响Seebeck系数。不同的A原子和B原子组合会导致电子结构的差异，从而使得Seebeck系数发生变化。对于CuCrO₂和CuFeO₂，由于Cr和Fe的电子结构不同，它们的Seebeck系数也有所不同。通过对晶体结构的调控，如掺杂、形成固溶体等，可以改变电子和热学性能，从而优化热电性能。1.5掺杂对铜铁矿结构氧化物材料性能影响研究进展掺杂是调控铜铁矿结构氧化物材料性能的重要手段，通过引入不同种类和浓度的杂质原子，可以改变材料的晶体结构、电子结构，进而对其电学、热学和光学等性能产生显著影响。在电学性能方面，掺杂对电导率和Seebeck系数有着重要影响。当在铜铁矿结构氧化物中引入施主或受主杂质时，会改变材料的载流子浓度。在CuFeO₂中掺入Zn²⁺，由于Zn²⁺的价态与Cu⁺不同，会导致电子的得失，从而改变载流子浓度。适量的Zn掺杂可以在一定程度上提高电导率。这是因为掺杂引入的额外载流子增加了载流子浓度，使得在电场作用下参与导电的载流子增多。当掺杂浓度过高时，杂质原子可能会形成散射中心，增加电子的散射几率，从而导致电导率下降。在Seebeck系数方面，掺杂会改变材料的电子结构，影响载流子的散射机制和能量分布。通过合理的掺杂，可以优化载流子的散射，使得Seebeck系数得到提高。研究发现，在CuCrO₂中进行特定元素的掺杂，能够调整费米能级附近的电子态密度，改变载流子的散射过程，进而提高Seebeck系数。热学性能也是掺杂影响的重要方面，其中热导率是关键参数。铜铁矿结构氧化物的热导率主要由晶格热导率和电子热导率组成。掺杂可以通过多种机制降低晶格热导率。一方面，掺杂原子与基质原子的尺寸差异会产生晶格畸变，这种畸变会散射声子，增加声子的散射几率，从而降低晶格热导率。在CuMnO₂中掺入尺寸不同的杂质原子，会导致晶格局部结构的改变，形成声子散射中心，有效降低晶格热导率。另一方面，掺杂还可以引入点缺陷、位错等微观缺陷，这些缺陷也会对声子产生散射作用。通过引入纳米级的第二相粒子，与基体形成界面，界面处的声子散射也能降低热导率。在制备过程中，控制掺杂元素的分布和微观结构，如形成纳米复合结构，可以进一步增强对声子的散射，从而更有效地降低热导率。在光学性能方面，掺杂对铜铁矿结构氧化物的光学带隙和光吸收特性有明显影响。由于掺杂原子的电子结构与基质原子不同，会改变材料的电子能级分布，进而影响光学带隙。在AgFeO₂中进行掺杂，研究发现掺杂可以调控其光学带隙。合适的掺杂可以使光学带隙变窄或变宽，这取决于掺杂原子的种类和浓度。当光学带隙发生变化时，材料对光的吸收特性也会相应改变。如果光学带隙变窄，材料对长波长光的吸收能力会增强；反之，光学带隙变宽则会使材料对短波长光的吸收更明显。这种光学性能的调控使得铜铁矿结构氧化物在光电器件领域具有潜在的应用价值，如在光探测器、发光二极管等器件中的应用。理论研究在理解掺杂对铜铁矿结构氧化物材料性能影响方面发挥了重要作用。第一性原理计算可以从原子尺度深入研究掺杂体系的电子结构、电荷分布和化学键性质等。通过计算不同掺杂体系的能带结构、态密度等，可以揭示掺杂对电子结构的影响机制，从而预测材料的电学、热学和光学性能。理论计算还可以研究掺杂原子与基质原子之间的相互作用，以及掺杂引起的晶格畸变对性能的影响。计算结果为实验研究提供了理论指导，帮助实验人员设计合理的掺杂方案，优化材料性能。例如，通过理论计算预测在CuFeO₂中进行某种特定元素的掺杂可以提高热电性能，实验人员可以根据这一预测进行实验验证，大大提高了研究效率。1.6研究内容与创新点1.6.1研究内容CuMO₂体系电子结构的第一性原理计算：运用基于密度泛函理论的第一性原理计算方法，选用平面波赝势方法或全电子线性化缀加平面波方法等，对不同M元素（如Cr、Fe、Mn等）的CuMO₂体系进行电子结构计算。分析其能带结构，明确价带和导带的分布、带隙大小及特性，判断材料是金属、半导体还是绝缘体；研究态密度，包括总态密度和分波态密度，了解不同原子轨道对电子态的贡献，以及费米能级附近电子态的分布特征。例如，对于CuCrO₂，通过计算其能带结构，观察Cr的3d轨道与O的2p轨道形成的杂化轨道对带隙的影响；分析态密度，确定费米能级处主要由哪些原子轨道的电子占据，为后续热电性质的研究提供电子结构基础。CuMO₂体系热电性质的理论预测：基于计算得到的电子结构，利用玻尔兹曼输运理论，结合第一性原理计算结果，预测CuMO₂体系的热电性质。计算电导率，考虑电子的散射机制，如声子散射、杂质散射等对电导率的影响，分析不同散射机制在不同温度下的主导作用；计算Seebeck系数，探究其与电子结构的关系，如载流子浓度、费米能级位置等对Seebeck系数的影响规律；计算热导率，包括晶格热导率和电子热导率，研究晶格振动模式对晶格热导率的影响，以及电子-声子相互作用对总热导率的贡献。通过理论预测，获得不同温度和载流子浓度下的热电优值ZT，评估CuMO₂体系的热电性能。掺杂对CuMO₂体系电子结构和热电性质的影响研究：选择合适的掺杂元素，如在CuFeO₂中掺入Zn、Mg等元素，通过第一性原理计算研究掺杂对CuMO₂体系晶体结构的影响，包括晶格参数的变化、原子位置的偏移等。分析掺杂后体系的电子结构变化，如杂质能级的引入、能带结构的改变、态密度的重新分布等。研究掺杂对热电性质的影响机制，探讨如何通过掺杂优化电导率、Seebeck系数和热导率，以提高热电优值ZT。例如，计算不同掺杂浓度下的电学和热学输运性质，找出最佳的掺杂方案，为实验制备高性能的CuMO₂基热电材料提供理论指导。与实验结果对比验证及分析：收集已有的CuMO₂体系热电材料的实验数据，包括电子结构和热电性质方面的数据。将理论计算结果与实验数据进行对比分析，验证理论模型和计算方法的准确性。对于计算结果与实验数据存在差异的情况，深入分析原因，如计算中未考虑的因素（如实验中的缺陷、杂质分布不均匀等），或者实验条件的不确定性等。通过对比验证，进一步完善理论模型，提高对CuMO₂体系电子结构和热电性质的理解，为后续的理论研究和实验制备提供更可靠的依据。1.6.2创新点研究方法创新：本研究综合运用多种先进的第一性原理计算方法，如平面波赝势方法结合全电子线性化缀加平面波方法的高精度计算，全面深入地研究CuMO₂体系的电子结构，这种多方法联用的方式能够更准确地描述体系的电子行为，克服单一方法的局限性。在研究掺杂对热电性质的影响时，采用了基于机器学习辅助的第一性原理计算方法，通过构建大量的掺杂体系数据集并进行机器学习训练，快速筛选出具有潜在高性能的掺杂方案，大大提高了研究效率，这在CuMO₂体系热电材料研究中是一种创新性的尝试。研究视角创新：从原子尺度和电子结构层面深入探究CuMO₂体系热电性能的内在物理机制，不仅关注常见的电学和热学输运性质，还研究了电子关联效应、轨道杂化等微观因素对热电性能的影响，为理解该体系热电性能提供了全新的视角。首次考虑了CuMO₂体系中晶体结构的动态变化对热电性质的影响，通过分子动力学模拟与第一性原理计算相结合的方法，研究温度变化过程中晶体结构的演变及其对电子和声子输运的影响，这在以往的研究中较少涉及。研究内容创新：对多种不同M元素的CuMO₂体系进行系统研究，全面对比分析它们的电子结构和热电性质，揭示不同元素对体系性能的影响规律，为开发新型高性能热电材料提供更丰富的选择和理论依据，这种系统的研究在同类研究中具有创新性。探索了新型的掺杂策略，如多元素协同掺杂和梯度掺杂，研究其对CuMO₂体系电子结构和热电性质的协同调控作用，为进一步提高热电材料的性能开辟了新的途径。二、理论基础与计算方法2.1密度泛函理论密度泛函理论（DensityFunctionalTheory，DFT）是一种研究多电子体系电子结构的量子力学方法，在凝聚态物理、材料科学、化学等领域有着广泛的应用。它将体系的基态能量表示为电子密度的泛函，通过求解电子密度来确定体系的基态性质，从而大大简化了多电子体系的计算。在研究CuMO₂体系的电子结构和热电性质时，密度泛函理论为我们提供了重要的理论框架和计算工具。2.1.1绝热近似绝热近似，又称Born-Oppenheimer近似，由麦克斯・玻恩（MaxBorn）和弗拉基米尔・福克（VladimirFock）于1928年提出。其基本概念基于电子与原子核的质量和运动特性差异。电子质量极小，约为原子核质量的千分之一甚至更小，这使得电子的运动速度远快于原子核。在量子力学的框架下，对于一个包含电子和原子核的多体体系，当原子核发生微小移动时，由于电子运动速度极快，电子能够迅速调整其状态，以适应原子核位置的变化。从时间尺度来看，电子的运动时间尺度约为10^{-16}秒，而原子核的运动时间尺度约为10^{-13}秒，二者相差约三个数量级。这意味着在原子核运动的一个微小时间间隔内，电子已经完成了多次的状态调整，仿佛电子始终处于与原子核瞬时位置相对应的基态。基于上述特性，绝热近似在密度泛函理论中起着关键作用。它允许我们将电子与原子核的运动进行解耦。在研究多电子体系时，我们可以先固定原子核的位置，将原子核视为提供一个固定的外部势场。在这个固定势场下，电子的哈密顿量仅包含电子的动能、电子与电子之间的相互作用能以及电子与原子核之间的势能。通过求解这个仅涉及电子的哈密顿量，我们可以得到电子的波函数和能量。在处理CuMO₂体系时，我们首先将Cu、M（如Cr、Fe、Mn等）和O原子核的位置固定，构建电子在这些原子核形成的势场中的哈密顿量，进而求解电子的状态。当需要考虑原子核的运动时，如在研究晶格动力学或热膨胀等问题时，我们可以在电子结构计算的基础上，将电子结构计算得到的能量作为原子核运动的势能面，再利用经典力学或量子力学的方法来研究原子核在这个势能面上的运动。这种解耦的方式极大地简化了多体体系的计算，使得我们能够分别对电子和原子核的行为进行深入研究。2.1.2Hohenberg-Kohn定理Hohenberg-Kohn定理是密度泛函理论的重要基石，包含两个核心内容。Hohenberg-Kohn第一定理表明，对于一个处于外部势场v(r)中的多电子体系，其基态的所有物理性质都由基态电子密度n(r)唯一确定。这里的外部势场v(r)主要由原子核产生的库仑势构成。该定理意味着，一旦我们确定了体系的基态电子密度，那么体系的基态能量、电荷分布、电子云形状等所有物理性质都随之确定。从数学角度来看，体系的哈密顿量H可以表示为H=T+V_{ee}+\sum_{i}v(r_i)，其中T是电子的动能算符，V_{ee}是电子-电子相互作用能算符，\sum_{i}v(r_i)是外部势场。而体系的基态能量E_0是哈密顿量H在基态波函数\Psi_0下的期望值，即E_0=\langle\Psi_0|H|\Psi_0\rangle。根据第一定理，由于基态电子密度n(r)与基态波函数\Psi_0存在一一对应的关系（尽管这种对应关系在实际中很难直接给出解析表达式），所以基态电子密度n(r)能够唯一确定基态能量E_0以及其他物理性质。Hohenberg-Kohn第二定理指出，存在一个以基态电子密度n(r)为变量的普适泛函E[n(r)]，使得当电子密度取真实的基态电子密度时，该泛函达到最小值，且这个最小值就是体系的基态能量。这个普适泛函E[n(r)]可以表示为E[n(r)]=T[n(r)]+E_{ee}[n(r)]+E_{ext}[n(r)]，其中T[n(r)]是电子的动能泛函，E_{ee}[n(r)]是电子-电子相互作用能泛函，E_{ext}[n(r)]是外部势场作用能泛函。在实际应用中，虽然我们无法精确知道这些泛函的具体形式，但可以通过各种近似方法来求解。例如，在局域密度近似（LDA）中，假设体系中每个小区域内的电子行为类似于均匀电子气，从而对交换关联能泛函进行近似。Hohenberg-Kohn定理为密度泛函理论提供了坚实的理论基础。它打破了传统量子力学中必须通过求解多电子波函数来确定体系性质的局限，将多电子体系的研究转化为对电子密度的研究。电子密度是一个三维空间的函数，相比多电子波函数（其变量数为电子数的三倍，是一个非常高维的函数），电子密度的维度大大降低，使得计算在数值上变得可行。这使得我们能够通过对电子密度的变分求解，来确定体系的基态性质，为后续的Kohn-Sham方程的建立和实际计算提供了理论依据。在研究CuMO₂体系时，我们基于该定理，通过构建合适的电子密度泛函，来计算体系的电子结构和相关物理性质。2.1.3Kohn-Sham方程Kohn-Sham方程是密度泛函理论的核心方程，其推导基于Hohenberg-Kohn定理和单电子近似。首先，根据Hohenberg-Kohn定理，体系的基态能量是基态电子密度的泛函。为了求解这个泛函，Kohn和Sham提出了一个重要的近似，将多电子体系看作是一组无相互作用的电子在一个有效势场中运动。我们从多电子体系的哈密顿量开始推导。多电子体系的哈密顿量H可以表示为：H=-\frac{1}{2}\sum_{i=1}^{N}\nabla_{i}^{2}+\sum_{i=1}^{N}v_{ext}(r_{i})+\frac{1}{2}\sum_{i\neqj}^{N}\frac{1}{|r_{i}-r_{j}|}其中，-\frac{1}{2}\sum_{i=1}^{N}\nabla_{i}^{2}是电子的动能项，\sum_{i=1}^{N}v_{ext}(r_{i})是电子与外部势场（主要是原子核产生的势场）的相互作用项，\frac{1}{2}\sum_{i\neqj}^{N}\frac{1}{|r_{i}-r_{j}|}是电子-电子相互作用项。Kohn-Sham引入了一个虚构的无相互作用电子体系，其哈密顿量H_{KS}为：H_{KS}=-\frac{1}{2}\sum_{i=1}^{N}\nabla_{i}^{2}+\sum_{i=1}^{N}v_{eff}(r_{i})其中，v_{eff}(r_{i})是有效势场，它包含了外部势场v_{ext}(r_{i})、Hartree势v_{H}(r_{i})以及交换关联势v_{xc}(r_{i})，即v_{eff}(r_{i})=v_{ext}(r_{i})+v_{H}(r_{i})+v_{xc}(r_{i})。Hartree势v_{H}(r_{i})描述了电子与其他电子的平均库仑相互作用，其表达式为：v_{H}(r_{i})=\int\frac{n(r')}{|r_{i}-r'|}dr'其中，n(r')是电子密度。交换关联势v_{xc}(r_{i})则包含了电子之间的交换能和关联能的贡献，它是一个复杂的泛函，目前没有精确的解析表达式，通常采用各种近似方法来处理。对于这个虚构的无相互作用电子体系，其波函数可以用一组单电子波函数\psi_{i}(r)的Slater行列式来表示。根据变分原理，体系的基态能量E可以通过对波函数的变分来求解。对H_{KS}应用变分原理，得到Kohn-Sham方程：\left[-\frac{1}{2}\nabla^{2}+v_{eff}(r)\right]\psi_{i}(r)=\epsilon_{i}\psi_{i}(r)其中，\epsilon_{i}是单电子本征能量。在实际求解Kohn-Sham方程时，通常采用自洽迭代的方法。首先，给定一个初始的电子密度分布n^{(0)}(r)，根据上述公式计算出有效势场v_{eff}^{(0)}(r)。然后，将v_{eff}^{(0)}(r)代入Kohn-Sham方程，求解得到一组新的单电子波函数\psi_{i}^{(1)}(r)。根据这组新的波函数，可以计算出新的电子密度分布n^{(1)}(r)=\sum_{i=1}^{N}|\psi_{i}^{(1)}(r)|^{2}。接着，用n^{(1)}(r)重新计算有效势场v_{eff}^{(1)}(r)，再代入Kohn-Sham方程求解新的波函数，如此循环迭代，直到前后两次迭代得到的电子密度或能量的差异小于某个设定的收敛阈值，此时认为计算达到自洽，得到的结果即为体系的基态电子结构。在计算CuMO₂体系的电子结构时，Kohn-Sham方程发挥着关键作用。通过精确求解Kohn-Sham方程，我们可以得到体系中每个原子周围的电子分布情况，进而分析其能带结构、态密度等电子结构信息，为研究热电性质提供基础。例如，通过分析能带结构，我们可以确定体系的带隙大小，判断其是金属、半导体还是绝缘体，这对于理解体系的电学性质至关重要；通过研究态密度，我们可以了解不同原子轨道对电子态的贡献，以及费米能级附近电子态的分布特征，从而进一步分析体系的电子输运性质。2.1.4交换关联能交换关联能是密度泛函理论中一个重要的概念，它描述了电子之间除了经典库仑相互作用之外的量子力学相互作用。在多电子体系中，电子是全同费米子，由于泡利不相容原理，具有相同自旋的电子不能占据同一量子态，这导致电子之间存在一种特殊的相互作用，即交换能。电子之间还存在由于电子的瞬时位置相关性而产生的关联能，它反映了电子之间的动态相互作用。交换关联能E_{xc}是交换能E_{x}和关联能E_{c}之和，即E_{xc}=E_{x}+E_{c}。目前，由于交换关联能的精确泛函形式未知，在实际计算中通常采用近似方法。常见的近似方法包括局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）。局域密度近似假设体系中某点的交换关联能只与该点的电子密度有关，且等于具有相同密度的均匀电子气的交换关联能。其交换关联能泛函E_{xc}^{LDA}[n]的表达式为：E_{xc}^{LDA}[n]=\intn(r)\epsilon_{xc}^{LDA}(n(r))dr其中，\epsilon_{xc}^{LDA}(n(r))是均匀电子气的交换关联能密度，它是电子密度n(r)的函数。LDA在处理一些简单体系时取得了较好的结果，例如金属体系，它能够较为准确地描述电子的行为和体系的性质。对于一些具有明显非均匀电子密度分布的体系，如表面体系和分子体系，LDA的近似效果相对较差。这是因为LDA没有考虑电子密度的梯度信息，而在这些体系中，电子密度的变化较为剧烈，梯度信息对交换关联能的影响不可忽略。广义梯度近似在LDA的基础上，考虑了电子密度的梯度对交换关联能的影响。其交换关联能泛函E_{xc}^{GGA}[n]不仅依赖于电子密度n(r)，还依赖于电子密度的梯度\nablan(r)。常见的GGA泛函有Perdew-Burke-Ernzerhof（PBE）泛函等。以PBE泛函为例，其交换关联能泛函的形式较为复杂，它通过引入一些与电子密度和梯度相关的参数和函数来描述交换关联能。GGA在处理具有非均匀电子密度分布的体系时，通常比LDA表现得更好。在研究半导体材料时，GGA能够更准确地预测带隙大小，这对于理解半导体的电学性质和应用具有重要意义。GGA也存在一定的局限性，对于一些强关联体系，如过渡金属氧化物，GGA可能无法准确描述电子之间的强相互作用，导致计算结果与实验值存在较大偏差。交换关联能的近似方法对计算结果有着显著的影响。不同的近似方法会导致交换关联能的计算值不同，进而影响体系的总能量、电子结构和其他物理性质的计算结果。在研究CuMO₂体系时，选择合适的交换关联能近似方法至关重要。如果选择的近似方法不合适，可能会导致对体系带隙的计算出现较大偏差，从而影响对其电学性质的判断；在计算热导率等热电性质时，交换关联能的近似误差也可能会传播到最终的计算结果中，影响对热电性能的评估。因此，在实际研究中，需要根据体系的特点和研究目的，合理选择交换关联能的近似方法，并通过与实验结果或更精确的理论计算进行对比，来验证和优化计算结果。2.2全电子线性化缀加平面波法全电子线性化缀加平面波（Full-potentialLinearizedAugmentedPlaneWave，FLAPW）方法是一种基于密度泛函理论的高精度电子结构计算方法，在研究CuMO₂体系的电子结构和热电性质时具有重要应用。它能够精确地处理电子与原子核之间的相互作用，以及电子之间的交换关联作用，为深入理解材料的微观结构和物理性质提供了有力工具。2.2.1Muffin-tin球Muffin-tin球是FLAPW方法中的一个重要概念。在固体中，电子的分布可以分为两个区域：以每个原子核为中心，半径为R的球形区域，以及所有球与球之间的区域。这里的球形区域被称为Muffin-tin球区，球外的区域称为间隙区。芯态电子波函数完全在Muffin-tin球内，只有价态电子可以延伸到球外的间隙区。在Muffin-tin球内，通常取球对称势，球外则取常数势，这种势场模型称为Muffin-tin势。通过这种划分，能够分别对不同区域的电子行为进行合理描述。在计算CuMO₂体系时，Muffin-tin球内的势场主要由Cu、M（如Cr、Fe、Mn等）和O原子核产生，球内的电子行为接近自由原子中的电子行为。而球外间隙区的电子则可以用平面波更好地近似。这种划分方式使得在描述电子波函数时，能够充分利用原子电子波函数在近核区域的准确性和平面波在远离核区域的有效性。在计算CuCrO₂的电子结构时，通过合理定义Muffin-tin球，可以精确描述Cr原子周围的电子分布，以及电子在间隙区的传播，从而为准确计算体系的能带结构和态密度提供基础。2.2.2线性化缀加平面波的基函数线性化缀加平面波的基函数是FLAPW方法的关键组成部分。它的构造基于Muffin-tin势。在球内，基函数由径向波函数和球谐函数的乘积组成，其中径向波函数满足径向薛定谔方程；在球外，势场为零，解具有平面波的形式。通过将平面波展开成球谐函数和贝塞尔函数，并利用球面上波函数的连续条件，可以确定基函数中的系数。与传统的缀加平面波（APW）基函数不同，线性化缀加平面波（LAPW）基函数在Muffin-tin球内对APW基函数中的径向函数做了泰勒展开，使得径向薛定谔方程的解不再是能量本征值的函数，而是某一个选定的能量参数。这样做的好处是，确定能量中心后，便可以方便地获得这个能量附近的径向波函数，大大提高了计算速度。例如，在计算CuFeO₂的电子结构时，通过构建合适的LAPW基函数，能够准确描述Fe原子的3d电子和O原子的2p电子的波函数，进而精确计算体系的电子态密度和能带结构。这种基函数的构造方式使得FLAPW方法能够更准确地描述电子波函数，对于研究CuMO₂体系中复杂的电子结构和电子相互作用具有重要意义。2.2.3线性缀加平面波方法FLAPW方法的计算步骤基于密度泛函理论和线性化缀加平面波基函数。首先，需要对体系进行几何结构优化，确定原子的坐标和晶格参数。在研究CuMO₂体系时，通过优化结构，可以得到最稳定的晶体结构，为后续的电子结构计算提供准确的几何模型。然后，基于优化后的结构，构建Muffin-tin势，并确定线性化缀加平面波基函数。接下来，利用这些基函数展开电子波函数，并根据Kohn-Sham方程构建哈密顿矩阵和重叠矩阵。通过求解本征方程，得到体系的电子能量本征值和波函数。在求解过程中，通常采用自洽迭代的方法，不断更新电子密度和势场，直到计算结果收敛。FLAPW方法的原理在于将多电子体系的薛定谔方程转化为一组单电子的Kohn-Sham方程进行求解。通过线性化缀加平面波基函数的构建，能够精确描述电子在Muffin-tin球内和间隙区的行为，从而准确计算体系的电子结构。这种方法的优势在于它能够处理全电子体系，精确考虑电子与原子核之间的相互作用，对于过渡金属氧化物等含有d电子的体系，能够准确描述d电子的行为，得到高精度的电子结构和物理性质。在研究CuMO₂体系时，由于过渡金属M的d电子对体系的电学和热电性质有重要影响，FLAPW方法能够很好地处理这些d电子，为深入理解体系的性质提供准确的理论计算结果。FLAPW方法也存在一定的局限性。由于其计算精度高，计算量相对较大，对于大规模体系的计算效率较低。在处理复杂体系时，计算过程可能会变得非常复杂，需要较高的计算资源和较长的计算时间。2.2.4Wien2k软件Wien2k软件是一款基于FLAPW方法的电子结构计算软件，在材料科学领域有着广泛的应用。它的功能十分强大，能够计算晶体材料的电子结构，包括能带结构、态密度、电荷密度等。通过计算能带结构，可以了解材料中电子的能量分布和运动状态，判断材料是金属、半导体还是绝缘体；计算态密度可以分析不同原子轨道对电子态的贡献，以及费米能级附近电子态的分布特征；计算电荷密度可以直观地看到电子在原子周围的分布情况，研究原子之间的化学键性质。在研究CuMO₂体系时，Wien2k软件可以准确计算其电子结构，为研究热电性质提供基础。Wien2k软件具有一些显著的特点。它采用全势方法，能够精确考虑原子势场的非球形部分，这对于准确描述电子与原子核之间的相互作用非常重要。软件的计算精度高，能够得到与实验结果较为吻合的理论计算结果。它还提供了丰富的输入输出选项，方便用户进行各种参数设置和结果分析。使用Wien2k软件进行计算时，通常需要遵循一定的流程。首先，需要准备输入文件，包括晶体结构信息、原子坐标、基函数设置等。在研究CuMO₂体系时，需要准确输入其晶体结构参数，如晶格常数、原子的坐标等。然后，运行软件进行计算，在计算过程中，软件会根据输入参数进行自洽迭代计算，直到满足收敛条件。计算完成后，需要对输出结果进行分析。用户可以通过软件自带的分析工具，或者将输出数据导入其他软件进行处理，得到所需的物理量和图像，如绘制能带结构和态密度图，分析电荷密度分布等。通过这些结果分析，可以深入了解CuMO₂体系的电子结构和物理性质。2.3玻尔兹曼理论在研究CuMO₂体系的热电性质时，玻尔兹曼理论起着至关重要的作用。它为我们理解电子在材料中的输运过程提供了坚实的理论基础，通过玻尔兹曼方程可以定量地描述电子的输运性质，如电导率、Seebeck系数和热导率等。同时，BoltzTraP软件作为基于玻尔兹曼理论的计算工具，能够高效地计算这些热电输运性质，为我们深入研究CuMO₂体系的热电性能提供了有力的支持。2.3.1玻尔兹曼方程玻尔兹曼方程是描述非平衡态统计物理中粒子输运过程的重要方程，它在研究材料的热电性质时起着关键作用。在固体中，电子的输运受到多种因素的影响，如声子散射、杂质散射等。为了推导玻尔兹曼方程，我们考虑一个处于外场（如电场、温度梯度等）中的多电子体系。假设电子的分布函数为f(\vec{r},\vec{k},t)，它表示在时刻t，位置为\vec{r}，波矢为\vec{k}的电子数密度。根据刘维尔定理，在没有散射的情况下，分布函数沿着相空间的轨道是守恒的。但在实际的固体中，电子会与声子、杂质等发生散射，从而改变其运动状态。考虑散射过程后，分布函数f(\vec{r},\vec{k},t)的变化率可以表示为：\frac{\partialf}{\partialt}=-\vec{v}\cdot\nabla_{\vec{r}}f-\frac{e\vec{E}}{\hbar}\cdot\nabla_{\vec{k}}f+\left(\frac{\partialf}{\partialt}\right)_{scatt}其中，\vec{v}是电子的速度，\vec{v}=\frac{1}{\hbar}\nabla_{\vec{k}}\epsilon(\vec{k})，\epsilon(\vec{k})是电子的能量；\vec{E}是外电场；\left(\frac{\partialf}{\partialt}\right)_{scatt}是散射项，表示由于散射导致的分布函数的变化率。对于稳态情况，即\frac{\partialf}{\partialt}=0，上式可简化为：\vec{v}\cdot\nabla_{\vec{r}}f+\frac{e\vec{E}}{\hbar}\cdot\nabla_{\vec{k}}f=\left(\frac{\partialf}{\partialt}\right)_{scatt}这就是玻尔兹曼方程的一般形式。在研究热电性质时，通常考虑线性响应情况，即外场较弱，分布函数f相对于平衡态分布函数f_0的偏离较小。此时，可将f表示为f=f_0+f_1，其中f_1是偏离平衡态的部分，且|f_1|\ll|f_0|。将f=f_0+f_1代入玻尔兹曼方程，并忽略高阶项，得到线性化的玻尔兹曼方程：\vec{v}\cdot\nabla_{\vec{r}}f_0+\frac{e\vec{E}}{\hbar}\cdot\nabla_{\vec{k}}f_0=-\left(\frac{\partialf_1}{\partialt}\right)_{scatt}玻尔兹曼方程的物理意义在于，它描述了在各种外场和散射机制作用下，电子分布函数的变化规律。方程左边第一项\vec{v}\cdot\nabla_{\vec{r}}f表示由于电子的漂移运动导致的分布函数在空间上的变化；第二项\frac{e\vec{E}}{\hbar}\cdot\nabla_{\vec{k}}f表示外电场对电子运动的影响，即电子在电场作用下的加速运动导致的分布函数在波矢空间上的变化。方程右边\left(\frac{\partialf}{\partialt}\right)_{scatt}则体现了散射过程对电子分布函数的改变，散射使得电子的运动状态发生随机变化，从而影响电子的输运。在输运性质计算中，玻尔兹曼方程是核心方程。通过求解玻尔兹曼方程，可以得到电子的分布函数f，进而计算出各种输运性质。对于电导率\sigma，可根据电流密度\vec{J}与电场\vec{E}的关系\vec{J}=\sigma\vec{E}，利用分布函数f计算出电流密度\vec{J}，从而得到电导率\sigma。在计算Seebeck系数时，考虑体系存在温度梯度\nablaT，通过玻尔兹曼方程求解分布函数，进而得到Seebeck系数。对于热导率\kappa，也可以基于玻尔兹曼方程，结合能量流密度的计算来得到。玻尔兹曼方程为我们从微观层面理解电子的输运过程提供了理论框架，是研究材料热电性质的重要工具。在研究CuMO₂体系时，通过求解玻尔兹曼方程，可以深入了解电子在该体系中的输运特性，为优化其热电性能提供理论依据。例如，通过分析散射项对分布函数的影响，可以探讨如何通过调控散射机制来提高电导率或降低热导率，从而提高热电优值。2.3.2BoltzTraP软件BoltzTraP软件是一款专门用于计算晶体材料热电输运性质的软件，它基于玻尔兹曼理论，结合第一性原理计算得到的电子结构信息，能够高效准确地计算材料的电导率、Seebeck系数、热导率等热电性质。BoltzTraP软件的功能十分强大。它可以读取多种第一性原理计算软件输出的电子结构数据，如VASP、WIEN2k等软件的计算结果。通过读取这些数据，BoltzTraP软件能够获取材料的能带结构、态密度等信息，为后续的热电性质计算提供基础。在计算电导率时，软件根据玻尔兹曼方程，考虑电子的散射机制，计算电子在外电场作用下的输运情况，从而得到电导率。在计算Seebeck系数时，软件考虑体系存在温度梯度时电子的输运特性，通过求解玻尔兹曼方程得到Seebeck系数。对于热导率的计算，BoltzTraP软件分别考虑晶格热导率和电子热导率，综合计算得到总热导率。使用BoltzTraP软件进行热电性质计算时，通常遵循一定的流程。首先，需要使用第一性原理计算软件对材料的电子结构进行计算，得到相关的输出文件。在研究CuMO₂体系时，我们可以使用VASP软件进行结构优化和电子结构计算，得到包含能带结构、态密度等信息的输出文件。然后，将这些输出文件转换为BoltzTraP软件能够读取的格式。有些第一性原理计算软件与BoltzTraP软件有直接的接口，可以方便地进行数据转换；对于没有直接接口的软件，可能需要使用一些中间工具进行格式转换。接下来，设置BoltzTraP软件的计算参数，如温度范围、载流子浓度等。根据研究的需求，合理设置这些参数，以得到准确的计算结果。运行BoltzTraP软件进行计算，计算完成后，软件会输出各种热电性质的计算结果，包括不同温度和载流子浓度下的电导率、Seebeck系数、热导率等。用户可以对这些结果进行分析和处理，例如绘制热电性质随温度或载流子浓度变化的曲线，分析热电性能的变化规律。在热电性质计算中，BoltzTraP软件具有重要的作用。它将复杂的理论计算过程进行了封装，使得研究人员可以方便地利用第一性原理计算结果来研究材料的热电性质。通过BoltzTraP软件，我们可以快速地得到不同材料或不同掺杂体系的热电性质，为材料的筛选和优化提供依据。在研究CuMO₂体系时，使用BoltzTraP软件可以系统地研究不同M元素的CuMO₂体系以及掺杂后的体系的热电性质，深入分析结构与性能之间的关系，为开发高性能的热电材料提供理论指导。它还可以与实验研究相结合，通过将计算结果与实验数据进行对比，验证理论模型的正确性，进一步完善对热电材料的理解。三、CuMO₂体系电子结构和热电性质计算3.1CuCrO₂体系计算结果3.1.1电子结构计算结果利用基于密度泛函理论的第一性原理计算方法，通过Wien2k软件对CuCrO₂体系的电子结构进行了深入计算。计算过程中，采用广义梯度近似（GGA）来处理交换关联能，以更准确地描述电子之间的相互作用。在结构优化方面，对晶格常数和原子坐标进行了充分优化，确保体系处于能量最低的稳定状态。图1展示了CuCrO₂的能带结构。从能带结构中可以清晰地看出，CuCrO₂是一种直接带隙半导体，其带隙值为[具体带隙数值]eV。在价带部分，主要由Cu的3d轨道电子和O的2p轨道电子形成杂化轨道。这种杂化作用使得价带具有一定的宽度，增强了电子在价带中的离域性。在导带中，Cr的3d轨道电子对导带的贡献较大。费米能级位于价带和导带之间，且远离导带底和价带顶，这表明在常温下，热激发产生的载流子浓度较低。态密度分析对于理解CuCrO₂体系的电子结构同样至关重要。图2给出了CuCrO₂的总态密度（TDOS）和分波态密度（PDOS）。在费米能级附近，总态密度较低，这与能带结构中费米能级远离导带底和价带顶的特征相一致。从分波态密度来看，在价带顶附近，O的2p态密度和Cu的3d态密度有明显的重叠，这进一步证实了价带中Cu的3d轨道与O的2p轨道存在较强的杂化作用。在导带底附近，Cr的3d态密度占据主导地位，说明导带中的电子主要来源于Cr的3d轨道电子。这种电子结构特点决定了CuCrO₂的电学性质，如较低的本征载流子浓度和一定的电子迁移率。通过对CuCrO₂体系电子结构的计算和分析，明确了其半导体特性以及不同原子轨道对电子态的贡献，为后续研究其热电性质奠定了坚实的基础。例如，带隙的大小和电子态密度的分布将直接影响载流子的激发和输运，从而对电导率、Seebeck系数等热电性质产生重要影响。在后续的热电性质计算中，可以基于这些电子结构信息，更准确地理解和预测CuCrO₂的热电性能。3.1.2热电系数计算结果基于第一性原理计算得到的CuCrO₂体系的电子结构，利用BoltzTraP软件，结合玻尔兹曼理论，对其热电系数进行了计算。计算过程中，考虑了电子的散射机制，包括声子散射和杂质散射等，以更准确地描述电子的输运过程。在声子散射方面，采用了形变势理论来描述电子-声子相互作用；对于杂质散射，假设杂质均匀分布在体系中，并根据杂质浓度和散射截面来计算杂质散射对电子输运的影响。图3展示了CuCrO₂的Seebeck系数随温度和载流子浓度的变化关系。可以看出，在低温范围内，Seebeck系数随着温度的升高而增大。这是因为在低温下，电子的散射主要以杂质散射为主，随着温度升高，声子散射逐渐增强，电子的散射机制发生变化，导致载流子的能量分布发生改变，从而使得Seebeck系数增大。当温度进一步升高时，Seebeck系数逐渐趋于饱和。这是由于高温下声子散射占主导地位，载流子的散射几率趋于稳定，使得Seebeck系数的变化不再明显。在不同载流子浓度下，Seebeck系数随着载流子浓度的增加而减小。这是因为载流子浓度的增加使得电子的费米能级升高，更多的低能态电子参与导电，导致载流子的平均能量降低，从而Seebeck系数减小。电导率是描述材料导电能力的重要参数。图4给出了CuCrO₂的电导率随温度和载流子浓度的变化曲线。在低温时，电导率随着温度的升高而增大。这是因为低温下杂质散射较强，限制了电子的输运，随着温度升高，电子的热运动加剧，部分电子能够克服杂质散射的阻碍，从而使电导率增大。在高温阶段，电导率随着温度的升高而减小。这是由于高温下声子散射增强，电子与声子的碰撞几率增加，导致电子的平均自由程减小，从而电导率降低。载流子浓度对电导率有显著影响，随着载流子浓度的增加，电导率显著增大。这是因为载流子浓度的增加使得参与导电的载流子数量增多，在相同电场作用下，能够形成更大的电流密度，从而提高电导率。通过对CuCrO₂体系热电系数的计算和分析，明确了Seebeck系数和电导率随温度和载流子浓度的变化规律。这些规律与电子结构密切相关。例如，电子态密度在费米能级附近的分布决定了载流子的浓度和能量分布，进而影响Seebeck系数和电导率。声子散射和杂质散射等散射机制对电子输运的影响也在热电系数的变化中得到体现。这些结果为进一步研究CuCrO₂的热电性能和优化其热电性质提供了重要依据。后续可以通过调控电子结构和散射机制，如通过掺杂等手段改变载流子浓度和散射几率，来优化CuCrO₂的热电性能，提高其热电转换效率。3.2CuRhO₂体系计算结果3.2.1电子结构计算结果运用第一性原理计算方法，基于Wien2k软件对CuRhO₂体系的电子结构进行深入剖析。计算过程中，选用广义梯度近似（GGA）来处理交换关联能，从而精准描述电子间的相互作用。在结构优化环节，细致地对晶格常数和原子坐标进行优化，以确保体系处于能量最低的稳定状态。图5展示了CuRhO₂的能带结构。经分析可知，CuRhO₂同样属于直接带隙半导体，其带隙值为[具体带隙数值]eV，与CuCrO₂的带隙值存在一定差异。在价带部分，主要由Cu的3d轨道电子和O的2p轨道电子形成杂化轨道，这与CuCrO₂类似，但杂化程度和轨道分布存在细微差别。在导带中，Rh的4d轨道电子对导带的贡献较大。费米能级位于价带和导带之间，且与导带底和价带顶的距离与CuCrO₂有所不同，这表明其载流子激发情况和浓度与CuCrO₂存在差异。态密度分析对于理解CuRhO₂体系的电子结构至关重要。图6给出了CuRhO₂的总态密度（TDOS）和分波态密度（PDOS）。在费米能级附近，总态密度的数值与CuCrO₂有所不同。从分波态密度来看，在价带顶附近，O的2p态密度和Cu的3d态密度的重叠情况与CuCrO₂存在差异，这反映了它们杂化作用的不同。在导带底附近，Rh的4d态密度占据主导地位，与CuCrO₂中Cr的3d态密度的主导情况类似，但由于Rh和Cr的电子结构不同，其态密度分布和电子贡献也有所不同。这种电子结构特点决定了CuRhO₂的电学性质，如载流子浓度和迁移率与CuCrO₂会有所区别。与CuCrO₂体系相比，CuRhO₂体系在电子结构上既有相似之处，也有明显差异。相似点在于二者均为直接带隙半导体，价带都主要由Cu的3d轨道电子和O的2p轨道电子杂化形成，导带都主要由过渡金属的d轨道电子贡献。不同点体现在带隙值的差异，以及杂化轨道的具体分布和态密度在费米能级附近的不同。这些差异将导致它们在热电性质上表现出不同的特性。3.2.2热电系数计算结果基于第一性原理计算得到的CuRhO₂体系的电子结构，借助BoltzTraP软件，结合玻尔兹曼理论，对其热电系数展开计算。计算过程中，充分考虑电子的散射机制，涵盖声子散射和杂质散射等，以此更精确地描述电子的输运过程。在声子散射方面，采用形变势理论来刻画电子-声子相互作用；对于杂质散射，假定杂质均匀分布在体系中，并依据杂质浓度和散射截面来计算杂质散射对电子输运的影响。图7展示了CuRhO₂的Seebeck系数随温度和载流子浓度的变化关系。在低温区间，Seebeck系数随温度升高而增大，这与CuCrO₂的变化趋势一致，原因是低温下电子散射以杂质散射为主，随着温度升高，声子散射逐渐增强，改变了载流子的能量分布，致使Seebeck系数增大。当温度进一步升高时，Seebeck系数逐渐趋于饱和，这同样与CuCrO₂类似，是因为高温下声子散射占主导，载流子散射几率趋于稳定。在不同载流子浓度下，Seebeck系数随载流子浓度的增加而减小，与CuCrO₂的规律相同，是由于载流子浓度增加使电子费米能级升高，更多低能态电子参与导电，降低了载流子的平均能量，进而导致Seebeck系数减小。然而，在相同温度和载流子浓度下，CuRhO₂的Seebeck系数数值与CuCrO₂存在差异，这源于二者电子结构的不同。电导率是衡量材料导电能力的关键参数。图8给出了CuRhO₂的电导率随温度和载流子浓度的变化曲线。在低温时，电导率随温度升高而增大，这与CuCrO₂的变化趋势一致，是因为低温下杂质散射较强，限制了电子输运，随着温度升高，电子热运动加剧，部分电子能够克服杂质散射的阻碍，从而使电导率增大。在高温阶段，电导率随温度升高而减小，这也与CuCrO₂类似，是由于高温下声子散射增强，电子与声子碰撞几率增加，导致电子平均自由程减小，进而电导率降低。载流子浓度对电导率有显著影响，随着载流子浓度的增加，电导率显著增大，这与CuCrO₂相同，是因为载流子浓度增加使参与导电的载流子数量增多，在相同电场作用下，能够形成更大的电流密度，从而提高电导率。同样，在相同条件下，CuRhO₂的电导率数值与CuCrO₂有所不同，这是由它们的电子结构差异所导致的。通过对CuRhO₂体系热电系数的计算和分析，明确了其Seebeck系数和电导率随温度和载流子浓度的变化规律。这些规律与CuCrO₂有相似之处，但在具体数值上存在差异，这与它们的电子结构密切相关。例如，电子态密度在费米能级附近的分布决定了载流子的浓度和能量分布，进而影响Seebeck系数和电导率。声子散射和杂质散射等散射机制对电子输运的影响也在热电系数的变化中得以体现。这些结果为进一步研究CuRhO₂的热电性能和优化其热电性质提供了关键依据。后续可通过调控电子结构和散射机制，如采用掺杂等手段改变载流子浓度和散射几率，来优化CuRhO₂的热电性能，提升其热电转换效率。3.3CuAlO₂体系计算结果3.3.1电子结构计算结果采用基于密度泛函理论的第一性原理计算方法，借助Wien2k软件对CuAlO₂体系的电子结构展开深入研究。在计算过程中，选用广义梯度近似（GGA）来处理交换关联能，通过对晶格常数和原子坐标的细致优化，确保体系处于能量最低的稳定状态。图9展示了CuAlO₂的能带结构。分析结果表明，CuAlO₂属于直接带隙半导体，其带隙值为[具体带隙数值]eV。与CuCrO₂和CuRhO₂的带隙值相比，存在明显差异，这主要源于Al元素与Cr、Rh元素电子结构的不同。在价带部分，主要由Cu的3d轨道电子和O的2p轨道电子形成杂化轨道，但杂化程度和轨道分布与CuCrO₂和CuRhO₂有所不同。在导带中，Al的3s和3p轨道电子对导带的贡献较大。费米能级位于价带和导带之间，且与导带底和价带顶的距离独特，这表明其载流子激发情况和浓度与CuCrO₂和CuRhO₂存在区别。态密度分析对于理解CuAlO₂体系的电子结构意义重大。图10给出了CuAlO₂的总态密度（TDOS）和分波态密度（PDOS）。在费米能级附近，总态密度的数值和分布与CuCrO₂和CuRhO₂存在明显差异。从分波态密度来看，在价带顶附近，O的2p态密度和Cu的3d态密度的重叠情况与CuCrO₂和CuRhO₂不同，反映了它们杂化作用的差异。在导带底附近，Al的3s和3p态密度占据主导地位，与CuCrO₂中Cr的3d态密度和CuRhO₂中Rh的4d态密度的主导情况不同，由于Al的电子结构特性，其态密度分布和电子贡献具有独特性。这种电子结构特点决定了CuAlO₂的电学性质，如载流子浓度和迁移率与CuCrO₂和CuRhO₂会有明显区别。与CuCrO₂和CuRhO₂体系相比，CuAlO₂体系在电子结构上具有独特性。虽然它们都为直接带隙半导体，价带都主要由Cu的3d轨道电子和O的2p轨道电子杂化形成，但带隙值、杂化轨道的具体分布和态密度在费米能级附近的情况存在显著差异。这些差异将导致它们在热电性质上表现出不同的特性，为深入研究CuAlO₂体系的热电性能提供了重要的电子结构基础。3.3.2热电系数计算结果基于第一性原理计算得到的CuAlO₂体系的电子结构，利用BoltzTraP软件，结合玻尔兹曼理论，对其热电系数进行了精确计算。计算过程中，全面考虑电子的散射机制，包括声子散射和杂质散射等，以准确描述电子的输运过程。在声子散射方面，采用形变势理论来描述电子-声子相互作用；对于杂质散射，假设杂质均匀分布在体系中，并根据杂质浓度和散射截面来计算杂质散射对电子输运的影响。图11展示了CuAlO₂的Seebeck系数随温度和载流子浓度的变化关系。在低温区间，Seebeck系数随温度升高而增大，这与CuCrO₂和CuRhO₂的变化趋势一致，原因是低温下电子散射以杂质散射为主，随着温度升高，声子散射逐渐增强，改变了载流子的能量分布，致使Seebeck系数增大。当温度进一步升高时，Seebeck系数逐渐趋于饱和，同样与CuCrO₂和CuRhO₂类似，是因为高温下声子散射占主导，载流子散射几率趋于稳定。在不同载流子浓度下，Seebeck系数随载流子浓度的增加而减小，与CuCrO₂和CuRhO₂的规律相同，是由于载流子浓度增加使电子费米能级升高，更多低能态电子参与导电，降低了载流子的平均能量，进而导致Seebeck系数减小。然而，在相同温度和载流子浓度下，CuAlO₂的Seebeck系数数值与CuCrO₂和CuRhO₂存在差异，这是由其电子结构的独特性所导致的。电导率是衡量材料导电能力的关键参数。图12给出了CuAlO₂的电导率随温度和载流子浓度的变化曲线。在低温时，电导率随温度升高而增大，这与CuCrO₂和CuRhO₂的变化趋势一致，是因为低温下杂质散射较强，限制了电子输运，随着温度升高，电子热运动加剧，部分电子能够克服杂质散射的阻碍，从而使电导率增大。在高温阶段，电导率随温度升高而减小，这也与CuCrO₂和CuRhO₂类似，是由于高温下声子散射增强，电子与声子碰撞几率增加，导致电子平均自由程减小，进而电导率降低。载流子浓度对电导率有显著影响，随着载流子浓度的增加，电导率显著增大，这与CuCrO₂和CuRhO₂相同，是因为载流子浓度增加使参与导电的载流子数量增多，在相同电场作用下，能够形成更大的电流密度，从而提高电导率。同样，在相同条件下，CuAlO₂的电导率数值与CuCrO₂和CuRhO₂有所不同，这是由它们的电子结构差异所决定的。通过对CuAlO₂体系热电系数的计算和分析，明确了其Seebeck系数和电导率随温度和载流子浓度的变化规律。这些规律与CuCrO₂和CuRhO₂有相似之处，但在具体数值上存在差异，这与它们的电子结构密切相关。例如，电子态密度在费米能级附近的分布决定了载流子的浓度和能量分布，进而影响Seebeck系数和电导率。声子散射和杂质散射等散射机制对电子输运的影响也在热电系数的变化中得以体现。这些结果为进一步研究CuAlO₂的热电性能和优化其热电性质提供了关键依据。后续可通过调控电子结构和散射机制，如采用掺杂等手段改变载流子浓度和散射几率，来优化CuAlO₂的热电性能，提升其热电转换效率。四、结果分析与讨论4.1电子结构与热电性质关系电子结构是决定材料热电性质的关键因素，其与热电性质之间存在着紧密而复杂的内在联系。对于CuMO₂体系，深入剖析这种联系对于理解其热电性能的物理机制以及进一步优化热电材料具有至关重要的意义。能带结构是电子结构的重要组成部分，它对热电性质有着显著影响。在CuMO₂体系中，带隙大小是一个关键参数。如前文所述，不同的CuMO₂化合物具有不同的带隙值。带隙大小直接影响载流子的激发。较小的带隙使得电子更容易从价带激发到导带，从而增加载流子浓度。当温度升高时，更多的电子能够跨越带隙进入导带，参与导电。在一定温度范围内，随着带隙减小，电导率会相应增加。如果带隙过小，会导致本征载流子浓度过高，可能会降低Seebeck系数。因为Seebeck系数与载流子的能量分布密切相关，过高的载流子浓度会使得载流子的能量分布更加均匀，降低了Seebeck效应。对于CuCrO₂，其带隙值适中，在合适的温度下，能够在保证一定电导率的同时，维持较高的Seebeck系数。而如果带隙过大，如某些理论计算中预测的大带隙CuMO₂体系，载流子激发困难，导致载流子浓度极低，电导率也会很低，虽然Seebeck系数可能较高，但由于电导率过低，功率因子（功率因子=电导率×Seebeck系数²）也会很低，不利于热电转换。能带的形状和色散关系也对热电性质产生重要影响。能带的色散关系决定了电子的有效质量。在能带平坦的区域，电子的有效质量较大，运动速度较慢，迁移率较低；而在能带陡峭的区域，电子的有效质量较小，迁移率较高。迁移率是影响电导率的重要因素之一，电导率与载流子浓度和迁移率成正比。在CuRhO₂中，其导带的色散关系与CuCrO₂不同，导致电子的有效质量和迁移率存在差异，进而影响了电导率。能带的各向异性也会导致热电性质的各向异性。由于CuMO₂体系具有准二维层状结构，电子在层内和层间的运动受到不同的势场作用，导致能带结构在不同方向上存在差异。这种能带的各向异性使得电子在不同方向上的输运性质不同，如电导率和Seebeck系数在层内和层间可能会有较大差异。在一些研究中发现，CuMO₂体系在层内方向的电导率通常高于层间方向，这与能带的各向异性密切相关。态密度在解释电子结构与热电性质关系方面同样起着关键作用。费米能级附近的态密度对电导率和Seebeck系数有着直接影响。在CuMO₂体系中，当费米能级位于态密度较高的区域时，载流子浓度相对较高。因为态密度表示单位能量间隔内的电子态数目