🔬 materials science

Temperature-Induced Reorganization of Supported Zn $_3$ Clusters on Cu(111): From Minimum-Energy Structures to Finite-Temperature Ensembles

本研究表明，有限温度下的熵效应从根本上改变了负载型 Zn $_3$ 集群在 Cu(111) 上的优选结构，使其从紧凑的三角形形状转变为延伸的线性构型，从而对仅基于静态 0 K 最低能量模型来识别催化活性位点的可靠性提出了挑战。

原作者： Jiayan Xu, Zheng Yu, Abhirup Patra, Amar Deep Pathak, Sharan Shetty, Detlef Hohl, Roberto Car

发布于 2026-08-19

📖 1 分钟阅读☕ 轻松阅读

原作者： Jiayan Xu, Zheng Yu, Abhirup Patra, Amar Deep Pathak, Sharan Shetty, Detlef Hohl, Roberto Car

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在工业化学领域，催化剂是让反应更快、更高效地发生的幕后英雄，且其自身并不会被消耗。可以将催化剂想象成一个专门的工作台，原材料在这里相遇、分解并重新组装成新的物质。几十年来，科学家们一直试图从原子层面理解这些工作台的具体形态，希望如果他们能看到原子的精确排列方式，就能设计出更好的催化剂。一个普遍的观点是，最稳定的原子排列方式——即维持结构所需能量最低的那种——才是执行工作的结构。这个想法在寒冷且静止的状态下运作良好，但在现实世界的化学反应中，温度很高，原子一直在不断地跳动和移动。长期以来困扰研究人员的问题是：在开启加热模式后，催化剂结构的静态、冰冷图像是否足以解释其行为。

一个研究小组通过研究一种用于制造甲醇（一种至关重要的燃料和化学构建模块）的特定类型催台剂来解决这一问题。这种催化剂是铜、氧化锌和氧化铝的混合物。虽然铜表面已被充分理解，但位于其上的微小锌簇所扮演的角色一直是一个谜。长期以来，科学家们假设这些锌簇会为了实现能量效率最大化而收缩成紧凑的三角形形状。然而，研究人员怀疑这种寒冷、静态的视角可能会产生误导。为了测试这一点，他们利用由机器学习驱动的高级计算机模拟，观察这些微小的锌簇在温度从冰点上升到工业反应器中的高温环境时的行为。

该研究重点关注了通常被氧和氢包围的三个锌原子的微小群体，它们坐落在铜表面上。研究人员首先证实了大家的预期：在零度时，这些簇倾向于聚拢在一起，形成一个紧凑的三角形。但当他们模拟高温环境时，故事发生了变化。随着温度升高，这些簇开始拉伸。它们不再保持紧凑的三角形，而是展开成了长长的直线。这种转变并非由原子寻找更低能量状态所驱动；事实上，拉伸后的形状实际上需要稍高的能量来维持。这种变化的原因完全是另一回事：熵。简单来说，拉伸后的形状为原子提供了更多的摆动和移动的空间。在高温下，这种运动的自由度变得比维持紧凑三角形所节省的能量更加珍贵，从而导致这些簇重新组织成这些线形结构。

研究人员还观察了簇下方的表面如何影响这种行为。他们测试了纯铜表面、仅混入单个锌原子的表面，以及顶层有四分之一为锌的表面。在纯铜和单锌表面上，这些簇可以自由移动，像桌上的弹珠一样在表面上滑动。这种流动性是一个问题，因为如果簇移动过多，它们会相互碰撞并合并成较大的、失去活性的团块，这一过程被称为烧结。然而，在锌浓度较高的表面上，这些簇停止了移动。表面的锌原子起到了锚定的作用，将这些簇锁定在原地。这一发现表明，虽然簇可能会因热量而改变形状，但表面的成分决定了它们是留在原地还是四处游荡。

