Point-group filtering in unitary coupled-cluster ansätze: exact criticality under the Abelian filter, freeness measured near equilibrium under the full group

技術要約：ユニタリ結合クラスター・アンサンブルにおける点群フィルタリング

問題提起
本論文は、分子の点群対称性によるユニタリ結合クラスター（UCC）アンサンブルのフィルタリングにおける、変分コストおよび理論的正当性を、特に既約表現（irreps）が縮退している非アーベル群の文脈で扱っている。対称性フィルタリングはアーベル群においては標準的な手法であるが、非アーベル群（例：アンモニアの$C_{3v}$、メタンの$T_d$）への適用は、近年の文献において「未解決」の問題となっている。既存の研究では相反する結果が報告されており、一部の研究はフィルタリングされたアンサンブルとフィルタリングされていないアンサンブルの間のエネルギー差は無視できると主張する一方で、他の研究では有意な誤差（最大40 mHartree）を報告している。さらに、群が次元1より大きい既約表現を持つ場合のフィルタリングの変分コストに関する確定した理論的言明や、標準的なアーベル部分群アプローチに対する全点群によるパラメータ圧縮の明確な測定値も存在しない。

手法
著者らは、表現論、変分微積分、および数値的裏付けを組み合わせた厳密な数学的枠組みを用いている。

1. 理論的枠組み:

   • 定義: 論文では2種類のフィルタを定義している。$G$-等変プール（$P_G$）は、全点群$G$に対して不変な演算子を含む。アーベル・プール（$P_{Sym}$）は、極大アーベル部分群$N \subseteq G$に対して不変な演算子を含む。
   • 仮説: 証明は3つの仮説に基づいている。(A1) ハミルトニアンは群と可換である。(A2) 参照決定式は群に対して不変である（閉殻適応軌道）。(A3) パラメータのタイイング（係数の結合）はパリティに対して一貫している。
   • 臨界部分多様体: 著者らは、対称性に適応したパラメータ空間の多様体が、エネルギーの臨界部分多様体であることを証明している。これは、エネルギーの勾配が、フィルタによって除去されたすべての方向において、単に参照決定式においてだけでなく、その部分多様体上のあらゆる点において消失することを意味する。
   • スパンの包含: 鍵となる理論的結果は、完全に不変な演算子の線形スパン（$P_G$）が、アーベル・フィルタリングされた演算子のスパン（$P_{Sym}$）に含まれることを確立している。その結果、アーベル・フィルタリングされたアンサンブルによって到達可能な状態の集合は、全点群による不変アンサンブルによって到達可能なすべての状態を含む。

2. 数値的裏付け:

   • 古典的CCSD: 6つの分子（HF, LiH, BeH$_2$, H$_2$O, NH$_3$, CH$_4$）に対し、STO-3Gおよび6-31G基底関数を用いて計算を行った。著者らは、標準的なパッケージ（例：PySCF）が最大のアーベル部分群をデフォルトとするため、全群対称性を強制するためのカスタムプロジェクターを実装した。
   • ユニタリ・アンサンブル: STO-3GにおいてUCCSDの計算を行った。著者らは、既約表現のラベルを固定したまま縮退軌道を回転させることで、軌道基底の不整合（misalignment）の影響を分析し、フィルタの「自由度」をテストした。
   • 幾何学的スキャン: アンモニアを平衡構造から静的相関領域まで引き伸ばし、全群フィルタが劣化する幾何学的境界を特定した。

主要な貢献

1. 厳密な臨界性（定理 4.3）: 著者らは、既約表現の次元に関わらず、対称適応部分多様体がエネルギーの臨界部分多様体であることを証明した。除去されたすべての方向における勾配は、この部分多様体上のあらゆる場所で消失する。

   • 帰結: 参照決定式で初期化され、厳密な勾配を持つ準ニュートン法を用いた場合、アンサンブルが全群によってフィルタリングされているかアーベル部分群によってフィルタリングされているかにかかわらず、最適化は全く同じ軌跡を辿り、全く同じエネルギーを返す。これまでに報告されているエネルギーの差は、実装エラー（例：非解析的な勾配や基底の不整合）の証拠であり、フィルタの変分コストではない。

