邓煜 Fields Medal 解读 —— 为物理学家准备的 Hilbert 第六问题与 Boltzmann 方程的数学视角研究导览

Hilbert 的野望

1900 年 8 月,巴黎召开第二届国际数学家大会,Hilbert 在会上作了那场后来广为人知的演讲,开列出一批他认为将牵引二十世纪数学走向的问题。这批问题——写进正式文本时共二十三个——绝大多数是纯数学的家事:连续统假设、Riemann 猜想、Diophantus 方程的可解性等。唯独排在第六位的那个显得格格不入,因为它谈的不是数学内部,而是物理:如何用公理化的方式处理物理学。

Hilbert 的底气来自前一年。1899 年他刚出版《几何基础》,把欧几里得两千年来靠直觉支撑的体系重新架在一组清晰的公理之上。他要问的是:既然几何能这样做,那些数学已经深度介入的物理理论——他点名了概率论与力学——为什么不能照办?第六问题字面上的主张,就是要给物理学也找出这样一套地基。

不过对我们要讲的故事而言,真正重要的是他埋在这段话里的一个具体得多的要求。Hilbert 特别提到 Boltzmann 关于力学原理的工作,说

其中「仅仅被示意」的那些极限过程值得从数学上认真发展起来——那是一些「从原子论的观点导出连续介质运动定律」的极限过程。

这句话对物理学家而言稍微有些晦涩难懂。如果翻译成物理学家喜欢的话,大意是:一团气体,底下是数不清的分子按 Newton 定律你推我搡;上面我们却用密度、速度、温度这些光滑的场,拿 Euler 方程或 Navier–Stokes 方程去描述它。两幅图像都对,也都在用,可从下面那幅怎么严格地走到上面这幅,Hilbert 认为还是一笔糊涂账,而这道坎该由数学来迈过去。

说它是「一笔糊涂账」并不夸张,因为这里并排站着的其实是同一团气体的三层描述。最底层是微观力学:约 个分子,每个都是一颗遵守 Newton 方程的质点或硬球,确定、可逆,原则上完备,实际上没法直接用。最上层是宏观流体:几个连续场满足的偏微分方程,就是工程师和气象学家每天在解的那一套。夹在中间的是 Boltzmann 在 1872 年给出的动理学描述——一个关于分布函数 的方程,用统计的语言回答「某处、某速度附近有多少分子」。物理学家在这三层之间穿梭得非常自然:谁都知道气体由分子组成,谁也都会写 Navier–Stokes。但「两头都知道」和「从一头严格推出另一头」是两回事,中间靠的是一串极限——分子数趋于无穷、平均自由程趋于零、时空尺度放大——把这些极限做成真正的数学定理,证明微观的解确实收敛到介观、宏观的解,才是 Hilbert 那句话真正的分量。

这样一份「把物理学界早已通行的推导补成数学定理」的清单,乍看近乎多此一举,实则不然:这条链子里从 Boltzmann 时代起就卡着一个死结——可逆的微观力学,怎么会生出不可逆的宏观规律?1872 年 H 定理刚给出熵只增不减,Loschmidt 就反问:把所有分子的速度反号,轨道原样倒退、熵反而下降,单调性又从何谈起?Boltzmann 的回答是「分子混沌假设」——认定相撞的两个分子在碰撞前互不关联——它巧妙地把时间之箭塞进了统计描述,却把一个数学问题原封留在那儿:这个假设本身凭什么成立、又从哪里来。

再往下追,物理学家和数学家就分道扬镳了。前者把答案一路推到时间之箭和宇宙的低熵初态,认为不可逆性是从边界条件继承来的;后者不碰宇宙学,只问能否从可逆的多体动力学本身、在恰当的极限下把分子混沌严格地推出来——这正是 Hilbert 第六问题的数学化身,也是本文的主线。个中曲折留待后文;这里只需记住,要看懂邓–Hani–马究竟做了什么,得先回到那个被反复追问的中间层:Boltzmann 方程本身。

时间会倒流吗?

志得意满的 Boltzmann

1872 年,Boltzmann 完成了他一生中最为满意的工作,也就是被后世称为 Boltzmann 积分微分方程的那个长长的等式。它描述了稀薄气体中分布函数 的演化:

这里头和位于,速度是 的气体分子碰撞的是带星号的那些(),而碰撞后的量都带上撇号;散射后,带星号的这些粒子进入立体角 , 表示的是 这些粒子的分布,其他以此类推;要考虑所有这些过程的贡献,就是要对所有的散射方向和碰撞例子积分。正因如此,这个方程左边是微分、右边是积分,还把三维位形空间和三维速度空间搅在一起,即便 Boltzmann 本人也求不出它的一般解。可真正让他志得意满的,从来不是解这个方程,而是不必解就能从它读出的两件事。

其一,它自己认得平衡态。方程右边那个积分——碰撞项——记的是碰撞如何持续改写速度分布。如果没有外力,而又密度均匀,要让分布不再变,只需让右边这个积分处处归零;而由被积式可见,这要求任何一次碰撞的前后都满足 。两边取对数,就是说 在碰撞中守恒。可一次弹性碰撞真正守恒的只有粒子数、动量和动能三样,于是 只能是它们的线性组合,分布被逼成一个速度的高斯:

平衡

这正是 Maxwell 几年前靠对称性论证猜出的那条速度分布。Boltzmann 不需要任何额外假设,方程就把它当作「碰撞项的零点」自动吐了出来——平衡态无非是动理学演化的不动点。

其二,它自带一支只朝一个方向的箭头。定义一个泛函

Ɦ

Boltzmann 证明,只要 按他的方程演化,就恒有

Ɦ

而等号当且仅当分布退化成上面那条 Maxwell 分布时才成立。换成物理学家熟悉的话:Ɦ 差一个负的玻尔兹曼常数因子就是熵,Ɦ,于是「Ɦ 只减不增」就是「熵只增不减」——热力学第二定律。一条原本只能凭经验断言,如今竟从分子的碰撞里被推了出来。这就是 H 定理。只是这份证明里悄悄嵌着一个核心假设,而那,正是日后所有麻烦的入口。

对 Boltzmann 而言,这几乎是热力学的盖棺定论:熵增不再是一条从天而降的经验律,而是分子在 Newton 力学下彼此碰撞的必然归宿。他有理由志得意满——事实上,若不是同时代还并肩站着一个 Maxwell,他几乎必然会被视为 19 世纪下半叶最伟大的物理学家。

然而这份辉煌并没能陪他走到最后。Boltzmann 的后半生几乎都耗在为原子论辩护上——与以 Mach、Ostwald 为首、坚信唯有能量才是实在的「唯能论」缠斗不休,直到 1906 年在抑郁中自尽。而写下这个方程的 1872 年,离那场悲剧尚有 34 年之遥;给这套志得意满的理论带来第一次劫难的,也不是他的论敌,而是他的密友 Loschmidt。

来自 Loschmidt 的质疑

1876 年,Loschmidt 给 Boltzmann 写了一封信,提出一个后来被称作「可逆性佯谬」的反驳。它简单却无懈可击:Newton 方程对时间反演是对称的——若 是一段合法的分子运动,那么把每个分子在某一刻的速度统统反号、让它们沿原路倒行, 同样合法。于是取任一段熵在升、Ɦ 在降的演化,到终点把所有速度一反,系统便精确地沿原轨迹退回,Ɦ 一路回升、熵不增反降。既然每条熵增的轨道都配着一条镜像的熵减轨道,「熵永远只增不减」就绝无可能是纯从可逆力学里推出的定理——它对那条反向轨道本该同样成立,可那里熵分明在掉。

Loschmidt 戳中的,正是上一节末尾埋下的那个「核心假设」。H 定理的证明里藏着一步看似无害的手脚:算碰撞项时,Boltzmann 默认即将相撞的两个分子在碰撞前彼此独立,成对的分布 1 可以拆成两个单分子分布的乘积,

这就是「分子混沌假设」(德文 Stoßzahlansatz,字面即「碰撞数假定」)。它看着像常识,实则极不对称:它只对碰撞前成立。一次碰撞过后,两个分子的速度已被这次碰撞绑在一起、带上了关联;你若此刻反演时间,这些出射关联就变成入射关联,分子混沌当场失效。也就是说,这个假设本身就挑好了一个时间方向——是它,而不是 Newton 力学,把那支箭头塞了进来。

佯谬于是化解:H 定理从不断言「任意微观态的熵都在增」,它只断言「满足分子混沌的态的熵在增」。而 Loschmidt 反演出的那个态,恰是分子混沌被破坏得最彻底的态——它每一对分子都被预先精细关联好,专等着「反着撞回去」把熵送低。这种态在所有微观态里是极罕见的例外,可以合法地排除在定理前提之外。重点从来不是它存在,而是它有多重要——这一点,要到后面用相空间的「体积」才算得清。

再往深一层,这一步暴露了 Boltzmann 方程的真身:它不是 Newton 方程的恒等变形,而是一个统计方程。分子混沌是一道粗粒度的、概率性的输入——它至多在「绝大多数」微观态上、在分子数趋于无穷的极限里近似为真,绝不是对每条力学轨道都成立的机械事实。Boltzmann 到底有没有说服 Loschmidt,已无从查考;但他显然没能说服自己给这道假设一个干净的交代。「凭什么分子混沌对典型的态成立」——这个被他一句带过的问题,将耗去此后一百多年,也正是本文最终要回到的地方。

Ehrenfest 的狗和跳蚤

Boltzmann 没能等到替他翻案的那一天。1906 年 9 月,他在的里雅斯特附近的度假地自缢;而就在前一年,Einstein 已用一纸关于布朗运动的论文,把「看不见的分子」折算成显微镜下悬浮颗粒可测的抖动——原子论最有力的那块实验拼图正在合拢,Perrin 的验证也只在两三年之后。假若他再多撑几年,或许一切都会不同;然而历史不容倒带,时代的车轮径直碾了过去。

接过这套理论的,是他的高足 Paul Ehrenfest。1907 年,Ehrenfest 夫妇写下一篇题为《关于 Boltzmann H 定理的两个著名反驳》的短文,用一个玩具模型把 Loschmidt 的诘难连同另一个更棘手的反驳一并讲得透亮;四年后,他们又为数学百科全书写就那篇统计力学基础的经典综述。可命运待这一支同样刻薄:1933 年,Ehrenfest 也在抑郁中结束了自己的生命。为这门学问奠基的两个人以同样的方式离场,或许并非全是偶然——关于「不可逆究竟从何而来」的追问,此后也一直没能挪到物理学舞台的正中央。

那个玩具模型,后来被昵称为「狗与跳蚤」。设想两条狗 A、B,身上共有 只编了号的跳蚤。每一步随机报出一个 到 的号码,被点到的跳蚤就从当前那条狗跳到另一条身上——规则至此说完。有两点要紧:其一,这套动力学严格可逆,把上一步的号码再报一次,跳蚤就跳回原位,整个过程倒放同样合法,Loschmidt 的时间对称性一分不差;其二,我们真正盯着的不是哪只跳蚤在哪儿,而是一个宏观量—— A 身上的跳蚤数 。

假设从这样的初态出发: 只跳蚤全趴在 A 上()。接下来不必模拟就能断定——绝大部分时候 都在往下掉,直到在半数 附近落定、然后微小起伏。统计上说,这是因为 A 上恰有 只跳蚤的排布数是 ,它在 处尖锐地堆成一座高峰。用 一读,贴近半分的宏观态「体积」大得压倒一切,远离半分的态则寥寥可数。于是从一个稀有的偏态出发,随机一跳几乎必然把系统推向态数更多的那一侧——这就是熵增,不多不少,全是概率。

那 Loschmidt 呢?那条熵减的倒放轨迹当然存在,但它要你恰好从一个「排布精巧、专等着跳回全满」的态起步,概率小到可以抹去。另有一个更凶的反驳来自 Zermelo:Poincaré 证明过,这样的有限系统迟早会回到任意接近初态之处——跳蚤总有一刻重新聚回 。这一条模型也答了:回归确会发生,但平均要等上约 步;对一撮真气体的 ,这个数远远盖过宇宙的寿命,于是「必然复现」在任何可观测的意义上都等同于「永不复现」。可逆又可复现的微观规则,就这样心安理得地托起了一条单向、不回头的宏观曲线——所谓不可逆,自始至终只是压倒性的概率。

高熵的海洋

Ehrenfest 的比方虽妙,却终归是一条人为设定的 Markov 链——每一步「随机报号」在真实气体里并无对应,分子不掷骰子,只死板地沿 Newton 方程走。从这样一个掷骰子的模型到真正的哈密顿动力学,中间隔着山河大海。但它点出的内核是对的:微观可逆与宏观不可逆并不矛盾,甚至能安放进同一个框架。

而这个内核根本不靠那点人造的随机性——把随机性彻底抽掉,结论照样成立。1956 年,Mark Kac 搭了一个如今叫「Kac 环」的模型: 个格点连成一圈,每个格点上坐着一颗非黑即白的球;圈上有 条边被预先钉了「记号」,位置一经指定便不再变。时钟每走一格,所有球一齐沿环前进一步,凡跨过带记号的边就翻色、跨过干净的边则不变——规则至此说完。它是确定的,一切都由固定的记号分布与初始颜色定死;又是严格可逆的,倒着放(逆行一格、遇记号照翻)便精确还原。可就在这么一个死板可逆的系统里,宏观量竟会「弛豫」:记白球比黑球多出的净数为 ,只要假定每颗球前方是否撞上记号彼此无关、比例恒为 (正是分子混沌的翻版),一步就给 乘上一个绝对值小于 的因子,到第 步的时候,平均而言,这个净数 指数式地衰向零——黑白各半的平衡。一个货真价实、可逆系统里的 H 定理。

当然,精确的动力学并不会真的一路衰减到底。转满 步,每颗球绕环两周、每条边恰好各跨两次,翻色成双成对地抵消,整个系统一丝不差地弹回初态——它不但复现,还快得离谱:区区 步,远不是跳蚤那 的天文数字。破绽正在此处:那条单调衰减从来不是精确动力学的定理,而是借了「彼此无关」这条分子混沌假设;它只在关联尚未累积的 窗口里成立,一旦球转够了圈、关联反噬,复现便赢了回去。于是连这个从头到尾都能验算的玩具,也只在复述同一句话:可逆的微观、加一条统计假设、再加压倒性的数目,就足以撑起一段看似不可逆的宏观弛豫。剩下的唯一疑问是——凭什么那条统计假设对「典型」的初态成立?

Kac 环的一步演化: 个格点各坐一颗球(白或黑),红杠标出钉了记号的边。每一步所有球顺时针挪一格,跨过记号即翻色。图中这一步, 个记号各把一颗经过的白球翻黑(橙圈标出),另两颗原有的黑球只随之平移、不翻色;净色差 由 降到 ,向黑白各半靠拢。整套规则确定又可逆,转满 步后还会精确复位。

