[論文レビュー] 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）法に対して、線形独立かつスピン完全な空間置換演算子を生成する、きめの細かい4段階のアルゴリズムを提示する。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 法におけるスピン不完全性が引き起こす偽のスピン混合と相関エネルギーの誤差を解消すること。
• 10次励起までおよび多重度 2S+1=11 まで、空間置換演算子の線形独立性とスピン完全性を保証すること。
• スピン軌道形式や切断された CSF 多様体に依存せず、任意の高スピン基底状態に一般化可能で自動化可能なフレームワークを提供すること。

提案手法

• スピン固有状態を生成しスピン完全性を保証するために、Löwdin の射影演算子技法を採用する。
• 基底 CSF に作用させたときに異なる空間関数を生成する、指定されたランクのすべてのプロトタイプ空間置換演算子（ˆE）を体系的に生成する。
• スピン粒子から空間軌道への割り当てをトポロジー的マッピングにより行い、スピン対称性を保持する。
• 占有、仮想、活性の軌道空間に対して標準化された順序付けを適用し、演算子定義の一貫性を確保し、重複を回避する。
• モジュラーかつブラックボックス型の入出力システムとしての手続きを採用：与えられた置換ランク m に対して、線形独立かつスピン完全な ˆE 演算子の集合を出力する。
• CSF 表現の行列ランク解析による妥当性評価と、Weyl-Robinson 維度公式による確認により、手法の正当性を裏付けた。

実験結果

リサーチクエスチョン

• RQ1 体系的かつ一般化可能な手法が、任意の高スピン開殻基底状態および励起ランクに対して、線形独立かつスピン完全な空間置換演算子を生成できるか？
• RQ2 スピン軌道形式や切断された CSF 多様体に依存せずに、カップルドクラスター理論におけるスピン完全性を厳密に強制する方法は何か？
• RQ3 スピン不完全性が開殻 CC 計算における相関エネルギーに与える影響は何か？ そして、その影響は定量的に評価可能か？
• RQ4 スピン適応的ではあるがスピン不完全な演算子を使用した場合、完全にスピン完全な実装と比較して顕著な誤差が生じるか？

主な発見

• 提案手法は、10次励起までおよび多重度 2S+1=11 まで、線形独立かつスピン完全な空間置換演算子を成功裏に生成し、行列ランクと次元公式による確認が得られた。
• 生成された演算子を用いた概念実証用 CC 実装は、ボロン原子の高スピン状態に成功裏に適用され、スピン射影誤差がゼロの純粋なスピン状態が得られた。
• スピン適応的かつスピン完全な CC 法による相関エネルギーは、スピン軌道 CC と比較して僅かな差異（四重項状態では約 0.75%）に留まり、高い精度を示した。
• 直接比較において、スピン適応的ではあるがスピン不完全な CC は、相関エネルギーに 0.4–0.05 mH の持続的誤差を示し、スピン完全性の重要性を裏付けた。
• 二重精度の範囲内で、直交および非直交な演算子集合間に相関エネルギーに顕著な差異は認められず、非直交性が精度を損なわないことを示唆した。