← 返回 PaperDaily 大模型与智能体

相变理论70年后的大一统:边界条件彻底沦为“工具人”

Lee-Yang零点理论诞生近八十年后,Rutgers团队给出了一个异常干净的定理:只要配分函数的零点在复数平面上不贴着实轴,粒子数密度和涨落就乖乖收敛到热力学极限,且边界条件说了不算。把相变与零点分布之间的因果链又拧紧了一圈。

原论文信息如下:
论文标题:
Lee-Yang Zeros And Particle Fluctuations

发表日期: 2026年7月

发表单位: Rutgers大学数学系(美国)

原文链接: https://arxiv.org/pdf/2607.26975v1.pdf


在统计力学的世界里,有一个让无数物理学家又爱又恨的名字——Lee-Yang零点。1952年,李政道和杨振宁把配分函数当成复变函数来研究,发现相变的秘密全部藏在配分函数零点的分布里。然而七十多年过去,零点与可观测物理量之间的严格数学关系,仍然有许多空白。
这篇文章做的事情,一句话概括:证明了只要Lee-Yang零点在复平面上离实轴足够远,那么有限体积的粒子数密度和粒子数方差在热力学极限下收敛到同一组值,且与边界条件选取无关。
本篇文章是理论数学物理方向,没有实验和代码实现,核心是严谨的定理证明与数学推导。因此文章重点放在:用通俗的语言解释Lee-Yang零点、Gibbs测度等核心概念,拆解定理含义与证明思路,梳理可应用的具体物理系统,阐明边界条件的推广逻辑,最后以龙迷三问与龙哥点评收尾。文中无表格可展示,但会在合适位置插入公式的HTML文本形式进行说明。

一杯水的启示:相变背后的零点密码

想象一下这样的场景:你站在厨房里盯着锅中的水,眼看着它从平静的液面开始翻涌,冒出无数气泡,最终变成翻滚的蒸汽。从液态到气态的转变,是一个典型的相变过程。但为什么偏偏是100℃?水分子之间到底发生了什么,让宏观性质突然改变?这些问题,统计力学家们已经追问了一个多世纪。
在平衡态统计力学里,描述一个系统的核心对象是巨正则配分函数(grand canonical partition function),记作 Ξ(μ; Λ, b)。其中 μ 是化学势,Λ 是系统所占的体积区域,b 表示边界条件。配分函数的对数除以体积得到压力 p(μ; Λ, b) = |Λ|⁻¹ log Ξ(μ; Λ, b),而物理量的期望值都可以通过压力对化学势的导数来获得。
具体来说,粒子数密度是压力的一阶导数除以 β,即 ρ(μ; Λ, b) = β⁻¹ ∂p/∂μ;而粒子数密度涨落(也就是粒子数方差除以体积)则对应压力的二阶导数,σ_Λ/|Λ| = β⁻² ∂²p/∂μ²。这个二阶导数跟系统的等温压缩率直接相关,是探测相变的关键信号。在液气临界点附近,压缩率发散,密度涨落变得异常剧烈,甚至能让光发生强烈散射——这就是著名的临界乳光(critical opalescence)现象。早在1910年,爱因斯坦就对此给出了理论解释。从那时起,通过密度涨落来研究相变就成了统计力学最经典的手段之一。
然而这里有一个十分微妙的问题:真实系统总是有限的,而我们关心的相变行为发生在体积趋向无穷的极限下。有限体积系统的粒子数方差 σ_Λ 和无限体积下的二阶导数 β⁻² ∂²p/∂μ² 之间,到底是什么关系?在不同的边界条件下,这个极限是否一致?举个例子,在某个温度下,同一个系统可以同时存在多个不同的Gibbs测度(Gibbs measures),即满足Dobrushin–Lanford–Ruelle(DLR)方程的无限体积概率分布。它们分别描述不同的宏观状态,比如密度高低不同的两个共存相。这些不同的 Gibbs 测度在局部性质上可以完全不同,但它们的压力是相同的。那么,它们各自的粒子数涨落在热力学极限下会收敛到同一个值吗?
这个问题还有一个更加奇特的背景:超均匀性(hyperuniformity)。超均匀系统指的是粒子数涨落增长慢于体积,例如只按表面积增长,导致压缩率在热力学极限下趋于零。零温下的周期晶体和许多准晶体就属于这种类型。但在正温条件下,经典平衡态系统真的能保持超均匀吗?早在1967年,Ginibre 就给出了一个著名下界:对于具有正对势或者硬核势的系统,有限体积的压缩率不会太小,涨落总是广延的,从而排除了正温超均匀性。Dereudre 和 Flimmel 后来把这一结果推广到了更广的 Gibbs 点过程类别,证明了对于超稳定、下正则、可积对势的系统,超均匀性不会出现。
OIP-C.jpeg
这个问题的难点在于:不同的边界条件可能挑选不同的无限体积 Gibbs 测度,不同 Gibbs 测度可以有不同的压缩率。更微妙的是,粒子数方差是一个"非局域"的量,它对边界条件的选择极其敏感。而本论文抓住了统计力学中最经典的工具——Lee-Yang 零点的分布,一举打通了"零点远离实轴"和"涨落收敛到共同极限"之间的逻辑链条。

