[論文レビュー] Inversion symmetry of singular values and a new orbital ordering method in tensor train approximations for quantum chemistry

Mi-Song Dupuy, Gero Friesecke|arXiv (Cornell University)|Feb 18, 2020

Tensor decomposition and applications参考文献 28被引用数 6

ひとこと要約

本稿では、非相互作用ハミルトニアンの基底状態における特異値の予期せぬ反転対称性を活用して、量子化学におけるテンソル積（TT）近似のための新しい軌道順序付け手法を提案する。この対称性から導かれる隠れた不変量を最適化することで、特異値の尾部減衰を4–5桁改善し、Fiedler順序付けなどの既存手法を著しく上回る。

ABSTRACT

The tensor train approximation of electronic wave functions lies at the core of the QC-DMRG (Quantum Chemistry Density Matrix Renormalization Group) method, a recent state-of-the-art method for numerically solving the $N$-electron Schr\\"odinger equation. It is well known that the accuracy of TT approximations is governed by the tail of the associated singular values, which in turn strongly depends on the ordering of the one-body basis. Here we find that the singular values $s_1\\ge s_2\\ge ... \\ge s_d$ of tensors representing ground states of noninteracting Hamiltonians possess a surprising inversion symmetry, $s_1s_d=s_2s_{d-1}=s_3s_{d-2}=...$, thus reducing the tail behaviour to a single hidden invariant, which moreover depends explicitly on the ordering of the basis. For correlated wavefunctions, we find that the tail is upper bounded by a suitable superposition of the invariants. Optimizing the invariants or their superposition thus provides a new ordering scheme for QC-DMRG. Numerical tests on simple examples, i.e. linear combinations of a few Slater determinants, show that the new scheme reduces the tail of the singular values by several orders of magnitudes over existing methods, including the widely used Fiedler order.

研究の動機と目的

• 量子化学におけるテンソル積（TT）近似のための基底順序付けという重要な課題に取り組み、QC-DMRG手法の精度を決定づける要因を解明する。
• 非相互作用ハミルトニアンにおけるTT近似の特異値の振る舞いを調査し、特に軌道順序付けに依存する特異値の尾部減衰の依存関係を明らかにする。
• 相関波動関数の特異値尾部減衰を最小化する理論的根拠に基づいた最適化された軌道順序付けスキームを開発する。
• Fiedler順序付けなどの既存手法を上回る特異値減衰とTT近似品質を達成する、実用的で数値的に検証可能な手法を提供する。

提案手法

• 非相互作用基底状態を表すテンソルの特異値に現れる新しい反転対称性を同定する：$ s_1 s_d = s_2 s_{d-1} = \text{const} $、ここで $ s_1 \to s_d $ は降順に並べた特異値である。
• 特異値の尾部挙動を完全に特徴付ける隠れた不変量 $ I = s_1 s_d $ を定義する。
• この不変量を最小化することで特異値減衰を抑制する「最良係数順序付け」——軌道再順序付け戦略を提案する。
• 個々のスレーター行列式からの不変量の重ね合わせを用いて、相関波動関数への拡張を図り、尾部減衰の上界を導出する。
• 典型的な量子化学状態を模倣する、少数のスレーター行列式の線形結合に対して、手法を実装・テストする。
• 数値実験を通じて、新手法がFiedler順序付けや他の標準的手法と比較して特異値減衰の観点でどのように優れているかを評価する。

実験結果

リサーチクエスチョン

• RQ1 非相互作用系における1体軌道の順序付けが、テンソル積近似における特異値の尾部減衰に与える依存関係は何か？
• RQ2 非相互作用基底状態の特異値スペクトルには、近似の改善に利用可能な隠れた対称性が存在するか？
• RQ3 スレーター行列式成分から導かれる不変量を用いて、相関波動関数の特異値の尾部挙動を上界で抑えられるか？
• RQ4 単一の不変量を最小化する軌道順序付けスキームは、Fiedler順序付けなどの既存手法と比較して、顕著に優れたTT近似品質を達成できるか？

主な発見

• 非相互作用基底状態（スレーター行列式）のテンソル積表現の特異値は、正確な反転対称性 $ s_1 s_d = s_2 s_{d-1} = \text{const} $ を示す。この対称性は、軌道順序付けに明示的に依存する。
• 隠れた不変量 $ I = s_1 s_d $ が特異値の尾部挙動を完全に決定づけ、尾部制御の複雑さを1つのスカラー量に簡略化する。
• 「最良係数順序付け」による不変量最適化により、標準的手法と比較して特異値の尾部減衰が4–5桁改善される。
• 本手法はFiedler順序付けを著しく上回り、特異値減衰の改善は1桁程度にとどまる。
• 相関波動関数においては、個々のスレーター行列式からの不変量の重ね合わせによって尾部減衰が上界で抑えられ、一般化された順序付けスキームが可能になる。
• 少数のスレーター行列式の線形結合に対する数値的検証により、本手法がテストケース全体で一貫して優れた特異値減衰を達成することが確認された。