[論文レビュー] Quantum algorithm for time-dependent differential equations using Dyson series

Dominic W. Berry, Pedro C. S. Costa|arXiv (Cornell University)|Dec 7, 2022

Quantum Computing Algorithms and Architecture参考文献 17被引用数 7

ひとこと要約

本論文は、時間に依存する線形微分方程式を解くための量子アルゴリズムを提示する。Dyson級数を用いて時間順序付きの発展演算子を表現し、ブロック符号化により解を量子状態として符号化し、最適な量子線形システムソルバーを用いて解く。誤差と微分係数への依存が対数的であり、複雑性は発展演時間に対して線形で、精度に対しては対数的である。高次元においては古典的手法に比べて指数的高速化を達成する。

ABSTRACT

Time-dependent linear differential equations are a common type of problem that needs to be solved in classical physics. Here we provide a quantum algorithm for solving time-dependent linear differential equations with logarithmic dependence of the complexity on the error and derivative. As usual, there is an exponential improvement over classical approaches in the scaling of the complexity with the dimension, with the caveat that the solution is encoded in the amplitudes of a quantum state. Our method is to encode the Dyson series in a system of linear equations, then solve via the optimal quantum linear equation solver. Our method also provides a simplified approach in the case of time-independent differential equations.

研究の動機と目的

• 時間に依存する線形微分方程式を、複雑性スケーリングを改善した形で効率的に解く量子アルゴリズムの開発を目的とする。
• 特に次元が増加する際に古典的手法のスケーリングが著しく悪化する状況において、このような方程式を量子的手法で解く課題に取り組むことを目的とする。
• 時間順序付きの発展演算子としてのDyson級数を活用し、それを量子線形システムに符号化して、効率的な解法を実現することを目的とする。
• 誤差と微分係数への依存を対数的とすることで、先行研究と比較して複雑性解析を単純化することを目的とする。
• 時間に依存する場合と時間に依存しない場合の両方に適用可能な統一的かつ簡素化されたフレームワークを提供することを目的とする。

提案手法

• 係数行列 A(t) の時間順序付き発展演算子としてのDyson級数を用いて、解を同次型と非同次型に分解する。
• Dyson級数をブロック符号化を用いてブロック行列構造に符号化し、ハミルトニアンシミュレーションで用いられる手法に類似させる。
• Dyson級数から導かれる線形方程式系を、Costaら（2021）の最適な量子線形システムソルバーを用いて解く。
• ユニタリ操作 U_A、U_b、U_x をそれぞれ行列 A(t)、非同次項 b(t)、初期状態 x₀ の符号化に用いる。
• 高精度の最終量子状態を準備するために振幅増幅を適用し、その呼び出し回数は逆誤差の対数的スケーリングを示す。
• A(t) が対角化可能である必要がないようにし、複雑性の上限を単純化するために対数的ノルムを用いる。

実験結果

リサーチクエスチョン

• RQ1 Dyson級数を量子線形システムに効率的に符号化することで、誤差の対数的スケーリングを達成する時間に依存する微分方程式を解くことは可能か？
• RQ2 時間に依存する微分方程式を解く際の複雑性は、発展演時間 T、誤差 ε、および A(t) や b(t) の微分係数にどのように依存するか？
• RQ3 ハミルトニアンシミュレーションで用いられるブロック符号化アプローチを、時間に依存する線形常微分方程式に適応させ、複雑性解析を改善できるか？
• RQ4 量子アルゴリズムにおける U_A、U_b、U_x の呼び出し回数のスケーリングは何か？ これにより先行研究と比較してどのように異なるか？
• RQ5 時間に依存しない場合の複雑性解析を、時間に依存する定式化の特殊ケースとして簡素化できるか？

主な発見

• U_b および U_x への呼び出し回数は O(ℛλT log(1/ε))、U_A への呼び出し回数は O(ℛλT log(1/ε) log(λT/ε)) であり、追加のゲート数は O(ℛλT log(1/ε) log²(λT/ε)) に比例する。
• 誤差依存が 1/ε の対数的であるため、これは量子アルゴリズムにおいて最適であり、多項式的または線形的スケーリングと比較して顕著な改善を示す。
• 本手法は、古典的ソルバーと比較して次元に指数的高速化を達成する一方で、誤差および A(t) や b(t) の微分係数への依存が対数的である。
• 複雑性スケーリングは、全発展演時間 T およびブロック符号化パラメータ λ に対して線形であり、精度に対しては対数的である。
• ブロック行列に余分な行を追加する必要がなく、ハミルトニアンシミュレーションに類似したDyson級数に基づくブロック符号化を用いることで、先行研究と比較して解析が単純化される。
• 時間に依存しない場合の複雑性式も、誤差に対して対数的スケーリングを維持しつつ、システムパラメータへの依存関係が明確に単純化された形で得られる。