核心定理:零点绕着走,导数就听话

1952年,李政道和杨振宁提出了一个革命性的想法:把配分函数看成复变量 z = e^(βμ) 的函数,研究其在复平面上的零点分布。他们证明了一条深刻的结论:在这个复数平面上,相变的发生必然伴随着零点在热力学极限下"逼近"正实轴。如果零点在某个实轴上点的邻域内完全不出现(对足够大的体积一致地保持有界距离),那么这个点附近就不会有相变,压力是解析函数。这就是著名的Lee-Yang零点理论
本篇论文的主要结果,可以看作Lee-Yang理论向更高阶导数的一个延伸。为了看清这一点,首先需要注意:有限体积压力 p(μ; Λ, b) 对 μ 是凸函数(实轴上 z > 0),因为巨正则配分函数的系数都是非负的。凸函数的导数在极限下有一个众所周知的特性:如果一列凸函数逐点收敛到某个极限函数,那么它们在极限函数可导的点处,一阶导数也收敛。这意味着密度 β⁻¹∂p/∂μ 的极限行为相对"仁慈"——只要极限压力可导,有限体积的密度就自动收敛,无需任何零点条件。
但二阶及以上的导数就没有这种好运了。凸性只能保证单调性之类的性质,完全无法控制二阶导数的极限行为。这就是为什么需要零点假设——正是在这里,零点条件开始发挥作用。
fine.JPEG
论文的核心定理(Theorem 1.1)可以这样表述:
设 B 是一个可容许边界条件集合,所有边界条件共享同一个热力学压力 p(x) = lim_{L→∞} p(x; Λ_L, b)(在实轴 x > 0 上一致,且与 b 无关)。固定实数 z₀ > 0,假设存在一个包含 z₀ 的复平面邻域 B_δ(z₀),配分函数在其中没有零点,且这个性质对所有足够大的 L 一致成立。那么压力函数 p 可以全纯延拓到 z₀ 的一个邻域,并且对任意正整数 n,都有 lim_{L→∞} ∂ⁿ_z p(z₀; Λ_L, b) = ∂ⁿ_z p(z₀)。也就是说,热力学极限和求导可以在 z₀ 处交换顺序。
定理的推论直接表达了物理意义:令 μ₀ = β⁻¹ log z₀,那么有限体积的密度 ρ(μ₀; Λ, b) 趋向 β⁻¹∂_μ p(μ₀),粒子数方差 σ_Λ/|Λ| 趋向 β⁻²∂²_μ p(μ₀),而且这个极限值与边界条件无关。
这个结论的深刻之处在于:即便系统存在多个 Gibbs 测度(例如硬核格点系统在高逸度下不同子格对应的共存相),它们的粒子数方差仍然会收敛到同一个值。粒子总数是一个由压力决定的体积量,它不会因为边界条件的不同而"分裂"。边界条件可以移动有限体积系统的零点,但无法移动它们定义的极限。这为解决"共存相之间压缩率是否一致"这一悬而未决的问题提供了第一个一般性的肯定答案。

