好久没写过长文了，上一次对Wang先辈的论文做了不少解读，在我的主页可以找到零零散散的汇总，我们废话不多说，这次来看看Deng的。

阅读提示：本文大量数学史私货（没说没有其他私货），所以旁边有目录，不喜欢看庞大数学史论述的可以看一眼第一节标出来的希尔伯特第六问的三层级，直接跳到后面看对工作的解读。

希尔伯特第六问题

希尔伯特之梦：希尔伯特及其事业、斗争

"In der Mathematik gibt es kein Ignorabimus."
——Hilbert

19世纪末，康托尔创立了朴素集合论，为整个现代数学提供了统一的底层语言。希尔伯特对其推崇备至，称之为“人类纯粹智力活动的最高成就”。

在1900年之前，罗巴切夫斯基、黎曼等人发展的非欧几何打破了欧几里得几何作为唯一物理现实的垄断地位。数学家们意识到，几何学实际上可以脱离直观的物理空间，成为一套基于公理的纯逻辑演绎系统，1899年，希尔伯特出版了划时代的著作《几何基础》，以前所未有的严密性将几何学彻底公理化。在希尔伯特的哲学观中，几何学本质上是物理学的一个分支——既然带有强烈经验色彩的几何学可以被剥离经验、转化为严密的纯数学公理系统，那么物理学的其他分支理应也能接受同样的数学处理。

在这种时代精神的感召下，1900年巴黎第二届国际数学家大会上，他提出了23个数学问题，深刻地塑造了20世纪数学的发展轨迹，第六问题核心诉求是将物理学公理化，即将物理理论剥离为三个独立部分：物理公理、分析工具以及物理学解释。

第一层级的目标是公理化概率论与力学，希尔伯特在原稿中明确指出：“在这些（需要公理化的）物理科学中，排在首位的是概率论和力学。” 他要求找到一套完备、独立且不矛盾的公理系统，作为这两个学科的绝对起点。这是因为19世纪下半叶，麦克斯韦和玻尔兹曼创立了气体动理论和统计力学，这一理论试图用大量微观粒子的牛顿力学运动来解释宏观的流体力学和热力学现象，然而，这一跨越面临着严峻的逻辑和数学困难：微观粒子遵循确定性的经典力学，而宏观演化却引入了概率和不可逆性。当时的概率论仍然被视为一种处理赌博和经验数据的边缘工具，缺乏严密的数学基础，甚至充满悖论。希尔伯特敏锐地察觉到，如果不从数学上严密定义概率，并理清从微观到宏观的极限过渡，物理学的根基将是不稳固的。

第二层级目标是从微观到宏观的数学极限证明，这是希尔伯特对玻尔兹曼理论的直接回应，他要求数学家从原子层面的微观牛顿定律出发，在粒子数 N \to \infty 且体积缩小的极限下，严格推导出介观的玻尔兹曼方程，并进一步推导出宏观的流体动力学方程，如Euler或Navier-Stokes方程，本质上是要求证明：物理学的多尺度现象（微观、介观、宏观）在数学上是同源且逻辑自洽的。

终极目标是消除物理学中的隐性假设，希尔伯特希望通过公理化，彻底查明物理理论中隐藏的前提条件。一旦物理学变成了纯粹的数学推导，科学家就可以通过检验公理的相容性，来判断物理理论是否可能存在自相矛盾。

当然，此时的希尔伯特还不知道自己要面对什么。

风云多变，1902年，伯特兰·罗素提出了著名的罗素悖论，罗素将这个悖论写信告诉了即将出版《算术基本法则》第二卷的弗雷格，直接摧毁了弗雷格毕生的逻辑主义心血，集合论中的悖论如布拉利-福尔蒂悖论、康托尔悖论接连浮出水面，标志着第三次数学危机的全面爆发。

天下大乱，形势大好。

要理解希尔伯特如何在这场纷乱的斗争中走上舞台，并完善自己的思想，就要理解希尔伯特后半生的死对头——布劳威尔。

布劳威尔是一位带有神秘主义色彩的思想家，1905年，也就是在他24岁时，他出版的短篇专著《生活、艺术与神秘主义》是他神秘主义思想的集中表现。他主张人类应当从物质世界和外在现象中撤离，退回到自我意识的封闭堡垒中，外部世界在本质上是一种幻觉或次要的衍生物；同时，在他的描述中，当个体闭上双眼，切断与外部世界的联系时，就能体验到一种神秘的、无时间性的宁静。

在布劳威尔的宇宙观中，人类历史是一场悲剧性的堕落史，而导致堕落的罪魁祸首正是理智和语言，他认为，人类为了在自然界中生存和统治其他客体，发展出了理智，而理智，就其本身而言，其并非为了探索真理而存在的官能，理智将世界客体化、割裂化，它关心的只是“因果序列”和“控制”。他认为语言纯粹是人类为了意志的交锋、为了驱使他人或达成社会协作而发明的工具，语言无法传达真实的内在体验，一旦我们试图用语言去描述那个神秘的内在真理，真理就已经被扭曲和贬低了。

这种思想深深渗透到了他后来的数学构造实践中，布劳威尔将康德的先验时间观推向了神秘主义的极致，他认为，当静止的内在意识发生第一次波动，即从现在流向过去，从统一分裂为二时，就产生了时间的直觉，所有的自然数、连续统以及整个直觉主义数学大厦，都是在这个无语言的、纯粹的内在时间直觉中一步步构造出来的，数学本质上是一场无语言的、内在的神秘冥想活动，数学公式和论文，只是数学家为了唤起其他数学家相似的内在体验而使用的、并不完美的暗示符号。

1907年，布劳威尔在阿姆斯特丹大学发表博士论文《论数学的基础》，首次系统性地提出了直觉主义的哲学框架，标志着布劳威尔正式将其神秘主义和直觉主义主张付诸数学实践。

最开始，希尔伯特和布劳威尔是亲密的朋友，1909年，两人在荷兰的斯赫弗宁根海滨初次相遇，这次会面给年轻的布劳威尔留下了极其深刻的印象，他在给朋友的信中深情地形容希尔伯特是“穿透我生命的一道美丽的新光芒”。

布劳威尔早期在拓扑学领域取得了震撼古典数学界的成就，如不动点定理、维数拓扑不变性等，这些至今都还是代数拓扑学的必学内容。希尔伯特对这位年轻人的数学才华赞叹不已，并在1912年鼎力支持布劳威尔获得了阿姆斯特丹大学的正式教职，希尔伯特甚至曾试图将布劳威尔挖到当时的“数学麦加”——哥廷根大学，不过被布劳威尔出于对远离喧嚣的隐居生活的向往而婉拒了。

但是，布劳威尔的野心越来越膨胀，直觉主义革命的他正在不断蚕食希尔伯特形式主义哥廷根派的空间，他开始要求数学界剔除所有不具备构造性的证明，这就要求放弃经典逻辑中的排中律，这意味着希尔伯特最为倚重的大量理论将被视为无意义的符号游戏。他宣称，大部分基于非构造性证明的现代数学（包括希尔伯特早年的成名作不变量定理）都是毫无意义的废话，必须被抛弃，二人开始互生罅隙。

最让希尔伯特难堪的是，外尔，这个希尔伯特最得意、最聪明的门生。他在1918年出版《连续统》，开始偏离古典数学轨道，试图在不使用非构造性方法的前提下，重新建立实数理论，哥廷根内部开始出现裂痕，三年后，1921年，外尔突然倒戈，发表了支持直觉主义的宣言《关于数学基础的新危机》，宣称“布劳威尔就是革命”，这对晚年患病的希尔伯特造成了极大的打击。

这期间，希尔伯特在1920年发表了希尔伯特纲领，正式把他从十九世纪末期至二十世纪初想要存在的稳定地基做了一个系统的论述，我们称之为希尔伯特之梦，他要建立一个绝对安全、严密且全知全能的数学乌托邦，其核心定义包含四大支柱：

1.彻底的形式化： 将所有数学命题、概念和证明完全转化为没有任何经验意义的符号组合与逻辑演算，剥离人类的直觉和模糊性。数学变成一场“按照既定规则操作符号的游戏”。

2.绝对的相容性： 必须在数学体系内部，严格证明这个形式系统是绝对不会推导出矛盾的。这是“梦”的根基。

3.逻辑的完备性： 任何一个真实的数学命题，都必须能够在这个形式系统内被推导和证明出来。即数学中不存在“绝对的不可知”。

4.普遍的可判定性： 存在一种机械的算法过程，能够在有限步内判定任何给定的数学命题是真还是假。

1922年，希尔伯特正式吹响了反击的号角，发表了题为《数学的重新基础：第一篇报告》的著名论文，他开始着手构建希尔伯特纲领的梦想，在《数学的重新基础》中，希尔伯特明确将矛头对准外尔和布劳威尔，他指出，如果按照直觉主义的要求去阉割数学，废除排中律和非构造性存在证明，那么数学家将失去微积分、康托尔的集合论以及绝大部分现代分析学。希尔伯特在这篇论文中确立了一个核心基调：不能退缩，必须一劳永逸地为古典数学建立一个绝对安全的基础，彻底平息关于悖论的恐惧。

这篇论文最重要的理论创举，是希尔伯特在数学内部划定了一道界限，提出了对象语言与元语言的分离。

希尔伯特提出，要彻底切断数学符号与直观经验之间的联系。数学命题（如微积分、集合论）不再是关于现实世界或抽象实体的陈述，而仅仅是一堆没有内在意义的符号串（如 \exists x, \forall y）。 数学定理的证明，不再是逻辑上的“推理”，而是按照一套严格规定的、机械的推导规则（如分离规则），对这些符号串进行物理上的排列组合与变换。

数学本身变成了没有意义的符号游戏，那么谁来赋予整个体系合法性？希尔伯特由此创立了一门全新的学科——元数学（或称证明论）。元数学的对象是数学证明本身，希尔伯特提出，我们要站在形式系统之外，把写在纸上的那些符号和推导步骤作为物理实体，用一种极其直观、极其安全的逻辑，来研究这些符号系统的性质。

