《Journal of Chemical Theory and Computation》：Periodic GFN2-xTB in CP2K: Multipolar Ewald Electrostatics, k-Point Sampling, and Transferability Benchmarks for Solids
编辑推荐：
扩展紧束缚（extended tight-binding，xTB）方法因其保留了显式电子哈密顿量，且计算成本远低于常规密度泛函理论（density functional theory，DFT），在凝聚相模拟中具有吸引力。然而，若无真正周期性的哈密顿量、布里渊区
扩展紧束缚（extended tight-binding，xTB）方法因其保留了显式电子哈密顿量，且计算成本远低于常规密度泛函理论（density functional theory，DFT），在凝聚相模拟中具有吸引力。然而，若无真正周期性的哈密顿量、布里渊区（Brillouin-zone）采样和对称性感知的晶体计算，其向周期性材料的可迁移性无法被可靠评估。研究人员在CP2K中实现了第二代几何、频率和非共价相互作用参数化的扩展紧束缚（GFN2-xTB）方法的周期性版本。核心理论要素是针对GFN2-xTB的自洽壳层电荷、原子偶极和四极的多极Ewald（multipolar Ewald）公式，以及用于几何优化、晶胞优化和分子动力学的解析力和应力贡献。在DMC-ICE13数据集的冰多晶型物排序中，周期GFN2-xTB在收敛的33k点网格上将相对多晶型物能量优于GFN1-xTB，平均绝对误差从8.01降至3.46 kJ mol?1。由于GFN1-xTB和GFN2-xTB在多个相互耦合、联合参数化的项上存在差异，这种整体模型比较无法孤立多极静电的贡献。对于LC10立方共价和离子固体集合，GFN2-xTB为全部十个体系生成了有效的物态方程最小值，晶格常数和结合能的平均绝对误差分别为0.062 ?和1.30 eV每原子。周期GFN2-xTB使用X23b数据集评估了分子晶格能量和晶胞参数。使用本征Bloch 23网格的全对称性约化晶胞优化，随后在优化晶胞上进行33网格能量评估，周期GFN2-xTB对完整X23b弛豫晶胞集合给出了14.09 kJ mol?1的平均绝对晶格能量误差和5.84%的平均绝对弛豫晶胞体积误差。这些基准测试将周期GFN2-xTB定位为用于筛选、预优化和大规模探索性模拟的高效电子结构层级。该工作为扩展紧束缚框架中的周期性各向异性静电及其在CP2K中的高效应用奠定了基础。
**周期性GFN2-xTB方法在CP2K中的实现与基准评估**
半经验量子化学（semiempirical quantum chemical, SQC）方法在可迁移的经典力场与第一性原理电子结构理论之间占据着重要的中间地位。这类方法保留了显式电子哈密顿量，因此能够提供电荷、轨道、键级、响应性质以及化学键断裂等关键信息，同时通过物理启发的近似和经验参数替代了最昂贵的电子结构计算部分。在分子量子化学领域，从忽略双原子微分重叠（neglect-of-diatomic-differential-overlap, NDDO）模型如MNDO和AM1，到后来广泛适用的PM6方法，这一思想有着悠久的历史。然而，这些方法主要针对孤立分子和团簇设计，其参数集主要由元素或原子价态组织，限制了其在凝聚相体系中的直接应用。经验紧束缚（empirical tight binding, ETB）或扩展Hückel理论则提供了另一种基于DFT的SQC途径，已在密度泛函紧束缚（density functional tight binding, DFTB）框架中形式化。扩展紧束缚（extended tight binding, xTB）方法则通过解析表达式替代数值表格积分，并聚焦于元素特异性而非元素对特异性的参数化，实现了广泛的元素覆盖。GFN家族方法旨在涵盖大部分化学空间，提供可靠的几何结构、振动频率和非共价相互作用。其中GFN1-xTB基于三阶DFTB方案，通过引入环境依赖性参数有效考虑多体效应；GFN2-xTB则通过包含原子偶极和四极的各向异性静电以及密度相关的色散贡献，进一步改善了哈密顿量的物理内容。
GFN1-xTB已主要应用于分子固体和大规模模拟，证明了SQC方法在凝聚相模拟中的价值。尽管tblite库通过软件接口提供了GFN2-xTB参数化，但这与对其自洽偶极和四极进行完整的周期性长程处理是不同的概念。此前已有多项研究尝试将GFN1-xTB扩展到周期性体系，但对GFN2-xTB的各向异性静电项的周期性处理仍缺乏完整的实现。从分子体系扩展到周期性GFN2-xTB远比添加常规的仅含电荷的Ewald求和困难得多，因为GFN2-xTB使用了自洽壳层电荷以及原子偶极和四极，必须将静电核从单极推广到多极Ewald公式，并一致地处理实空间和倒空间项、短程阻尼、自洽势、解析力和应力贡献。此外，分子参数化的半经验哈密顿量在有限体系中可能表现优异，但在测试结合能、晶体体积和布里渊区收敛的多晶型能量学时，可能会显示出系统性误差。因此，对一个真正周期性的GFN2-xTB实现进行定量评估势在必行。
研究人员在CP2K中实现了基于多极Ewald求和的周期性GFN2-xTB方法。该实现集成了tblite库提供模型成分，而CP2K负责周期性模拟机制、电子最小化、力累积和大规模并行执行。实现支持GFN2-xTB多极变量的周期性静电、解析力、用于几何和晶胞优化的应力张量、显式k点采样以及对称性感知的晶体计算。方法评估采用了三个互补的基准数据集：DMC-ICE13数据集提供13种冰多晶型物的扩散蒙特卡罗（diffusion Monte Carlo, DMC）参考晶格能；LC10数据集提供零点点校正的实验和后Hartree-Fock晶格常数及结合能；X23b分子晶体数据集则探测一般的分子间堆积、弛豫晶胞体积和结合能。GFN1-xTB和GFN2-xTB作为完整发布的参数化版本在整个评估过程中进行了比较，二者的差异不仅在于各向异性静电的存在与否，还涉及耦合的哈密顿量、排斥和色散项。
主要关键的技术方法包括：第一，多极Ewald静电公式，这是将GFN2-xTB静电模型引入周期性边界条件的核心步骤，通过实空间短程阻尼贡献与倒空间结构因子的结合，一致地处理壳层电荷的标量、原子偶极的矢量和原子四极的张量相互作用；第二，本征Bloch k点采样机制，通过将图像分辨的哈密顿量和重叠矩阵转换为Bloch形式，在每个采样的晶体动量处独立求解广义本征值问题，并重构图像分辨的密度矩阵，支持Monkhorst-Pack网格和MacDonald偏移网格，并利用spglib库进行完全的对称性约化；第三，解析力与应力张量的推导，力项包含Hellmann-Feynman部分和Pulay修正，应力张量遵循Nielsen-Martin形式的第一性原理应力形式，支持KEEP_SPACE_GROUP过程以保持晶体的初始对称性。
**4.1 基准测试策略**：DMC-ICE13提供了氢键分子晶体的严格首项测试，研究人员以冰Ih为参考计算相对多晶型物能量；LC10扩展至十种立方共价和离子固体，通过本征Bloch 43网格物态方程确定晶格常数，并以匹配的孤立原子参考评估结合能；X23b组合了结合能与完全弛豫的晶胞体积。这些基准区分了数值收敛与模型误差，以及结构可迁移性与结合能准确性。
**4.2 DMC-ICE13冰多晶型物**：研究人员对13种DMC-ICE13参考几何进行了周期GFN1-xTB和GFN2-xTB单点计算。Γ点仅计算给出MAE分别为6.69和5.58 kJ mol?1，但与收敛的布里渊区采样结果存在显著差异，尤其对于冰VII，其GFN2-xTB相对能量误差从Γ点的+33.46降至收敛33网格上的+3.41 kJ mol?1。33网格已有效收敛，更大网格的误差变化小于0.01 kJ mol?1。在33k点水平上，GFN2-xTB将相对能量MAE从8.01降至3.46 kJ mol?1，RMSE从8.65降至3.81 kJ mol?1，但非Ih相仍平均过稳定，最大偏差为冰XIII（-6.35 kJ mol?1）。与最佳DFT近似相比，周期GFN2-xTB尚未达到最精确的冰多晶型物能量描述水平。
**4.3 LC10共价和离子固体**：所有十种体系在两种GFN方法下均产生有效的最小值。完整的GFN2-xTB模型将晶格常数MAE从0.145降至0.062 ?，结合能MAE从1.544降至1.299 eV每原子，相对于GFN1-xTB均有改善。但结合能误差仍远大于SCS-MP2和SOS-MP2参考值，表明分子参数化在定量预测周期性结合能方面仍存在不足。
**4.4 X23b分子晶体**：本征Bloch 23网格晶胞优化配合完全对称性约化和保角晶胞优化方案，给出GFN2-xTB的平均绝对体积误差为5.84%，优于GFN1-xTB的7.51%。平均符号体积误差为-1.66%，但残留误差具有体系依赖性，最大偏差出现在金刚烷、蒽、咪唑和萘。晶格能量比较更为严格：在23网格优化结构上进行33网格能量重评估，X23b晶格能量MAE为14.09 kJ mol?1，对应RMSE为21.34 kJ mol?1。主要的离群体系为咪唑、萘、金刚烷、六亚甲基四胺和吡唑，其中多环芳烃晶体的强过束缚表明当前可迁移分子参数化在无周期性重拟合的情况下无法提供色散主导的分子晶体结合能的平衡描述。
周期性GFN2-xTB在当前实现中的优势在于结构筛选、预优化和大型探索性模拟；对于最终定量周期性结合能，GFN2-xTB应被视为高精度参考的补充而非替代。最直接的未来路线是构建一个结合分子晶体晶格能、弛豫晶胞指标、冰多晶型物相对能量以及周期力和应力的周期训练集，检验能否通过重新平衡扩展体系中的排斥、色散和多极静电场来减少当前偏差。
