[論文レビュー] Physical insights from imaginary-time density--density correlation functions

Tobias Dornheim, Zhandos A. Moldabekov+3 more|arXiv (Cornell University)|Sep 6, 2022

Advanced Chemical Physics Studies被引用数 5

ひとこと要約

本稿は、虚時間密度-密度相関関数 $ F(\boldsymbol{q},\tau) $ が動的構造因子 $ S(\boldsymbol{q},\tilde{\tau}) $ と同等の物理的情報を含み、解析的続行を経ずに温度や準粒子励起エネルギーといった主要な物理量を直接抽出可能であることを示している。一貫した電子気体の第一原理経路積分モンテカルロデータを用いて、ロトンモードのような非自明な特徴が $ F(\boldsymbol{q},\tau) $ において明確に現れることを示し、$ S(\boldsymbol{q},\tilde{\tau}) $ の周波数モーメントと $ \tau = 0 $ 近傍における $ F(\boldsymbol{q},\tau) $ の曲率との間の直接的な関係を確立した。これにより、実周波数解析に対する安定で強力な代替手法が得られる。

ABSTRACT

The accurate theoretical description of the dynamic properties of correlated quantum many-body systems such as the dynamic structure factor $S(\mathbf{q},ω)$ constitutes an important task in many fields. Unfortunately, highly accurate quantum Monte Carlo methods are usually restricted to the imaginary time domain, and the analytic continuation of the imaginary time density--density correlation function $F(\mathbf{q},τ)$ to real frequencies is a notoriously hard problem. In this work, we argue that no such analytic continuation is required as $F(\mathbf{q},τ)$ contains, by definition, the same physical information as $S(\mathbf{q},ω)$, only in an unfamiliar representation. Specifically, we show how we can directly extract key information such as the temperature or quasi-particle excitation energies from the $τ$-domain, which is highly relevant for equation-of-state measurements of matter under extreme conditions. As a practical example, we consider \emph{ab initio} path integral Monte Carlo results for the uniform electron gas (UEG), and demonstrate that even nontrivial processes such as the \emph{roton feature} of the UEG at low density straightforwardly manifest in $F(\mathbf{q},τ)$. In fact, directly working in the $τ$-domain is advantageous for many reasons and holds the enticing promise for unprecedented agreement between theory and experiment.

研究の動機と目的

• 虚時間密度-密度相関関数 $ F(\boldsymbol{q},\tau) $ が動的構造因子 $ S(\boldsymbol{q},\tilde{\tau}) $ と同等の物理的情報を含み、よりアクセスしやすい表現形であることを示すこと。
• 温度や準粒子励起エネルギーといった主要な物理量が、実周波数への解析的続行を必要とせず、$ F(\boldsymbol{q},\tau) $ から直接抽出可能であることを示すこと。
• $ S(\boldsymbol{q},\tilde{\tau}) $ の周波数モーメントと $ \tau = 0 $ における $ F(\boldsymbol{q},\tau) $ の曲率との間の定量的関係を確立し、新たな診断ツールの開発を可能とすること。
• 特に低密度におけるロトン特徴の捉え込みに優れる、一貫した電子気体の第一原理経路積分モンテカルロシミュレーションを用いて、この手法の妥当性を検証すること。
• 状態方程式や温い密度物質の診断において、$ F(\boldsymbol{q},\tau) $ を実周波数心度関数の実用的で正確な代替手段として推奨すること。

提案手法

• 著者らは虚時間のフラクチュエーション-分散定理を導出し、静的密度応答関数 $ \chi(\boldsymbol{q},0) $ が $ F(\boldsymbol{q},\tau) $ の虚時間区間 $ \tau \in [0,\beta] $ における積分と直接関係することを示した。式は $ \chi(\boldsymbol{q},0) = -n \int_0^\beta F(\boldsymbol{q},\tau) d\tau $ で与えられる。
• $ S(\boldsymbol{q},\omega) $ の周波数モーメントと $ \tau = 0 $ の周りにおける $ F(\boldsymbol{q},\tau) $ のテイラー展開係数、特に2階微分（曲率）との間の数学的関係を確立した。この曲率は第2周波数モーメントを符号化している。
• 有限温度および密度の変化する一様電子気体の $ F(\boldsymbol{q},\tau) $ を生成するために、第一原理経路積分モンテカルロ（PIMC）シミュレーションを用いた。
• 特に $ \tau $-領域における曲率と減衰挙動の関数的性質を分析し、動的構造因子におけるロトン最小値のような物理的特徴を推定した。
• 異なる波数ベクトルおよび温度における $ F(\boldsymbol{q},\tau) $ の振る舞いを比較し、ロトンのような非自明な特徴が $ \tau $-ドメインにおいて直接可視化されることを示した。
• 虚時間表現で直接作業することにより、解析的続行問題の数値的不安定性を回避し、量子モンテカルロ法が虚時間において安定かつ高精度に機能することを活用した。

実験結果

リサーチクエスチョン

• RQ1 解析的続行を経ずに、虚時間密度-密度相関関数 $ F(\boldsymbol{q},\tau) $ から温度や準粒子励起エネルギーといった物理的情報を直接抽出可能か？
• RQ2 動的構造因子 $ S(\boldsymbol{q},\omega) $ の周波数モーメントと、$ \tau = 0 $ における $ F(\boldsymbol{q},\tau) $ の曲率との間にはどのような関係があるか？
• RQ3 一様電子気体におけるロトンモードのような非自明な特徴が、実周波数に変換しないまま $ F(\boldsymbol{q},\tau) $ において明確に識別可能か？
• RQ4 $ F(\boldsymbol{q},\tau) $ が温い密度物質や状態方程式測定の診断において、$ S(\boldsymbol{q},\omega) $ と同等に信頼的かつ実用的な代替手段として用いられる可能性はどの程度か？
• RQ5 高次虚時間相関関数は、非線形応答関数や3体・4体相関を捉えるために果たす役割は何か？

主な発見

• 虚時間密度-密度相関関数 $ F(\boldsymbol{q},\tau) $ は動的構造因子 $ S(\boldsymbol{q},\omega) $ と同等の物理的情報を含み、多くの物理的洞察を得るには解析的続行が不要である。
• $ \tau = 0 $ における $ F(\boldsymbol{q},\tau) $ の曲率は、$ S(\boldsymbol{q},\omega) $ の第2周波数モーメントと直接関係しており、$ \tau $-ドメインの挙動と実周波数スペクトルモーメントとの間の定量的リンクを提供する。
• 低密度における一様電子気体のロトン特徴が、$ F(\boldsymbol{q},\tau) $ において明確に可視化されることを示した。非自明なスペクトル特徴が $ \tau $-表現そのものに直接現れることを示した。
• 静的密度応答関数 $ \chi(\boldsymbol{q},0) $ は、$ \tau \in [0,\beta] $ における $ F(\boldsymbol{q},\tau) $ の積分から直接計算可能であり、虚時間のフラクチュエーション-分散定理の妥当性が裏付けられた。
• 虚時間ドメインで直接作業することで、解析的続行の数値的不安定性を回避でき、電子-電子静的構造因子のような周波数積分済み物理量の安定的抽出が可能となった。
• PIMCシミュレーションで入手可能な高次虚時間相関関数を用いることで、非線形応答関数や多体相関を調べる手段が得られ、量子多体系における診断の新たな道筋が開かれた。