[論文レビュー] Thermodynamic coherence of the Variational Average-Atom in Quantum Plasmas (VAAQP) approach

R. Piron, T. Błeński+1 more|arXiv (Cornell University)|Feb 17, 2009

High-pressure geophysics and materials参考文献 14被引用数 7

ひとこと要約

本稿では、ヴァイナル・ウィルス・スフィアの電荷中性を仮定せず、Virialの定理を満たすことで熱力学的整合性を確保する、完全変分的平均原子モデルであるVAAQPコードを提示する。電子圧力の簡潔な解析的表現を数値微分を用いず導出しており、非変分的モデルやInferno型モデルと比較して、アルミニウムの10 eV等温線において一貫した平均イオン化状態と圧力が得られることを示している。

ABSTRACT

A new code called VAAQP (Variational Average-Atom in Quantum Plasmas) is reported. The model as well as main results of previous studies are briefly recalled. The code is based on a new fully variational model of dense plasmas at equilibrium with quantum treatment of all electrons. The code can calculate the Average Atom structure and the mean ionization from the variational equations respecting the virial theorem and without imposing the neutrality of the Wigner-Seitz sphere. The formula obtained for the electronic pressure is simple and does not require any numerical differentiation. A description of the principal features of the code is given. The thermodynamic consistency of the results obtained with VAAQP is shown by a comparison with another approach on the example of the aluminium 10 eV isotherm EOS curve. A first comparison to an INFERNO-type model is also presented.

研究の動機と目的

• 熱力学的整合性を持つ量子平均原子モデルを、密度の高いプラズマに適用し、Virialの定理を満たすように開発すること。
• 非変分的モデル（Inferno や NWS）の限界を是正すること。これらは人工的な中性化条件と自由エネルギーの数値微分に依存しており、圧力計算に問題を生じる。
• Wigner-Seitz球の中性化を仮定せず、平均イオン化状態と電子圧力を正確に計算可能にする。
• Wigner-Seitz半径を超えたFriedel振動を捉える変分的フレームワークを提供し、温い密度物質領域での精度を向上させること。
• Inferno や非変分的 NWS 手法と比較可能な基準となる、熱力学的整合性を持つ結果を確立すること。

提案手法

• 熱力学的平衡状態にある密度の高いプラズマに対して、完全変分的モデルを構築し、すべての電子を量子力学的に取り扱う。
• 変分原理を用いて、電子密度およびポテンシャルの自己無撞着方程式を導出し、Virialの定理を満たすようにする。
• Wigner-Seitz球の中性化仮定を排除し、平均イオン化状態が変分解から自然に導かれるようにする。
• 自由エネルギーの数値微分を避けるために、電子圧力の閉形式解析的表現（式17）を導出する。
• VAAQPコードを実装し、アルミニウムの10 eVにおける平均原子構造、平均イオン化状態、状態方程式（EOS）を計算する。
• 同じ熱力学的条件下で、非変分的NWSおよびInferno型モデルと比較し、一貫性と正確性を評価する。

実験結果

リサーチクエスチョン

• RQ1 完全変分的平均原子モデルは、Wigner-Seitz球の中性化を仮定せず、熱力学的整合性を達成できるか？
• RQ2 変分原理から導かれた解析的圧力式は、非変分的モデルにおける数値微分による圧力と比べてどのように異なるか？
• RQ3 Friedel振動がWigner-Seitz半径を超えて及ぼす影響は、温い密度物質領域における電子密度および平均イオン化状態にどの程度寄与するか？
• RQ4 同じ熱力学的状態において、変分的VAAQPモデルと非変分的またはInferno型モデルとの間で、電子密度および平均イオン化状態にどのような差が生じるか？
• RQ5 長波長の電荷振動を捉える変分的アプローチにより、温い密度物質領域で熱力学的整合性を持つ結果を得られるか？

主な発見

• VAAQPモデルは、熱力学的整合性を持つ結果を生み出し、アルミニウムの10 eV等温線に沿って、解析的圧力（式17）と自由エネルギーの数値微分による圧力が完全に一致することを確認した。
• これに対して、非変分的NWSモデルでは、解析的圧力と数値的圧力の間に顕著な不一致が見られ、熱力学的整合性の欠如を示している。
• 電子圧力は類似しているものの、VAAQPとNWSの平均イオン化状態には最大10%の差が生じており、圧力だけではモデルの妥当性を検証できないことを示している。
• Friedel振動領域、特に温い密度物質領域において、VAAQPとNWSの電子密度プロファイルには顕著な差が生じている。
• Inferno型モデルで中性化を仮定しているにもかかわらず、その電子密度は非変分的NWS解よりも変分的VAAQP解に近い値を示しており、予期しない安定化効果が生じている可能性を示唆している。
• 変分的アプローチにより、Inferno や NWS のようなモデルが無視するWigner-Seitz半径を超えた長波長電荷振動を的確に捉えることができ、より正確な原子構造予測が可能になった。