QALB(Quantum Algorithm for Lattice Boltzmann,Itani-Sreenivasan-Succi,Physics of Fluids 36, 017112 (2024))是第一个不丢弃非线性碰撞项的完整量子 LBM 算法:它借 Kowalski (1997) 的映射把非线性动力学变成玻色模式(bosonic mode)的演化,Carleman 线性化阶数被吸收进 Fock 空间截断,从而绕开经典 Carleman 变量的跨格点爆炸Itani et al. 2024 (PoF)。它与显式线路派(carleman-qlbm)同属 Succi 学派,但走的是二次量子化编码(见bosonic-encoding-lbm)。

Kowalski 框架:非线性 ODE → 玻色模式

Kowalski 的一般构造(QALB 引文 [55],J. Math. Phys. 38, 2483 (1997)):对多项式向量场 的非线性 ODE ,把经典变量提升为玻色位置算符 ( 湮灭、 产生算符),动力学由生成函数给出

其中 是 的联合本征态(QALB 式 (67)–(72))。精确解仍是原非线性解的编码;截断每个模式的 Fock 空间到 个能级后,演化变成有限维的——且截断阶正好对应 Carleman 展开阶: 型乘积活在 个量子激发的扇区里。这就是”Carleman 阶数 ↔ 玻色占据数截断”的对应。

编码红利: 个 Carleman 阶只需 个量子比特(二进制占据数表示,QALB 式 (40)、(50);各阶还可按 重标度以改善数值条件)。对比显式 Carleman 线路按组合数增长的变量寄存器,这是二次量子化路线的核心卖点。

关键构件

  • 平衡分布与 Hermite 多项式:LBM 平衡分布本就是 Hermite 展开的截断,其生成函数与谐振子相干态结构天然契合(QALB 式 (74)–(78))——“格子 Boltzmann ↔ 玻色谐振子”的数学纽带。
  • 非厄米哈密顿量: 一般非厄米(耗散方向),QALB 用 LCU(幺正的线性组合)实现其演化;文中另给出厄米化对照版本(附后处理耗散项)。
  • 流送:作用在格点位置寄存器上的幺正平移 (式 (125))——代价是牺牲流送精确性换取幺正性:截断 Fock 空间中的平移不再精确(对照carleman-convergence-analysis中经典侧”流送破坏 Carleman 封闭性”的对称困难)。

复杂度与数值验证

比特账本(QALB 式 (122)、(126)–(128)):格点位置编码 ;位置+离散密度合计 。作者还用湍流自由度估计 给出”表示湍流所需量子比特数” (式 (1))——几百度 的流动只要几十个理想比特。按 Table III 的方案对比,碰撞与二进制位置嵌入的组合是唯一在格点总体积 上不线性缩放的方案;但碰撞成本随时间步数 与速度数 幂律增长,作者明言”尚未达到量子优势”Itani et al. 2024 (PoF)。

数值验证(D1Q3 单格点逐次碰撞,图 9–10):非厄米哈密顿量解在截断 qc=2、 时与经典碰撞解的相对误差保持在 20% 以内;厄米化版本给出系统性的对照。误差依赖初值相对平衡值的距离,而非初值本身Itani et al. 2024 (PoF):

D1Q3 单格点量子解与经典碰撞解对照
图:D1Q3 单格点逐次碰撞——量子(非厄米 H)解与经典解的对照,截断 qc=2、N=3 时误差约 20% 以内(来源:QC040)

定位与批判

QALB 的意义在于证明”非线性碰撞可以不丢”:相比丢弃 的线性化方案(budinski-qlbm)与维度爆炸的显式 Carleman 线路(carleman-qlbm),玻色编码给出第三条路。遗留问题同样明确:Fock 截断误差与非幺正误差的长时间行为尚缺系统控制(综述列为该路线的公开问题Sanavio & Succi 2024 (IntechOpen));20% 级精度距工程需求尚远;且全为理想化数值演示,未涉及真机噪声。同族后续(Itani-Succi 的显式碰撞/流送线路 2301.05762)把同一思想的碰撞算子写成玻色升降算符多项式,利用截断算符的幂零性 使平衡分布指数自动截断(见bosonic-encoding-lbm)。

延伸阅读