没有外部热浴，孤立量子系统怎样呈现热平衡？

没有周围的热水，温度从哪里来

把一个小物体放进热水里，它与周围交换能量，最后接近水温。这种热平衡很容易理解。但如果把一团原子或一组量子比特尽量与外界隔开，先制备一个特殊状态，再让它自行演化，局部测量也可能逐渐符合某个温度下的统计分布。1 2

问题在于，理想封闭系统的量子演化是可逆的。已知完整初态和哈密顿量，原则上可以一路算下去；纯态也不会在这种演化中变成整体混合态。为什么局部却像是忘记了许多制备细节，只剩能量、温度等少量参数？

回答需要追踪三个不同环节：整体中的关联怎样改变局部可见的信息，局部波动怎样变得不再明显，以及最后留下的局部统计为什么恰好接近热分布。 少掉任何一环，都可能把“看起来变乱”错当成热化。

本文主要比较有限多体量子系统的实验与理论。制备初态、快速改变参数或缓慢扫描哈密顿量时，外部控制可以对系统做功；我们说的封闭、定能量演化，指随后固定哈密顿量并忽略环境作用的阶段。不能把制备过程中的能量变化与演化阶段的能量守恒混在一起。

整体纯态，为什么局部可以是混合态

“纯态”不是没有涨落，也不是所有粒子各自状态都确定。它表示完整系统可以由一个量子态向量描述。设系统分成我们观察的 $A$ 和暂时不观察的 $B$，局部测量由约化密度矩阵描述：

这个操作不是把 $B$ 物理切掉，也不是让它停止作用，而是在计算 $A$ 的观测概率时不读取 $B$ 的结果。

最简单的例子是两个量子比特处于

整体是一个确定的纯态。单看任何一个量子比特，测到0或1的概率各为一半，而且任何只作用于单个比特的测量都只能看到相应混合态。若联合观察两者，关联则会显露出来。整体纯态与局部混合，并不互相矛盾。1 2

纯度可以把区别写成一个数：

纯态的纯度为1；混合态小于1。上面的整体纯度为1，每个单比特的纯度却为 $1/2$。2016年的原子实验正是同时测整体与局部纯度，直接检查这个区别，而没有只凭一条局部响应曲线推测整体发生了什么。2

产生纠缠，还不够叫作热化

再看一个真正随时间变化的例子。两个量子比特之间只允许一个激发来回交换，从 $|10\rangle$ 开始，在把交换耦合和 $\hbar$ 设为1的时间单位下：

单个比特的纯度随时间变成 $\cos^4t+\sin^4t$。它会降低，说明产生了纠缠；随后又升回1，整个过程周期性重复。下面的曲线由这个模型实际计算，并与矩阵指数演化核对。

这个小反例说明三件事要分开：

“系统的其余部分像一个热浴”可以帮助理解第一步，但不能拿这句话跳过后两步。其余部分本来也未必处在平衡中，必须说明它们如何共同演化。1 4

从许多相位，到少量热学参数

把初态按哈密顿量的能量本征态展开，理想封闭演化只让各部分积累不同相位：

这里的 $|c_n|^2$ 不随时间变化。系统并没有在能量本征态之间不断重新抽签；完整能量分布保留着初态信息。

但一个局部观测量的结果还包含不同本征态之间的交叉项。它们以不同频率摆动，叠加时可能逐渐难以保持同一方向。退相位说的就是这种相位关系在所观察的量上互相抵消的过程，而不是每一个量子相位被环境擦掉。1

不过，求出时间平均并不等于证明所有时刻都接近平均。两比特例子的时间平均也存在，后来的大幅振荡却不会消失。对较大的多体系统，还需要限制相位重合、参与的能量态数量以及观测矩阵元，才能说明长期波动足够小。1 4

假如相关能级没有简并，某个观测量的无限时间平均具有如下形式：

