Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

本文提出了一种针对无限周期体系的高效、线性缩放的从头算多体量子嵌入与局部关联方法实现，该方法利用插值可分密度拟合技术，能够为弱关联和强关联固体提供精确的热力学极限耦合簇基态能量估计。

要点

想象一下，你正试图通过研究一间单体房屋来理解一座庞大、无限的城市是如何运作的。在量子化学的世界里，这座“城市”是一个晶体（如钻石或金属氧化物），而这间“房屋”则是被称为“晶胞”的微小重复单元。科学家们想要知道这座无限城市的精确能量，以预测其性质，但计算这个无限网格中每一个电子之间的相互作用，就像试图同时计算地球上每一片沙滩上的每一粒沙子一样。使用传统方法，这在计算上是不可能的，因为随着网格变大，工作量增长得太快。

本文介绍了一种巧妙的新型“捷径”，称为插值可分密度拟合 (Interpolative Separable Density Fitting, ISDF)，并结合了 FFTISDF 技术来解决这个问题。以下是其工作原理的简单类比：

问题：“邻居太多”的问题

在晶体中，电子不仅与直接相邻的邻居发生相互作用，还能感受到远处电子的引力。为了获得准确的答案，你需要用一个点阵（称为 k 点）对晶体进行采样。

• 旧方法： 想象一下，试图通过询问体育场里的每一个人他们听到了什么，来计算整个体育场的噪音水平。如果你增加人数（k 点），你需要追踪的对话数量就会爆炸式增长。这就是为什么以前的方法在模拟大型无限晶体时会遇到瓶颈。
• 目标： 作者希望使用高达 1,000 个 k 点（一个非常密集的网格）来计算这些无限晶体的能量，从而得到代表“热力学极限”（即材料真实无限大小）的结果。

解决方案：“智能总结者”

作者开发了一种类似于智能总结者或翻译官的方法。

1. 插值点（“关键证人”）：
   该方法并没有询问晶体中每一个电子关于其相互作用的情况，而是挑选了一组精选的小规模“证人”（称为插值点）。可以把它们想象成新闻编辑室里的关键记者。与其采访每一位市民来了解城市的情绪，不如采访一组精心挑选的 100 个人，他们就能准确地代表整座城市的情绪。

   • 论文表明，通过使用这些“证人”，他们可以高精度地重建整个电子云的行为，但所花费的工作量却极小。

2. 线性缩放（“魔法电梯”）：
   在旧方法中，如果你将模拟规模增加一倍（更多的 k 点），所需的时间会增加四倍甚至更多（就像是在爬一座陡峭且无尽的山坡）。

   • 使用这种新方法后，计算所需的时间呈线性增长。如果你将 k 点数量增加一倍，计算时间也仅仅增加一倍。这就像拥有一部魔法电梯，无论山有多高，你都能轻松上升而不会感到疲惫。这使得他们能够运行高达 1,000 个 k 点的模拟，这在以前是不可能的。

工具：“嵌入”与“局部相关性”

为了获得最准确的能量数值，论文使用了两种特定策略：

• 密度矩阵嵌入（“焦点小组”）： 这就像是选取一小群具有代表性的人（晶体的一个片段）进行深度研究，同时将剩余的城市视为简化的背景。这使得能够对“局部”相互作用进行非常精确的观察。
• 局部自然轨道相关性（“高效分类”）： 这种方法对电子进行排序，使得只有那些对特定相互作用真正重要的电子才会被进行高细节计算，而忽略那些距离太远而无关紧要的电子。

测试内容

团队在四种不同类型的材料上测试了这种新的“智能总结者”：

1. 钻石： 一种坚硬的宽禁带半导体。
2. 二氧化碳 (CO2)： 一种分子晶体（类似于干冰）。
3. 氧化镍 (NiO)： 一种电子“强相关”的材料（电子表现得像混乱的群体而非独立的个体）。
4. CaCuO2： 一种具有层状结构的铜氧化物超导体。

