[論文レビュー] Quantum algorithm for nonlinear differential equations

Seth Lloyd, Giacomo De Palma+6 more|arXiv (Cornell University)|Nov 12, 2020

Quantum Computing Algorithms and Architecture参考文献 4被引用数 24

ひとこと要約

この論文では、非線形微分方程式を解くための量子アルゴリズムを提示している。量子レジスタにシステム状態を符号化し、複数のエンタングルドコピーを用いて多項式非線形性を非線形シュレーディンガー方程式の近似によってシミュレートすることにより、古典的手法と比較して指数的スピードアップを達成している。アルゴリズムは積分時間に対して2次的にスケーリングし、状態空間次元に対して対数的にスケーリングするため、ナビエ＝ストークス方程式やプラズマ流体力学といった高次元問題の効率的解法が可能である。

ABSTRACT

Quantum computers are known to provide an exponential advantage over classical computers for the solution of linear differential equations in high-dimensional spaces. Here, we present a quantum algorithm for the solution of nonlinear differential equations. The quantum algorithm provides an exponential advantage over classical algorithms for solving nonlinear differential equations. Potential applications include the Navier-Stokes equation, plasma hydrodynamics, epidemiology, and more.

研究の動機と目的

• 非線形微分方程式を解く量子アルゴリズムを構築し、古典的ソルバーに対して指数的優位性を達成すること。
• 非線形方程式に対する従来の量子手法が示す指数的リソーススケーリングの課題を、安定で複数コピーのアプローチにより克服すること。
• 近い将来の量子コンピュータを用いて、流体力学やプラズマモデルのような高次元非線形系の効率的シミュレーションを可能にすること。
• スパースで計算可能なテンソルと安定な数値積分技術を活用することで、アルゴリズムのスケーラビリティと正確性を保証すること。
• 線形系に限らない物理的に重要な非線形ダイナミクスを含む、量子微分方程式ソルバーの適用範囲を拡張すること。

提案手法

• 非線形微分方程式の状態ベクトルを、次元dのヒルバート空間内の量子状態として符号化する。
• 初期量子状態の複数コピーを用いて、非線形シュレーディンガー方程式の近似により多項式非線形項をシミュレートする。
• n個のシステムコピーに作用する時間発展演算子であるハミルトニアンを実装し、非線形ダイナミクスを表すスパーステンソルFによって相互作用を定義する。
• Trotter-Suzuki分解を適用して、多コピー系の時間発展演算子をシミュレートし、有効な単一系ダイナミクスを近似する。
• 数値的精度を保証するため、コピー数nを積分時間Tの2乗に比例してスケーリングし、誤差抑制を1/nの比例で維持する。
• 多コピー系上で量子線形微分方程式ソルバーを実装し、非線形ダイナミクスを効率的にシミュレートする。

実験結果

リサーチクエスチョン

• RQ1 非線形性が標準的量子線形ソルバーの挑戦要因であるにもかかわらず、非線形微分方程式を解く量子アルゴリズムが指数的スピードアップを達成できるか。
• RQ2 量子フレームワーク内で、線形操作のみを用いて非線形微分方程式の非線形性を効率的に符号化・シミュレートする方法は何か。
• RQ3 本アルゴリズムのリソーススケーリングは、積分時間と状態空間次元に関してどのようになるか。また、指数的スケーリングを超えて改善可能か。
• RQ4 本アルゴリズムで用いられる平均場近似が、長時間シミュレーションにおいても有効に保たれる条件は何か。
• RQ5 ナビエ＝ストークス方程式やボルツマン方程式のような物理的に重要な非線形方程式に、本手法をスパースで計算可能な相互作用テンソルを用いて適用可能か。

主な発見

• 本量子アルゴリズムは、状態空間次元dに対して対数的にスケーリングし、積分時間tに対して2次的にスケーリングすることで、古典的ソルバーと比較して指数的スピードアップを達成している。
• 必要とされる量子リソース（コピー数を含む）は、積分時間Tに対して2次的にスケーリングし、従来の方法が示す指数的スケーリングを回避している。
• コピー数nがトレーターステップ数Tよりも著しく大きい場合、誤差抑制が1/nに比例して維持され、精度が保たれる。
• ナビエ＝ストークス方程式、プラズマ流体力学、ボルツマン方程式など、局所的相互作用と保存則に起因するスパースで計算可能なテンソルFを有する方程式に本手法は適用可能である。
• 量子フーリエ変換や量子特異値変換などの技術を用いた量子後処理により、パワースペクトルや主成分といった特徴を抽出できる。
• 各ステップでコピーを破棄する必要があった従来の非線形方程式用量子アルゴリズムに比べ、本手法はリソーススケーリングの点で優れている。