[論文レビュー] On the pure state $v$-representability of density matrix embedding theory

Fabian M. Faulstich, Raehyun Kim+4 more|arXiv (Cornell University)|Oct 18, 2021

Physics of Superconductivity and Magnetism参考文献 57被引用数 4

ひとこと要約

本稿では、密度行列埋め込み理論（DMET）における非相互作用純状態v-可代表性（NI-PS-V）仮定の緩和を目的として、増大ラグランジュ法（alm-DMET）を提案する。これにより、Aufbau原理が成立しない場合でも、高レベルと低レベルの密度行列ブロックの正確な一致が可能となり、低レベルハミルトニアンの偽のギャップなしの状態を実現する。任意の軌道占有状態を許容しつつも、射影性（占有数が0または1に制限）を維持することで、標準DMETが抱える偽のギャップなし状態を回避し、強い相関系においても高い精度を達成する。

ABSTRACT

Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition and the vanishing low-level gap are related to the assumption that the high-level density matrix blocks are non-interacting pure-state $v$-representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy.

研究の動機と目的

• 非相互作用純状態v-可代表性（NI-PS-V）仮定が現実的でないために、DMETにおける正確な密度行列一致が失敗する問題に対処すること。
• Aufbau原理に違反した場合に生じる低レベルバンドギャップの消失問題を解消すること。
• Aufbau原理に依存しない正確な密度行列ブロック一致を実現する強固な手法の開発。
• 標準DMETが失敗する強い相関系において、数値的安定性と精度を向上させること。
• 2次元 Hubbardモデルおよび水素鎖モデルにおいて、困難な条件下で提案手法の有効性を示すこと。

提案手法

• 正確な密度行列一致を実現するため、相関ポテンシャルのフィッティングを制約付き最適化問題として定式化し、増大ラグランジュ法を用いて制約を強制する。
• 低レベル1-RDMにおける任意の占有状態を許容することで、NI-PS-V条件を緩和するが、同時に射影性（占有数が0または1に制限）を維持する。
• 制約付き最適化を解くために、投影勾配法を用い、ペナルティパラメータとラインサーチによる収束を実現する。
• 収束を制御し、計算コストを削減するため、ハイパーパrameter（例：初期ペナルティ、最大内部反復回数）を導入する。
• 全システムの対角化回数を制限するため、内部反復回数を上限で制御する。
• 収束性を改善するため、10回ごとの外反復でハイパーパrameterを適応的に更新する戦略を導入する。

実験結果

リサーチクエスチョン

• RQ1 高レベル1-RDMブロックがNI-PS-Vでない場合でも、DMETで正確な密度行列一致を達成できるか？
• RQ2 低レベル1-RDMの構築においてAufbau原理を緩和することで、低レベルハミルトニアンにおける偽のギャップなし状態が解消されるか？
• RQ3 増大ラグランジュ法が、強い相関系におけるDMETの数値的ロバスト性と収束性を向上させられるか？
• RQ4 標準の最小二乗法DMETと比較して、提案手法alm-DMETの精度と収束性はどのように異なるか？
• RQ5 ハイパーパrameterの設定が、alm-DMETの収束性と計算コストに与える影響は何か？

主な発見

• 穴ドーピングされた2次元 Hubbardモデル（n = 32/36）のような困難なケースにおいて、Aufbau原理の緩和により、密度行列ブロックの正確な一致が達成された。
• alm-DMET手法は、NI-PS-V条件が満たされない場合に標準DMETが直面する低レベルバンドギャップの消失問題を効果的に回避した。
• 数値結果から、特に強い相関領域において、エネルギーと密度行列一致の両方で、標準の最小二乗法DMETよりも高い精度が得られた。
• 異なる初期化条件を用いたH36鎖系の実験により、alm-DMETの収束が初期RHF推定値に依存しないことが実証された。
• 内部反復回数の制限により、計算コストが削減され、収束性が安定化したが、最適なハイパーパラメータチューニングが性能発揮に不可欠であることが判明した。
• 異なるハイパーパラメータ設定下でも同一のエネルギー軌道が得られたため、エネルギー収束のロバスト性が確認された。