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

基于第一性原理探究Sm-CeO₂(111)表面CO分子催化氧化的微观机制一、引言1.1研究背景与意义在现代工业与科技迅猛发展的进程中，催化领域始终是化学、材料科学等多学科交叉研究的核心前沿。CeO₂作为一种极为关键的稀土氧化物，凭借其独特的物理化学性质，在催化领域占据着举足轻重的地位。CeO₂的晶体结构为立方萤石型，这种结构赋予了它特殊的电子特性，尤其是Ce元素具备+3和+4两种稳定氧化态，使得CeO₂在不同环境条件下能够便捷地发生Ce³⁺与Ce⁴⁺之间的氧化还原反应，进而展现出卓越的氧化还原性能。CeO₂在催化领域的应用极为广泛，在汽车尾气净化催化剂中，它是不可或缺的关键组分。随着汽车保有量的持续攀升，汽车尾气已成为大气污染的主要来源之一，其中包含的一氧化碳（CO）、碳氢化合物（HC）和氮氧化物（NOx）等污染物对环境和人体健康危害巨大。CeO₂能够利用自身的氧化还原性能，储存和释放氧，有效调节尾气中氧气的浓度，从而有力地促进CO、HC和NOx等污染物发生氧化还原反应，将其转化为无害的二氧化碳（CO₂）、水（H₂O）和氮气（N₂），显著降低汽车尾气对环境的污染。在燃料电池中，CeO₂既可以充当催化剂的载体，为活性组分提供高比表面积的支撑，增强活性组分的分散性和稳定性；也能够直接作为催化剂的活性组分，积极参与电极反应中的氧还原反应和氢氧化反应，大幅提高燃料电池的性能和效率，推动清洁能源的高效利用。此外，在工业废气处理、化工合成等诸多领域，CeO₂同样发挥着至关重要的作用，展现出巨大的应用潜力和价值。CO催化氧化反应在环境保护、能源领域以及工业生产等方面都具有至关重要的意义，一直是催化领域的研究热点。在环境保护层面，CO是一种对人体健康和生态环境危害严重的有毒气体。它主要来源于冶金、化工等行业化石燃料的不完全燃烧、垃圾焚烧以及汽车尾气排放等。据相关数据显示，我国工业排放的CO中，钢铁行业的贡献占比较大，其中烧结工序的排放占比高达40%。CO与人体血液中的血红蛋白具有极强的亲和力，其结合能力远高于氧气与血红蛋白的结合能力，一旦CO进入人体，会迅速与血红蛋白结合，形成碳氧血红蛋白，阻碍氧气的输送，导致人体组织缺氧，引发头晕、恶心、呕吐等中毒症状，严重时甚至会危及生命。此外，CO在大气中还会参与光化学反应，产生一系列二次污染物，如臭氧（O₃）等，进一步加剧大气污染，破坏生态平衡。通过高效的CO催化氧化技术，将CO转化为CO₂，可以显著减少CO的排放，有效改善空气质量，保护生态环境，保障人类健康。从能源角度来看，CO催化氧化在能源转化和利用过程中扮演着关键角色。在固体氧化物燃料电池（SOFC）中，CO是一种重要的燃料来源。然而，未经催化氧化的CO直接参与电池反应时，反应效率较低，且容易在电极表面产生积碳，导致电池性能下降和寿命缩短。通过CO催化氧化反应，可以将CO高效地转化为CO₂，并释放出电能，提高燃料电池的能量转换效率，促进清洁能源的开发和利用，缓解能源危机。在一些化工合成过程中，CO也是重要的原料之一。例如，在合成甲醇等有机化学品的反应中，CO催化氧化反应的高效进行能够提高原料的利用率，降低生产成本，提升化工生产的经济效益和可持续性。在众多CeO₂基催化剂的研究中，Sm-CeO₂体系因其独特的优势受到了广泛关注。Sm作为一种稀土元素，其离子半径与Ce离子半径相近，能够较为容易地进入CeO₂晶格中，形成固溶体。这种固溶体结构的形成会导致CeO₂晶格发生畸变，产生更多的氧空位。氧空位在催化反应中具有至关重要的作用，它可以作为活性位点，增强对反应物分子的吸附和活化能力，从而显著提高催化剂的活性和选择性。研究Sm-CeO₂(111)表面CO分子催化氧化机理，对于深入理解该催化剂的催化过程具有不可替代的重要性。通过揭示CO分子在Sm-CeO₂(111)表面的吸附、活化以及与氧气发生反应的详细步骤和微观机制，可以明确催化剂表面活性位点的性质和作用，了解反应过程中电子转移、化学键的形成与断裂等关键信息，为从原子和分子层面认识催化反应提供理论依据。这有助于我们深入理解催化过程的本质，为开发新型高效的催化剂提供坚实的理论基础。研究Sm-CeO₂(111)表面CO分子催化氧化机理对于优化催化剂性能具有重要的指导意义。基于对催化机理的深入认识，可以有针对性地对催化剂进行设计和改性。例如，通过调控Sm的掺杂量、优化催化剂的制备工艺等手段，进一步优化催化剂的结构和性能，提高其对CO催化氧化的活性、选择性和稳定性，降低催化剂的成本，推动CO催化氧化技术在实际应用中的广泛推广和发展。在实际应用中，高效的CO催化氧化催化剂可以应用于汽车尾气净化装置、工业废气处理设备以及燃料电池等领域，为解决环境污染问题和促进能源的高效利用提供有力的技术支持，具有重要的现实意义和应用价值。1.2CeO₂的性质及相关应用1.2.1CeO₂的结构与性质CeO₂属于立方晶系，具有典型的萤石型结构。在这种结构中，Ce原子位于面心立方的晶格节点位置，每个Ce原子被周围8个O原子以立方体的方式配位，形成Ce-O₈多面体；而O原子则填充在四面体空隙中，每个O原子与周围4个Ce原子相连。这种有序的原子排列方式赋予了CeO₂独特的物理化学性质。从晶体结构参数来看，CeO₂的晶格常数约为0.541nm，这种特定的晶格常数对其内部原子间的相互作用和电子云分布产生影响，进而影响材料的整体性能。在物理性质方面，CeO₂通常呈现为淡黄色或白色的粉末状固体，具有较高的熔点，约为2600℃，这使得它在高温环境下仍能保持稳定的固态结构，为其在高温催化等领域的应用提供了基础。CeO₂还具有适中的硬度，莫氏硬度约为6-7，这种硬度特性使其在一些需要耐磨性能的应用场景中，如抛光材料等方面具有一定的优势。在电学性能上，CeO₂具有一定的离子导电性，尤其是在高温和存在氧空位的情况下，氧离子可以通过氧空位在晶格中进行迁移，从而实现离子导电，这一特性在固体氧化物燃料电池等领域有着重要的应用价值。从化学性质角度，Ce元素存在+3和+4两种稳定的氧化态，这使得CeO₂在不同的化学环境下能够发生Ce³⁺与Ce⁴⁺之间的氧化还原反应。在氧化性氛围中，Ce³⁺可以失去电子被氧化为Ce⁴⁺；而在还原性氛围中，Ce⁴⁺又能够获得电子被还原为Ce³⁺。这种可逆的氧化还原特性使得CeO₂在催化领域展现出独特的性能，例如在氧化还原反应中可以作为电子传递的媒介，促进反应物的氧化或还原过程。CeO₂具有较好的化学稳定性，在一般的酸碱环境中，不易与常见的酸碱发生化学反应，能够保持自身的结构和化学组成稳定，这为其在复杂化学体系中的应用提供了保障。1.2.2CeO₂在汽车尾气催化方面的应用及原理随着全球汽车保有量的持续增长，汽车尾气排放已成为大气污染的主要来源之一。汽车尾气中含有一氧化碳（CO）、碳氢化合物（HC）、氮氧化物（NOx）等多种有害污染物，这些污染物对环境和人体健康造成了严重的危害。CO是一种无色、无味的有毒气体，它与人体血液中的血红蛋白具有极强的亲和力，会导致人体缺氧中毒；HC是形成光化学烟雾的重要前体物，会引发呼吸道疾病等问题；NOx则会导致酸雨、臭氧层破坏等环境问题。因此，汽车尾气的净化对于环境保护和人类健康至关重要。CeO₂在汽车尾气净化催化剂中扮演着核心角色，是三元催化剂（TWC）的关键组成部分。三元催化剂能够同时对CO、HC和NOx三种污染物进行催化转化，将它们转化为无害的CO₂、H₂O和N₂。CeO₂在其中主要发挥了储氧和调节氧浓度的作用，以及利用其氧化还原特性促进反应进行。在汽车发动机的工作过程中，尾气中的氧含量会随着工况的变化而发生波动。当发动机处于富燃状态时，尾气中氧气含量较低，此时CeO₂中的Ce⁴⁺会捕获尾气中的氧原子，被还原为Ce³⁺，从而储存氧：CeO₂+2e⁻+2H⁺→Ce₂O₃+H₂O；当发动机处于贫燃状态时，尾气中氧气含量较高，Ce³⁺又会释放出储存的氧原子，被氧化为Ce⁴⁺：Ce₂O₃+1/2O₂→2CeO₂。通过这种储氧和释氧的循环过程，CeO₂能够有效地调节尾气中的氧浓度，使其始终保持在一个合适的范围内，为CO、HC和NOx的催化转化提供稳定的反应环境。在CO的催化氧化反应中，CeO₂表面的活性位点可以吸附CO分子，使其发生活化。同时，CeO₂中的晶格氧也能够参与反应，将CO氧化为CO₂：CO+O→CO₂（其中O表示CeO₂表面的晶格氧）。在这个过程中，CeO₂的氧化还原循环起到了关键作用。当晶格氧参与反应后，Ce⁴⁺被还原为Ce³⁺，形成氧空位；随后，气相中的氧气分子会吸附在氧空位上，发生解离并填充氧空位，使Ce³⁺重新被氧化为Ce⁴⁺，从而完成一个氧化还原循环。这个循环过程不断进行，实现了CO的持续催化氧化。在HC的催化氧化反应中，CeO₂同样发挥了重要作用。HC分子首先吸附在CeO₂表面的活性位点上，被活化后的HC分子与CeO₂表面的晶格氧或气相中的氧气发生反应，逐步被氧化为CO₂和H₂O。CeO₂的储氧能力可以在反应过程中提供额外的氧源，促进HC的完全氧化，提高净化效率。在NOx的还原反应中，CeO₂可以通过与贵金属（如Pt、Pd、Rh等）协同作用，促进NOx的还原。CeO₂表面的活性位点能够吸附NOx分子，使其发生解离和活化；同时，CeO₂提供的氧空位和调节氧浓度的能力，有助于维持反应过程中的氧化还原平衡，促进NOx与CO、HC等还原剂之间的反应，将NOx还原为N₂。1.2.3CO的催化氧化机理概述CO的催化氧化反应是一个复杂的多步骤过程，涉及到反应物在催化剂表面的吸附、活化、表面物种的相互作用以及产物的脱附等多个环节。目前，关于CO催化氧化机理的研究主要基于几种经典的反应模型，包括Langmuir-Hinshelwood（L-H）机理、Eley-Rideal（E-R）机理和Mars-vanKrevelen（M-vK）机理。在L-H机理中，CO和O₂分子首先分别吸附在催化剂表面的不同活性位点上，形成吸附态的CO和O₂（表示吸附态）。随后，吸附态的O₂发生解离，形成两个吸附态的O原子。吸附态的CO和O在催化剂表面发生反应，生成吸附态的CO₂，最后CO₂从催化剂表面脱附，进入气相。整个反应过程可以表示为：CO+⇌CO，O₂+2⇌2O*，CO*+O→CO₂，CO₂*⇌CO₂+*。在这个机理中，反应速率主要取决于表面吸附物种之间的反应步骤，并且要求催化剂表面具有足够多的相邻活性位点来吸附CO和O原子。E-R机理则认为，气相中的CO分子直接与吸附在催化剂表面的O原子发生反应，生成CO₂。即CO+O*→CO₂。在这个机理中，反应不需要CO分子在催化剂表面进行吸附和活化，而是直接与表面的活性氧物种发生反应，反应速率主要取决于气相CO分子与表面O原子的碰撞频率和反应活性。M-vK机理强调催化剂晶格氧的参与。在反应过程中，催化剂表面的晶格氧与CO分子发生反应，将CO氧化为CO₂，同时催化剂晶格中产生氧空位。随后，气相中的O₂分子吸附在氧空位上，发生解离并填充氧空位，使催化剂恢复到初始状态。反应过程可以表示为：CO+Oₗ→CO₂+Vₒ（Oₗ表示晶格氧，Vₒ表示氧空位），1/2O₂+Vₒ+2e⁻→Oₗ。在这个机理中，催化剂的氧化还原性能和氧空位的形成与迁移能力对反应速率起着关键作用。实际的CO催化氧化过程可能是多种机理共同作用的结果，具体取决于催化剂的种类、表面结构、反应条件（如温度、压力、气体组成等）。不同的催化剂表面具有不同的活性位点和电子结构，会影响反应物的吸附和活化方式，从而导致不同的反应机理占主导地位。在低温条件下，L-H机理可能更为常见，因为此时反应物分子的扩散和反应活性较低，需要在催化剂表面进行充分的吸附和活化；而在高温条件下，M-vK机理可能更容易发生，因为高温有利于晶格氧的迁移和氧空位的形成。1.2.4CeO₂在其他方面的应用除了在汽车尾气催化领域的重要应用外，CeO₂凭借其独特的物理化学性质，在燃料电池、化学传感器等多个领域也展现出了广泛的应用前景。在燃料电池领域，CeO₂作为固体氧化物燃料电池（SOFC）的关键材料，发挥着至关重要的作用。在SOFC中，CeO₂既可以作为电解质材料，又可以作为电极催化剂的组成部分。作为电解质，CeO₂具有较高的氧离子电导率，尤其是在中低温范围内（500-700℃）。在这个温度区间，CeO₂晶格中的氧离子能够在电场的作用下通过氧空位进行迁移，实现电荷的传导，从而保证电池的正常工作。与传统的氧化钇稳定的氧化锆（YSZ）电解质相比，CeO₂基电解质具有更低的工作温度要求，这有助于降低电池的制造成本和提高电池的启动速度。在电极催化剂方面，CeO₂可以与贵金属（如Pt、Pd等）或过渡金属氧化物复合，形成高性能的电极催化剂。CeO₂的氧化还原性能和储氧能力能够促进电极反应中的氧还原反应（ORR）和氢氧化反应（HOR）。在阴极，CeO₂能够吸附和活化氧气分子，促进氧分子的解离和氧离子的形成，加速氧还原反应的进行；在阳极，CeO₂可以提高燃料（如氢气、一氧化碳等）的氧化活性，增强电极的催化性能，提高燃料电池的整体效率。在化学传感器领域，CeO₂常被用于制备气体传感器，用于检测环境中的有害气体，如CO、NOx、H₂S等。CeO₂基气体传感器的工作原理主要基于其表面与气体分子之间的相互作用引起的电学性能变化。当目标气体分子吸附在CeO₂表面时，会与表面的活性位点发生化学反应，导致CeO₂表面的电子结构发生改变，从而引起其电阻、电容等电学参数的变化。通过检测这些电学参数的变化，就可以实现对目标气体浓度的检测。在检测CO气体时，CO分子在CeO₂表面被氧化为CO₂，同时CeO₂表面的电子会转移给CO分子，导致CeO₂的电阻发生变化。通过测量电阻的变化量，就可以定量地确定CO的浓度。CeO₂的氧化还原性能和对气体分子的吸附选择性，使得它在气体传感器领域具有较高的灵敏度和选择性，能够快速、准确地检测出低浓度的有害气体，为环境监测和工业安全生产提供了重要的技术支持。CeO₂在光催化领域也有一定的应用。在光催化反应中，CeO₂能够吸收特定波长的光，产生电子-空穴对。这些电子和空穴具有较强的氧化还原能力，可以参与到有机污染物的降解反应中。在光照条件下，CeO₂表面产生的空穴可以将吸附在其表面的有机污染物氧化为二氧化碳和水等无害物质，而电子则可以与氧气分子反应生成具有强氧化性的超氧自由基等活性氧物种，进一步促进有机污染物的降解。CeO₂的光催化性能与其晶体结构、表面缺陷、氧空位等因素密切相关，通过对这些因素的调控，可以提高CeO₂的光催化活性，实现对有机污染物的高效降解，在污水处理、空气净化等环境治理领域具有潜在的应用价值。1.3国内外研究现状及发展趋势1.3.1负载或掺杂金属原子提高CeO₂催化活性为了进一步提升CeO₂的催化活性，科研人员在负载或掺杂金属原子方面开展了大量研究。常见的负载金属有贵金属（如Pt、Pd、Rh等）和非贵金属（如Fe、Co、Ni等），掺杂金属则包括稀土金属（如Sm、Gd、Y等）以及过渡金属（如Ti、Zr、Mn等）。贵金属负载在CeO₂上能够显著提高其催化活性。Pt/CeO₂催化剂在CO催化氧化反应中表现出优异的性能，Pt的高催化活性与CeO₂的储氧能力相结合，促进了CO和O₂在催化剂表面的吸附和活化。Pt原子可以作为活性位点，吸附CO分子，使其发生活化；同时，CeO₂表面的氧空位和晶格氧能够与Pt协同作用，加速O₂的解离和参与反应的过程，从而提高CO催化氧化的反应速率。研究表明，在低温条件下，Pt/CeO₂催化剂能够实现CO的高效转化，其起燃温度明显低于纯CeO₂催化剂。非贵金属负载在CeO₂上也能提升其催化性能，且具有成本优势。Fe负载在CeO₂上形成的Fe-CeO₂催化剂，在CO催化氧化反应中，Fe的引入改变了CeO₂的电子结构，增加了表面氧空位的浓度，从而提高了催化剂对CO和O₂的吸附能力和反应活性。Fe的存在可以促进CeO₂表面Ce³⁺/Ce⁴⁺氧化还原循环的进行，使得晶格氧的迁移和参与反应更加容易，进而提高CO的氧化效率。在一些研究中，通过优化Fe的负载量和制备工艺，Fe-CeO₂催化剂在中高温条件下展现出良好的CO催化氧化性能，具有较高的转化率和稳定性。在掺杂金属方面，稀土金属Sm掺杂CeO₂形成的Sm-CeO₂固溶体是研究的热点之一。Sm的掺杂导致CeO₂晶格发生畸变，产生更多的氧空位，增强了催化剂对反应物的吸附和活化能力。氧空位可以作为活性位点，吸附和活化O₂分子，使其更容易解离为活性氧原子，同时也有利于CO分子的吸附和反应。研究发现，适量的Sm掺杂能够显著提高CeO₂在CO催化氧化反应中的活性，降低反应的活化能。通过第一性原理计算和实验表征相结合的方法，深入研究了Sm-CeO₂体系中氧空位的形成机制、电子结构变化以及对CO催化氧化反应路径的影响，为进一步优化催化剂性能提供了理论依据。过渡金属Zr掺杂CeO₂也能有效改善其催化性能。Zr的掺杂可以稳定CeO₂的晶格结构，提高其热稳定性；同时，Zr的存在会改变CeO₂表面的电子云分布，影响表面氧物种的性质和活性，从而对CO催化氧化反应产生积极影响。Zr掺杂的CeO₂催化剂在高温条件下仍能保持较好的催化活性和稳定性，在一些工业催化过程中具有潜在的应用价值。通过控制Zr的掺杂量和制备条件，可以调控催化剂的结构和性能，满足不同反应体系的需求。1.3.2催化反应过程机理的研究方法与进展在研究CO催化氧化机理方面，实验和理论计算都发挥着不可或缺的作用，并且取得了一系列重要进展。实验研究方法主要包括原位表征技术和动力学研究。原位表征技术能够在反应条件下对催化剂的结构、表面物种和电子状态等进行实时监测，为揭示催化反应机理提供直接的实验证据。X射线光电子能谱（XPS）可以分析催化剂表面元素的化学态和电子结合能，从而确定反应过程中表面物种的变化。在CO催化氧化反应中，通过XPS可以监测到CeO₂表面Ce³⁺/Ce⁴⁺的比例变化，以及CO、O₂等反应物和CO₂等产物在催化剂表面的吸附和反应情况，了解反应过程中电子转移和化学键的形成与断裂。傅里叶变换红外光谱（FT-IR）则可以检测催化剂表面吸附物种的振动频率，从而确定表面物种的种类和结构。利用FT-IR可以观察到CO在催化剂表面的吸附形态，如线性吸附、桥式吸附等，以及吸附态CO与表面氧物种之间的反应过程，为研究反应机理提供详细的信息。动力学研究通过测量反应速率与反应物浓度、温度等因素之间的关系，建立反应动力学模型，从而推断反应的速率控制步骤和反应机理。通过改变CO和O₂的浓度，测量不同条件下的反应速率，利用动力学方程拟合得到反应的速率常数和活化能，进而确定反应的速率控制步骤是CO的吸附、O₂的解离还是表面反应等过程。结合不同催化剂体系和反应条件下的动力学研究结果，可以总结出CO催化氧化反应的一般规律和影响因素，为优化催化剂性能和反应条件提供指导。理论计算方法主要包括密度泛函理论（DFT）计算和分子动力学（MD）模拟。DFT计算能够从原子和分子层面研究催化剂表面的电子结构、吸附能、反应路径和活化能等，为深入理解催化反应机理提供微观层面的理论支持。通过构建CeO₂表面模型，利用DFT计算CO和O₂在CeO₂表面的吸附能，确定最稳定的吸附位点和吸附构型；计算反应过程中各基元步骤的活化能和反应热，从而确定反应的最优路径和速率控制步骤。在研究Sm-CeO₂(111)表面CO催化氧化机理时，DFT计算可以揭示Sm掺杂对CeO₂表面电子结构的影响，以及这种影响如何改变CO和O₂的吸附和反应行为，为解释实验现象提供理论依据。MD模拟则可以在原子尺度上模拟催化反应的动态过程，研究反应体系中原子的运动轨迹、分子间的相互作用以及反应过程中的能量变化。通过MD模拟可以观察到CO和O₂分子在催化剂表面的扩散、吸附和反应过程，以及反应过程中催化剂表面结构的动态变化，为理解催化反应的微观动力学提供直观的图像。将DFT计算和MD模拟相结合，可以更全面、深入地研究CO催化氧化反应机理，从静态和动态两个方面揭示催化过程的本质。1.3.3国内外研究CeO₂趋势目前，国内外在CeO₂基催化剂研究方面呈现出多个前沿方向和热点问题。在催化剂的设计与制备方面，研究重点逐渐转向精准调控催化剂的微观结构和组成，以实现更高的催化活性、选择性和稳定性。通过纳米技术制备具有特定形貌和尺寸的CeO₂纳米材料，如纳米棒、纳米立方体、纳米多孔结构等，这些特殊形貌的纳米材料能够暴露更多的高活性晶面，增加表面活性位点的数量，从而提高催化剂的活性。研究不同形貌的CeO₂纳米棒，发现其{110}晶面暴露比例较高时，在CO催化氧化反应中表现出更高的活性，因为{110}晶面具有特殊的原子排列和电子结构，有利于CO和O₂的吸附和活化。开发新型的复合催化剂也是研究热点之一。将CeO₂与其他功能性材料复合，如与具有高导电性的碳材料复合形成CeO₂-C复合材料，利用碳材料的高导电性和大比表面积，提高电子传输效率，增强催化剂对反应物的吸附能力；与具有特殊光学性能的半导体材料复合形成CeO₂-半导体复合材料，实现光催化与热催化的协同作用，拓展催化剂的应用范围。在光热协同催化领域，CeO₂与TiO₂复合形成的CeO₂-TiO₂催化剂，在光照和加热条件下，能够同时利用光生载流子和热活化作用，实现对CO和其他有机污染物的高效降解。在催化反应机理的深入研究方面，随着实验技术和理论计算方法的不断发展，研究人员致力于从更微观、更动态的角度揭示催化反应的本质。结合原位表征技术和理论计算，实时跟踪催化反应过程中催化剂表面的结构变化、电子转移和化学反应路径，深入理解活性位点的形成和作用机制。利用原位X射线吸收精细结构（XAFS）技术和DFT计算，研究CeO₂基催化剂在CO催化氧化反应过程中Ce原子的配位环境变化和电子结构演变，揭示活性位点的本质和反应过程中的关键步骤。探索CeO₂基催化剂在新领域的应用也是当前的研究趋势之一。除了传统的汽车尾气净化、燃料电池等领域，CeO₂基催化剂在新兴的能源存储与转化、环境修复等领域展现出潜在的应用价值。在可充电电池中，CeO₂作为电极材料或添加剂，能够改善电池的充放电性能和循环稳定性；在土壤修复领域，利用CeO₂的氧化还原性能和吸附能力，去除土壤中的重金属污染物和有机污染物，为环境保护提供新的解决方案。1.4课题提出和研究内容1.4.1课题提出及重要性在CeO₂基催化剂的研究体系中，Sm-CeO₂体系凭借其独特的结构和性能优势，成为研究CO催化氧化反应的理想模型体系之一。Sm的掺杂使CeO₂晶格产生畸变并形成大量氧空位，显著影响了催化剂表面的电子结构和化学活性，然而，目前对于Sm-CeO₂体系中Sm掺杂对CeO₂电子结构的具体影响机制，以及这些微观层面的变化如何精准调控CO和O₂分子在催化剂表面的吸附、活化及反应过程，仍缺乏深入且系统的认识。尽管已有研究涉及Sm-CeO₂的催化性能，但多集中于实验宏观表征，对于原子和分子尺度下的催化反应细节及电子转移过程的研究还不够透彻，这在一定程度上限制了对催化剂活性和选择性本质的理解，也阻碍了基于微观机制的高效催化剂的理性设计和开发。选择Sm-CeO₂(111)表面作为研究对象，主要基于晶体表面晶面效应的考量。晶体的不同晶面具有不同的原子排列和电子云分布，从而表现出各异的表面能和化学活性。CeO₂的(111)晶面是其最稳定且暴露比例较高的晶面之一，在催化反应中扮演着重要角色。研究Sm-CeO₂(111)表面，能够深入探究Sm掺杂对CeO₂最具代表性晶面的影响，为全面理解Sm-CeO₂催化剂的催化性能提供关键信息。通过第一性原理计算研究Sm-CeO₂(111)表面CO分子催化氧化机理，具有多方面的创新点和重要意义。从研究方法上，第一性原理计算能够从原子和电子层面深入剖析催化反应过程，避免了传统实验研究难以触及微观细节的局限性，为揭示催化反应本质提供了原子和分子水平的理论视角。通过计算可以精确获取CO和O₂在Sm-CeO₂(111)表面的吸附能、吸附构型以及反应过程中的活化能和反应热等关键参数，为深入理解催化反应路径和速率控制步骤提供定量依据。从研究内容上，本研究聚焦于Sm-CeO₂(111)表面这一特定体系，系统研究Sm掺杂量、氧空位浓度等因素对CO催化氧化性能的影响规律，以及这些因素与催化剂表面电子结构和反应活性之间的内在联系。这种对特定体系和关键影响因素的深入研究，有助于填补该领域在微观机制研究方面的空白，为CeO₂基催化剂的设计和优化提供更为精准的理论指导。研究Sm-CeO₂(111)表面CO分子催化氧化机理，对于拓展CeO₂基催化剂的应用领域也具有重要的推动作用。通过深入理解催化反应机理，可以为开发新型高效的CO催化氧化催化剂提供理论基础，促进其在汽车尾气净化、工业废气处理、燃料电池等领域的广泛应用，为解决环境污染和能源危机等实际问题提供有力的技术支持。1.4.2研究主要内容本研究借助第一性原理计算方法，深入系统地研究Sm-CeO₂(111)表面CO分子催化氧化机理，具体研究内容涵盖以下多个关键方面：Sm-CeO₂(111)表面结构优化与稳定性分析：构建合理的Sm-CeO₂(111)表面模型，细致考虑不同Sm掺杂量和氧空位浓度对表面结构的影响。运用第一性原理计算对表面模型进行全面优化，获取最稳定的表面原子结构和电子结构。通过计算表面形成能、结合能等关键参数，深入分析Sm-CeO₂(111)表面的稳定性，明确Sm掺杂和氧空位引入对表面稳定性的影响规律，为后续研究提供坚实的结构基础。CO和O₂在Sm-CeO₂(111)表面的吸附行为研究：精确计算CO和O₂分子在Sm-CeO₂(111)表面不同吸附位点的吸附能，确定最稳定的吸附构型。结合电子结构分析，深入探究吸附过程中分子与表面之间的电荷转移、化学键的形成与变化等微观机制。通过对比不同Sm掺杂量和氧空位浓度下的吸附能和吸附构型，揭示Sm掺杂和氧空位对CO和O₂吸附性能的影响规律，明确表面活性位点的性质和分布。CO催化氧化反应路径与机理研究：基于对表面吸附行为的深入理解，运用爬坡nudgedelasticband（CI-NEB）等方法，系统搜索和计算CO催化氧化反应在Sm-CeO₂(111)表面的可能反应路径和对应的活化能。通过分析反应过程中各基元步骤的能量变化和原子结构变化，确定反应的最优路径和速率控制步骤。结合实验研究结果，深入探讨反应机理，明确反应过程中涉及的表面物种、反应中间体以及电子转移过程，全面揭示Sm-CeO₂(111)表面CO催化氧化的微观机制。温度对CO催化氧化反应的影响研究：运用分子动力学（MD）模拟方法，在不同温度条件下对Sm-CeO₂(111)表面CO催化氧化反应进行动态模拟。通过模拟观察反应体系中原子的运动轨迹、分子的扩散和反应过程，研究温度对反应速率、产物选择性以及表面结构稳定性的影响。结合阿伦尼乌斯方程等动力学理论，建立温度与反应速率之间的定量关系，深入分析温度对反应动力学的影响机制，为实际反应条件的优化提供理论依据。与实验结果对比与验证：广泛收集相关实验数据，将本研究的理论计算结果与实验结果进行细致对比和深入分析。通过对比验证，进一步完善和修正理论模型，提高理论计算的准确性和可靠性。同时，基于理论计算结果，为实验研究提供有针对性的指导和建议，如优化催化剂的制备工艺、调整反应条件等，实现理论与实验的有机结合，共同推动对Sm-CeO₂(111)表面CO分子催化氧化机理的深入理解和高效催化剂的开发。二、第一性原理的理论基础2.1绝热近似在多粒子体系中，由于原子核的质量远大于电子的质量，其运动速度相对电子而言极为缓慢。基于此，绝热近似作出了一个关键假设：在研究电子的运动时，可以认为原子核固定在其瞬时位置上，仿佛处于静止状态。这一假设的合理性源于电子与原子核质量的巨大差异，电子质量约为9.11\times10^{-31}kg，而原子核质量通常是电子质量的数千倍甚至更高。这种质量上的悬殊导致电子的运动速度比原子核快得多，使得在电子运动的时间尺度内，原子核的位置变化可以忽略不计。从物理图像来看，电子在原子核形成的相对稳定的势场中快速运动，就像在一个固定的框架内活动。以氢原子为例，电子围绕着原子核高速旋转，而原子核几乎不动，电子的运动主要受到原子核库仑引力的作用。在多原子分子或固体材料中，虽然情况更为复杂，但绝热近似同样适用，电子在由多个原子核构成的固定势场中进行运动。基于上述假设，绝热近似将多粒子体系的哈密顿量进行了简化。在不考虑外场作用时，多粒子体系的哈密顿量H包含电子的动能项T_e、电子-电子相互作用项V_{ee}、原子核的动能项T_N、原子核-原子核相互作用项V_{NN}以及电子-原子核相互作用项V_{eN}，即H=T_e+V_{ee}+T_N+V_{NN}+V_{eN}。在绝热近似下，由于认为原子核固定不动，原子核的动能项T_N可以忽略不计，此时哈密顿量简化为H=T_e+V_{ee}+V_{NN}+V_{eN}。这种简化使得我们可以将电子的运动和原子核的运动分开考虑，极大地降低了问题的复杂性。原本需要同时处理电子和原子核的多体运动问题，现在可以转化为在固定原子核势场下研究电子的运动，将多粒子问题简化为多电子问题。在研究金属晶体中电子的运动时，通过绝热近似，我们可以将金属原子的原子核视为固定在晶格节点上，电子在这些固定原子核形成的势场中运动。这样就可以专注于电子的行为，而不必同时考虑原子核的复杂运动，从而能够运用量子力学方法对电子的能量状态、波函数等进行求解，为进一步研究金属的电学、热学等性质奠定基础。2.2单电子近似在解决多电子体系问题时，由于电子间存在着复杂的相互作用，精确求解多电子体系的薛定谔方程面临着巨大的困难。为了简化这一复杂问题，单电子近似应运而生。单电子近似的核心思想是将多电子问题转化为多个独立的单电子问题进行处理。在多电子体系中，每个电子都受到原子核的吸引以及其他电子的排斥作用，其运动状态相互关联，使得问题极为复杂。单电子近似假设每个电子都在原子核和其他电子所产生的平均势场中独立运动，这样就可以将多电子体系的哈密顿量分解为多个单电子哈密顿量之和。以晶体中的电子为例，晶体中包含大量的电子，它们之间存在着库仑相互作用。在单电子近似下，我们可以将每个电子看作是在由所有原子核和其他电子形成的平均势场中运动。这种平均势场可以看作是一个等效的周期性势场，其形式与晶体的晶格结构密切相关。通过这种近似，我们可以为每个电子构建一个独立的薛定谔方程，即：\left(-\frac{\hbar^{2}}{2m}\nabla^{2}+V_{eff}(r)\right)\psi_{i}(r)=E_{i}\psi_{i}(r)其中，\frac{\hbar^{2}}{2m}\nabla^{2}是电子的动能算符，V_{eff}(r)是平均势场，\psi_{i}(r)是第i个电子的波函数，E_{i}是第i个电子的能量。在实际应用中，实现单电子近似的一种常见方法是Hartree-Fock近似。Hartree-Fock近似首先假设多电子体系的波函数可以表示为单电子波函数的乘积形式，即：\Psi(r_{1},r_{2},\cdots,r_{N})=\psi_{1}(r_{1})\psi_{2}(r_{2})\cdots\psi_{N}(r_{N})然后，通过变分原理来确定这些单电子波函数。在这个过程中，需要考虑电子的交换反对称性，即当两个电子交换位置时，波函数会改变符号。为了满足这一性质，引入了Slater行列式来描述多电子波函数。通过求解Hartree-Fock方程，可以得到每个单电子的波函数和能量，从而实现多电子问题向单电子问题的转化。另一种实现单电子近似的重要方法是基于密度泛函理论（DFT）。DFT认为，体系的基态性质可以由基态电子密度唯一确定。通过引入交换关联泛函来描述电子之间的复杂相互作用，将多电子体系的问题转化为一个等效的单电子问题。在DFT中，Kohn-Sham方程是核心方程，它的形式与上述单电子薛定谔方程类似，但其中的势场包含了交换关联势等因素，以更准确地描述电子间的相互作用。通过求解Kohn-Sham方程，可以得到体系的电子密度和能量等性质，实现对多电子体系的有效处理。2.3密度泛函理论(Densityfunctionaltheory)2.3.1Hohenberg-Kohn定理密度泛函理论（DFT）的基础是Hohenberg-Kohn定理，这一定理建立了系统基态电子密度和外势场之间的深刻对应关系。Hohenberg-Kohn第一定理指出，对于一个处于外部势场V_{ext}(r)中的多电子体系，其基态性质可以由基态电子密度n(r)唯一确定。这里的外部势场主要来源于原子核的库仑吸引势，而电子密度则描述了电子在空间中的分布情况。这意味着，只要知道了基态电子密度，就可以确定体系的所有基态性质，如能量、电荷分布、电子波函数等。从数学角度来看，体系的哈密顿量H可以表示为电子动能项T、电子-电子相互作用项V_{ee}以及外部势场项V_{ext}的总和，即H=T+V_{ee}+V_{ext}。而根据Hohenberg-Kohn第一定理，体系的基态能量E_0是基态电子密度n(r)的泛函，记为E_0[n(r)]。Hohenberg-Kohn第二定理进一步表明，基态能量是一个泛函，它依赖于电子密度函数。对于给定的外部势场V_{ext}(r)，可以通过对电子密度函数进行泛函变分来求解体系的基态能量。具体来说，通过寻找一个合适的电子密度函数n(r)，使得泛函E[n(r)]取最小值，这个最小值就是体系的基态能量。在实际计算中，通常采用变分原理来实现这一过程。假设一个试探电子密度函数n'(r)，计算对应的能量泛函E[n'(r)]，然后通过不断调整n'(r)，使得E[n'(r)]逐渐减小，直到达到最小值，此时的n'(r)就是体系的基态电子密度，对应的能量就是基态能量。Hohenberg-Kohn定理的重要意义在于它将复杂的多体问题转化为一个只需要处理三个维度（x、y、z方向）的电子密度问题。在传统的量子力学方法中，处理多电子体系需要考虑电子之间复杂的相互作用，其波函数是一个包含多个电子坐标的高维函数，计算量极其庞大。而根据Hohenberg-Kohn定理，我们可以将注意力集中在电子密度这一三维函数上，大大简化了电子结构的计算。这使得密度泛函理论在处理大分子、固体材料等多电子体系时具有显著的优势，能够在可接受的计算资源下得到较为准确的结果。2.3.2Ko-Sham方程Kohn-Sham方程是密度泛函理论的核心方程，它将多电子体系的问题巧妙地转化为一个等效的单电子问题。其推导过程基于Hohenberg-Kohn定理，旨在寻找一种有效的方法来求解体系的电子密度。在多电子体系中，体系的总能量E可以表示为电子动能项T、电子-电子相互作用项（库仑相互作用项V_{ee}和交换关联项E_{xc}）以及外部势场项V_{ext}的总和，即E=T+V_{ee}+E_{xc}+V_{ext}。为了简化问题，Kohn和Sham引入了一个虚构的无相互作用电子体系，这个体系的电子在一个有效势场V_{eff}(r)中运动。这个有效势场包含了外部势场V_{ext}(r)、Hartree势V_{H}(r)（描述电子-电子之间的经典库仑相互作用）以及交换关联势V_{xc}(r)，即V_{eff}(r)=V_{ext}(r)+V_{H}(r)+V_{xc}(r)。对于这个虚构的无相互作用电子体系，每个电子的运动可以用一个单电子波函数\psi_{i}(r)来描述，满足以下的Kohn-Sham方程：\left(-\frac{\hbar^{2}}{2m}\nabla^{2}+V_{eff}(r)\right)\psi_{i}(r)=\epsilon_{i}\psi_{i}(r)其中，\frac{\hbar^{2}}{2m}\nabla^{2}是电子的动能算符，\epsilon_{i}是第i个电子的本征能量。体系的电子密度n(r)可以通过单电子波函数表示为：n(r)=\sum_{i=1}^{N}\vert\psi_{i}(r)\vert^{2}其中N是电子总数。通过求解Kohn-Sham方程，可以得到单电子波函数\psi_{i}(r)，进而计算出电子密度n(r)。在实际计算中，通常采用自洽迭代的方法来求解Kohn-Sham方程。首先，假设一个初始的电子密度n^{(0)}(r)，根据这个电子密度计算出有效势场V_{eff}^{(0)}(r)，然后将其代入Kohn-Sham方程求解得到一组新的单电子波函数\psi_{i}^{(1)}(r)，再根据这些波函数计算出新的电子密度n^{(1)}(r)。接着，用新的电子密度重新计算有效势场，再次代入Kohn-Sham方程求解，如此反复迭代，直到前后两次计算得到的电子密度或能量的变化小于某个设定的收敛标准，此时认为计算达到自洽，得到的电子密度和波函数即为体系的基态解。Kohn-Sham方程的物理意义在于，它将多电子体系中复杂的电子-电子相互作用通过有效势场的方式进行了近似处理，将多电子问题转化为一系列单电子问题。每个单电子在这个有效势场中独立运动，使得计算过程得到了极大的简化。通过求解Kohn-Sham方程，我们能够得到体系的电子结构信息，如电子密度分布、能级结构等，这些信息对于理解材料的物理化学性质至关重要。在研究材料的导电性时，电子密度和能级结构的信息可以帮助我们判断材料中电子的移动能力和能带结构，从而确定材料是导体、半导体还是绝缘体。2.3.3局域密度近似(LDA)与广义梯度近似(GGA)在密度泛函理论的实际计算中，交换关联泛函E_{xc}[n]的准确描述是一个关键而又极具挑战性的问题。由于电子之间的交换关联作用非常复杂，目前尚无法精确求解，因此通常采用近似方法来处理。局域密度近似（LDA）和广义梯度近似（GGA）是两种常用的近似方法，它们在描述电子相互作用时各有特点和适用范围。LDA是一种较为简单的近似方法，它将电子密度在空间中各点的变化近似为均匀。具体来说，LDA假设在空间某一点r处的交换关联能密度\epsilon_{xc}(n(r))只与该点的电子密度n(r)有关，且等于均匀电子气在相同密度下的交换关联能密度。基于这一假设，体系的交换关联能E_{xc}^{LDA}可以表示为：E_{xc}^{LDA}=\int\epsilon_{xc}(n(r))n(r)drLDA的优点是计算相对简单，计算量较小，在处理一些电子密度变化较为平缓的体系时，能够给出较为合理的结果。在研究简单金属的电子结构时，LDA能够较好地描述电子之间的相互作用，计算得到的能带结构和态密度等与实验结果有一定的吻合度。LDA也存在明显的局限性。由于它完全忽略了电子密度的梯度变化对交换关联能的影响，当体系中电子密度变化较大时，如在分子和固体表面等区域，LDA的计算结果往往会出现较大偏差。在计算分子的键长和键能时，LDA通常会高估分子的稳定性，导致计算得到的键长比实验值短，键能比实验值大。GGA则考虑了电子密度在空间中的梯度变化，相对于LDA更为精确。GGA认为交换关联能不仅与电子密度n(r)有关，还与电子密度的梯度\nablan(r)有关。其交换关联能E_{xc}^{GGA}的表达式通常可以写成：E_{xc}^{GGA}=\int\epsilon_{xc}(n(r),\nablan(r))n(r)dr通过引入电子密度的梯度项，GGA能够更好地描述电子之间的非局域相互作用，在处理分子、半导体和过渡金属等体系时表现出更好的性能。在计算半导体的能带结构时，GGA能够更准确地预测能带间隙，与实验值更为接近；在研究过渡金属的磁性时，GGA也能更合理地描述电子的自旋-轨道耦合等相互作用，得到更符合实际的磁性性质。GGA的计算相对复杂，计算量比LDA大。在处理大规模体系时，计算时间和计算资源的消耗会显著增加。而且GGA在某些情况下也存在一定的局限性，对于一些具有强电子关联效应的体系，GGA的计算结果仍然不够准确。2.3.4Hubbard参数U的引入(DFT+U计算)在一些体系中，如过渡金属氧化物和稀土化合物等，电子之间存在着强关联效应。传统的密度泛函理论（DFT），无论是采用局域密度近似（LDA）还是广义梯度近似（GGA），在处理这类体系时都存在一定的局限性。这是因为LDA和GGA主要适用于弱关联体系，它们对于电子之间复杂的强相互作用描述不够准确。在过渡金属氧化物中，过渡金属离子的d电子具有较强的局域性，电子之间的库仑相互作用和交换关联作用非常显著，传统DFT方法难以准确描述这些相互作用，导致计算得到的能带结构、磁性等性质与实验结果存在较大偏差。为了更准确地描述强关联体系，引入Hubbard参数U的DFT+U计算方法应运而生。Hubbard参数U主要用于描述原子内电子之间的强库仑相互作用。在DFT+U方法中，通过在哈密顿量中添加一个与U相关的项，来修正电子之间的相互作用。具体来说，对于具有强关联的原子轨道（如过渡金属的d轨道或稀土元素的f轨道），DFT+U方法将电子的能量分为两部分：一部分是传统DFT计算得到的能量，另一部分是由Hubbard项描述的强关联能量。Hubbard项的形式通常为：E_{U}=\frac{1}{2}\sum_{i\neqj}\sum_{\sigma,\sigma'}U_{ij\sigma\sigma'}n_{i\sigma}n_{j\sigma'}其中，U_{ij\sigma\sigma'}是Hubbard参数，表示不同轨道i和j上自旋为\sigma和\sigma'的电子之间的库仑相互作用强度，n_{i\sigma}和n_{j\sigma'}分别是轨道i和j上自旋为\sigma和\sigma'的电子占据数。引入Hubbard参数U后，DFT+U方法能够更准确地描述强关联体系中电子的行为。在计算过渡金属氧化物的能带结构时，DFT+U方法可以有效地修正传统DFT方法对能带间隙的低估问题，使得计算结果与实验值更加接近。在研究磁性材料的磁性时，DFT+U方法能够更好地描述电子的自旋-轨道耦合和磁相互作用，从而得到更准确的磁性性质。DFT+U方法也存在一些问题，其中最主要的是Hubbard参数U的取值具有一定的经验性。目前，U值的确定通常需要参考实验数据或通过一些理论模型进行估算，不同的取值方法可能会导致计算结果存在一定的差异。如何准确地确定Hubbard参数U，仍然是DFT+U方法研究中的一个重要课题。2.3.5赝势和PAW(ProjectorAugmentedWave)方法在第一性原理计算中，由于电子与原子核之间的相互作用非常复杂，特别是内层电子与原子核的紧密束缚，使得直接求解包含所有电子的多体问题计算量巨大且困难重重。为了简化计算并提高计算效率，赝势和投影缀加波（PAW）方法被广泛应用。赝势方法的基本原理是将原子核和内层电子看作一个整体，即离子实，用一个相对简单的赝势来代替离子实与价电子之间的复杂相互作用。在真实的原子体系中，内层电子紧紧围绕原子核，其波函数在原子核附近会发生剧烈振荡，这使得计算变得极为复杂。赝势通过构建一个平滑的、相对简单的势函数，使得价电子在这个赝势场中的运动能够近似地反映真实体系中价电子的行为。这个赝势的构建需要满足一定的条件，如在远离原子核的区域，赝波函数与真实波函数具有相同的渐近行为；赝势能够正确地描述价电子的能量和波函数等。通过使用赝势，我们可以将计算重点放在价电子上，大大减少了计算量，提高了计算效率。在计算固体材料的电子结构时，使用赝势可以快速得到价电子的能量和波函数，进而计算出材料的能带结构和态密度等性质。PAW方法是在赝势方法的基础上发展起来的一种更精确的方法。PAW方法的核心思想是在全电子波函数和赝波函数之间建立一种映射关系。它定义了一组投影算子，通过这些投影算子将全电子波函数和赝波函数联系起来。在PAW方法中，体系的总能量和其他物理量的计算可以在赝波函数空间中进行，同时通过投影算子的作用，能够准确地考虑全电子波函数的信息。与传统的赝势方法相比，PAW方法在保持计算效率的同时，能够更准确地描述电子的行为，特别是在处理原子内层电子的贡献时具有明显的优势。在研究过渡金属和稀土元素等具有复杂电子结构的体系时，PAW方法能够更精确地计算电子的能量和波函数，得到更准确的电子结构信息。PAW方法还具有很好的可转移性，即对于不同的化学环境和原子构型，PAW方法能够保持较好的计算精度。这使得PAW方法在处理各种材料体系时都具有较高的可靠性和适用性，成为现代第一性原理计算中广泛使用的方法之一。在计算不同晶体结构的过渡金属化合物时，PAW方法都能够准确地描述电子结构和化学键的性质，为研究材料的物理化学性质提供了有力的工具。2.4第一性原理计算的求解过程基于密度泛函理论的第一性原理计算的求解过程是一个复杂且严谨的过程，通常包含多个关键步骤。以研究Sm-CeO₂(111)表面CO分子催化氧化机理为例，首先需要构建合理的晶体结构模型。对于Sm-CeO₂(111)体系，需要确定CeO₂的晶体结构参数，如晶格常数等。CeO₂具有立方萤石结构，其晶格常数约为0.541nm。在此基础上，考虑Sm的掺杂情况，确定Sm的掺杂位置和掺杂量。假设掺杂量为x（x表示Sm原子在Ce原子位置的替代比例），通过构建超晶胞模型来模拟Sm掺杂的情况。超晶胞的大小和对称性会影响计算结果的准确性和计算量，一般根据研究目的和计算资源选择合适的超晶胞，如2×2×1的超晶胞，这样可以在保证一定计算精度的同时，控制计算量在可接受范围内。还需要考虑氧空位的影响，在模型中引入氧空位，通过移除特定位置的O原子来实现。确定氧空位的浓度和位置，例如在CeO₂(111)表面的特定原子层和原子位置上引入氧空位，以研究氧空位对CO分子催化氧化的影响。模型构建完成后，进行结构优化。结构优化的目的是使体系的总能量达到最小，从而得到最稳定的原子结构。在优化过程中，主要调整原子的位置和晶胞的形状。使用共轭梯度法、BFGS（Broyden-Fletcher-Goldfarb-Shanno）算法等优化算法。共轭梯度法通过迭代搜索能量梯度的反方向来逐步调整原子位置，使体系能量降低；BFGS算法则是一种拟牛顿算法，通过近似海森矩阵来加速收敛，更快地找到能量最小值。在优化过程中，设置收敛标准，如能量收敛精度为10⁻⁶eV/atom，力收敛精度为0.01eV/Å。当体系的能量和原子受力满足这些收敛标准时，认为结构优化完成。通过结构优化，可以得到Sm-CeO₂(111)表面原子的稳定排列方式，为后续研究CO和O₂分子的吸附以及催化氧化反应提供准确的结构基础。在完成结构优化后，进行电子结构计算。这一步骤的核心是求解Kohn-Sham方程，以获得体系的电子密度分布和能级结构。在计算过程中，需要选择合适的交换关联泛函。常用的交换关联泛函有局域密度近似（LDA）和广义梯度近似（GGA）。LDA适用于电子密度变化较为平缓的体系，计算相对简单，但对于一些电子密度变化较大的体系，如分子和固体表面等，计算结果可能存在偏差。GGA考虑了电子密度的梯度变化，在处理分子、半导体和过渡金属等体系时表现更优，但计算相对复杂，计算量比LDA大。对于Sm-CeO₂(111)体系，由于涉及到过渡金属Ce和稀土金属Sm，电子结构较为复杂，一般选择GGA中的PBE（Perdew-Burke-Ernzerhof）泛函，它在描述这类体系时具有较好的准确性。确定平面波基组的截断能量，截断能量的大小会影响计算精度和计算量。截断能量过低，计算精度会受到影响；截断能量过高，计算量会大幅增加。通常通过测试不同截断能量下的计算结果，如体系总能量、电子密度分布等，来选择合适的截断能量。对于Sm-CeO₂(111)体系，经过测试，选择截断能量为400eV，在此截断能量下，计算结果能够较好地收敛，且计算量在可接受范围内。通过求解Kohn-Sham方程，可以得到体系的电子密度分布、能级结构等信息，这些信息对于理解Sm-CeO₂(111)表面的电子性质和化学反应活性至关重要。2.5本课题选用的软件和硬件条件2.5.1VASP软件包简介VASP（ViennaAb-initioSimulationPackage）是一款在材料计算领域具有广泛应用和卓越影响力的软件包。它基于密度泛函理论（DFT），能够高效地进行电子结构计算和量子力学-分子动力学模拟。VASP在处理材料体系时，展现出诸多显著的功能和特点。在功能方面，VASP能够精确计算材料的电子结构，包括能级、电荷密度分布、能带、电子态密度等关键信息。通过计算能带结构，可以清晰地了解材料中电子的能量分布情况，判断材料是导体、半导体还是绝缘体。在研究半导体材料时，VASP计算得到的能带间隙数据对于理解材料的电学性能和光电器件的应用具有重要指导意义。VASP还可以计算材料的光学性质，如介电函数、吸收光谱等。这些光学性质的计算结果有助于研究材料在光电子学领域的应用，如发光二极管、光电探测器等器件的设计和优化。在计算材料的磁学性质方面，VASP能够分析材料的磁性起源、磁矩大小以及磁相互作用等。对于磁性材料的研究，这些计算结果对于开发高性能的磁性存储器件和磁性传感器至关重要。从特点来看，VASP具有高精度和高效率的优势。它采用平面波基组和赝势或投影缀加波（PAW）方法对多电子系统进行描述，能够准确地处理电子与离子间的相互作用。PAW方法在保持计算精度的同时，有效提高了计算效率，使得VASP能够处理较大规模的材料体系。在研究复杂的合金材料时，PAW方法能够准确描述合金中不同原子的电子结构和相互作用，为研究合金的性能提供可靠的计算结果。VASP在迭代求解过程中采用了Broyden和Pulay密度混合方案，加速了自洽循环的收敛速度。这使得VASP在进行大规模计算时，能够更快地达到收敛状态，节省计算时间和资源。VASP在材料计算领域的应用极为广泛。在凝聚态物理研究中，它被用于探索新型量子材料的电子结构和物理性质，如拓扑绝缘体、高温超导体等。通过VASP的计算，研究人员能够深入了解这些材料的奇特量子特性，为新材料的开发和应用提供理论支持。在材料科学领域，VASP可用于设计和优化新型材料，如高强度合金、高性能催化剂等。通过计算材料的结构参数、力学性质和电子结构等，研究人员可以预测材料的性能，并指导实验合成具有理想性能的材料。在化学领域，VASP可用于研究化学反应机理、分子结构和光谱性质等。在研究有机化学反应时，VASP能够计算反应过程中的能量变化和电子转移情况，帮助研究人员理解反应机理，优化反应条件。2.5.2MaterialsStudio软件简介MaterialsStudio是一款功能强大且全面的材料模拟软件，在材料研究领域发挥着重要作用，尤其是在构建模型和可视化分析方面展现出独特的优势。在构建模型方面，MaterialsStudio提供了丰富多样的工具和功能，能够满足不同类型材料模型的构建需求。对于晶体材料，它可以根据晶体的空间群信息和原子坐标，快速准确地构建出晶体结构模型。在构建CeO₂晶体模型时，只需输入CeO₂的立方萤石结构的空间群信息以及Ce和O原子的坐标，MaterialsStudio就能自动生成精确的CeO₂晶体结构。通过调整原子的位置、晶格常数等参数，可以进一步优化晶体结构，以满足不同研究目的。对于表面模型，MaterialsStudio能够方便地创建各种晶体表面，如CeO₂(111)表面。通过对晶体结构进行切割和表面弛豫处理，可以得到稳定的表面模型，用于研究表面吸附、催化反应等过程。还可以构建复合材料模型，将不同材料的组分进行组合，模拟复合材料的结构和性能。在构建CeO₂与金属纳米颗粒复合的催化剂模型时，MaterialsStudio可以精确地控制金属纳米颗粒的尺寸、形状和在CeO₂表面的分布，为研究复合材料的协同催化作用提供基础。在可视化分析方面，MaterialsStudio具有直观、便捷的可视化界面，能够将复杂的材料结构和计算结果以直观的图形方式展示出来。它可以清晰地展示晶体结构的三维模型，通过旋转、缩放等操作，研究人员可以从不同角度观察晶体结构，了解原子的排列方式和相互作用。在分析CeO₂晶体结构时，通过可视化界面可以直观地看到Ce原子和O原子的配位情况，以及晶体中的晶格缺陷和氧空位等信息。MaterialsStudio还能够对电子结构计算结果进行可视化分析，如电荷密度分布、能带结构、态密度等。将电荷密度分布以等值面的形式展示出来，可以直观地了解电子在材料中的分布情况，分析化学键的形成和电子转移过程。通过绘制能带结构和态密度图，可以直观地观察到材料的电子能量分布和态密度变化，为研究材料的电学和光学性质提供直观的依据。MaterialsStudio还支持对分子动力学模拟结果进行可视化分析，展示分子的运动轨迹和相互作用过程，帮助研究人员深入理解材料的动态行为。2.5.3课题组的硬件条件本课题组配备了高性能的计算服务器，其硬件配置对于本课题的第一性原理计算研究提供了坚实的支撑。服务器采用了多核心的中央处理器（CPU），例如IntelXeonPlatinum系列处理器，该系列处理器具有较高的核心数和主频，能够并行处理大量的计算任务。在进行Sm-CeO₂(111)表面结构优化和电子结构计算时，多核心CPU可以同时计算不同原子的受力和电子波函数，大大提高了计算效率，缩短了计算时间。服务器还配备了大容量的内存，一般为128GB甚至更高。充足的内存能够保证在计算过程中，大量的中间数据和计算结果能够及时存储和读取，避免了因内存不足导致的计算中断或性能下降。在处理大规模的超晶胞模型时，大容量内存可以存储更多的原子坐标和电子密度信息，确保计算的顺利进行。为了进一步提高计算速度，服务器还配备了高性能的图形处理器（GPU），如NVIDIATesla系列GPU。GPU具有强大的并行计算能力，在进行复杂的矩阵运算和电子结构计算时，能够显著加速计算过程。在求解Kohn-Sham方程时，GPU可以快速计算哈密顿矩阵的本征值和本征向量，从而加速电子结构的收敛。课题组还拥有高速的存储设备，采用了固态硬盘（SSD）阵列，提供了快速的数据读写速度。在存储和读取大量的计算数据和结果文件时，SSD阵列能够大大缩短数据传输时间，提高工作效率。在保存和读取Sm-CeO₂(111)表面结构优化和电子结构计算的结果文件时，高速存储设备可以快速完成操作，方便研究人员及时分析和处理数据。这些硬件条件的合理配置，使得本课题组在进行Sm-CeO₂(111)表面CO分子催化氧化机理的第一性原理计算研究时，能够高效、准确地完成各项计算任务，为深入研究提供了有力的保障。2.6本章小结本章深入阐述了第一性原理的理论基础及计算方法，这是研究Sm-CeO₂(111)表面CO分子催化氧化机理的核心理论工具。绝热近似通过将电子和原子核的运动分离，极大地简化了多粒子体系问题，使我们能够专注于电子的运动，为后续研究奠定了基础。单电子近似进一步将多电子问题转化为单电子问题，通过假设每个电子在平均势场中独立运动，为求解多电子体系的薛定谔方程提供了可行的途径。密度泛函理论是第一性原理计算的核心，Hohenberg-Kohn定理确立了基态电子密度与体系性质的紧密联系，为理论计算提供了坚实的基础。Kohn-Sham方程作为DFT的关键方程，通过引入有效势场，将复杂的多电子相互作用简化为单电子问题求解，使得计算过程更加高效和可行。在实际计算中，局域密度近似（LDA）和广义梯度近似（GGA）等方法用于处理交换关联泛函，它们各有优劣，适用于不同的体系。对于电子密度变化较为平缓的体系，LDA计算简单且能给出合理结果；而对于电子密度变化较大的体系，GGA考虑了电子密度的梯度变化，计算结果更为准确，但计算量相对较大。在处理具有强电子关联效应的体系时，引入Hubbard参数U的DFT+U方法能够更准确地描述电子之间的相互作用，弥补了传统DFT方法的不足。赝势和PAW方法则是为了简化电子与原子核相互作用的计算，提高计算效率。赝势方法将原子核和内层电子看作离子实，用简单的赝势代替其与价电子的复杂相互作用；PAW方法在赝势的基础上，通过建立全电子波函数和赝波函数的映射关系，在保持计算效率的同时，更准确地描述了电子行为，尤其适用于处理具有复杂电子结构的体系。第一性原理计算的求解过程包括构建晶体结构模型、结构优化和电子结构计算等关键步骤。构建模型时，需综合考虑晶体结构参数、原子掺杂和缺陷等因素，以确保模型的准确性和代表性。结构优化通过调整原子位置和晶胞形状，使体系总能量最小，得到稳定的原子结构。电子结构计算则通过求解Kohn-Sham方程，获得体系的电子密度分布和能级结构，为深入理解材料的物理化学性质提供关键信息。本课题选用VASP软件包进行计算，它基于密度泛函理论，具备强大的功能，能够高效计算材料的电子结构和性质，在材料模拟和计算领域应用广泛。MaterialsStudio软件则用于构建模型和可视化分析，为研究提供了便捷的工具。课题组配备的高性能计算服务器，拥有多核心CPU、大容量内存、高性能GPU和高速存储设备等硬件设施，为第一性原理计算提供了有力的硬件支持，确保计算任务能够高效、准确地完成。这些理论基础、计算方法、软件工具和硬件条件的有机结合，为深入研究Sm-CeO₂(111)表面CO分子催化氧化机理提供了坚实的保障。三、CeO₂模型建立、计算及验证3.1CeO₂(111)表面计算模型在构建CeO₂(111)表面的原子模型时，首先需要深入了解CeO₂的晶体结构特征。CeO₂具有立方萤石型结构，其空间群为Fm-3m，晶格常数约为0.541nm。在这种结构中，Ce原子位于面心立方晶格的顶点和面心位置，每个Ce原子被8个O原子以立方体的形式配位，形成Ce-O₈多面体；而O原子则填充在四面体空隙中，每个O原子与4个Ce原子相连。为了准确模拟CeO₂(111)表面的性质，我们采用了超晶胞模型。超晶胞模型的选取需要综合考虑计算精度和计算资源的限制。经过一系列的测试和分析，我们选择了2×2×3的超晶胞结构。这种结构包含了16个Ce原子和32个O原子，在保证计算精度的同时，能够有效控制计算量，使计算过程在合理的时间和资源范围内完成。在确定超晶胞结构后，需要对其进行表面构建和优化。我们沿着[111]方向对超晶胞进行切割，以获得CeO₂(111)表面。在切割过程中，需要考虑表面原子的弛豫现象。由于表面原子的配位环境与体内原子不同，表面原子会发生一定程度的位移，以降低表面能，达到更稳定的状态。为了模拟这种弛豫现象，我们在计算中允许表面原子进行充分的弛豫，而固定体内原子的位置。通过多次结构优化计算，使用共轭梯度法等优化算法，不断调整表面原子的位置，直到体系的总能量和原子受力达到收敛标准。能量收敛精度设置为10⁻⁶eV/atom，力收敛精度设置为0.01eV/Å。经过优化后，得到了稳定的CeO₂(111)表面结构。在优化后的CeO₂(111)表面结构中，我们对表面原子的几何参数进行了详细分析。通过计算发现，表面Ce-O键长与体内Ce-O键长存在一定差异。表面Ce-O键长略微伸长，这是由于表面原子的配位不饱和，受到的束缚力相对较小，导致键长伸长。具体来说，体内Ce-O键长约为0.234nm，而表面Ce-O键长在优化后约为0.238nm。这种键长的变化对表面的电子结构和化学活性产生了重要影响。键长的伸长使得表面Ce原子的电子云分布发生改变，表面的电子密度降低，从而增强了表面对反应物分子的吸附能力，为后续的CO和O₂分子吸附以及催化氧化反应提供了更有利的条件。我们还对表面的原子排列和对称性进行了分析。CeO₂(111)表面呈现出六边形的原子排列结构，具有良好的对称性。这种原子排列和对称性对表面的物理化学性质具有重要影响，例如影响反应物分子在表面的吸附位点和吸附构型，进而影响催化反应的活性和选择性。3.2计算方法和相关公式在对Sm-CeO₂(111)表面CO分子催化氧化机理的研究中，采用基于密度泛函理论（DFT）的第一性原理计算方法，运用VASP软件包进行具体计算。在计算过程中，选用广义梯度近似（GGA）中的PBE泛函来描述电子之间的交换关联相互作用。相较于局域密度近似（LDA），PBE泛函考虑了电子密度的梯度变化，在处理具有复杂电子结构的体系时表现更为出色，能够更准确地描述Sm-CeO₂体系中电子之间的相互作用。对于Sm-CeO₂体系中Ce和Sm原子的内层电子，采用投影缀加波（PAW）方法进行处理。PAW方法能够在保证计算精度的同时，有效提高计算效率。它通过定义一组投影算子，将全电子波函数和赝波函数联系起来，在赝波函数空间中进行计算的同时，准确考虑全电子波函数的信息，尤其适用于处理具有复杂电子结构的体系，如包含过渡金属和稀土金属的Sm-CeO₂体系。平面波基组被用于描述电子的波函数。平面波基组具有数学形式简单、易于处理的优点，能够在倒易空间中对电子的波函数进行精确展开。在计算中，设置平面波截断能为400eV。截断能的选择至关重要，它直接影响计算的精度和效率。截断能过低，会导致平面波基组对电子波函数的描述不够精确，从而影响计算结果的准确性；截断能过高，则会大幅增加计算量，耗费更多的计算资源和时间。通过多次测试不同截断能下体系的总能量、电子密度分布等关键物理量，发现当截断能为400eV时，计算结果能够较好地收敛，且计算量在可接受范围内，因此选择400eV作为平面波截断能。为了模拟真实的表面体系，采用