铍团簇吸附一氧化碳分子的密度泛函理论研究.pdf

第 卷第 期 年 月西南科技大学学报 ：收稿日期：；修回日期：基金项目：四川省自然科学基金项目（）作者简介：第一作者，林莉（），女，硕士研究生；通信作者，邝向军（），男，博士，教授，研究方向为材料结构和性质的设计，：铍团簇吸附一氧化碳分子的密度泛函理论研究林莉邝向军（西南科技大学数理学院四川绵阳 ）摘要：采用密度泛函理论和广义梯度近似方法，在 水平上对 （）团簇进行了全电子计算。结果表明：分子倾向于平行吸附在纯 团簇的表面并占据铍团簇的外部位置，吸附后主团簇结构基本没有发生较大的突变，但是 分子的吸附会使 团簇的对称性下降，时铍团簇对 分子的吸附相对较强，时铍团簇对 分子的吸附相对较弱；吸附 分子后，多数 键长以及 的键长变长，即一氧化碳的吸附会使 以及 的相互作用减弱；吸附之后 团簇的能隙、垂直电离势、垂直亲和势均发生了较为明显的变化，的能隙出现了较为明显的奇偶振荡并在 出现了最大值，吸附前后团簇的垂直电离势、垂直亲和势曲线变化大趋势基本相同，随着原子数的增加团簇的垂直电离势逐渐降低而垂直亲和势逐渐升高。分子吸附后对纯 团簇的稳定性产生了较大影响，的键长、键长及团簇的能隙均发生了较大变化，吸附后仍存在幻数结构 。关键词：铍团簇一氧化碳分子吸附密度泛函理论中图分类号：文献标志码：文章编号：（），（，）：（），：；近年来，随着团簇学科的迅猛发展，研究者对团簇的各种物理和化学性质以及它们随着团簇分子大小的演变规律已经有了较为系统的认知。人们希望通过在团簇表面吸附不同的原子、分子来合成具有独特性质的新型团簇材料，这将对催化科学、表面科学、材料科学、凝聚态物理等领域产生深刻的影响。金属团簇作为团簇的重要组成部分，其在理论及实验上是人们关注的热点。二价碱土金属铍是一种稀有的金属材料，也是最轻的金属之一。在宏观上，金属铍是航天、航空、原子能及冶金工业必不可少的材料 。从微观上看，铍也有着一些有趣的特性，二聚体铍是弱结合，而较大的铍团簇表现为自由电子体系并显示金属性质 。等指出，小分子铍原子之间的相互作用类似于范德瓦耳斯相互作用，随着团簇尺寸的增加，铍原子之间的成键会转换成共价键，再进一步转变成金属键 。此外，等发现当 ，时铍团簇中出现了幻数簇 ，这与电子壳模型的推论是一致的 。铍原子除了本身的特殊性之外，其另一重要用途是以铍原子为主团簇掺杂其他金属，进而形成二元合金。葛桂贤 对铍团簇掺杂 ，等金属进行了研究，结果表明掺杂后团簇的结构及稳定性都会发生显著的变化；雷雪玲等 分析了 与 的基态结构、稳定性以及极化率，结果表明锂原子的掺杂增加了团簇的极化率。这些研究表明，金属原子的掺杂对铍团簇结构、能量、稳定性都产生了较大的影响。在过去的研究中，人们更多讨论了金属掺杂对铍团簇的影响，而对于铍团簇吸附非金属的研究较少。事实上，铍原子可以和许多非金属相互作用形成多元化合物，例如氧化铍（），尽管它具有较大的毒性，但它可以用于制造铍氧化物、陶瓷、玻璃，也可用作核反应器燃料和减速剂，还可用作有机反应催化剂 。李淑萍等 发现铍团簇表面存在大量悬挂键，能够更好吸附一氧化氮分子，进而减轻一氧化氮排放对空气的污染。董兰等 对铍团簇吸附氢气进行了研究，发现随着铍团簇原子数量增加，其与氢原子的杂化程度逐渐变大。尽管越来越多的研究者投入于对该领域的研究，但有关 在铍团簇表面的吸附行为的研究未见报道。小分子铍团簇不仅具有复杂的表面，且原子之间的相互作用较弱，这与其他金属有着明显的区别，因此研究铍团簇吸附一氧化碳的结构、能量、吸附位置、吸附强度具有重要意义。本文从第一性原理出发，在 （）水平上运用广义梯度近似方法对一氧化碳分子在铍团簇上的吸附进行计算，并对以下问题进行分析：（）一氧化碳的吸附是否会对纯 团簇的结构和稳定性产生影响；（）团簇的结构和稳定性是如何演化的；（）吸附后 分子以及 团簇的键长、能隙如何变化；（）吸附后团簇是否仍存在幻数结构。计算方法在计 算 过 程 中，采 用 中 的 模块，选择广义度近似（）方法进行拟合，选用电子的交换关联函数（）进行全电子计算，容限为 。为了加速计算，使用迭代子空间中的直接反演（）方法，并且将拖尾值设置为 （，下同）。在结构优化过程中，自旋不受约束，结构的对称性不受约束。最大力、能量和最大位移的收敛公差分别为 ，和 。在进行 团簇吸附一氧化碳分子的相关计算时，最关键的是确定最低能量结构。首先需要得到纯 团簇的最稳定结构。结合前人对纯 团簇的研究，发现密度泛函理论中的 在纯 团簇的结构和稳定性上表现得非常好，因此笔者采用相同的方法构建纯 团簇并对其进行几何优化，优化结构如图 所示，与前人的研究成果一致 。在此基础上，为了得到 团簇的初始结构，让铍团簇的每个可能的非等价位置与 分子成键，得到了 团簇的初始结构，包括 种可能的成键结构：原子与铍团簇成键、原子与铍团簇成键以第 期林莉，等：铍团簇吸附一氧化碳分子的密度泛函理论研究及 ，原子同时与铍团簇成键。之后计算了吸附前后团簇的总能量（）、平均结合能（）、振动频率、键长、最高已占据轨道能级、最低未占据轨道能级、能隙（）、垂直亲和势（）以及垂直电离势（）等参数。团簇的平均结合能 ，计算公式如下：（）（）（）（）（）（）（）（）式中：，是单个铍的能量；是 团簇的总能量；（）表示得到一个电子后团簇的总能量；（）表示失去一个电子后团簇的总能量。（）团 簇 的 平 均 结 合 能（），的计算公式如下：（）（）（）（）（）（）（）（）（）（）式中：表示 （）团簇的总能量；表示 （）团簇的吸附能；（）表示得到一个电子后团簇的总能量；（）表示失去一个电子后团簇的总能量。为了检测该方法的可靠性，比较了二聚体铍的相关数据。由表 可知，的键长为 ，这与文献 的实验结果（）及之前报道的理论数据（和 ）一致，表明我们选择的方法能够较准确研究 团簇对一氧化碳的吸附行为。结果与讨论 团簇的几何结构优化与电子性质采用密度泛函理论对 （）团簇的几何结构进行了优化，获得的最低能量结构如图 所示。（）团簇的对称性、总能 、平均原子结合能 、振动频率、最高已占据轨道能级 、最低未占据轨道能级 、能隙 、垂直电离势 以及垂直亲和势 如表 所示。表 结果显示，随着团簇 原子数的增加，团簇的总能量降低，平均结合能增加，团簇越发稳定。此外，我们知道 和 的差值 反应的是团簇的电子结构稳定性，小，电子就容易从 跃迁到 ，电子结构的稳定性就比较低，反之，则电子就不容易从 跃迁到 ，电子结构的稳定性就比较高。由表 可知，和 的 较高，电子不易跃迁，结构稳定性较高，而 和 的 小于 ，结构不稳定。（）的基态结构为了得到 （）团簇的最低能量结构，模拟了两种可能的情况：一种是 平行吸附在 团簇上；第二种是以 分子的 端或者 端吸附在 团簇上。在对各种可能的情况进行 结构优化及频率分析后，得到的 （）团簇的最低能量结构如图 所示，（）基态团簇的对称性、总能量 、吸附能 、振动频率 、键长、最高已占据轨道能级 、最低未占据轨道能级 、能隙 、垂直电离势 及垂直亲和势 如表 所示；图给出了 团簇及 团簇中 平均键长；图 为 团簇中 ，原子间距。的基态几何结构是具有 （）对称性的四边形，是体系中仅有的平面结构，相比于纯 团簇，吸附使得 的平均键长缩短了 ，而 键长则伸长了 ；的基态是具有（）对称性的六面体结构，这是 体系中的第一个三维结构，此时 的平均键长缩短了 ，而 键长则伸长了 ；仍是具有 （）对称性的空间结构。值得注意的是，从图可以看出，开始，吸附后的 平均键长图 （）团簇的基态结构 （）西南科技大学学报第 卷表 （）能级参数 （）团簇分子 对称性 总能量 平均结合能 振动频率 团簇分子 对称性 （）（）总能量 平均结合能 振动频率 表 （）团簇的能级参数 （）团簇分子对称性 （）（）（）（）（）（）（）（）（）（）（）（）都高于吸附前的 平均键长，此时 ，原子间距仍被拉长；的基态结构是在 的基础上添加一个铍原子得到的，具有 对称性，此时 键长变为 ，伸长了 ，而铍原子平均键长只缩短了 ；是正八面体结构，具有较高的对称性，吸附后对称性降为，键长也发生了较为明显的伸长，相比于 端吸附以及 端吸附，此时能量为 ，是几种可能的结构中最稳定的结构。是具有 对称性的空间结构，分子平行吸附于铍团簇的表面，吸附第 期林莉，等：铍团簇吸附一氧化碳分子的密度泛函理论研究注：表示 基态结构，表示 吸附 时的初始构型，表示 驰豫后的构型图 （）和 （）团簇的基态结构 （）（）图 团簇及 团簇中 的平均键长 图 基态团簇中 ，原子间距 后 键长伸长了 。团簇可以看作具有 对称性的双帽三角反棱柱，分子的吸附使团簇的对称性降为 （），同时减弱了 与 之间的相互作用。从 到 ，每种结构都可以由前一种结构加上额外的 原子产生。可以看作是在 的边上添加了一个 原子，此时 中 键长为 ，相比吸附前伸长了 ，说明吸附明显减弱了 分子间的相互作用。则可以看作是在 的顶端四面体底部添加了一个 原子，具有与 相同的 对称性。端四棱锥底部连接一个 原子形成 最低能西南科技大学学报第 卷量结构，仍具有 对称性，但是 键的伸长量显著小于相邻团簇的伸长量。以此类推，的最稳定结构可以通过将 原子附着到 的底部来获得，此时团簇对称性由原来的 降为，这表明 分子的吸附极大地降低了 团簇的对称性。的基态结构则是在 的中心添加一个 原子得到的，与 相同，此时的对称性也降到了。通过图 可以看出 ，的 键长相对较长，而 ，的 键长较短，这表明在 ，时 团簇对一氧化碳分子的吸附强度高于 ，时的吸附强度。由上述分析可得：（）分子倾向于平行吸附在 团簇的表面并占据铍团簇的外部位置，且吸附后主团簇结构没有发生较大突变。（）的吸附会使 团簇的对称性下降。（）分子的吸附会使 的平均键长、键长发生不同程度的伸长，被吸附后活性提高。能量性质和相对稳定性团簇的结合能是团簇热力学稳定性的度量。从表 和图 可以发现，所有 团簇的结合能都大于相应的纯 团簇的结合能，并且随着铍原子数的增加，纯铍团簇的结合能逐渐增加，在 时达到最大值 。同时，团簇的结合能逐渐减小，在 时达到最小值 ，然后在 附近振荡。这表明所有 团簇在能量上都比相应的铍团簇更稳定。从表 和图 中 团簇与一氧化碳之间的吸附曲线可以看出，吸附能在 ，时出现了峰值，这表明此时铍团簇和一氧化碳的相互作用较强，吸附后的结构较为稳定，且在 之后团簇的吸附能出现了持续下降，即 团簇相对于其他团簇更稳定；在 ，时出现了最小值，表明 ，团簇的稳定性较低。此外，垂直电离势（）常被用来研究小团簇的化学稳定性，越大，能级越深，反应活性越低，化学稳定性越高。能隙 是用于考查团簇电子稳定性的另一个重要参数。越大，从价带到导带激发电子所需的能量越高，对应的稳定性越高。图 分别给出了吸附前后团簇的 随团簇尺寸的变化规律。对于纯 团簇，其 的峰值出现在 ，这与已有的研究成果是一致的。与吸附前相比，吸附 后团簇的 出现了 种情况：，时，吸附 后团簇的吸附能显著提升，特别是 ，其 为 ，相比于吸附前提升了 ，表明这 个团簇具有高的热力学稳定性；，时，吸附后团簇的 显著减小，尤其是 ，此时 减小了 ，约为吸附前的一半，这说明 分子的吸附使得 团簇的化学活性得到了显著提升；其他情况下，吸附前后的 变化比较小。图 给出了吸附前后团簇的 随团簇尺寸的变化关系，当 ，时，的吸附显著提升了团簇的 ；此外，在 和 处出现的小峰值以及随后的下降表明这两个团簇比它们相邻的团簇更不容易失去电子。计算了 的垂直电子亲合势 ，结果如图所示。图表明，除了 分子的吸附极大地降低了 的 值之外，两条曲线的趋势几乎相同。此外，在 ，时 出现了两个极小值，这表明 和 团簇在获得电子方面的相对稳定性。综上所述，可得出 是吸附后体系中最稳定的结构。图 和 团簇结合能和吸附能的变化 图 ，（）基态团簇的 ，（）第 期林莉，等：铍团簇吸附一氧化碳分子的密度泛函理论研究图 ，（）基态团簇的 ，（）图 ，（）基态团簇的 ，（）结论采用密度泛函理论和广义梯度近似方法在 水平上对 （）团簇进行了全电子计算，结论如下：（）分子倾向于平行吸附在纯 团簇的表面并占据铍团簇的外部位置，吸附后主团簇结构没有发生较大的突变，但是 的吸附会使 团簇的对称性下降。（）分子的吸附会使 的平均键长、键长发生不同程度的伸长，的化学活性显著提升。（）的结合能高于纯 团簇，即 能量稳定性高于纯 团簇，在 ，时 团簇的吸附能最大。当 ，时 团簇的能隙相比于纯 团簇有显著提升；，和 时，的吸附显著提升了团簇的垂直电离势；，时出现了垂直亲和势的两个极小值。是吸附后该体系最稳定的结构。参考文献 ，：，：，（）：?，（）：张文献，刘 磊，李郁芬 铍团簇 （）的基态能量、结构及其拓扑性质 原子与分子物理学报，（），：，：，：，（）：，（）：，（）：葛桂贤，井群，杨致，等第一性原理对 （）团簇最低能量结构及其电子性质的研究 物理学报，（）：，（），（）：葛桂贤过渡金属（，）掺杂纯铍团簇的第一性原理研究 郑州：河南大学，（下转第 页）西南科技大学学报第 卷参考文献 方桂花，孙鹏博，于孟欢，等石蜡相变材料热性能提升研究进展 应用化工，（）：，李国俭相变储能材料开发与封装技术研究进展 热力发电，（）：黄港，邱玮，黄伟颖，等 相变储能材料的研究与发展 材料科学与工艺，（）：林文珠，凌子夜，方晓明，等 相变储热的传热强化技术研究进展 化工进展，（）：贺斌，何光进，孙彩云，等相变储能材料研究进展及应用 信息记录材料，（）：王若钰，梁斌，朱英明，等石蜡基高热导率相变储能材料的制备 化学试剂，（）：陈璞，蒋达华，徐玉珍，等 相变材料选用及储能墙体节能研究进展 化工新型材料，（）：史琛，王平，杨柳建筑用石蜡类相变储能材料的改性研究进展 中国材料进展，（）：张静，何伟，王浩，等石蜡基新型定形相变蓄热材料的热 物 性 分 析 建 筑 热 能 通 风 空 调，（）：，王宇，李琳，马静 石墨基定形相变材料研究进展 化工新型材料，（）：陈雪芳，门明燃，李强，等石蜡 聚丙烯定型相变储能材料的泄漏率和定型能力的研究 化工新型材料，（）：，（）：，：，：，：，（）：徐众，侯静，李军，等膨胀石墨 有机质复合相变材料的制备及性能 化工进展，（）：闫嘉森，韩现英，党兆涵，等石蜡 膨胀石墨 石墨烯复合相变储热材料的制备及性能 高等学校化学学报，（）：翟天尧，李廷贤，仵斯，等 高导热膨胀石墨 硬脂酸定形相变储能复合材料的制备及储 放热特性 科学通报，（）：邓杨冲，陈华，金妍 膨胀石墨 石蜡复合相变材料的制备与物性研究 制冷技术，（）：张佳利，丁宇，曲丽洁，等石蜡 膨胀石墨复合相变储热单元的放热性能 储能科学与技术，（）：杨永梅，葛凤华，王松 低质水相变墙体的性能模拟 暖通空调，（）：檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪 （上接第 页），：，（）：，（，）（），李淑萍，孟江，王继刚 在金属 （）团簇表面的平行吸附 原子与分子物理学报，（）：董兰，蒋树斌，郑申声，等（）团簇吸附的密度泛函研究 原子能科学技术，（）：，?，：，（）：，（）：，（）：，：，（）：第 期杨琨，等：石蜡 膨胀石墨复合定形相变材料的制备及季节适应性研究