[论文解读] Are there really phase transitions in 1-d heat conduction models?

Lei Yang, Peter Grassberger|arXiv (Cornell University)|Jun 6, 2003

Advanced Thermodynamics and Statistical Mechanics被引用 3

一句话总结

本文挑战了关于一维热传导模型中有限温度相变的既有主张，特别是Frenkel-Kontorova（FK）模型与耦合转子模型，此前研究曾声称在特定温度下热导率趋于无穷。通过采用高精度、大规模的Langevin动力学模拟，积分时间长、系统尺寸最大达N=8192，作者表明，看似导率发散的现象实为在T→0和T→∞时的谐振极限所导致的缓慢交叉过渡，并非真正的相变。

ABSTRACT

Recently, it has been claimed (O. V. Gendelman and A. V. Savin, Phys. Rev. Lett. {\bf 84}, 2381 (2000); A.V.Savin and O.V.Gendelman, arXiv: cond-mat/0204631 (2002)) that two nonlinear classical 1-d lattice models show transitions, at finite temperatures, where the heat conduction changes from being finite to being infinite. These are the well known Frenkel-Kontorova (FK) model and a model for coupled rotators. For the FK model we give strong theoretical arguments why such a phase transition is not to be expected. For both models we show numerically that the effects observed by Gendelman {\it et al.} are not true phase transitions but are rather the expected cross-overs associated to the conductivity divergence as $T o 0$ and (for the FK model) $T o\infty$.

研究动机与目标

• 检验Frenkel-Kontorova（FK）模型在有限温度下是否表现出相变，使热导率从有限值变为无穷大。
• 研究耦合转子模型是否也表现出类似相变，如Gendelman和Savin先前所报告。
• 解决早期模拟与一维非线性晶格热传导理论预期之间的矛盾。
• 确定观察到的导率发散是否为真实相变，还是有限尺寸效应与非平衡动力学导致的人为现象。
• 阐明温度依赖的声子散射与局域化在决定一维系统热导率中的作用。

提出的方法

• 对长达8192个粒子的一维链执行大规模Langevin动力学模拟，其中N₀ = 40个粒子为热浴边界层。
• 使用时间步长为0.05的跃迁积分器，积分时间长达5×10⁸个单位，以确保稳态收敛。
• 施加白高斯噪声热浴，温度分别为T_high与T_low，以产生恒定的温度梯度（ΔT ≈ 10–20%）。
• 测量热流J，并根据系统尺寸N与温度T计算有限尺寸导率κ = J·N/ΔT。
• 分析κ的尺寸依赖性，以区分真正的发散（指示相变）与缓慢的交叉行为。
• 将结果与Gendelman和Savin的早期主张进行对比，重点关注低温与高温下κ的行为。

实验结果

研究问题

• RQ1 Frenkel-Kontorova模型与耦合转子模型是否表现出真正的有限温度相变，使热导率从有限值变为无穷大？
• RQ2 在低温与高温下观察到的导率发散是否为真实相变，还是由于谐振极限导致的交叉过渡？
• RQ3 在固定温度下，导率κ(N)的尺寸依赖性如何区分有限导率与发散行为？
• RQ4 在转子模型中，局域激发与声子散射在低温下如何阻碍热输运？
• RQ5 为何早期模拟在某些温度范围内报告无限导率？这是否由于系统尺寸或积分时间不足所致？

主要发现

• 对于ε=1.0的Frenkel-Kontorova模型，所有研究温度（T=0.125至T=128）下导率κ始终保持有限，未见发散或相变迹象。
• 在T=0.125与T=128时，κ随N缓慢增加，但当N>1000时趋于饱和，表明为交叉过渡而非真正发散。
• 对于ε=3.0与ε=10.0，所有有限温度下κ均保持有限，且在低温下亦无发散迹象，与早期主张相矛盾。
• 在转子模型中，T=0.2、0.3、0.45与0.6时，κ随N增加，但大N下曲线最终趋于平坦，表明所有有限温度下导率均为有限值。
• 早期研究中观察到的导率发散归因于系统尺寸与积分时间不足，而非真实相变。
• 在T→0与T→∞时观察到的行为与系统趋于有效谐振极限一致，导致缓慢交叉过渡行为，而非临界相变。