[论文解读] 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

一句话总结

本文提出了VAAQP代码，这是一种用于致密量子等离子体的全变分平均原子模型，通过满足virial定理并消除对Wigner-Seitz球体电中性性的假设，确保了热力学一致性。该模型推导出电子压强的简单解析表达式，无需数值微分，与非变分模型及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 在温稠密物质中，超出Wigner-Seitz半径的Friedel振荡对电子密度和平均电离度的影响程度如何？
• RQ4 在相同热力学状态下，变分VAAQP模型与非变分或Inferno型模型在电子密度和平均电离度上的差异有多大？
• RQ5 是否能够通过采用能捕捉长波长电荷振荡的变分方法，在温稠密物质区域实现热力学一致的结果？

主要发现

• VAAQP模型产生了热力学一致的结果，这通过铝10 eV等温线上解析压强（公式17）与数值微分自由能结果的完全一致得到了证实。
• 相比之下，非变分NWS模型在解析压强与数值压强之间表现出显著差异，表明其存在热力学不一致性。
• 尽管电子压强相似，VAAQP与NWS的平均电离度差异可达10%，这表明仅靠压强不足以验证模型的可靠性。
• 在Friedel振荡区域，VAAQP与NWS的电子密度分布存在显著差异，尤其是在温稠密物质区域。
• 尽管Inferno型模型假设电中性，其电子密度分布却比非变分NWS模型更接近变分VAAQP解，表明其存在非预期的稳定化效应。
• 变分方法成功捕捉了超出Wigner-Seitz半径的长波长电荷振荡，而这些效应在Inferno和NWS等模型中被忽略，从而实现了更精确的原子结构预测。