实验结果

• 准确度： 他们证明了该方法可以极其精确地预测这些材料的能量，其精度能与速度慢得多的旧方法相媲美，但完成时间仅为后者的一小部分。
• “热力学极限”： 通过使用高达 1,000 个 k 点并进行数学上的“外推”（预测向无穷大的趋势），他们能够为这些无限晶体的基态能量提供迄今为止最准确的估计。
• 磁性质： 对于氧化镍和 CaCuO2，他们计算了原子之间如何进行磁性相互作用（具体为“交换常数”）。他们的结果比以往的计算更接近真实的实验值，这证明了包含这些“强相关性”对于理解这些材料至关重要。

核心结论

本文提出了一种新的计算引擎，使得模拟无限晶体变得可行，其细节程度可以达到以往仅限于微小分子研究的水平。通过使用“智能总结者”（ISDF）来降低电子相互作用的复杂性，他们将一项在计算上几乎不可能完成的任务转变为高效且可扩展的任务。这使得科学家们终于能够获得关于固体材料真实、无限本质的可靠答案，而无需依赖一颗行星大小的超级计算机。

技术摘要

技术摘要：利用插值可分密度拟合实现晶体材料的从头算多体量子嵌入与局部相关计算

问题陈述
关联波函数方法通过系统可改进的多体展开，为实现计算材料科学中的可控精度提供了路径。然而，由于其对系统规模具有陡峭的计算缩放特性，将其应用于无限周期系统受到了阻碍，这成为了获得热力学极限（TDL）结果的重要障碍。虽然现有策略利用平移对称性（布洛赫轨道）将成本降低了 $N_k^{-1}$ 到 $N_k^{-2}$ 倍，并采用局部性或低秩近似（如 Cholesky 分解或密度拟合/分辨率恒等式）来管理相关性，但仍缺乏一种能够实现高阶多体方法随 k 点数量 ($N_k$) 线性缩放的统一方法。具体而言，标准的密度拟合技术在应用于实现 TDL 收敛所需的庞大 k 点网格时，往往会面临二次方缩放或难以承受的存储需求。

方法论
作者提出了一种针对无限周期系统的高效从头算多体量子嵌入与局部相关方法实现。该方法的核心是将**平移对称性适配的插值可分密度拟合（FFTISDF）**集成到高斯-平面波框架（在 PySCF 中实现）中。

关键技术组件包括：

1. FFTISDF 表示法： 库仑相互作用通过张量超收缩（THC）形式进行表示。通过基于旋转 Cholesky 分解确定的选定插值点（IPs），对实空间网格上的基函数乘积进行近似。这使得二电子积分可以表示为因子化形式：
   $$(\mu k_\mu \nu k_\nu | \lambda k_\lambda \sigma k_\sigma) = \sum_{IJ} X^{k_\mu}_{I\mu} X^{-k_\nu}_{I\nu} X^{k_\lambda}_{J\lambda} X^{-k_\sigma}_{J\sigma} W^q_{IJ}$$
   其中 $X$ 是插值向量，$W$ 是库仑核。
2. 线性缩放算法： 通过利用卷积定理和快速傅里叶变换（FFT），加速了度规张量、库仑核和交换矩阵的评估。这使得积分评估、库仑矩阵构建和交换矩阵构建的缩放严格降低至关于 $N_k$ 的线性水平。
3. 量子嵌入与局部相关： FFTISDF 形式被应用于：
   • 密度矩阵嵌入理论（DMET）： 通过使用因子化积分将 k 空间晶体轨道哈密顿量转化为片段空间，从而加速片段哈密顿量的构建。
   • 局部自然轨道（LNO）相关： 促进全局低阶相关计算（如 k-SOS-MP2, k-dRPA）以及局部高阶修正（LNO-CCSD, LNO-CCSD(T)）。
