🔬 physics

Rotation Collision based Quantum Lattice Boltzmann Methods

本文提出了一种新型幺正量子格子玻尔兹曼方法，通过引入一种基于旋转的碰撞算子，解决了非幺正 BGK 碰撞算子的根本挑战，并为对流扩散方程和不可压缩纳维-斯托克斯方程构建并验证了门级量子电路。

原作者： Kangyang Zeng, Changshen Huang, Xi Liu, Xiaodong Niu, Zhenhua Chai

发布于 2026-09-16

📖 1 分钟阅读☕ 轻松阅读

原作者： Kangyang Zeng, Changshen Huang, Xi Liu, Xiaodong Niu, Zhenhua Chai

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

几十年来，科学家们一直依赖一种被称为格子玻尔兹曼方法（lattice Boltzmann method）的强大数学工具来模拟流体运动和热量传播。想象一下，试图预测烟雾如何在建筑物周围旋绕，或者血液如何在微小的毛细血管中流动；这种方法将流体分解为由微小点组成的网格，追踪每个点上的粒子如何与邻居碰撞并移动到下一个位置。它是一种折中方案，介于单个分子的混沌运动与河流的平滑流动之间，使其成为解决复杂工程问题的理想选择。然而，随着计算机变得越来越强大，研究人员开始转向量子计算机，这种机器利用量子力学的奇特规则来比今天的超级计算机更快地解决问题。挑战在于，标准的模拟处理粒子碰撞的方式并不符合量子力学的严格规则，因为量子力学要求每一步都必须是完全可逆的，并且保持某种特定的数学平衡。

来自华中科技大学和汕头大学的中国研究团队现在开发了一种弥合这一差距的新方法。他们创建了一个能在量子计算机上自然运行的格子玻尔兹曼方法版本，通过将标准的碰撞步骤替换为“旋转”。在他们的新框架中，不再是强迫流体粒子向一个平静状态弛豫（这会导致信息丢失），而是让量子计算机轻轻地将粒子的状态旋转向那个平静状态。这种旋转与量子力学完美契合，因为它是一种可逆的操作，能够保持系统的平衡。研究人员使用两种不同的模型测试了这个想法：一种用于更简单的热量和质量传递问题，另一种用于更复杂的不可压缩流体流动（例如水在管道中流动）。他们的工作表明，这种新方法可以准确地模拟真实流体的行为，同时遵循严格的量子物理定律。

该团队解决的核心问题是流体行为与量子计算机运作方式之间的根本不匹配。在传统方法中，当粒子发生碰撞时，它们会损失一部分能量并趋于稳定，这是一个不可逆的过程。然而，量子计算机无法执行会丢失信息的运算；每一步都必须是一个幺正变换（unitary transformation），这意味着总信息量保持不变，就像洗一副完美的扑克牌，没有任何牌会丢失一样。以前的修复尝试涉及复杂的变通方法，例如使用额外的辅助比特，这些方法往往会失败，或者需要完全移除碰撞步骤，从而改变了模拟的物理特性。由中国团队提出的新方法通过将碰撞视为数学空间中的一次旋转，而非能量的损失，从而避免了这些陷阱。他们证明，通过将流体粒子的状态旋转到恰当的角度，系统表现得几乎与传统方法完全一致，但又不会破坏量子力学的规则。

为了证明其想法的可行性，研究人员在经典计算机上进行了详细的数值实验，以模拟其新量子方法的表现。他们首先测试了一个经典的流体力学问题——泰勒-格林涡（Taylor–Green vortex），这是一种由于摩擦而随时间衰减的流体旋涡模式。他们将这种基于旋转的方法与标准经典方法进行了对比，发现结果几乎完全相同。流体的旋涡模式（包括涡流的位置和强度）都完美匹配。他们还测试了一种被称为泊肃叶流（Poiseuille flow）的流动，这种流动描述了流体如何在通道中运动，结果同样显示新方法产生的结果与已知的数学解高度一致。团队发现，新方法与旧方法之间的微小差异非常之小，只有在计算使用的步长非常大时才会出现，这证实了该新方法具有极高的准确性。

这种新方法的一个特别巧妙之处在于它如何处理不同类型的流体问题。对于用于研究热量和扩散的较简单模型，研究人员发现，流体试图达到的“目标”状态是固定的，且不会随着模拟过程而改变。这意味着，负责准备该目标状态的量子计算机电路可以构建一次并重复使用，从而节省大量的计算资源。然而，对于用于模拟水流的更复杂的模型，目标状态会根据流体在瞬间的运动速度而改变。在这种情况下，研究人员展示了系统仍然可以运作，即通过在每一步测量流速并相应地更新电路，尽管这需要一个稍微复杂一点的过程。

研究结论指出，这种基于旋转的框架为在量子计算机上运行流体模拟提供了一种坚实的、在数学上严谨的方法。它解决了如何使碰撞步骤在不牺牲精度或不需要不可能的计算量的情况下实现可逆性的长期难题。虽然目前的成果是在经典计算机上进行的模拟，但其底层逻辑是专门为量子硬件设计的。研究人员指出，未来的工作需要解决如何完全在量子框架内实现粒子的运动和模拟的边界，但这种新方法已经将碰撞步骤确立为一个可行的、幺正的过程。通过将耗散碰撞转化为精确的旋转，该团队为将流体力学最重要的工具之一带入量子时代开辟了一条清晰的路径。

技术摘要：基于旋转碰撞的量子格子玻尔兹曼方法

问题陈述
格子玻尔兹曼方法（LBM）是一种广泛使用的介观计算流体力学（CFD）方法，其碰撞-迁移结构非常适合并行计算。然而，将 LBM 转化为量子计算面临两个根本障碍：