希尔伯特清楚，如果在元数学中继续使用复杂的古典逻辑，就会陷入循环论证，因此，他在此文中规定：元数学的推理方法，必须是极其严格的“有限主义”的。在元数学中，绝对不允许使用任何关于“无限”的概念。所有的元数学推理，必须是直观的、构造性的、在有限步骤内可以完成的验证。这实际上是希尔伯特向布劳威尔做出的巨大让步，他等同于同意直觉主义者关于真理必须是有限构造的哲学观点。但把这种有限主义限制在元数学层面，只要用他们认可的有限构造方法，证明了包含无限的古典形式数学系统是安全的，他们就再也无话可说了。

在确立了形式化和元数学框架后，希尔伯特在论文中提出了纲领的终极目标，也是整个系统合法性的唯一判据：相容性证明。在《数学的重新基础》中，他将这一目标转化为一个纯粹的符号组合学难题：

1.假设存在一个形式系统（包含了所有的古典数学逻辑和算术公理）。

2.我们需要在“有限主义”的元数学框架内证明：在这个系统中，绝对不可能同时推导出公式 A 和它的否定 \neg A。

只要完成了这个有限主义的相容性证明，排中律和实无限就可以作为一种理想元素，合法地保留在数学中。悖论的幽灵将被彻底驱散，康托尔的乐园将得到永恒的保卫。

1922年的《数学的重新基础》标志着希尔伯特从防守转向了战略反攻，然而，希尔伯特的身体也在逐渐恶化，到了1928年，希尔伯特患上了严重的恶性贫血，病情危机。

《数学年刊》是当时世界上最具权威的数学期刊，希尔伯特是主编，而布劳威尔也是核心编委之一。他在病榻上产生了一个极其强烈的恐惧：一旦自己离世，凭借布劳威尔在学术界的巨大威望和强硬霸道的性格，哥廷根学派布卢门塔尔等人将无法压制他。布劳威尔极有可能会实质性掌控《数学年刊》，并利用主编的权力，将直觉主义确立为期刊的审稿标准——这意味着，使用了排中律和非构造性证明的古典数学论文将被大量拒稿。因此，他决定在自己死前，动用一切政治手段，彻底清除布劳威尔。

那一夜他思索良久。他想到了千千万万在经典数学奋斗的数学家，难道他们为之奋斗的一切，就换来这样一个数学界吗？他想到广大的数学学生，他们怎么办？他很不甘心。他做出了一个决定，这件事哪怕是他摔得粉碎，即使失败了，他也不惧怕。因为他知道：“理性的风终究会吹走尘埃”。

希尔伯特从哥廷根出来了。在召集学者们开会时，说了一句惊世骇俗的话：“以前我带你们抗击直觉主义，现在我又要带你们抗击直觉主义！”

他发动了一场轰轰烈烈的《数学年刊》编辑部政变，1928年10月，希尔伯特在未通知布劳威尔的情况下，单方面向出版商斯普林格以及编辑部其他成员发出了一封强硬的信件。他宣称，由于布劳威尔的直觉主义立场与期刊的宗旨完全不相容，他拒绝继续与布劳威尔共事。希尔伯特下达了最后通牒：要么布劳威尔走人，要么我希尔伯特退出《数学年刊》。

作为共同主编的爱因斯坦收到了这份要求罢免布劳威尔的信件。爱因斯坦对此感到极其厌恶，他拒绝在罢免信上签字。爱因斯坦给波恩写信，留下了著名的评价：“这简直是一场青蛙与老鼠的战争！” 随后，爱因斯坦为了表达对这种学阀作风的抗议，主动辞去了《数学年刊》编委的职务，抽身退步。1928年11月2日，另一位著名数学家卡拉西奥多里试图在两人之间斡旋，建议布劳威尔主动辞职以保具体面。但生性骄傲、自认没有做错任何事的布劳威尔断然拒绝，并发表长信痛斥希尔伯特的行为是“不公正的政变”。

1928年11月27日，在希尔伯特的强大压力和出版商斯普林格的默许下，执行主编布卢门塔尔被迫充当了刽子手。他给布劳威尔写了一封极其简短冷酷的信件，正式通知他：为了期刊的生存，经过希尔伯特的决定，你已被解除了编辑职务。为了在程序上使这次驱逐显得合法，希尔伯特和斯普林格出版社甚至采取了一个极端的手段，他们宣布直接解散整个《数学年刊》编辑部，然后立刻重组。 在重新公布的编辑名单中，希尔伯特依然在列，而布劳威尔的名字则被永远剔除了。

遭到无理驱逐的布劳威尔陷入了极度的抑郁，他感到被整个数学界背叛，在此后的长达十多年的时间里，即所谓的“沉默的十年”，他停止了所有数学论文的发表，撤回了自己在拓扑学和直觉主义上的前沿阵地，转而研究哲学，甚至变得孤僻和偏执。

斗争以一种最滑稽的两败俱伤走向了结束

1930年9月8日，希尔伯特在家乡发表了那场著名的退休广播演讲，以“我们必须知道，我们必将知道”作为唯理主义的最高昂宣誓；而就在此前一天的圆桌会议上，24岁的库尔特·哥德尔首次宣布了他的第一不完备性定理。

希尔伯特并没有参加哥德尔发言的那场圆桌会议，将这个噩耗带给哥廷根学派的，是当时在场的约翰·冯·诺依曼，以及后来仔细研读了哥德尔论文的保罗·伯奈斯。

据伯奈斯回忆，当希尔伯特最初听到哥德尔的结论时，他的第一反应是“愤怒”与深深的挫败感，在很长一段时间里，希尔伯特拒绝完全接受哥德尔的结论，试图让伯奈斯去寻找证明中的逻辑漏洞。

第一定理：哥德尔证明了，在任何包含基本算术的无矛盾形式系统中，必然存在一个命题，既不能被证明为真，也不能被证明为假。

第二定理：哥德尔证明了，如果一个系统是相容的，那么这种相容性绝对无法在系统内部被证明。

这意味着，希尔伯特试图建立的“绝对安全的基础”在逻辑上是一个悖论。

面对无可辩驳的数学证明，希尔伯特最终表现出了一代宗师的体面，他不再寻求“绝对的有限主义证明”。1936年，哥廷根学派的新锐格哈德·根岑通过引入“超限归纳法”证明了算术系统的相容性，这虽然偏离了希尔伯特最初的严苛标准，但被视为希尔伯特纲领在哥德尔打击下所能达成的最好妥协。

布劳威尔对此态度冷淡，他一直主张，数学是人类心智的自由直觉活动，语言和符号只是对这种心智活动的拙劣模仿。对布劳威尔而言，哥德尔的不完备性定理本质上仍然是一个关于形式系统的定理。既然布劳威尔打心底里认为形式系统只是没有生命的纸面游戏，那么去证明一个毫无意义的游戏是不完备的，又有什么哲学价值呢？

哥德尔既不是希尔伯特的形式主义者，也不是布劳威尔的直觉主义者，他是一个坚定的柏拉图主义者，哥德尔的搅局，让这场纷争以一种极不干净的结尾收场了。

余梦难消

1933年的哥廷根，空气中弥漫着政治的狂热与学术的死寂。后来，当纳粹教育部长向晚年的希尔伯特询问哥廷根数学研究所的现状时，这位曾经的统帅只留下了一句凄凉的回答：“哥廷根已经没有数学了。”肉体上的哥廷根在纳粹的铁蹄下化为焦土，精神上的哥廷根在哥德尔的定调中轰然倒塌。对于纯数学而言，希尔伯特的那个“绝对无矛盾”的古典之梦，似乎已经随着那个时代一起被连根拔起。

这一年，是哥廷根学派覆灭的元年。

但许多同时代的人并未察觉，同样是在1933年，远在莫斯科的苏联数学家柯尔莫哥洛夫出版了薄薄的《概率论基础》；而就在此前一年的1932年，流亡前夕的冯·诺依曼出版了《量子力学的数学基础》。

尽管希尔伯特企图用公理化驯服纯粹算术的战役以一种极不干净的方式收场了，但他试图用公理化驯服“经验世界”的宏大构想，却恰恰在哥廷根覆灭的前夜迎来了大爆发。

柯尔莫哥洛夫的《概率论基础》天才地将概率论与勒贝格的测度论结合，他提出了概率空间的三元组 (\Omega, \mathcal{F}, P)，其中 \Omega 为样本空间，\mathcal{F} 为 \sigma-代数（事件），P 为概率测度。通过五条极简的公理，柯尔莫哥洛夫将概率论完全纳入了现代纯粹数学的范畴。这一工作被数学界公认为完美解决了希尔伯特第六问题在概率论方向的诉求。

1900年时量子力学尚未诞生，但当1920年代海森堡（矩阵力学）和薛定谔（波动力学）建立量子理论后，狄拉克给出了一个物理上优美但数学上不严密的形式体系（如引入了当时的数学无法定义狄拉克 \delta 函数，虽然后来在泛函分析的广义函数论解决了这个问题）。量子力学的公理化自然成为了第六问题的新核心。

《量子力学的数学基础》在希尔伯特等人的谱理论基础上引入了可分希尔伯特空间 \mathcal{H} 的概念。他将量子力学的状态公理化为 \mathcal{H} 中的单位射线，将物理可观测量公理化为 \mathcal{H} 上的自伴算子，并利用谱定理严格统一了薛定谔的波动方程和海森堡的矩阵力学。这一工作在非相对论量子力学领域，不仅解决了公理化问题，也直接推动了泛函分析和算子代数的繁荣。这是希尔伯特第六问的一个意外收获。

其实从这里，你可以发现希尔伯特第六问并不是一个可以被一劳永逸解决的问题，而是一个“程序性的召唤”，有趣的是，我们都知道数学是以问题导向为基础的，来自于这一个问题的导向为数学带来了很多，物理学公理化的失败或重构，反哺了整个20世纪的纯粹数学，催生了测度论、泛函分析、广义函数论、算子代数体系以及现代微分几何，这条从物理直觉走向数学严密的无尽之路，本身构成了现代数学物理的一部协同进化的历史。

在现代

因此，我想在这里指出，邓煜解决了希尔伯特第六问题这一说法完全不准确，因为这并不是一个可以被解决的问题，每次当物理走向新的领域，第六问题就永远在膨胀，这是一个open problem。

