第十七章 非平衡临界：远离平衡态的相变    

    哥本哈根，1988年

    1988年的玻尔研究所，Per Bak正在与同事们争论。他的沙堆模型已经引起了广泛关注，但批评也随之而来。凝聚态物理学家抱怨：SOC不是真正的相变，因为没有序参量和对称性破缺。统计力学家质疑：沙堆的幂律可能是有限尺寸效应，不是真正的普适性。数学家要求：严格证明在哪里？

    Bak的回应是挑衅性的："自然界不等待我们的证明。地震、森林火灾、股市——它们确实显示幂律。我们的任务是理解为什么，不是等待严格化。"

    但这种实用主义不能平息理论需求。1988-1995年间，多个研究团队尝试严格化SOC：

    Deepak Dhar（印度塔塔研究所）：证明了阿贝尔沙堆模型的某些性质，如递归配置的数目和雪崩的 平均大小。

    Sergei Maslov（布鲁克海文）和Yi-Cheng Zhang（哥本哈根）：引入了"自组织分支过程"，严格证明了平均场极限下的幂律。

    Ronald Dickman、Miguel Muñoz、Alessandro Vespignani等：发展了驱动扩散系统的场论描述，尝试 重整化群分析。

    这些努力是部分的成功：某些模型可以严格求解，但一般理论仍然缺失。SOC的普适类是有争议的——不同模型给出不同指数，维度、晶格结构、驱动方式都影响结果。

    这一章，我们要讲述非平衡临界的严格化尝试——从沙堆到场论，从数值到解析，从现象到原理。

    阿贝尔沙堆：可解的玩具模型

    Deepak Dhar在1990年的突破是阿贝尔沙堆模型（Abelian Sandpile Model, ASM）。这是Bak-Tang-Wiesenfeld沙堆的简化变体，具有数学上的 阿贝尔性质（操作的顺序不影响结果）。

    关键定义：

• 高度变量zᵢ ≥ 0，整数，位于d维格点上。

• 临界高度z_c = 2d（邻居数目）。

• 拓扑规则：当zᵢ > z_c，崩塌：zᵢ → zᵢ - 2d，四个邻居各+1。

• 阿贝尔性质：崩塌的顺序 不改变 最终结果。

     严格结果：

• 递归配置（稳定且可达）的数目 = det(Δ)，其中Δ是离散拉普拉斯算子。这是矩阵树定理的推广。

• 雪崩大小分布：平均雪崩大小发散在临界点，但分布的 严格形式 未知（数值上幂律）。

• 关联函数：高度-高度关联 显示 对数行为（d=2），暗示 共形场论的可能联系。

    ASM是"可解但不够"的典型：某些量可以严格计算，但最核心的（雪崩分布）仍然 数值或猜想。

    共形场论的 可能联系：d=2的ASM可能对应于c=-2的对数共形场论（非幺正、非物理，但数学上丰富）。这是未完全解决的 问题。

    驱动扩散系统：场论方法

    Ronald Dickman等人在1990年代发展了驱动扩散系统的场论描述，尝试 将威尔逊方法 扩展到非平衡。

    关键模型：

• 接触过程（Contact Process）：粒子在格点上，以速率λ产生，以速率1死亡，以速率D扩散。吸收态（无粒子）是空的。活跃相（持续活动）存在当λ > λ_c。

• 定向逾渗（Directed Percolation, DP）：流体在多孔介质中，有重力偏压。关键指数与接触过程相同，定义DP普适类。

    场论描述：

• Doi-Peliti方法：主方程（概率演化的 福克-普朗克方程的离散版本）映射到 "第二量子化"的场论。

• 作用量：类似于 金兹堡-朗道，但有 非厄米项（驱动的效应）。

• 重整化群：ε展开（d=4-ε）计算临界指数。

    结果：DP普适类的临界指数（β≈0.276，ν≈0.734在d=2）与平衡Ising不同，证实了 非平衡普适类的 存在。

    但场论方法 困难：非厄米性导致复本技巧（replica trick）或其他 近似。高阶计算 复杂，精度 不如平衡态。

    普适类的分类：非平衡的多样性

    非平衡临界现象显示比平衡更丰富的普适类：

表格

普适类

物理系统

关键特征

