[論文レビュー] The full configuration interaction quantum Monte Carlo method in the lens of inexact power iteration

Jianfeng Lu, Zhe Wang|arXiv (Cornell University)|Nov 24, 2017

Quantum many-body systems参考文献 12被引用数 10

ひとこと要約

本稿は、量子化学における大規模固有値問題を解くために、FCIQMC や FRI といったランダム化アルゴリズムに適用可能な、不正確なパワー反復法の一般化された収束解析フレームワークを提案する。理論的収束速度を確立し、不正確な行列・ベクトル積における分散低減が精度と効率を向上させることを示し、Hubbard モデルや分子における数値実験で iFCIQMC と FRI が決定的手法を上回ることを確認した。

ABSTRACT

In this paper, we propose a general analysis framework for inexact power iteration, which can be used to efficiently solve high dimensional eigenvalue problems arising from quantum many-body problems. Under the proposed framework, we establish the convergence theorems for several recently proposed randomized algorithms, including the full configuration interaction quantum Monte Carlo (FCIQMC) and the fast randomized iteration (FRI). The analysis is consistent with numerical experiments for physical systems such as Hubbard model and small chemical molecules. We also compare the algorithms both in convergence analysis and numerical results.

研究の動機と目的

• 量子多体問題に用いられる不正確なパワー反復法を統一的に分析する理論的枠組みを構築すること。
• 一般化された不正確な行列・ベクトル積の設定下で、FCIQMC や FRI などのランダム化アルゴリズムの収束定理を確立すること。
• 収束速度、圧縮誤差、計算コストの観点から、FCIQMC、FRI、iFCIQMC、HT（ヘルミート＝テイラー）法の性能を比較すること。
• 不正確な行列・ベクトル積の誤差が固有ベクトル近似およびエネルギー推定子の精度に与える影響を定量化すること。
• 不正確な行列・ベクトル積におけるメモリ使用量と分散の主なボトル neck を特定することで、今後のアルゴリズム改善を導くこと。

提案手法

• 行列・ベクトル積を近似写像 $ F_m(A, x_t) $ で置き換える一般化された不正確なパワー反復フレームワークを提案し、確率的または切断された行列・ベクトル積をモデル化する。
• 反復列と真の最大固有ベクトルの間の角度 $ \theta(v_t, u_1) $ を用いて収束を分析し、累積的な不正確さ $ \xi_t $ によって誤差が上限付けられることを示す。
• 固有値比 $ \lambda_2 / \lambda_1 $ に基づく理論的収束速度を導出し、誤差伝搬解析を用いて不正確な設定に拡張する。
• 部分空間収束の指標として、角度の正接 $ \tan \theta(v_t, u_1) $ を用い、これと不正確な行列・ベクトル積の誤差を直接関連付ける。
• Hubbard モデルおよび小規模な分子（Ne、H₂O）に対して、標準基底関数集合を用いた数値実験を実施し、理論的予測の妥当性を検証する。
• 平均二乗誤差（MSE）、平均誤差、標準偏差、反復あたりのウォールクロック時間などの指標を用いてアルゴリズムを比較する。

実験結果

リサーチクエスチョン

• RQ1 行列・ベクトル積における不正確さは、FCIQMC や FRI などのランダム化固有値ソルバーの収束にどのように影響するか？
• RQ2 行列・ベクトル積が確率的または切断的に近似された場合の、不正確なパワー反復の理論的収束速度は何か？
• RQ3 分散低減手法（iFCIQMC や FRI）が実際には決定的手法（HT）よりも高い精度を達成する理由は何か？
• RQ4 不正確な行列・ベクトル積における圧縮誤差は、反復列と真の固有ベクトルとの距離にどのように影響するか？
• RQ5 Hubbard モデルと分子系の両方において FCIQMC と FRI の性能に観察される差異は、どのように説明できるか？

主な発見

• 不正確なパワー反復フレームワークは、FCIQMC、FRI 及びその変種の収束行動を統一的な理論的視点でうまく説明できた。
• Hubbard モデルにおいて、iFCIQMC と FRI は FCIQMC や HT よりも低い平均二乗誤差（MSE）を達成し、iFCIQMC は $ 9.6 \times 10^{-9} $ という最小の MSE を示した。
• aug-cc-pVDZ 基底を用いた Ne 原子において、iFCIQMC は平均誤差 $ 7.7 \times 10^{-5} $ を達成し、反復あたり 1.1 秒の時間を要したが、精度と効率の両面で FCIQMC を上回った。
• cc-pVDZ 基底を用いた H₂O において、iFCIQMC は MSE $ 2.7 \times 10^{-10} $ を達成し、FCIQMC（$ 2.1 \times 10^{-9} $）および HT（$ 1.2 \times 10^{-6} $）よりも顕著に低い値を示した。
• 分子系（Ne および H₂O）において、FRI は FCIQMC と同等の MSE を示したが、Hubbard モデルとは異なり、系依存の性能差が顕在した。
• FRI と HT の時間コストは、FCIQMC や iFCIQMC よりもはるかに高かった。これは、ベクトルがスパースであっても、正確な行列・ベクトル積が高価であるためである。