Reaching the thermodynamic limit of periodic CCSD cohesive energies and band gaps with denser Brillouin zone sampling

本文提出了一种可扩展的分布式内存周期性 CCSD 理论实现方案，该方案能够实现密集的布里渊区采样（高达 216 个 k 点），从而可靠地将内聚能和带隙外推至热力学极限，并为八种半导体和绝缘体提供了确定的基准值，其与实验数据的误差分别约为 0.1–0.2 eV 和 0.4 eV。

要点

想象一下，你正试图理解一个巨大、完美的晶体（如一颗巨大的钻石或一块食盐）的行为。在现实世界中，这些晶体非常庞大，包含数以万亿计的原子。然而，为了在计算机上进行研究，科学家通常必须将问题缩小到一个微小的、易于处理的“微型晶体”（超胞），并假定它会无限重复。

这种微型晶体方法的问题在于，这就像试图通过只观察一个后院来了解整个大陆的天气。你会错过大局，从而产生“有限尺寸误差”——即由于样本过小而导致的错误。

挑战：“像素化”的晶体
在量子化学的世界里，你样本的“大小”由你在被称为**布里洛因区（Brillouin zone）**的数学空间中采样的点数决定。可以将这些点视为数字图像中的像素。

• 低分辨率（像素较少）： 你得到的是一张模糊、不准确的晶体属性图像（例如原子结合得有多紧，或者电子跨越能隙所需的能量）。
• 高分辨率（像素较多）： 你得到的是一张清晰、准确的图像。

问题在于，使用被称为 CCSD（含单激发和双激发的耦合集群理论）的金标准方法来计算这些属性，其计算成本极其昂贵。这就像是用计算器去渲染一部 4K 电影。由于这种方法计算量巨大，以往的研究只能负担得起“低分辨率”的图像（较小的像素网格），从而被迫去猜测高分辨率图像的样子。这些猜测往往会导致显著的误差。

解决方案：超级强大的计算机团队
本文作者构建了一个全新的、超高效的软件程序，它就像一个庞大的工人团队（可以在多达 12 个计算节点、每个节点 96 个核心的配置下运行）。这个团队可以协同工作，处理比以往任何时候都更大的“图像”。

他们不再仅仅观察微小的像素网格，而是能够采样高达 216 个点（一个 $6 \times 6 \times 6$ 的网格）在布里洛因区内。这就像是从一张模糊的缩略图升级到了高清 4K 图像。凭借这种全新的清晰度，他们终于能够直接观测到真实的、“热力学极限”下的完美无限晶体行为，而无需进行任何猜测。

发现： “金标准”基准测试
利用这种高清晰度的方法，团队计算了八种常见材料（如氧化镁、硅和钻石）的两种关键属性：

1. 结合能（Cohesive Energy）： 将晶体拆解为单个原子所需的能量。
2. 带隙（Band Gap）： 电子导电所需的能量“间隙”（本质上决定了材料是绝缘体还是半导体）。

他们将预测结果与现实世界的实验进行了对比：

• 结合能： 他们的预测结果非常接近现实，误差通常仅为 0.1 到 0.2 eV。他们倾向于稍微低估了原子结合的紧密程度。
• 带隙： 他们的预测也相当不错，误差约为 0.4 eV，但他们倾向于稍微高估了能隙（即预测该材料是比实际情况更佳的绝缘体）。

“间接”之谜
某些材料具有“间接”带隙，这类带隙的计算更为棘手。这就像是试图测量两个彼此无法直接观察到的点之间的距离。作者发现，标准的猜测方法（外推法）在这里经常失效，会导致对能隙的低估。他们开发了一种聪明的“复合”策略——先测量直接路径，然后加上针对间接路径的修正——从而获得更准确的结果。

二氧化钛测试
为了证明其方法在更复杂材料上的有效性，他们将该方法应用于金红石型二氧化钛（一种常见的白色颜料和光催化剂）。他们的计算预测带隙为 4.17 eV。这比实验值（经修正后约为 3.9 eV）略高，但作者指出，这种微小的误差与其已知方法的局限性是一致的，并表明可能需要更复杂的物理过程（如三激发）才能达到完美的准确度。

底线
本文不仅仅提供了一些新的数据；它提供了一个决定性的基准。通过使用庞大的计算机团队来生成“高清晰度”的数据，作者为 CCSD 方法的效果建立了一个新的、可靠的标准。他们表明，虽然该方法非常优秀，但在与现实世界对比时，仍然存在大约 0.3–0.4 eV 的可预测“模糊度”（误差）。这有助于其他科学家明确了解在未来可以对类似计算产生多少信任。

技术摘要

技术摘要：实现周期性 CCSD 内聚能与带隙的热力学极限

问题陈述
包含单激发与双激发的周期性耦合簇理论（CCSD）是一种计算固体电子结构的极具前景的非微扰方法，相比于密度泛函理论（DFT）等平均场理论，它具有系统性的可改进性。然而，其应用在热力学极限（TDL）——即消除有限尺寸误差的物理相关极限——方面受到了计算成本的严重限制。周期性 CCSD 的存储和执行时间分别随 $O(n^4 N_k^3)$ 和 $O(n^k N_k^4)$ 缩放（其中 $n$ 是每个晶胞中的原子数，$N_k$ 是 $k$ 点数）。因此，以往的研究被限制在 $k$ 点网格不超过 $4^3$ 的三维固体研究中。这种限制阻碍了可靠向热力学极限的外推，特别是对于收敛速度（$N_k^{-1/3}$）比基态能量（$N_k^{-1}$）更慢的带隙而言。由于无法获得足够密集的 $k$ 点采样，内聚能和带隙的预测会受到未量化有限尺寸误差的影响，从而阻碍了建立 CCSD 层级理论的确定性基准。

