[論文レビュー] Identifying the maximum entropy method as a special limit of stochastic analytic continuation

K. S. D. Beach|arXiv (Cornell University)|Mar 1, 2004

Advanced Thermodynamics and Statistical Mechanics被引用数 81

ひとこと要約

本稿では、最大エントロピー法（MEM）が、スペクトル関数を再構築するために場の配置をサンプリングする動的アプローチである確率的解析接続（SAC）の平均場極限として確立される。データフィットとエントロピーから導かれるハミルトニアンを持つ古典的場系に解析接続を写像することで、著者らはMEMがゼロ温度極限として現れることを示し、SACはMEMの平滑化バイアスを超えて鋭いスペクトル特徴をよりよく解像する熱揺らぎを含むことを明らかにする。

ABSTRACT

The maximum entropy method is shown to be a special limit of the stochastic analytic continuation method introduced by Sandvik [Phys. Rev. B 57, 10287 (1998)]. We employ a mapping between the analytic continuation problem and a system of interacting classical fields. The Hamiltonian of this system is chosen such that the determination of its ground state field configuration corresponds to an unregularized inversion of the analytic continuation input data. The regularization is effected by performing a thermal average over the field configurations at a small fictitious temperature using Monte Carlo sampling. We prove that the maximum entropy method, the currently accepted state of the art, is simply the mean field limit of this fully dynamical procedure. We also describe a technical innovation: we suggest that a parallel tempering algorithm leads to better traversal of the phase space and makes it easy to identify the critical value of the regularization temperature.

研究の動機と目的

• 広く用いられている最大エントロピー法（MEM）を、より一般的な確率的解析接続（SAC）フレームワークと厳密に結びつけること。
• MEMが固有の平滑化バイアスにより鋭いスペクトル特徴を損なうという限界を克服すること。
• 恣意的平均化手法を越えた、体系的かつ統計的に根拠のある解析接続手法を構築すること。
• SACにおけるサンプリング効率の向上と最適正則化温度の特定を目的とした並列温度法アルゴリズムの導入。

提案手法

• 解析接続問題を、場の配置がスペクトル関数に対応する単位区間上の相互作用を有する古典的場系に写像する。
• 基底状態が虚時間データの正則化なし逆問題に対応するハミルトニアンを定義する。
• 小さな擬似温度における場配置の熱平均を用いて正則化を導入し、モンテカルロサンプリングを実行する。
• MEMがSAC手順の平均場（ゼロ温度）極限に一致することを示すことにより、MEMとSACの関係を導出する。
• 位相空間の到達性向上と臨界正則化温度の特定を目的とした並列温度法アルゴリズムを実装する。
• スペクトル関数の広い動的範囲に対応し、数値的安定性を向上させる対数メッシュ離散化を用いる。

実験結果

リサーチクエスチョン

• RQ1 量子モンテカルロデータの解析接続の文脈において、最大エントロピー法（MEM）と確率的解析接続（SAC）とはどのように関係しているか？
• RQ2 確率的解析接続法は、統計的場理論の枠組みから体系的に導出可能か？
• RQ3 SACがMEMの動的一般化として、平滑解へのバイアスを低減する熱揺らぎを含むことで、より優れたスペクトル再構築を達成するか？
• RQ4 正則化温度はSACにおいてどのような役割を果たし、最適に決定できるか？
• RQ5 標準のマルコフ連鎖モンテカルロ法と比較して、並列温度法はSACの収束性と精度を向上させるか？

主な発見

• 最大エントロピー法が、確率的解析接続の平均場極限として厳密に示され、MEMはSAC手順のゼロ温度極限に対応することが明らかになった。
• SACは熱揺らぎを含むため、MEMがそれらを洗練させてしまうのとは対照的に、鋭い特徴をよりよく解像するスペクトルを生成する。
• 並列温度法アルゴリズムの使用により、サンプリング効率が著しく向上し、臨界正則化温度の信頼性のある同定が可能になった。
• 本手法は、恣意的平均化技術の代替として体系的かつベイズ的根拠を持つものであり、古典的場動力学の明確な物理的解釈を有する。
• BCS超伝導体に対する数値的テストでは、SACがMEMよりも正確な固有状態対角化結果に近いスペクトルを生成し、特に超伝導ギャップとコherenceピークの解像度が向上した。
• ラグランジュ乗数を用いた正規化とエントロピー制約の一貫した取り扱いにより、解の安定性と一意性が保証された。