早在 1880 年代,Boltzmann 本人就给了一幅干净的图。把 个分子每一种可能的微观状态——每个分子的位置连同速度——都设想成一个点,它们撑起一个 维的辽阔空间,学过统计力学的人会脱口而出它的名字:Γ 空间。系统此刻的状态是其中一个点,随时间演化便扫出一条轨迹;这条轨迹本身当然可逆,把所有速度轴反号不过是给 Γ 空间照一次镜子,物理内容分毫不改。真正决定命运的是另一件事——延续跳蚤模型的算法,只需把那里的「排布数」换成连续相空间里的「体积」,Boltzmann 熵便是 2。显而易见,「平衡、高熵」那一大类宏观态占的体积大得没有对手,「远离平衡、低熵」的态则挤在一个可怜的角落;两者之比往往是 的 次方量级——对 ,这是一个大到失去日常意义的悬殊。

于是结论近乎几何式的必然:你若不巧被搁在一座低熵的小岛上,下一刻脚下几乎必是茫茫大海——不是有什么力推着你离岸,纯粹因为海太大、岛太小,随便迈一步都落进海里。(相流还保体积,这是刘维尔定理,所以这并非谁把小岛压扁了,而只是「可去之处」几乎全是高熵。)而一旦入海,你要在这片浩瀚里游荡到另一座小岛,所需的时间长得穷尽宇宙之寿也远远不够。这,就是我们看到的熵增:它不是一条写在运动方程里的定律,而是一桩关于「体积」的压倒性事实,与微观可逆没有半分冲突——你尽可以逆着轨迹走回小岛,只是这样的初态稀有到可以当作不存在。

高熵之海:把 个分子所有可能的微观态摊开,得到 维的 Γ 空间;「远离平衡、低熵」的那一大类宏观态在其中不过是一座孤岛,四周尽是辽阔得多的「平衡、高熵」之海。系统的状态沿一条轨迹漂流——轨迹本身完全可逆,也没有谁在推着它离岸;只是可去之处几乎全在海里,于是从岛上迈出一步便几乎必然入海,而要在浩瀚中游荡到重新踏上小岛,所需的时间远远超过宇宙的寿命。

时间的箭头

高熵的海洋那幅图足够有说服力,物理学家于是顺势追问:那整个宇宙呢?我们之所以总能亲眼见到熵由低走高,是因为我们总能亲手制备一个低熵的孤立系统,再看它自发地爬向高熵。可这桩买卖有代价——它在消耗低熵:要让一角变得有序,就得把更大的一片搅乱,每一次「亲眼见到第二定律」,其实都在支取一笔低熵的存款。而低熵是这世上最结实的不可再生资源,比环保主义者口中的石油煤炭还要「不可再生」得多——那是字面意义上的不可再生,无论科技怎么发展都造不出净的低熵,你至多把它从一处挪到另一处、且每挪一次都净亏。于是真正的问题浮出水面:我们正大手大脚花着的这笔庞大低熵,最初究竟从哪儿来?答案只有一个——宇宙诞生于一个熵低到离谱的初态,而所谓「时间之箭」,无非是这笔遗产被单向花掉的漫长过程。

以上这套推理,就是大名鼎鼎的「宇宙学箭头」,绝大多数物理学家对它相当受用。妙处有二。其一,它把不可逆彻底交还给了我们自己:微观世界一丝一毫都可逆,没有哪条方程偏爱哪个时间方向;我们之所以看见不可逆,不过是因为一副仅有千亿量级神经元的大脑,根本盯不住 量级的粒子,只好退而用粗粒度的眼光去看——是这层近似、而非世界本身,生出了「不可逆」。其二,正因为宇宙的开端足够低熵,我们才有本钱一次次制备低熵、再看它奔向高熵;若非如此,我们打出生起就泡在一片一览无余的热寂里,什么箭头、什么故事,都无从谈起。

可这套宇宙学箭头有个薄弱处。回过去看那片高熵之海——就在它与低熵小岛的交界线上,坐着一批特殊的高熵态——它们的粒子之间编码着复杂的关联;顺着动力学往下走,它们会轻而易举地滑回岛上、滑向低熵。这并非空想:想象一桌子的台球——既然你能一杆把码成三角形的 15 粒球打得四散奔逃,那反过来,这 15 粒球只要位置和速度都恰到好处,当然会一颗颗倒退着重新排回那个三角形。其实我们刚才已经见过了:Kac 环里那 步之后的精确复位,正是关联一点点攒足、终于把整个系统拽回初态的现场。同样的逻辑还能再狠一层——哪怕我们此刻就待在低熵态里,只要凑够高阶关联,照样能滑向更低的熵。可现实中谁也没见过熵自个儿往下掉,为什么?宇宙学箭头再一次赶来救场:你根本没有那么多关联可花。这和低熵的道理如出一辙,只不过这次见底的是另一种资源——关联;宇宙给的初始存货里能拿来逆天而行的关联无比稀薄,不够你用来构建这个初态。反过来说,如果你出生在一个充满了精巧的关联的宇宙中,那,那边的时间箭头也许是反过来的!

把一杆开球倒着放:15 粒台球从满桌乱撞中一颗颗收拢,末了严丝合缝地扣回那个三角形,母球退回杆头线。为「归位」并没有另行设计什么——这段画面是先正着算出一次开球,再把每一帧倒序播放,故它逐帧都严格满足 Newton 方程,是一个货真价实的合法解(图中不计摩擦,故球永不停歇)。它之所以从未在哪家球房里上演,只因要把 15 组位置与速度凑得这般恰到好处,需要一份精巧到离谱的关联——而分子不是 15 个,是 个。

对绝大部分物理学家,整个故事就已经可以结束了,可是总是有人依然不服气儿:你只说了关联不是来自初态的,但凭什么关联不会在演化里自己悄悄长出来?——别忘了,Kac 环里它偏偏就长出来了!

退一步看,这一切的正中央,始终立着 Boltzmann 那一步「分子混沌」。推导方程时,他每算一次碰撞,就默认相撞的两个分子在碰撞前互不关联、把成对分布重新拆回单分子分布的乘积——说白了,就是每一拍都把刚刚攒起的关联清零、重新注入一次「无关联」。这跟宏观上那层粗粒度近似根本是同一件事:都在把关联扔掉。于是所有的问号最终拧成三个:这样一次次扔掉关联,凭什么是对的?为什么如此粗暴的一步近似,换来的 Boltzmann 方程却精确到足以管住真实气体、误差小得可以忽略?而这被忽略的误差,又与那些被扔掉的关联究竟是什么关系?要把这三问问得清清楚楚、进而答得严严实实,物理学家那点直觉已经不够使了,得请出一套专为「关联」记账的机器——一串把单分子分布、双分子分布、三分子分布……层层咬合起来的方程,一个由五个人姓氏首字母拼成的名字:BBGKY。

什么是关联?

万恶之源:Liouville

要真正解开「关联是什么」这个谜,得把时钟往回拨上大约一百年——早在统计力学诞生以前,一个叫 Liouville 的法国数学家,就在纯粹的经典力学里写下了那条日后大名鼎鼎的 Liouville 方程。他问的是个再朴素不过的问题:在相空间里撒下一团概率密度,任它随力学演化,这团密度会怎么变?对我们这团 个分子的气体,这密度就是活在 维 Γ 空间里的联合分布

它一次记下所有分子的位置和速度,意思是说,对这样一个系统,它的第一个粒子在 ,速度是 ,第二个粒子在 ,速度是 ,等等等等,的概率密度,在时刻 是多少呢?这是一个有着 个自变量的极高维的函数,下文为省笔墨,宗量不再写全,但下标 一直保留。至于给 头上加的这一飘(),是一句关键的提醒:它是归一化到 的概率密度——和 Boltzmann 方程里那个「数分子」、整团积出来是 的 并不是同一个东西;两者究竟差多少,到下一节自会算清。

Liouville 的推理,反过来从守恒出发,只有两步。第一步,写下连续性方程。相点既不生也不灭,密度 就像一种守恒的流体——它随时间的变化,全由「流进流出」的净通量定。对一般的流体,守恒律方程说:,而流密度恰为 。在相空间中,流速沿第 个粒子第 根速度和位置轴的分量,正是 (位置的变化率)与 (速度的变化率),因而概率流便是 与 ;对全部 个分量求散度,连续性方程便是

这一步还没碰任何力学,只是概率守恒。

第二步,把 Hamilton 正则方程代进去。拿乘积求导,对任意一对 、 上式每一项都裂成两半(现在可以省略每一项自己的标记 了):一半是 、(密度被流带着走),另一半是 (流本身的胀缩)。关键在后一半,而它恒等于零:这是由于,正则方程逐分量地说

质量是常数, 对 的偏导与 对 的偏导分子分母同差一个 、恰好约掉,于是胀缩项里的每一对都是

两个二阶混合偏导恒等相消——对任何 Hamiltonian 都成立,压根不用写出它的具体形式。这就是相流保容,也是全部玄机所在。

于是刘维尔定理自己冒了出来。后一半既已消失,连续性方程只剩「密度被流带着走」那一半,正好凑成密度沿轨迹的全导数:

密度沿每一条轨迹一分不改——这,就是 Liouville 方程。

请务必看清:它什么新东西都没做。它逐字逐句就是那 条 Newton 方程(、),只不过改用「一整团初始条件」的密度来讲;没添半条假设、没丢一丝信息——换句话说,到此为止我们一丁点粗粒化都没有做,一根关联都没有扔。而这份诚实的代价立竿见影:正因为什么都没扔,Liouville 方程也就给不出任何不同于 Newton 的东西。它和 Newton 方程一样完全可逆——「相流保容」正是我们在高熵之海里用过的那条 Liouville 定理;密度沿轨迹一分不改,靠 整体算出的那份「精细的熵」便一步也不挪。这里面没有熵增,没有热力学,自然也没有时间箭头、没有不可逆,一切如同 Newton 力学本身一样干净美好。

动理学五剑客:BBGKY

麻烦只有一个: 是 维的怪物,跟直接解 Newton 方程一样动不了。于是退而求其次,只留 个分子的信息、把其余的积分掉——定义 粒子约化分布

也就是把第 到第 个分子的位置速度统统积掉,只剩前 个的联合分布。留意这一步只是「边缘化」,归一化纹丝不动: 积出来照样是 。比如全部积到只剩 1 号,得到的 就是「1 号那个特定分子在哪儿」的概率分布;当然也可以积掉除 2 号之外的所有分子,得到 2 号自己的 ;以此类推, 个分子人手一份。而 Boltzmann 那个「数分子」的分布,问的是「那儿有多少个分子」——谁坐在那儿都算数——把这 份统统加起来就是它:。再简单一点:若假设每个分子自己的分布都长得一样(注意,在经典力学里这不是白给的——量子力学有粒子全同性兜底,经典力学没有,这是实打实的一条假设),便可干脆写成 3。无论如何,要做的事情已经一目了然:从 Liouville 方程出发,推出 Boltzmann 方程。

思路已经定了,就是把 2 到 号分子全部积分,看看剩下来的 听谁指挥。先把家伙事儿摆齐:设无外场,分子只靠两两之间的力互相作用;这里考虑的典型情形,是多体相互作用能量写成 ,于是 只依赖位置、不依赖速度;本段把 暂记作 ,把「积掉 2 到 号」的测度暂记作 。把 Liouville 方程按粒子归堆,1 号一堆、其余一堆:

其中 是 1 号受的合力( 是 号施给它的那一份)。现在对整条方程做 ,逐项看它们的下场。

头两项直接穿过积分号。时间导数与积分次序可换, 也只碰 1 号自己的变量:

求和里的项全军覆没。对任何 :飞行项的系数 不含 ,力 不含 4——于是都能整个塞进导数号里,变成「对被积变量的全导数」:

全导数的积分只剩边界值,而密度在无穷远处衰减为零——这一大片求和一笔勾销。

只有 1 号受的力赖着不走。把合力拆开逐项看:

每一项里, 可以提出积分号,但 依赖被积的 ,这一重积分动不了;能动的是其余 重——把除 号之外的分子全积掉,按定义剩下的恰是 1 号与 号的联合分布:

果然——不知道二体联合分布,1 号的演化就没法闭合:光有 ,不够用。这时上一段那条「人人一样」的假设再次上岗(初值对编号交换对称,而对称的动力学会把这份对称一路保持下去): 个施力者谁来都一个样,这 项一模一样,并成 倍的一项。把幸存者们凑到一起:

左边是 1 号分子孤零零的自由飞行;右边那个积分,俨然一副「碰撞项」的长相——连那个 都在提醒你,它数的是「可能来撞我的分子」——可它根本算不动,除非你先知道 。

你不甘心,转头去给 写方程:把 3 到 号统统积掉,同样的检点再来一遍。这一回的结构更有看头:1 号与 2 号彼此之间的相互作用完整地留在了左边——两人小组按精确的二体力学演化,方才右边那位「碰撞」的本尊,如今就坐在组内;而右边挂出的,是组内每人与「组外闯进来的一个分子」的相互作用,因子变成 ,积分里坐着 。到这里规律已经一目了然:任何 人小组都按组内力学精确演化,唯一的修正来自组外闯入的那一个分子——于是每一层都向高一层索取信息,环环相扣,直到顶层那个谁也解不动的 为止:

这一整列首尾咬合的方程,就是 BBGKY 层级,名字取自上世纪四五十年代先后独立写下它的五个人——Bogoliubov、Born、Green、Kirkwood、Yvon5。它仍旧只是 Liouville 方程的另一种写法,一样精确、一样可逆;关联,就明明白白地住在 , (及更高阶)里,随层级一层层往上耦合。

快刀斩乱麻:Boltzmann 归来

要从这一整列里榨出一个只关于 的封闭方程,就必须在最低一层斩断层级。上一节末尾那条第一层方程,意思直白得很:1 号分子若不和别人碰,就只是自由飞行;真正改写它速度分布的,是右边那一坨 。所以问题并没有消失,只是换了名字——要关上第一层的门,就得知道两个分子凑近时的联合分布长什么样。(下面这一路都把外场关掉:没有外场,分子只靠两两相互作用散射;外场要加回去,也只会回到左边的输运项里,不改碰撞项这场戏的本质。)

第一件要看清的事,是「远处的分解」帮不上忙。若 远大于相互作用的力程 𝒹,两个分子当然可以视作素不相识6:

𝒹

可第一层方程右边偏偏乘着 ,而力只在 𝒹 的碰撞区里才不为零。真正要紧的,从来不是两分子相隔很远时能不能分解,而是它们即将进碰撞区、正在碰撞区、刚刚出碰撞区时, 由谁决定——恰恰是两人纠缠最深、乘积近似最不成立的地方。这把刀要砍,也只能砍在碰撞区的边界上。

