Thermodynamic Invariants of Coupled Channels: A Many-Channel Tolman-Ehrenfest Effect

本文通过将托尔曼-埃伦费斯特效应推广至耦合热力学通道，推导出一个涉及鲁平纳联络和乐的独特不变量，该不变量从几何上解决了颗粒物质物理学中长期存在的难题，并预测了剪切带的一个可检验关系。

要点

想象一下，你试图理解人群如何移动，或者一堆沙子在压力下如何发生位移。在旧有的思维方式（经典热力学）中，科学家将系统的不同部分视为彼此独立的“房间”。如果一个房间的温度发生变化，下一个房间发生了什么其实无关紧要；一切最终都会趋于单一、均匀的温度。

本文论证道，对于像致密沙子、湿土或活跃人群这样的复杂材料而言，“独立房间”这一概念是错误的。相反，万物交织成一张纠缠的网。如果你推挤沙子（应力），它会改变沙子的堆积方式（体积），而这两者相互影响之强烈，以至于你再也无法将它们分开描述。

以下是他们发现的拆解，使用了简单的类比：

1. “扭曲”的温度

在普通房间里，热量会流动，直到温度处处相同。但在这些复杂的耦合系统中，“温度”（对于沙子而言，它是衡量颗粒是处于抖动状态还是紧密堆积状态的指标）并不会保持均匀。

作者发现了一条隐藏的规则。这就像你在攀登一座山峰。在平坦的世界里，你只需直线行走。但在有强风（即“耦合”）的山上，为了保持在同一海拔高度，你必须沿着曲线行走。

他们发现了一个新的“不变量”（一条永不改变的规则）。该规则指出，如果你将局部的“温度”乘以一个特殊的“修正因子”（他们称之为 $\zeta$），无论你在系统中的哪个位置，结果总是一个相同的数值。

• 类比：想象货币兑换。如果你在一个国家持有美元，在另一个国家持有欧元，汇率会根据地点而变化。你不能简单地说“ everywhere 1 美元 = 1 欧元”。但是，如果你将你的美元乘以当地的汇率，你总是能得到相同的“真实价值”。在这篇论文中，“汇率”就是修正因子 $\zeta$，而“真实价值”则是系统的真实平衡态。

2. “隐藏的扭曲”（Holonomy）

为什么存在这个修正因子？论文使用了来自几何学的概念，称为“霍洛诺米”（holonomy）。

• 类比：想象你在平坦田野上的圆形跑道上行走。当你回到起点时，你面对的方向与出发时相同。现在，想象你在一个球体（如地球）上的跑道上行走。如果你从北极走到赤道，横跨赤道，再走回北极，当你回到起点时，你面对的方向与出发时不同。你被世界的形状“扭曲”了。

在这篇论文中，“世界的形状”就是材料的熵表面。由于不同的通道（体积和应力）是耦合的，在系统中绕一圈行走会“扭曲”温度。这种扭曲由 $\zeta$ 来衡量。如果通道没有耦合，就不会有扭曲，温度也会是均匀的（即旧有的、简单的观点）。

3. 解决困扰 60 年的沙子谜题

本文将这一理论应用于颗粒材料（如沙子）。60 年来，科学家们一直知道一条称为**罗氏定律（Rowe's Law）*的规则，它描述了沙子在剪切时如何膨胀（剪胀）与施加在其上的应力之间的关系。然而，一直存在一个令人烦恼的问题：该定律中的一个特定数值（称为 $K_\mu$）会随着沙子堆积密度的不同而改变。科学家们无法解释为什么*它会改变；他们只能每次重新测量它。

作者表明，这个变化的数值并非谜团；它只是修正因子 $\zeta$ 在起作用。

• 结果：当沙子松散时，通道未耦合，扭曲为零，旧规则完美适用。但当沙子变得非常紧密（接近“阻塞”状态）时，耦合变得巨大。修正因子 $\zeta$ 急剧增大，这正好解释了为什么数值 $K_\mu$ 似乎发生了变化。它并没有改变；我们只是忘记了将其乘以“汇率” $\zeta$。

4. 这对实验意味着什么

这篇论文不仅仅是在做数学推导；它提供了两种在现实世界中进行测试的具体方法：

1. 均匀性测试：如果你观察剪切带（沙子正在滑动的区域），“温度”（可动性，compactivity）和“应力温度”（角动性，angoricity）看起来会杂乱无章且不均匀。但是，如果你将它们乘以各自的修正因子，结果在整个剪切带上应该是完美平滑且均匀的。
2. 尺度测试：沙子开始表现异常（即修正因子急剧上升）的点，应该发生在非常特定的尺度上，这与沙子内部结构重新排列的速度有关。

总结

该论文声称，当复杂系统相互作用时，你不能将其各部分视为独立的。系统中存在一种几何上的“扭曲”，迫使你调整测量值。通过应用这种调整（即 $\zeta$ 因子），他们解决了一个困扰 60 年的谜题，即为什么沙子在阻塞状态下表现不同，表明这种“怪异”行为实际上是系统形状带来的可预测的几何后果。

技术摘要

技术摘要：耦合通道的热力学不变量

