基于第一性原理探究Sm-CeO₂(111)表面CO分子催化氧化微观机制

基于第一性原理探究Sm-CeO₂(111)表面CO分子催化氧化微观机制一、引言1.1研究背景与意义在现代社会，随着工业化进程的加速和汽车保有量的不断增加，环境污染和能源问题日益严峻，CO作为一种常见的有害气体，其排放对环境和人类健康造成了严重威胁。在环保领域，CO主要来源于汽车尾气、工业废气排放等。汽车尾气中含有大量未完全燃烧的CO，当这些尾气排放到大气中，会与人体血液中的血红蛋白结合，降低血液的输氧能力，导致人体缺氧，引发头痛、头晕、恶心等症状，严重时甚至会危及生命。在工业生产中，如钢铁冶炼、化工合成等过程，也会产生大量的CO废气。这些废气如果未经处理直接排放，会对周边空气质量造成严重污染，影响生态平衡。因此，有效治理CO污染成为环保领域的重要任务之一。在能源领域，CO的催化氧化也具有重要意义。在燃料电池中，CO的存在会毒化电极催化剂，降低电池的性能和效率。通过催化氧化去除CO，可以提高燃料电池的稳定性和使用寿命，推动清洁能源的发展。在一些化工合成过程中，CO的氧化反应可以为反应提供能量，实现能量的有效利用。CO催化氧化反应作为消除CO污染的关键技术，一直是研究的热点。传统的CO催化氧化技术主要依赖贵金属催化剂，如铂（Pt）、钯（Pd）等。然而，贵金属催化剂存在成本高、资源稀缺、易中毒等缺点，限制了其大规模应用。因此，开发高效、低成本、稳定的非贵金属催化剂成为CO催化氧化领域的研究重点。氧化铈（CeO₂）作为一种重要的稀土氧化物，具有独特的物理化学性质，在催化领域展现出巨大的应用潜力。CeO₂具有萤石型立方晶体结构，其中每个铈离子被八个氧离子紧密环绕，构筑出三维立方的氧化物骨架。这种结构赋予CeO₂优异的氧化还原性能，Ce元素具有+3和+4两种稳定的氧化态，在不同的环境条件下能够很容易地发生Ce³⁺与Ce⁴⁺之间的氧化还原反应，使其在催化、环境保护等领域表现出独特的性能。在汽车尾气净化催化剂中，CeO₂可以利用其氧化还原性能，存储和释放氧，调节尾气中氧气的浓度，从而促进一氧化碳、碳氢化合物和氮氧化物等污染物的氧化还原反应，使其转化为无害的二氧化碳、水和氮气，达到净化尾气的目的。为了进一步提高CeO₂的催化性能，研究人员通过掺杂其他元素对其进行改性。钐（Sm）作为一种稀土元素，其掺杂可以有效改变CeO₂的晶体结构和电子性质，从而提高其催化活性。Sm-CeO₂复合氧化物在CO催化氧化反应中表现出比纯CeO₂更高的催化活性和稳定性。其原因在于Sm的掺杂引入了晶格缺陷，增加了氧空位的浓度，促进了氧的迁移和活化，从而提高了CO的氧化效率。在众多的晶体表面中，Sm-CeO₂(111)表面由于其特定的原子排列和电子结构，对CO分子具有独特的吸附和催化活性。研究Sm-CeO₂(111)表面CO分子催化氧化机理，对于深入理解CO催化氧化过程、开发高效的CO催化剂具有重要的理论意义。通过第一性原理计算，可以从原子和电子层面揭示CO分子在Sm-CeO₂(111)表面的吸附、反应和脱附过程，为催化剂的设计和优化提供理论指导。在实际应用中，基于对Sm-CeO₂(111)表面CO分子催化氧化机理的研究，可以有针对性地调整催化剂的制备工艺和组成，提高催化剂的性能，降低成本，推动CO催化氧化技术在环保和能源领域的广泛应用。1.2CeO₂的性质及相关应用1.2.1CeO₂的结构与性质氧化铈（CeO₂）在材料科学领域中占据着重要地位，其晶体结构属于立方晶系，具体为萤石型结构。在这种结构中，每个铈（Ce）原子被八个氧（O）原子以立方体的顶点位置包围，形成了一个面心立方的点阵结构。这种紧密堆积的结构赋予了CeO₂较高的稳定性和一定的硬度，使其在高温、高压等极端条件下仍能保持结构的完整性。从电子特性角度来看，Ce元素具有独特的电子构型，其4f电子轨道的存在使得CeO₂展现出特殊的物理化学性质。在CeO₂中，Ce主要以+4价态存在，但由于4f电子的可参与性，Ce也可以呈现+3价态。这种变价特性使得CeO₂在氧化还原反应中表现出卓越的性能，能够在不同的环境条件下快速地进行Ce³⁺与Ce⁴⁺之间的相互转化。当CeO₂处于还原性气氛中时，部分Ce⁴⁺会获得电子被还原为Ce³⁺，同时晶格中会产生氧空位，以维持电荷平衡；而在氧化性气氛中，Ce³⁺又会失去电子被氧化为Ce⁴⁺，氧空位则会被填充。这种氧化还原的可逆性使得CeO₂在催化、电池等领域具有广泛的应用前景。在催化领域，CeO₂的氧化还原性质使其能够有效地促进氧化还原反应的进行。在CO的催化氧化反应中，CeO₂可以通过提供活性氧物种，将CO氧化为CO₂，同时自身被还原。随后，在氧气的存在下，CeO₂又能够重新被氧化，恢复到初始状态，从而实现催化循环。在汽车尾气净化催化剂中，CeO₂利用其储氧和释氧能力，调节尾气中氧气的浓度，促进一氧化碳、碳氢化合物和氮氧化物等污染物的氧化还原反应，使其转化为无害的二氧化碳、水和氮气。在能源领域，CeO₂的氧化还原性质也发挥着重要作用。在固体氧化物燃料电池中，CeO₂作为电解质材料，不仅需要具备良好的离子导电性，还需要在电池的工作条件下能够稳定地进行氧化还原反应，以维持电池的正常运行。CeO₂的氧化还原特性使其能够在电池的阳极和阴极之间传递氧离子，促进电极反应的进行，提高电池的性能和效率。1.2.2CeO₂在汽车尾气催化方面的应用及原理随着汽车保有量的持续增长，汽车尾气排放对环境造成的污染日益严重。汽车尾气中含有一氧化碳（CO）、碳氢化合物（HC）和氮氧化物（NOx）等多种有害污染物，这些污染物不仅会对空气质量产生负面影响，还会危害人体健康。CeO₂作为一种重要的汽车尾气催化材料，在尾气净化过程中发挥着关键作用。CeO₂在汽车尾气催化中的主要作用是利用其独特的储氧能力。在汽车发动机的工作过程中，尾气的组成会随着工况的变化而发生波动，时而处于富燃状态（氧气不足），时而处于贫燃状态（氧气过量）。CeO₂能够在贫燃条件下储存氧气，将其晶格中的氧原子与尾气中的氧气结合，形成高价态的CeO₂；而在富燃条件下，CeO₂又能够释放储存的氧气，将其提供给尾气中的CO、HC等还原性污染物，促进它们的氧化反应。这种储氧和释氧的过程可以有效地调节尾气中氧气的浓度，使催化剂始终处于最佳的工作状态，提高尾气净化效率。在尾气催化反应中，CeO₂与其他催化剂组分（如贵金属催化剂）协同作用，共同促进CO、HC和NOx的转化。对于CO的氧化反应，在贵金属催化剂（如Pt、Pd等）的作用下，CO分子被吸附在催化剂表面，同时CeO₂释放出的活性氧物种也吸附在催化剂表面。CO与活性氧发生反应，生成CO₂，从而实现CO的去除。反应方程式如下：2CO+O_2\stackrel{催化剂}{\longrightarrow}2CO_2，在这个过程中，CeO₂不仅提供了活性氧，还通过与贵金属的相互作用，改变了贵金属的电子结构和表面性质，提高了贵金属的催化活性和稳定性。对于HC的氧化反应，CeO₂同样发挥着重要作用。HC分子在催化剂表面被吸附后，与CeO₂释放的活性氧发生反应，逐步被氧化为CO₂和H₂O。不同类型的HC（如烷烃、烯烃、芳烃等）在催化氧化过程中的反应路径和速率有所不同，但CeO₂的存在都能够有效地促进这些反应的进行。在NOx的还原反应中，CeO₂可以通过与还原剂（如CO、HC等）的相互作用，将NOx还原为N₂。在这个过程中，CeO₂参与了反应中间物种的形成和转化，促进了NOx的还原过程。例如，在选择性催化还原（SCR）反应中，CeO₂可以作为助催化剂，提高SCR催化剂的活性和选择性，使NOx在较低的温度下就能被有效地还原。1.2.3CO的催化氧化机理CO的催化氧化是一个复杂的过程，涉及多个反应步骤和表面物种。目前，关于CO催化氧化的机理主要有两种被广泛接受的理论：Langmuir-Hinshelwood（L-H）机理和Eley-Rideal（E-R）机理。在L-H机理中，CO和O₂分子首先分别吸附在催化剂表面的不同活性位点上。CO分子通过C原子与催化剂表面的金属原子或氧原子形成化学键，而O₂分子则通过解离吸附的方式，在催化剂表面形成两个活性氧原子。吸附的CO和活性氧原子之间发生反应，生成CO₂分子，然后CO₂分子从催化剂表面脱附。整个反应过程可以表示为：O_2+2*\longrightarrow2O*（O₂分子在催化剂表面解离吸附，*表示催化剂表面的活性位点）CO+*\longrightarrowCO*（CO分子在催化剂表面吸附）CO*+O*\longrightarrowCO₂*（吸附的CO和活性氧原子反应生成CO₂）CO₂*\longrightarrowCO₂+*（CO₂分子从催化剂表面脱附）在E-R机理中，O₂分子首先在催化剂表面解离吸附，形成活性氧原子。然后，气相中的CO分子直接与催化剂表面吸附的活性氧原子发生反应，生成CO₂分子，并从催化剂表面脱附。反应过程如下：O_2+2*\longrightarrow2O*（O₂分子在催化剂表面解离吸附）CO+O*\longrightarrowCO₂*（气相CO与吸附的活性氧反应生成CO₂）CO₂*\longrightarrowCO₂+*（CO₂分子从催化剂表面脱附）除了这两种基本机理外，实际的CO催化氧化过程可能还涉及其他复杂的因素，如催化剂表面的缺陷、杂质、反应物和产物的扩散等。在一些具有氧空位的催化剂表面，氧空位可以作为活性位点，促进O₂分子的吸附和解离，同时也可以增强CO分子的吸附能力，从而加速CO的催化氧化反应。催化剂表面的杂质可能会影响活性位点的性质和分布，进而影响反应的速率和选择性。反应物和产物在催化剂表面和体相中的扩散速率也会对反应速率产生重要影响，如果扩散速率较慢，可能会导致反应物在催化剂表面的积累或产物的脱附困难，从而降低反应效率。在CeO₂基催化剂上，CO的催化氧化机理还与CeO₂的氧化还原性质密切相关。CeO₂可以通过自身的氧化还原循环，提供活性氧物种，参与CO的氧化反应。在反应过程中，CeO₂中的Ce⁴⁺会被还原为Ce³⁺，同时产生氧空位，这些氧空位可以吸附O₂分子并使其解离，生成的活性氧原子再与CO反应。随着反应的进行，Ce³⁺又会被氧化为Ce⁴⁺，氧空位被填充，完成一个氧化还原循环。这种氧化还原循环机制使得CeO₂在CO催化氧化中具有独特的性能。1.2.4CeO₂在其他方面的应用CeO₂凭借其优异的物理化学性质，在众多领域展现出广泛的应用前景。在燃料电池领域，CeO₂作为固体氧化物燃料电池（SOFC）的关键材料之一，发挥着重要作用。在SOFC中，CeO₂既可以作为电解质材料，又可以作为电极催化剂的组成部分。作为电解质，CeO₂具有较高的氧离子导电性，能够在高温下有效地传导氧离子，促进电池的电极反应。与传统的氧化钇稳定的氧化锆（YSZ）电解质相比，CeO₂基电解质在较低温度下具有更好的离子导电性，这有助于降低SOFC的工作温度，提高电池的效率和稳定性，减少电池组件的材料成本和制造难度。在有机合成领域，CeO₂也被广泛用作催化剂或催化剂载体。在一些氧化反应中，CeO₂可以作为活性催化剂，直接参与反应过程，促进反应物的氧化。在醇类的氧化反应中，CeO₂可以将醇氧化为相应的醛或酮，其催化活性和选择性与CeO₂的晶体结构、表面性质以及反应条件等因素密切相关。CeO₂还可以作为催化剂载体，负载其他活性金属或金属氧化物，提高催化剂的分散性和稳定性，从而增强催化剂的性能。负载在CeO₂上的贵金属催化剂在一些有机合成反应中表现出更高的活性和选择性，能够有效地促进反应的进行，提高目标产物的收率。在环境保护领域，除了汽车尾气净化外，CeO₂还可用于处理工业废气中的其他污染物，如挥发性有机化合物（VOCs）。CeO₂基催化剂能够有效地催化氧化VOCs，将其转化为无害的CO₂和H₂O，从而减少VOCs对大气环境的污染。在光催化领域，CeO₂也具有一定的应用潜力。通过与其他半导体材料复合，CeO₂可以提高光催化反应的效率，用于降解有机污染物、分解水制氢等。在电子材料领域，CeO₂可用于制备电阻器、电容器等电子元件。由于其具有一定的导电性和介电性能，CeO₂可以作为电子材料的组成部分，调节材料的电学性能，满足不同电子元件的需求。在抛光材料领域，CeO₂因其硬度适中、化学稳定性好等特点，被广泛应用于光学镜片、半导体芯片等精密器件的抛光过程，能够有效地去除材料表面的划痕和杂质，提高表面的平整度和光洁度。1.3国内外研究现状及发展趋势1.3.1负载或掺杂金属原子提高CeO₂催化活性为了提升CeO₂的催化活性，国内外学者广泛开展了负载或掺杂金属原子的研究。在负载金属原子方面，贵金属如Pt、Pd、Au等常被负载于CeO₂表面，以增强其催化性能。研究表明，负载Pt的CeO₂催化剂在CO氧化反应中展现出优异的低温活性，Pt原子的存在促进了CO分子在催化剂表面的吸附与活化，降低了反应的活化能。这是因为Pt具有良好的电子传导能力，能够快速地将电子传递给CO分子，使其更容易被氧化。对于非贵金属负载，过渡金属如Mn、Fe、Co等也受到了极大关注。负载Mn的CeO₂催化剂在挥发性有机化合物（VOCs）催化氧化中表现出良好的活性，Mn的负载改变了CeO₂的表面氧物种分布，增加了表面活性氧的含量，从而提高了催化剂的氧化能力。Mn的电子结构与CeO₂相互作用，使得表面氧物种的活性增强，更易于与VOCs分子发生反应。在掺杂金属原子方面，稀土元素如La、Sm、Gd等常被用于掺杂CeO₂。La掺杂的CeO₂在汽车尾气净化中具有重要应用，La的掺杂稳定了CeO₂的晶体结构，提高了其储氧能力和热稳定性，进而提升了催化剂在复杂工况下的性能。La的离子半径与Ce相近，掺杂后能够较好地融入CeO₂晶格，减少晶格缺陷的产生，增强晶体结构的稳定性。过渡金属离子掺杂同样能够显著影响CeO₂的催化性能。Zr掺杂的CeO₂在固体氧化物燃料电池（SOFC）中作为电解质材料，Zr的掺杂提高了CeO₂的氧离子导电性，降低了电池的工作温度，提高了电池的效率和稳定性。Zr的掺杂改变了CeO₂晶格的氧空位浓度和分布，使得氧离子在晶格中的迁移更加容易，从而提高了氧离子导电性。不同金属原子的负载或掺杂对CeO₂催化活性和选择性的影响机制各有不同。负载金属原子主要通过改变催化剂表面的活性位点和电子结构，影响反应物的吸附和活化；而掺杂金属原子则通过改变CeO₂的晶体结构、电子性质和氧空位浓度等，调控催化剂的氧化还原性能和表面活性。负载Au的CeO₂催化剂中，Au原子作为活性位点，能够特异性地吸附CO分子，使其在较低温度下就能够被氧化；而Sm掺杂的CeO₂中，Sm的存在增加了氧空位的浓度，促进了氧的迁移和活化，提高了CO氧化的效率。1.3.2催化反应过程机理的研究催化反应过程机理的研究一直是催化领域的核心内容。实验研究是探究催化反应机理的重要手段之一。通过原位光谱技术，如原位红外光谱（in-situIR）、原位拉曼光谱（in-situRaman）等，能够实时监测催化反应过程中反应物、中间体和产物在催化剂表面的吸附、反应和脱附过程。在CO催化氧化反应中，原位红外光谱可以检测到CO分子在催化剂表面的吸附态，以及反应过程中生成的CO₂的特征吸收峰，从而推断反应的中间步骤和反应路径。程序升温技术，如程序升温脱附（TPD）、程序升温还原（TPR）等，能够研究反应物在催化剂表面的吸附强度、脱附行为以及催化剂的氧化还原性能。TPD可以测定CO在催化剂表面的脱附温度，从而了解CO与催化剂表面的相互作用强弱；TPR则可以表征催化剂被还原的难易程度，反映其氧化还原活性。理论计算方法，如第一性原理计算、分子动力学模拟等，为深入理解催化反应机理提供了微观层面的信息。第一性原理计算基于量子力学理论，能够精确计算催化剂的电子结构、吸附能、反应能垒等参数，从而揭示催化反应的本质。通过第一性原理计算，可以研究CO分子在催化剂表面的吸附构型和吸附能，分析反应过程中化学键的断裂和形成，确定反应的决速步骤和反应机理。分子动力学模拟则可以在原子尺度上模拟催化反应过程，研究反应物和产物在催化剂表面的扩散行为、反应动力学等。在研究CO催化氧化反应时，分子动力学模拟可以模拟CO和O₂分子在催化剂表面的扩散过程，以及它们之间的反应动力学，为优化催化剂的结构和性能提供理论依据。实验和理论计算相结合的方法，能够更全面、深入地研究催化反应机理。通过实验观察到的现象，利用理论计算进行微观层面的解释和验证；而理论计算预测的结果，又可以通过实验进行验证和优化。这种相互验证和补充的研究方法，推动了催化反应机理研究的不断深入。1.3.3国内外研究CeO₂趋势目前，国内外对CeO₂的研究呈现出多方面的发展趋势。在新型催化剂设计方面，研究者致力于开发具有更高活性、选择性和稳定性的CeO₂基催化剂。通过调控催化剂的微观结构，如纳米结构、多孔结构等，增加活性位点的数量和暴露程度，提高催化剂的性能。设计具有纳米多孔结构的CeO₂催化剂，能够提供更大的比表面积和更多的活性位点，促进反应物的吸附和反应。探索新型的掺杂或负载策略，引入更多种类的金属或非金属元素，实现催化剂性能的协同优化。研究发现，同时掺杂两种或多种金属元素的CeO₂催化剂，在某些反应中表现出比单一掺杂催化剂更优异的性能，这是由于多种元素之间的协同作用，优化了催化剂的电子结构和表面性质。在多尺度研究方面，结合宏观实验、微观理论计算和介观模拟，从不同尺度全面理解CeO₂的催化性能和反应机理。通过宏观实验获得催化剂的整体性能数据，利用微观理论计算揭示原子和电子层面的反应机制，借助介观模拟研究催化剂的结构与性能之间的关系，实现对CeO₂催化过程的深入认识。利用宏观实验测定CeO₂催化剂在CO氧化反应中的转化率和选择性，通过微观理论计算分析CO分子在催化剂表面的吸附和反应过程，再通过介观模拟研究催化剂颗粒的大小、形状等因素对反应性能的影响。在工业应用方面，CeO₂催化剂的大规模制备技术和应用工艺的优化成为研究重点。开发高效、低成本的制备方法，提高催化剂的生产效率和质量稳定性，以满足工业生产的需求。探索CeO₂催化剂在更多工业领域的应用，如化工合成、废气处理等，拓展其应用范围。在化工合成中，CeO₂催化剂可以用于一些重要的有机合成反应，提高反应的选择性和产率；在废气处理中，CeO₂催化剂可以用于处理工业废气中的多种污染物，实现废气的达标排放和资源的回收利用。1.4课题提出和研究内容1.4.1课题提出及重要性在当前环保和能源领域的研究中，CO的催化氧化一直是备受关注的焦点。CO作为一种有害气体，其排放不仅对环境造成严重污染，如形成酸雨、破坏臭氧层等，还对人类健康构成直接威胁，引发呼吸系统疾病、心血管疾病等。传统的CO催化氧化技术依赖贵金属催化剂，然而，这些贵金属资源稀缺、成本高昂，且易受到杂质的影响而中毒失活，这在很大程度上限制了其大规模的实际应用。因此，开发高效、低成本且稳定的非贵金属催化剂成为解决CO污染问题的关键。氧化铈（CeO₂）因其独特的萤石型立方晶体结构和优异的氧化还原性能，在催化领域展现出巨大的潜力。CeO₂中的Ce元素具有+3和+4两种稳定的氧化态，能够在不同的环境条件下快速进行氧化还原反应，这使得CeO₂在CO催化氧化中具有重要的应用价值。为了进一步提升CeO₂的催化性能，通过掺杂其他元素进行改性成为一种有效的策略。钐（Sm）作为一种稀土元素，其掺杂可以显著改变CeO₂的晶体结构和电子性质，从而提高其催化活性。研究表明，Sm的掺杂能够引入晶格缺陷，增加氧空位的浓度，促进氧的迁移和活化，进而提高CO的氧化效率。在CeO₂的众多晶体表面中，Sm-CeO₂(111)表面由于其特定的原子排列和电子结构，对CO分子具有独特的吸附和催化活性。研究Sm-CeO₂(111)表面CO分子催化氧化机理，对于深入理解CO催化氧化过程具有重要的理论意义。通过探究这一表面的催化氧化机理，可以从原子和电子层面揭示CO分子在催化剂表面的吸附、反应和脱附过程，为CO催化氧化反应提供微观层面的理论依据，填补该领域在微观机理研究方面的空白。从实际应用角度来看，基于对Sm-CeO₂(111)表面CO分子催化氧化机理的研究，可以有针对性地设计和优化催化剂的结构和组成。通过调整Sm的掺杂量、制备工艺等参数，可以提高催化剂的活性、选择性和稳定性，降低催化剂的成本，从而推动CO催化氧化技术在环保和能源领域的广泛应用。在汽车尾气净化中，可以开发出更高效的催化剂，降低尾气中CO的排放，减少对环境的污染；在燃料电池中，可以提高电池的性能和稳定性，促进清洁能源的发展。因此，本研究对于解决环境污染和能源问题具有重要的现实意义。1.4.2研究主要内容本研究运用第一性原理计算方法，对Sm-CeO₂(111)表面CO分子催化氧化机理展开深入研究，具体内容如下：Sm-CeO₂(111)表面结构的优化与分析：采用第一性原理计算，对Sm-CeO₂(111)表面的原子结构进行优化。通过计算不同原子位置的能量，确定最稳定的表面结构，分析表面原子的排列方式、键长、键角等结构参数。研究Sm原子在CeO₂(111)表面的掺杂位置对表面结构的影响，以及由此产生的晶格畸变和应力分布。通过对表面结构的分析，为后续的CO吸附和氧化反应研究提供基础。CO在Sm-CeO₂(111)表面的吸附特性研究：计算CO分子在Sm-CeO₂(111)表面不同吸附位点的吸附能，确定最稳定的吸附构型。分析CO分子与表面原子之间的相互作用，包括化学键的形成、电子转移等。研究吸附过程中CO分子的电子结构变化，以及表面电子云分布的改变。通过吸附特性的研究，了解CO分子在Sm-CeO₂(111)表面的吸附行为，为氧化反应的起始阶段提供理论依据。CO在Sm-CeO₂(111)表面的氧化过程研究：构建CO在Sm-CeO₂(111)表面氧化的反应路径，计算反应过程中各中间体和过渡态的能量。确定反应的决速步骤，分析反应的活化能和反应热。研究氧化过程中表面氧物种的参与方式和作用，以及CeO₂的氧化还原循环机制。通过对氧化过程的研究，揭示CO在Sm-CeO₂(111)表面催化氧化的微观机理。反应条件对CO催化氧化性能的影响研究：考虑温度、压力等反应条件对CO在Sm-CeO₂(111)表面催化氧化性能的影响。通过热力学和动力学分析，计算不同反应条件下的反应速率常数和平衡常数。研究反应条件对反应路径和产物选择性的影响，为实际应用中优化反应条件提供理论指导。二、第一性原理的理论基础2.1绝热近似在研究多粒子体系的过程中，由于原子核与电子的质量存在显著差异，绝热近似应运而生。原子核的质量相较于电子质量大得多，约为电子质量的1836倍。这种巨大的质量差异导致电子的运动速度远快于原子核，使得在电子运动的时间尺度内，原子核的位置变化可以忽略不计。形象地说，电子如同高速运动的粒子，而原子核则像是静止不动的物体。基于这一特性，绝热近似将原子核和电子的运动进行分离处理。在处理电子的运动时，可假定原子核固定在其瞬时位置上，不考虑原子核的动能，仅关注电子在原子核所产生的固定势场中的运动。此时，多粒子体系的哈密顿量可分为电子哈密顿量和原子核哈密顿量两部分。电子哈密顿量描述了电子的动能、电子之间的相互作用以及电子与原子核之间的相互作用；而原子核哈密顿量则主要涉及原子核的动能以及原子核之间的相互作用。通过这种分离，原本复杂的多粒子体系问题被简化为相对简单的多电子问题。在实际计算中，绝热近似极大地降低了计算的复杂度。在研究晶体的电子结构时，无需同时考虑原子核和电子的动态变化，只需专注于电子在固定原子核势场中的行为。这使得计算量大幅减少，计算效率显著提高，为后续运用其他理论和方法进行深入计算奠定了基础。从量子力学的角度来看，绝热近似是合理且有效的。电子的波函数主要由其与原子核之间的相互作用决定，而在短时间内，原子核的缓慢运动对电子波函数的影响较小。因此，将原子核视为固定不动，能够准确地描述电子的基态和低激发态性质。在研究分子的电子结构和化学反应过程中，绝热近似能够提供较为准确的结果，帮助研究者理解分子的成键、反应机理等重要信息。2.2单电子近似在多电子体系中，电子之间存在着复杂的相互作用，包括电子与电子之间的库仑排斥力、交换作用以及关联作用等，这使得直接求解多电子体系的薛定谔方程变得极为困难。单电子近似作为一种重要的理论方法，旨在将多电子体系简化为单电子问题，从而使计算得以实现。单电子近似的核心思想是，将每个电子视为在固定的原子核势场以及其他电子的平均场中独立运动。在晶体中，原子核构成了周期性的晶格结构，电子在这个晶格势场中运动。同时，由于电子之间的相互作用较为复杂，难以精确处理，单电子近似采用平均场的概念，将其他电子对某个电子的作用平均化，视为一个等效的势场。形象地说，就好像每个电子在一个由原子核和其他电子的平均影响所构成的“平均势场”中独自运动，而忽略了电子之间的瞬时相互作用。以金属晶体为例，在金属中，价电子在整个晶体中自由移动。单电子近似认为，每个价电子所感受到的势场，是由金属离子形成的周期性晶格势场，以及其他价电子的平均分布所产生的势场之和。在这种近似下，每个价电子的运动可以用一个单电子薛定谔方程来描述，方程中的势能项包含了原子核的吸引势以及其他电子的平均排斥势。通过单电子近似，多电子体系的薛定谔方程被简化为一组单电子的薛定谔方程。每个单电子的波函数描述了该电子在平均势场中的运动状态，而体系的总波函数可以近似表示为这些单电子波函数的乘积。这种简化极大地降低了计算的复杂度，使得我们能够利用已有的数学方法和计算技术来求解电子的能量和波函数。在实际计算中，我们可以通过迭代的方法，不断调整平均势场，使得计算得到的电子波函数和能量能够自洽，即满足一定的收敛条件。单电子近似在许多领域都有着广泛的应用。在材料科学中，它是计算材料电子结构的基础，通过单电子近似可以计算材料的能带结构、态密度等重要性质，从而深入理解材料的电学、光学、磁性等物理性质。在化学中，单电子近似可以用于研究分子的电子结构和化学反应机理，帮助我们理解分子的成键方式、反应活性等。2.3密度泛函理论2.3.1Hohenberg-Kohn定理密度泛函理论（DFT）的基石是Hohenberg-Kohn定理，该定理为用电子密度来描述多电子体系的基态性质提供了坚实的理论依据。Hohenberg-Kohn第一定理表明，对于一个处在外部势场V_{ext}(r)中的多电子体系，其基态电子密度\rho(r)是唯一确定的，并且反过来，外部势场V_{ext}(r)（除了一个常数项）也由基态电子密度\rho(r)唯一确定。这意味着体系的所有性质都可以由电子密度\rho(r)来确定，因为外部势场决定了体系的哈密顿量，进而决定了体系的所有量子力学性质。从物理意义上讲，电子密度是电子在空间中的分布情况，它反映了电子与原子核以及电子之间的相互作用，因此通过电子密度可以确定体系的能量、电荷分布等重要性质。Hohenberg-Kohn第二定理指出，体系的基态能量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]。这一定理从数学上证明了可以通过寻找合适的电子密度泛函来计算多电子体系的基态能量和其他性质，而无需像传统的量子力学方法那样直接处理多电子波函数。由于电子密度是一个三维空间的函数，相比多电子波函数（其变量个数随电子数增加而急剧增加），处理起来要简单得多，这为发展高效的计算方法提供了可能，使得密度泛函理论在处理大量原子的复杂体系时具有很大的优势，在物理、化学、材料科学等领域得到了广泛的应用。例如，在研究材料的电子结构时，可以通过计算电子密度的分布来确定材料的能带结构、态密度等性质，从而深入理解材料的电学、光学、磁性等物理性质。2.3.2Ko-Sham方程基于Hohenberg-Kohn定理，Kohn和Sham提出了著名的Ko-Sham方程，为计算多电子体系的基态性质提供了一种有效的方法。Ko-Sham方程的核心思想是将复杂的多体问题转化为单电子有效势问题。具体来说，引入一个无相互作用的参考体系，其电子密度与真实多电子体系的电子密度相同。通过变分原理，对能量泛函求极值，得到Ko-Sham方程，其具体形式为：[-\frac{\hbar^2}{2m}\nabla^2+V_{eff}(r)]\psi_i(r)=\epsilon_i\psi_i(r)，其中\psi_i(r)是第i个Ko-Sham轨道波函数，\epsilon_i是对应的轨道能量，V_{eff}(r)是有效势能。有效势能V_{eff}(r)包括几个重要部分：外部势能V_{ext}(r)，例如原子核产生的势能，它描述了电子与原子核之间的静电吸引作用；电子间的库仑势能V_H(r)，即Hartree势能，描述电子间的直接库仑相互作用，其表达式为V_H(r)=\int\frac{\rho(r')}{|r-r'|}dr'，其中\rho(r')是电子密度；交换关联势能V_{xc}(r)，描述电子间的交换和关联作用，这是电子之间复杂相互作用的体现，由于其精确形式难以确定，通常需要通过近似来计算。Ko-Sham方程将多电子问题转化为一组单粒子方程，使得我们可以通过求解这些单粒子方程来获得多电子体系的基态能量和电子密度等性质。在实际计算中，通过迭代的方法，不断调整电子密度和有效势能，使得计算结果达到自洽，即前后两次计算得到的电子密度和能量的差异在一定的误差范围内。这种迭代过程使得我们能够逐步逼近多电子体系的真实基态性质。通过求解Ko-Sham方程，可以得到材料中电子的波函数和能量，进而计算出材料的各种物理性质，如能带结构、电荷密度分布等。2.3.3局域密度近似(LDA)与广义梯度近似(GGA)在密度泛函理论的实际应用中，交换关联能的计算是一个关键问题，由于其精确形式难以确定，通常需要采用近似方法。局域密度近似（LDA）和广义梯度近似（GGA）是两种常用的近似方法，它们在处理交换关联能时采用了不同的策略，具有各自的特点和适用范围。LDA假设交换关联能只与局域的电子密度有关，即认为体系中某一点的交换关联能等于相同密度的均匀电子气的交换关联能。其交换关联能泛函E_{xc}^{LDA}[\rho]可以表示为E_{xc}^{LDA}[\rho]=\int\rho(r)\epsilon_{xc}(\rho(r))dr，其中\epsilon_{xc}(\rho(r))是均匀电子气的交换关联能密度，它是电子密度\rho(r)的函数。LDA在处理一些简单体系时，如金属和半导体，能够给出较为合理的结果，计算效率较高，因为它只依赖于局域的电子密度信息，计算过程相对简单。但LDA也存在一定的局限性，它往往会高估体系的结合能，低估体系的带隙。这是因为LDA忽略了电子密度的梯度变化对交换关联能的影响，在实际体系中，电子密度并非完全均匀分布，其梯度变化会对交换关联能产生重要影响，而LDA没有考虑这一因素，导致计算结果与实际情况存在偏差。GGA则考虑了电子密度梯度的修正，认为交换关联能不仅与局域电子密度有关，还与电子密度的梯度有关。GGA的交换关联能泛函E_{xc}^{GGA}[\rho]通常表示为E_{xc}^{GGA}[\rho]=\int\rho(r)\epsilon_{xc}(\rho(r),\nabla\rho(r))dr，其中\epsilon_{xc}(\rho(r),\nabla\rho(r))是依赖于电子密度和密度梯度的交换关联能密度。通过考虑电子密度梯度，GGA能够更准确地描述电子之间的相互作用，从而得到更合理的结合能。在一些复杂体系的计算中，GGA的表现优于LDA，能够给出更接近实验值的结果。GGA泛函可以运用拟合方法提高泛函精度，如BLYP泛函，也可以通过物理限制得到较好的泛函形式，如PBE泛函。然而，GGA计算的带隙仍偏小，在一些对带隙要求较高的应用中，可能需要使用杂化泛函（如HSE03、HSE06）来得到与实验更符合的带隙数值。2.3.4Hubbard参数U的引入(DFT+U计算)在处理一些强关联体系时，如过渡金属氧化物、高温超导体等，传统的密度泛函理论（DFT）存在一定的局限性。这些体系中存在着局域的d、f轨道，电子之间的强关联作用使得电子的行为不能简单地用DFT中的平均场近似来描述。在过渡金属氧化物中，过渡金属离子的d电子之间存在着强烈的库仑相互作用和交换作用，导致电子的局域化程度较高，电子之间的关联效应显著。在这种情况下，传统DFT方法往往会低估电子之间的相互作用，无法准确描述体系的电子结构和物理性质。为了解决这一问题，引入了Hubbard参数U，采用DFT+U方法进行计算。DFT+U方法通过在能量泛函中添加一项与Hubbard参数U相关的修正项，来描述电子之间的强关联作用。这一修正项主要考虑了同一原子上不同轨道电子之间的库仑排斥能，使得计算能够更准确地反映强关联体系中电子的行为。Hubbard参数U的引入具有重要作用。它能够有效地修正DFT方法在处理强关联体系时的误差，使得计算结果更接近实验值。在研究过渡金属氧化物的磁性时，DFT+U方法能够准确地预测出材料的磁矩和磁有序状态，而传统DFT方法往往会得到错误的结果。通过调整Hubbard参数U的值，可以优化计算结果，使其与实验数据更好地吻合。Hubbard参数U的取值通常需要通过实验数据或其他高精度的理论计算来确定，以确保计算的准确性。在一些研究中，通过与角分辨光电子能谱（ARPES）等实验技术的结合，来确定合适的Hubbard参数U，从而提高计算结果的可靠性。2.3.5赝势和PAW(ProjectorAugmentedWave)方法在第一性原理计算中，准确描述原子核与电子之间的相互作用是至关重要的，但同时也是计算中的一个难点。赝势和PAW（ProjectorAugmentedWave）方法作为两种重要的技术手段，在简化原子核与电子相互作用计算方面发挥了关键作用。赝势方法的基本思想是将电子波函数分为外层电子和芯电子两部分，并使它们正交，从而得到比全电子波函数更平缓的赝波函数。赝波函数和全电子波函数对应的本征能级相同，在截断半径外，两类波函数重合。模守恒赝势要求两类波函数在截断半径内对应电荷积分相同，以保证在处理化学成键等问题时的准确性；超软赝势则取消了这一约束，得到了更平滑的赝波函数，虽然在一定程度上牺牲了电荷守恒的严格性，但大大提高了计算效率，使得在处理较大体系时成为可能。赝势方法通过这种方式，有效地简化了原子核与电子之间复杂的相互作用，减少了计算量，使得大规模的第一性原理计算成为现实。在计算晶体的电子结构时，使用赝势方法可以快速地得到体系的能带结构和态密度等信息，为研究晶体的物理性质提供了便利。PAW方法则是在赝势方法的基础上发展起来的一种更为先进的技术。PAW方法通过相对简单的操作，将赝波函数变换得到全电子波函数，它综合了赝势方法计算效率高和全电子方法精度高的优点，逐渐成为主流的计算方法之一。PAW方法在处理重元素和复杂体系时表现出了独特的优势，能够更准确地描述原子核与电子之间的相互作用，从而得到更精确的计算结果。在研究含有过渡金属元素的材料时，PAW方法能够准确地处理过渡金属原子的电子结构，为深入理解材料的物理性质提供了有力的支持。PAW方法在计算过程中能够更好地保留电子的局域信息，使得计算结果在描述化学键的形成和断裂等过程时更加准确，这对于研究化学反应机理等问题具有重要意义。2.4第一性原理计算的求解过程第一性原理计算的求解过程是一个系统而严谨的过程，它基于量子力学的基本原理，通过一系列的数学方法和近似处理，来获得多电子体系的电子结构和能量等重要信息。首先，需要建立合适的模型体系。对于Sm-CeO₂(111)表面CO分子催化氧化体系，要构建包含Sm、Ce、O等原子的晶体结构模型，并考虑表面的原子排列和CO分子的吸附位置。在构建晶体结构模型时，需依据CeO₂的萤石型立方晶体结构，确定Ce和O原子的位置，再根据掺杂比例引入Sm原子，同时考虑(111)表面的原子层和周期性边界条件，以准确描述表面的性质。对于CO分子的吸附，要考虑不同的吸附位点，如顶位、桥位、穴位等，以确定最稳定的吸附构型。基于绝热近似，将原子核和电子的运动分离，把多粒子体系简化为多电子体系。在此基础上，运用单电子近似，将多电子问题转化为单电子问题，即假设每个电子在固定的原子核势场以及其他电子的平均场中独立运动。通过构造单电子哈密顿量，将多电子体系的哈密顿量简化为单电子哈密顿量的总和，从而可以分别求解每个单电子的薛定谔方程。引入密度泛函理论（DFT），将多电子体系的能量表示为电子密度的泛函。根据Hohenberg-Kohn定理，体系的基态能量和其他性质可以通过电子密度来确定。通过变分原理，对能量泛函求极值，得到Ko-Sham方程。在求解Ko-Sham方程时，需要选择合适的交换关联泛函，如局域密度近似（LDA）或广义梯度近似（GGA），以描述电子间的交换和关联作用。对于Sm-CeO₂体系，由于Ce原子存在4f电子，可能需要考虑Hubbard参数U的引入，采用DFT+U方法进行计算，以更准确地描述电子之间的强关联作用。在实际计算中，通常采用平面波赝势方法（PWPM）。将电子波函数用平面波展开，通过赝势来描述原子核与电子之间的相互作用，以减少计算量。在平面波展开中，需要确定合适的截断能，截断能过小会导致计算精度不足，截断能过大则会增加计算量。通过测试不同的截断能，选择既能保证计算精度又能控制计算量的截断能值。利用赝势方法，可以将复杂的原子核与电子相互作用简化为赝势与电子的相互作用，从而提高计算效率。利用迭代的方法求解Ko-Sham方程，使电子密度和能量达到自洽。在迭代过程中，首先猜测一个初始的电子密度，然后计算有效势能，进而求解Ko-Sham方程得到新的电子波函数和能量。根据新的电子波函数计算新的电子密度，与上一次的电子密度进行比较，若两者的差异在设定的收敛标准内，则认为计算达到自洽，得到的电子结构和能量即为体系的基态性质；若未达到收敛标准，则用新的电子密度重新计算有效势能，继续迭代，直到满足收敛条件为止。在收敛过程中，需要设置合适的收敛标准，如电子密度的变化量、能量的变化量等，以确保计算结果的准确性和可靠性。2.5本课题选用的软件和硬件条件2.5.1VASP软件包简介VASP（ViennaAb-initioSimulationPackage）软件包是一款基于密度泛函理论（DFT）的平面波赝势方法的第一性原理计算软件，在材料科学和凝聚态物理等领域得到了广泛应用。其核心优势在于采用平面波基组展开电子波函数，这种方式具有计算精度高、通用性强等特点，能够精确地描述晶体和表面体系的电子结构。在计算晶体的能带结构时，平面波基组能够准确地反映电子在周期性势场中的运动状态，从而得到精确的能带信息。VASP软件利用赝势方法来描述原子核与电子之间的相互作用，极大地简化了计算过程，提高了计算效率。赝势方法将电子波函数分为外层电子和芯电子两部分，并使它们正交，得到比全电子波函数更平缓的赝波函数，从而减少了计算量。VASP支持多种赝势，如投影缀加波（PAW）赝势和超软赝势（USPP），用户可以根据具体的研究需求选择合适的赝势。在研究含有过渡金属元素的材料时，PAW赝势能够更准确地描述过渡金属原子的电子结构，从而得到更精确的计算结果。在处理交换关联能时，VASP软件提供了多种近似方法，包括局域密度近似（LDA）、广义梯度近似（GGA）以及杂化泛函等。这些方法能够根据不同体系的特点，合理地描述电子间的交换和关联作用，为研究材料的各种性质提供了有力的支持。对于金属体系，LDA方法通常能够给出较为准确的结果；而对于一些具有强关联作用的体系，如过渡金属氧化物，GGA方法或杂化泛函能够更好地描述电子间的相互作用，得到更符合实际情况的结果。VASP软件还具备强大的功能模块，能够进行结构优化、分子动力学模拟、态密度计算等多种计算任务。在结构优化方面，VASP软件通过迭代计算，不断调整原子的位置和晶格参数，使体系的能量达到最小值，从而得到稳定的结构。在分子动力学模拟中，VASP软件能够模拟原子在一定温度和压力下的运动轨迹，研究材料的动力学性质。通过态密度计算，VASP软件可以得到材料中电子态在能量空间的分布情况，为分析材料的电学、光学等性质提供重要依据。2.5.2MaterialsStudio软件简介MaterialsStudio软件是一款功能强大的材料模拟软件，在构建和可视化晶体结构方面具有独特的优势。它提供了直观、便捷的操作界面，使用户能够轻松地构建各种复杂的晶体结构模型。用户可以通过简单的鼠标操作，在软件中定义原子的种类、位置和晶格参数，快速搭建出所需的晶体结构。在构建Sm-CeO₂(111)表面结构时，用户可以利用MaterialsStudio软件的晶格构建工具，准确地设置CeO₂的萤石型立方晶格参数，然后根据掺杂比例引入Sm原子，并确定其在(111)表面的位置，从而构建出精确的Sm-CeO₂(111)表面结构模型。MaterialsStudio软件具备丰富的晶体结构数据库，包含了大量常见材料的晶体结构信息，用户可以直接从数据库中获取所需的结构，在此基础上进行修改和完善，大大节省了构建结构的时间和精力。如果用户需要研究某种已知材料的性质，可以直接从数据库中调用其晶体结构，然后进行后续的模拟和分析。该软件还具有出色的可视化功能，能够以多种方式展示晶体结构，帮助用户更直观地理解结构特征。用户可以通过旋转、缩放、剖切等操作，从不同角度观察晶体结构，清晰地看到原子的排列方式、键长、键角等结构信息。MaterialsStudio软件还可以将晶体结构以三维动画的形式展示出来，动态地呈现原子的运动和相互作用过程，为研究晶体的性质和反应机理提供了更直观的手段。在构建Sm-CeO₂(111)表面CO分子吸附模型时，MaterialsStudio软件能够清晰地展示CO分子在表面的吸附位置和吸附构型，以及CO分子与表面原子之间的相互作用，帮助用户更好地理解吸附过程。MaterialsStudio软件还可以与其他计算软件（如VASP）相结合，将构建好的结构模型直接导入到计算软件中进行后续的计算分析，实现了从结构构建到计算模拟的无缝对接，提高了研究效率。2.5.3课题组的硬件条件课题组拥有高性能的计算机集群，为第一性原理计算提供了强大的硬件支持。计算机集群由多台计算节点组成，每个计算节点配备了高性能的中央处理器（CPU）和图形处理器（GPU）。CPU采用了英特尔至强系列处理器，具有多核心、高主频的特点，能够快速处理复杂的计算任务。每个计算节点配备了32核心的英特尔至强铂金8380处理器，主频可达2.3GHz，能够在短时间内完成大规模的矩阵运算和数值迭代。GPU则选用了英伟达的A100系列显卡，其强大的并行计算能力能够显著加速计算过程，特别是在处理大规模的电子结构计算时，能够大大缩短计算时间。A100显卡拥有10752个CUDA核心，能够同时处理大量的计算任务，相比传统的CPU计算，能够将计算时间缩短数倍甚至数十倍。计算机集群还配备了大容量的内存和高速的存储设备。每个计算节点配备了256GB的DDR4内存，能够满足大规模计算对内存的需求。存储设备采用了高速的固态硬盘（SSD）和并行文件系统，能够快速读写数据，保证计算过程的高效性。并行文件系统的带宽可达数GB/s，能够支持多个计算节点同时读写数据，提高了数据传输的效率。为了保证计算机集群的稳定运行，课题组还配备了专业的散热系统和不间断电源（UPS）。散热系统采用了液冷和风冷相结合的方式，能够有效地降低计算节点的温度，保证CPU和GPU在高负载运行时的稳定性。UPS则能够在市电中断时，为计算机集群提供一定时间的电力支持，确保计算任务的正常进行，避免数据丢失。通过这些硬件设备的支持，课题组能够高效地进行第一性原理计算，为研究Sm-CeO₂(111)表面CO分子催化氧化机理提供了坚实的保障。2.6本章小结本章系统阐述了第一性原理的理论基础，包括绝热近似、单电子近似和密度泛函理论等关键内容。绝热近似基于原子核与电子质量的显著差异，将原子核和电子的运动分离，把多粒子体系简化为多电子体系，大大降低了计算的复杂度。单电子近似则将多电子问题转化为单电子问题，假设每个电子在固定的原子核势场以及其他电子的平均场中独立运动，为后续的计算提供了基础。密度泛函理论是第一性原理计算的核心，其基石Hohenberg-Kohn定理证明了多电子体系的基态性质可由电子密度唯一确定。在此基础上发展的Ko-Sham方程，通过引入无相互作用的参考体系，将多电子问题转化为一组单粒子方程，使得计算多电子体系的基态能量和电子密度等性质成为可能。在实际计算中，常用的交换关联能近似方法有局域密度近似（LDA）和广义梯度近似（GGA），它们分别考虑了电子密度的局域和梯度信息，各有其特点和适用范围。对于强关联体系，引入Hubbard参数U的DFT+U方法能够更准确地描述电子之间的强关联作用。赝势和PAW方法则在简化原子核与电子相互作用计算方面发挥了重要作用，赝势方法通过构建比全电子波函数更平缓的赝波函数，减少了计算量；PAW方法则综合了赝势方法计算效率高和全电子方法精度高的优点，成为主流的计算方法之一。本章还详细介绍了第一性原理计算的求解过程，包括建立模型体系、运用各种近似理论和方法进行计算，以及通过迭代求解Ko-Sham方程使电子密度和能量达到自洽。在实际应用中，选用了VASP软件包和MaterialsStudio软件，VASP软件基于密度泛函理论，采用平面波赝势方法，具备强大的计算功能；MaterialsStudio软件则在构建和可视化晶体结构方面具有独特优势。课题组拥有的高性能计算机集群，为第一性原理计算提供了强大的硬件支持，确保了计算的高效性和准确性。这些理论、方法和工具的综合运用，为研究Sm-CeO₂(111)表面CO分子催化氧化机理奠定了坚实的基础。三、CeO₂模型建立、计算及验证3.1CeO₂(111)表面计算模型为了深入研究Sm-CeO₂(111)表面CO分子催化氧化机理，采用超晶胞模型构建CeO₂(111)表面结构。CeO₂具有萤石型立方晶体结构，空间群为Fm-3m，晶格常数约为5.41Å。在构建超晶胞时，考虑到计算精度和计算量的平衡，选择了4×4×1的超晶胞，该超晶胞包含64个Ce原子和128个O原子。在MaterialsStudio软件中，首先构建CeO₂的立方晶胞，根据其晶体结构特点，确定Ce原子位于面心立方的顶点和面心位置，O原子位于四面体间隙位置。通过调整原子坐标和晶格参数，使其符合CeO₂的晶体结构特征。然后，利用软件的超晶胞构建功能，将立方晶胞扩展为4×4×1的超晶胞，以模拟较大尺寸的CeO₂(111)表面。为了引入Sm原子进行掺杂，选择在CeO₂(111)表面的特定Ce原子位置进行替换。在确定Sm原子的掺杂位置时，考虑了表面原子的配位环境和对称性，选择了表面层中具有代表性的Ce原子位置进行替换，以确保Sm原子的掺杂能够对表面结构和电子性质产生显著影响。通过多次测试不同的掺杂位置，最终确定了在表面层中靠近中心位置的Ce原子被Sm原子替换，此时体系的能量相对较低，结构较为稳定。在构建Sm-CeO₂(111)表面结构时，还考虑了表面的弛豫效应。由于表面原子的配位环境与体相原子不同，表面原子会发生一定程度的弛豫，以降低体系的能量。在VASP软件中，通过设置原子位置的优化参数，对表面原子的位置进行了充分的弛豫计算。在优化过程中，采用共轭梯度法对原子的位置进行迭代优化，使体系的能量和受力达到收敛标准。经过多次迭代计算，表面原子的位置逐渐稳定，体系的能量也达到了最小值。优化后的表面结构中，表面原子的键长和键角发生了一定的变化，与体相结构相比，表面原子的间距略有增大，键角也有所调整，这表明表面原子发生了明显的弛豫，形成了更加稳定的表面结构。3.2计算方法和相关公式本研究采用基于密度泛函理论（DFT）的平面波赝势方法（PWPM），运用VASP软件进行计算。平面波赝势方法将电子波函数用平面波展开，通过赝势来描述原子核与电子之间的相互作用，从而有效减少计算量。在平面波展开中，单粒子轨道波函数\psi(r)可以表示为：\psi(r)=\frac{1}{\sqrt{N}}\sum_{G}\varphi(G)\exp(i(K+G)\cdotr)其中，N为归一化因子，\varphi(G)是展开系数，K是第一Brillouin区的波矢，G是原胞的倒格矢。为了在实际计算中取有限个平面波数，引入截断能E_{cut}，满足：E_{cut}=\frac{\hbar^2(G+K)^2}{2m}对G的求和限制在(G+K)^2/2\leqE_{cut}的范围内，以保证用于展开的波函数能量小于E_{cut}，确保计算精度与效率的平衡。在本研究中，通过测试不同的截断能，确定了合适的截断能值，以保证计算结果的准确性，经过多次测试，最终选取截断能为500eV，此时体系总能量的变化量小于0.01eV/atom，满足收敛精度要求。在处理交换关联能时，采用广义梯度近似（GGA）中的PBE泛函，其交换关联能泛函E_{xc}^{GGA}[\rho]表示为：E_{xc}^{GGA}[\rho]=\int\rho(r)\epsilon_{xc}(\rho(r),\nabla\rho(r))dr其中，\epsilon_{xc}(\rho(r),\nabla\rho(r))是依赖于电子密度和密度梯度的交换关联能密度，考虑了电子密度梯度的修正，能够更准确地描述电子之间的相互作用。对于Sm-CeO₂体系，由于Ce原子存在4f电子，存在较强的电子关联效应，因此采用DFT+U方法进行计算，在能量泛函中添加与Hubbard参数U相关的修正项E_{U}：E_{total}=E_{DFT}+E_{U}其中，E_{DFT}是传统DFT计算得到的能量，E_{U}用于描述电子之间的强关联作用，通过合理设置Hubbard参数U，可以更准确地描述体系的电子结构和物理性质。在本研究中，通过参考相关文献和测试计算，确定了Hubbard参数U的值为3.5eV，使得计算结果与实验值和其他高精度理论计算结果更为吻合。在结构优化过程中，以体系的总能量和原子受力作为收敛判据。总能量的收敛标准设置为小于1×10⁻⁵eV/atom，原子受力的收敛标准设置为小于0.01eV/Å。在优化过程中，采用共轭梯度法对原子的位置进行迭代优化，使体系的能量和受力达到收敛标准。通过多次迭代计算，表面原子的位置逐渐稳定，体系的能量也达到了最小值，从而得到稳定的Sm-CeO₂(111)表面结构。3.3CeO₂(111)表面相关计算结果3.3.1氧空位形成能在CeO₂(111)表面，氧空位的形成能是衡量其稳定性的关键指标。通过第一性原理计算，得到了不同位置氧空位的形成能。在表面第一层的O原子位置形成氧空位时，计算得到的形成能为2.56eV；而在第二层O原子位置形成氧空位时，形成能为3.12eV。这表明表面第一层的氧空位相对更容易形成，稳定性较高。从晶体结构的角度来看，表面第一层的O原子与Ce原子的配位环境与体相中的O原子不同，其配位不饱和程度更高，与周围Ce原子的化学键相对较弱。当形成氧空位时，表面第一层的O原子脱离所需克服的能量较小，因此形成能较低，氧空位更易形成且相对稳定。而第二层的O原子受到周围原子的束缚更强，形成氧空位时需要更高的能量来破坏其与周围原子的相互作用，所以形成能较高，氧空位的稳定性相对较低。从电子结构的角度分析，氧空位的形成会导致周围电子云分布的改变。在表面第一层形成氧空位时，由于其配位不饱和，周围Ce原子的电子云更容易发生重排，以补偿氧空位带来的电荷变化，从而降低体系的能量，使得氧空位相对稳定。而在第二层形成氧空位时，周围原子的电子云重排受到更多限制，难以有效补偿电荷变化，导致体系能量升高，氧空位的稳定性降低。不同位置氧空位形成能的差异，对CeO₂(111)表面的化学反应活性有着重要影响。表面第一层的氧空位由于其较低的形成能和较高的稳定性，更容易成为化学反应的活性位点。在CO催化氧化反应中，CO分子更容易吸附在表面第一层的氧空位附近，与周围的活性氧物种发生反应，从而促进CO的氧化。而第二层的氧空位由于其稳定性较低，在反应过程中可能需要更高的能量才能参与反应，对CO催化氧化反应的贡献相对较小。3.3.2几何结构变化当CeO₂(111)表面形成氧空位后，其几何结构发生了显著变化。在氧空位形成前，CeO₂(111)表面的Ce-O键长约为2.34Å，O-O键长约为3.26Å。形成氧空位后，以表面第一层氧空位为例，其周围的Ce-O键长发生了明显的收缩，与氧空位相邻的Ce-O键长缩短至2.28Å，这是因为氧空位的出现使得Ce原子周围的电子云密度发生变化，Ce原子对剩余O原子的吸引力增强，导致键长缩短。同时，与氧空位不直接相邻但处于同一层的Ce-O键长则略有伸长，达到2.38Å，这是由于电子云的重新分布以及原子间的相互作用调整所引起的。从表面原子的位移情况来看，与氧空位直接相邻的Ce原子向氧空位方向发生了位移，位移距离约为0.06Å，这种位移使得Ce原子与氧空位的距离更近，增强了Ce原子与氧空位之间的相互作用。而表面其他位置的Ce原子也发生了微小的位移，以适应氧空位形成后体系的结构变化。这些原子的位移和键长的变化，导致表面原子的配位环境发生改变，原本与氧原子形成配位的Ce原子，由于氧空位的存在，配位原子数目减少，配位结构变得更加不饱和。几何结构的变化对CeO₂(111)表面的电子结构和化学反应活性产生了重要影响。键长的变化改变了原子间的电子云重叠程度，进而影响了电子的分布和能级结构。Ce-O键长的收缩使得Ce原子与O原子之间的电子云重叠增加，电子云更加集中在Ce-O键区域，导致Ce原子的电子云密度降低，O原子的电子云密度升高。这种电子结构的变化使得表面原子的化学活性发生改变，Ce原子由于电子云密度降低，更容易失去电子，表现出更强的氧化性；而O原子由于电子云密度升高，更容易接受电子，在化学反应中更易参与氧化还原过程。原子的位移和配位环境的改变也会影响表面对反应物分子的吸附能力和吸附构型。表面原子的配位不饱和使得表面更容易吸附反应物分子，且吸附构型会根据表面原子的位移和配位环境的变化而调整，从而影响化学反应的路径和速率。在CO吸附过程中，由于表面几何结构的变化，CO分子可能会以不同的吸附构型吸附在表面，不同的吸附构型会导致CO分子与表面原子之间的相互作用强度不同，进而影响CO的活化和氧化反应的进行。3.3.3电子结构变化氧空位的形成对CeO₂(111)表面的电子结构产生了显著影响，这可以通过态密度（DOS）和电荷分布的变化来体现。在未形成氧空位的CeO₂(111)表面，Ce的4f电子态主要分布在费米能级以下，呈现出一定的局域性。O的2p电子态与Ce的5d电子态在费米能级附近存在一定的杂化，形成了较为复杂的电子结构。当表面形成氧空位后，在态密度图中，费米能级附近出现了新的电子态。这是由于氧空位的存在，使得周围Ce原子的电子云分布发生变化，Ce原子的4f电子态向费米能级移动，部分4f电子的能量升高，进入到费米能级附近，形成了新的电子态。这些新的电子态具有较高的活性，为表面化学反应提供了更多的电子参与途径。从电荷分布来看，氧空位形成后，周围Ce原子的电荷发生了明显的重新分布。与氧空位相邻的Ce原子失去部分电荷，电荷数从原来的+3.65降低到+3.48，这表明Ce原子的氧化态发生了变化，部分Ce⁴⁺被还原为Ce³⁺。这是因为氧空位的出现导致电荷不平衡，为了保持体系的电中性，周围的Ce原子会失去电子，从而使自身的氧化态降低。而氧空位周围的O原子则获得部分电荷，电荷数从原来的-1.82增加到-1.75，这使得O原子的电子云密度增加，化学活性增强。这种电子结构的变化对CeO₂(111)表面的化学反应活性有着重要影响。费米能级附近新电子态的出现，使得表面更容易与反应物分子发生电子转移，促进化学反应的进行。在CO催化氧化反应中，CO分子可以与表面费米能级附近的电子态相互作用，接受电子而被活化，从而更容易与表面的活性氧物种发生反应生成CO₂。Ce原子氧化态的变化以及O原子电荷的增加，改变了表面原子的化学活性和反应活性位点的分布。Ce³⁺的出现增加了表面的还原性，使得表面更容易提供电子，促进氧化还原反应的进行；而O原子电荷的增加使其更容易参与亲电反应，与CO分子中的C原子发生反应，推动CO的氧化过程。3.4Sm负载CeO₂(111)表面的第一性原理研究3.4.1Sm/CeO₂(111)负载模型的建立及计算方法为深入探究Sm对CeO₂(111)表面CO分子催化氧化性能的影响，构建了Sm负载在CeO₂(111)表面的模型。在MaterialsStudio软件中，基于已优化的CeO₂(111)表面4×4×1超晶胞模型，选择表面层的一个Ce原子位置，将其替换为Sm原子，以此构建Sm/CeO₂(111)负载模型。在选择Sm原子替换的Ce原子位置时，考虑了表面原子的配位环境和对称性，优先选择了表面层中配位不饱和程度较高且处于对称位置的Ce原子，这样的选择能够使Sm原子的引入对表面电子结构和化学活性产生更显著的影响。替换后，对模型进行结构优化，确保模型的稳定性。在VASP软件的计算设置方面，采用平面波赝势方法，平面波截断能设置为500eV，这是通过前期的收敛性测试确定的，在此截断能下，体系总能量的变化量小于0.01eV/atom，能够保证计算结果的准确性。采用广义梯度近似（GGA）中的PBE泛函来描述电子间的交换关联能，以准确描述电子之间的相互作用。考虑到Ce原子的4f电子存在强关联效应，采用DFT+U方法进行计算，设置Hubbard参数U为3.5eV，该值是参考相关文献并结合本体系的特点，通过多次测试计算确定的，能够较好地描述Ce原子的电子结构。在结构优化过程中，以体系的总能量和原子受力作为收敛判据，总能量的收敛标准设置为小于1×10⁻⁵eV/atom，原子受力的收敛标准设置为小于0.01eV/Å。采用共轭梯度法对原子的位置进行迭代优化，使体系的能量和受力达到收敛标准。在优化过程中，对Sm原子、与Sm原子相邻的Ce原子和O原子的位置进行了充分的优化，确保这些原子的位置达到最稳定状态，以准确描述Sm负载后表面结构的变化。3.4.2氧空位形成能的变化及分析通过第一性原理计算，对比了Sm负载前后CeO₂(111)表面氧空位形成能的变化。在未负载Sm的CeO₂(111)表面，表面第一层氧空位的形成能为2.56eV。当在CeO₂(111)表面负载Sm后，相同位置氧空位的形成能降低至2.12eV。这表明Sm的负载显著降低了氧空位的形成能，使氧空位更容易形成。从晶体结构角度分析，Sm的离子半径（1.08Å）大于Ce的离子半径（0.97Å），Sm原子替换Ce原子后，会导致周围晶格发生畸变，使周围原子的配位环境发生变化。这种晶格畸变和配位环境的改变，削弱了Sm周围O原子与其他原子之间的相互作用，使得O原子更容易脱离晶格，从而降低了氧空位的形成能。从电子结构角度来看，Sm的掺杂会引起电子云分布的变化。Sm原子的电子结构与Ce原子不同，其外层电子的分布和能级与Ce原子存在差异。Sm原子的引入使得表面电子云重新分布，部分电子从Sm原子周围转移到O原子周围，导致O原子与周围原子之间的化学键减弱，从而降低了氧空位的形成能。Sm原子的电子转移还会影响表面的电荷分布，形成局部的电荷不平衡，进一步促进了氧空位的形成。氧空位形成能的降低对Sm-CeO₂(111)表面的催化性能具有重要影响。氧空位作为表面的活性位点，其浓度的增加和形成能的降低，使得表面更容易吸附和活化反应物分子。在CO催化氧化反应中，氧空位能够吸附O₂分子并使其解离，生成活性氧物种，这些活性氧物种能够与CO分子发生反应，促进CO的氧化。较低的氧空位形成能意味着在相同的反应条件下，表面能够形成更多的氧空位，从而提供更多的活性位点，提高CO催化氧化的反应速率和效率。3.5本章小结本章运用第一性原理计算方法，对CeO₂(111)表面及Sm负载后表面的结构和电子性质进行了深入研究。通过构建CeO₂(111)表面的超晶胞模型，并对其进行结构优化，得到了稳定的表面结构。计算结果表明，CeO₂(111)表面的氧空位形成能与位置密切相关，表面第一层氧空位的形成能为2.56eV，相对较低，更易形成且稳定性较高。氧空位的形成导致表面几何结构发生显著变化，Ce-O键长和表面原子的配位环境改变，这对表面的电子结构产生了重要影响，使费米能级附近出现新的电子态，Ce原子的氧化态发生变化，O原子的电荷增加，从而改变了表面的化学反应活性。在研究Sm负载CeO₂(111)表面时，构建了Sm/CeO₂(111)负载模型，并采用平面波赝势方法和DFT+U计算，设置了合适的计算参数。计算发现，Sm的负载显著降低了氧空位的形成能，从2.56eV降至2.12eV，这是由于Sm的离子半径较大，引起晶格畸变和电子云分布变化，削弱了O原子与其他原子之间的相互作用。氧空位形成能的降低使得表面更容易形成氧空位，增加了表面的活性位点，这对Sm-CeO₂(111)表面的催化性能具有重要意义，为后续研究CO在该表面的吸附和氧化反应奠定了基础。四、Sm掺杂改性CeO₂(111)表面的第一性原理研究4.1Sm₀.₀₈Ce₀.₉₂O₂(111)及其还原表面模型的建立及特点4.1.1Sm掺杂的Sm₀.₀₈Ce₀.₉₂O₂(111)表面为深入探究Sm掺杂对CeO₂(111)表面结构和性能的影响，构建了Sm₀.₀₈Ce₀.₉₂O₂(111)表面模型。