第二件事,是稀薄气体递来一个小参数。把第二层 BBGKY 方程写出来:左边是 1、2 两分子彼此之间的精确二体力学,右边挂着第三个分子闯入的三体项。三体项什么时候有贡献?只有第 3 个分子也恰好钻进那个大小约 𝒹 的碰撞区时才有。于是它相对二体项被一个无量纲因子压低:

𝒹

是气体体积。对寻常稀薄气体,这个数远小于 ;物理图像也很清楚——两颗分子正贴在一起短程相互作用时,第三颗刚好也挤进同一小团空间,机会实在渺茫。于是第二层方程在领头阶上可以就地截断,只留 1、2 两人自己的散射:

这条式子已经把关键说出来了:在一次二体碰撞的时间尺度上,二体分布并非任意变化,而是被二体散射轨迹拖着走。

第三件事,是把快慢分开。换成质心与相对变量,这句话能说得更干净:

质心 与总速度只在平均自由程、平均自由时间那类慢尺度上晃动;真正在碰撞区里急剧变化的,是相对变量 与 。而碰撞历时不过

碰𝒹

( 是典型相对速度)。只要 碰 远小于平均自由时间 自由,那么在碰撞那一瞬,、质心输运、以及一切更慢的宏观变化都可以先冻结,剩下的只是相对坐标上的定态二体散射:

这可不是说两个分子没有相互作用;恰恰相反,它说的是相互作用主宰了碰撞区内部的演化——只不过这演化是精确的二体散射,于是 沿散射轨迹守恒:入口处给它什么值,二体力学就原封不动地把这个值搬到出口。碰撞区中间的细节从此不必逐点去求,只需盯住「流进多少」和「流出多少」。

于是回到第一层方程的右边。碰撞区很小, 在空间上又已被粗粒化到不分辨 𝒹 的尺度,索性把所有一体分布都看作在同一个 处取值。沿相对速度方向取一根碰撞轴 ,垂直平面上用瞄准参数 𝒷 标记入射束,碰撞区的入口、出口分别记作 、。相对速度穿过碰撞管的通量给出

碰𝒷𝒷

这是整个推导里最实在的一步:原先那个「力乘速度梯度」的体积分,被二体散射方程改写成了碰撞管两端的通量差。入口少了多少、出口多了多少,就是碰撞项里的损失与增益。

接下来只差把 𝒷𝒷 翻译成散射角的语言。固定入射相对速度 之后,二体力学自会把瞄准参数 𝒷 与散射角 钩连起来,于是

𝒷𝒷𝒷𝒷𝒷

这就是微分截面。它的妙处在于把具体相互作用势的全部细节都封装了起来:你若知道两分子之间的力,就能从经典散射算出 𝒷 与 的关系;你若只关心 Boltzmann 方程长什么样,那就只需把这份几何因子记成 便可。譬如半径为 𝓇 的硬球——两球心一旦靠到 𝒹𝓇 便算接触,于是

𝓇总𝓇𝒹

正好是半径 𝒹 的那个圆的面积——两球心擦肩而过的最大偏移恰是 𝒹,多一分就撞不上。所以「截面」这名字不是比喻,它真的在数:入射束里哪一片面积会被散射到某个立体角中去。留意这里的长度是两个半径之和而非单个分子的大小:碰撞是两颗球的事,不是一颗球的事。

现在,才轮到分子混沌正式登场。入口 处,两颗分子尚未相遇,取它们互不关联:

出口 处,两颗分子已换上出射速度 、,也重新视作互不关联:

这里的撇号不是另一个时刻,而是同一次弹性二体散射把入射速度对 映成的出射速度对 ;粒子数、总动量、总动能在这一映射中守恒。入口那一刀,正是 Loschmidt 会死死盯住的地方;出口那一刀,则是在说:一次碰撞制造出的二体关联,会在飞离碰撞区后迅速摊薄,不再影响下一次碰撞的领头阶演化。

最后把 换回 Boltzmann 那个「数分子」的分布 ,质心位置 改回寻常空间点 ,1 号的标签抹掉、2 号的标签改成星号(、),然后注意到对立体角总要积分,在带 的参考系和原系的积分结果总是相同,便得到

一字不差——这正是本文开头那条长长的等式(无外场版本)。左边仍是自由飞行;右边第一项数的是「别的速度对碰完之后撞进 」的增益,第二项数的是「本来速度为 的分子被撞走」的损失。至于 、 与 之间那些看着烦人的 因子,在这一步与截面标度一道,恰好约成一个有限的碰撞率;真正留下来的不是分子编号,而是相对速度乘微分截面。Boltzmann 当年凭物理直觉写下的增益与损失,如今从 Liouville 方程出发,一路积分、截断、代入,堂堂正正地走了出来。

关联在哪里?

到这里,Boltzmann 方程算是从 BBGKY 层级里「推」出来了——但请留意这个「推」字的引号。我们并没有从可逆力学里无损地变出不可逆,而是一路上把关联往外扔了三次:其一,稀薄气体让三体同时挤在一起的机会被 𝒹 压低,于是三分子之间的关联整个不予理会;其二,碰撞时长 碰 远小于平均自由时间 自由,于是碰撞区里那点纠缠被当作瞬时完成、来不及回头影响慢变量;其三,也是最要命的一条——分子混沌假设,把入口与出口的二体分布重新拆回一体分布的乘积,把两分子之间的关联干脆抹平。

三刀砍下去,方程封闭了、不可逆也有了。可被砍掉的东西并没有凭空消失——Liouville 方程从头到尾一分信息都没丢,我们只是不再去看它罢了。那么问题就成了:关联究竟躲到哪里去了?

物理上头一个给出交代的是 Bogoliubov。他指出稀薄气体里藏着三个悬殊的时间尺度——

  • 碰:一次碰撞的时长,两分子真正「贴在一起」的那一瞬,约等于力程 𝒹 除以热速度,短得可怜;
  • 自由:平均自由时间,一个分子两次碰撞之间的等待,约等于平均自由程除以热速度——这恰是 Boltzmann 方程演化的节拍;
  • 流:流体时间,宏观密度、温度梯度松弛所需的漫长时光。

Bogoliubov 的三阶段:横轴是时间,每跨过一道坎,描述就往下塌一级——阶梯的下降本身,就是「描述在坍缩」:自变量的数目一路从 落到 、再落到 。碰 时关联正在生成,非得整条 BBGKY 层级、Γ 空间那 个变量不可;跨过 碰,扔掉高阶关联(分子混沌),只剩单分子相空间的 个变量,也就是 Boltzmann 方程;再跨过 自由,连速度分布的细节也扔掉(局域平衡),只剩密度、流速、温度这 个守恒场——而它们各自只依赖 个空间坐标,便是 Euler 与 Navier–Stokes。本文走的是第一级台阶;第二级不展开,只在最后一节交代一句。

三者一个比一个大得多,碰自由流,每道鸿沟都宽如气体的稀薄程度。关联的身世于是清楚了:它生在碰撞里——历时 碰 的那一瞬,两个分子的速度被牢牢绑在一起;可碰完两人便各奔东西,这点瓜葛被摊薄到一整段自由程上,等它们下一回撞上别人(已是 自由 之后),迎面来的几乎必是个素不相识的新搭档。也就是说,关联在快尺度上生成,又在快尺度上被稀释,根本来不及在慢尺度上攒起来。Bogoliubov 把这幅图像凝成一句「泛函假设」:过了几个 碰,所有高阶分布都不再各自记忆,而是被同一时刻的 奴役、沦为它的泛函;其中最低阶的那个,恰恰就是

下文为了方便,就把这个式子简写为

所以关联并没有被销毁,只是往层级上头跑了。两个分子撞过一次,二体的那点关联很快被稀释;可它并未归零,而是顺着 BBGKY 的链条爬进 、……躲进越来越高阶、越来越细的那些分布里,也躲进「谁曾经跟谁撞过」的那部碰撞史里。我们那三刀,砍的从来不是关联本身,而是不再追问它跑去了哪儿。

而这正是全部风险所在。Bogoliubov 那套说辞终究是一幅物理图像、一个在「稀薄」里取领头阶的启发式论证:它假定了关联一散便不回头,却没有证明。Newton 力学偏偏不禁止两个刚撞过的分子日后再度相遇——一次「重碰撞」,就足以把当年那点旧关联从高阶分布里重新拽回场中,让分子混沌当场失效。而这类反噬一旦不再是零星几例、而是成规模地串通起来,我们在时间之箭那一节设想过的奇观就会当真上演:那 15 粒四散奔逃的台球,一颗颗倒退着重新排回那个规规矩矩的三角形。谁也没见过这一幕,并不是因为力学禁止它,而是因为要把那么多高阶关联凑得恰到好处实在太难——特别是你的粒子不是 15 个,而是 ——可「太难」终究不是证明。短时间里重碰撞稀罕到可以不计;可它们会不会一直稀罕下去?那三刀会不会始终砍得干净?换句话说:一开始几乎无关联的气体,在自身的演化下,能不能一直保持几乎无关联?这个问题此刻已是一道字面精确的数学题,也有了它正式的名字——propagation of chaos(混沌传播)。

混沌会传播吗?

稀薄到刚刚好

到这里为止,「稀薄」一直只是个形容词。我们截断 BBGKY 第二层时用过一个小参数 𝒹,可那句话至多是说「这个数很小」;要把整套推导变成一条能证的定理,就必须交代得更死: 究竟怎样趋于无穷、𝒹 怎样趋于零,以及两者以什么样的步调结伴而行。单说一句 是空话——分子越添越多而个头不变,气体只会越来越稠,最后挤成液体,那与 Boltzmann 方程再无干系。

有两个要求在互相拉扯。

其一,碰撞不能消失。分子若越来越难碰上,方程右边的碰撞项就趋于零,剩下的不过是一群互不理睬的分子各自沿直线飞行——那是自由流动,不是气体的物理,更谈不上什么熵增。

其二,第三者不能出场。我们之所以敢把第二层就地截断,靠的正是「两颗分子正贴在一起时,第三颗恰好也挤进那一小团空间」的机会渺茫。这要求气体足够稀。

一头要碰得足够勤,一头要碰得足够干净——每一次碰撞都只有两个人在场。Grad 看出,这两个要求之间恰好留着一条独木桥,宽度不多不少正好一条线7。

算一算就明白。让一个分子在气体中穿行,它扫出一根碰撞管:管子的截面就是「多近算撞上」的那块横截面,对直径 𝒹 的硬球即 𝒹;管子伸长的速度是典型相对速度 。于是它每单位时间扫过的体积是 𝒹,再乘上数密度 ,便是它每秒撞上别人的次数:

自由𝒹

要让平均自由时间在极限里既不塌成零、又不涨到无穷,只剩一条路可走:

𝒹有限常数

体积固定时,这就是 𝒹 定住不动。这便是 Boltzmann–Grad 标度:分子数按 𝒹 涨上去,直径按 𝒹 落下来,两者的乘积岿然不动,于是碰撞率、平均自由程、平均自由时间统统保持有限——气体还是那团气体。

妙处在下一步。回头看那个稀薄参数:

𝒹𝒹𝒹𝒹

分子实际占据的体积份额自动归零。也就是说,我们当初为截断第二层而手工塞进去的那个小参数,在这个标度下根本不是一条额外假设,而是极限本身的推论——你只要求「碰撞率有限」,稀薄性就自己送上门来。前面留下的因子疑问也在这里解开:第一层方程右边那个 与截面里的 𝒹,单看一个在发散、一个在塌缩,唯有合成 𝒹 才是一个有限的碰撞率。这两个看着碍眼的因子,从来就该被当作一个整体。

独木桥的两侧则都是悬崖。若 𝒹,每个分子每秒要撞上无穷多次,平均自由程缩得比分子间距还短,三体、四体的照面不再稀罕——那是稠密流体的地界,Boltzmann 方程当场失效。若 𝒹,分子彼此擦肩而过,碰撞项归零,又退回自由流动。所谓稀薄到刚刚好,说的正是这条独木桥:稀薄到第三者可以不计,却又没稀薄到无人相撞。

而这条桥并非数学家的臆造。常温常压下的空气,分子直径约 ,数密度约 ,于是平均自由程约 、平均自由时间约 。自由程是分子直径的将近 倍,而稀薄参数 𝒹——小得足以让三体项无足轻重,却又远远不是零。我们每天呼吸的这团东西,本来就坐在那条独木桥上。

标度既已钉死,「从 Newton 方程推出 Boltzmann 方程」这句话,才第一次有了可以下笔去写的形状。

混沌的初始

万事俱备,只欠东风。剩下的只有一件事:在 Boltzmann–Grad 极限下,把上文那个物理问题翻译成一句数学上说得清、也驳得倒的话。

怎么翻译?统计力学的骨子里就是概率论,那自然该用概率论自己的语言——Kolmogorov 在 1933 年立下的那一套:把概率架在测度之上。这套语言格外对本文的胃口,因为它天生就分得清「几乎所有」与「例外」,而我们一路上所有的纠结,归根结底都卡在这个分寸上。

翻译只需三步:说清楚样本空间是什么、随机变量是什么、收敛又是什么意思。这一节先看样本空间。

既然要谈概率,那微观那一侧交出来的,就不再是某一个具体的微观态——不说「1 号分子此刻在这儿、速度是那么多;21 号分子此刻在那儿、速度又是多少」——而是整个 Γ 空间,连同其上的一个概率测度:初始条件落在何处,各有多大的可能。请留意这里的一个样本点是什么:它不是一个分子,而是整团气体的一整套初始条件—— 个数,全体分子的位置与速度一次性给齐8。抽一次样,抽到的就是这么一整套。

而 在概率论的话语体系下,就被描述为测度的密度:对 Γ 空间里任意一块微观态 ,对这个集合积分就得到它的概率:

诸君即便不熟悉测度论也不必紧张:这里其实没有什么新东西,在 Liouville 那里我们已经跟它打过交道,「概率测度」说人话就是概率,这一整套无非是虚张声势的统计物理。唯一真正新添的本事只有一样:测度是给「集合」称重的,而不只是给「点」赋值。后文那句「例外稀有到可以忽略」,正是靠这杆秤才说得出口。

然而在 维空间上定义测度的方式有千千万万,数学在这里什么话也没有。Γ 空间上任何一个概率密度都是合法的:你可以让它集中成一个尖峰,那等于干脆指定了某一个具体的微观态;也可以让它专挑那些「排布精巧、专等着熵减」的初始条件去偏爱。Newton 力学与 Kolmogorov 的公理,对此都不会有半句异议。换句话说,Loschmidt 那个反演态、15 粒台球复位的那个态,都是合法的微观态,Newton 力学一条也没有禁止;你若要求结论对每一个初始条件都成立,H 定理当场就被它们推翻。为了得到热力学,我们必须在测度上对他们加以限制。

所以接下来这一步不是推导,是约定:本文自始至终只考虑「分子之间事先毫无串通」的那一类测度——写成式子,就是单分子分布的直积。

