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

1. 核心物理图象

• 核心思想：传统格子玻尔兹曼方法（LBM）中的BGK碰撞项是耗散的、非酉的，无法直接在量子计算机上实现。本文的核心物理图象是：将碰撞过程重新解释为量子态在希尔伯特空间中的纯旋转——当前分布函数对应的量子态 $|\psi\rangle$ 向平衡态 $|\psi_{\rm eq}\rangle$ 旋转一个角度 $\theta = \phi/\tau$（$\phi$ 是两态之间的夹角，$\tau$ 是弛豫时间）。由于旋转天然是酉操作，这就绕开了BGK碰撞的非酉性障碍。

• 贡献概述：作者提出了“旋转碰撞算子”这一酉替代方案，并给出了D2Q5（对流扩散方程）和D2Q9（不可压Navier–Stokes方程）两种模型的门级量子线路。理论上证明了旋转碰撞与经典BGK碰撞的差异仅为 $O(\delta t^2)$，数值实验验证了其精度。

2. 关键术语解释

• 旋转碰撞算子（Rotation Collision Operator）：本文提出的核心新概念。它不是直接在分布函数空间做线性插值（那是非酉的），而是在量子态的振幅空间中，将当前态 $|\psi\rangle$ 绕其与平衡态 $|\psi_{\rm eq}\rangle$ 张成的二维子空间旋转一个角度 $\theta = \phi/\tau$。由于旋转保持态矢的模长不变，该操作天然满足量子力学对酉性的要求。

• 振幅编码（Amplitude Encoding）：将分布函数 $f_i$ 编码为量子态振幅 $\sqrt{f_i}$ 的方式，即 $|\psi\rangle = \frac{1}{\sqrt{\rho}}\sum_i \sqrt{f_i}|i\rangle$。这是连接经典LBM与量子计算的桥梁，也是理解为何BGK碰撞在量子框架下“非线性且非酉”的关键。

• 顺序线路拼接（Sequential Circuit Concatenation）：指将不同时间步的碰撞–流动酉块按顺序连接成一个完整的量子线路，而不需要在中间步骤做投影测量。对于D2Q5模型，平衡态制备线路可预编译并复用，是实现多步酉演化的关键。

3. 主要贡献 (Key Contributions)

• 首次从第一性原理构造了酉碰撞算子：不同于此前将非酉算子分解为酉算子线性组合（需辅助比特和后选择，成功率随系统指数衰减）或线性化碰撞算子（引入近似误差）的方案，本文的旋转碰撞算子在构造上就是严格酉的，不需要辅助比特，也不需要丢弃碰撞项。

• 理论证明了旋转碰撞与BGK碰撞的等价性精度：在连续极限 $\delta t \to 0$ 下，旋转碰撞与经典BGK碰撞的差异为 $O(\delta t^2)$，与标准LBM的截断误差同阶。特别地，当 $\tau = 1$ 时两者精确相等。

• D2Q5模型的平衡态线路可复用性：由于D2Q5的平衡态分布 $f_i^{\rm eq} = \lambda_i \phi$ 是宏观量 $\phi$ 的线性函数，归一化后的平衡态方向 $|\hat{e}_{\rm eq}\rangle$ 与瞬时分布无关。因此平衡态制备线路可一次编译、全程复用，只需更新与瞬时态相关的旋转参数。

• 完整的门级线路构造与数值验证：基于Shende–Bullock–Markov态制备方案和相位反射算子，给出了D2Q5和D2Q9模型的显式量子线路，并通过Taylor–Green涡和Poiseuille流两个基准算例验证了方法的精度和有效性。

4. 研究方法 (Methodology)

• 理论框架：以LBM的碰撞–流动结构为基础，将分布函数编码为量子态振幅 $|\psi\rangle = \sum_i \sqrt{f_i/\rho}|i\rangle$。定义当前态与平衡态之间的夹角 $\phi = \arccos\langle\psi|\psi_{\rm eq}\rangle$，构造在 $\text{span}\{|\psi\rangle, |\psi_{\rm eq}\rangle\}$ 二维子空间中的旋转算子 $U_{\rm col}$（见原文Eq. 2.13），旋转角 $\theta = \phi/\tau$。

• 线路实现：旋转算子通过广义相位反射 $S_\phi(\alpha) = I + (e^{i\alpha}-1)|\phi\rangle\langle\phi|$ 实现，而相位反射又借助SBM态制备线路 $S_\phi(\alpha) = \text{SBM}_\phi P_0(\alpha)\text{SBM}_\phi^\dagger$ 转化为多控相位门。整个碰撞线路为 $\text{SBM}_{\psi_t} \to \text{SBM}_{\rm eq}^\dagger \to P_0(\beta_t) \to \text{SBM}_{\rm eq} \to \cdots$（见原文Eq. B.7）。流动步则可用量子行走或量子傅里叶变换实现。

• 模型差异处理：D2Q5的平衡态方向固定，线路可复用；D2Q9的平衡态依赖于宏观速度 $u$（含 $u^2$ 项），需在每个时间步测量当前态、经典计算 $u$ 和 $f_i^{\rm eq}$、再更新碰撞算子。

5. 实验结果与结论 (Results and Conclusion)

• 主要结论：

1. 旋转碰撞算子与经典BGK算子在Taylor–Green涡和Poiseuille流中给出几乎相同的流场（涡量分布、速度剖面均吻合），且与解析解一致。
2. 两种碰撞算子的差异 $E_u^{Q,\rm BGK} = O(\delta t^2)$，与理论预测一致；当 $\tau = 1$ 时差异降至浮点舍入水平。
3. D2Q5模型的门级线路资源（transpile至 $\{U_3, CX\}$ 基）为：4量子比特、95个CX门、180个U3门、线路深度189（不含态制备）。

• 对领域的意义：本文首次为LBM的碰撞步骤提供了严格酉、无需辅助比特、无需后选择的量子实现方案，为量子计算流体力学（量子CFD）扫除了一个基础性障碍。D2Q5模型的可复用平衡态线路结构尤其有利于近期量子硬件的实现。

• 开放性问题：

1. 态制备与终端测量的成本：复杂度分析未包含态制备和测量，其 $O(N)$ 的标度可能部分抵消量子优势。
2. 流动步和边界条件的量子化：目前流动步和半反弹边界条件仍以经典方式处理，如何完全量子化仍是未解决问题（可能可借助SVD框架）。
3. 全相干的多步演化：D2Q9模型中状态相关参数的“相干确定”需要额外的量子子程序；目前采用每步测量+经典反馈的混合方案。
4. 硬件噪声影响：所有结果均来自经典模拟，真实量子硬件上的噪声效应有待研究。

6. 论文标签 (Tags)

量子算法, 模拟, 编译与优化, 量子信息

📄 点击此处展开/折叠原文 PDF

原文链接： Rotation Collision based Quantum Lattice Boltzmann Methods