Simulating generalised fluids via interacting wave packets evolution

本文介绍了一种高效的模拟框架，该框架将广义流体力学建模为相互作用的半经典波包气体，从而能够对带有非积分扰动的拟可积系统进行快速的大规模研究，并揭示了即使在局部可观测物理量表现出热化时，长程相关性仍可以无限期持续存在。

要点

想象一个拥挤的舞池，每个人都随着一种非常特定且复杂的节奏在律动。在物理学世界中，这就像是一个一维粒子系统（比如细管中的原子），它们是“可积的”（integrable）。这意味着它们遵循严格且可预测的规则，彼此碰撞时永远不会真正丢失各自的能量，也不会变得“混乱”。

长期以来，科学家们一直有一种描述这群人平均运动的好方法，叫做广义流体力学（Generalized Hydrodynamics, GHD）。你可以把 GHD 想象成对舞池进行的天气预报：它告诉你哪里人群密集，哪里人群稀疏，以及运动的“风”是如何流动的。

问题所在：
现实生活并不完美。有时，舞池并不完全平坦（外部势阱），或者舞者会撞到一些不该撞的东西（破坏可积性的扰动）。当这些微小的缺陷发生时，旧的“天气预报”（GHD）就会失效。计算变得异常困难，并且它无法预测当系统试图趋于平衡（热化）时所产生的那些微小的、混沌的涨落。这就像是试图用一张忽略了阵风的地图来预测一场风暴。

新的解决方案：“幽灵”舞者
论文作者提出了一种聪明的新方法来模拟这些系统。他们不再试图求解整个人群复杂的数学方程，而是将系统想象成一团由半经典波包组成的粒子气体。

这里有一个创意类比：
想象真实的、相互作用的舞者很难追踪，因为他们不断地推搡和拉扯。作者建议，我们假装这些舞者实际上是**“幽灵”舞者**（称为“裸粒子”），他们沿着直线行走，彼此从不接触。

然而，这里有一个魔术技巧：

1. 我们追踪这些沿着直线移动的幽灵舞者。
2. 然后，我们应用一个数学“透镜”或映射，将他们的直线位置转换为真实舞者的扭动位置。
3. 这个映射考虑到了当真实的舞者靠近时，他们实际上会如何“移动”彼此的位置（就像硬棒相互碰撞一样）。

为什么这很酷？

• 速度快： 追踪直线对于计算机来说非常容易。复杂的“碰撞”过程由最后的数学透镜来处理，而不是在实时模拟中处理每一次碰撞。
• 能够处理混沌： 如果你在舞池里加了一个凸起（外部势能）或稍微改变了规则，你只需要改变幽灵舞者的运动方式。数学透镜会自动调整，以展示真实的人群如何做出反应。
• 捕捉“细节”： 旧的方法忽略了微小的、随机的抖动（涨落）。这种新方法自然地包含了这些细节，就像真实的群体中人们会有脚步挪动，而不仅仅是齐步走。

大惊喜：“长程宿醉”
研究人员使用这个新工具研究了当舞池弯曲（像个碗或陷阱）时会发生什么。他们原本预期人群最终会平静下来，变成一种随机的热力学混乱状态（平衡态）。

他们发现了一些令人惊讶的现象：

• “表面”看起来很平静： 如果你从远处观察人群（仅仅检查平均速度或密度），它看起来似乎已经稳定下来，达到了一个平和的热平衡状态。
• “记忆”依然存在： 然而，如果你仔细观察人群中不同部分是如何相互连接的（相关性），你会发现它们在极长的距离内仍然保持着联系。这就像是人群虽然看起来已经放松了，但仍然记得很久以前跳过的一个特定舞步。

结论：
这篇论文表明，即使一个系统看起来已经“热化”（达到了某种稳态或随机状态），它实际上可能仍处于一种由于这些隐藏的长程连接而导致的、远离平衡的长期状态中。这种“幽灵舞者”模拟法证明了，真正的弛豫过程比我们之前认为的要长得多，尤其是在受限空间内。

简而言之，他们通过追踪“幽灵”而非“真实”粒子，构建了一种更快、更智能的方法来模拟拥挤的量子系统，并发现这些系统保持记忆的时间比我们想象的要长得多。

技术摘要

技术摘要：通过相互作用波包演化模拟广义流体

问题陈述
广义流体力学（GHD）为描述一维可积及拟可积系统的输运提供了一个通用框架。虽然标准的欧拉尺度 GHD 能够准确捕捉平均密度和电流的演化，但在引入破坏可积性的扰动（如外部捕获势或多体相互作用）时，它无法描述涨落，且在数值计算上变得难以处理。现有的扩展方法（如涉及纳维-斯托克斯修正的扩散 GHD）依赖于长程相关性在流体单元之间可以忽略不计的假设。然而，最近的理论研究表明，中间时间尺度的弹道输运会产生持久的长程相关性，这使得作为标准扩散流体力学基础的局部平衡假设失效。求解需要捕捉这些相关性的完整方程组（涉及二点函数及更高阶函数）在计算上是极其昂贵的，对于拟可积系统而言尤其如此。