或许最令人惊讶的发现是，特定化学中间体（甲酸盐分子）的存在并未改变游戏的基本规则。即使连接着这样一个额外的分子，这些簇在温度升高时仍然会从三角形转变为线形，且富锌表面依然能阻止它们移动。研究表明，仅观察低温下最低能量结构的传统方法忽略了活性位点的真实本质。在实际运行条件下，催化剂不是一个单一、刚性的形状，而是一个动态的集合体，在紧凑构型与拉伸构型之间不断切换。

这些结果挑战了科学家建模催化剂的标准方式。多年来，假设一直是零温下最稳定的结构就是最重要的结构。本文证明了这一假设是有缺陷的。活性位点并非静态的雕塑，而是随温度变化的流体群体。研究人员发现，维持形状所需的能量与通过移动获得的熵之间的竞争，决定了其结构。通过使用机器学习来加速模拟，他们能够观察到这些由于温度驱动而产生的变化，而这在以前较慢的旧方法中是无法捕捉到的。

这些发现对于设计更好的催化剂具有重要意义。如果工程师想要创造一种寿命更长、性能更好的催化剂，他们不能仅仅寻找在寒冷实验室中最稳定的原子排列。他们必须考虑原子在反应器高温且充满移动分子的环境下会如何表现。研究表明，在铜表面添加更多的锌是一个关键策略，这不仅是为了改变簇的形状，更是为了将它们钉在原地，防止它们因游荡而失去活性。这项工作提供了一个更清晰、更真实的图景，展示了催化剂内部究竟发生了什么，将该领域从对静态快照的观察转向了对原子生命动态电影的记录。

技术摘要：温度诱导的负载型 Cu(111) 上 Zn₃ 集群的重组

问题陈述
在多相催化领域，特别是用于甲醇合成的工业级铜/氧化锌/氧化铝（CZA）催化剂中，活性位点的原子级结构及其稳定性仍存在争议。虽然研究认为在 Cu 界面处存在微小的 Zn 基物种作为活性位点，但对其识别高度依赖于零温（0 K）密度泛函理论（DFT）优化。这种方法假设 0 K 下的最小能量构型代表了反应条件下的活性位点。然而，催化剂在有限温度下通常表现出显著的结构灵活性、原子迁移以及与吸附物的动态相互作用。因此，热力学相关的活性位点可能是在操作温度下分布的构型系综，而非单一的静态结构。由于所需的时间和长度尺度限制，利用第一性原理方法显式采样有限温度下的自由能景观在计算上是难以实现的。

方法论
为了应对这些挑战，作者采用了一个结合机器学习原子间势函数（MLIP）、全局优化、分子动力学（MD）以及增强采样自由能计算的多尺度计算框架。

MLIP 开发： 开发了一个深度势能（Deep Potential, DP）模型，该模型基于包含纯金属（Cu, Zn）、合金、氧化物及负载型集群的 44,269 个结构的数据集进行训练。训练数据通过 PBE+D3 水平的 DFT 生成。所得模型在能量上的均方根误差为 0.005 eV/atom，在力上的均方根误差为 0.056 eV/Å，从而在实现接近 DFT 精度的同时，显著提升了模拟速度。
系统模型： 研究调查了两个与甲醇合成相关的集群模型：羟基化的 $Zn_3(OH)_3$ $Z n_{3} (O H)_{3}$ 和含甲酸盐的 $Zn_3(OH)_2CHOO$ $Z n_{3} (O H)_{2} C H O O$ 。这些模型被负载在三种不同的基于 Cu(111) 的表面上，以模拟不同程度的 Zn 合金化程度：
- 原始 Cu(111)（0% Zn）。
- 单 Zn 位点表面（稀释合金化，约 1.56% Zn）。
- 合金化 $Cu_3Zn$ 表面（25% Zn）。
全局优化： 使用带有 DP 模型的遗传算法来识别各表面上低能量的三角形和线性构型，随后进行 DFT 验证。
有限温度采样：
- 无偏 MD： 在 300 K 和 450 K 下进行了 20 ns 的模拟，以观察自发的结构转变和集群迁移率。
- 增强采样： 使用在线概率增强采样（OPES）方法重建了 150、300 和 450 K 下的自由能景观。采用了两个集体变量（CV）：最大 Zn-Zn-Zn 夹角（ $\theta_{max}$ ）用于区分三角形（ $\approx 60^\circ$ ）与线性（ $\approx 170^\circ$ ）几何结构；以及一个基于距离的描述符（ $\tilde{d}_{max}$ ）以加速转变过程。
迁移率分析： 计算了集群质心的均方位移（MSD）和位点相关函数，以评估空间稳定性和抗烧结能力。