证明拆解:三个经典工具撑起一个优雅结论

这套定理的证明只用了三个经典复分析工具:Borel–Carathéodory不等式、Vitali定理和Cauchy积分公式。没有复杂的现代数学机器,没有繁琐的技术引理堆砌,完全是一个"大巧若拙"的证明。这也是本文最让人感到舒服的地方,龙哥看到这里都想拍一下大腿。
先来看第一个关键步骤:配分函数的一致上界。对于稳定、温控、下正则的对势,加上均匀有界密度的边界条件,可以得到一个非常干净的不等式,对任意复数 z、任一边界条件 b 和任意边长 L:
|Ξ(z; Λ_L, b)| ≤ Ξ(|z|; Λ_L, b) ≤ exp( |z| |Λ_L| e^(β(B+K)) )
其中 B 是稳定性常数(满足 φ(x₁, ..., x_N) ≥ -BN),K 是由边界密度上限和势函数下正则包络共同决定的常数。这个"一致上界"有两个作用:一是保证配分函数是整函数(对任意 z 都收敛),二是为后续应用 Borel–Carathéodory 不等式提供所需要的增长控制。
接着是核心的复分析论证。假设在圆盘 B_δ(z₀) 内没有零点,那么在单连通区域上就可以定义 log Ξ 的解析分支。这个分支在实轴上与有限体积压力一致(因为虚部在实轴上是连续的且取值于 2πZ,在 z₀ 处为零,故恒为零)。记之为 p(z; Λ_L, b)。然后应用 Borel–Carathéodory 不等式来获得一致局部有界性。
Borel–Carathéodory 不等式是复分析中的一个经典工具:它用解析函数的实部来控制系统整体的模。具体来说,如果 f 在圆盘 B_R(z₀) 上解析,那么对 0 < r < R,有 sup_{B_r(z₀)} |f| ≤ (2r/(R-r)) sup_{B_R(z₀)} Re f + (R+r)/(R-r)|f(z₀)|。这里用 Re p ≤ c₁|z| + c₂(来自可容许性条件),再加上 p(z₀) 的非负有界性,就得到有限体积压力族在 B_δ(z₀) 内局部一致有界。
有了局部一致有界性,就可以引出第二个工具——Vitali定理。Vitali定理说的是:如果一列解析函数在某个区域内局部一致有界,并且在某个具有聚点的子集上逐点收敛,那么这一列函数在整个区域内内闭一致收敛到某个解析函数。这里,"逐点收敛"来自于热力学极限的定义,在实轴区间 I 上成立;"局部一致有界"则来自上面的 Borel–Carathéodory 估计。Vitali定理于是给出了在整个圆盘 B_δ(z₀) 上的内闭一致收敛性。由于在实轴上的极限是压力 p(与边界条件无关),这个解析延拓也就自动与边界条件无关。
最后一步,把一致收敛"升级"为导数的收敛,这只需要第三个工具——Cauchy积分公式。对于任意 z ∈ B_r(z₀),取 ρ > 0 使得圆周 |ζ - z| = ρ 完全落在 B_δ(z₀) 内,则有:
∂ⁿ_z(p(z; Λ_L, b) - p(z)) = n!/(2πi) ∮ (p(ζ; Λ_L, b) - p(ζ))/(ζ - z)^(n+1) dζ
对圆周上的积分路径进行估计,立刻得到 |∂ⁿ_z(p(z; Λ_L, b) - p(z))| ≤ n!/ρⁿ sup_{B_{r+ρ}(z₀)} |p(·; Λ_L, b) - p|。右边随 L → ∞ 趋于零(因为一致收敛),从而证明了任意阶导数在 z₀ 附近一致收敛到对应的极限导数。取 r = δ/2、ρ = δ/4,就获得更强的结论:在 B_{δ/2}(z₀) 上一致收敛。
meng.jpeg
整个证明干净利落,没有任何多余的枝蔓。读完让人感叹:好的数学不需要堆砌工具,而在于找到正确的视角。

