[论文解读] Generation of Spin-Adapted and Spin-Complete Substitution Operators for (High-Spin) Open-Shell Coupled Cluster of Arbitrary Order

Nils Herrmann, Michael Hanrath|arXiv (Cornell University)|Aug 23, 2020

Advanced Chemical Physics Studies参考文献 38被引用 9

一句话总结

本文提出了一种严格的四步算法，用于生成任意激发阶数和自旋量子数 S 的高自旋开壳层耦域簇（CC）方法的线性无关且自旋完备的空间取代算符。通过利用 Löwdin 的投影算符方法，该方法确保了自旋完备性并避免了自旋污染，从而实现了对高自旋态（如硼原子中的态）的精确 CC 计算，结果表现出可忽略的自旋投影误差以及与自旋轨道 CC 方法相比微小的关联能差异。

ABSTRACT

A rigorous generation of spin-adapted (spin-free) substitution operators for high spin ($S=S_z$) references of arbitrary substitution order and spin quantum number $S$ is presented. The generated operators lead to linearly independent but non-orthogonal CSFs when applied to the reference and span the complete spin space. To incorporate spin completeness, spectating substitutions (as e.g. $\hat{E}_{iv}^{va}$) are introduced. The presented procedure utilizes L\"owdin's projection operator method of spin eigenfunction generation to ensure spin completeness. The generated operators are explicitly checked for (i) their linear independence and (ii) their spin completeness for up to 10-fold substitutions and up to a multiplicity of $2S+1 = 11$. A proof of concept implementation utilizing the generated operators in a coupled cluster (CC) calculation was successfully applied to the high spin states of the Boron atom. The results show pure spin states as well as small effects on the correlation energy compared to spinorbital CC. A comparison to spin-adapted but spin-incomplete CC shows a significant spin incompleteness error.

研究动机与目标

• 开发一种系统化、严谨的程序，用于生成高自旋开壳层耦域簇方法的自旋自洽且自旋完备的取代算符。
• 解决开壳层 CC 中长期存在的自旋不完备性问题，该问题会导致虚假的自旋污染和关联能误差。
• 确保最多十阶激发和多重度 2S+1=11 的空间取代算符的线性无关性和自旋完备性。
• 提供一种可推广的、自动化的框架，适用于任意高自旋参考态，避免依赖自旋轨道形式体系。

提出的方法

• 该方法采用 Löwdin 的投影算符技术生成自旋本征态并确保自旋完备性。
• 系统地生成所有给定阶数的原型空间取代算符（ˆE），这些算符在作用于参考 CSF 时可产生不同的空间函数。
• 采用拓扑映射方法分配自旋粒子到空间轨道的排列，以保持自旋对称性。
• 对轨道空间（占据、虚拟、活性）应用规范索引排序，以标准化算符定义并避免冗余。
• 该过程模块化，作为黑箱输入输出系统运行：给定取代阶数 m，返回一组线性无关且自旋完备的 ˆE 算符。
• 通过在确定子基中对 CSF 表示的矩阵秩分析进行验证，并使用 Weyl-Robinson 维数公式确认。

实验结果

研究问题

• RQ1 是否能够通过一种系统化、可推广的方法，为任意高自旋开壳层参考态和激发阶数生成线性无关且自旋完备的空间取代算符？
• RQ2 在不依赖自旋轨道形式体系或截断的 CSF 流形的前提下，如何严格强制实现耦域簇理论中的自旋完备性？
• RQ3 自旋不完备性对开壳层 CC 计算中关联能的影响是什么？是否可以对其进行量化？
• RQ4 与完全自旋完备的实现相比，使用自旋自洽但自旋不完备的算符是否会引入显著误差？

主要发现

• 所提出的方法成功生成了最多十阶取代和多重度 2S+1=11 的线性无关且自旋完备的空间取代算符，结果通过矩阵秩和维数公式检查得到确认。
• 使用生成算符的原理性 CC 实现被成功应用于硼原子的高自旋态，获得了自旋纯态，自旋投影误差为零。
• 自旋自洽且自旋完备的 CC 方法所得关联能与自旋轨道 CC 方法相比仅存在微小偏差（四重态情况下约为 0.75%），表明其具有高精度。
• 直接对比显示，自旋自洽但自旋不完备的 CC 方法在关联能上存在持续 0.4–0.05 mH 的误差，证明了自旋完备性的重要性。
• 在双精度范围内，正交与非正交算符集之间未观察到关联能的显著差异，表明非正交性不会损害精度。