问题陈述
经典热力学依赖于熵表面有效对角化的假设，从而允许热力学通道独立达到平衡。然而，对于致密颗粒介质、反应性地球物理流体和活性物质等系统而言，热力学通道之间的耦合是主导物理机制。包括 Blumenfeld–Edwards 颗粒统计力学、规范理论热力学和交叉扩散约化在内的现有框架，虽然承认耦合修正的必要性，但未能推导出耦合通道在平衡态下所共享的特定不变量。因此，均匀强度变量（例如恒定温度）的标准条件在这些耦合机制下失效。

方法论
作者将最初为相对论引力场 formulated 的 Tolman–Ehrenfest (TE) 效应，推广至耦合热力学系统中的熵流形。

1. 几何表述：研究利用 Ruppeiner 度规 $g_{ij} = -\partial^2 S / \partial q_i \partial q_j$ 来量化对角近似的失效。非对角分量（$i \neq j$ 时的 $g_{ij}$）代表广延变量 $q_i$ 与其共轭强度 $\beta_i$ 之间的耦合。
2. 不变量的推导：通过分析熵流形上的水平集流（$dS = \sum \beta_i dq_i = 0$），作者推导出了支配强度演化的微分方程。他们寻求一个沿该流所有轨迹保持恒定的第一积分 $F(\beta_1, \dots, \beta_n)$。
3. 和乐与联络：该解涉及定义一个源自 Ruppeiner 联络的耦合系数 $\zeta_i$。作者证明，该系数的对数微分 $d \log \zeta_i$ 对应于一个由度规分量和强度构造的一形式 $\omega_i$。该一形式沿闭合回路的积分代表了 Ruppeiner 联络的和乐，仅当通道解耦时该和乐才为零。

主要贡献与结果
主要结果是推导出了多通道 Tolman–Ehrenfest (MCTE) 关系：
$$\zeta_i T_i = C$$
其中 $T_i = \beta_i^{-1}$ 是强度变量，$C$ 是跨越所有通道和点的单一标量不变量，$\zeta_i$ 是耦合系数（Ruppeiner 联络的和乐）。

• 颗粒实现：该框架被应用于 Edwards 统计力学中的双通道（体积 $V$ 和应力 $\sigma$）颗粒系综。在此，温度被替换为压缩性 $\chi$，应力共轭于角动量 $A$。
• 耦合的几何起源：非对角 Ruppeiner 曲率 $g_{V\sigma}$ 被确定为耦合的几何驱动力。论文表明，$g_{V\sigma}$ 源于从三通道系统中绝热消除一个快速中间通道（机械织构自由度），从而将耦合印记在约化后的双通道熵表面上。
• 定量基准：使用玩具熵表面模型，作者证明，虽然裸压缩性 $\chi$ 沿熵水平集显著变化（变化因子为 3），但乘积 $\zeta_V \chi$ 在机器精度范围内保持恒定（$\sim 10^{-13}$）。在接近阻塞相变时，修正因子 $\zeta_V$ 偏离单位值 $O(1)$，表明单通道近似会引入此量级的系统误差。

关于颗粒塑性的三个具体预测
该论文在不引入自由参数的情况下，推导出了颗粒系统的三个可检验推论：

1. 剪胀比：剪胀比被证明是加权了强度变量比率的非对角 Ruppeiner 曲率的直接函数（$D \approx \frac{g_{V\sigma}}{g_{VV}} \frac{\chi}{A}$）。
2. Rowe 能量限制：Rowe 定理中摩擦耗散的非负性被证明等价于 $2 \times 2$ Ruppeiner 度规的半正定性。因此，Rowe 界限是热力学稳定性（$\delta^2 S \leq 0$）的结果，而非独立的唯象假设。
3. 状态依赖性 $K_\mu$ 的解决：解决了关于阻塞附近应力 - 剪胀系数 $K_\mu$ 状态依赖性的 60 年谜题。正确的关系是 $\sigma'_1 / \sigma'_3 = \zeta_V(\sigma) K_\mu R$。$K_\mu$ 偏离常数的现象由非对角熵曲率几何决定，其中 $\zeta_V$ 在接近阻塞时产生 $O(1)$ 的修正。

意义与主张
该论文声称提供了当热力学通道耦合时替代裸强度变量的唯一不变量。它将经典的均匀温度条件和相对论 Tolman 温度推广至统一的几何框架。

作者断言，MCTE 框架通过表明经典结果（Rowe 剪胀、能量限制以及 $K_\mu$ 的状态依赖性）并非独立的经验定律，而是非对角 Ruppeiner 曲率的几何后果，从而解决了颗粒塑性领域的长期问题。

实验与计算测试
该论文提出了两个直接测试以验证 MCTE 图景：

1. C-均匀性测试：在发展的剪切带中，单个强度（$\chi, A$）应在空间上不均匀，而乘积 $\zeta_V \chi$ 应保持空间均匀。这将 MCTE 模型与经典的多温度替代方案区分开来。
2. 长度尺度测试：修正后的 $K_\mu$ 出现 $O(1)$ 偏离的起始点应发生在快速织构通道的扩散长度尺度 $\ell \sim D_w^{1/2}$ 处。该尺度被预测与对偶反应 - 扩散系统中准孤子波的波数重合，这一巧合可通过离散单元法 (DEM) 模拟进行检验。

该工作得出结论：$N=2$ 框架描述了平衡背景态，而涉及奇宇称系统和非关联流动的 $N=3$ 动力学将在未来的工作中予以探讨。