哪些系统能用?无零点区域的“白名单”

既然定理的假设是"零点保持有界距离",那么自然的问题就是:哪些系统确实具有这种性质?论文的第三节列出四个已知的无零点区域。这实际上是一个"应用白名单":只要系统落入其中任一情形,定理的全部结论自动成立。
第一,低逸度区域。当 |z| 足够小时,Mayer 团展开(cluster expansion)和 Kirkwood–Salsburg 方程给出收敛的级数表示。事实上,Ruelle、Penrose 和 Groeneveld 早在1960年代就严格证明了这一区域的收敛性,压力可以写成 βp(z) = Σ b_n zⁿ 的形式,其中 b_n 是 Mayer 团积分。在这个区域内,有限体积的粒子数方差自动收敛到热力学极限值,没有任何零点方面的担忧。
第二,铁磁自旋系统。对于 Ising 铁磁模型等系统,Lee-Yang 圆定理(Lee-Yang circle theorem)保证所有零点都落在单位圆 |z| = 1 上。因此,当 z₀ ≠ 1(等价于外磁场不为零)时,零点和实轴上的任意点都保持正距离。这个经典结论已经被 Griffiths、Simon、Lieb、Sokal 和 Newman 等人推广到更一般的铁磁相互作用。最近,Jiang 和 Newman 还证明了第一个 Lee-Yang 零点的模随系统尺寸增大而减小,并收敛到零点的分布端点,说明零点不会泄漏到更靠近原点的位置。
第三,大逸度下的硬核格点系统。这部分是最有趣的,因为它同时涉及"存在相变"和"压力解析"这两个看似矛盾的性质。以 Z² 上的硬方块模型、硬六边形模型以及二部格点上的最近邻排斥模型为代表,当逸度 z 足够大时,压力对 z⁻¹ 是解析的,但系统同时存在多个 Gibbs 测度——例如偶子格和奇子格上的相。这正是 Theorem 1.1 发挥威力的舞台:压力解析性保证了零点不贴实轴,而零点不贴实轴保证了所有导数的收敛。Pirogov–Sinai 理论为这类系统提供了系统化的处理框架。Jauslin 和 Lebowitz 的近期工作进一步把结果扩展到非滑动(non-sliding)模型,证明了在逸度平面上的一个环形邻域内不存在零点。
第四,精确可解模型。二维 Ising 模型、硬六边形模型和六顶点模型的配分函数有闭式表达式,可以在复逸度平面上完全确定奇点结构。这些模型构成了验证一般定理的良好的测试平台。

边界条件的“紧箍咒”:从均匀密度到无界密度的推广

论文对边界条件的处理非常细致,这是理论工作中容易被忽视但极为重要的环节。热力学极限的独立性是否成立,往往取决于允许哪些边界条件参与竞争。
正文部分处理的是"均匀有界密度"的边界条件集合。定义 B_ρ̄:对每个单位立方体 Δ_q,边界构型 b 在其中的粒子数 m_q(b) 不超过 ρ̄。这类边界条件在物理上相当自然——实际系统中不会出现某个位置堆积异常多的粒子。对于这类边界条件,论文证明了可容许性条件(9)成立,即配分函数沿立方体序列被一致指数上界控制。
均匀有界密度听起来已经挺合理了,但有些物理上有效的边界条件并不满足这个约束。于是论文在附录中推广到两类更宽泛的情形。
第一类是 Procacci 和 Yuhjtman 的无界密度边界条件。这类边界条件允许局部密度随位置增长,但增长速度受到一个函数 g 的控制:m_q(b) ≤ ρ̄(1 + g(|q|_∞))。如果势函数满足超稳定性(super-stability),增长条件又足够温和,那么配分函数依然有可控的上界。这里可以看到势函数和边界条件之间的权衡:想要放宽边界条件的限制,就必须加强势函数的约束。从"稳定"升级到"超稳定",相当于用更强的势条件换取更弱的边界条件限制。
第二类是 Ruelle 的温控边界条件。这类边界条件不要求局部密度有界,而是要求其平方和随体积增长的速度满足一定的限制。同样在超稳定的条件下,论文验证了可容许性条件。由于两类边界条件都覆盖了,它们的叠加情形也被隐式地包含了。
这种"势条件与边界条件之间的权衡"是经典统计力学中反复出现的主题。可容许性条件的本质,是确保配分函数的增长速度可控,从而使得有限体积压力序列足够"规矩",不会出现失控的涨落。边界条件的泛化推广不仅增加了定理的适用面,更揭示了零点条件与热力学极限之间关系的某种稳健性。

