[論文レビュー] Maximum Entropy and the Variational Method in Statistical Mechanics: an Application to Simple Fluids

Chih‐Yuan Tseng, Ariel Caticha|ArXiv.org|Nov 25, 2004

Statistical Mechanics and Entropy参考文献 28被引用数 3

ひとこと要約

本稿は、統計力学における複雑な確率分布の最適近似を生成するための最大エントロピー（ME）フレームワークを導入する。特に単純な流体に対して有効である。相対エントロピーを用いて、硬球の直径を変化させた試行分布の族の中から最良の分布を選択することで、従来の変分法を改善し、Lennard-Jones流体に対してはBHおよびWCAといった最先端の摂動理論と同等の結果を達成する。

ABSTRACT

We develop the method of Maximum Entropy (ME) as a technique to generate approximations to probability distributions. The central results consist in (a) justifying the use of relative entropy as the uniquely natural criterion to select a "best" approximation from within a family of trial distributions, and (b) to quantify the extent to which non-optimal trial distributions are ruled out. The Bogoliuvob variational method is shown to be included as a special case. As an illustration we apply our method to simple fluids. In a first use of the ME method the "exact" canonical distribution is approximated by that of a fluid of hard spheres and ME is used to select the optimal value of the hard-sphere diameter. A second, more refined application of the ME method approximates the "exact" distribution by a suitably weighed average over different hard-sphere diameters and leads to a considerable improvement in accounting for the soft-core nature of the interatomic potential. As a specific example, the radial distribution function and the equation of state for a Lennard-Jones fluid (Argon) are compared with results from molecular dynamics simulations.

研究の動機と目的

• 統計力学における不確実な確率分布近似を生成する一般化された手法を開発すること。
• 相対エントロピーが試行分布族から最良の近似を選ぶために、一意的かつ自然な基準であることを正当化すること。
• 単一の分布近似にとどまらず、試行分布の重み付き平均を組み込むことで変分法を拡張すること。
• 古典的単純流体、特にLennard-Jones流体に対してMEフレームワークをテストし、精度と頑健性のベンチマークを提供すること。
• MEに基づく近似が、BHおよびWCAといった既存の摂動理論と同等またはそれを上回る、熱力学的および構造的性質の予測を可能にすることを示すこと。

提案手法

• 相対エントロピーで定量化される情報損失を最小化しながら、新たな制約を満たすように事前分布を事後分布に更新する最大エントロピー（ME）の原則を用いる。
• 異なる直径を持つ硬球試行分布の族を用いて、流体の正確な正準分布を近似するME法を適用する。
• 異なる硬球直径における重み分布をMEの原則に従って最適化し、柔ららかくコアの相互作用をよりよく捉える重み付き平均分布を得る。
• ME最適化分布から径方向分布関数（RDF）および状態方程式を導出し、分子動力学シミュレーションおよび摂動理論と比較する。
• 単一の硬球直径のみを用いる場合に、Bogoliubov変分原理がME法の特殊な場合であることを示す。
• 相対エントロピーの分析を通じて、非最適な分布がどれほど排除されるかを評価し、ME法が標準的変分法をどの程度改善するかを定量化する。

実験結果

リサーチクエスチョン

• RQ1 統計力学における扱いにくい確率分布の最適近似を生成するため、最大エントロピー法を体系的に適用する方法は何か？
• RQ2 なぜ相対エントロピーが試行分布族から最良の近似を選ぶために、一意的かつ自然な基準となるのか？
• RQ3 試行分布の重み付き平均を許容することで、単一の分布にとどまらない標準的変分法に比べ、ME法はどの程度改善されるのか？
• RQ4 ME法は、Lennard-Jones流体の径方向分布関数および状態方程式について、既存の摂動理論と同等の精度で予測を可能にするか？
• RQ5 ME法は、分子動力学シミュレーションおよびBHやWCAといった他の先進的理論と比較して、異なる密度および温度における流体の挙動をどの程度的確に捉えられるか？

主な発見

• ME法は、異なる直径を持つ硬球分布の重み付き平均としての確率分布を成功裏に生成し、柔ららかくコアの相互作用を有する原子間ポテンシャルの記述を改善した。
• ME法から導出された径方向分布関数（RDF）は、硬球直径を調整しないままでも、Verletシミュレーション結果の関数的形を非常によく再現した。
• ME法を用いて計算された状態方程式は、単純な変分法よりも明確な改善を示し、特に中間密度領域において分子動力学シミュレーションと良好な一致を示した。
• 低密度領域では、すべての手法（ME、変分法、BH、WCA）がシミュレーションと一致するが、長距離2粒子相関の取り扱いの限界により、乖離が生じる。
• 高密度領域ではWCA摂動理論が最も優れた性能を示すが、ME法も競争力があり、広範な条件下で優れた性能を示した。
• 低温度における圧力-密度等温曲線において、ME法は液体-気体相転移の正しい定性的な挙動を示すヴァン・デル・ワールスのループを生成した。これは、シミュレーションや実験でも観察されるものと一致する。