Quantum Simulation of Stokes Flow via Schrödingerisation and Artificial Compressibility

技術要約：シュレディンガー化と人工圧縮性を用いたストークス流の量子シミュレーション

問題提起
非圧縮ストークス流のシミュレーションは、マイクロ流体力学、低レイノルズ数流体力学、および多孔質媒体内の輸送などのアプリケーションにおいて極めて重要である。しかし、古典的な数値手法は、ストークス系の鞍点構造（saddle-point structure）という根本的なボトルネックに直面している。これは、離散レベルにおいて非圧縮条件（$\nabla \cdot u = 0$）を厳密に満たす必要があるためである。この結合は、大規模で条件数の悪い線形システムをもたらし、問題の次元が増大するにつれて計算リソースの要求量が指数関数的に増大するという、計算コストの高い課題を引き起こす。人工圧縮性のような古典的な正則化手法も存在するが、高次元の設定においては、必要なリソースの指数関数的なスケーリングにより、依然として計算上の制約を受ける。

手法
著者らは、人工圧縮性と**シュレディンガー化（Schrödingerisation）**技術を統合することで、これらの課題に対処する統一された量子アルゴリズムの枠組みを提案している。

1. 人工圧縮性による正則化: 厳密な非圧縮条件を、擬似圧力の発展方程式を導入することで緩和する。具体的には、連続式 $\nabla \cdot u = 0$ を $\nabla \cdot u = \varepsilon p$ に置き換える。ここで、$\varepsilon$ は小さな人工圧縮性パラメータである。これにより、元の鞍点問題を、運動量方程式と放物型の圧力発展方程式を含む結合された初期値問題へと変換する。著者らは、正則化された解が真の非圧縮解に対して線形収束（$O(\varepsilon)$）することを示す誤差評価を提供している。

2. スタガード格子（Staggered-Grid）による離散化: 数値的安定性と適切な速度・圧力の結合を確保するために、スタガード格子有限体積法を用いてシステムを離散化する。この手法では、速度成分と圧力を異なるノード点に配置することで、コロケート格子（collocated grids）で一般的に見られる偽の振動を回避している。

3. シュレディンガー化マッピング: 本研究の核心となる革新は、非エルミートかつ散逸的な偏微分方程式（PDE）を、ユニタリな量子系へとマッピングすることにある。 「歪んだ位相変換（warped phase transformation）」 ($w = e^{-q}u$) と補助変数 $q$ に関するフーリエ変換を用いることで、元の動力学を高次元のヒルベルト空間へと埋め込む。この埋め込みにより、散逸的なシステムを、エルミート・ハミルトニアン $H_{Stokes}$ を持つシュレディンガー型の方程式（$\partial_t \hat{v} = i H_{Stokes} \hat{v}$）へと変換する。

4. 量子回路の構築: 著者らは、時間発展演算子 $U_{Stokes}(\tau) = \exp(i H_{Stokes} \tau)$ をシミュレートするための明示的な量子回路を設計している。

   • ハミルトニアンは、第1次Lie-Trotter-Suzuki公式を用いて基本要素へと分解される。
   • 差分演算子（シフトおよびラプラシアン）は、バイナリ符号化とシフト演算子（$S_+, S_-$）を用いて表現される。
   • 回路は、補助変数 $q$-空間と運動量 $\eta$-空間の間を切り替えるための量子フーリエ変換（QFT）を利用し、非ユニタリな元の動力学をユニタリな量子発展によってシミュレートすることを可能にする。
   • 解は、逆QFTと補助変数部分空間（$q > 0$）における射影測定を通じて回収される。

主な貢献

• 統一された枠組み: 本論文は、人工圧縮性とシュレディンガー化を組み合わせた、ストークス方程式に対する初の明示的な量子アルゴリズムを提示しており、結合された速度・圧力システムに対して一貫した定式化を提供している。
• 明示的な回路設計: 著者らは、ハミルトニアンを基本ゲート（CNOT、アダマール、位相、および多制御回転ゲート）へと分解する詳細な量子回路実装を導出し、厳密なゲート数解析を提供している。
• 複雑性解析: 厳密な複雑性解析により、理論的な量子優位性を確立している。著者らは、高次元の問題において、提案手法が古典的手法と比較して問題の次元数に対して指数関数的な加速を達成することを実証している。具体的には、古典的な複雑性が特定の仮定の下で $O(\delta^{-(d+4)})$ とスケールするのに対し、提案手法は次元 $d$ への指数関数的な依存を回避している。
• 数値的検証: Qiskitプラットフォーム上で行われた数値シミュレーションにより、アルゴリズムの妥当性とスケーラビリティが裏付けられている。これには、ベンチマーク問題に対する精度テストやパラメータ収束解析が含まれる。

結果
Qiskitによる数値シミュレーションは、アルゴリズムの精度を検証している。著者らは、解析解を持つ結合ストークス系に対して本手法をテストし、離散化パラメータ（$n_x, n_q$）が増加し、時間ステップ（$dt$）および人工圧縮性パラメータ（$\varepsilon$）が減少するにつれて、点別誤差（pointwise error）が減少することを実証した。これらの結果は、微分演算子の離散化の正当性と、量子回路設計全体の安定性を裏付けている。

意義と主張
本論文は、古典的なストークス流シミュレーションに固有の次元のボトルネックを克服することにより、計算流体力学における重要なブレイクスルーを提供すると主張している。非ユニタリな鞍点問題をユニタリな量子発展へと変換することで、本手法は高次元の流体力学問題を解くための量子コンピューティングの可能性を活用している。著者らは、提案手法が問題の次元数に対して「証明可能な指数関数的加速」を提供することを強調しており、複雑な幾何学的形状や高次元の流体力学における将来的なアプリケーションへの有力な候補としている。本研究は、このようなアルゴリズムを近未来の量子ハードウェアへ移植し、最終的には非線形ナビエ・ストークス方程式へと拡張するための基礎的な一歩として位置づけられている。