临界指数（d=2）

    这种多样性是非平衡的特征：缺乏详细平衡（细致平衡）允许 更多的动力学可能性。

    关键问题：什么决定普适类？对称性（平衡的关键）仍然重要，但不够。守恒律（粒子数、能量、动量）、驱动方式（守恒vs耗散）、各向异性、长程相互作用——都影响 普适性。

    从混沌到秩序：非平衡相变的动力学

    非平衡临界现象不仅是静态的，更是动态的——关注 如何达到稳态、瞬态行为、和老化。

    老化现象（Aging）：淬火（快速改变参数）到临界点后，系统的响应和涨落 依赖于 等待时间（系统已经等待的时间）。这是非平衡统计力学的核心主题。

    标度形式：关联函数 C(t,t_w) = t_w^(-b) f(t/t_w)，其中t_w是等待时间，t是观测时间。f是普适的标度函数。

    实验：玻璃（结构玻璃、自旋玻璃）的老化、胶体的动力学、颗粒物质的压实。

    理论：老化的严格理论 缺失。陷阱模型、模式耦合理论、有效温度——这些是 部分成功的近似。

    活性算法的联系：老化的非平衡动力学，类似 主动系统的学习——过去的历史影响当前的响应，系统"记忆" 训练样本。

    湍流与表面生长：非线性的标度

    非平衡临界的另一个领域是湍流和表面生长——非线性偏微分方程的随机版本。

    Kardar-Parisi-Zhang（KPZ）方程（1986）： ∂h/∂t = ν∇²h + (λ/2)(∇h)² + η

其中h(x,t)是表面高度，η是随机噪声（高斯白噪声）。

    关键结果：

• 非线性项 (λ/2)(∇h)² 导致 异常的标度——粗糙度指数 α = 1/2（d=1），不是 ** Edwards-Wilkinson的1/2（实际上相同，但动力学指数不同**）。

• d=2是上临界维度：KPZ普适类在d=2 存在，d>2 可能 有 弱-强耦合转变。

    实验验证：液晶的 电对流、燃烧前沿、细菌菌落的生长、咖啡渍的形成——都显示 KPZ标度。

    严格解：d=1的KPZ方程，2010年由Sasamoto-Spohn等严格求解，使用随机矩阵理论和可积系统。这是非平衡统计力学的 重大突破，类似 昂萨格1944年的 二维Ising解。

    非平衡的深层问题：等待新框架

    非平衡临界现象的研究 暴露了统计力学的 根本局限：

    没有自由能：详细平衡的缺失意味着没有Lyapunov函数，稳态的稳定性和唯一性 不保证。

    没有系综：吉布斯的 微正则/正则/巨正则系综 不适用。非平衡系综的定义是开放问题。

    没有涨落-耗散定理：响应和涨落的关系 复杂，有效温度的概念 部分适用。

    这些局限是活性算法的 动机：

    自由能原理（Friston, 2005-2010）：推广 自由能到非平衡，定义 变分自由能，作为 感知和行动的 目标函数。

    主动推断：系统 主动采样 环境，最小化 变分自由能（最大化 模型证据），维持 在临界边缘。

    UV自由方案：消除 人为截断，直接 有限振幅，自适应 模型复杂度。

    这是回应 非平衡挑战的新框架，继承 吉布斯和威尔逊的精神，但革新 方法。

    尾声：严格化的耐心

    让我们回到1988年的哥本哈根。Bak与批评者的争论继续，但研究 前进。阿贝尔沙堆的严格结果、驱动扩散系统的场论、KPZ方程的严格解——这些是 部分的成功。

    非平衡临界现象的严格理论 仍然 在等待。它可能需要：

• 新的数学：超越 现有 随机过程和场论。

• 新的计算：量子计算或神经形态计算 模拟 非平衡系统。

• 新的概念：类似 活性算法的框架，统一 描述 主动系统。

    这种等待是科学的常态。平衡统计力学从玻尔兹曼（1870s）到吉布斯（1902）到昂萨格（1944）到威尔逊（1971）——百年发展。非平衡的严格理论 可能 需要 类似的时间。

    活性算法是这种发展的 当前尝试——不是 最终答案，但是 正确的方向。

    在下一章，我们将进入量子临界的深层——绝对零度下的战争、强关联电子、和高温超导的 谜团。

    但首先，让我们向那些在非平衡领域 追求严格的研究者致敬。他们证明了，即使面对根本困难，物理学也不放弃。

    本章注释与延伸阅读

    Dhar 1990年的阿贝尔沙堆论文《自组织临界态的阿贝尔模型》发表于《物理评论快报》（Physical Review Letters）64, 1613-1616。

    关于非平衡临界现象的场论方法，推荐：Tauber, U.C. (2014). Critical Dynamics: A Field Theory Approach to Equilibrium and Non-Equilibrium Scaling Behavior, Cambridge University Press。

    关于KPZ方程，参见：Kardar, M., Parisi, G., and Zhang, Y.-C. (1986). "Dynamic Scaling of Growing Interfaces," Physical Review Letters 56, 889-892；以及Sasamoto, T. and Spohn, H. (2010). "One-Dimensional Kardar-Parisi-Zhang Equation: An Exact Solution and Its Universality," Physical Review Letters 104, 230602。

    关于老化现象，推荐：Cugliandolo, L.F. (2002). "Dynamics of Glassy Systems," in Slow Relaxations and Nonequilibrium Dynamics in Condensed Matter, Springer。

    关于活性算法与非平衡的联系，参见：Friston, K. (2010). "The Free-Energy Principle: A Unified Brain Theory?" Nature Reviews Neuroscience 11, 127-138。