关键结果

温度诱导的构型重组：
- 在 0 K 时，对于原始 Cu(111) 上的 $Zn_3(OH)_3$ ，紧凑的三角形构型通常是最小能量结构。
- 有限温度自由能计算显示，随着温度升高，构型表现出向扩展线性构型的显著转变。在 450 K 时，线性构型在热力学上变得与三角形构型具有竞争性。
- 这种转变主要由熵稳定化驱动。虽然不同构型之间的势能差在不同温度下保持相对恒定，但线性构型具有更高的熵，随着温度上升，高熵效应日益抵消其焓值上的劣势。
表面合金化的影响：
- 原始及稀释表面： 在原始 Cu(111) 和单 Zn 位点表面上，随着温度升高，线性系综逐渐趋于稳定，导致在 450 K 时三角形和线性结构的分布比例相当。
- 合金化表面： 在 $Cu_3Zn$ 表面，能量排序甚至在 0 K 时就发生了反转，线性构型比三角形构型更稳定 0.19 eV。广泛的 Zn 合金化在整个温度范围内稳定并硬化了线性构型。
- 吸附物的角色： 甲酸盐（ $Zn_3(OH)_2CHOO$ ）的存在进一步向线性构型偏移了能量平衡，但保留了关于温度和合金化程度的相同定性趋势。
空间稳定性与迁移率：
- 在原始 Cu(111) 上，负载型 $Zn_3$ 集群表现出显著的横向迁移率，表明其具有高度的迁移和凝聚（烧结）倾向。
- Zn 合金化强烈抑制了扩散。在单 Zn 位点表面，尽管其构型平衡与原始表面相似，但孤立表面 Zn 原子나的存在创造了稳定的 $Zn_c\text{-adsorbate-}Zn_s$ 桥联基元。这些基元充当了锚定中心，显著延长了集群的停留时间，并与不稳定的 $Zn_c\text{-Cu}$ 相互作用相比，极大地抑制了迁移。
- 在高度合金化的 $Cu_3Zn$ 表面，扩散几乎可以忽略不计，有效地使孤立的界面 Zn 物种免受烧结影响。

意义与主张
本文证明，负载型 Zn 基集群的热力学相关结构与基于静态 0 K 优化的预测存在根本差异。活性位点系综受焓与熵竞争的影响，其中熵可以在反应条件下驱动从紧凑几何向扩展几何的转变。

作者主张，仅基于最小能量结构来识别催化活性位点是不充分的。相反，活性位点应被视为动态系综，其分布和表面迁移率会随温度、吸附物覆盖度以及局部组成而演变。本研究强调，虽然静态优化可能预测三角形集群占主导地位，但有限温度采样揭示了线性结构存在显著的群体，特别是在操作温度下。此外，这项工作强调了表面组成不仅控制构型平衡，还控制着集群的空间稳定性，其中 Zn 合金化通过创建稳定的锚定基元在防止烧结方面发挥了关键作用。

这些发现为解释实验观察（例如显示三角形结构的低温 STM 数据）提供了一个框架，并建议在操作温度附近，应将表征重点放在涌现出的扩展、动态系综上。该研究验证了使用 MLIP 加速的自由能采样作为弥合理想化静态模型与实际工作中复杂动态结构之间差距的有效途径。

技术摘要：温度诱导的负载型 Cu(111) 上 Zn₃ 集群的重组

类似论文