技術要約：時間依存および時間独立なリンドブラディアン・シミュレーションのための高速ランダム・コンパイラ

問題設定

リンドブラディアン（Lindblad master equations）に従う開いた系の量子シミュレーションは、ノイズや設計された散逸を伴う現実的な物理シナリオをモデル化するために不可欠である。ランダム化された積公式（qDRIFTなどの）は、一次のシミュレーションにおいてリンドブラディアン動力学への拡張に成功しているが、高次のランダム化コンパイラは、主に閉じた系（ハミルトニアン）のシミュレーションに対して開発されてきた。さらに、時間依存のリンドブラディアンのシミュレーションは、非ユニタリな進化と、完全正値トレース保存（CPTP）ダイナミクスを維持する必要性 due to の追加的な課題を提示する。

本研究は、時間独立および時間依存のリンドブラディアン・シミュレーションのための、一次および二次のランダム化サンプリング・アルゴリズムの開発に取り組むものである。主な目的は、標準的な一次手法の $\varepsilon^{-1}$ 依存性を超え、目標精度（$\varepsilon$）を達成するために必要なタイムスライス数（$r$）のスケールを改善することである。

手法

1. 時間独立リンドブラディアン・シミュレーション

著者らは、リンドブラディアン $\mathcal{L}$ を局所的なハミルトニアンおよび散逸生成子に分解する、既存の qDRIFT 型のアプローチに基づいている。

• 一次（qDRIFT型）: 標準的な平均化チャネル $Q_\tau$ は、係数に従って生成子をサンプリングすることによって構築される。これは $O(1/r)$ の誤差スケーリングをもたらす。
• 二次（qSWIFT型）: ハミルトニアンに対する qSWIFT アルゴリズムに着想を得て、著者らは、厳密な進化 $e^{\tau \mathcal{L}}$ のテイラー展開と、平均化された qDRIFT チャネル $Q_\tau$ との比較により、二次補正を導出している。
  • 彼らは、厳密な進化の二次項の差を表す補正マップ $\mathcal{L}^{(2)}$ を定義する。
  • 二次マップ $Q_t^{(2)}$ は、この補正項を $r$ 個のタイムスライスのうちの1つに挿入し、可能なすべての挿入位置について和をとることによって構築される。
  • 実装上の課題: 得られる二次マップは必ずしも CPTP ではない。そのため、著者らは出力状態の準備ではなく、観測量の期待値の推定を提案している。彼らは、qSWIFT コントロール量子ビットと**ヘラルディッド・ユニタリ拡張（heralded unitary dilations）**を利用して、非物理的な補正マップを密度行列の非対角ブロック上にコヒーレントに実装し、期待値の不偏推定を可能にする。

2. 時間依存リンドブラディアン・シミュレーション

時間依存の系については、進化は時間順序指数によって記述される。

• 一次（連続時間サンプリング）: 著者らは、連続的な qDRIFT アプローチを拡張している。彼らは、ある時刻におけるリンドブラディアンのダイヤモンドノルムに比例する確率密度に従って時刻 $\tau$ をサンプリングし、対応する生成子を適用するランダム化チャネルを定義する。
• 二次補正: 厳密な時間順序プロパゲータと、ランダム化された一段階チャネルとの間の二次微分から、補正マップが導出される。時間独立の場合と同様に、この補正はタイムスライスのシーケンス内に挿入される。
• 線形結合入力モデル: フレームワークは、リンドブラディアンが局所的な時間依存生成子（ハミルトニアンおよび散逸）の和であるシナリオに拡張される。アルゴリズムは、各区間内の局所的な項のインデックスと時間を共にサンプリングし、これによりシミュレーションをフル・グローバル・リンドブラディアンではなく、局所的な進化として定式化することを可能にする。

3. 観測量の推定

二次補正されたマップは CPTP が保証されないため、本論文は観測量 $G$ の期待値の推定に焦点を当てている。