打开的问题:零点逼近速率与超均匀性

论文在最后一节明确指出了当前结果的局限性和未来方向。作者诚实地承认:定理本身并不提供零点的具体位置信息,而是把零点无区域当作一个假设来使用。在一般情况下,并没有先验的控制办法来保证零点不逼近实轴。这决定了定理的适用范围只能限定在那些已经有独立方法(如团展开、圆定理、Pirogov–Sinai理论)确认了无零点区域的系统。
论文预期,在压力解析的每一个点上,所有导数都应该收敛且与边界条件无关,不管有限体积零点以何种方式逼近实轴。少数零点漂移到一个压力仍然光滑的实轴点附近,理论上不应该造成任何影响。但要证明这一点,需要定量地控制零点逼近实轴的速率和密度。这正是目前技术无法做到的——除了少数特殊模型(例如铁磁系统第一个零点就有精确控制)。
更根本的开放性问题是:真正的超均匀态在正温平衡系统中是否存在?论文的定理明确排除了零点的作用,但超均匀性往往与长程相互作用相关(例如库仑系统),或与临界点上的零点贴实轴密切相关(此时压缩率发散归零意味着二阶导数消失,必然伴随零点抵达实轴)。在这些情形下,热力学极限和求导的交换需要完全不同的机制。这篇论文只是把边界往前推了一大步,但尚未抵达终点。

龙迷三问

下面是龙哥对大家可能疑问的解答:

这篇论文到底解决了一个什么问题?它证明了:在Lee-Yang零点与正实轴保持正距离的条件下,有限体积压力的任意阶导数在热力学极限下收敛到相应的极限导数,并且极限值与边界条件无关。换句话说,只要零点不捣乱,粒子数密度和粒子数涨落就会收敛到由热力学压力唯一确定的数值——即使系统存在多个共存相也是如此。这为"共存相之间的压缩率是否一致"这个问题给出了第一个一般性答案。

为什么一阶导数不需要零点条件,二阶却需要?因为有限体积压力在化学势 μ 上是凸函数。凸函数序列如果逐点收敛,那么它们的一阶导数在极限函数可导的地方自动收敛(这是凸分析的标准结论)。一阶导数就是粒子数密度,所以密度在远离一级相变处自然收敛。但凸性对二阶导数没有同样的控制能力——它只能保证二阶导数为正的测度,无法约束其极限行为。二阶导数对应粒子的涨落,这正是需要零点条件介入的地方。

超均匀性和刚性定理之间是什么关系?Ginibre(1967)对有限体积建立了压缩率的下界,所以任何具有稳定相互作用且对势非负的系统,在正温下涨落必然是广延的,不可能超均匀。Dereudre 和 Flimmel 把这个下界推广到了超稳定的 Gibbs 点过程。本篇论文则关注另一个方向:不仅证明涨落是广延的,还要证明涨落·每单位体积的极限值是唯一的——不依赖于边界条件。这两个方向互补,合在一起描绘出更完整的图景。

如果你还有哪些想要了解的,欢迎在评论区留言或者讨论~

龙哥点评

论文创新性:★★★★☆ 新定理把Lee-Yang零点与粒子数涨落严格联系起来,方向新颖,与历史文献的对话充分。

