非阿贝尔点群情况下多体方法中对称性利用的直接乘积分解方法

《The Journal of Physical Chemistry A》：The Direct-Product Decomposition Approach for Symmetry Exploitation in Many-Body Methods in Case of Non-Abelian Point Groups

　　

研究人员针对C_3v特例，展示了用于处理耦合簇（CC）计算中对称性的直接乘积分解方案如何扩展到非阿贝尔点群。研究人员表明，通过将两个不可约表示的可约乘积分解为其不可约表示，对于双电子积分和CC振幅可以获得块结构。为了处理必要的双电子积分和振幅排序重排、自旋适应以及CC计算中的??(M⁵)收缩（M为基函数数量），研究人员提出了一种策略，该策略同时使用相应量的约化表示和非约化表示，并在两者之间来回切换。约化表示用于??(M⁶)收缩，而其他步骤则更适合在非约化表示中进行。研究人员对采用不同基组进行NH₃和PH₃测试计算时实施的CC单双激发（CCSD）方法的试点实现证实，可实现显著节省（与无对称性处理相比超过20倍，与使用C_s对称性处理相比约5倍），这些发现表明，利用非阿贝尔对称性将使大型高对称分子的CC计算成为可能。

以下是对论文《非阿贝尔点群情况下多体方法中对称性利用的直接乘积分解方法》的解读，该论文发表在《The Journal of Physical Chemistry A》。

**研究背景与问题**
对称性是自然科学各领域的重要概念。在化学中，对称性可用于解释化学键（分子轨道理论、配体场理论）、通过选择规则解释光谱以及通过Woodward–Hoffmann规则解释化学反应性等。在量子化学计算中，利用对称性可以降低计算成本和内存需求。对称性在量子化学中的应用包括：积分评估和自洽场（SCF）步骤、原子轨道（AO）到分子轨道（MO）的积分变换以及电子相关处理。对于电子相关处理，Stanton等人提出了直接乘积分解（DPD）方案，该方案在阿贝尔点群（如D_2h及其子群）中非常有效，能实现所有量的对称分块且无显著开销。然而，原始DPD方案仅限于阿贝尔点群，因为它依赖于两个不可约表示的直积仍是另一个不可约表示这一性质，而这对于非阿贝尔点群不成立。因此，现有研究无法直接处理非阿贝尔点群的对称性，限制了大型高对称分子（如C₆₀）的耦合簇（CC）计算。本研究旨在将DPD方案扩展到非阿贝尔点群，以C_3v为例展示该方法，并评估其计算节省潜力。

**研究内容与结论**
研究人员针对C_3v点群（具有6个对称操作和3个不可约表示A₁、A₂、E），提出了DPD方案的扩展。关键思路是通过对两个不可约表示的可约乘积（如E?E）进行约化，分解为不可约表示的直和，从而恢复四维量（如双电子积分、CC振幅）的对角块结构。研究人员同时使用约化表示和非约化表示，对于??(M⁶)收缩使用约化表示，而对于重排、自旋适应和??(M⁵)收缩等步骤，则先在非约化表示中处理，再通过定义的“reduce”例程来回切换。该方法解决了非阿贝尔点群下“缺失积分/振幅”问题（即非约化表示中缺少的部分可通过已有量重构）。研究人员在试点程序QUENA中实现了CCSD方法（CC单双激发），利用CFOUR计算的轨道能量和双电子积分，对NH₃和PH₃采用不同基组（cc-pVDZ、cc-pVTZ、cc-pVQZ）进行测试。结果表明，与无对称性（C₁）处理相比，操作数节省约21–22倍；与仅采用最大阿贝尔子群（C_s）处理相比，节省约5.5–6倍。节省因子在涉及更多虚拟轨道的项（如环项和粒子-粒子梯项）中最大。研究结论认为，利用非阿贝尔对称性能显著降低CC计算成本，使大型高对称分子的计算成为可能。

