[論文レビュー] The cost of solving linear differential equations on a quantum computer: fast-forwarding to explicit resource counts

David Jennings, Matteo Lostaglio+3 more|arXiv (Cornell University)|Sep 14, 2023

Quantum Computing Algorithms and ArchitectureComputer Science参考文献 45被引用数 3

ひとこと要約

この論文は、量子コンピュータ上で線形常微分方程式（ODE）をシミュレートするための、初めての非漸近的リソース推定を提供する。安定性の性質が高速前進（fast-forwarding）を可能にし、負の対数ノルムを有する系では、クエリ複雑性がO(√T)にまで低減されることを示している。特定の状況では、最先端の手法と比較して最大4桁の改善が達成されるが、ハミルトニアン系のダイナミクスは、高速前進の定理に従い、その対象とならない。

ABSTRACT

How well can quantum computers simulate classical dynamical systems? There is increasing effort in developing quantum algorithms to efficiently simulate dynamics beyond Hamiltonian simulation, but so far exact resource estimates are not known. In this work, we provide two significant contributions. First, we give the first non-asymptotic computation of the cost of encoding the solution to general linear ordinary differential equations into quantum states -- either the solution at a final time, or an encoding of the whole history within a time interval. Second, we show that the stability properties of a large class of classical dynamics allow their fast-forwarding, making their quantum simulation much more time-efficient. From this point of view, quantum Hamiltonian dynamics is a boundary case that does not allow this form of stability-induced fast-forwarding. In particular, we find that the history state can always be output with complexity $O(T^{1/2})$ for any stable linear system. We present a range of asymptotic improvements over state-of-the-art in various regimes. We illustrate our results with a family of dynamics including linearized collisional plasma problems, coupled, damped, forced harmonic oscillators and dissipative nonlinear problems. In this case the scaling is quadratically improved, and leads to significant reductions in the query counts after inclusion of all relevant constant prefactors.

研究の動機と目的

• 線形ODEを解くための量子アルゴリズムの明示的・非漸近的リソースカウントを提供し、漸近的解析を超えること。
• 古典的力学系の安定性特性が、量子シミュレーションにおける高速前進を可能にし、シミュレーション時間スケーリングを低減する仕組みを調査すること。
• 量子線形方程式系ソルバー（QLSA）のリソース推定を、履歴状態と最終状態出力も含めた動的系に拡張すること。
• リソース推定の複雑性解析に、リャプノフに基づく安定性（対数ノルム）を統合するフレームワークを確立し、効率の向上を図ること。
• 安定性を持つ系（負の対数ノルム）が、ハミルトニアン系とは異なり、シミュレーション時間Tに対して非線形スケーリングを示すことを実証すること。

提案手法

• 線形ODEをL y = cという線形方程式系に埋め込み、状態準備にための量子線形ソルバー（QLS）サブルーチンを活用する。
• ブロック符号化技術を用いて、系行列Aと時間発展演算子を表現し、QLSAを介した効率的量子シミュレーションを実現する。
• 新しい系ẋ′ = A′x′に変換することで、負の対数ノルムを有する系を構築し、安定性に起因する時間発展演算の圧縮により高速前進を可能にする。
• 時間積分法と履歴状態符号化を用い、離散時間tにおけるx(t)に比例する量子状態を出力し、1ノルムで誤差を有界に保つ。
• 参考文献[3]の非漸近的QLSA境界を用い、前因子と失敗確率補正を含む明示的なクエリ複雑性式を導出する。
• 線形ソルバーへの還元により、複雑性の下界を導出し、最悪ケース仮定下で、安定系はΩ(min(log^γ(1/ε), T^γ))のリソースを要することを示す。

Figure 1: Query counts for negative log-norm systems. We consider homogeneous ODE systems $\dot{\boldsymbol{x}}=A\boldsymbol{x}$ for which $A$ has a negative log-norm $\mu=\lambda_{\textrm{max}}(A/2+A^{\dagger}/2)<0$ , $\|A\|=1$ and set $\epsilon=10^{-10}$ . (Main figure): For all such systems, the

実験結果

リサーチクエスチョン

• RQ1 古典的力学系の安定性特性が、シミュレーション時間Tに対して非線形スケーリングを達成する量子シミュレーションを可能にするか？
• RQ2 故障耐性のある量子コンピュータ上で線形ODEをシミュレートする際の、明示的・非漸近的リソースコスト（クエリ数、キュービット数）は何か？
• RQ3 負の対数ノルムを持つ系による高速前進は、標準的な量子ODEシミュレーション手法と比較して、クエリ複雑性にどのような影響を与えるか？
• RQ4 定数の前因子と非漸近的境界は、従来の漸近的解析と比較して、量子ODEシミュレーションの実用的関連性をどの程度向上させるか？
• RQ5 なぜハミルトニアン系のダイナミクスは高速前進の範囲外にあり、これはノーファストフォワーディング定理と整合的であるのか？

主な発見

• 負の対数ノルム（安定系）を有する線形ODEに対しては、クエリ複雑性がO(√T)にスケーリングされ、最先端の漸近的解析と比較して最大4桁の改善が達成される。
• 履歴状態出力の明示的リソースカウントは、T ∈ [10⁶, 10¹⁵] および ε = 10⁻¹⁰ の範囇で、11,900√T log(T) から 10,300T log(T) の間で有界である。
• 負の対数ノルムを持つ系への変換により高速前進が可能になるが、誤差制約ε′_TD = e^(-zT)ε_TDが元に戻ると、Ω(T)の複雑性下限が回復され、ノーファストフォワーディング定理と整合的であることが保たれる。
• A⁻¹bが平衡点である安定系は、長時間極限において線形ソルバーに還元され、最悪ケース仮定下でΩ(min(log^γ(1/ε), T^γ))の複雑性下限が得られる。
• 本フレームワークは、従来のQLSA境界を動的系に拡張し、失敗確率補正を含む期待クエリ数の明示的式を提供する。
• ハミルトニアン系（A = -iH）は境界ケースである：安定性に起因する高速前進は不可能であり、ノーファストフォワーディング定理と整合的である。

Figure 2: Algorithm schematics. ODE algorithm schematics. (Red, dashed): classical precomputations. (Blue, solid): quantum operations. The algorithm reduces solving the ODE to solving a linear system of equations. O1, O2, O3 are the three oracles given on the previous page.