右边仍有初态权重，为什么能变成热学答案？本征态热化假说，简称ETH，提供一种有条件的机制：对适当的局部观测，能量相近的许多本征态给出的结果也相近。只要初态主要覆盖的能量密度范围足够窄，改变这些态的具体权重，对局部观测影响就不大；与此相配的非对角矩阵元性质则限制时间波动。1

所以这里不是“细节已经不存在”，而是所选观测对这些细节不再敏感。若初态覆盖很宽的能量范围、有影响局部统计的额外守恒量，或考察一个精心设计的整体观测量，这种简化都可能失败。

常用的热参照是Gibbs态：

其中 $\beta=1/(k_BT)$。能量守恒帮助选择温度；存在粒子数等约束时，还要在对应约束下建立参照。局部比较应当先构造完整系统的参照，再取约化态。尤其对很小或边界耦合不可忽略的系统，不能未经检查就把局部状态写成孤立小块自身的 $e^{-\beta H_A}/Z_A$。

2016年实验：整体与局部可以同时测

Kaufman等研究的是六个原子在六个格点上的运动。初始每格放一个原子，再允许原子在相邻格点间隧穿。同格点多粒子的相互作用和隧穿共同产生多体演化。2

他们使用两份按相同方式制备的系统，通过两份系统的干涉及粒子数奇偶统计测量纯度。这样可以同时检查两类结果：局部粒子数分布是否接近热参照，以及完整六格系统是否仍保持较高纯度。实验发现，局部统计逐渐呈现热态特征，整体纯度却没有变得像相应混合热态那样低。

这里的实验语言应比理想方程谨慎。测得的整体纯度不等于严格的1，干涉操作、初态制备及两份系统的微小差异都会降低可见纯度。补充材料还说明了残余熵背景的处理，以及按粒子总数后选择数据；在相应数据处理中保留下来的原始运行约占30%。这些校准和选择是解释实验的组成部分，不能被“整体纯态”四个字抹去。2

实验直接读出的是与纯度相连的二阶Rényi熵：

它与常用于热力学的von Neumann熵有关，但不是同一个量。两者都能反映局部混合程度，不能因为都叫熵，就把只对其中一种成立的热力学关系搬到另一种上。

把六格系统算一遍，哪些地方接近，哪些仍不同

本库建立一个对应现象的理想化模型：六个玻色子、六个开放边界格点，初态严格取每格一个粒子的乘积态。用 $J$ 表示相邻格点的隧穿能量，$U$ 表示同格点相互作用，哈密顿量为：

我们采用实验中一组参数的比例 $U/J=103/66$，完整保留固定六粒子空间的462种构型。数值直接求解幺正演化，不加入热浴；但没有复现原实验的有限时长参数切换、初始基态的小修正、测量背景和后选择。因此下图是有原文参数依据的理想化机制计算，不是重拟合实验数据。2

这里的迹距离采用0到1的规范：0表示这部分系统的全部测量统计相同，1表示可以用某种测量将两态完全区分。因而它检查的要求比一条粒子数分布曲线相近更强。

计算中的总能量和能量方差保持不变；整体纯度始终为1。与此同时，局部纯度明显降低。完整六格链在同总粒子数、同平均能量下的Gibbs参照，整体纯度却只有约0.0081。两个整体状态相差很大，并不妨碍一些局部统计相近。

下面左图比较中间单格的粒子数分布。有限窗口平均、谱计算得到的无限时间平均与Gibbs参照已经接近。右图逐步增大所观察的子系统：对演化中的整体纯态，子系统与其补集的纠缠熵相同，观察范围扩展到全部六格时，熵必须回到零；混合热态没有这个要求。

这张图不能用来宣告“小系统已经完全热化”。单格谱时间平均态与Gibbs态的迹距离约0.017；扩大到中间两格，约0.066。而在 $Jt/\hbar\in[10,20]$ 的窗口里，两格瞬时态到Gibbs参照的距离平均约0.181，仍有明显差别。把状态先取时间平均再计算距离，会得到约0.089，不能拿这个较小的数冒充每一时刻都很接近。

