[论文解读] Dynamical and thermodynamical stability of isothermal distributions in the HMF model

Pierre-Henri Chavanis|arXiv (Cornell University)|Jul 28, 2010

Statistical Mechanics and Entropy参考文献 4被引用 5

一句话总结

本文基于变分原理，推导了哈密顿平均场（HMF）模型中等温分布的动力学与热力学稳定性条件。分析严格表明：能量高于临界值 $E_c$ 时，均匀相是熵的最大值（稳定）；非均匀相始终为稳定的熵最大值，且在正则系综与微正则系综中磁化率方差不同，证实了相变附近系综非等价性。

ABSTRACT

We provide a new derivation of the conditions of dynamical and thermodynamical stability of homogeneous and inhomogeneous isothermal distributions in the Hamiltonian Mean Field (HMF) model. This proof completes the original thermodynamical approach of Inagaki [Prog. Theor. Phys. 90, 557 (1993)]. Our formalism, based on variational principles, is simple and the method can be applied to more general situations. For example, it can be used to settle the dynamical stability of polytropic distributions with respect to the Vlasov equation [Chavanis & Campa, arXiv:1001.2109]. For isothermal distributions, the calculations can be performed fully analytically, providing therefore a clear illustration of the method.

研究动机与目标

• 推导HMF模型中等温分布的动力学与热力学稳定性的严格条件。
• 通过适用于更广泛系统的变分形式，扩展并完善Igaki的热力学方法。
• 通过在固定能量与质量下进行熵最大化，阐明均匀相与非均匀相的稳定性。
• 证明非均匀相始终为稳定的熵最大值，而均匀相在 $E_c$ 以下变得不稳定。
• 量化正则系综与微正则系综中磁化率方差的差异，展示系综非等价性。

提出的方法

• 在固定质量 $M$ 与能量 $E$ 的约束下，将变分原理应用于熵泛函 $S[f]$，推导稳定性条件，推广了Poincaré关于平衡解线性系列的研究方法。
• 通过计算熵的二阶变分 $S''(B)$ 确定稳定性：$S''(B) < 0$ 表示局部熵最大值（稳定），$S''(B) > 0$ 表示最小值（不稳定）。
• 利用关系式 $S''(B) = -\beta F''(B)$ 将热力学稳定性与自由能曲率关联，实现解析计算。
• 在 $N \to \infty$ 极限下，应用鞍点近似与Hubbard-Stratonovich变换，计算配分函数与热容曲线。
• 推导出两种系综中磁化率的显式表达式：$\langle B^2 \rangle_{MCE} = -1/S''(B)$ 与 $\langle B^2 \rangle_{CE} = -1/F''(B)$。
• 在临界点（$\epsilon \to 1^-$，$\eta \to 1^+$）与基态（$\epsilon \to -2^+$）附近进行渐近展开，分析涨落发散行为。

实验结果

研究问题

• RQ1 HMF模型中等温分布的动力学与热力学稳定性的充要条件是什么？
• RQ2 熵的二阶导数 $S''(B)$ 如何决定某一相（均匀或非均匀）是稳定最大值还是不稳定最小值？
• RQ3 为何在非均匀相中，正则系综与微正则系综的磁化率方差不同？这对系综等价性有何含义？
• RQ4 临界能量 $E_c$ 与临界温度 $T_c$ 如何标记均匀相由稳定转为不稳定的相变点？
• RQ5 相变附近与基态极限下磁化率涨落的定量行为如何？它们如何随系统参数变化？

主要发现

• 当 $E > E_c$ 时，均匀相是熵的最大值（稳定），由 $S''(0) = -\frac{8\pi^2}{k^2M}(1 - \frac{E_c}{E}) < 0$ 严格证明；当 $E < E_c$ 时，该相变得不稳定。
• 非均匀相（$B \neq 0$）始终为稳定的熵最大值，因为 $S''(B) < 0$ 对所有非均匀区域内的 $B$ 成立。
• 在均匀相中，微正则系综的磁化率方差为 $\langle B^2 \rangle_{MCE} = \frac{k^2M}{4\pi^2} \frac{1}{1 - E_c/E}$，适用于 $E > E_c$。
• 在非均匀相中，正则系综与微正则系综的方差不同：$\langle ( abla B)^2 \rangle_{CE} \neq \langle ( abla B)^2 \rangle_{MCE}$，证实了系综非等价性。
• 在临界点附近（$\epsilon \to 1^-$），正则系综方差标度为 $\sim \frac{5}{4(1 - \epsilon)}$，微正则系综为 $\sim \frac{1}{4(1 - \epsilon)}$，表明涨落发散但不等价。
• 在基态极限（$\epsilon \to -2^+$）下，两种方差分别标度为 $\sim \frac{( abla + 2)^2}{32}$ 与 $\sim \frac{( abla + 2)^2}{64}$，表明低能量下涨落被强烈抑制。