[论文解读] Thermodynamic Stability at the Two-Particle Level

Alexander Kowalski, M. Reitner+7 more|arXiv (Cornell University)|Sep 20, 2023

Advanced Chemical Physics StudiesPhysics and Astronomy被引用 3

一句话总结

本文通过将强关联费米子系统的热力学稳定性条件重新表述为局域单体和二体关联函数（特别是广义局域电荷感应率的本征值）的形式，而非配分函数的导数，实现了重新表述。关键结果是一个基于感应率最低本征值的单一、紧凑判据：当该本征值低于 $-\beta/t^2$ 时即发生不稳定性，该判据同时捕捉了负压缩性和此前被忽视的正压缩性不稳定性，尤其在莫特转变附近。

ABSTRACT

We show how the stability conditions for a system of interacting fermions that conventionally involve variations of thermodynamic potentials can be rewritten in terms of one- and two-particle correlators. We illustrate the applicability of this alternative formulation in a multi-orbital model of strongly correlated electrons at finite temperatures, inspecting the lowest eigenvalues of the generalized local charge susceptibility in proximity of the phase-separation region. Additionally to the conventional unstable branches, we address unstable solutions possessing a positive, rather than negative compressibility. Our stability conditions require no derivative of free energy functions with conceptual and practical advantages for actual calculations and offer a clear-cut criterion for analyzing the thermodynamics of correlated complex systems.

研究动机与目标

• 解决通过配分函数的二阶导数评估热力学稳定性的概念与计算挑战。
• 提出一种基于局域可观测量的判据，用于诊断不稳定性，而无需计算自由能函数的导数。
• 将稳定性分析扩展至具有正压缩性的体系，这些体系在传统方法中曾被忽视。
• 展示新判据在近洪特金属与莫特绝缘体转变的多轨道 Hubbard 模型中的适用性。

提出的方法

• 通过配分函数 $\Omega$ 对 $\mu$ 和 $U$ 的 Hessian 矩阵推导传统稳定性条件。
• 将这些条件重新表达为单体格林函数 $G$ 以及广义二体感应率 $\chi$ 的本征值与本征矢量的形式。
• 建立 $\chi$ 的最低本征值 $\lambda_I$ 与热力学不稳定性之间的直接联系，临界阈值为 $-\beta/t^2$。
• 通过在双轨道 Hubbard 模型中对广义电荷感应率进行数值对角化，探测相图中各区域的稳定性。
• 分析压缩率 $\kappa$ 随本征值贡献的变化，以识别主导的不稳定性机制。
• 通过与半满时 $\partial D/\partial U$ 的已知发散行为对比，验证该判据，确认其与已有结果的一致性。

实验结果

研究问题

• RQ1 热力学稳定性条件能否在不依赖配分函数二阶导数的前提下重新表述？
• RQ2 广义电荷感应率的最低本征值在表征热力学不稳定性中起什么作用？
• RQ3 是否存在传统判据所遗漏的、具有正压缩性的不稳定解？
• RQ4 感应率的本征值如何与双占据数等物理可观测量的发散相关联？
• RQ5 能否将 $\chi$ 的本征值谱作为关联费米子系统中不稳定性的一个普适指标？

主要发现

• 广义电荷感应率的最低本征值 $\lambda_I$ 作为热力学不稳定的普适指标：当 $\lambda_I < -\beta/t^2$ 时即发生不稳定性。
• 该判据成功识别出负压缩性不稳定性以及此前未被识别的莫特转变附近的正压缩性不稳定性。
• 在半满临界终点附近，$\lambda_I \to -\beta/t^2$ 与 $\partial D/\partial U$ 的发散行为一致，将本征值与已知奇异性行为联系起来。
• 当系统接近不稳定性时，压缩率 $\kappa$ 主要由最低本征值 $\lambda_I$ 的贡献主导。
• 第二小的本征值 $\lambda_{II}$ 始终降低 $\kappa$，但其本身不足以引发不稳定性。
• 该方法避免了显式计算热力学导数，为强关联体系的数值研究提供了在概念和实践上均更优的替代方案。