这是有限系统中一条有用的提醒：粒子数分布看着相近、时间平均相近和完整局部状态在大多数时刻相近，是递进的要求。我们的计算展示了前两者及后者的残余波动，没有证明热力学极限中的ETH。

理想计算还允许把演化算子反向作用，数值返回初态的误差约 $10^{-13}$。这只是对理想方程及数值演算的核验，并不意味着实验可以轻易精确反转所有自由度。它说明：局部看不见的细节，在这种封闭模型中确实没有被整体演化删除。

2025年实验：同一个能量，也可能留下不同的整体状态

2025年，Andersen等把类似问题推进到二维超导量子模拟平台。平台有69个量子比特，其中一些能量与相关测量使用65比特子系统；实验结合数字操作制备和测量状态，再让模拟的相互作用同时运行。3

对本题最有区分力的部分，并不是器件规模本身，而是用不同方式制备最终平均能量相近的状态，再检查哪些观测会相同。研究者一方面改变参数扫描时间，另一方面改变初始激发数。若最后得到的相关长度、涡旋密度等都只依赖最终能量，而对制备路线不再敏感，就支持这些量已经受到热描述约束。

实验确实观察到：把来自不同初始激发和扫描时间的数据按最终能量重画，相关长度与涡旋密度代理量呈现共同关系。这比只在一种初态上看到“曲线稳定了”更有信息，因为它比较了不同制备历史。3

但能量涨落没有落到同一条曲线上。对固定哈密顿量的封闭演化，$H$ 的整个概率分布都不变，$\langle H^2\rangle-\langle H\rangle^2$ 当然也不必被不同初态重置成相同数值。实验通过二比特和四比特关联重建相应量，看到了仍可区分的制备差异；部分未直接测齐的关联使用近各向同性假设下的插值，并与额外测量和模拟比较。这里的严格守恒论证针对固定的理想XY哈密顿量，实际器件的附加项与噪声需另行处理。3

这使“信息去了哪里”有了很具体的回答：某些信息在关心的少量观测中变得难以察觉，却仍能通过更完整的统计组合被分辨。也正因此，ETH讨论局部或适当物理观测时的适用范围不能省略；它不是所有可构造算符都只由平均能量决定的命题。

这套处理器也不是完全无噪声的封闭系统。论文单列哈密顿量校准、读出错误、激发衰减和后选择处理；有些修正依赖激发数守恒，有些退相干影响仍然保留。参数扫描阶段还存在外部做功。本文采用的是其中“比较制备路线与最终观测”的证据，不把经过误差处理的结果说成完美幺正演化，也不据一个有限平台确认无限系统的全部热化性质。

2026年的严格结果，解决了哪部分

Pilatowsky-Cameo与Choi刊于2026年卷期的定理（2025年12月已在线发表），为“容易制备的非平衡纯态怎样局部变得像Gibbs态”提供了一个严格答案。4 但要先读清它回答的量词。

在其主要设定中，量子比特放在周期晶格上，相互作用具有几何局域性，哈密顿量具有平移不变性、保持不随时间变化，并满足特定的能谱无共振条件。温度需要高于与相互作用有关的阈值。这里的“高温”是相对于这个量子模型的能量尺度，不是一个适用于所有材料的摄氏度界线。

所说无共振，要求除平凡的相同配对外，不同能级之间的能量差不能精确相同；这比仅要求能级本身不同更强。它限制的是不同演化频率精确重合的方式。

对初态也有要求。每个纯态可以由乘积态经过深度受限的电路制备；初态集合的平均能量确定有效温度，而且这个集合在相应复杂度和能量约束下满足文中规定的最大熵条件。定理不是只写“高温、均匀”，就保证任意挑出的纯态都热化。

