[论文解读] Thermal response of a Fermi-Pasta-Ulam chain with Andersen thermostats

Federico D’Ambrosio, Marco Baiesi|arXiv (Cornell University)|Jun 21, 2017

Advanced Thermodynamics and Statistical Mechanics参考文献 46被引用 2

一句话总结

本文利用安德森热浴，为惯性非平衡系统推导出一种涨落-响应关系，实现了无需微小扰动即可精确计算热响应率（如比热和热膨胀系数）。该方法在强非平衡区域中优于直接扰动模拟，且在边界热浴作用下的费米-帕斯塔-乌拉姆链中，能获得平滑且低误差的响应率估计。

ABSTRACT

The linear response to temperature variations is well characterised for equilibrium systems but a similar theory is not available, for example, for inertial heat conducting systems, whose paradigm is the Fermi-Pasta-Ulam (FPU) model driven by two different boundary temperatures. For models of inertial systems out of equilibrium, including relaxing systems, we show that Andersen thermostats are a natural tool for studying the thermal response. We derive a fluctuation-response relation that allows to predict thermal expansion coefficients or the heat capacitance in nonequilibrium regimes. Simulations of the FPU chain of oscillators suggest that estimates of susceptibilities obtained with our relation are better than those obtained via a small perturbation.

研究动机与目标

• 为惯性非平衡系统建立线性响应理论，以计算热响应率，其中标准平衡涨落-响应关系失效。
• 解决在存在热流的系统（如受驱的费米-帕斯塔-乌拉姆（FPU）链）中，缺乏可靠框架来计算热响应系数（如比热和热膨胀系数）的问题。
• 证明安德森热浴可作为拉格朗日热浴的数学上可处理且物理解释清晰的替代方案，用于推导响应关系。
• 通过数值模拟表明，所提出的涨落-响应关系在估计响应率方面，比直接小扰动方法更准确且更稳定。

提出的方法

• 基于随机轨迹的路径权重比，推导出使用安德森热浴驱动的惯性系统中热响应率的涨落-响应关系。
• 利用安德森热浴中速度重置的随机性，应用标准响应函数技术，避免了拉格朗日方法中常见的奇异性。
• 将响应率表示为可观测量和熵产生率的泛函，关键结果由公式（12）–（14）给出，涉及路径权重比的对称与反对称分量。
• 将该理论框架应用于计算FPU链长度对边界温度变化的响应率，使用平衡态与非平衡稳态两种情况。
• 对具有两个不同温度热浴的FPU链进行数值模拟，通过所推导公式计算响应率，并与直接扰动结果进行比较。
• 采用1/t的线性外推法，估计t → ∞时的静态响应率χs_L，以改善收敛性并减少统计噪声。

实验结果

研究问题

• RQ1 能否为具有边界热浴的惯性非平衡系统（如FPU链）推导出热响应率的涨落-响应关系？
• RQ2 所推导的涨落-响应公式的性能与直接小扰动方法相比，在估计热响应率方面如何？
• RQ3 系统长度对边界温度变化的响应率在平衡态与非平衡稳态下是否不同？这揭示了何种非平衡效应？
• RQ4 即使在非平衡状态下对称与反对称路径权重比分量（C+ 和 C−）发散，能否仍用其准确重构响应率？
• RQ5 在非平衡稳态下，静态响应率χs_L是否依赖于平均温度与热流方向？

主要发现

• 所推导的涨落-响应关系（公式12–14）为惯性非平衡系统提供了一种稳健且精确的热响应率计算方法，无需依赖小扰动。
• 数值结果表明，通过涨落-响应公式估计的响应率比直接扰动模拟获得的结果更平滑、噪声更小。
• 在平衡态下，响应的对称与反对称分量（C+ 与 C−）相等，符合预期，且响应率与标准库伯公式χ = 2C−一致。
• 在非平衡稳态下，C+ 与 C− 显著发散，表明存在强烈的非平衡效应，但完整公式仍能给出正确的响应率估计。
• 在平衡态下，静态响应率χs_L随温度升高而增加：T = 0.1时χs_L = 0.524(3)，T = 0.3时χs_L = 0.601(2)，T = 0.5时χs_L = 0.625(2)。
• 在非平衡态下，χs_L依赖于热流方向：当T1 = 0.1且TN = 0.5时，χs_L = 0.625(2)；当T1 = 0.5且TN = 0.1时，χs_L = 0.553(3)，表明存在真实非平衡效应，无法由平衡平均捕捉。