• 全体の期待値は、一次の項（標準的なサンプリングを介して推定される）と、二次補正項（$\delta g$）に分解される。
• 補正項は、コントロール量子ビットを用いた アダマール・テスト型の構成 を用いて推定される。
• Swift マップ: 非ユニタリな生成子（ハミルトニアンの交換子およびリンドブラッド散逸の両方）をコヒーレントに実装するために、著者らはシステム・コントロール状態の非対角ブロック上に作用する「swift マップ」を定義する。これらは、ヘラルディング量子ビットを伴うユニタリ拡張を通じて実現される。ヘラルディングの成功した結果はサンプル平均に寄与し、失敗はゼロとして記録され、明示的な正規化なしに不偏推定量となることを保証する。

主な結果

誤差境界と複雑性

• 時間独立の場合: 全体の強度 $\lambda$ と進化時間 $t$ を持つ時間独立リンドブラディアンに対して、二次 qSWIFT 型近似は以下の正規化ダイヤモンドノルム誤差を達成する：
  $$\frac{1}{2} \| e^{t\mathcal{L}} - Q_t^{(2)} \|_\diamond = O\left( \frac{\lambda^3 t^3 + \lambda^4 t^4}{r^2} \right)$$
  $\lambda t \geq 1$ のとき、誤差 $\varepsilon$ を達成するために必要なタイムスライス数は：
  $$r = O\left( \frac{\lambda^2 t^2}{\sqrt{\varepsilon}} \right)$$
  これは、一次の qDRIFT 法の $\varepsilon^{-1}$ 依存性から $\varepsilon^{-1/2}$ への改善を表している。

• 時間依存の場合: 一般的な時間依存リンドブラディアンに対して、二次サンプリング・アルゴリズムは、リンドブラディアンの積分ダイヤモンドノルムに依存する、正規化ダイヤモンドノルムにおける $O(1/r^2)$ の誤差スケーリングを達成する。その結果、目標精度 $\varepsilon$ に対する依存性も $\varepsilon^{-1/2}$ に改善される。

• 線形結合モデル: 結果は線形結合入力モデルに拡張され、ゲート複雑性は局所的な生成子の積分ノルムの総和によって決定される。

実装の詳細

本論文は、以下に関する具体的なアルゴリズムを提供している：

1. 導出された確率分布に従ったインデックスと時間のサンプリング。
2. 非 CPTP 補正を実現するために必要な「swift マップ」のためのユニタリ拡張の構築。
3. ヘラルディング失敗時のスケーリング因子とサンプリング確率を含む、観測量の期待値の不偏推定量の計算。

意義と主張

本論文は、時間独立および時間依存の両方の設定において、リンドブラディアン動力学のために特別に設計された最初および二次のランダム化サンプリング・アルゴリズムを提供することを主張している。

• 精度の向上: 主要な意義は、タイムスライス数に関する精度の依存性が、一次（qDRIFT）の $\varepsilon^{-1}$ から二次の $\varepsilon^{-1/2}$ へと改善されたことにある。これは、高次のハミルトニアン・シミュレーション手法（qSWIFTなど）で見られるスケーリングの改善と一致するが、これを開いた量子系に適用したものである。
• 非 CPTP マップの取り扱い: 本研究は、高次の補正が開いた系においては物理的なチャネルではないという技術的障害に対処している。状態準備から観測量推定へと目標をシフトさせ、ヘラルディッド・ユニタリ拡張を利用することで、著者らはこれらの高次補正を利用するための実用的な操作的手法を実証している。
• 汎用性: フレームワークは一般的な時間依存リンドブラディアンを受け入れ、線形結合入力モデルに拡張されており、広範なクラスの局所的な開いた量子系への適用を可能にしている。

著者らは、二次構成は厳密であるが、より高次（$p > 2$）への形式的な拡張はスケッチされているのみであり、さらなる厳密な誤差解析が必要であることを述べている。同様に、提案された推定量の統計的オーバーヘッドと分散も、将来の最適化領域として特定されている。