每个分子按同一个 独立撒下去——别被「同一个」误导: 是任意的单分子分布,可以是热平衡的 Maxwell 分布,也可以是「一大团热浴外加一丛确定能量的粒子」这样的双峰9;「独立同分布」约束的是分子之间不串通,而不是要它们处境相同。

不过刚才那个例子里藏着一处值得当心的地方。「往热浴里打一丛束流」若照字面写成分布,其实是

束浴束束

——它把「哪些编号属于束流」钉死了,因而对编号交换并不对称,与上面那个式子严格说来并非同一个测度。补救的办法是对称化:对「究竟哪些编号属于束流」求平均。这样得到的测度恢复了对编号的对称性,而它的固定阶边缘,恰好就是那个双峰分布的直积。

可这道手续凭什么是合法的?凭什么把编号搅匀之后,物理分毫不改?这要等下一节把「我们究竟在观测什么」讲清楚了才好回答——先记下,姑且称作悬念一:它的答案,在于那个观测量里。

至于那个「」,它的来历值得说清,因为含糊不得:硬球终归不能重叠,而这件事根本不是加在测度上的修正,是样本空间本身就少了一块——球心靠得比 𝒹 更近的那些点压根不在 Γ 空间里,谈不上概率。于是上式精确的意思是「把 个独立同分布的分子撒下去,扔掉重叠的那些样本,再重新归一」,那个「」不多不少就等于这道手续。好在它小得无害:被扔掉的份额正是上一节那个自动归零的 𝒹。

这个出发点有个名字,叫初始混沌。这里务必看清一件事:这条约定并没有禁止任何微观状态。15 粒台球倒退着排回三角形的那个初态,在这个独立同分布的测度下照样分到一份不为零的概率密度——你真要去抽,是抽得到的;它既没有被剔出样本空间,也没有被赋值为零。我们做的全部事情,自始至终只有一样:通过测度让它变得「很不可能」。不靠禁止,只靠称重——这也正是前面说「测度是给集合称重的」那句话的用武之处。至于为何如此一来这些态便「很不可能」,此处暂且按下不表,留作悬念二,待到下节,自然清楚。

不过这里还有件事值得说道说道:这条约定凭什么立得住?

先把最难听的一句说出来:这确实是一条约定,不是推导。Newton 力学里没有任何一条方程指定初始测度;定理是一句蕴含式,而它就是前件。所谓「从 Newton 方程推出 Boltzmann 方程」,严格说永远拖着一个前提从句——这一点不必遮掩。

但约定也有高下之分,而这一条并不是随手挑出来的窄类。「每个分子独立地按 初 撒下」这句话,正是「我只控制了宏观廓线,此外什么也没安排」的数学写照——而这恰恰就是实验室里「制备一团气体」的全部含义。反过来,要制备一个带关联的系综,你得逐对分子去安排它们的相对位置与速度;那不是拧个旋钮能办到的事。

「可是,能写成直积的分布,在全体分布里只占极小的一部分。」这话对。物理学家对此再熟悉不过:量子力学里,不纠缠的可分离态在全体量子态中同样是非典型的稀客10。

那凭什么可以这么设?答案我们其实早就写下过,只是当时用的是另一套词——宇宙学箭头。精巧的关联和低熵一样,是一种不可再生资源;宇宙给的初始存货里,根本没有那么多关联供你去搭一个反常的初态。「初始混沌」这条假设,就是那段话的数学编码。于是整条链子清清楚楚:可逆的 Newton 力学,加上一个来自宇宙初条件的前件,才生出不可逆的 Boltzmann 方程——箭头是随着前件一道进来的,一步也没有凭空产生。

「那演化之后呢?总不至于一直待在那个稀有的集合里吧——那岂不是更离谱的精细调节?」问得好,它一秒也没待过。

事实上连 都不在上面。还记得那条排除条件么——若真按独立同分布把分子撒进去,重叠的对数期望是 𝒹;代进常温常压的 空气,这个数是 对。任意指定的一对分子重叠的概率不过 量级,可全体加起来却是个天文数字。所以真实的测度与那个理想的直积测度,从一开始就相去甚远,谈不上「待在里面」;碰撞一旦发生,更是把它推得越来越远。

于是问题变得尖锐:既然那个 维的测度一秒也没有停在直积的形状上,我们说「初始的混沌被传播下去」,到底还剩下什么意思?

从定律到定理

要明白这点,我们得看看翻译还剩的那两步:说清楚随机变量是什么,以及收敛是什么意思。

随机变量。我们拿什么去比对微观与动理学两幅图像?答案朴素得很——划格子数数。在单分子相空间(三维位置加三维速度)里圈定一块区域 ,去数有多少分子落在里面、占了总数的几分之几:

⟦时刻落在中的分子⟧

而要证明的,无非是说:在 Boltzmann–Grad 极限下,这个随机变量 依概率收敛到 Boltzmann 方程的解(经过归一化后) 在 上的积分。

读到这里多半有人要抗议: 本来不就定义成「某处、某速度附近有多少分子」吗?那这个 岂不就是 本人,说它收敛到 ,岂不是同义反复?这个抗议正中要害。而答案是:不是同义反复——因为「有多少分子」这句话,藏着两种截然不同的读法。

第一种,对一次具体的微观实现去数。抽定一组满足初始混沌的初始条件,让它严格按 Newton 方程演化, 时刻数一数格子里有几个球。这样得到的 有三个要命的特征:它是随机的——换一组初始条件就换一个值;它是尖刺状的——不过 个孤零零的点,谈不上什么光滑的密度;而且它不满足任何偏微分方程,它只是 条 Newton 轨道的一个副产品。在这个基础上, 是一个结结实实的随机变量,而且只能取整数。

第二种,把它当作一个光滑的函数,并且要求它自己满足一条封闭的演化方程。那才是 Boltzmann 的 ,你即便把他对区域 积分,也只能得到一个平均值,而不是某一次具体的读数。

很明显,这两者在有限的 上根本不是一回事。更要命的是:即便你把第一种按概率测度去积分(其实就是取平均)、把随机性抹掉,得到的 虽然光滑了,却依然不满足 Boltzmann 方程——它满足的是 BBGKY 的第一层,而那一层不封闭,右边挂着 ,正是我们费尽周折才砍断的那根链条。

所以定理要断言的,是两件都不平凡的事捆在一起:

  • 涨落消失:单次实现数出来的 ,向它的系综平均集中;
  • 平均满足封闭方程:而那个系综平均,在极限下由一条自足的六维方程支配。

Boltzmann 当年一步就跨过了这两件事,把「数出来的分布」与「满足我这条方程的分布」直接当成同一个东西。要证的,正是这个等同——它只在极限下才成立。

观测量既已摆明,上一节挂着的悬念一就可以还了。而说清楚它,只需给「重排编号」这件事一个算符:记 为把全体分子的标签洗一次牌的操作(交换两个、三个或者任意多个,下标应该不同,但结论一样,此处省略)。于是要点化成两条「对易关系」。

其一,观测量不看编号。 数的是「有几个分子落在 里」;把标签洗过牌,落进 的还是那几个相点,读数一个不差。
其二,动力学也不看编号。硬球彼此全同,哈密顿量对重排编号不变,于是幺正演化算符 与它对易11:

说人话就是:先洗牌再演化,与先演化再洗牌,给出同一个位形。

两条一摆,剩下的就是物理学家最熟的那套逻辑了——与演化算符对易的算符,其本征值是守恒量。初始测度若是 的本征态(本征值 ,也就是对编号交换对称),那它永远是;于是一个测度与它对称化之后的版本,给出的 的分布律逐时刻完全相同。那道对称化手续确实免费,而且不是「物理上差不多」,是字面意义上的一模一样。

归一化上还有一处要当心。 是个占比,归一到 ;所以它只能收敛到同样归一到 的那个分布,也就是 ,而不是 Boltzmann 方程里那个「数分子」的 ——后者随 一道发散,本身谈不上收敛。两者只差一个常数因子,而这个因子连同截面一起,恰是前面钉住的 𝒹:正因为那个乘积保持有限,除掉 之后剩下的方程才不多不少、仍旧是 Boltzmann 方程。

收敛。现在,要证的话可以说出口了。在 Boltzmann–Grad 极限下,

依概率

而右边那个 ,正是 Boltzmann 方程以 初 为初值的解。

这是一条大数定律。就像抛 次硬币、正面的比例依概率收敛到 那样:分子数趋于无穷时,那个随机的格子读数集中到一个确定的值;而支配这个确定值随时间如何演化的,正是 Boltzmann 方程。

于是「Boltzmann 方程是 Newton 力学的一个近似」这句流传已久的话,其实说轻了,也说偏了。它不是近似——它是极限:不是某个误差还可以再估一估的替代品,而是随机涨落被 碾平之后,剩下的那条确定性法则。热力学第二定律与掷骰子的大数定律,在这里原是同一件事的两种说法。

大数定律只管一体的读数。要回答上一节结尾那个尖锐的问题,还得看更高阶的分布。更强的命题是:对每一个固定的阶数 ,

这就是上一节那个问题的答案。真正被动力学一路带下去的,不是那个 维测度的直积形状——它一秒也没有停在上面——而是每一个固定阶的边缘的可分解性。态早已飞离那张薄片,可它投在任何低维屏幕上的影子,看上去仍旧是乘积的。初始的无关联不是被演化摧毁,而是以这种方式被它传播下去——propagation of chaos 这个名字,字面讲的就是这件事。

这里的量词顺序不可颠倒:先固定 ,再让 。任意有限个分子之间渐近独立,绝不等于全体 个分子彼此独立。关联并没有消失,它只是被摊薄到越来越高的阶上去了。

顺带厘清一个用词。前文一路说「某个态编码着复杂的关联」「凑够高阶关联」,可概率论里的「关联」严格讲是测度的属性——联合分布不等于边缘之积。单独一个微观态不过是 个数,既谈不上相关不相关,也给不出什么光滑的密度:它在六维空间给出的只有 个尖刺。那么「这个态编码着关联」,到底在说什么?

答案是:判断有没有关联,根本不需要光滑密度,数数就够了。

单体的数法,我们已经在用。挑一块区域 ,数落在里面的分子、除以总数,得到的就是 ——「有多大比例的分子落在 里」。

二体的数法,把同一个动作照做一遍,只是改数成对的分子。挑两块区域 、,去数有序分子对(两个编号不同)里「前一个落在 、后一个落在 」的有多少,再除以有序对的总数 :

⟦时刻一个落在、另一个落在的有序分子对⟧

而关联,就是这两个数对不上。分子之间若毫无串通,「一个在 、一个在 」的对所占的比例,理应等于两个占比相乘:

对不上,就说这个态带着二体关联;三体照办,数三元组、跟三个占比之积相比。从头到尾没有出现任何光滑函数——等号两边都是从那 个数直接数出来的两个比例12。

眼尖的人一定已经发现,这与之前所说的关联,严丝合缝地对应:取系综平均, 给出的正是 , 给出的正是 。经验的那一份是单次实现数出来的,系综的那一份是它的平均。两重意思于是接上了——一个测度「混沌」,正意味着从它抽出的典型微观态,这些数对得上。

既然如此,上文挂着的悬念二的答案也就呼之欲出了。

那些反常的初态,恰恰就是「这两个数对不上」的态。15 粒台球要倒退着排回三角形,意味着此刻成千上万对分子的相对位置与速度都被调得恰到好处;用刚才的话说,就是 在「即将迎面相撞」那一类区域上,远远偏离 。

而独立同分布的抽样,恰恰就是让这两个数高度吻合的抽样方式。这不是额外加的什么规定,这就是「独立」二字的全部内容:一个一个、互不商量地撒下去,数出来的对子比例自然趋近于占比之积。要抽出一个显著偏离的样本,需要的正是这些分子集体串通——而串通得越广,概率掉得越狠。

所以「很不可能」的来历,说穿了平淡无奇:我们选的那把秤,本身就是按「互不串通」刻度的;凡是非靠串通才能成立的态,在这把秤下自然轻得可以忽略。它一点也没有被禁止——只是称出来的分量太轻了罢了。用 Kolmogorov 的语言:依概率收敛,说的正是这个集合的概率测度趋于零。

要留意是「趋于零」而非「等于零」:台球复位的那段倒放画面,逐帧都严格满足力学;它之所以不会出现在你眼前,不是因为力学不许,而是因为你抽不中。「不可逆」自始至终是一句关于概率的话——而现在,它是一句可以证明的话。

不过要当心,「看不见」不等于「不存在」。关联并没有蒸发,它只是躲进了更高的阶里。你去算二体,对得上;算三体,对得上;一路算到十体、百体,还是对得上——可谁也没担保它不藏在第一千阶、第一万阶上。定理许诺的是「每一个固定的阶」,而从来不是「所有的阶」;那道缝隙,正是关联的藏身之处。

这幅图景物理学家同样不陌生:一个热化了的多体量子态,少数几个格点的约化密度矩阵干净得像乘积,全局却纠缠得一塌糊涂。纠缠是真实且巨大的,只是少体观测量探不到它。

于是「那些关联不重要」这句话,也终于有了准确的意思——它不是说「离直积态不远」。全局意义上,真实的测度离直积远得很,关联的总量非但不小,反而大得惊人;它只是分散在多得数不清的高阶自由度里,被任何固定的低阶窗口漏了过去。

翻译至此完成。最后,把定理没有声称的东西也摆明。它说的是硬球,彼此只在接触的一瞬发生作用;是稀薄气体,坐在那条独木桥上;是经典力学,没有量子。它给的是一个极限,不是有限 下的误差估计。它更没有说 Boltzmann 方程对稠密气体、对长程力、对任意长的时间都成立。

于是整句口号,终于有了它完整的形状:在 Boltzmann–Grad 极限下,从几乎无关联的初始测度出发,在 𝓉 这段时间里,硬球气体的格子读数依概率收敛到 Boltzmann 方程的解。

只剩一个空格没有填:𝓉 能有多大。

这个空格分量极重。若 𝓉 短得只够分子碰上一两次,定理虽然干净,却离「气体的宏观行为」还差得远——熵增、弛豫、趋向 Maxwell 分布,哪一样都要熬过无数次碰撞才谈得上。而 𝓉 还受着另一重牵制,来自一个也许出人意料的方向:Boltzmann 方程自己的解存不存在、光不光滑,本身就是一桩没解决的难题。你没法证明微观的读数收敛到一个你尚且不确定是否存在的东西。所以这类定理的形状注定是有条件的——「在 Boltzmann 方程拥有一个足够光滑的解的那段时间上,微观的读数收敛到它」。方程解的寿命有多长,严格化的战线就只能推到多远。

展开一棵碰撞树

𝓉 那个空格,头一个把数字填进去的是 Oscar Lanford——1975 年,在一次讲习班的讲义里。他填进去的那个数很难看,可那是人类第一次真的把它填上。