4. 外推策略： 作者采用了两种外推技术以达到热力学极限（TDL）：
   • k 点外推： 对能量相对于 $N_k^{-1}$ 进行线性拟合以进行 k 点采样。
   • 相关域外推： 使用 k-SOS-MP2 作为基准，通过改变截断阈值，将 LNO-CCSD 和 LNO-CCSD(T) 能量外推至完整的相关域。

主要贡献

• 算法实现： 实现了一种适用于周期系统的稳健 FFTISDF 算法，能够实现多体计算中随 $N_k$ 线性缩放，克服了传统高斯密度拟合（GDF）和 FFT 密度拟合（FFTDF）在交换和哈密顿量构建中存在的二次方缩放瓶颈。
• 可扩展性演示： 在包含高达 1000 个 k 点的晶体系统上，成功执行了关联基态计算（DMET 和基于 LNO 的 CCSD/CCSD(T)）。
• 存储效率： 与 GDF 的 $N_k^2 N_{AO}^2 N_{aux}$ 缩放相比，存储需求大幅降低，缩放比例为 $N_k N_{IP}^2 + N_k N_{IP} N_{AO}$。对于一个 1000 k 点的金刚石计算，存储量从约 881 GB 减少到了约 2.1 GB。
• 热力学极限估计： 具备将 CCSD(T) 基态能量外推至热力学极限的能力，适用于包括弱相关和强相关固体在内的系统。

结果
该方法在四个代表性系统上进行了基准测试：

1. 金刚石（宽带隙半导体）：
   • 使用 $c_{IP}=14$ 的 FFTISDF 相对于高截断 FFTDF 参考值，实现了每个原子的能量误差低于 $10^{-4}$ a.u.。
   • 墙钟时间分析证实了 FFTISDF 在所有步骤（积分、库仑、交换、哈密顿量）中均符合线性缩放，而 GDF 和 FFTDF 在交换和哈密顿量构建中表现出二次方缩放。
   • TDL 外推得到了收敛的 CCSD(T) 能量，证明了该方法能够采样由于成本限制而无法通过标准方法实现的 k 点收敛过程中能量增加的阶段。
2. 二氧化碳（分子晶体）：
   • 验证了一种外推程序，即通过将 LNO 能量（截断域）相对于全局 k-SOS-MP2 能量的差值（$\Delta_{LNO}$）进行线性拟合。
   • 外推得到的 TDL 能量（MP2, CCSD 和 CCSD(T)）与全局参考值表现出极佳的一致性（每个原子的误差 $\sim 10^{-5}$ 至 $10^{-4}$ a.u.），证实了域外推策略的有效性。
3. 氧化镍（NiO）和铜酸钙（CaCuO2，强相关铜氧化物）：
   • 应用 FFTISDF-DMET 结合 CCSD 解算器来确定磁交换常数（$J$）。
   • 对于 NiO，DMET 得出的次近邻交换常数 $J_2 = -14.05$ meV，较 k-HF 的低估值（$-7.59$ meV）有所改善，但仍略低于实验范围。
   • 对于 CaCuO2，DMET 得出 $J = -165.44$ meV，其量级显著大于 k-HF 值（$-38.97$ meV），并且与源自自旋波色散的实验值（$-142$至$-158$ meV）高度吻合。

意义与主张
本文声称，这项工作建立了一个稳健的计算框架，能够为晶体在热力学极限下的高阶关联方法（DMET, LNO-MP2, LNO-CCSD, LNO-CCSD(T)）获取收敛的基态能量。其主要意义在于相对于 k 点数量的线性缩放，这使得获取此前在计算上难以实现的 TDL 结果变得切实可行。

作者指出，尽管该方法取得了成功，但目前仍需要相对较多的插值点（通常为原子轨道的 10 倍）才能达到所需的精度，这使得为大型单元格确定这些点在计算上非常繁重。此外，对于含有过渡金属的系统，实空间网格尺寸仍然较大。作者总结道，虽然该框架是一个重要的进步，但未来的工作应侧重于开发更高效的协议来生成紧凑的插值点，并探索 LNO 方法的其他相关空间。