方法论
为了解决这些局限性，作者开发并应用了一种基于 PySCF 框架的新型、优化的分布式内存周期性 CCSD 和方程运动 CCSD（EOM-CCSD）软件实现。关键方法论特征包括：

• 并行化与对称性： 该代码利用分布式和共享内存并行技术，并利用空间群对称性在不可约布里渊区进行计算，从而实现了在多达 12 个节点（每个节点 96 个核心）及每节点 1.5 TB 内存上的计算。
• 采样密度： 该实现允许在双 $\zeta$（DZ）基组中进行高达 $N_k = 6^3 = 216$ 个 $k$ 点的布里渊区采样，在三 $\zeta$（TZ）基组中为 $5^3$，在四 $\zeta$（QZ）基组中为 $4^3$，且无需近似。
• 外推方案： 本研究采用了特定的外推模型以达到热力学极限：
  • 对于内聚能，Hartree-Fock（HF）能量使用三参数模型（$N_k^{-1} + N_k^{-2}$）进行拟合，而相关能（MP2 和 CCSD）使用二参数模型（$N_k^{-1}$）。
  • 对于带隙，作者评估了三种模型：模型 A（领先阶 $N_k^{-1/3}$）、模型 AB（包含至二阶）以及模型 ABC（包含至三阶，包括 $N_k^{-1}$ 项）。
  • 还测试了一种“GW-EOM”方案，该方案假设 EOM-CCSD 与 $G_0W_0@HF$ 计算之间的有限尺寸误差系数成比例，以此预测较小 $k$ 网格下的热力学极限带隙。
• 间接带隙的复合方法： 考虑到间接带隙（电离能与电子亲和能之和）由于导带底（CBM）需要 $k$ 网格偏移而难以收敛，本文提出了一个复合方法。该方法将计算分为价带顶（VBM）处的直接带隙以及 CBM-VBM 偏移，从而实现更稳定的外推。
• 基组与势函数： 计算使用相关一致基组（cc-pVXZ，$X=D, T, Q$）和 GTH 伪势。基组不完备误差通过外推（例如 $X^{-3}$ 形式）或通过比较 DZ、TZ 和 QZ 的结果来进行修正。

主要贡献与结果
作者报告了八种简单半导体和绝缘体（C, Si, BN, BP, MgO, LiH, LiF, LiCl）以及一种新应用——金红石型 TiO$_2$ 的收敛 CCSD 内聚能和带隙。

• 内聚能： 本研究提供了收敛至约 0.1 eV 的 CCSD 内聚能的确定性基准数值。

  • CCSD 往往比实验值（经零点能修正后）低估内聚能，平均低估 0.1–0.2 eV，而 MP2 则倾向于高估。
  • 研究发现 CCSD 的准确性与标准 DFT 泛函（PBE, PBEsol, SCAN）相当或更好。
  • 对于大多数固体，伪势误差通常很小（<0.05 eV），但在 BN 和 C 中略大。

• 带隙： 本文确立了可靠的带隙热力学极限外推需要比以往可及的规模显著更大的 $k$ 网格（$N_k \ge 3^3$）。

  • 外推模型： 当应用于小网格（$N_k \le 4^3$）时，简单的单项外推（模型 A）一致地低估带隙 0.3–0.5 eV。经过 $N_k=3^3$–$6^3$ 拟合的三参数模型 ABC 被确定为最稳健的参考模型。
  • GW-EOM 方案： 使用 $G_0W_0$ 来修正 CCSD 有限尺寸误差的 GW-EOM 方法表现出可靠性差异。它对某些材料（MgO, LiCl, LiF）效果良好，但对其他材料（BN, BP, Si, C）则高估带隙 0.4–0.5 eV，这表明误差系数的比例假设并不普遍成立。
  • 间接带隙： 复合方法成功解决了间接带隙（如 BN）的收敛问题，揭示了 VBM 处的直接带隙外推是稳定的，而 CBM 偏移的收敛速度很快。
  • 准确度： 最终的 TDL EOM-CCSD 带隙相对于经零点能修正后的实验值，其平均绝对误差（MAE）在采用伪势时为 0.39 eV，在全电子计算时为 0.42 eV。该方法通常会高估带隙。
  • TiO$_2$： 对于金红石型 TiO$_2$，TDL EOM-CCSD 的预测值为 4.17 eV，比经零点能修正后的实验值（~3.9 eV）高出约 0.3 eV。

• 单激发特性： 作者指出，与实验偏差较大（特别是岩盐结构如 MgO 和 LiH）的材料，其单激发特性较低（$n_{1}^{IP} \approx 93-94\%$），而误差较小的材料（$>95\%$），这表明对于这些系统，三激发对于提高精度可能是必要的。

意义
这项工作是迄今为止基于正则轨道开展的最大规模周期性 CCSD 计算，成功解决了有限尺寸和基组误差，为热力学极限下的 CCSD 内聚能和带隙提供了首个可靠的基准。研究表明，虽然 CCSD 在内聚能方面具有高精度，但其在带隙方面的应用对外推方案和 $k$ 点密度非常敏感。作者得出结论，对于受限于 $N_k \le 4^3$ 的实际计算，AB 外推模型是目前最准确的方法，尽管相对于完全收敛的 TDL 仍存在约 0.25 eV 的残余误差。结果强调，EOM-CCSD 与实验之间剩余的差异（0.3–0.4 eV）可能源于电子相关处理不足（特别是对三激发的处理需求）或电子-声子相互作用，而非有限尺寸误差。这为更经济高效的方法设定了新的基准标准，并阐明了 CCSD 在固体电子结构研究中的当前局限性。