招数说来平淡。回想 BBGKY 那条链子的窘境:第一层告诉我们 由 决定,第二层说 又由 决定……我们当初的对策是斩断它——在第二层就地截断,硬把 拆成乘积。Lanford 换了个法子:不斩,往回代——直接写下 Liouville 方程的一般解13:

「怎会如此!」物理学家这时候一定脱口而出——这不就是 Dyson 级数吗!

正是如此。要弄明白怎么回事,请暂时把目光从经典力学移开——我们要看的是量子力学里的一位老朋友:编时乘积算符。含时微扰论里,要算一个受扰体系怎么演化,办法是把所有的相互作用按发生的时刻排好队:

读法诸君再熟悉不过:第 项就是「一共受了 次微扰」,那 个时刻排好队、各自积分——Feynman 图数的正是这些项。指数展开本该带的 ,被「只取排好序的那一块区域」吸收掉了。

现在把整套东西搬到碰撞上来。要展开的不再是波函数,而是 ;扮演「微扰」的不再是势,而是「外人闯进来撞一下」。于是需要两个算符: 是上一节那个幺正演化算符,只不过如今只管 个硬球,把它们按精确动力学演化一段时间; 则「添一个新分子,让它撞上来」。

后面这个算符不是单次碰撞,而是一整包。真实世界里还有 位外人不在碰撞族谱上,每一位都可能是闯进来的那个,严格写下来要对他们逐个求和;而每一项里挂着的,是把「族谱那 位 + 这个特定外人」之外统统积掉所得的边缘分布。好在这些项根本就是同一个函数,而理由上一节已经备好了。那里我们请出过一个算符 ,它把编号洗一次牌,并且与演化算符对易;于是「对编号交换对称」是个守恒的本征值——初始测度是 的本征态,此后一直是。既然如此,那位外人究竟原名 还是第 号,积出来分毫不差。再加上碰撞算符是线性的14,求和号便整个提到前面,只留下一个计数因子:

可是且慢——「撞的是谁」那个求和,凭什么不能照样塌掉?族谱里那 位,在 里明明也是彼此对称的。分野只有一条:谁被积掉了。把第二层的碰撞项摊开写最清楚——族谱里已有 1、2 两位,外人叫 3:

碰

外人的坐标 、 是被积掉的哑变量。积分一做完,它在结果里连一丝痕迹都不剩——那位外人原名 3 号、17 号还是第 号,给出的是同一个函数。于是那 项合并成一个计数因子: 的来历就在这里。而族内的 、 是结果的自由宗量。方括号里那两项,一项踢 1 号的速度,一项踢 2 号的速度——是两个不同的函数。结果记得踢的是谁,而这个信息没有被任何积分抹掉。

有了这个算符,BBGKY 的第 层就能写成一条 Duhamel 公式。等号左边写完整的变量表达式;右边的时间则作为函数族的外部参数,使演化算符与碰撞算符只作用于相空间函数:

读法很朴素:若没有外人闯入,这 个分子就从 到 精确演化;若有外人闯入,就追溯到最后一次有外人闯入的时刻 ,先让族谱里原有的 位从 精确演化到 ,再在 插入一次 ,把问题交给一个 粒子的分布。可这个 粒子的分布也未知——同一招再来一遍便是。如此一路往回代,每代一次就多出一个更早的碰撞时刻,也多出一位新加入的分子。

本节开头那条总公式的第一个等号把零次、一次、二次碰撞逐项摊开,第二个等号再按「一共插了 次碰撞算符」收拢。一模一样的形状:演化一段、插一次、再演化一段、再插一次……时序积分也一样。第 项数的是这个分子一共被外人撞了几次;而那串编号 数的不是真实气体里的身份证号——那些早被 抹平了——而是「第几个被纳入族谱」。

不过还有三处得说清楚。

其一,这里的「自由段」并不自由: 装着这一组分子彼此之间的全部碰撞——他们撞多少次都在这个算符肚子里,一次也不往外冒;真正被当作「微扰」剥出来的,只有组外闯入这一件事。所以这个展开不是按相互作用的强度排序,而是按牵扯进来的分子数排序。这句话我们其实早就说过,只是当时是另一副面孔——BBGKY 那节的结论:任何 人小组都按组内力学精确演化,唯一的修正来自组外闯入的那一个分子。当初是层级的结构,如今成了展开的规则。

其二,每插一次,体系就大一圈。量子那边插入 不改变体系;这里每插一次就添一位新分子,所以那串下标 一路往上涨。

其三,最右端落在 。那个 ,正是我们唯一确切知道测度的地方;按初始混沌,它近似于 。级数里的演化算符实际作用时,再把这个裸分布量展开成前式所写的相空间函数。

于是整件事有了清清楚楚的形状——把「 时刻怎么样」这个未知,换算成「所有可能的来路各自贡献多少」,而每条来路的贡献,都由已知的初条件决定。到这一步为止,一丝近似都没有做;这仍是 Liouville 方程的原话,只不过换了一种组织方式。

真正的关键在第二条上头——正是它会长出一棵树。那个塌不掉的求和 正是分岔的来处:每往回退一步,都要在「撞的是谁」上分一次岔。量子那边一阶一阶展开出来的是一条链,这里却是一棵树——树根是我们盯着的那个分子,在 时刻;每退到一个 就分出一条新枝(那一位从此加入);枝上还能再分枝;一直分到 ,枝梢落地。一条碰撞史,就是把每一步的选择都定死:第 次退回时新来者撞的是谁,定出树的形状;那次碰撞发生在哪个 、以什么样的瞄准参数与方位角,定出树的尺寸。形状离散、要求和,尺寸连续、要积分,上式两样都写着。一棵碰撞树记录的,恰恰就是那些没被积掉、因而留下了痕迹的选择——「哪个外人」既已被积分抹平,它连一根枝都不占。而链变成树这个看似无害的差别,正是后面全部困难的源头。


倒着展开出来的碰撞史。每一次碰撞都是一个节点:树上的节点是三叉的——一条线从上方进来,两条线往下走(原来那位继续、新来那位加入)。左:一棵树,每退到一个 就撞开一个节点,四条边落到 ,初始混沌正是在那里求值。边上的编号就是「撞的是谁」那个求和选出来的东西: 那次撞的是 2 号,而不是 1 号或 3 号。中:同一棵树,外加 处 1 号与 2 号再撞一次——两条边绕成的环。留意这个环非绕不可:1 号在 被 3 号打偏、2 号在 被 4 号打偏,它们才可能重逢;纯自由飞行下两球一旦分开就再也回不来,所以环上必然挂着别人。右:换一种,2 号与 3 号从没直接撞过,只是都当过 1 号的对手,却在 撞上了——三条边绕成的环,没有哪一对撞过两次。两种都算重碰撞,也由同一步一并排除。

这个展开要用得成,还得先过两关。第一关是收敛:树的数目随着退回次数疯长,每一棵都要算、都要加起来,加得起来吗?第二关是环——而这一关的对手,我们早就见过。

什么叫「环」?树规定了哪几次碰撞;可它没规定两条分叉不许再连起来。一旦有哪一次碰撞把两条已经分开的枝重新连上,那就不再是树,而是带环的网。(这里的「已经分开」是顺着展开的方向说的——展开倒着走,所以在物理时间里,那是「将要分开」。)这类碰撞统称重碰撞15:它藏在 那段组内精确演化里,不会长成新的枝。

「重」字容易把人引偏,以为非得是同一对分子撞两次。环确实有这一种:1 号与 2 号在 撞过,后来又撞上了,两条边绕成一个环。可它远不是唯一的一种——譬如 2 号与 3 号从没直接撞过,只是都当过 1 号的对手;若这两位在半路上撞了一下,三条边便绕成一个环。判据是闭合了环,不是重逢。

带环恰恰就是分子混沌失效的形状,而它失效在什么地方,两种环讲的是同一件事。回想那把刀砍在哪儿:我们假定即将相撞的两个分子彼此素不相识——说得准确些,是假定这两位此前的履历涉及完全不同的分子。第一种环里,两位的履历直接重叠:上回就是彼此的对手。第二种环里,履历间接重叠:都撞过 1 号,而 1 号把 2 号的消息捎给了 3 号,这条链客观存在。哪一种都让「素不相识」当场作废。树 ↔ 履历互不相交,环 ↔ 关联反噬——两幅图像在这里严丝合缝地对上了。

Lanford 证明的第一件大事,是环可以忽略:在 Boltzmann–Grad 极限下,带环的那些碰撞史,总测度趋于零。道理不难体会,而且对两种环是同一个道理——一个尺度估计就能看明白。

先看树上的碰撞怎么算。新加入的那位可以是 位外人中的任何一个,而每一位每单位时间扫出的碰撞管体积是 𝒹;两者一乘,正是前面钉死的那个有限碰撞率

𝒹

Boltzmann–Grad 标度保的就是它:𝒹 在塌缩, 在发散,合起来岿然不动。

可闭环那一次没有这个 。撞的是谁、被撞的是谁,图上早已点了名,「可以是任何一位」这份自由度用不上;剩下的只有孤零零的一根碰撞管,于是

环𝒹

一句话:树上的碰撞有 位候选人替那块小截面撑腰,闭环的那一次孤立无援。

几何上,这件事长成这样:其余的边一旦定下,图里各条轨迹就都被钉死了;此时要再闭上一个环,等于要求某一条尚且自由的轨迹恰好瞄准另一条已经定死的轨迹——速度与偏折角得凑得恰到好处,差之毫厘便擦肩而过16。至于这个环是两条边绕成还是三条边绕成、那两位从前认不认识,分类上并不重要:要付的代价来自闭环本身。

于是每一段 都退化成自由飞行,组内不再有额外的碰撞,只剩树上规定的那些,级数里也就只剩下树。

而剩下的这些树,逐项恰好凑成 Boltzmann 方程的迭代解。回头看一眼 Boltzmann 在 1872 年写下的方程,右边那个东西我们叫它碰撞算符,它吃掉两个分布函数、吐出来一个:

请留意它长在 上,不是 ——Boltzmann 方程本来就是「数分子」的方程。把碰撞项挪去当源项,左边只留自由飞行,就得到它的 Duhamel 形式:

读法与层级那边一模一样:第一项管「一次也不碰」,纯是自由飞行;积分那一项则把最后一次碰撞发生的时刻 穷尽一遍,于是右边冒出那一刻的 ——把整条式子代回它自己便是,这就是 Picard 迭代。碰撞算符的两个宗量,正是分岔的来处:每代一次,一个 裂成两个;再代一次,四个。于是迭代出来的东西自动是一棵二叉树,第 项有 个排好序的碰撞时刻,末端挂着 个初始分布。而这棵二叉树,想必诸君都已经猜到了——这不就是碰撞树吗!

于是逐项对起来即可。环既已砍尽,碰撞树剩下的清一色是二叉,与 Picard 迭代长出来的那些棵棵对得上。剩下三处核对都不费事:自由飞行那一段,两边本是同一个算符;末端的初态,就是那个直积;至于碰撞算符,把 Boltzmann 的换成 Liouville 的,中间只有一点琐碎的归一化要走17。这几样核完——其实不必真去核,摆在那儿已经看得出来——两套语言逐项对上,一棵树对一棵树,一个杈对一个杈。不是「近似地像」,是同一组项逐项取了极限。

而最值得回味的是这一刀砍在了哪儿。当初推 Boltzmann 方程,我们是硬拆:在第二层就地截断,把 按手换成 ——那是一刀砍在方程上。如今这一刀不必砍了:什么也没有拆,只是把一撮总测度趋于零的碰撞史扔掉,剩下的东西自己长成了 Boltzmann 方程。Boltzmann 当年凭直觉写下的那个碰撞项,就这么从一堆倒着生长的家谱树里,一项一项地长了出来。

至此,propagation of chaos 头一回不再是物理学家的信念,而是一条被证明了的命题;可逆的微观力学与不可逆的 Boltzmann 方程之间,不存在矛盾——1872 年写下的那句话,1975 年终于画上句号18。

但坏消息在第一关:那个级数,只在很短的时间里收敛。而这个增长,正好落在「链变成树」那个差别上。先数一数树有多少棵。诀窍只有一句:场上有几位,新来的就有几种撞法。起初场上只有树根那一位,所以第一个加入的分子别无选择,只能撞他—— 种。这下场上是两位了,第二个加入的可以挑其中任何一位—— 种。场上又多一位,第三个便有 种。轮到第 个加入时,场上恰好已有 位,于是 种。一棵树的形状,就是把这一路上每一步的选择全部定死;把各步的选法数乘起来,

——退 次,形状共有 棵。

留意这 是从哪儿冒出来的:它正是那个塌不掉的求和累乘的结果。外人的编号被积分抹平了,「谁闯进来」那个求和于是塌成一个因子;族内的编号却留在宗量里,「撞的是谁」那个求和塌不掉,第 步有 项。族谱之内编号算数,族谱之外不算数——同一条分野,先送来 凑成有限的碰撞率,如今又送来 把收敛半径压死。这正是量子那边没有的因子:那里每一阶只有一条链,这里有 棵树。

救兵还是编时那块区域。 个碰撞时刻必须排好队,,而这块区域的体积正是

——就是刚才被指数展开吸收掉的那个 ,此刻它长在分母上。两边一乘恰好抵消,只剩 。而每退一次还得为那次碰撞付一份代价,正是刚才那个碰撞率 。于是级数第 项的大小约为 ——一个几何级数, 时收敛,一旦越过 就发散。而 本就是碰撞率,也就是 自由,于是收敛的条件赤裸裸地摆在那里:

自由

Lanford 能证到的时间,是平均自由时间的一个零头。代进常温常压的空气:平均自由时间约 ,定理管辖的还不到这个数。也就是说,它覆盖的是一个分子还没来得及撞上第二下的那段光阴。而熵增、弛豫、趋向 Maxwell 分布——哪一样都要熬过无数次碰撞才谈得上。我们费尽周折想严格化的那个「气体的宏观行为」,在 𝓉 的这个取值下,一样都还没有开始。

面对这样一个结果,第一反应总是「证明手法太笨,换个聪明点的就好了」。可事情恐怕没这么便宜——在那个时刻附近,被展开的对象本身变了性质。想想那棵树在长什么:时间越长,牵扯进来的同伙越多,树也就越大。数值与启发式的论证提示,当 接近一个平均自由时间,碰撞图的尺寸会突然爆炸——从一个个互不相干的小团簇,骤然长出一个囊括几乎所有分子的巨分量19。这个转变物理学家眼熟得很:它就是渗流,与随机图越过临界点时长出巨分量是同一回事。而一旦所有分子被同一张网连在一起,「把碰撞史一条条列出来分别求和」这个想法本身就垮掉了——不是加不动,是根本没有「一条条」可言。

所以那个短得可怜的 𝓉,并不单是证明技术的缺口。Lanford 之后,这条路上又走了三十年,把论证的每一处细节补齐、把结论推广、把估计做得更精细20——可没有人能把 𝓉 推过那道坎。

