技術要約：量子斜交固有射影 (Quantum Oblique Eigenprojection)

問題設定

本論文は、量子コンピュータ上で斜交固有射影 (oblique eigenprojections) を実行するという課題に取り組んでいる。エルミート行列（固有空間が直交する）に対する量子アルゴリズムは、量子特異値変換 (QSVT) などのフレームワークを用いて確立されているが、一般的な非正規行列は大きな障壁となる。非正規行列の場合、ヒルベルト空間は一般化された非直交の固有部分空間へと分解される。課題は、ターゲットとなる複素スペクトルの領域（例：単位円板内部の固有値）に関連する特定の不変部分空間への入力状態の射影を構築し、補完的な領域に関連する成分を抑制することである。

この問題は、スペクトルが複素数である場合に根本的に阻害される。実数スペクトルの場合とは異なり、多項式は複素平面上において指示関数を近似することが最大加法原理（最大値原理）によって制限されるため、近似を行うことができない。その結果、標準的なQSVTや量子固有値変換 (QEVT) の手法は、解析関数に対する多項式近似に依存しているため、この問題を直接解くことはできない。

レゾルベント積分（コーシーの積分公式を使用）のような従来のアプローチは、非解析関数を扱うために、ターゲットとなるスペクトルを囲む輪郭に沿って行列のレゾルベント $(\zeta I - A)^{-1}$ を積分することができる。しかし、これらの手法には決定的な非効率性が存在する。すなわち、正規化係数 $\alpha$ が $\Theta(\|\Pi\|/\delta)$ （ここで $\|\Pi\|$ は射影のノルム、$\delta$ はスペクトルギャップ）としてスケールするという問題である。この膨張した正規化係数は、射影をサブルーチンとして使用するあらゆるアルゴリズムのクエリ複雑性に大きな影響を与える。なぜなら、ブロックエンコーディングの成功確率は $\alpha^2$ に反比例するためである。

手法

著者らは、離散化されたレゾルベント積分の生成されるブロックエンコーディングの全構造を利用することで、標準的なレゾルベント積分法を改善する新しいアプローチを提案している。

1. 離散化されたレゾルベント積分: 斜交固有射影 $\Pi$ は、レゾルベントの輪郭積分を台形則を用いて離散化することにより近似される。これにより、シフトされたレゾルベントの和が得られる：$\Pi \approx \frac{1}{n} \sum_{j=0}^{n-1} (I - \omega^{-j}_n A)^{-1}$。
2. ブロックエンコーディングの構築: 標準的な手法では、これらのシフトされたレゾルベントを含むブロック対角行列 $M$ を構築し、QSVTを用いて逆行列を求める。随伴レジスタにアダマール変換を適用することで、結果として得られる逆行列の上左ブロックが $\Pi$ を近似する。しかし、残りのブロックは構造化されていないため、正規化係数は行列全体のワーストケースのノルム、すなわち $\Theta(1/\delta)$ によって決定されてしまう。
3. 離散フーリエ変換 (DFT) 解析: 本研究の鍵となる洞察は、逆行列 $M^{-1}$ のブロックは任意ではなく、巡回構造 (circulant structure) を持っているという点である。$M^{-1}$ をアダマール変換の代わりに量子フーリエ変換 (QFT) で共役変換することにより、すべてのブロック $(k, l)$ が、シフトされたレゾルベントの特定の離散フーリエモード $S_m$ に対応していることが明らかになる。
   • 第0モード $S_0$ は、固有射影 $\Pi$ を近似する。
   • 高次モード $S_m$ は、スペクトル集合条件（例：行列が対角化可能であるか、または有界な数値範囲を持つ場合）の下で幾何級数的に減衰する。
4. 二側ブロック前処理 (Two-Sided Block Preconditioning): 著者らは、フーリエ基底における二側ブロック前処理技術を導入している。これは、逆転させる前に行列を修正するスケーリング演算子 $T$（計算基底において対角）を適用するものである。
   • スケーリングパラメータ $t$ が減少するにつれ、前処理された行列の第1ブロック行および列は増大する。
   • しかし、この前処理は、抽出される上左ブロックの正規化係数の逆ギャップ成長を抑制し、それを $\Theta(\|\Pi\|/\delta)$ から $O(c \|\Pi\|)$ へと減少させる（ここで $c$ はスペクトル集合条件に依存する定数、例えば数値範囲の場合は $c = 1+\sqrt{2}$、対角化可能な場合は条件数である）。この正規化はほぼ最適であり、$\|\Pi\|$ は $\Pi$ の任意のブロックエンコーディングの理論的下限である。