它得到的是：从这类初态集合中取一个典型初态，经过演化，在典型时间观察固定大小的局部区域，结果与Gibbs态越来越难区分；误差随完整系统增大而减小。 初态平均、时间平均、局部区域大小和系统大小的极限，都是结论的一部分。

最大熵初态集合的限制也不等于“每个初态已经是热态”。集合的平均可以是Gibbs态，集合里的每个成员却仍是容易制备、低复杂度的纯态。定理真正连接的是这些成员的局部动力学与热统计；但它确实限制了允许怎样抽取初态。补充材料另外讨论了看似自然的其他乘积态采样方式，明确保留它们是否热化的未决问题。4 4补

证明思路也能用大白话分成两步。第一步，说明一个典型成员的局部观测接近自己的时间平均。第二步，说明这些时间平均在局部接近热参照。把两项误差相加，才得到局部热化。单独证明第一步，并没有完成第二步。

它还没有保证所有符合条件的系统都很快热化。主要定理使用长时平均；补充材料给出了有限时间版本，但其界依赖能谱中“能隙之差”的最小间距，可能随系统大小迅速恶化，甚至给出指数大的时间界。能够证明最终的典型行为，与能预测实验窗口内多久出现，是不同成果。4补

本库六格模型是开放边界玻色链，并不满足该量子比特晶格定理的全部设定；它承担解释机制的任务，不能当作定理的数值验证。2025年的实际处理器也有边界、噪声与具体制备流程，需要逐项判断与理想定理之间的关系。

哪些信息会使普通热描述失效

能量之外若还有大量独立、会约束局部观测的守恒量，仅靠一个温度就可能不够。某些可积系统可以退相位、使局部观测稳定，却仍保留这些守恒量所携带的初态差异；这时需要包含相应约束的广义Gibbs描述。1

这里要避免一个形式陷阱：任何固定哈密顿量都有 $H^2,H^3$ 等守恒，并不意味着所有系统都因此可积。真正相关的是额外、非平凡且影响局部统计的独立守恒结构。2025实验中整体能量涨落不必统一，与一些局部量可由热描述概括，完全可以同时成立。

2006年的“量子牛顿摆”实验是一个早期提醒：近一维玻色气体经历大量碰撞后，在可观察时间内仍未明显接近通常的热平衡。原论文摘要把近可积性作为解释，同时明确实际俘获势和非理想相互作用会影响这个条件。本文只核到这一原始摘要，用它作有限时间的现象入口，不宣称这些实验气体永远不会热化。5

另外，有限小系统可能有明显的返回与振荡；近似守恒量可以产生长暂态；不同观测的松弛时间也可能相差很大。某个量在一段时间内稳定，只能先回答它的平衡问题。要称为热化，还须改变初态、比较合适的热参照、检查守恒约束，并把内部退相位和环境退相干的影响分开。

从这条机制链可以理解什么

理想封闭量子系统无需先与一个已经热化的外界交换能量，才有可能出现局部热描述。相互作用产生关联和纠缠，许多能量分量的相位演化让局部观测逐渐不再敏感于部分初始相位；在合适的能谱、初态与守恒条件下，剩余局部统计又能由少量热学参数概括。

这条链的每一步都有独立边界：纠缠可以出现而没有平衡，平衡可以出现而不服从普通Gibbs分布，局部热化可以出现而整体仍然保留信息。2016实验、2025不同制备路线的比较和2026有明确量词的定理，从不同角度支持这些区别。本库小系统计算则把它们放进同一套可回查数字里，并保留了有限尺寸下不够接近的部分。

它没有从无条件的量子定律推出整个宇宙的时间方向：我们仍选择了特定制备的非平衡初态，也仍限定所观察的区域与时间。这里解释清楚的是，同一个整体量子演化，为什么能够保留完整细节，同时让局部观测服从简单得多的热学描述。

本库计算代码位于 research/calculations/isolated_quantum_thermalization.py；可下载参数与核验结果、六格演化时间序列、能量本征态权重、局部占据分布及子系统熵。