[論文レビュー] Variational determination of the two-particle density matrix as a quantum many-body technique

Brecht Verstichel|arXiv (Cornell University)|Mar 26, 2012

Quantum many-body systems参考文献 79被引用数 4

ひとこと要約

本博士論文では、半定値計画法を用いて、量子多体系における2粒子縮約密度行列（2RDM）を直接的に決定する変分的手法を提示する。N-実現可能性を満たす密度行列に対する制約付き最適化問題として定式化し、対称性およびスピン適応型リフト条件を組み込むことで、標準的な2インデックス制約が失敗する1次元 Hubbard モデルにおける強い電子相関効果を正確に捉えることができる。

ABSTRACT

In this thesis the variational optimisation of the density matrix is discussed as a method in many-body quantum mechanics. This is a relatively unknown technique in which one tries to obtain the two-particle reduced density matrix directly in a variational approach. In the first Chapter we introduce the subject in the broader context of many-body quantum mechanics, and briefly sketch the history of the field. The second Chapter tries to summarise what is known about N-representability of reduced density matrices, and derives these results in a unified framework. The optimisation problem can be formulated as a semidefinite program. Three different algorithms are discussed in Chapter three, and their performance is compared. In Chapter four we show how symmetry can be exploited to reduce the computational cost considerably. Several applications of the method are discussed in Chapter five. We show that the standard two-index constraints fail in the strong-correlation limit of the one-dimensional Hubbard model. In Chapter six we identify the spin-adapted lifting conditions as the most compact constraints that can correctly describe this limit.

研究の動機と目的

• 多体量子系における2粒子縮約密度行列（2RDM）を直接的に決定するための変分的手法の開発。
• 縮約密度行列理論におけるN-実現可能性の課題に取り組み、統一的な枠組み内で制約を導出する。
• 特に強い相関極限における1次元 Hubbard モデルを対象とし、強相関系における精度の向上。
• 強相関を正確に記述できる、具体的にはスピン適応型リフト条件である、コンactかつ物理的に意味のある制約の特定と実装。
• 2RDMの変分最適化における対称性の活用による計算コストの低減。

提案手法

• 2RDMの変分最適化を半定値計画問題（SDP）として定式化し、数値的効率な解法を可能にする。
• 2RDMが有効な多体量子状態に対応することを保証するため、N-実現可能性条件を導出し適用する。
• SDPを解くための3つの異なるアルゴリズムを採用し、収束性およびスケーラビリティの観点から性能を比較する。
• 最適化問題の次元を低減し、計算コストを削減するための対称性活用技術を導入する。
• 強相関をコンactかつ物理的に正確に記述できる制約として、スピン適応型リフト条件を用いる。
• 1次元 Hubbard モデルへの応用を通じて手法の妥当性を検証し、正確な解と比較する。

実験結果

リサーチクエスチョン

• RQ1 強相関を示す多体量子系において、2粒子密度行列に対する変分的手法は、従来の手法を上回る性能を発揮できるか？
• RQ2 1次元 Hubbard モデルの強い相関極限において、標準的な2インデックスN-実現可能性条件がなぜ失敗するのか？
• RQ3 2RDM形式において、強い電子相関を正確に記述できる最小かつ物理的に意味のある制約の集合は何か？
• RQ4 2RDM最適化の計算コストを低減するために、どのようにして対称性を体系的に活用できるか？
• RQ5 スピン適応型リフト条件は、標準的な制約よりもよりコンactかつ正確なN-実現可能性の表現を提供できるか？

主な発見

• 標準的な2インデックスN-実現可能性制約は、1次元 Hubbard モデルの強い相関極限を正しく記述できず、物理的に意味のない結果を生じる。
• スピン適応型リフト条件は、2RDM形式において強い相関効果を正確に捉えることのできる、最もコンactかつ正確な制約集合であると特定された。
• 対称性低減を組み合わせた変分的2RDM手法は、精度を損なわずに顕著な計算コスト削減を達成した。
• 3つのSDPを解くためのアルゴリズムをベンチマーク化した結果、収束速度および数値的安定性に差が認められた。
• 本手法は、強い相関領域におけるHubbardモデルの主要な物理的特徴を正確に再現でき、強相関系への応用の有効性を裏付けた。
• N-実現可能性条件の統一的枠組みにより、さまざまな物理系に適用可能な一貫性のある制約の導出が可能となった。