五十年弹指一挥间,少年子弟江湖老,但这道坎还是卡在那里,不知多少英雄折腰。

二次元也可以拿菲尔茨奖!

先把那道坎重新量一量。上一节我们算出 Lanford 的 𝓉 是平均自由时间的一个零头,可有一句话没说:这个数,恰好也是「我们会解 Boltzmann 方程多久」的那个数。

这不是巧合,上一节其实已经把证据摆在眼前了。微观那条 Dyson 级数,第 项的大小是 ;宏观那边解 Boltzmann 方程用的 Picard 迭代,第 项的大小也是 。两条级数逐项对齐,那么收敛半径当然也得一模一样。换句话说:Lanford 的证明之所以在一个平均自由时间处失效,与「不动点法解 Boltzmann 方程只能解这么久」是同一件事的两个说法。于是那道坎的性质变了:五十年来卡住微观推导的,其实是我们解宏观方程的本事。这不是气体的性质,是方法的性质。

邓煜–Hani–马做的事,就是把这根绳子剪断。他们的定理不是「把 𝓉 往前推了一截」,而是换了一种句式:

你给我一个光滑的 Boltzmann 解,它活多久,我就把微观极限推多久——不管这个解是怎么造出来的:能量方法、稳定性理论,还是别人二十年后才证出来的存在定理。

写下来是这样。这里的 𝓉 还是上一节那个角色——「证明管得到多远」——只是换了位主人:Lanford 的 𝓉 由他那条级数的收敛半径定死,邓煜–Hani–马的 𝓉 则由 Boltzmann 解的寿命定。他们证明的核心是下面这条定理:在 Boltzmann–Grad 低密度标度下,对时间区间 𝓉 里的每个 、以及每个阶数 𝒹,

𝒹固定(在意义下)

唯一的条件是:Boltzmann 方程的解 在整个 𝓉 上连同它的空间梯度一致有界21。

请留意这不止是 。式子说的是:随便挑 个分子,它们的联合分布收敛到单分子分布的直积——这就是混沌传播(propagation of chaos) 本身,是本节标题那个问句的正面回答。混沌不但传播,而且一直传到 Boltzmann 解存在的最后一刻。与 Lanford 相比,有三处不同。头一处是时间:从「平均自由时间的零头」变成「Boltzmann 解存在的整段」。第二处是代价:收敛只在 里成立,比 Lanford 那边弱一些。而第三处才是全部意义所在——限制的来源变了。Lanford 的 𝓉 说的是「我们的手艺只到这儿」,邓煜–Hani–马的 𝓉 说的是「方程自己只到这儿」。前者是我们的短处,后者是问题的形状。而只要 Boltzmann 解全局光滑、一直一致有界,𝓉 甚至可以随 𝒹 一起发散——现有证明给到 𝒹 的量级,作者明说不认为这是最优的,期望能推到 𝒹。而把这个速度定量化,正是通往流体极限的那扇门。

那么,怎么做到的?说穿了只有一句话:把 Lanford 那一段首尾相接地用 次,并在每个接口处重新制备一个直积态。请留意:他们没有把 Lanford 的射程推远一分一毫。射程还是那么短,短得可怜。变的只是反复用——既然级数只在一个平均自由时间之内收敛,那就把 𝓉 切成 段,每段都短过那个收敛半径,一段一段地走。段长记作 层。任谁都会想到这一招,而且它在每一段之内确实管用——层内就是一次完整的 Lanford,上一节那套原封不动地跑得通。

出事的是接口。Lanford 那套展开的最右端落在 ,靠的是那里我们确切知道测度:初始混沌,一个干干净净的直积。可第二层的「初态」是第一层末尾的状态,而那时候分子已经互相认识了——撞过的两位之间留着关联,没直接撞过、却有过共同对手的,也留着。请务必回顾前文的强调——关联只是躲到高阶了,并不是演化后依然离直积态不远。那直积态没有了,地基就没有了。于是接口处要做一件事:把地基重新铺上。做法是个三步循环:

  1. 读出——把这一层末尾的 体态积掉其中 个分子,留下一个单粒子边缘分布,记作 主;
  2. 重制备——拿这个 主 去张量幂,制备一个新的直积态。地基就这么补回来了;
  3. 再跑——让这群「重新素不相识」的分子,按真实的动力学跑完下一层。

跑到下一个层末,再从第 1 步来一遍。请留意每一层末态是循环自己跑出来的,不是真实气体的末态——它的初态本就是上一轮制备出来的那个假直积。把各层串起来,主 就是一条离散的链——只在 层层𝓉 这 个接口上取值,一步一个格点,每一步就是上面那三步22。

到这儿有个岔口,务必看清:从此以后有两条链在跑。

  • 真实那条:气体老老实实按 Newton 方程演化,谁也没碰它。它每时每刻都有定义,处处连续。
  • 辅助链:每到接口就挨一刀, 体态被一把换成直积。它在接口上是跳的——而这个跳变不是物理,是我们亲手做的手术。

主 是辅助那条链在格点上的读数。层内它不需要有定义;所有接口数据都只在格点上更新。至于每一刀切下去丢掉的东西——那正是下面要追踪的误差。

顺带把这个新面孔认准,免得会错意。它是个单粒子函数,一个 、一个 、一个 ,就这么多;第 个接口上的「直积态」是 主层,不是它自己。 时真实状态本就是直积,所以 主,种子是精确的;从第一个接口起两者才分家——因为第 2 步制备出来的直积是假的,真实末态并不是直积。至于它跟 Boltzmann 解的关系:也不相等,它比 Boltzmann 解更贴近真相。三者由外而内是:Boltzmann 解、主、真实的 。

这套动作物理上是什么,值得单说一句:每过 层,就把这罐气体重新制备一次。不是换一罐气体——密度、速度分布、动量、能量统统原样保留;被抹掉的只有「谁跟谁有过交情」这一件事。而它的名字我们前面早就见过:这就是 Bogoliubov 的泛函假设。他说,过了几个 碰,所有高阶分布都不再各自记忆,而是被同一时刻的 奴役、沦为它的泛函。三步循环正是把这句话当真执行一遍:系统的记忆在每个接口上被清零,只留下单分子分布这一个状态变量。换成今天的说法,这一步就是马尔可夫化——把一个有记忆的、可逆的 Hamilton 系统,硬按成一个无记忆的过程。上面那条闭合的链,走的正是马尔可夫过程该走的步子。连层长都落在该落的地方。碰层自由——正是 Bogoliubov 三阶段里动理学那一格。关联在 碰 的尺度上生成、又在同一尺度上被稀释,熬到层末已经所剩无几;这时候清零,清掉的本来就不多。

到此为止,我们做的还只是一个近似,不是一条定理。因为「每层重制备一次直积」这个动作本身一点也不新——它就是分子混沌,只不过换了个节奏。把这篇文章里的三代人按「这一刀砍几次」排一排,一眼就看得明白:

  • Boltzmann 1872:每算一次碰撞就砍一次,把刚攒起的关联当场清零。
  • Lanford 1975:只在 砍一次,此后一刀不砍。所以严格,也所以撑不过一个平均自由时间——不砍,就得把整棵树列到底,而树在疯长。
  • 邓煜–Hani–马 2024:每层砍一次。

单看这三行,第三行不过是把节奏调到了前两者中间:比 Boltzmann 稀疏,比 Lanford 频繁。而这一点新意也没有——Bogoliubov 那句泛函假设讲的正是同一件事,都在假定关联散了就不回来,照这样做下去,我们只是把一个老的启发式论证换个节奏重演一遍,一分严格性都没多出来。所以第三行还少着半句,而全部的新意,恰恰就在这半句里:每层砍一次——而每一刀砍掉的东西,一律不许扔。Boltzmann 砍完就扔,Lanford 干脆不砍,这里是砍完还要把误差留下来——下一层继续追踪,一分都不能丢。

误差要追踪,就得有分解方式。分解的对象是真实的 体分布——不是我们制备出来的那个假直积;假直积按定义完全独立,没有偏离独立的部分。取某个接口 层,把真实分布拆成「完全独立的那部分」加上「偏离独立的那部分」,并按是哪几位偏离了分门别类:

层主层主层𝓎层

读法很直白。第一项是完全因子化——这 位互不连通,各自只跟不相干的旁人打交道。求和那一项里, 是这 位分子的任一非空子集,𝓎 装的正是「 里这几位经由碰撞彼此连通」的那份关联,而 之外的人照旧独立地乘在旁边23。一句话:等号左边是真相,右边第一项是我们制备出来、马上要喂进下一层的那个假直积,求和那一堆就是两者之差。(严格说右边还挂着一个余项,原文说它小得可以忽略、后续分析里径直略去,这里从简。)

这个东西物理学家熟得不能再熟:它就是累积量(cumulant),就是连通关联函数,就是 Feynman 图里那些连通图。统计场论里我们天天做同一件事——把 点函数减去所有能拆成低阶乘积的部分,剩下的连通部分才携带新信息。这里一字不改,只不过「连通」的含义变得极其具体:存在一条碰撞链把这几位串起来。上一节里我们盯着的那些树和那些环,画的正是这个连通关系。𝓎 的定义性质也就一句话:当且仅当 里这几位彼此完全独立时,它为零。

目标随之明确,要证的是两件事:主项 主 一路贴着 Boltzmann 解;而所有的 𝓎 都小,并且牵扯的人越多越小,具体说,希望有 𝓎𝒹。这也正是「分子混沌」这句口号在数学上的精确形状:不是「关联为零」,而是「关联小,且随牵扯的分子数指数地更小」。

万事俱备,剩下的看起来只是例行公事。证明这个小量能从这一层传到下一层:第 层末尾小,则第 层末尾也小。层层递推,一路推到 𝓉,收工。然而这一步做不到。而且——这才是要紧的——它不是「暂时做不到」。

最能说明问题的是 Pulvirenti 与 Simonella 的工作:他们确实做出了一个传递,但入口要的是 的小,出口给的只有 的小。拓扑掉了一级,第二轮就接不上了。看上去像个技术缺口,等一个更会做估计的人来补。但邓煜–Hani–马指出:没有人补得上。对任何 Banach 空间都补不上。

反证漂亮得不像话,而且从头到尾是物理的。Newton 力学可逆:从第一层末尾那个时刻出发,往前走一层到 层,与往后走一层回到 ,这两件事只差一个速度反号。而任何自然的范数——、、,随便你造一个——在速度反号下都不变。于是,如果真有一个「范数小 范数小」的递推,它必须对前后两个方向一视同仁,于是同时给出

层层层

——荒唐,第二式的主项该是 才对。

现在把这件事翻译成物理,它的分量才显出来:范数是一个数。它只知道关联有多大,不知道关联长什么样。而「往前」和「往后」这两幅图景里,关联的大小一模一样,形状却南辕北辙。比如一对即将相撞的分子:往前看,这两位各自从自己的来路上飞来,履历互不相交是常态;要他们的履历相交(撞过同一个人,或者干脆就撞过彼此),得靠一次额外的巧合——小概率事件。分子混沌在这个方向上成立。往后看,把同一幅图倒过来放,这两位是刚刚撞完分开的。他们的履历必然相交——那次碰撞就明明白白写在数据里,概率是 ,谈不上什么小概率。分子混沌在这个方向上当场作废。同一个构型,同一份关联,同一个范数值;朝前朝后,一个成立一个不成立。所以想用一个数字把这份不对称传下去,从一开始就是不可能的。

于是不可逆性在这里露了一次脸,而这一次它不是作为哲学露脸,是作为一个具体的、指着能看的技术障碍:要传下去的不是大小,是结构。这大概是整项工作里最该讲给物理学家听的一段。一百五十年里,「不可逆从哪儿来」一直是个可以争到天亮的问题;而在这里,它第一次变得清晰起来:任何在时间反演下不变的量,都传不动这件事。一路走来我们说过好几次:Boltzmann 方程的不可逆,来自我们把关联扔掉。这话对了一半——扔关联是必要的,可它远远不够。其实物理学家早在上世纪三十年代就知道了:把稀薄气体换成另一种极限——每个分子感受到的不再是偶尔一次剧烈的撞击,而是所有旁人加在一起的一份微弱的平均力,。这是等离子体,是自引力星系,是所谓的平均场极限。推导它用的也是混沌传播(propagation of chaos,一字不差就是同一个词),做的也是把 维压成 维、只留一个单分子分布。样样都和 Boltzmann 那边一样。但得到的是

而这个方程在 、 之下严格不变24。没有 H 定理,没有熵增,没有箭头。它和 Newton 力学一样可逆。

同样地粗粒化,同样地扔关联,一个生出了箭头,一个没有。差在哪儿?差在被扔掉的那部分长什么样。平均场里相互作用是弱而集体的,根本没有「一次碰撞」这种事件可言,关联是均匀地、无结构地小,前后一个样子。稀薄气体里相互作用稀而剧烈,是一颗颗离散的碰撞事件;正是这份离散性,让关联带上了组合结构——就是上一节我们数过的那些树和那些环。

马骁把这件事说得最直接,而且他把它标为整项工作在物理上最重要的新观察:

系统的时间不可逆性,是累积量(cumulant)以及高阶 Feynman 图的特殊结构的推论。

于是前面那个说法可以修正得更准了。不是粗粒化生出了不可逆,是被粗粒化扔掉的那部分东西的形状生出了不可逆。扔,是必要的;扔掉的是什么形状,才是决定性的。

既然逐层归纳这条路被堵死了——那就绕开它。绕法说穿了倒是出人意料地朴素。回想那个 的发散是从哪儿来的:来自把整个分布一路往回展到 ,每退一次多一位分子、多攒一份 。可现在我们手上的东西已经分成两块了,那就分开处置:

  • 主项就地停住,一步也不往回展。它不需要展——我们要它干的事只是「贴着 Boltzmann 解」,而那是另一套办法,后面再谈。
  • 只把累积量往回展,而且一展到底,一路展回 。

妙处在于两边各得其所。发散只跟着主项走,主项不展,发散就无从发生;累积量倒是一路展到了底,可它自带 𝒹 那份小量,展得再狠也压得住。作者自评这种部分时间展开在此前的文献里没有出现过,并且明说它是「从可逆动力学得到不可逆方程」的关键一步。

于是先前那句话在这里兑现了:Lanford 头一回让混沌传播从信念变成命题;到这里,它更进一步,从「只在开头成立一次」变成「每一层都重新成立一次,且每一次的误差都记在 𝓎 上、可以一路追下去」。

现在,大道理讲完了,该铺开纸、掏出笔,看组合估计怎么落地了。

把一段碰撞史的几何全部抹掉,只留拓扑:画一张图,节点是一次碰撞(原文称 atom),边是分子的轨迹段,每个节点四条边(两进两出);同一个分子先后的几次碰撞按时间串起来,并且记着各自属于哪一层。这样一张图,原文称为 molecule,这里不妨叫它分子图。