我们知道希尔伯特第六问题想要解决的是从微观推导出宏观，这其中又分为两个部分，从微观到介观、从介观到宏观。

邓煜解决的部分是从大量微观粒子的牛顿动力学系统推导出介观的统计分布方程即玻尔兹曼方程在时间尺度上的核心瓶颈，这个是从微观到介观的部分，在严格的数学意义上彻底打通了从牛顿力学推导流体方程的通道。

从历史到尝试

玻尔兹曼方程

要理解邓煜的核心工作——从牛顿动力学系统推导出玻尔兹曼方程，我们得先知道玻尔兹曼方程是什么，在传统意义下如何推导。

考虑一个由大量全同单原子分子组成的气体系统。定义单体分布函数 f(\mathbf{r}, \mathbf{v}, t)。在时刻 t，处于位置 \mathbf{r} 到 \mathbf{r} + \mathrm d\mathbf{r}，速度在 \mathbf{v} 到 \mathbf{v} + \mathrm d\mathbf{v} 范围内的相空间体积元 \mathrm d^3\mathbf{r} \mathrm d^3\mathbf{v} 中的粒子数为： \mathrm dN = f(\mathbf{r}, \mathbf{v}, t) \mathrm d^3\mathbf{r} \mathrm d^3\mathbf{v} 如果没有粒子间的碰撞，相空间中的粒子仅在外部保守力场 \mathbf{F}(\mathbf{r}, t) 的作用下运动。根据哈密顿力学刘维尔定理，相空间流体是不可压缩的。在时间 \mathrm dt 内，位于 (\mathbf{r}, \mathbf{v}) 的粒子将运动到 (\mathbf{r} + \mathbf{v}\mathrm dt, \mathbf{v} + \dfrac{\mathbf{F}}{m}\mathrm dt)。基于粒子数守恒： f\left(\mathbf{r} + \mathbf{v}\mathrm dt, \mathbf{v} + \frac{\mathbf{F}}{m}\mathrm dt, t +\mathrm dt\right)\mathrm d^3\mathbf{r}'\mathrm d^3\mathbf{v}' = f(\mathbf{r}, \mathbf{v}, t)\mathrm d^3\mathbf{r}\mathrm d^3\mathbf{v} 相体积不变，因此 \mathrm d^3\mathbf{r}'\mathrm d^3\mathbf{v}' = \mathrm d^3\mathbf{r}\mathrm d^3\mathbf{v} ，将左式进行泰勒展开并取一阶项，可得全导数为零： \frac{\mathrm df}{\mathrm dt} = \frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{r}} f + \frac{\mathbf{F}}{m} \cdot \nabla_{\mathbf{v}} f = 0 在实际气体中，粒子间存在相互作用（碰撞），导致粒子会跳入或跳出我们所关注的相空间体积元 \mathrm d^3\mathbf{r} \mathrm d^3\mathbf{v}。因此，分布函数的全导数等于碰撞导致的净变化率：

\frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{r}} f + \frac{\mathbf{F}}{m} \cdot \nabla_{\mathbf{v}} f = \left(\frac{\partial f}{\partial t}\right)_{\text{coll}} 为了推导右侧的碰撞积分，需要引入以下假设，即稀薄气体极限：

1.只考虑二体碰撞

2.碰撞是局域的且瞬时的

3.分子混沌假设，在碰撞前，参与碰撞的两个粒子的速度是不相关的，即二体分布函数可以分解为单体分布函数的乘积：f^{(2)}(\mathbf{r}, \mathbf{v}, \mathbf{v}_1, t) = f(\mathbf{r}, \mathbf{v}, t)f(\mathbf{r}, \mathbf{v}_1, t)。

注意这里，分子混沌假设从牛顿力学的严格角度来看，这个等式是错误的，当两个粒子发生碰撞后，它们的速度和位置由于牛顿相互作用已经产生了强烈的关联，牛顿力学是时间可逆的，因此牛顿力学要求这种关联在未来的演化中被精确记忆，但分子混沌假设强行抹除了这种关联，假设粒子在参加下一次碰撞前，又忘记了彼此的历史，变成了统计独立的，正是这种人为的信息截断和概率化处理，打破了牛顿系统的时间可逆性。

考虑两个质量为 m 的粒子发生弹性碰撞，碰撞前速度为 (\mathbf{v}, \mathbf{v}_1)，碰撞后速度为 (\mathbf{v}', \mathbf{v}'_1)。 由动量守恒和能量守恒： \mathbf{v} + \mathbf{v}_1 = \mathbf{v}' + \mathbf{v}'_1\mathbf{v}^2 + \mathbf{v}_1^2 = (\mathbf{v}')^2 + (\mathbf{v}'_1)^2 定义相对速度 \mathbf{g} = \mathbf{v} - \mathbf{v}_1 及碰撞后的相对速度 \mathbf{g}' = \mathbf{v}' - \mathbf{v}'_1。能量守恒表明 |\mathbf{g}| = |\mathbf{g}'|。