幺正性（Unitarity）： 量子演化要求幺正算符，而标准的 BGK 碰撞算符是耗散且非幺正的。
非线性（Nonlinearity）： 在 D2Q9 等模型中，平衡分布函数包含速度的二次项。如何在不违反不可克隆定理的情况下，将这些项直接编码进量子电路，是一个挑战。

现有的应对策略存在显著局限性。一些方法将非幺正算符分解为幺正算符的线性组合，但这导致成功概率随系统规模呈指数级下降。另一些方法则完全避开碰撞项，从而限制了物理范围；或者对算符进行线性化处理，这会引入偏离原始物理模型的近似误差。此外，将碰撞算符表示为状态相关幺正矩阵的方法，通常需要在每个时间步重新计算，导致门复杂度呈指数增长。至关重要的是，目前尚未建立一种能够从第一原理出发、系统性构建并能精确还原 BGK 松弛结构的幺正替代方案。

方法论
本研究提出了一种基于旋转碰撞算符的新型量子格子玻尔兹曼方法（QLBM）。其核心方法是将非幺正的 BGK 松弛映射为振幅空间内分布函数的幺正旋转。

振幅编码： 分布函数 $f_i$ 被编码为量子态 $|\psi\rangle = \frac{1}{\sqrt{\rho}} \sum_i \sqrt{f_i} |i\rangle$ 。在此空间中，BGK 碰撞通过因子 $(1 - 1/\tau)$ 减少与平衡态的偏差。
旋转碰撞算符： 作者定义了当前态 $|\psi\rangle$ 与平衡态 $|\psi_{eq}\rangle$ 之间的夹角 $\phi$ 。BGK 松弛通过一个幺正旋转 $U_{col}$ 来模拟，该旋转将 $|\psi\rangle$ 向 $|\psi_{eq}\rangle$ 旋转角度 $\theta = \phi/\tau$ 。该算符利用广义相位反射和态准备电路（具体采用 Shende–Bullock–Markov 方案）构建。
理论一致性： 本文证明，虽然旋转在希尔伯特空间中是精确的，但将其映射回分布函数空间所得的结果在 $\delta t$ 的一阶项上与经典 BGK 碰撞一致。这种差异被证明是 $O(\delta t^2)$ 阶的，这与标准 LBM 的截断误差一致。
模型实现：
- D2Q5（对流扩散模型）： 对于该模型，平衡分布关于宏观变量 $\phi$ 是线性的。因此，归一化后的平衡态是一个独立的固定向量。这使得平衡态准备电路可以被预编译并在整个模拟过程中重复使用。
- D2Q9（不可压缩纳维-斯托克斯方程）： 在此模型中，平衡分布依赖于局部速度的非线性项。平衡态和碰撞参数随状态而变化。在提出的实现中，通过每个时间步的测量来重建宏观变量，以更新碰撞算符，尽管该框架也支持在指定参数后进行幺正块的顺序串联。

主要贡献

幺正公式化： 本文建立了一种基于振幅空间旋转的幺正碰撞公式，从第一原理出发解决了 BGK 算符的非幺正性问题，无需借助概率性后选择或无碰撞近似。
门级电路： 作者为 D2Q5 和 D2Q9 模型构建了显式的门级量子电路。
- 对于 D2Q5，由于平衡态的固定特性，可以实现显著的电路复用，其中平衡态准备电路仅需编译一次。
- 对于 D2Q9，框架详细说明了如何处理状态相关更新，并比较了不同的电路合成策略（例如，直接幺正矩阵应用与特征基编码）。
顺序串联： 本文概述了一种“顺序电路串联”策略，即在不进行中间投影测量的情况下连接连续时间步的碰撞-迁移块，从而在幺正演化过程中保持量子相干性（假设状态相关参数是通过相干方式确定或通过混合经典-量子更新确定的）。

数值结果
开发的 QLBM 通过经典模拟器进行了数值验证：

泰勒-格林涡（D2Q9）： 对 $Re=100$ 的二维泰勒-格林涡进行的模拟表明，通过旋转碰撞算符获得的涡度场与经典 BGK 算符得到的空间模式几乎完全一致。微小的量级差异归因于振幅演化与分布重建之间的非线性映射。
泊肃流（Poiseille Flow）（D2Q9）： 在 $Re=10, 40,$ 和 $100$ 下的模拟显示，QLBM 速度剖面与解析解及经典格子 BGK 结果吻合良好。
误差缩放： 参数扫描证实，旋转碰撞与经典 BGK 碰撞之间的差异随 $\delta t^2$ 阶缩放，验证了理论分析。其误差行为在时间步长和松弛时间的变化方面与经典方法一致。

意义与主张
本文声称提供了一种幺正且具备门级可实现性的格子玻尔兹曼方法公式，解决了 BGK 算符带来的根本性的幺正性挑战。其主要理论贡献之一是识别出了 D2Q5 模型中可复用的平衡态准备电路结构，这通过减少动态重编译的需求，为量子实现提供了显著优势。

作者谦虚地指出，虽然碰撞算符是幺正的，但目前的 D2Q9 实现依赖于一种混合方法，即通过重建宏观变量来更新碰撞参数。他们承认，如何完全在纯量子框架内相干地确定这些状态相关量仍是一个开放的算法问题。此外，该研究局限于经典模拟；硬件噪声对实际量子设备的影响，以及态准备与测量的完整成本（其缩放为 $O(N)$ ）被确定为未来的研究方向。本文并不声称在不进行测量的情况下完全在纯量子框架内解决 D2Q9 的非线性问题，而是提供了一种系统性的幺正替代方案，用以保留物理上的松弛结构。

类似论文