[论文解读] Solution of the Metropolis dynamics on a complete graph with application to the Markov chain Mpemba effect

Israel Klich, Marija Vucelja|arXiv (Cornell University)|Dec 31, 2018

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

一句话总结

本文通过推导完全图上梅特罗波利斯动力学的所有本征值和本征态，解析求解了该系统的动力学，揭示了由于单调弛豫，系统不表现出弱或强梅特罗波利斯效应。然而，当跃迁受阻时（例如在稀疏图上），梅特罗波利斯效应——即更热的系统冷却得更快——可以出现，如一个缺少一条边的三状态模型所示。

ABSTRACT

We find analytically the complete set of eigenvalues and eigenvectors associated with Metropolis dynamics on a complete graph. As an application, we use this information to study a counter-intuitive relaxation phenomenon, called the Mpemba effect. This effect describes situations when upon performing a thermal quench, a system prepared in equilibrium at high temperatures relaxes faster to the bath temperature than a system prepared at a temperature closer to that of the bath. We show that Metropolis dynamics on a complete graph does not support weak nor strong Mpemba effect, however, when the graph is not complete, the effect is possible.

研究动机与目标

• 解析确定完全图上梅特罗波利斯动力学的完整本征系统，尽管其依赖于能量层级具有非解析性。
• 研究梅特罗波利斯效应——即从更高初始温度更快弛豫——是否可能在该系统中出现。
• 探索在梅特罗波利斯动力学控制的马尔可夫跳跃过程中，梅特罗波利斯效应出现的条件。
• 对比完全图与不完全图的行为，以识别梅特罗波利斯效应产生的结构要求。

提出的方法

• 使用置换矩阵 $ S_r $ 将稳态概率 $ \pi_i $ 按降序重排，以简化谱分析。
• 在重排基下推导梅特罗波利斯速率矩阵 $ R $ 的本征值 $ \lambda_k $ 和本征态 $ |v_k\rangle $，其中 $ \lambda_1 = 0 $，且对 $ k > 1 $，有 $ \lambda_k = -\left(\pi_{r_{k-1}}^{-1}Z_{k-1} + k - 2\right) $。
• 用累积概率 $ Z_k = \sum_{l=k}^n \pi_{r_l} $ 表达衰减最慢的本征态 $ |v_2\rangle $，从而实现对弛豫动力学的分析。
• 应用卢与拉兹（2017）的公式计算控制弛豫速率的系数 $ a_2(\beta) $，使用本征态 $ |v_2\rangle $ 和浴的能级分布。
• 分析 $ a_2(\beta) $ 的单调性，以排除弱或强梅特罗波利斯效应所需的非单调行为。
• 构建一个三状态模型，其中一条跃迁被阻断，以证明在稀疏图中梅特罗波利斯效应的出现，显示 $ a_2(\beta) $ 的非单调性。

实验结果

研究问题

• RQ1 完全图上的梅特罗波利斯动力学是否支持梅特罗波利斯效应，即非单调弛豫行为？
• RQ2 跃迁图的哪些结构特征（如完全性）是梅特罗波利斯效应发生的必要条件？
• RQ3 完全图上梅特罗波利斯速率矩阵的谱分解如何影响弛豫动力学？
• RQ4 当图中移除直接跃迁时，梅特罗波利斯效应是否可能在马尔可夫跳跃过程中出现？
• RQ5 弛豫系数 $ a_2(\beta) $ 的单调性在决定梅特罗波利斯效应是否存在方面起什么作用？

主要发现

• 完全图上梅特罗波利斯动力学的完整本征系统可被解析求解，本征值为 $ \lambda_k = -\left(\pi_{r_{k-1}}^{-1}Z_{k-1} + k - 2\right) $，本征态 $ |v_k\rangle $ 以累积概率 $ Z_k $ 表示。
• 弛豫系数 $ a_2(\beta) $ 在逆温度 $ \beta $ 上是单调的，这意味着在完全图情况下，弱或强梅特罗波利斯效应均不会发生。
• 完全图中梅特罗波利斯效应的缺失是由于 $ e^{\beta E_1}Z_1(\beta) $ 的严格单调性，这确保了 $ a_2(\beta) $ 仅在 $ \beta = \beta_b $ 处为零，且无中间零点。
• 相比之下，当一条跃迁被阻断时（例如在三状态系统中），弛豫系数 $ a_2(\beta) $ 变为非单调，满足 $ a_2(\beta_b) = 0 $ 且 $ \lim_{\beta \to \infty} a_2(\beta) = 0 $，从而实现间接的梅特罗波利斯效应。
• 梅特罗波利斯效应在稀疏图中出现，是由于对称性的破坏和弛豫动力学的非单调性，尽管在完全连接的情况下该效应并不存在。
• 结果表明，梅特罗波利斯效应并非梅特罗波利斯动力学的普遍特征，而是关键依赖于底层跃迁图的拓扑结构，稀疏图可实现非平凡的弛豫现象。