方法论：波包气（WPG）
作者提出了一种被称为“波包气”（Wave Packet Gas, WPG）的替代性粒子表示法，用于描述广义流体。该方法将广义流体建模为一组半经典、相互作用的波包（类孤子粒子）组成的“气体”，其轨迹被映射到非相互作用的“裸”粒子轨迹上。

1. 裸-相互作用映射（Bare-Interacting Mapping）： 该方法的核心是一个将相互作用坐标 $x_i$ 与裸坐标 $X_i$ 联系起来的变换。对于一般的可积模型，该变换定义为：
   $$X_i = x_i + \frac{1}{2} \sum_{k} a_{ik} \text{sign}(x_i - x_k)$$
   其中 $a_{ik}$ 是模型相关的散射位移（例如，对于硬棒模型为常数，对于 Lieb-Liniger 模型则随动量变化）。裸粒子以速度 $v^{\text{bare}}(\theta_i)$ 进行弹道运动，而相互作用粒子则经历由映射导出的有效力。
2. 处理破坏可积性： 当存在破坏可积性的项（例如外部势 $V(x)$）时，该映射将相互作用哈密顿量转化为受状态相关、非局部有效势作用的裸粒子系统。裸粒子的运动方程变为：
   $$\dot{X}_i = v^{\text{bare}}(\theta_i), \quad \dot{\theta}_i = -\left. \partial_x V(x) \right|_{x=x_i}$$
   这要求在每个时间步长内进行裸坐标与相互作用坐标之间的转换，以评估力，从而有效地将多体复杂性编码进单粒子动力学中。
3. 统计采样： 为了恢复流体力学行为，该方法对初始态系综进行采样。对于经典系统，通过泊松点过程初始化裸位置，并从目标分布中采样快速度。对于量子系统，通过相空间分箱和占据概率引入费米-狄拉克或玻色-爱因斯坦统计。随后通过反转映射（通常通过对凸曲面进行梯度下降最小化）获得相互作用坐标。
4. 涨落与相关性： 不同于以往失败的“跳蚤气”（flea gas）算法，WPG 自然地嵌入了正确的初始涨落。通过对系综求平均，该方法不仅能重现欧拉尺度的演化，还能重现高阶流体力学修正，包括扩散扩散和二点相关函数，而无需显式求解耦合的积分-微分方程。

主要贡献与结果

• 高效模拟拟可积系统： WPG 算法成功模拟了存在破坏可积性扰动（外部势和二体相互作用）的系统，在这些系统中，由于强梯度和“湍流”相的形成，求解扩散 GHD 的偏微分方程在数值上是不稳定的。
• 重现已知极限： 该方法定量地重现了标准 GHD 在可积极限下的结果，并在短时间内（$t \ll \ell$）与扩散纳维-斯托克斯预测相匹配，这构成了一个非平凡的基准。
• 长程相关性的持续性： 作者证明，在受外部势（谐振、四次或余弦阱）约束的系统中，局部一点观测值（如快速度分布）在扩散时间尺度（$t \sim \ell^2$）上可能表现出热化特征（收敛于高斯分布）。然而，二点相关函数（密度-密度和能量-能量相关函数）保持有限并表现出长程行为（相关长度 $\sim O(10\ell)$），即使在这些时间尺度下也是如此。
• 热化时间尺度： 结果表明，由所有相关性衰减定义的真正热化并不发生在扩散时间尺度。相反，二点函数的弛豫需要更长的时间尺度（$t \sim \ell^3$），这涉及到三点相关性，而标准扩散流体力学无法捕捉到这一点。
• 对初始条件的鲁棒性： 研究表明，长程相关性的产生和持续性对初始态非常敏感。诸如“切割”热云（创造牛顿摇篮装置）之类的方案会放大这些非热相关性，使其在最先进的实验平台中变得可探测。

意义与主张
本文声称 WPG 框架为模拟广义流体提供了一条透明且高效的路径，它是一种自然的硬棒粒子推广，能够自动嵌入 GHD 的精确涨落流体力学扩展。

其核心意义在于证明了在存在弱破坏可积性的情况下，真正的热化并未在扩散时间尺度上达到。作者认为，虽然局部观测值可能暗示系统处于热态，但系统会长时间维持远离平衡态的长程相关性。这挑战了以往关于受限可积系统热化的假设，并表明相关性层级（涉及三点及更高阶函数）在弛豫过程中起着至石之作用。

该方法被呈现为一个通用的工具，用于探索一维流体的现象，在严格可积模型与受现实扰动支配的机制之间进行插值，包括与冷原子实验相关的复杂的破坏可积性的二体相互作用势（例如偶极相互作用）。作者的结论指出，我们的发现使得有必要重新评估拟可积系统的热化问题，并强调了使用定制初始条件来探测这些长程相关性动力学的潜力。