主な結果

本論文は、主に3つの定理を確立している。

• 定理 1 (行列レゾルベントのDFT): レゾルベント積分の離散化のフーリエ基底におけるブロックの閉形式の表現を提供する。各ブロックは、シフトの累乗に伴い、行列の内部および外部の不変部分空間への制限によって決定されることを示している。
• 定理 2 (単位円板射影): スペクトルが単位円と離れている行列 $A$ に対して、内部部分空間への斜交固有射影 $\Pi$ を、正規化係数 $\alpha_\Pi = O(c \|\Pi\|)$ とクエリ複雑性 $O\left(\frac{c \alpha_A \|\Pi\|}{\delta} \log \frac{c \|\Pi\|}{\varepsilon} \log \frac{1}{\varepsilon}\right)$ でブロックエンコードできることを証明している。これは、標準的なアプローチと同じクエリ複雑性を維持しつつ、正規化係数を劇的に改善している。
• 定理 3 (一般領域): 任意の複素平面内のターゲット領域への拡張を行う。解析的な分離写像 $f$（例：ターゲット領域を単位円板に写す有理関数）を用いることで、問題を $f(A)$ に対する単位円板のケースに帰着させる。コストは $f(A)$ の単位円に対するギャップに依存する。

用途

著者らは、この改良されたプリミティブが以下の3つの具体的な用途において大幅な高速化をもたらすことを示している。

1. 複素スペクトルに対する固有状態準備: 非正規行列（複素固有値を持つ）に対してQEVTアルゴリズムを拡張する。新しい手法は、固有状態を $O(\|\Pi\|^2 / (\eta \delta))$ のクエリ複雑性で準備する。これは、標準的なレゾルベント積分アプローチの $O(\|\Pi\|^2 / (\eta \delta^2))$ というスケーリングと比較して、二次的な改善である。この改善は、振幅増幅のラウンド数がギャップに依存せず $O(\|\Pi\|/\eta)$ となることに起因する。
2. 連続時間代数リッカチ方程式: 最適制御や量子化学において現れるこれらの方程式を解くことは、非正規ハミルトニアンの不変部分空間への射影に帰着する。新手法は、Rodenas-Ruiz, Zhao, および Lee による以前のソルバーと比較して、ギャップ依存性を三次的に改善し、総複雑性から $1/\delta$ の因子を2つ取り除いている。
3. 常微分シルベスター方程式: ブロック行列の不変部分空間への射影による $AX + XB = C$ の解決。この手法は、Wang と Liu による直接的な拡張手法のギャップ依存性を二次的に改善する。

重要性と主張

本論文は、本研究が非正規行列に関する量子アルゴリズムの基礎的なプリミティブを提供することを主張している。正規化係数を最適な $\|\Pi\|$ に近づけつつ、逆ギャップに対してほぼ線形のクエリ複雑性を維持することで、著者らは、エルミート行列に対する QSVT の効率性と、非エルミート問題に必要な能力との間の溝を埋めている。

著者らは、以下の理由から本研究が重要であると強調している。

• 非解析関数に対してこれまで問題を抱えてきた「膨張した正規化係数」の問題を解決したこと。
• 非解析的な行列関数（具体的には斜交固有射影）が、多項式近似を使用できないにもかかわらず、QSVT と同様の複雑性で量子コンピュータ上に実装可能であることを示したこと。
• DFT を通じた全ブロック構造の解析と二側前処理の適用という手法が、ブロックエンコーディングを最適化するための新しいパラダイムを提供しており、他の非解析関数にも適用可能であること。

著者らは、その範囲について慎重な姿勢を崩しておらず、前処理はレゾルベントの特定の閉形式の DFT に依存しており、すべての構成（例：Wang と Liu の積レゾルベント法）に即座に適用できるわけではないと述べている。また、QSCT においてコストがターゲット関数の幾何学的形状のみによって決定されるのに対し、本手法のコストは入力行列のスペクトルギャップによって支配されることも指摘している。