Topological Foundations of Multi-Field Instabilities in Continua: Part 1: Foundations Part 2:Analytical Formulation for 1-D Spin Chains Part 3: Numerical Upscaling

技术摘要：连续介质中多场不稳定性之拓扑基础：第二部分 —— 一维自旋链的解析公式化

问题陈述
尽管单个颗粒间的相互作用具有时间可逆性，但颗粒材料仍表现出不可逆的宏观变形（剪胀与失效）。传统力学将这种不可逆性主要归因于颗粒接触处的摩擦耗散。本文认为这种归因是不完整的，提出几何耦合机制在颗粒接触尺度上发挥着基本作用，且该作用独立于摩擦。本文旨在解决的核心问题是：识别一维颗粒链中构型不可逆性的参数无关的、拓扑学起源，特别是隔离出驱动剪胀驱动矩阵不稳定性（dilatancy-driven matrix instability）的机制，而不依赖于经验能量景观或摩擦系数。

方法论
本文基于增广 Onsager 框架（$A = D + L$）开发了一个解析公式，其中 $D$ 是对称耗散块，$L$ 是反对称保守块。研究方法通过以下步骤进行：

1. 拓扑分解： 将颗粒接触状态分解为最小的体积—力学—构型（VMC）三元组（$N=3$）。
   • $V$：体积通道（重叠/法向力）。
   • $M$：力学通道（切向/剪切传递）。
   • $C$：构型通道（织构取向/接触法线旋转）。
2. 图论过滤： 将一维颗粒链建模为路径图 $P_N$。本文利用 Hodge–Helmholtz 分解证明，在一维情况下，第一贝蒂数 $\beta_1 = 0$。这一拓扑约束充当了“上循环过滤器”（cocycle filter），消除了所有旋度驱动的环流模式（solenoidal cycles），并将系统的演化完全限制在上循环（梯度）子空间内。
3. 李代数映射： 将 $N=3$ 三元组的复可逆耦合块 $L$ 识别为李代数 $\mathfrak{so}(3)$ 的一个元素。构型通道 $C$ 被重新解释为经典自旋自由度。算符 $L$ 作为旋转生成元，使状态向量绕耦合轴 $\omega$ 发生进动。
4. 解析推导： 本文针对 $N=3$（奇数）和 $N=4$（偶数）通道数的线性演化方程 $\dot{q} = Lq + f_0$（在无摩擦极限下 $D=0$），推导出了闭式解析解。通过与直接数值积分进行对比，验证了代数预测的准确性。

核心贡献

• VMC 三元组作为最小门户： 本文确立了三通道 VMC 接触是能够产生“门户层”（Gateway Layer）的最小拓扑单元。不同于形成“稳定层”（Stable Layers）并具有闭合辛轨道（symplectic orbits）的偶数维系统（$N=2, 4$），奇数维系统（$N=3$）由于宇称定理（Parity Theorem, $\det(L)=0$）的存在，拥有一个结构零特征值。
• 不可逆性的自旋重释： 通过将接触动力学映射到 $\mathfrak{so}(3)$，本文证明了由宇称决定的零模（null mode）恰好对应于耦合生成元的旋转轴。由于旋转算符在其自身轴向上不产生力矩，任何投影在该轴上的边界强迫都无法通过保守块进行耗散或重定向。这导致了沿零模方向的久期漂移（secular drift，即无界增长），为即便在无摩擦情况下依然存在的构型不可逆性提供了一种几何学解释。
• 上循环过滤器： 本文严格定义了一维笛卡尔限制并非一种简化，而是一种精确的拓扑过滤器。本文证明了一维链的无环性质（$\beta_1=0$）在数学上禁止了 Schnakenberg 环流和旋度驱动的剪切局部化。因此，在一维空间中，唯一的不可逆路径是由上循环携带的零模漂移，这在宏观上表现为剪胀。
• 无摩擦不可逆性： 分析表明，不稳定性的产生是接触拓扑（即奇数通道数 $N=3$）的结果，而非接触摩擦。即使在存在耗散（$D > 0$）的情况下，门户数 $G = \|L\|_2^2 / (\sigma_{min}(D)\sigma_{max}(D))$ 决定了弛豫前的瞬态过冲幅度，但漂移机制的存在是拓扑性的，且独立于 $G$。

结果

• 宇称驱动动力学：
  • 奇数 $N$（门户层）： 系统沿零特征向量 $v_0 \propto (L_{MC}, 0, L_{VM})^T$ 表现出久期漂移。在持续的轴向边界强迫下，零模振幅随时间线性增长（$v_0^T q(t) \propto t$），代表了从初始状态的无界逃逸。
  • 偶数 $N$（稳定层）： 系统表现出谐波约束。动力学被限制在一个紧凑的不变环面上，表现为围绕偏移平衡点的有界准周期振荡。
• 闭式解： 本文提供了 $N=3$ 和 $N=4$ 状态演化的显式解析表达式。这些解证实了奇数 $N$ 系统会发生久期漂移，而偶数 $N$ 系统保持有界，从而验证了两者之间的理论差异，且无需数值近似。
• 激活机制： 本文阐明了在无摩擦极限（$D=0$）下，漂移是精确且无条件的。在存在耗散（$D > 0$）时，门户数 $G$ 决定了弛豫前的瞬态过冲程度，但漂移机制的存在是拓扑性的，独立于 $G$。

意义与主张

本文声称为一维颗粒不可逆性提供了一个无参数的、拓扑的基础。其主要意义在于将剪胀和构型漂移的起源与摩擦耗散解耦。通过将 VMC 三元组识别为不稳定的不可约生成元，作者认为颗粒不可逆性是接触空间连通性的一种涌现属性，类似于经典陀螺仪中被保留的自旋轴。

作者断言，该框架提供了一种严谨的替代方案，用来自代数必然性取代经验性的本构建模，而这种代数必然性源于热力学通道数的宇称。这项工作构成了连接第一部分建立的通用拓扑原则与第三部分详细描述的数值验证及 DEM 对应关系的分析桥梁。本文明确限定了其研究范围仅限于线性算符演化，指出非线性饱和以及向非线性波动方程的映射属于未来的研究工作。