而核心计数单位只有一个:图里独立环的个数,也就是图论里的圈秩(circuit rank),记作 。正是上一节那个环。全部的胜负,落在两个都以 为指数的量的赛跑上。他们用了一系列令人眼花缭乱的技术,算出来带 个环的图,数目至多是

𝒹

而每一张这样的图,贡献至多是

⟦图的那一项⟧层𝒹

一句话:每多一个环,要算的项数多一个对数因子,每一项要付的代价多一个幂次。幂次打败对数,级数收敛。整个全局结果的支点,说穿了就是这一行。

而这里发生了一次很值得玩味的地位翻转。上一节里,环是要排除的坏事——带环的碰撞史总测度趋于零,砍掉,剩下的树才凑成 Boltzmann 方程。如今换了个角度:正因为环要付一个 𝒹 的代价,环多的图才被压得住。环从「碍事的东西」变成了控制估计的参数。同一件事,短时里用来排除坏历史,长时里用来压住复杂历史。

至于难在哪儿——难在这份代价不是白拿的。环可以极其纠缠:一串重碰撞连锁发生时,逼出一个环会连带触发一堆环,而你不能指望每个环都独立地送你一个 𝒹。真正要命的是跨层的环:两层之间的枝叶粘在一起绕成的环,两头的信息不在同一层里,估计没法在单层内部闭合。

对付它的是一整套切割算法(cutting algorithm):按精心设计的顺序把一张大分子图逐级拆成若干种基本块,每种基本块按它能提供多少几何约束分成三类——有的不赚不亏;有的实打实地送来一个 𝒹,而这一类编码的恰恰就是一次重碰撞的代价(一条尚且自由的边被迫瞄准一条已经定死的轨迹,正是上一节那句话);还有的是亏本的,得靠前一类补回来。拆解的顺序至关重要——顺序选错,本该正比于环数的增益会退化成孤零零一个 𝒹,全盘皆输。跨层的环则要另立几个模型分别对付。

这一块是论文里最长最硬的部分,好在叙事上我们只需要知道它在干什么:保证「环数」这个组合数字,真的能一个不落地兑成 𝒹 的幂次。

要让这套机器跑得起来,还得先给动力学装上几道闸:每一层里团簇的尺寸、重碰撞的次数、擦身而过的次数统统强行截断,超出的部分整体丢进一个余项,再证明这个余项在 里随 𝒹 趋零。而处理「重碰撞次数超标」的那一步,用了一个粗糙得离谱却极其普适的老结果:在无界的自由空间里,给定数目的硬球,无论怎么摆、怎么给初速,总的碰撞次数都有一个只依赖于球数的上界——而且这个上界大得荒唐,是重指数量级的25。这就够用了——若少数几个球之间重碰撞太多,就必然有额外的分子被卷进环里,于是又多出一份自由度可赚。不过这是有代价的:那个结果吃的是自由空间的弥散性——球飞散开去就不再回来。换到环面上它就不成立了,绕一圈还会撞回来。所以整个证明搬不到周期边界条件上去,这个问题后面还要另算。

还剩一件事:主项凭什么一路贴着 Boltzmann 解?

框架就是 Lanford 的。一层之内,主项的初态是直积、时长短过收敛半径——这恰恰是 Lanford 定理的设定,而 Lanford 说的就是:这样跑出来的单粒子边缘贴着 Boltzmann 解。原文对这一半的全部评价也只有一句:还是短时理论那一套。可没这么便宜,有两处差异得处理清楚。

其一,输入本来就不一样了。Lanford 那边两条线同从 的直积出发,起点重合;这里不然——第 层开头喂进去的那个假直积,跟 Boltzmann 解在同一时刻的值早就分了家(上一层留下的差异)。所以每层的对照不是「相等」,而是「差多少」,且这个差还得带着往下一层走。

其二,碰撞算符差着一个 𝒹。硬球相撞时两个球心相距一个直径,Boltzmann 方程却把碰撞当作发生在同一点。这一处上一节其实打过照面,当时归在「一点琐碎的归一化」里带过了——在 Lanford 那儿它确实琐碎,一取极限就没了;可这里要迭代 次,就得知道它消失得多快,于是冒出一个 𝒹。定理条件里那个「空间梯度有界」,就是从这儿来的——它不是装饰。

这两处各自怎么估、误差又怎么逐层传下去,是纯粹的技术活,这里不展开。

还有一句值得单独抄下来:除了「Boltzmann 方程有正则解」这一条之外,整个证明没有用到 Boltzmann 方程的任何其它相消机制。这话的分量在于,这套办法对方程本身的结构要求低得出奇——它不关心你的解是怎么来的,也不指望方程自带什么好性质。

最后得交代一句:这套东西是从哪儿来的。答案是——从波那边来的。

邓煜与 Hani 此前解决的是另一个长期公开问题:从弱非线性的 Schrödinger 方程出发,长时间地推出波动理方程(wave kinetic equation)。那是波湍流理论的基本方程,在物理上的地位与 Boltzmann 方程在气体里一模一样——同样是物理学家凭直觉写下的,同样卡在短时间上。邓煜–Hani–马的摘要里自己就点明,这项工作「沿用前两位作者引入的范式」。

对应关系整齐得让人吃惊:Wigner 函数 ↔ 硬球气体的经验密度(两边都是随机场,都满足一条大数定律);波动理方程 ↔ Boltzmann 方程;而分子图(molecule)那套组合分类、「哪些变量才算真自由度」的挑法,全部源自波那边。差别也有:波那边的级数只条件收敛,得靠相位相消才活得下来;粒子这边绝对收敛,可几何约束硬得多。同一套方法论,在两个看似毫不相干的领域里各解决一个搁置了数十年的公开问题——这大概是 Fields 奖最典型的形状。

多年以后,面对评奖委员会,邓煜准会想起抱着 iPad 追《天才麻将少女》的那个遥远的下午。

世纪回响

说完好话,接下来的话必须说得同样清楚:Hilbert 第六问题没有被解决。被走通的,是「稀薄硬球气体」这条链路上的一段,而且是有条件的。清单如下:

  1. Boltzmann 方程本身的全局光滑性,仍然是个公开问题。这是条件性结果:方程的解活多久,定理就管多久;解要是在某个时刻爆掉,定理也到此为止。
  2. 接不上流体极限。定理要求的那个函数空间要分布在空间方向上衰减、总质量有限,这恰恰把流体极限最需要的设定(Maxwell 分布附近的扰动)排除在外。2025 年的姊妹工作处理了环面情形并一路接到 Euler 与不可压 Navier–Stokes,但要付出很强的正则性假设,而激波那类场景是已知会违反它的26。
  3. 环面还没拿下:上面那条弥散性用不上,得另造一套 cutting algorithm。
  4. 长程势完全开放。没有截断的长程相互作用里截面不可积,增益项与损失项各自都不收敛,非得靠两者相消才有意义——这一步现在还完全不知道怎么严格地做。
  5. 长时间的涨落与大偏差仍缺。短时间的涨落方程与大偏差泛函已经有了,长时间还没有。

从 Newton 到流体那条完整的链子,这项工作填上的是中间最硬的一环,而且填得很结实。两头都还开着口子。

最后说一件我认为最容易被误读的事。 这项工作在物理学界的第一反应,多半是这样的:Boltzmann 方程本来就是对的,实验验了一百多年,工程上算了一百多年,从来没出过岔子;数学家如今补上一个证明,那是数学家的事,物理上什么也没变。

这个反应错过了全部的要点。

因为这个证明真正回答的,从来就不是「Boltzmann 方程对不对」。它回答的是「不可逆是从哪一步、以什么形式进来的」——而这是个纯粹的物理问题,物理学家为它吵了一百五十年,从 Loschmidt 一直吵到今天,吵出来的多半是些不太好验证的说法:粗粒化、信息丢失、观察者的局限、宇宙的初条件。

而这里给出的答案,具体得出人意料。它不落在任何一个哲学词儿上,它落在关联的组合结构上,而且落到了可以数的地步:每一次碰撞,都把信息从低阶的累积量推向更高阶的累积量。低阶的那些——也正是我们真正观测的那些——因此越来越接近独立,Boltzmann 方程也就越来越准;而被推上去的那些高阶关联,以 𝒹 的速度躲进我们的分辨率以下。

于是那三个悬了一路的问号,在这里一并有了着落:

  • 「一次次扔掉关联,凭什么是对的?」——因为被扔掉的那部分确实小,小到 𝒹,而且牵扯的分子越多越小。
  • 「为什么如此粗暴的一步近似,换来的方程却精确得可以忽略误差?」——因为误差是按环数计价的,而每一个环都要付一个幂次的代价。
  • 「被忽略的误差,与被扔掉的关联究竟是什么关系?」——它们根本就是同一样东西,就是那些 𝓎,一个都没少。

信息一分也没有丢。Liouville 方程从头到尾可逆,这一点从来没变过。变的只是:信息一路往越来越高阶的关联里跑,而那里我们读不出来——不是原则上读不出来,是它躲得比我们的分辨率还快。

那段台球倒放的画面,现在可以给个更锋利的说法了。那 15 粒球倒着归位不违反任何定律,它需要的是一份恰好长成那个形状的关联。而这项工作做的,正是把「那个形状」量化了:要让气体倒着走,你得让那些本该按 𝒹 衰减下去的高阶关联反过来精确地长起来——而这正是宇宙的初条件里没有给你的那笔存款。

可惜的是,出于同样的理由,我们没法回到一百五十年前的奥地利,把这一切亲口告诉 Boltzmann。

革命尚未成功

上一节说的是这项工作没有做到什么。这一节往前看:接下来还有哪些路可以走,哪些坎还横在那里。

先说两条已经看得见路的。

其一,硬球只是个模型,真实分子之间当然不是硬球。硬球这个模型在数学上相当奇异:力全部集中在接触的那一瞬,之前之后一概为零。真实分子之间更像光滑势:近处强烈排斥,远处微弱吸引,Lennard-Jones 那副样子。所以问题很清楚——怎么把这套方法从硬球势推广成真实的势,而且无论对长时间还是短时间都成立。短时间那一层,这个推广早就完成了:Lanford 之后的三十年里,紧支撑的短程排斥势也被做了出来,与硬球并列。所以「Boltzmann 方程只对硬球成立」从来不是问题。

可长时间这一层还没有。邓–Hani–马目前把短程势列在「未来方向」里,说这是够得着的问题,同时也明说还需要新的技术零件。障碍在哪儿,正好照见前面那套记账机器:光滑势下,三个以上的分子有可能同时相互作用;而「一次碰撞=一对分子」恰恰是 molecule 的地基。图上每个节点四条边,两进两出,说的就是「两个分子进来、两个分子出去」。如果三体同时纠缠在一起,这个节点该画成什么样,就得重新定义。

所以准确的说法是:从硬球到光滑短程势,短时间那一层早已走通;长时间这一层还空着,只是作者认为它属于可以推进的方向。

其二,流体那一级台阶。

这是 Hilbert 第六问题的第二个箭头。严格说,Hilbert 当年真正想要的正是流体方程;Boltzmann 方程只是他点名的中间站。

怎么走,物理学家早就知道:在局域 Maxwell 分布附近做展开,也就是密度、流速、温度都随时空变化的 Maxwell 分布。领头阶给出可压缩 Euler 方程(、、 三个场),再往下一阶给出不可压 Euler 与 Navier–Stokes。这就是 Hilbert 展开与 Chapman–Enskog 展开,也正是那张 Bogoliubov 三阶段图里的第二级台阶:连速度分布的细节也扔掉,只留五个守恒场。

这一支的严格化,其实比第一个箭头做得好得多。邓–Hani–马自己的原话是:与 Boltzmann 方程的推导相比,流体方程的推导「要清楚得多」。这里有几十年的积累。但它们大多是从 Boltzmann 方程出发的,脚下并不直接接微观。Hilbert 第六问题要的不是两段各自严格,而是一条不断的链子:从 Newton 力学一路走到 Navier–Stokes,中间不换乘。所以,他们的姊妹论文做的正是这个接缝:在环面上做一个本文定理的类比,再与已有的流体极限结果拼起来。这样一来,Hilbert 纲领的两步第一次被连着走完了。代价是一组相当强的正则性假设;而激波那类场景,是已知会违反它的。那扇门开了一条缝,还没推开。

说完这两条看得见路的,剩下两道坎就没那么客气了:力伸得太远,和分子是量子的。

先说长程力。技术上的困难说来直白:长程势下的散射截面不可积。两个分子哪怕隔着老远擦身而过,也会被彼此轻轻偏折一下,而这种「掠射」的碰撞多到发散。于是 Boltzmann 碰撞项里的增益与损失各自都是无穷大,只有两者相减才有限。物理学家的对策是碰撞角截断:把小角度那一段一刀切掉。库仑力更惨——它连衰减都谈不上,得靠 Debye 屏蔽把有效力程掐在一个长度上,再拿库仑对数 兜住剩下的发散。这些做法在物理上管用,但还不是从微观力学推出的定理。要想在数学上严格化,麻烦还要更深一层。回想我们赖以记账的那套机器:把碰撞史抹成 molecule,节点是一次碰撞,环数 是计价单位,cutting algorithm 一刀一刀把它拆开。这整套东西,都架在「一次碰撞=一个离散事件」这个前提上。长程力干的事,恰恰是把这个前提取消掉。每个分子时时刻刻被无数个远处的分子轻轻推着,没有哪一瞬可以指着说「这里发生了一次碰撞」。没有 atom 可数,环也就无从数起。

而这件事的另一半,我们其实早就见过了。上一节那个 Vlasov 反例说的正是:把相互作用摊成弱而集体的平均力,得到的方程可逆,一点箭头也没有。真实的等离子体恰恰卡在两者中间——它满足的是 Vlasov–Landau 方程:一个集体的平均场项,加上一个稀有的近距离碰撞项。箭头只藏在后面那一半里。怎么把它从前面那一半的包围中严格地择出来,今天还没有人知道。

再说量子。方程本身早就有了,而且物理学家用了快一百年:把 Boltzmann 碰撞项里的四个分布函数按末态占据数修正一遍,就得到 Uehling–Uhlenbeck 方程(数学文献里多叫 Boltzmann–Nordheim,核物理里习惯叫 BUU):

碰

上取号玻色、下取号费米。多出来的那几个 就是量子统计的全部:末态已经有粒子在那儿,会反过来影响这一次碰撞跳不跳得过去。玻色是受激增强,越挤越想去;费米是 Pauli 阻塞,占满了就不许去。经典 Boltzmann 方程里没有这一项,因为经典粒子不在乎末态有没有人。

留意这里换了个符号。括号里那个 是个纯数,所以 也必须是纯数。它不再是我们一路用的那个「每单位相空间体积有几个分子」,而是每个量子态上平均住着几个分子,也就是占据数。两者相差的,正是「一个量子态占多大相空间」:

其中 是自旋简并度;至于那个 ,只因为本文的相空间一直按速度而非动量来量27。

请留意,当初我们提到说 不过是「某个用来调整量纲的相空间体积单位」,取什么都行;反正外面套着 ,换个取法只平移熵的零点。可到了这里,它不能随便取了:量子力学亲手把相空间切成了格子,一格 ,不多不少。当初那个可有可无的摆设,如今是硬的。 还顺手给了一个能直接估的判据。它的量级正是简并参数:数密度乘上热 de Broglie 波长的立方。于是「什么时候必须用量子版」不再是一句定性的话,而是一个能算的数: 是稀薄或高温,那几个括号退成 ,BUU 当场退回 Boltzmann; 逼近 ,量子简并开始,括号就非留不可。

这个方程管用得很。重离子碰撞的输运模型靠它,半导体里的载流子输运靠它,超冷原子气体的弛豫也靠它。它的平衡解自动是 Bose–Einstein 与 Fermi–Dirac 分布,而不是 Maxwell;H 定理照样成立,只是熵的形式换了。

但从多体 Schrödinger 方程严格地推出它,基本上还是一片空白。难在哪儿?在下面三条:

  • 轨迹没了。「一条尚且自由的边被迫瞄准一条已经定死的轨迹」——这句话是整套 cutting algorithm 的地基,而量子力学里根本没有轨迹可言,这个句子连意思都说不通。
  • 全同性从记号方便升格成了物理本体。我们当初请出 ,只是为了把编号洗匀、好让求和塌掉;在量子里,交换对称与反对称是波函数的硬性质,不是记账的方便。
  • 「初始混沌」该怎么写?经典这边是一个干净的直积态。量子这边直积态叫 Hartree 或 Fock 态,看着像,可纠缠并不是关联的简单推广:两个粒子可以毫无经典关联却纠缠得一塌糊涂。那么「初始无纠缠」到底该假设成什么,本身就是个问题。
    不过量子这一支并不是白纸一张。邓煜自己早在另一头留了半个脚印。

波湍流那条线本来就是量子的:他与 Hani 出发的那个弱非线性 Schrödinger 方程,正是玻色气体的平均场描述。那边没有轨迹,也就不可能有几何约束;他们找到的替代品是相位相消——级数只条件收敛,靠相位的相互抵消才活下来。所以「量子里没有轨迹怎么办」这个问题,其实已经有过一个答案,只不过那个答案是给波准备的,还没人知道怎么搬到粒子这边来。

玻色那边还埋着一个漂亮的意外。这个方程的玻色解会在有限时间里塌掉:空间均匀、各向同性,而且粒子数与能量密度落在超临界区间时,分布会在零动量处堆起一个 峰,方程本身在那一刻失效。这不是猜测,是已经证明了的28。而邓–Hani–马的定理,恰恰是「解活多久,我管多久」的句式。于是这里出现了一件很微妙的事:那个「爆掉」不是证明的技术缺陷,它就是玻色–爱因斯坦凝聚本身。动理学方程在它面前失效,是物理,不是数学的失手。

最后,把最大的那个类比亮出来。我们这一路证明的核心,压成一句话是:关联一路往越来越高阶跑,跑得比我们的分辨率还快。 信息一分没丢,只是躲进了没人读得出来的地方。

量子那边,同一句话有另一个名字。先说本征态热化假说(ETH)。它常被概括成「混沌的量子系统会热化」,可这句话什么信息也没给。它真正断言的是一件相当激进的事:对一个混沌的多体系统,取它任意一个能量本征态,随便挑一个局域观测量 ,那么

微正则⟦随体系增大而指数衰减的涨落⟧

也就是说,单独一个本征态,自己就已经「看起来是热的」了。热力学系综不是从很多个态平均出来的,它早就藏在每一个本征态里面。于是热化根本不需要什么额外机制:从任意初态出发,各个本征态之间的相位随时间乱掉,局域观测量自然就趋向那个人人共有的系综值。

这与分子混沌是同一类东西:一条关于「典型性」的假设,有压倒性的数值支持,却几乎没有证明。它们在各自领域的地位,乃至各自面临的怀疑,都惊人地相似。

而它们说的其实是同一件事。ETH 里,初态的全部信息一分不少地保存在各个本征态的相位里;演化是幺正的,Liouville 定理的量子版一步没让。可你要把它读出来,需要指数精细的分辨率。信息跑到了读不出来的地方去,跟本文这一路讲的一模一样。

这句话还有一个可以直接量出来的版本:纠缠熵的线性增长。取系统的一小块,算它与其余部分的纠缠熵,你会看到它随时间线性地涨,直到饱和在热力学熵那个值上。「信息往越来越大的子系统里跑」——这就是它的字面量化,是我们那个 𝒹 在量子世界里的对应物。

这些话如今不再只是纸上的推演,因为有了一个干净的实验台:超冷原子气体。

它几乎是为这类问题定做的:Feshbach 共振可以把相互作用强度从零一路调到无穷,也就是把「稀薄参数」握在手里;飞行时间成像可以直接把动量分布拍下来,也就是直接测 ;一次淬火之后系统怎么弛豫,也可以一帧一帧看着它发生。Boltzmann 方程管不管用,不再只是算出来的,也是看出来的。

而它也立刻交出了一个反例。把气体压进一维、让它近乎可积,那么它不会热化。著名的「量子 Newton 摆」实验里,两团原子在一维管里来回对撞了成千上万次,动量分布始终保持着初态的记忆,就是不肯弛豫到热平衡29。后来的理论工作把这类系统的归宿讲清楚了:它们停在一个由无穷多守恒量共同决定的广义 Gibbs 系综上,而不是通常的热平衡。可积性在量子这边扮演的角色,正是 Kac 环在经典那边扮演的角色:一个明明白白的提醒——「会热化」从来不是定理,它是一句需要前提的话,而前提可以被违反。

Hilbert 当年提的那个问题,字面上要的是「物理学的公理化」。一百二十六年过去,我们刚刚把最简单的那一罐气体做完:硬球、稀薄、经典、无外场。

革命尚未成功。


Footnotes

  1. 此处使用了略微不同的记号: 是 Boltzmann 方程里那个「数分子」的分布, 是 Liouville 方程里那个「1 号分子」的分布;两者相差一个 的因子,详见下文。 ↩

  2. 这里的 是某个用来调整量纲的「相空间体积单位」,在外边取 时,其实没有任何影响(只平移了零点),但一般而言常见的取法是 . ↩

  3. 同理,「数对」的分布是 。这些 的因子眼下只是归一化因子;等后文把 的极限当真时,它们将与碰撞截面合流成 Boltzmann–Grad 标度,正是严格化的要害之一。 ↩

  4. 阻尼力那种依赖速度的「力」不是微观粒子之间的基本相互作用,而是已经做过平均之后冒出来的有效场;若坚持逐个粒子的语言,它本来就不该出现在这里。唯一值得认真单独看一眼的是 Lorentz 力:。这时把速度导数写成散度,可以证明多出来的那项仍然为零;所以它并不破坏这一处积分分部。只是「Lorentz 力」本身已经是把电磁场当作给定背景的简写;若完整地把带电粒子与电磁场一起纳入统计力学,相空间还要包含场自由度,例如矢势 与势函数 (以及相应共轭变量),问题立刻变成无穷多自由度的统计场论。严格处理要进入泛函空间,连天然测度都很难自然给出,复杂度已经远远超出本文讨论的 BBGKY 硬球图像。 ↩

  5. 这里默认分子间是处处有限的光滑对势。Lanford 与邓–Hani–马真正处理的硬球气体里,「力」集中在接触的一瞬、是奇异的,层级得改用碰撞面上的边界项来写——形式更繁琐,精神完全一致。 ↩

  6. 往下一路推导,物理学家的「当然可以」还会反复出现;就当诸君都不是挑剔的数学家,不再一一唠叨。 ↩

  7. 在 维空间里,碰撞管的截面是 𝒹,故标度写作 𝒹 保持有限;正文取 。这条低密度标度出自 Grad 1958 年为《物理学手册》所写的动理学理论长文。 ↩

  8. 技术上还要再松一道:连分子的总数都不必钉死,而是让它服从一个 Poisson 型的分布(所谓巨正则系综)。这在低密度下处理关联时要方便得多。 ↩

  9. 束流的能量不能真取成 函数:定理要求 初 是连续 Lipschitz 的概率密度、且速度方向带 Gauss 尾,把束流抹开一丁点即可。 ↩

  10. 不过不能简单说一句「直积态测度为零」就了事。这里的舞台是 个自变量的函数所成的空间,而无穷维空间上并不存在天然的平移不变测度(Lebesgue 测度在此没有对应物),「零测」得先人为挑一个先验才有意义,不同的先验给出不同的零测集。有限维的类比倒是干净: 个粒子的联合分布住在一个维数随 指数增长的单纯形里,而直积态只由一个单粒子分布决定,在其中不过是一个低维的代数子簇。 ↩

  11. 顺带一提,这个演化算符在相空间的平方可积函数上是幺正的,而幺正性的来源恰是相流保容——Liouville 定理换上算符的行头,就是 Koopman 与 von Neumann 那套经典力学的算符表述,其生成元正是 Liouville 算符。 ↩

  12. 两处细节。其一,这两个比例天然就差一个 : 里混进了「同一个分子既算进 又算进 」的项,而 把两个编号相同的排除在外;所以「有关联」指的是差得远比 大。其二,差别未必在随便哪两块区域上都显眼——对 Loschmidt 那个反演态,你得挑对区域,比如「两球正迎面即将相撞」那一类;在那儿差别是压倒性的,而那正是重碰撞的种子。 ↩

  13. 这里约定 、。当 时,多重积分、碰撞算符以及中间的演化算符都视为空,算符链就退化为 ,微分乘积则视为 ;省略号中,演化算符与碰撞算符交替出现,每段演化的时长都是相邻两个碰撞时刻之差。 ↩

  14. 整条 BBGKY 层级都是线性的——它本就是 Liouville 方程换个写法,而 Liouville 方程对 线性。Boltzmann 方程那个出名的二次非线性,要等到闭合那一步、把 换成 时才诞生。 ↩

  15. 「重碰撞」(recollision)这个词在文献里的范围比字面宽得多:判据是树没有规定的碰撞,任意一对分子都算,不要求它们此前撞过。Gallagher–Saint-Raymond–Texier 给的是「存在任意两个不同的分子与某个 ,使它们的两条轨迹相距不超过 𝒹」;Bourbaki 讲稿则直接写作「除碰撞树规定者之外不发生任何碰撞」。之所以以环计数而不以「重逢」计数,是因为代价正来自闭环本身:每闭一个环,就有一条自由的边被迫瞄准一条已被图定死的轨迹。这个计数法一直沿用到邓–Hani–马那里——他们对每张碰撞图给出的小量因子是 ,指数上的 数的正是环数。 ↩

  16. 这里说的是机制,不是证明。真做起来相当微妙:刚在上一个节点分开的两位本就相距 𝒹,并不处在一般位置,不能简单地「把坏集合加厚一圈」了事。Gallagher–Saint-Raymond–Texier 明确警告过这一点,并用一整章专讲它——先证「好构型在新粒子加入后仍是好构型」,再用几何引理把导致重碰撞的那些速度与偏折角的集合估出来。 ↩

  17. 所谓「琐碎的归一化」,可以这样理解:两侧的分布本就差 倍(),于是两个碰撞算符也差一个 。 ↩

  18. 说「落了地」要稍加限定:Lanford 1975 的讲义给的是纲领与关键估计,此后经一系列工作直到 2000 年代初,论证的细节才补齐、结论才被推广和精细化。 ↩

  19. 这一条尚非定理:讲稿只说「数值与形式论证表明可以预期」,并指出它与超临界随机图的强类比。 ↩

  20. 这条路上的名字很长一串,硬球之外还推广到了短程光滑势。 ↩

  21. 确切的条件是:存在 与常数 ,使 对一切 𝓉 成立,其中该范数在速度方向要求 Gauss 尾、在空间方向对单位胞逐个取上确界后求和(于是要求总质量有限、分布在空间上衰减)。正是这最后一条把「Maxwell 分布附近的均匀扰动」排除在外,见后文边界那一节的第 2 条。 ↩

  22. 留意这一步是非线性的——非线性正是在第 2 步「张量幂」上诞生的。这与前文那处完全同源:整条 BBGKY 层级本是线性的,Boltzmann 方程那个出名的二次非线性,也是在把 换成 的那一刻才诞生。每一次「重捏直积」,都是一次闭合,而闭合就是非线性的来处。 ↩

  23. 留意 那几项并不为零:𝓎主。它量的当然不是「一个分子的关联」(那说不通),而是参考廓型自身的偏差——上一层留下的关联照样在演化、照样反过来改动单粒子分布,加上截断丢掉的那些大团簇,都记在这一项上。若把参考廓型取成真实的 ,这一项确实归零(累积量的标准约定),可那样一来递推就不闭合了:你得先知道 才写得下分解。𝓎 正是为「主项必须能往前算」付的价钱。 ↩

  24. 验一下即可:令 逆(这里不能用波浪号——波浪号在本文是「Liouville 侧」的标记),则数密度 逆 在 处等于原来的在 处,于是平均力也照样翻转,三项各自变号,方程原样成立。这与 Newton 方程的时间反演对称是同一件事——本来就该如此,因为平均场极限根本没有丢掉任何「方向」。 ↩

  25. Burago–Ferleger–Kononenko 1998:无界空间中 个硬球的总碰撞次数不超过 量级。这个界不依赖球的大小、质量与初始构型,粗糙但普适,正合此处所需——这里要的不是好的界,是一个有限的界。 ↩

  26. 这条与正文主线关系密切,值得单说:从 Boltzmann 方程到 Euler/Navier–Stokes 是 Hilbert 第六问题的第三个箭头,而它在激波处必然出问题——激波是解的间断,任何要求强正则性的推导都在那里失效。所以「Newton → Boltzmann → 流体」这条链子,眼下是中间一环最结实,两头都还软着。 ↩

  27. 本文的相空间体积 一直按 来量(通篇写作 ),而 是 的体积,两者差一个 ;换到动量变量上,因子就干净地是 。开篇脚注里说 「常见的取法是 」,隐含的正是动量那一版——那里因为外面套着 、常数只平移零点,这点出入无伤大雅。 ↩

  28. Escobedo 与 Velázquez 的一系列工作。确切的条件是:空间均匀、各向同性、初值有界,且粒子数密度与能量密度落在「对应平衡态含有 Dirac 质量」的那个区间里;此时解在有限时间内于 范数下爆破,弱解则在原点处生成一个 Dirac 测度。 ↩

  29. Kinoshita、Wenger 与 Weiss 2006 年的实验。留意原始论文的论点只是不热化;「归宿是广义 Gibbs 系综」是此后一批理论工作(Rigol 等)建立的,两件事不宜混作一谈。 ↩