Efficient thermalization of one-dimensional noncommuting Hamiltonians via system-bath interactions with $\Theta(n)$ ancilla qubits

=== 要約 ===

技術的要約：システム・バス相互作用を介した1次元非可換ハミルトニアンの効率的な熱化

問題提起
量子多体系におけるギブス状態（熱状態）の準備は、物性物理学や量子化学への応用において極めて重要な課題である。リンドブラディアン（Lindbladian）に基づく手法は、近年、1次元（1D）の非可換ハミルトニアンに対して高速な混合（多項式対数または線形なシステムサイズ依存性）を実現しているが、その回路実装は複雑であり、複数のアンシラレジスタによって制御される順方向および逆方向のハミルトニアン発展を必要とする。

対照的に、システム・バス相互作用アルゴリズムは、より単純な実装モデルを提供する。これは、システムがリセット可能なアンシラ量子ビット（バス）と有限時間相互作用し、その後バスを破棄するというものである。しかし、1D非可換鎖に対するシステム・バスサンプラーの既存の保証は、リンドブラディアン手法の$\tilde{O}(n)$のスケーリングに対し、$\tilde{O}(n^6/\varepsilon)$のスケーリングを持つなど、著しく弱いものであった。中心となる課題は、システム・バスのラウンドが「ラムシフト（Lamb shift）」（コヒーレントなユニタリ誤差）および有限相互作用によるバイアスにより、ギブス状態を近似的にしか保存できない点にある。さらに、並列バスアーキテクチャ（$\Theta(n)$個のアンシラを結合）へと相互作用をスケールアップさせると、以前は制御が困難であった広範な誤差項が生じる。

手法
著者らは、$n$サイトの1D鎖の各サイトに対して、$4q-1$個のリセット可能なバス量子ビットを並列に結合するシステム・バス相互作用アルゴリズムを提案する。このアルゴリズムは、ランダム化された相互作用時間$T$の下で、システムとバスが時間依存ハミルトニアンの下で発展するラウンドで動作する。

急速な収束の証明は、主に3つの技術的ステップに基づいている。

1. チャネル近似: 物理的チャネル $\Phi_{\alpha, T}$ は、理想的なチャネル $\tilde{\Phi}_{\alpha, T} = U(T) e^{\alpha^2 L} U(T)$ （ここで $L$ は2次のダイソン展開から導出された有効なリンドブラディアン）によって近似される。著者らは、物理的チャネルと理想的チャネルの間の誤差が、KMS重み付きノルムにおいて $O(n\alpha^4)$ でスケールすることを証明した。これは、物理的チャネルのダイソン級数の項を、リンドブラディアンの指数展開と注意深く一致させ、1次元の複素時間局所性を利用して交換子誤差と時間切断誤差を抑えることによって達成される。決定的なのは、誤差の解析がトレースノルムではなく、$\chi^2$ダイバージェンスノルム（KMSノルム）において行われ、これにより誤差のスケーリングにおける $n$ の追加因子を回避している点である。

2. 近似定常状態の構成: 理想的な生成子 $L$ は、KMS詳細釣合い部分 $L_{KMS}$ とラムシフト項 $-i[H_{Lamb}, \cdot]$ に分解される。$H_{Lamb}$ はギブス状態 $\rho_\beta$ と可換ではないため、$L$ は $\rho_\beta$ を保存しない。著者らは、ラムシフトによる一次のドリフトを補正する近似定常状態 $\rho^*$ を構築する。これには、補正カーネルのキャンセル方程式を解くプロセスが含まれる。主要な技術的貢献は、「定常性欠陥（stationary defect）」（$\rho^*$ が固定点であることからの距離）を $O(n\alpha^4)$ に抑えることである。これは、ラムシフト項が $n$ 個の準局所項の和であり、ギブス中心の和に作用する際に打ち消し特性を示し、標準的な $O(n^2)$ ではなく $O(n)$ のモーメント境界をもたらすことを証明することによって達成される。

3. 収縮および混合時間の解析: 著者らは、平均化された物理チャネルの収縮を解析する。相互作用時間 $T$ のランダム化を利用してデフェージング（脱位相）を誘起させる。このデフェージングは、ラムシフトの高周波成分を抑制し、一方で散逸部分のダイナミクス（$L_{KMS}$ によって制御される）はシステムを平衡へと駆動する。1D系における強化された散逸推定（先行研究 [CD26] より）とデフェージング機構を組み合わせることで、チャネルが1ラウンドあたり $(1 - \Theta(\alpha^2))$ の割合で定常状態への距離を収縮させることを証明する。

主な貢献

• 並列バスアーキテクチャ: 本論文は、$\Theta(n)$ 個のアンシラ量子ビットがシステムと同時に相互作用する並列バスアーキテクチャを導入し、解析する。これにより、シングルアンシラ方式と比較して更新率が向上する。
• 改善された複雑度境界: 誤差 $\varepsilon$ 内でギブス状態を準備するための期待総ハミルトニアン発展時間を $\tilde{O}(n/\varepsilon)$ であると証明した。これは、既知の最適なリンドブラディアンベースのアルゴリズムと同じシステムサイズ依存性を持ち、従来のシステム・バスサンプラーの最良値である $\tilde{O}(n^6/\varepsilon)$ を $n^5$ の係数で改善している。
• 広範なラムシフトの制御: 並列アーキテクチャにおける広範なラムシフトを、単なる小さな摂動として扱うのではなく、制御する厳密な手法を提供する。これは、ランダムな順方向発展を利用してコヒーレント項を抑制し、ギブス中心の準局所演算子の和に対する新しいモーメント境界を確立することによって行われる。
• ノルム特異的な解析: 解析はKMS重み付き $L_2$ ノルム（$\chi^2$ ダイバージレンスに関連）において行われる。これは、システム・バスサンプラーの従来のトレースノルム解析では達成できなかった、誤差境界における $n$ への線形依存性を維持するために不可欠である。

結果
固定された逆温度 $\beta > 0$ における、有界な近接相互作用を持つ $n$ クディットの鎖に対し、本アルゴリズムは任意の初期状態からトレースノルム誤差 $\varepsilon$ 以内で、期待総ハミルトニアン発展時間：
$$T_{total} = \tilde{O}\left(\frac{n}{\varepsilon}\right)$$
を用いてギブス状態 $\rho_\beta$ を準備する。具体的には、必要なラウンド数は $m = \tilde{O}(n/\varepsilon)$ であり、各ラウンドの期待持続時間は $\sqrt{\log(n/\varepsilon)}$ に比例する。パラメータは $\alpha^2 \propto \varepsilon/n$ と設定される。

意義
本論文は、システム・バス相互作用アルゴリズムが（結合ハミルトニアン発展とリセット可能な量子ビットのみを必要とするため）近未来の量子ハードウェアへの実装が本質的に容易でありながら、複雑なリンドブラディアンベースの回路と同等の漸近的なシステムサイズ効率を実現できることを示している。並列相互作用と広範なラムシフトに関連する技術的課題を解決することで、本研究は、理論的な熱化の効率性と量子ハードウェアの実際的な制約との間の溝を埋めるものである。本結果は、より単純な実装のシステム・バスサンプラーが、スケーリング性能を犠牲にする必要はないことを裏付けている。