用经典复分析工具证明一个此前未知的结论,可以说是在老树上开出了新花。

实验合理度:★★★★☆ 纯理论证明没有实验环节,理论的内部自洽性和逻辑链条的严密性是评判依据。

从引理到定理的推导过程严谨,三个附录分别处理了不同边界条件,结构完整。

学术研究价值:★★★★☆ 对统计力学的严格化有直接推动作用,特别是对共存相涨落唯一性问题。

理论上为后续研究零点分布与导数收敛的关系提供了新框架。

稳定性:★★★★☆ 数学证明在逻辑上完全稳定。但在实际物理系统中的适用性依赖于零点条件的验证,这是暂时无法回避的局限。

适应性以及泛化能力:★★★★☆ 适用于低逸度、铁磁自旋、硬核格点大逸度、精确可解模型四类系统。但超出这些类别时零点条件难以验证。边界条件方面覆盖了三类:均匀有界密度、Procacci–Yuhjtman 无界密度、Ruelle 温控边界条件。

硬件需求及成本:★★★★★ 纯数学推导,无需任何硬件,理论研究零成本。

复现难度:★★★★☆ 证明步骤清晰完整,使用的工具都是复分析中的经典定理,只需扎实的数学基础即可逐步验证。论文中的附录提供了完整的推导。

产品化成熟度:★☆☆☆☆ 纯基础理论研究,无直接产品化目标。其价值体现在为统计力学的数学基础添砖加瓦。

可能的问题:定理的核心假设——零点保持有界距离——本身是难以验证的。与一些可以自动验证条件的现代统计物理工作相比,这篇论文把零点位置当作输入而不是结果,在一定程度上限制了结论的适用范围。此外,证明依赖于凸性和解析延拓,这些方法不太可能直接推广到动态过程或非平衡系统。


主要参考文献

[1] D. Ruelle. Statistical Mechanics: Rigorous Results. W. A. Benjamin, New York, 1969.
[2] T. D. Lee and C. N. Yang. Statistical theory of equations of state and phase transitions. II. Lattice gas and Ising model. Physical Review, 87(3):410–419, 1952.
[3] J. Ginibre. Rigorous lower bound on the compressibility of a classical system. Physics Letters A, 24:223–224, 1967.
[4] D. Dereudre and D. Flimmel. Non-hyperuniformity of Gibbs point processes with short-range interactions. Journal of Applied Probability, 2024.
[5] A. Procacci and S. A. Yuhjtman. Classical particles in the continuum subjected to high density boundary conditions. Annales Henri Poincaré, 23:799–827, 2022.
[6] D. Ruelle. Superstable interactions in classical statistical mechanics. Communications in Mathematical Physics, 18(2):127–159, 1970.
[7] H.-O. Georgii. Gibbs Measures and Phase Transitions, volume 9 of de Gruyter Studies in Mathematics. Walter de Gruyter, Berlin, 2nd edition, 2011.
[8] S. Torquato and F. H. Stillinger. Local density fluctuations, hyperuniformity, and order metrics. Physical Review E, 68(4):041113, 2003.
[9] 原论文链接:https://arxiv.org/pdf/2607.26975v1.pdf

*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击"阅读原文",查看更多原论文细节哦!       

end
Lee-Yang零点离实轴有多远,相变的秘密就藏在里面。想和龙哥一起追踪统计物理与数学物理的前沿进展?欢迎加入龙哥读论文粉丝群,扫描下方二维码或者添加龙哥助手微信号加群:kangjinlonghelper。一定要备注:研究方向+地点+学校/公司+昵称,根据格式备注,可更快被通过且邀请进群。
『龙哥读论文』微信群目前包含:图像处理、大模型及智能体、自动驾驶及机器人、AI医疗及AI金融5个群
wechat_helper dianzan
转发文章 微博 X LinkedIn Facebook
龙哥读论文 · PaperDaily

本文基于龙哥读论文 PaperDaily 数据库整理,结合论文原文与工程视角进行解读。