Float Lattice Gas Automata: A connection between Molecular Dynamics and Lattice Boltzmann Method for quantum computers

想象一个这样的世界：你无需依赖庞大且昂贵的超级计算机，就能预测一滴墨水在玻璃杯中如何旋转，或者一场风暴如何在城市上空肆虐。这就是计算流体力学（CFD）的梦想，它是致力于模拟液体和气体流动的科学分支。几十年来，科学家们一直在使用两套主要的“规则书”来参与这场游戏。第一种是格子气自动机（LGA），它就像是在网格上进行的跳棋游戏，微小的、坚硬的粒子在其中相互碰撞。它简单且快速，但这个游戏充满了噪声；粒子表现得过于像骰子，产生了一种统计上的模糊性，使得精确预测变得困难。第二种规则书是格子玻尔兹曼方法（LBM），它更加复杂。它不再追踪单个蹦跳的小球，而是追踪在特定位置发现粒子的概率，从而平滑了噪声并提供了更清晰的结果。然而，LBM 的计算量非常大，尤其是在尝试用量子计算机模拟复杂的旋流时，因为量子计算机在处理这些模拟所需的复杂非线性数学方面表现得很糟糕。一个核心问题一直是：我们能否在简单的、充满噪声的跳棋游戏与平滑的、沉重的概率数学之间架起一座桥梁，从而创造出一种既快速又精确，甚至能与下一代量子计算机兼容的工具？

本文介绍了一位领域内的新玩家，名为浮点格子气自动机（FLGA）。可以将 FLGA 视为一种“两全其美”的混合体。作者采用了旧有的、充满噪声的跳棋游戏（Losa），并将坚硬的整数粒子替换成了浮点数——这本质上允许粒子成为整体的一部分，比如 0.5 个粒子。通过这种方式，他们将旧游戏中混乱的随机碰撞转变为一种平滑的、确定性的流动，其表现几乎与高级的 LBM 完全一致，却无需承担沉重的数学负担。团队证明了这种新方法可以准确模拟现实世界的流体行为，例如冲击波和旋涡，其结果与公认的 LBM 高度吻合。此外，他们还构建了一个量子版本的算法（QFLGA），并在一个简单的一维冲击波模拟上成功进行了测试。虽然目前的量子版本在连续运行多个时间步长方面仍面临障碍，但结果表明，这种简化的方法可能是解锁未来量子机器高效流体模拟的关键，有望绕过目前阻碍进展的瓶颈。

新的“浮点”游戏

为了理解为什么这很重要，让我们看看旧的玩法。原始的 LGA 就像一个数字沙盒，每个盒子里都有固定数量的弹珠。如果两个弹珠相撞，它们就会弹开。但由于你只能拥有整数个弹珠（1, 2, 3...），结果是“多噪”的，就像旧收音机里的静电声。LBM 通过使用概率解决了这个问题，但它需要复杂的计算来确定这些概率如何变化，这对量子计算机来说很难快速完成。

作者 Antonio David Bastida Zamora 及其团队决定尝试一条中间道路。他们保留了旧有 LGA 简单的“碰撞”规则，但将弹珠改成了浮点数。想象一下，你不再拥有 3 个完整的弹珠，而是拥有 3.14 个弹珠。这一微小的改变使得系统可以直接计算碰撞的“平均”结果，而不是通过模拟成千上万次随机碰撞来寻找平均值。这消除了噪声，使数学过程变得更加平滑，有效地将 LGA 转变为一种表现得就像 LBM 一样的介观模型。

平滑粗糙的边缘

本文的一个巧妙之处在于他们如何处理粘度，即流体的“厚度”或“黏性”（例如蜂蜜与水的区别）。在旧有的 LGA 中，你必须多次运行模拟才能获得正确的厚度，这非常缓慢。在 FLGA 中，作者引入了一个“旋钮”（他们称之为参数 $C$），可以通过它来调节粒子碰撞的频率。通过转动这个旋钮，他们可以在不改变基本游戏规则的情况下，模拟极薄、极快的流动或浓稠、缓慢的流动。他们甚至发现了一种使用“过度松弛”（over-relaxation）的方法，即让系统稍微超过其平衡状态以加速，从而能够模拟以前难以建模的极低粘度流体。

团队在二维空间测试了这种新方法，模拟了一个泰勒-格林涡（流体的旋转模式）和一个顶盖驱动流腔（一个顶部盖子在滑动从而带动内部流体的盒子）。在这两种情况下，FLGA 的结果与金标准 LBM 的模拟几乎完全匹配。他们还比较了每种方法所需的计算机内存和时间。虽然 LBM 很快，但它需要大量的内存来存储所有中间计算过程。相比之下，FLGA 则更加精简；它直接从粒子的局部分布进行计算，需要的内存更少，完成任务所需的步骤也更少。

量子飞跃

论文中最令人兴奋的部分是尝试在量子计算机上运行它。量子计算机擅长同时处理许多可能性（叠加态），但在处理描述流体碰撞的非线性数学方面却很吃力。作者设计了一个量子电路（一套给量子计算机的指令），使用了振幅编码。他们不再将数字存储在比特中，而是将其存储为量子系统中波的“高度”。

他们构建了一个量子电路，该电路获取两个流体状态的副本，将它们相乘（以模拟碰撞），然后将它们移动到下一个位置。当他们在简单的 1D 冲击波（压力和密度的突然变化）上进行测试时，量子模拟与经典的 FLGA 结果完美契合。这证明了该概念是可行的。然而，这里有一个限制：目前的量子算法一次只能运行一个时间步长。要运行完整的模拟，你必须测量结果，记录下来，然后再将其加载回量子计算机进行下一步。这种“走走停停”的过程消耗了量子计算的加速优势。

作者指出，虽然这是一个局限性，但他们模型的简洁性可能提供一种解决办法。由于 FLGA 模型如此精简，或许可以在不进行完整测量的情况下近似计算下一步，或者使用诸如卡雷曼线性化（一种将非线性问题转化为线性问题的数学方法）之类的巧妙技巧，以让量子计算机运行得更久。

未来何去何从？

论文明确了其研究成果与局限。他们已经模拟并证明了 FLGA 在经典环境下有效，并且其量子版本在单步运行上可行。他们提出，这种方法可以解决量子流体力学的瓶颈，但尚未构建一个用于复杂二维流动的完整多步量子模拟。目前，进行完整二维模拟所需的量子比特（qubits）数量对于现有的机器来说过高。

然而，前行的道路看起来充满希望。通过简化游戏规则并消除噪声，作者在微观粒子世界与宏观流体流动世界之间架起了一座桥梁。如果未来的量子计算机能够处理“走走停停”的问题，FLGA 可能会成为模拟从人体血管中的血流到未来航天器空气动力学等一切内容的必备工具，而且这一切都将在未来的量子硬件上运行。目前来看，它是一个聪明且无噪声的替代方案，让我们离完美的流体模拟梦想又近了一步。