2. スパンの包含（定理 5.3 および 系数 5.4）: 本論文は、全群$G$に対して不変なすべての演算子が、アーベル・フィルタ$N$によって保持される演算子のスパン内に存在することを証明している。

   • 帰結: アーベル・フィルタリングされたプールは、あらゆる幾何構造において、全群不変なプールが到達可能なあらゆる状態に到達するために十分である。これは、アーベル・フィルタリングされた集合におけるエネルギーの下限が、全群不変な集合における下限以下であることを意味する（$\inf_{RSym} E \le \inf_{RG} E$）。論文内では、どちらの集合も厳密な基底状態を含むとは主張していないことを明記しており、両者の下限の等価性は証明されていない。この結果は、不変な生成子に対する十分性を確立するものであり、無制限のアンサンブルに対する十分性を証明するものではない。

3. 軌道基底の証明（系数 6.6）: 著者らは診断ツールを導出した。適応基底においては、ハミルトニアン情報に基づくプールは、対称性フィルタリングされたプールよりも多くのパラメータを含むことはない。もし公開された研究が、ハミルトニアン情報に基づくアンサンブルの方がパラメータ数が多いと報告しているならば、それは軌道基底が対称適応されていないことを証明している。これは、パラメータ数が不一致を示していたHeら[9]の先行研究における不一致を説明するものである。

4. 定量的圧縮と限界:

   • 圧縮: 全点群によるフィルタリングは、アーベル・フィルタと比較してパラメータ数を大幅に削減する（例：NH$_3$: 75 $\to$ 30、CH$_4$: 65 $\to$ 21）。
   • エネルギーコスト: 古典的CCSDの極限において、全群フィルタは6-31G基底において相関エネルギーを最大で$3.4 \times 10^{-12}$ mHartreeしか変化させない。
   • 幾何学的境界: 全群フィルタの「自由度」は平衡付近では保持されるが、静的相関領域では劣化する。引き伸ばされたアンモニアにおいて、全群フィルタによる誤差は、1.60 Åから1.80 Åの間で化学精度（1.6 mHartree）の閾値を越えるが、アーベル・フィルタは全スキャンを通じて自由な状態を維持する。

結果と意義
本論文は、参照が閉殻であり、かつ軌道が対称適応されている限り、全点群によるフィルタリングの「変分コスト」はゼロであると主張している。すなわち、アーベル・フィルタ（現在の標準的な手法）は、全群不変アンサンブルと同じエネルギー下限に到達するために理論的に十分である。先行文献における非アーベルケースでの失敗は、以下の要因に起因するとされている：

1. 基底の不整合: ある計算からの既約表現ラベルを別の計算の軌道に使用すること（例：異なるSCF収束パスや線形代数ライブラリのスレッド数によるもの）により、対称性の仮定が崩れる。
2. 実装上のアーティファクト: 最適化プロトコルにおける非解析的な勾配、または不完全な励起プール。

本研究の意義は、UCCアンサンブルにおける全点群対称性の使用に対する厳密な変分的正当性を提示した点にある。アーベル・フィルタは、参照が閉殻で軌道が対称適応されている場合、不変な生成子の到達可能な状態空間を犠牲にすることなく、探索空間（次元）を削減する「安全な」下限であることを証明した。また、全群による圧縮は、不変な生成子の到達可能な状態空間を損なうことなく、最適化の探索空間を削減するものである。著者らは、これらの結果は閉殻参照にのみ適用され、開殻または軌道縮退のある基底状態には拡張されないこと、および勾配消失の性質が保持されない可能性があること、そして制限された集合と無制限の集合の等価性は証明されていないことを謙虚に付記している。