**主要关键技术方法**
（1）直接乘积分解（DPD）方案：将两个不可约表示的可约乘积（如E?E）通过投影算子分解为不可约表示（A₁、A₂、E）的直和，从而恢复四维量的对角块结构并消除冗余。
（2）约化与非约化表示切换策略：对于??(M⁶)收缩使用约化表示，而对于重排、自旋适应和??(M⁵)收缩，先回退到非约化表示（基于最大阿贝尔子群C_s的排序）处理，再通过reduce例程转换为约化表示。
（3）缺失积分/振幅重构：利用已有量通过线性组合（如式12）重构非约化表示中缺失的元素。

**研究结果**
1. **2.1 一般概念**：对称选择规则决定哪些量需要计算和存储。分子轨道（MO）按不可约表示分类，可形成块结构。
2. **2.2 阿贝尔点群的DPD**：对于阿贝尔点群，所有不可约表示一维，两个不可约表示的直积仍是一维，因此四维量可通过超指标构造二维块矩阵，实现块状收缩。
3. **2.3 扩展到C_3v对称性**：非阿贝尔点群（如C_3v）中直积E?E = A₁⊕A₂⊕E，导致直接块结构含非对角块。通过投影算子进行约化（公式5–8），可恢复对角块结构（图5），去除冗余。
4. **2.4 非约化与约化表示**：重排和自旋适应在非约化表示中简单，而??(M⁶)收缩在约化表示中高效。通过reduce例程来回切换（图2）。
5. **2.5 四指标量重排**：在非约化表示中进行重排，若遇缺失元素，通过重构公式（式12）补齐，再约化。
6. **2.6 自旋适应**：闭壳层自旋适应需先回退非约化表示，进行自旋适应，再约化；缺失元素同样通过式12重构。
7. **2.7 ??(M⁵)收缩**：两种类型收缩（四指标→两指标和两指标→四指标）均在非约化表示中进行，需处理缺失元素和E表示二维性带来的额外因子。
8. **3. 实现与计算细节**：在QUENA程序中实现CCSD方程（式16–33），基于CFOUR的C_s对称性结果，通过轨道选择和约化步骤转换为C_3v。
9. **4. 探索性计算**：对NH₃和PH₃采用cc-pVQZ基组的操作数计数表明，C_3v对称性相比C₁节省约21.9倍（NH₃）和21.6倍（PH₃），相比C_s节省约5.5–6倍。节省因子在环项和粒子-粒子梯项中最大（约22.3和20.9）。
10. **5. 结论与展望**：本研究展示了如何将阿贝尔点群的DPD方案扩展到非阿贝尔点群（以C_3v为例）。通过使用约化和非约化表示并来回切换，解决了重排、自旋适应和??(M⁵)收缩中的问题。试点程序QUENA的测试表明，充分利用C_3v对称性可带来显著节省。预计该概念可推广到其他非阿贝尔点群，并扩展到三重激发的处理（如CCSD(T)）。

**结论翻译**：
在本工作中，研究人员证明了Stanton等人用于处理CC计算中阿贝尔点群对称性的直接乘积分解方法如何扩展到非阿贝尔点群。关键发现是，通过将两个不可约表示的可约乘积（在C_3v情况下为E?E）分解为其不可约表示，可以获得高维量（如双电子积分和振幅）的块结构。这也允许消除冗余（即仅保留一个E块），并使得??(M⁶)收缩能够以与阿贝尔情况相同的方式执行，而不会因对称性利用产生显著开销。然而，在重排、自旋适应以及??(M⁵)收缩中出现了复杂性。为克服这些复杂性，研究人员提出了一种策略，同时使用所涉及高维量的约化表示和非约化表示，因为所有问题步骤在非约化表示中处理起来相当直接（除了需要处理“缺失积分/振幅”问题，这是由于计算中并非所有积分均可用，但可通过可用量重构）。无需同时存储两种表示，因为可以在它们之间轻松切换，这一策略在QUENA新程序包中得以采用。该程序包含针对C_3v对称性的试点实现，并且对称性处理的部分目前针对该点群进行了硬编码。关于NH₃和PH₃的探索性计算表明，充分利用C_3v点群对称性相比在C₁（节省约21–22倍）和C_s（节省约5.5–6倍）对称性下的对应计算可带来显著节省。研究人员预计，将该概念推广到其他非阿贝尔点群是直接的。相关研究工作正在进行中，同时也在进行将当前对称性利用概念扩展到CCSD(T)框架中三重激发微扰处理的工作。