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

インクの一滴が水の中でどのように渦巻くか、あるいは嵐が街をどのように吹き荒れるかを、巨大で高価なスーパーコンピュータを一切使わずに予測できる世界を想像してみてください。これは、液体や気体の流れをシミュレートすることに捧げられた科学の一分野である、**計算流体力学（CFD）**の夢です。数十年にわたり、科学者たちはこのゲームを進めるために2つの主要な「ルールブック」を使用してきました。第一のものは、**格子ガスオートマトン（LGA）であり、これはグリッド上で小さな硬い粒子が互いに跳ね返り合うチェッカーのようなものです。これは単純で高速ですが、ゲームには「ノイズ」があります。粒子が少しダイス（サイコロ）のように振る舞いすぎるため、統計的な不確かさが生まれ、精密な予測を困難にします。第二のルールブックは、より洗練された格子ボルツマン法（LBM）**です。これは個々の跳ね返るボールを追跡する代わりに、特定の場所に粒子が存在する「確率」を追跡することで、ノイズを滑らかにし、よりクリアな結果をもたらします。しかし、LBMは計算負荷が高く、特に量子コンピュータ上で複雑に渦巻く流れをシミュレートしようとする際、量子コンピュータが苦手とする複雑で非線形な数学的处理に直面します。大きな疑問はこうでした。「単純でノイズの多いチェッカーのゲームと、確率という滑らかで重厚な数学との間に架け橋を築き、高速かつ精密で、おそらく次世代の量子コンピュータにも適したツールを作れるだろうか？」ということです。

本論文では、この分野に新たなプレイヤーとして、**フロート格子ガスオートマトン（FLGA）を紹介します。FLGAを「両方の良いとこ取り」のハイブリッドと考えてください。著者たちは、古いノイズの多いチェッカーのゲーム（LGA）を取り上げ、その硬い整数粒子を浮動小数点数（フロート）に置き換えました。これは、粒子を「0.5個の粒子」のように、全体の一部としての分数として扱うことを意味します。これによって、彼らは古いゲームの混沌としたランダムな跳ね返りを、洗練されたLBMとほぼ同じように見え、振る舞う、滑らかで決定論的な流れへと変貌させました。チームは、この新しい手法が衝撃波や渦巻きといった現実世界の流体挙動を正確にシミュレートでき、確立されたLBMの結果と一致することを証明しました。さらに、彼らはこのアルゴリズムの量子版（QFLGA）**を構築し、単純な1次元衝撃波シミュレーションでテストすることに成功しました。現在、量子版は連続して多くのタイムステップを実行するという課題に直面していますが、その結果は、この簡略化されたアプローチが、将来の量子マシンにおける効率的な流体シミュレーションを解き放つ鍵となる可能性を示唆しており、これまで進展を遅らせてきたボトルネックを回避できる可能性があることを示しています。

新しい「フロート」ゲーム

なぜこれが大きな出来事なのかを理解するために、従来の遊び方を見てみましょう。オリジナルのLGAは、各ボックスに決まった数のビー玉が入ったデジタル・サンドボックスのようなものでした。もし2つのビー玉が衝突すれば、それらは跳ね返りました。しかし、ビー玉は整数（1、2、3...）でしか存在できなかったため、結果には古いラジオの静電気のような「ノイズ」が含まれていました。LBMは確率を用いることでこれを解決しましたが、その確率の変化を計算するために複雑な計算が必要であり、それが量子コンピュータにとって高速に行うのが難しい部分でした。

著者であるアントニオ・ダビド・バスティダ・サモラとそのチームは、中間的な道を選びました。彼らは古いLGAの単純な「衝突」ルールは維持したまま、ビー玉を**フロート（浮動小数点数）**に変更したのです。例えば、3個の丸ごとのビー玉を持つ代わりに、3.14個のビー玉を持つ状況を想像してください。この小さな変更により、システムは衝突の平均的な結果を求めるために何千回ものランダムな跳ね返りをシミュレートする必要がなくなり、直接「平均的な」結果を計算できるようになります。これによりノイズが排除され、数学的に非常に滑らかになり、実質的にLGAをLBMと同じように振る舞うメゾスコピック（中間規模）モデルへと変貌させたのです。

粗いエッジを滑らかにする