考虑目标粒子具有速度 \mathbf{v}，它与速度在 \mathrm d^3\mathbf{v}_1 范围内的背景粒子发生碰撞。设微分散射截面为 \sigma(|\mathbf{g}|, \Omega)，其中 \Omega 是散射立体角。 在单位时间、单位空间体积内，发生这种特定碰撞的次数为： |\mathbf{v} - \mathbf{v}_1| \sigma(|\mathbf{g}|, \Omega) \mathrm d\Omega f^{(2)}(\mathbf{r}, \mathbf{v}, \mathbf{v}_1, t) \mathrm d^3\mathbf{v}\mathrm d^3\mathbf{v}_1 应用分子混沌假设，体积元内因碰撞导致粒子跳出的速率积分： \Gamma_{\text{loss}} = \mathrm d^3\mathbf{r} \mathrm d^3\mathbf{v} \int_{\mathbb{R}^3} \mathrm d^3\mathbf{v}_1 \int_{S^2} \mathrm d\Omega |\mathbf{v} - \mathbf{v}_1| \sigma(|\mathbf{g}|, \Omega) f(\mathbf{r}, \mathbf{v}, t) f(\mathbf{r}, \mathbf{v}_1, t) 同理，考虑逆碰撞，即碰撞前速度为 (\mathbf{v}', \mathbf{v}'_1) 的粒子发生碰撞，碰撞后速度变为 (\mathbf{v}, \mathbf{v}_1)，从而跳入我们的目标体积元。 逆碰撞的速率正比于 f(\mathbf{r}, \mathbf{v}', t) f(\mathbf{r}, \mathbf{v}'_1, t)。逆碰撞的截面与正碰撞相等，且相空间体积元满足关系： \mathrm d^3\mathbf{v} \mathrm d^3\mathbf{v}_1 = \mathrm d^3\mathbf{v}' \mathrm d^3\mathbf{v}'_1 因此，增益项可以写为： \Gamma_{\text{gain}} = \mathrm d^3\mathbf{r} \mathrm d^3\mathbf{v} \int_{\mathbb{R}^3} \mathrm d^3\mathbf{v}_1 \int_{S^2} \mathrm d\Omega |\mathbf{v}' - \mathbf{v}'_1| \sigma(|\mathbf{g}'|, \Omega) f(\mathbf{r}, \mathbf{v}', t) f(\mathbf{r}, \mathbf{v}'_1, t) 由于 |\mathbf{v}' - \mathbf{v}'_1| = |\mathbf{v} - \mathbf{v}_1| ，将增益项减去损失项并约去体积元即得到碰撞算子。引入简写记号 f \equiv f(\mathbf{r}, \mathbf{v}, t)， f_1 \equiv f(\mathbf{r}, \mathbf{v}_1, t)，f' \equiv f(\mathbf{r}, \mathbf{v}', t)，f'_1 \equiv f(\mathbf{r}, \mathbf{v}'_1, t)，将上述两项结合，得到最终的玻尔兹曼方程： \frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{r}} f + \frac{\mathbf{F}}{m} \cdot \nabla_{\mathbf{v}} f = \int_{\mathbb{R}^3} \mathrm d^3\mathbf{v}_1 \int_{S^2} \mathrm d\Omega |\mathbf{v} - \mathbf{v}_1| \frac{\mathrm d\sigma}{\mathrm d\Omega} \left( f' f'_1 - f f_1 \right) 至此完成。

Lanford的野望

（Lanford的论文是打字机打的，极度野蛮）

Landford首先尝试完成这个艰巨的任务，他如此考虑：

考虑一个包含 n 个粒子的系统，它们位于体积为 \Lambda 的容器中。对于具有光滑势能或硬球相互作用的粒子系统，系统初始的统计状态（概率测度）可以通过一个关于所有粒子坐标和动量对称的概率密度分布 \mu(x_1, \dots, x_n) 来描述，其中相空间坐标 x_i = (\vec{q}_i, \vec{p}_i)。随着时间 t 演化，分布函数 \mu_t 遵循刘维尔方程： \frac{\partial}{\partial t}\mu_t = H_n\mu_t 其中 H_n 是 n-粒子刘维尔算子。对于光滑势能，它是一个一阶线性微分算子： H_n = -\sum_{i=1}^n \frac{\vec{p}_i}{m} \frac{\partial}{\partial \vec{q}_i} + \sum_{i \neq j} \frac{\partial \Phi}{\partial \vec{q}_i}(\vec{q}_i - \vec{q}_j) \frac{\partial}{\partial \vec{p}_i} 对于硬球模型（粒子直径为 d），尽管演化在碰撞瞬间是不连续的，但通过在相空间中识别“入射”与“出射”碰撞构型，并在剔除零测集（如多体碰撞、擦边碰撞）后，系统依然存在光滑的流，且满足刘维尔方程。

\mu_t 包含了 6n 个变量，信息量过于庞大。为了建立宏观与微观的联系，我们将其对 n-j 个粒子的自由度进行积分，定义 j-粒子相关函数 \rho_j(t; x_1, \dots, x_j)： \rho_j(t; x_1, \dots, x_j) = n(n-1)\dots(n-j+1) \int \mathrm dx_{j+1} \dots \mathrm dx_n \mu_t(x_1, \dots, x_n) \rho_1(t; x_1) 描述了在相空间某处找到单个粒子的期望密度，而 \rho_2(t; x_1, x_2) 描述了同时在两个位置找到两个粒子的期望密度，反映了粒子间的关联。

通过对刘维尔方程两侧进行积分，并利用分布函数的对称性，微观演化被转化为一组耦合的无穷阶方程组，即BBGKY层次方程： \frac{\partial}{\partial t}\rho_j(t) = H_j\rho_j(t) + C_{j,j+1}\rho_{j+1}(t) \quad (j=1, 2, \dots, n) 这里的核心在于方程的非闭合性：演化 \rho_j 需要知道 \rho_{j+1} 的信息。H_j 是 j-粒子刘维尔算子，C_{j,j+1} 是碰撞算子，描述了额外的第 j+1 个粒子与前 j 个粒子发生碰撞而产生的改变。对于硬球系统，Lanford 给出了碰撞算子 C_{j,j+1} 的积分形式： (C_{j,j+1}\rho_{j+1})(x_1, \dots, x_j) = d^2 \sum_{i=1}^j \int \mathrm d\vec{p}_{j+1} \mathrm d\hat{\omega} \, \hat{\omega} \cdot \left(\frac{\vec{p}_{j+1} - \vec{p}_i}{m}\right) \rho_{j+1}(x_1, \dots, x_j, \vec{q}_i + \mathrm d\cdot\hat{\omega}, \vec{p}_{j+1}) 其中，\hat{\omega} 是单位球面上的立体角积分变量，代表碰撞时的相对方位；积分前的 d^2 来源于硬球碰撞的截面积参数。此即BBGKY方程组

为了推导玻尔兹曼方程，我们考虑的是一种空间不均匀的稀薄气体。Lanford采用了Grad提出的特殊极限（即Boltzmann-Grad极限）：

1.让粒子数 n \to \infty 且硬球直径 d \to 0。

2.在取极限的过程中，严格保持 nd^2 为一个有限的非零常数。

粒子的平均自由程 \lambda 大致可以估算为 \lambda \approx \dfrac{V(\Lambda)}{n\pi d^2}。保持 nd^2 固定，等价于保持宏观平均自由程固定。与此同时，气体的体积分数（正比于 nd^3）将趋于 0。这种极限数学上极好，它让粒子在空间中变为几何上的“点”（即排除了多体同时碰撞和有限体积效应），但同时保留了它们之间发生二体碰撞的有限概率。

随着 n \to \infty，之前定义的微观相关函数 \rho_j^{(d)} 的空间积分 \int \rho_j^{(d)} \mathrm dx_1 \dots \mathrm dx_j = n(n-1)\dots(n-j+1) 会发散趋于无穷大。为了在极限下得到具有良好数学行为的物理量，Lanford 引入了重标度相关函数 f_j^{(d)}： f_j^{(d)}(x_1, \dots, x_j) = n^{-j}\rho_j^{(d)}(x_1, \dots, x_j) 通过除以 n^j，我们保证了在 d \to 0 时，其积分的极限具有良好定义，即 \lim_{d \to 0} \int f_j^{(d)} \mathrm dx_1 \dots \mathrm dx_j = 1 将这个变量替换代入第一部分中的 BBGKY 层次方程，会发生一个关键的代数变化。原来的碰撞算子 C_{j,j+1} 包含参数 d^2： \frac{\partial}{\partial t}\rho_j(t) = H_j\rho_j(t) + C_{j,j+1}\rho_{j+1}(t) 代入 f_j 和 f_{j+1} 后，方程变为： \frac{\partial}{\partial t} f_j^{(d)}(t) = H_j f_j^{(d)}(t) + (n d^2) \tilde{C}_{j,j+1} f_{j+1}^{(d)}(t) 此时，方程中原本的截面因子 d^2 自然地与系统粒子数 n 结合，变成了有限常数 nd^2。

现在，我们对重标度的BBGKY方程形式上取 d \to 0 的极限，并假设 f_j^{(d)} 在极限下收敛于某函数 f_j^{(0)}。观察方程的各项如何演化：

当硬球直径 d \to 0 时，原本 j 个粒子系统内部发生的碰撞概率趋于零（因为它们占用的截面在缩小，且没有任何宏观尺度 n 的放大来补偿内部这有限的 j 个粒子的碰撞率）。因此，刘维尔算子 H_j 退化为一个纯粹的、无碰撞的自由流动算子 H_j^{(0)}： H_j^{(0)} = -\sum_{i=1}^j \frac{\vec{p}_i}{m} \cdot \frac{\partial}{\partial \vec{q}_i} 在微观硬球的碰撞积分中，第 j+1 个粒子与第 i 个粒子碰撞时，其空间坐标是 \vec{q}_{j+1} = \vec{q}_i + d\cdot\hat{\omega}。当 d \to 0 时，这种空间上的细微非定域性消失，碰撞变成完全定域的（\vec{q}_{j+1} \to \vec{q}_i）。我们将积分域分为两部分（区分碰撞前和碰撞后，利用相空间测度不变性），可以得到极限界的碰撞算子 C_{j,j+1}^{(0)} 的显式表达： (C_{j,j+1}^{(0)} f_{j+1})(x_1, \dots, x_j) = nd^2 \sum_{i=1}^j \int_{\hat{\omega} \cdot (\vec{p}_i - \vec{p}_{j+1}) \ge 0} \mathrm d\hat{\omega} \mathrm d\vec{p}_{j+1} \, \hat{\omega} \cdot \left(\frac{\vec{p}_i - \vec{p}_{j+1}}{m}\right)\times \{f_{j+1}(x_1, \dots, \vec{q}_i, \vec{p}_i', \dots, \vec{q}_i, \vec{p}_{j+1}') - f_{j+1}(x_1, \dots, \vec{q}_i, \vec{p}_i, \dots, \vec{q}_i, \vec{p}_{j+1})\} 这里，(\vec{p}_i', \vec{p}_{j+1}') 是使得出射动量为 (\vec{p}_i, \vec{p}_{j+1}) 时的入射动量，该形式反映了我们在前面推导过的碰撞增益和碰撞损失。

由此，我们得到了极限下描述系统的方程组——玻尔兹曼层次方程 \frac{\partial}{\partial t}f_j(t) = H_j^{(0)}f_j(t) + C_{j,j+1}^{(0)}f_{j+1}(t) 这一步的深刻之处在于，玻尔兹曼层次方程完全是一个线性方程组。

在得到线性化的玻尔兹曼层次方程后，直接求解偏微分方程是非常困难的。Lanford采用了理论物理中处理相互作用的标准手段：将非微扰项和微扰项分离，利用杜哈梅原理将微分方程转化为时间积分方程，并通过不断迭代得到一个无穷级数，类似于量子场论中的Dyson级数。

为了让数学表述更简洁，Lanford 将所有的相关函数序列 (f_1, f_2, \dots, f_j, \dots) 视为一个抽象的向量空间中的单一向量 \underline{f}。这样，无论是有限直径 d 的重标度BBGKY层次方程，还是极限下的玻尔兹曼层次方程，都可以抽象地写为同一个演化方程： \frac{\mathrm d}{\mathrm dt}\underline{f}(t) = H\underline{f}(t) + C\underline{f}(t) 其中，H 代表系统没有发生新粒子加入碰撞时的内部流动算子，C 代表发生外部碰撞的算子。

如果我们暂时忽略碰撞算子 C（即假设气体极其稀薄以至于碰撞不再发生），演化方程变为 \dfrac{\mathrm d}{\mathrm dt}\underline{f}(t) = H\underline{f}(t)。这个方程的解可以通过自由流动传播子 S(t) 直接写出，形式为 S(t)\underline{f}(0)。但必须注意微观BBGKY和极限玻尔兹曼在自由传播子上的本质区别：

微观BBGKY传播子 S(t) ：对于真实存在直径 d 的硬球，传播子作用下的解定义为： (S(t)\underline{f})_j(x_1, \dots, x_j) = f_j(T^{-t}(x_1, \dots, x_j)) 这里的 T^{-t} 是沿着真实微观动力学轨迹在时间上倒退 t。这意味着这 j 个粒子在倒退过程中，它们彼此之间仍然可能发生有限体积的弹性碰撞。

极限玻尔兹曼传播子 S^{(0)}(t)：当 d \to 0 时，这 j 个粒子在空间中变为没有体积的点，它们彼此之间发生碰撞的概率严格为 0。因此，H_j^{(0)} 退化为纯粹的直线运动。其自由传播子完全解耦： (S^{(0)}(t)\underline{f})_j(x_1, \dots, x_j) = f_j\left(\vec{q}_1 - t\frac{\vec{p}_1}{m}, \vec{p}_1, \dots, \vec{q}_j - t\frac{\vec{p}_j}{m}, \vec{p}_j\right) 现在我们将被忽略的碰撞算子 C 作为微扰加回来。通过对积分方程的不断迭代，微观BBGKY层次的精确形式解被展开为以下微扰级数： \underline{f}(t) = S(t)\underline{f}(0) + \sum_{m=1}^{\infty} \int_0^t \mathrm dt_1 \int_0^{t_1} \mathrm dt_2 \dots \int_0^{t_{m-1}} \mathrm dt_m S(t-t_1) C S(t_1-t_2) \dots C S(t_m)\underline{f}(0) 同理，极限玻尔兹曼层次的形式解也展开为类似的无穷级数： \underline{f}^{(0)}(t) = S^{(0)}(t)\underline{f}^{(0)}(0) + \sum_{m=1}^{\infty} \int_0^t \mathrm dt_1 \int_0^{t_1} \mathrm dt_2 \dots \int_0^{t_{m-1}} \mathrm dt_m S^{(0)}(t-t_1) C^{(0)} S^{(0)}(t_1-t_2) \dots C^{(0)} S^{(0)}(t_m)\underline{f}^{(0)}(0) 这个式子当然具有很清晰的物理图像，我们要计算 t 时刻 j 个粒子的相关函数，等价于我们把这 j 个粒子按自由轨迹往回倒推。在倒推的过程中：

第0项，粒子在整个 t 时间内完全独立、直线倒退，没有任何新粒子撞上它们。

第1项（m=1），在倒退到 t_1 时刻时，碰巧有第 j+1 个粒子撞上了这 j 个粒子中的某一个（算子 C 的作用）。于是带着这额外的第 j+1 个粒子，一共 j+1 个粒子继续直线倒推回初始时刻。

第 m 项，在历史时间 t_1 > t_2 > \dots > t_m 中，依次凭空创生（倒推视角）了 m 个与现有粒子发生碰撞的外部粒子。最终我们需要计算在初始时刻 t=0 时，这 j+m 个粒子的初始联合概率分布。

对于微观 BBGKY 级数，由于系统总共只有 n 个粒子，当追踪的相关粒子数超过 n 时 f_j \equiv 0。因此，BBGKY 的级数实际上是一个有限和，代数操作是合法的。

但是，对于玻尔兹曼极限下的无穷级数，由于 n \to \infty，我们不仅不能假定 f_j^{(0)} 对于很大的 j 为零（特别是如果我们要求初始状态满足因子化时，f_j^{(0)} 在 j 很大时绝对不为 0），这意味着玻尔兹曼级数是一个真正的无穷级数。我们完全不知道这个级数是否收敛，如果它发散，那么玻尔兹曼层次方程在数学上就不存在这样形式的解。

至此整个问题随之转化为了一个纯粹的泛函分析问题：如何证明这个无穷级数是收敛的？它能收敛多久？

在泛函分析中，要证明这样一个级数有意义，即收敛到一个有限的函数，必须定义一个合适的范数，并证明在这个范数下，级数的每一项都不会无限膨胀。

为了控制相关函数的增长，Lanford巧妙地引入了一个带有物理背景的泛函空间，其权重函数是麦克斯韦动量分布。首先，定义一个温度参数为 \beta 的理想麦克斯韦动量分布 \sigma_\beta(\vec{p})： \sigma_\beta(\vec{p}) = \left(\frac{\beta}{2\pi m}\right)^{3/2} \exp\left[-\frac{\beta}{2m}\vec{p}^2\right] 接着，对于定义在 j-粒子相空间上的函数 f_j，Lanford 引入了范数 \|f_j\|_\beta： \|f_j\|_\beta = \sup_{x_1, \dots, x_j} \frac{|f_j(x_1, \dots, x_j)|}{\prod_{i=1}^j \sigma_\beta(\vec{p}_i)} 物理上，这要求系统的 j-粒子关联函数在动量空间上的衰减速度不能慢于麦克斯韦分布。最后，对于整个无穷序列 \underline{f} = (f_1, f_2, \dots)，为了控制随着粒子数 j 增加带来的组合爆炸，引入了一个类比于统计力学中逸出度的参数 z > 0，定义总范数 \|\underline{f}\|_{z,\beta}： \|\underline{f}\|_{z,\beta} = \sup_j \frac{\|f_j\|_\beta}{z^j} 所有使得该范数有限的序列 \underline{f}，构成了一个巴拿赫空间，记为 Y_{z,\beta}。

我们评估自由流动传播子 S(t) 和碰撞算子 C 在这个空间中的表现，由于自由演化或纯弹性碰撞严格遵守能量守恒，系统总动能 H(x) = \sum \vec{p}_i^2 / 2m 保持不变，而我们的空间权重恰好是 e^{-\beta H(x)} 的形式。因此，演化不改变范数大小：\|S_j(t)f\|_\beta = \|f\|_\beta。这说明纯传播过程是非常可控的。

麻烦在于碰撞算子 C_{j,j+1} ，微观计算中，它涉及到对动量 \vec{p}_{j+1} 乘以相对速度 |\vec{p}_i - \vec{p}_{j+1}| 的积分。由于这个相对速度因子的存在，积分会产生无界的动量增长，导致算子 C_{j,j+1} 无法在同一个空间 Y_\beta 中保持有界。通过略微牺牲一点温度参数，即选择 \beta' < \beta（相当于允许目标空间的动量尾部略宽一点），可以吸收掉相对速度带来的多项式增长。此时，C_{j,j+1} 从 Y_\beta 映射到 Y_{\beta'} 是有界的。同时，为了压制随着粒子数 j 带来的组合增长，还需要调整参数 z' 使得 z' > (\beta/\beta')^{3/2} z。

现在我们来考察级数中第 m 项对应的多重积分： \int_0^t \mathrm dt_1 \int_0^{t_1} \mathrm dt_2 \dots \int_0^{t_{m-1}} \mathrm dt_m S(t-t_1) C S(t_1-t_2) \dots C S(t_m)\underline{f}(0) 被积函数中有 m 个碰撞算子 C 连乘。根据泛函估计，连续作用 m 次算子会产生一个正比于 m! 的巨大系数，因为第 1 次碰撞有 j 个目标，第 2 次有 j+1 个目标，以此类推，产生组合数累乘。如果直接估计，这个级数似乎必然发散。但是，多重积分 \int_0^t \mathrm dt_1 \dots \int_0^{t_{m-1}} \mathrm dt_m 的相空间体积恰好是 \dfrac{t^m}{m!} 。时间积分贡献的 1/m! 抵消了碰撞组合带来的 m! 增长。最终，第 m 项作为一个从 Y_{z,\beta} 到 Y_{z',\beta'} 的算子，其范数的量级被严格界定为： (|t| / t_0)^m 其中，常数 t_0 的表达式为： t_0 = \left[A \pi n d^2 z \left(\frac{m\beta}{3}\right)^{-1/2}\right]^{-1} 这个 t_0 具有极其清晰的物理意义：n z 代表密度，(3/m\beta)^{1/2} 代表热运动均方根速度 V_{rms}，那么 n d^2 z V_{rms} 正是碰撞频率，因此 t_0 的量级正好是系统平均自由时间的一个分数倍。

有了 (|t|/t_0)^m 的估计，这变成了一个极其简单的几何级数。只要 \vert{}t\vert{} < t_0，这个无穷级数就是绝对且一致收敛的。

最关键的是，Lanford证明了这种收敛性仅仅依赖于时间 t 与平均自由时间 t_0 的比值，而严格一致地独立于粒子数 n、直径 d 等参数（只要 nd^2 固定）。正是由于拥有了这种关于 d 的一致收敛性，我们在泛函空间中取极限操作 \lim_{d \to 0} 才是数学上合法的。于是，我们可以将微观BBGKY的无穷级数逐项取极限，证明它在 d \to 0 时精确地收敛于玻尔兹曼层次方程的无穷级数解。

我们最终的目标是那个著名的非线性玻尔兹曼偏微分方程，在 t=0 时刻，我们假设宏观气体的制备是“纯粹”的，即微观粒子之间没有预先存在的关联。这在数学上体现为相关函数的因子化。对于极限下的初始相关函数 f_j^{(0)}，因子化意味着任意 j 个粒子的联合分布，等于单个粒子分布的简单乘积： f_j^{(0)}(x_1, \dots, x_j) = \prod_{i=1}^j f_1^{(0)}(x_i) 物理上，Lanford 指出，这等价于大数定律的表现：对于空间中任意宏观区域 \Delta，粒子数比例 n_\Delta/n 依概率收敛于单粒子分布的期望值 \int_\Delta f_1^{(0)}(x_1)\mathrm dx_1这意味着在 t=0 时刻，粒子是完全统计独立的。

如果系统初始是没有关联的，那么在经历了一段时间的相互碰撞后，粒子之间还会保持独立吗？在微观层面，有限的直径 d ，碰撞不可避免地会产生关联，碰过的粒子之间显然有了动量纠缠。但Lanford证明了，在Boltzmann-Grad极限（d \to 0）下，这种因碰撞产生的关联在宏观单粒子分布的尺度上被稀释到了零。由于我们在前面已经证明了玻尔兹曼层次方程在泛函空间 Y_{z,\beta} 中的解是由微扰级数唯一确定的，我们只需做一个代数验证： 假设单粒子分布 f_1(t, x) 是玻尔兹曼方程的解，我们构造一个乘积形式的序列 f_j(t, x_1, \dots, x_j) = \prod_{i=1}^j f_1(t, x_i)。可以严格验证，这个因子化的序列形式上直接满足线性的玻尔兹曼层次方程。基于级数解在 0 < t < t_0 范围内的唯一性，由于真实的极限演化序列与我们构造的因子化序列具有相同的初始条件，因此它们必然在整个 0 < t < t_0 期间完全重合。如果在Boltzmann-Grad极限下初始满足分子混沌，那么在级数收敛的寿命 t_0 内，分子混沌假设严格成立，因子化性质被动力学严格保持。

一旦证明了 0 < t < t_0 时刻相关函数保持因子化，推导非线性玻尔兹曼方程就水到渠成了。我们回到极限下的玻尔兹曼层次方程的第一阶 \frac{\partial}{\partial t}f_1(t; x_1) = H_1^{(0)}f_1(t; x_1) + C_{1,2}^{(0)}f_2(t; x_1, x_2) 这里，单粒子演化依赖于双粒子相关函数 f_2(t; x_1, x_2)。现在，我们把刚刚证明的因子化结果 f_2(t; x_1, x_2) = f_1(t; x_1)f_1(t; x_2) 直接代入上述方程中的碰撞算子积分中： \frac{\partial}{\partial t}f_1(t; \vec{q}_1, \vec{p}_1) = -\frac{\vec{p}_1}{m}\cdot \frac{\partial}{\partial \vec{q}_1}f_1(t; \vec{q}_1, \vec{p}_1) + nd^2 \int_{\hat{\omega} \cdot (\vec{p}_1 - \vec{p}_2) \ge 0} d\hat{\omega} d\vec{p}_2 \, \hat{\omega} \cdot \left(\frac{\vec{p}_1 - \vec{p}_2}{m}\right)\times \left\{ f_1(t; \vec{q}_1, \vec{p}_1')f_1(t; \vec{q}_1, \vec{p}_2') - f_1(t; \vec{q}_1, \vec{p}_1)f_1(t; \vec{q}_1, \vec{p}_2) \right\} 这实际上就是 \frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{r}} f + \frac{\mathbf{F}}{m} \cdot \nabla_{\mathbf{v}} f = \int_{\mathbb{R}^3} \mathrm d^3\mathbf{v}_1 \int_{S^2} \mathrm d\Omega |\mathbf{v} - \mathbf{v}_1| \frac{\mathrm d\sigma}{\mathrm d\Omega} \left( f' f'_1 - f f_1 \right) 至此，我们看到了传统的工作如何解决这个问题。

笨蛋！问题是时间！

玻尔兹曼方程最著名的推论是H定理， \dfrac{\mathrm dH}{\mathrm dt} \le 0 ，即系统总是趋向于熵增的平衡态，这打破了时间反演对称性。

也就是说，微观动力学是时间可逆的，而玻尔兹曼方程是时间不可逆的，不可逆性究竟是从哪里钻出来的？Lanford的证明给出了数学答案：不可逆性不来源于人为引入的分子混沌假设或者说因子化，而是来源于 d \to 0 的玻尔兹曼-格拉德极限操作本身。

在证明前面的级数收敛以及构造极限碰撞算子 C_{j,j+1}^{(0)} 时，我们做了一个看似不起眼但打破时间对称性的操作：我们仅仅将即将发生碰撞的构型等效替换成了对应出射动量的函数值，从而消去了微观尺寸 d。在正向时间演化中，这对应于忽略碰撞后粒子分离时的微小空间偏差；但如果我们试图对时间反演（t < 0），同样的极限操作会导致碰撞算子的符号反转。

正如Lanford原文所述：“玻尔兹曼层次方程（与玻尔兹曼方程一样）在时间反演下是不具有不变性的。也就是说，不可逆性出现在取极限 d \to 0 的过程中，而不是出现在极限重标度相关函数因子化的假设中。”

但是，为什么说Lanford的证明完全没有给出一条美好的道路呢？因为从微观到介观，目的地必然是再走向宏观，然而，Lanford 只能证明玻尔兹曼方程在极短的时间（t < t_0）内有效，这个 t_0 仅仅是系统平均自由时间的一个极小分数倍。然而，诸如流体达到热平衡、或者气体表现出宏观Navier-Stokes方程行为等物理现象，都需要经历无数次的粒子碰撞和宏观时间累积。Lanford的“瞬间”有效性并未完整建立微观到宏观的链条。产生这个问题的根源在于Lanford依赖于杜哈梅原理展开的无穷微扰级数，产生了一个幂级数，这就注定了一旦 t 稍大，微扰级数就会在数学上彻底发散失效。在长时间演化中，气体系统不可避免地会发生大量环形碰撞或重复碰撞，比如粒子A撞击B，B撞击C，然后C在经历一系列演化后又撞回A。这种重复碰撞使得粒子在碰撞前就不再是统计独立的，从根本上污染了核心的分子混沌假设。Lanford通过限制时间在极短范围绕开了这个问题，却也使得该方法在长时演化面前无能为力。

一般而言，到了这种危难的时刻，我们的主角就要登场了，接下来有请邓煜登台。

邓煜和其工作

我们的时代是一个急功近利的时代，很多人着急最终问题的解决，但是

“每一新阶段的前进，即进一步的规定，同时也是向内部的退隐，更大的伸展同时也是更高的密集。”
——黑格尔《逻辑学》（下卷）绝对理念

我主张理解一个人的工作需要放在其历史工作底色下来看，自然，邓煜也不例外。

早期的工作和底色

希尔伯特六第六问题所涉及的PDEs，在色散偏微分方程领域表现为波湍流理论(Wave Turbulence Theory，下文简称WTT)，即从弱非线性色散波严格推导出波动力学方程。

他早期一篇相当出色的纯PDEs论文是Global solutions of the gravity-capillary water-wave system in three dimensions^[1]。在此之前，三维空间中纯重力水波（g>0, \sigma=0）或纯毛细水波（g=0, \sigma>0）的全局适定性已被解决，但这篇论文首次解决了重力-毛细水波系统（同时存在重力和表面张力，即 g>0, \sigma>0）在三维空间中小初值下的全局正则性。

在拟线性偏微分方程的全局存在性证明中，经典的Klainerman向量场方法和弱色散衰减是核心工具。但这篇论文面临两个极其致命的数学困难的叠加，这也是以往方法失效的原因。系统的线性色散关系为 \Lambda(\xi) = \sqrt{g|\xi| + \sigma|\xi|^3}。其群速度的导数 \Lambda''(\xi) 在特定频率 |\xi| = \gamma_0 处为零（退化）。这意味着线性解的最佳点态衰减速度仅为 |t|^{-5/6}，严格慢于三维通常的 \vert{}t\vert{}^{-1} 衰减。另一方面，对于 |t|^{-5/6} 的衰减，通常需要利用 Shatah 的正规型方法来克服。然而，该系统的二次非线性项存在一个余维为1的庞大时间共振集，即存在一个完整的共振球面，且缺乏可以抵消这些共振的零结构，导致传统正规型方法完全失效。

对于这个问题本身是如何解决的并不是主题，重点在于说他在里面展现出了什么样的思想和方法，这里面有几个技术我认为值得称道，一是共振流形的精细几何解析，系统同时存在导致点态衰减极慢的球面（半径为 \gamma_0）和引发时空共振的球面（半径为 \gamma_1），且它们彼此分离。面对庞大且无零结构抵消的时间共振集，引入了相函数在时间共振超曲面上的受限非退化条件，从而在 L^2 估计中通过 TT^* 论证获得了傅里叶积分算子的最优界。在WTT中，对“共振流形”的控制从连续几何走向了丢番图逼近与离散几何，水波中对 \nabla \Phi = 0 和 \Phi = 0 的非退化性分析深刻影响了后来WTT中在极大超立方体上对离散共振格点计数的数论边界控制，这是实现无限体积极限的几何基石。二是守恒律下的代数消去与强半线性结构，通过对Dirichlet-Neumann算子进行仿线性化，并构造经过二次修正的复未知量，将原有的拟线性偏微分方程转化为包含特殊二次结构的方程，揭示了与系统哈密顿能量守恒相关的强半线性结构，空间-时间积分乘子的大小与调制之间存在代数相关性，从而在共振区域获得了额外的导数增益，避免了能量估计中的导数损失。在推导波动力学方程时，薛定谔方程的非线性项会产生海量的相位振荡，我是认为这种思路和WTT中对威克配对机制中精确消去现象的捕捉是一脉相承的，通过代数结构的对称性来压制高阶发散是推导碰撞积分的关键。三是敢于放弃边角修补，放弃了传统的基于纯频率或纯空间的范数，构造了高度定制化的原子 Z-范数，该范数在空间和频率上同时进行极其精细的局部化分解（利用 Q_{jk} 算子），并且其结构严重依赖于时空共振集的几何位置与形状。传统的索伯列夫空间无法刻画波的统计独立性和相位随机性，邓煜在WTT中引入了类似于Besov空间与维纳代数结合的定制化随机张量空间，这种基于物理问题内禀几何构造空间的哲学保证了在动力学时间尺度下非线性相互作用不会破坏波包的统计特性。四是能够解决长时动力学展开的关键，本文中因为庞大共振集导致标准正规型方法失效，转而在傅里叶空间中对杜哈梅公式进行探索，发现了迭代共振的缓慢传播现象，即在三次和四次多重线性积分中，如果嵌套的子相位驻点存在，整体相位的频率导数也必然为零，以此将误差项的支撑集严格限制在特定频段。水波中的三次、四次Duhamel迭代，在WTT中爆发成了极其庞大和复杂的费曼图，水波问题中对嵌套驻点的估计技术升级为WTT中对带环树图相消干涉的积分界限控制，这是数学上第一次能够严格控制动力学时间尺度上无穷级数展开收敛性的技术。

在偏微分方程领域，许多学者习惯于套用成熟的框架，邓煜难得可贵的地方就是当标准工具因本质困难失效时，不试图强行修补外在框架，而试着做一些底部为之适配的构造工作，我很拜服这种方法论。另一方面，分析学中的发散或衰减不足，往往是因为放缩过度，邓煜非常擅长在看似杂乱无章的内容物中找到守恒律保证的对称性或代数零点，将硬分析上的一些界限问题，以一种优雅的手段转化为了寻找代数多项式在特定几何面上的代数消去这样的问题。总的来说，邓煜具有一种极其强悍且细腻的系统工程能力。

范式转移

参见：

【“老师，请解算我未来的人生”【中国计算机史填词】】 【精准空降到 02:29】

在早期的工作中，邓煜尚且还停留在一种范式：水波问题中，初始条件是确定性的。为了证明全局正则性，数学家必须针对最坏可能的波包分布进行估计，通过寻找系统深层的代数结构（如前文提过的强半线性结构）来抵消共振。

然而就像我们前面看到的，希尔伯特第六问题的核心是从微观的可逆动力学推导出宏观的不可逆统计方程，这必须依赖于大系统的随机性。这就不得不提理解邓煜后期所有工作的核心方法论文献Random tensors, propagation of randomness, and nonlinear dispersive equations^[2]

在这篇论文中，初始数据变成了由独立高斯随机变量给出的随机数据。邓煜与合作者将研究的核心问题转向了随机性如何在非线性色散流下传播，理解随机性结构（高频和精细尺度上的独立性）如何在非线性演化中得以保持，是解决WTT中随机相位近似能否在长时间内成立的根本前提。这篇论文放弃了传统的确定性函数空间，建立了一套全新的随机张量理论，这其实也是前文所提的方法论的某种呈现，

来仔细说说这篇论文吧，在奇异随机偏微分方程领域，Hairer的正则性结构和Gubinelli等人的仿控制微积分在热传导方程（抛物型）上取得了巨大成功。然而，邓煜等人指出，这些方法在色散方程（如薛定谔方程）面前遭遇了根本性的障碍，原因有三：

1.热传导方程满足最大值原理且热核有指数衰减，极其契合基于物理空间局部泰勒展开的正则性结构理论。但薛定谔方程是非局部的，不具备这些性质。

2.抛物理论严重依赖于 C^\alpha （Hölder）空间。然而薛定谔流在 C^\alpha 空间上是无界的，它只在 L^2 框架的Sobolev空间（H^s）上有界。

3.抛物方程的基本解 (\partial_t - \Delta)^{-1} 能直接提供两个空间导数的平滑提升，而薛定谔方程在空间上没有平滑效应，其平滑性只能通过时空耦合的 X^{s,b} 空间来体现，这与现有的SPDE理论完全不兼容。

因为这种困境，邓煜等放弃了物理空间，转向全局的傅里叶空间，把解展开为带系数的多线性高斯函数，这些系数就是所谓的随机张量。分析的重点也从估计 C^\alpha 范数，变成了估计这些张量的 L^2 \to L^2 算子范数。

为什么随机化能降低方程的正则性要求？有一个极其漂亮的启发式推导，为了简化，我们考察带有无配对非线性项 \mathcal{N}_{np} 的薛定谔方程： (\mathcal{N}_{np}(u))_k := \sum_{k_1 - \dots + k_p = k} u_{k_1}\overline{u_{k_2}}\dots u_{k_p} 并设初值 u(0) \in H^s，我们想要控制第一次非线性迭代 u^{(1)}(t) 在时间 |t|\sim 1 时的 H^s 范数。假设我们有一个集中在频域 |k| \sim N 上的确定性初值： u(0) = N^{-\alpha} \sum_{\vert{}k\vert{} \sim N} e^{ik\cdot x}, \quad \text{where } \alpha = s + \frac{d}{2} 此时 \|u(0)\|_{H^s} \sim 1 ，在傅里叶端，第一次迭代的系数 u_k^{(1)}(t) 大致为： u_k^{(1)}(t) \sim N^{-p\alpha} \sum_{k_j \in \mathbb{Z}^d, |k_j| \sim N} \frac{1}{\langle \Omega \rangle} 其中共振相位 \Omega := |k|^2 - |k_2|^2 + \dots - |k_p|^2。如果我们只限制在完全共振面（\Omega = 0）上，通过维度计数，内部求和的项数（即自由度）大小为 N^{pd-d-2}。因此： \|u^{(1)}(t)\|_{H^s} \sim N^{(p-1)(\frac{d}{2}-s)-2} 为了让迭代不发散，必须满足： s \ge \frac{d}{2} - \frac{2}{p-1} := s_{cr} 这就是经典的确定性临界指数。

现在，引入随机化。初值变为独立同分布的标准高斯随机变量 g_k(\omega) 的叠加： u(0) = N^{-\alpha} \sum_{|k| \sim N} g_k(\omega) e^{ik\cdot x} 此时的傅里叶系数迭代变为： u_k^{(1)}(t) \sim N^{-p\alpha} \sum_{k_j \in \mathbb{Z}^d, |k_j| \sim N} \frac{1}{\langle \Omega \rangle} g_{k_1}\overline{g_{k_2}}\dots g_{k_p} 由于 g_k 之间是相互独立的，中心极限定理告诉我们，大量独立随机变量求和时，其期望大小不是项数本身，而是项数的平方根。因此，同样在 \Omega = 0 面上，内部求和在极大概率下的真实大小不再是 N^{pd-d-2}，而是 (N^{pd-d-2})^{1/2} = N^{(pd-d-2)/2}。我们重新计算 H^s 范数： \|u^{(1)}(t)\|_{H^s} \sim N^{-(p-1)s-1} 正则性要求骤降为： s \ge -\frac{1}{p-1} := s_{pr} 这就是论文要攻克的全亚临界区域（s > s_{pr}）的根本来源。

由于在 t>0 后，解的傅里叶系数不再保持纯粹的独立性，上述平方根抵消的启发式论证会失效。为了严格追踪和控制这些被破坏的独立性，必须把解写成多线性高斯分布的展开形式 \Psi_k = \sum_{k_1,\dots,k_r} h_{kk_1\dots k_r} \prod_{j=1}^r \langle k_j \rangle^{-\alpha} g_{k_j}^{\pm}(\omega) 这里的系数 h_{kk_1\dots k_r} 就是一个 (r, 1) 随机张量，它包含了低频分量的所有信息。随着方程的迭代，这些张量会经历两种核心的代数操作：

1.合并，当两个多线性表达式相乘时，例如非线性项 \vert{}u\vert{}^{p-1}u 会导致相乘，对应的张量会发生半张量积式的收缩。文内给出了一个直观的模型：假设有两个表达式相乘 \Psi \cdot \Psi'。 \Psi_k = \sum_{a,b,c,d} h_{kabcd} (\text{高斯项}_1) 和 \Psi'_{k'} = \sum_{u,v,w} h'_{k'uvw} (\text{高斯项}_2) 如果在相乘时发生了频率配对，那么合并后的新张量 H 会对这些配对指标进行求和收缩： H_{mcdv} \sim \sum_{k+k'=m} \sum_{a,b} \langle a \rangle^{-2\alpha} \langle b \rangle^{-2\alpha} h_{kabcd} h'_{k'avb} 每次合并都会使张量的阶数和结构变得极度复杂，在后文必须证明，张量的某种算子范数在经过上述合并操作后，依然能保持有界。

2.修剪，在证明大偏差估计，保证平方根抵消成立时，张量 h 自身必须与它所乘的高斯变量 g_{k_j} 保持统计独立。但是在多次合并后，张量内部可能混入了与 g_{k_j} 相同频段的变量，破坏了独立性。为此，文内引入了修剪操作，在指定的截断频率 R 下，强制将张量 h 与自由高斯变量进行收缩求和，把高频部分消除： h'_{kk_{A'}} = \sum_{k_{A \setminus A'}} h_{kk_A} \prod_{j \in A \setminus A'} \langle k_j \rangle^{-\alpha} g_{k_j}^{\pm} 这就生成了一个低阶的新张量 h'，它严格依赖于低频，从而恢复了与高频高斯变量的独立性。而文中的下降法（Method of descent，Prop 2.8）证明了：修剪后的张量范数，被极大概率地控制在修剪前张量的各个分割范数的最大值之内。

为了实现宏大的随机张量拟设，邓煜必须发明一套极度精密的记号系统来追踪无穷无尽的迭代，随着方程的不断迭代，高斯项和余项会发生巨量的交叉相乘与配对，如果仅用图论中的树来记录，信息的层级会迅速丢失。为此，文章在第3.2节创造了“植物 \mathcal{S}”结构，它由三个不相交的有限集合组成 \mathcal{S} = (\mathcal{L}, \mathcal{V}, \mathcal{Y})：

1.树叶（Leaves, \mathcal{L}）：代表拟设中没有被收缩的高频线性高斯项，即方程的初值部分 \Delta_{N_l}f。每个叶子 l \in \mathcal{L} 带有频率 N_l 和符号 \zeta_l \in \{\pm\}。

2.花朵（ \mathcal{V}）：由于方程有余项 z_{N_f}，具有稍高的正则性但频率较低，花朵代表迭代树中碰到了余项的节点。

（请听：）

3.记忆/过去（ \mathcal{Y} ）：这是我觉得最巧妙的设计，当两个张量合并时，原本的树结构会坍缩，在合并时向集合 \mathcal{Y} 中塞入一个记忆点 \mathfrak{p}，其频率 N_{\mathfrak{p}} 设为参与合并的各分支频率的次大值，这样，极其复杂的树结构被压平成了一个集合，大大降低了代数复杂性，但衰减信息被完美保留。

这样我们可以回顾之前的那两个操作实际上干了什么

1.修剪操作 \mathcal{S}' = \text{Trim}(\mathcal{S}, R)：给定截断频率 R，把所有频率 N_l < R 的叶子从 \mathcal{L} 中砍掉。在张量层面，就是把张量 h 沿着这些低频叶子对应的随机高斯变量 g_{k_l} 缩并（求和）掉，生成一个阶数更低的新张量 h'。

2.合并操作 H = \text{Merge}(h^{(1)}, \dots, h^{(r)}, h, \mathscr{P}, \mathscr{O})：这是非线性项相互作用的结果。将多个低阶张量 h^{(j)} 通过底张量 h 连接起来。这里最难对付的是过度配对（即同一个频率 k 被强制配对了三次及以上）。合并定义通过集合划分 \mathscr{O} 严格枚举了所有可能的配对拓扑。

在合并张量时，底张量 h 是由非线性薛定谔方程的相位结构决定的。我们需要控制这个底张量作为算子的 L^2 \to L^2 范数，这直接转化为在共振流形上的整数格点计数问题。

具体来说，薛定谔方程的非线性项包含两个核心守恒律：

1.动量守恒：k = \sum_{j=1}^p \zeta_j k_j。

2.能量（完全共振）限制：|k|^2 - \sum_{j=1}^p \zeta_j |k_j|^2 = \Gamma。

假设我们要估算一个由上述条件支撑的张量 h_{k k_1\dots k_p}。给定一个指标的分割（Partition） P_0 和 Q_0，我们需要求算子范数 \|h\|_{k k_{P_0} \to k_{Q_0}}。

Schur's Lemma在这里派上了用场，要控制矩阵/张量算子范数，只需控制它的行和与列和的上确界： \|h\|_{k k_{P_0} \to k_{Q_0}} \le \left( \sup_{k_{Q_0}} \sum_{k, k_{P_0}} 1 \right)^{1/2} \left( \sup_{k, k_{P_0}} \sum_{k_{Q_0}} 1 \right)^{1/2} 这就意味着算子范数估计被完全等价地转化为了纯粹的组合计数问题，我们只需要计算：给定 k_{Q_0} 的取值时，有多少组 (k, k_{P_0}) 能满足动量和能量守恒的二次丢番图方程。

然后可以归纳法， p=3 时，方程化为求解 \zeta_1 |k_1|^2 + \zeta_2 |k_2|^2 + \zeta_3 |k_3|^2 = \Gamma。代入动量守恒消去一元，本质上是一个二次超曲面上的整数点计数。利用除数除数函数边界，证明了无配对情况下的解的数量上限为 (M_1 M_2)^{d-1+\theta}。

对一般奇数 p 的归纳是这样的，通过提取两个变量 (k_n, k_{n'})，利用动量守恒确定局部活动范围，将 p 个变量的问题递归降解为 p-2 个变量的问题。插值计算得出边界：\#S \le \prod_{j \in A} (M_j)^{2\alpha_0 + \theta}，其中 \alpha_0 = \dfrac{d}{2} - \dfrac{1}{p-1}。

最终，Prop 4.9整合了上述结果，证明了对于非过度配对的张量，其范数被控制在： \|h\|_{k k_{P_0} \to k_{Q_0}} \le \prod_{j=2}^p (M_j)^{\alpha_0 + \theta} \prod (M_{\text{局部}})^{\text{微小修正}} 然后来到下降法，我们知道修剪操作本质上是将一个高阶的确定性张量 h_{bck_A} 与一组独立的高斯自由变量 \eta_{k_j} 进行收缩求和，产生一个降阶的随机张量 H_{bc}： H_{bc} = \sum_{k_A} h_{bck_A} \prod_{j \in A} \eta_{k_j}^{\zeta_j} 我们需要极大概率下控制随机张量 H 的算子范数 \|H\|_{b \to c}。如果只用普通的大偏差估计计算期望 \mathbb{E}\|H\|^2，会发现各种交叉项的爆炸导致边界极度松散，无法提取到矩阵结构的对消性。

大概是借鉴了Bourgain的 gg^* 参数论证，引入了随机矩阵理论中的高阶迹方法，要控制矩阵最大奇异值，去计算巨大整数次幂 m 下的核 (HH^*)^m： \|H_{bc}\|_{b \to c}^{4m} = \|(HH^*)^m\|_{OP}^2 \le \sum_{b,b'} |(\mathcal{R}_{2m})_{bb'}|^2 其中 \mathcal{R}_n 是迭代相乘 n 次后的核。利用合并定理中的代数展开，\mathcal{R}_n 会产生多达 n 个独立副本的高斯变量 \eta。 (\mathcal{R}_n)_{bc} = \sum_{b'} (\mathcal{R}_{n-1})_{bb'} H_{b'c} = \sum_{b'} \sum_{k_Z, k_A} y_{bb'k_Z} h_{b'ck_A} \prod_{j \in Z} \eta_{k_j}^{\zeta_j} \prod_{j \in A} \eta_{k_j}^{\zeta_j} 因为 \eta 是均值为0的独立随机变量（确切地说是单位圆上的均匀分布），对这个巨大的表达式取极大概率界限时（或求期望），只有那些 \eta 正负配对抵消为 1 的项（即 1_{k_j = k_{j'}}）能够留下来。

所有存留下来的配对项实际上就是把原来的张量 h 从不同的维度进行扭结和自我收缩。文中指出，这种扭结恰好对应于将张量的输入集 A 划分为各种可能的二分图 (B, C)，并求普通半张量积。绝对上界可以严格地被限制为： \|y\| \le \left( \sup_{(B,C)} \|h\|_{bk_B \to ck_C} \right)^{n-1} \|h\|_{bck_A} 将上述界限代回 2m 阶矩不等式，得到： \|H_{bc}\|_{b \to c} \le (\tau^{-1}M)^\theta \left( \sup_{(B,C)} \|h\|_{bk_B \to ck_C} \right)^{1 - \frac{1}{2m}} \|h\|_{bck_A}^{\frac{1}{2m}} 让 m 趋于无穷大，末尾项会直接变为 1 ，可以得到 \|H_{bc}\|_{b \to c} \le \text{小对数损失} \cdot \max_{(B,C)} \|h\|_{bk_B \to ck_C} 这意味着将一个张量与随机高斯场收缩其产生的低阶随机张量的算子范数被极大概率地控制在原高阶确定性张量所有可能的切分范数的最大值之内。

到此为止，我们拥有了足够的手段和工具去处理真正的薛定谔方程，当然，直接处理原始的Wick排序薛定谔方程是灾难性的，去掉了无穷大的期望，残余的质量项在共振面上依然会引发相位偏移，因此需要先化简。

文中引入了一个依赖于时间 t 的相位旋转 \tilde{u_N}(t) = u_N(t) \cdot \exp\left(\frac{(p+1)i}{2} \int_0^t \int_{\mathbb{T}^d} W_N^{p-1}(u_N) \mathrm dt'\right) Wick 排序后的非线性项 W_N^p(u) 包含大量的形如 |u|^2 u 的共振项。通过上述规范变换，方程的非线性项被重整为一阶层级更低的非线性项组合： \sum_{3 \le q \le p} a_{pq} (m_N - \sigma_N)^{(p-q)/2} \mathcal{N}_q(\tilde{u_N}) 其中 m_N 是系统的守恒质量（L^2 范数），\mathcal{N}_q 是纯粹的简单多线性算子，它去除了那些会引发共振相移的恶劣配对。

利用大偏差估计，随机初始质量偏差 m_N^* := m_N - \sigma_N 会以极大概率收敛于一个确定的常数。因此，直接切除极小概率的异常集，将 m_N^* 视为有界的确定性参数。提取线性薛定谔演化流，令 v_N(t) := e^{it\Delta} \tilde{u_N}(t)。此时拉普拉斯算子 \Delta 从方程左侧消失，转移到了非线性项的傅里叶相位 \Omega 中： \Omega = |k|^2 - |k_1|^2 + \dots - |k_q|^2 方程彻底变成了一个积分方程。

为了使用能够捕捉平滑效应的全局时空 X^{s,b} 空间，文中强行引入了一个光滑时间截断函数 \chi(t) 和截断Duhamel算子 \mathcal{I}_\chi。我们将追踪相邻频率截断解的差值 y_N := v_N^\dagger - v_{N/2}^\dagger，只要控制住每一个 y_N，全解 v^\dagger = \sum y_N 就能求和收敛。

然后开始为 y_M 构造张量展开式，我们面对的是一个高阶的多线性方程，最难处理的交互是高低频交互，即一个高频输入的线性演化与多个低频解（如 v_{M^{[\delta]}}^\dagger）相乘。这部分不具备足够的平滑性，无法归入余项，直接展开又会导致无穷无尽的树结构。

定义一个纯粹的随机线性算子 \mathscr{Z}^M： (\mathscr{Z}^M w)_k(t) = -i \sum a_{pq} (m_M^*)^{\dots} \mathcal{I}_\chi \Pi_M \sum_{sym} \mathcal{M}_q(w, v_{M^{[\delta]}}^\dagger, \dots, v_{M^{[\delta]}}^\dagger)_k(t) 算子 \mathscr{Z}^M 的作用就是把高频 w 映射为下一次高低频交互的结果，方程被改写为： y_M = \text{初值} + \mathscr{Z}^M y_M + \text{其他非线性交互} 把 \mathscr{Z}^M y_M 移到左边并求逆，得到算子 \mathscr{R}^M = (1 - \mathscr{Z}^M)^{-1} 。

然后关键在于我们不显式地对这个高低频交互展开成无穷级数，而是利用前面定义的张量运算，证明 \mathscr{Z}^M 在 X^{s,b} 空间中的算子范数严格小于1，从而使得 \mathscr{R}^M 成为一个良定义的有界算子。现在，拟设被正式写出： y_M = \sum_{n \in \{0,1\}} \sum_{\mathcal{S}} \Psi_k[\mathcal{S}, h^{(\mathcal{S},n)}] + z_M 解被剖解为三部分：确定性基础树 h^{(\mathcal{S},0)}，随机迭代树 h^{(\mathcal{S},1)}，和高正则性余项 z_M。在 \mathscr{R}^M 的作用下，作者利用前文的Merge,Trim函数给出了这些张量的严格递归定义：

1.确定性张量 h^{(\mathcal{S},0)} ，它是直接来自初始数据高斯分布的纯线性或完全共振非线性演化，它是确定性的，只针对无花朵（\mathcal{V}=\emptyset）的“素植物”存在。通过积分算子 \mathcal{I}_\chi 直接作用在上一步合并出的底张量 H 上。

2.随机张量 h^{(\mathcal{S},1)} ，它是高度随机的，包含了从 \mathscr{Y}^{M} = \mathscr{R}^M - 1 算子中衍生出的高低频交互项。上一层的张量在输入本次Merge之前，被强制执行了 M^\delta 级别的Trim， \mathcal{S} = \text{Trim}(\text{Merge}(\text{Trim}(\mathcal{S}_1, M^\delta), \dots), M^\delta) 每次非线性交互后，如果不切掉低频部分，新生成的随机张量就会与环境中的低频高斯白噪声耦合。通过 \text{Trim} 强制截断，确保了张量 h^{(\mathcal{S},1)} 严格可测于 \mathcal{B}_{M^{[\delta]}}，从而保证了其与高频高斯场的统计独立，这就为之后运用下降法大偏差估计扫清了根本障碍。

3.余项 z_M ，在迭代中，总有些项不符合高频独立的要求，或者树已经深到不需要进一步展开了。这些项被扔进了 z_M ：

3.1.如果相互作用中包含太多低频余项 z_{N_j}

3.2.如果相互作用导致生成的植物体积 |\mathcal{S}| > D ，z_M 被证明是某个从 X^{b_0} \to X^{b_0} 压缩映射的唯一不动点，因而天然具备足够的正则性衰减 M^{-D_1}。

在定义了上述拟设后，提出了撑起全篇定理的拱顶石，命题 5.1。它宣布：只要存在一个宏大的归纳假设集Local(M)，这个解构体系就能无限进行下去。Local(M)在这里有五个边界

5.32.控制 h^{(\mathcal{S},0)} 在特定共振面 \Gamma 上的加权 L^2 算子范数。

5.34-46.控制 h^{(\mathcal{S},1)} 的双重时空（主调制 \lambda 与花朵调制 \lambda_{\mathcal{V}}）X_{\mathcal{V}}^{b, -b_0} 范数，以及极难处理的傅里叶局部化加权界限。

5.38.将张量塞回高斯变量中生成的解析函数 \Psi^{(\mathcal{S},n)} 必须严格属于次临界的 X^{s,b_0} 空间，且拥有随尺度 N 衰减的特性。

证明 z_M 在 X^{b_0} 具有强衰减 M^{-D_1}。

5.40.线性映射 \mathscr{L}^\zeta 的范数衰减，为了保证前面的 \mathscr{Z}^M 以及收缩映射能够成功的底层算子界限。

（先写到这，熬不动了，我后面再补）

本文少量小节部分内容由我和我的社团集体创作。

参考

1. ^Global solutions of the gravity-capillary water-wave system in three dimensions https://arxiv.org/abs/1601.05685
2. ^Random tensors, propagation of randomness, and nonlinear dispersive equations https://arxiv.org/abs/2006.09285