Enabling Ab-Initio Molecular Dynamics under Bias: The CP2K+SMEAGOL Interface for Integrating Density Functional Theory and Non-Equilibrium Green Functions

技术摘要：通过 CP2K+SMEAGOL 接口实现偏置下的从头算分子动力学

问题陈述
虽然结合非平衡格林函数（NEGF）的密度泛函理论（DFT）是模拟外部偏置电位下量子输运的标准方法，但传统的计算方法仅限于静态系统。现有的实现通常在固定原子位置的正则框架内运行，计算输运性质。这在理论理解与宏观实验结果之间造成了显著差距，特别是在电化学（EC）系统中，其转化过程本质上是巨正则的。为了准确描述电化学转化的自然环境，模拟需要偏置下的分子动力学（MD），以及能够描述电子进入和离开计算单元的显式开放边界描述。此外，现有的 DFT-NEGF 实现通常缺乏计算几何优化和动力学所需的电流诱导力，或者受限于特定的 DFT 代码（例如 SIESTA），而这些代码缺乏在凝聚相模拟优化的先进软件包中所具备的功能。

方法论
作者展示了一个连接流行的混合高斯/平面波电子结构软件包 CP2K 与 NEGF 代码 SMEAGOL 的新接口。该接口被称为 CP2K+SMEAGOL，旨在作为一个可与任何 DFT 软件链接的独立库，并尽可能地复用原始 SMEAGOL 代码。

关键技术实现包括：

• 电流诱导力： 该接口实现了偏置下作用于离子的非保守力的计算。这些力源自原子动量的时变导数，利用了全哈密顿量对原子位置的梯度。至关重要的是，该实现包含了源自密度矩阵（$\Omega_{ij}$）导数的力项，这解释了“电子风”（电子向离子传递动量）成分。这一项在其他实现中常被忽略，但在本文中被证明对于定性正确性至关重要。
• 密度矩阵构建： 为了正确占据偏置窗口内的束缚态，该接口更新了密度矩阵的计算方式。它没有使用单一轮廓评估或左右贡献的简单平均，而是采用了加权双轮廓方案（遵循 Brandbyge 等人的方法），通过由非平衡贡献决定的权重来对来自两侧电极的积分进行求和。
• 模拟设置： 该方法允许在偏置下进行 Born-Oppenheimer 分子动力学。系统由连接到半无限电极的中心“扩展分子”区域建模。电极的自能被视为静态且独立于中心区域的电荷密度，这一假设在中心区域足以屏蔽电荷变化时是成立的。

关键结果与验证
作者通过针对 SIESTA+SMEAGOL 在多个系统中的参考计算，对 CP2K+SMEAGOL 接口进行了验证：

1. 零偏置力（无限长金线）： 该接口成功重现了零偏置下的原子力，与仅使用 CP2K 的计算结果一致。它还表明，忽略能量密度矩阵项 ($\Omega$) 会导致定性错误的力。
2. 平行板电容器： 该接口重现了电容器中预期的电势降，显示出电极侧恒定的电势以及真空中的线性下降。它证实了加权双轮廓方案对于产生对称电荷密度和在高偏置下产生正确的哈特里电势是必要的，而单轮廓评估则会产生不对称且错误的结果。
3. Au-H$_2$-Au 接点几何优化： 在施加偏置（高达 1.5 V）的情况下，该接口正确预测了 H-H 键的伸长和原子位移，这与之前的 SIESTA+SMEAGOL 结果一致。这些变化归因于电流而非电场。
4. 溶剂化金线的分子动力学： 作者进行了首次在大规模凝聚相系统（由 166 个水分子溶剂化的单原子金线）上进行的 DFT-NEGF 分子动力学模拟。
   • 稳定性： 尽管具有开放系统的特性，动力学过程表现稳定，具有极小的长期能量漂移。
   • 电子迁移： 在 1 V 偏置下的 200 fs 轨迹中，模拟观察到了电子迁移效应，即局部电子密度极小值从中心原子移动到了最右侧的原子，这与屏蔽右侧电极正电荷的现象一致。
   • 传输： 溶剂化金线的传输谱显示出与理想无限长金线定性的相似性，差异归因于水分子开启了额外的隧穿通道。

性能
主要的计算瓶颈是评估作为格林函数积分的电子密度，这在每个 SCF 步中需要 96 个能量点计算（64 个实轴点，32 个复轴点），取代了单一的哈密顿量对角化。

• 并行化： 代码利用 MPI 进行实空间和 k 点计算，并利用 OpenMP 进行 DO 循环和 LAPACK 例程计算。性能在 32 个 OpenMP 线程下表现良好。
• 成本： 对于溶剂化金线系统，尽管使用了 8 倍更多的核心（2048 对 256），CP2K+SMEAGOL 的每个 SCF 步计算速度仅比标准 CP2K 慢约 3.2 倍。作者指出，这代表了数量级的成本增加，使得在没有进一步加速（如通过 ScaLAPACK 或机器学习力场）的情况下进行长时间 MD 轨迹变得具有挑战性。

意义与主张
论文声称，CP2K+SMEAGOL 是目前针对 CP2K 最功能丰富且最完整的 DFT-NEGF 实现，超越了原生 CP2K 实现（后者仅限于不带力的 $\Gamma$ 点输运）。

这项工作的主要意义在于证明了 DFT-NEGF 可用于在现实运行条件下对大规模凝聚相系统进行偏置下的分子动力学模拟。据作者所知，这是首次完成此类计算的工作。该接口使得研究电催化反应和电化学系统成为可能，在这些系统中，施加的电势、电流流动与原子动力学之间的相互作用至关重要，从而弥合了微观理论模型与宏观实验条件之间的鸿沟。