この論文における巧妙なトリックの一つは、粘性（流体の「粘り気」や「粘り強さ」、例えばハチミツと水の差）をどのように扱うかです。古いLGAでは、正しい粘性を得るためにシミュレーションを何度も実行する必要があり、時間がかかりました。FLGAでは、著者らは「つまみ」（$C$ と呼ばれるパラメータ）を導入し、粒子の衝突頻度を調整できるようにしました。このつまみを回すことで、根本的なルールを変えることなく、非常に薄い高速の流れや、粘り気の強い遅い流れをシミュレートできます。彼らはさらに、「オーバーリラクゼーション（過緩和）」という、システムを平衡状態からわずかに押し出すテクニックを見出し、これにより、従来はモデリングが困難だった極めて低粘性の流れのシミュレーションを可能にしました。

チームはこの新手法を2次元でテストし、テイラー・グリーン渦（流体の渦巻くパターン）と、リッド駆動キャビティ（上部の蓋が滑ることで内部の流体を引きずる箱）をシミュレートしました。どちらの場合も、FLGAの結果はゴールドスタンダードであるLBMのシミュレーションとほぼ完璧に一致しました。彼らはまた、各手法がどれだけのコンピュータメモリと時間を必要としたかも比較しました。LBMは高速ですが、すべての中間計算を保存するために多くのメモリを必要とします。対照的に、FLGAはより軽量です。粒子の局所的な分布から直接すべてを計算するため、メモリ消費が少なく、作業を完了させるためのステップ数も少なくなります。

量子への飛躍

この論文で最もエキサイティングな部分は、これを量子コンピュータで実行しようとする試みです。量子コンピュータは多くの可能性を同時に扱うこと（重ね合わせ）には長けていますが、流体が衝突する様子を記述する非線形な数学には苦戦します。著者らは、振幅エンコーディングを用いた量子回路（量子コンピュータへの一連の指示）を設計しました。数値をビットとして保存する代わりに、量子システムにおける波の「高さ」として保存するのです。

彼らは、流体の状態の2つのコピーを取り込み、それらを掛け合わせ（衝突をシミュレートするため）、次にそれらを次の場所へ移動させる量子回路を構築しました。これを単純な1次元衝撃波（圧力と密度の急激な変化）でテストしたところ、量子シミュレーションは古典的なFLGAの結果と完璧に一致しました。これはコンセプトが機能することを証明しています。しかし、一つ注意点があります。現在のところ、量子アルゴリズムは一度に1タイムステップしか実行できません。完全なシミュレーションを実行するには、結果を測定し、書き出し、次のステップのために再び量子コンピュータにロードする必要があります。この「ストップ・アンド・ゴー」のプロセスが、量子コンピューティングの速度上の利点を削いでしまいます。

著者らは、これが制限ではあるものの、彼らのモデルの単純さがこの問題への解決策を提供できる可能性があると示唆しています。FLGAモデルは非常に合理化されているため、完全な測定を行わずに次のステップを近似したり、あるいはカルレマン線形化（非線形問題を線形問題に変換する数学的手法）のような巧妙なテクニックを用いて、量子コンピュータをより長く稼働させ続けたりすることが可能かもしれません。

次なるステップは？

論文は、自分たちが成し遂げたことと、まだ成し遂げていないことを明確にしています。彼らは、FLGAが古典的に機能すること、そしてその量子版が単一のステップにおいて機能することをシミュレートし、証明しました。彼らは、このアプローチが量子流体力学のボトルネックを解決できる可能性があると示唆しましたが、複雑な2次元フローに対する完全なマルチステップの量子シミュレーションを構築したわけではありません。完全な2次元シミュレーションに必要な量子ビット（qubit）の数は、現在のマシンにとっては高すぎます。

しかし、前途は明るいです。ゲームのルールを簡略化し、ノイズを取り除くことで、著者らは微視的な粒子の世界と巨視的な流体の流れの世界の間の架け橋を作りました。将来の量子コンピュータがこの「ストップ・アンド・ゴー」の問題に対処できるようになれば、FLGAは、血管内の血流から将来の宇宙船の空気力学に至るまで、あらゆるものをシミュレートするための主要なツールとなり、将来の量子ハードウェア上で動作することになるでしょう。現時点では、それはノイズのない巧妙な代替案であり、完璧な流体シミュレーションという夢に一歩近づいたものです。