Near-optimal synthesis of non-Gaussian phase gates via qubit-oscillator Rabi control

技術要約：量子ビット・振動子ラビ制御による非ガウス位相ゲートの準最適合成

問題提起
ユニバーサルな連続変数（CV）量子計算には、ガウス演算を補完するために、多項式位相ゲート $e^{iT p(\hat{X})}$ のような非ガウス演算が必要である。しかし、必要な非線形性を設計することは大きなボトルネックとなっている。ハイブリッドな量子ビット・振動子プラットフォーム（捕捉イオンや超伝導回路など）は、線形なラビ相互作用を通じて効果的な非線形ダイナミクスを実現する経路を提供するが、既存のプロトコルは数値最適化に依存することが多く、目標誤差 $\varepsilon$ が減少する際の必要な相互作用時間の厳密なスケーリング則が欠けている。具体的には、線形ラビ制御が、非周期的な多項式位相ゲートを、基本限界に迫るような誤差 $\varepsilon$ に対する対数的な（polylogarithmic）相互作用時間で合成できるかどうかが不明である。さらに、有限個のパルス列によってこれらのゲートを正確に実現できるのか、あるいは本質的に近似的なものなのかも分かっていない。

手法
著者らは、次数 $R \ge 2$ の実多項式 $p(x)$ に対する多項式位相ゲート $U_p(T) = e^{iT p(\hat{X})}$ のための、効率的かつ解析的にコンパイル可能な合成スキームを開発した。このアプローチは、補助量子ビットがラビパルス $X_{\omega}(t) = \exp[it\sigma_{\omega} \otimes \hat{X}/2]$ を介して振動子の量子理 $\hat{X}$ に線形に結合するハイブリッド系を利用している。

1. 厳密な合成の不可能性： 本論文では、同一の量子化（quadrature）を持つ有限個のラビパルス列では、非線形な多項式位相ゲートを正確に実現できないことを最初に証明している。これは、任意の有限シーケンスのポストセレクション関数が定数係数の線形微分方程式を満たすのに対し、対象となる関数 $e^{iT p(x)}$ （非線形な $p$ の場合）はそれを満たさないことを示すことで証明される。
2. 一般化量子信号処理 (GQSP)： 高忠実度な近似を実現するために、著者らはGQSPを採用している。彼らは、有限の窓 $[-L, L]$ 上でターゲットとなる位相関数 $e^{iT p(x)}$ を近似する有界なローラン多項式 $P_N(z)$ を構築する。
3. 窓付き近似： 固定された有限フォック部分空間 $\mathcal{H}_n = \text{span}\{|0\rangle, \dots, |n\rangle\}$ 内の状態における一様ガウス裾境界を利用することで、窓の外側での近似誤差が指数関数的に小さいことが示される。窓のサイズ $L$ は $\Theta(\sqrt{n + 1 + \log(1/\varepsilon)})$ として選択される。
4. 解析的構成： パルスのパラメータ（軸および持続時間）は、変分最適化を行うことなく、ローラン係数から解析的に導出される。このスキームは、GQSP再帰法を通じてローラン多項式を実現するために、等しい持続時間を持つ偶数個のパルスを使用する。

主要な貢献と結果

• 対数的なスケーリング（Polylogarithmic Scaling）： 著者らは、次数 $R$ の位相ゲートに対して、目標誤差 $\varepsilon$ を達成するために必要な総ラビ相互作用時間 $T_{\text{tot}}$ が以下のようにスケールすることを証明した：
  $$T_{\text{tot}} = O\left(\log^{(R-1)/2 + o(1)}(1/\varepsilon)\right)$$
  これは、フーリエ・トロッター法のような多項式スケーリング（例：$O(1/\varepsilon)$）に対する指数関数的な改善を表している。
• 準最適性： 同一量子化ラビシーケンスにおいて、$T_{\text{tot}} = \Omega\left(\log^{(R-1)/2}(1/\varepsilon)\right)$ であるという下限が確立された。これは、提案された合成が準最適であり、基本限界とは $\log^{o(1)}(1/\varepsilon)$ という劣多項式因子のみで異なることを示している。
• ユニバーサル・ゲートセット： このスキームは、多項式位相ゲートの有限の合成へと拡張されている。著者らは、対数的なスケーリングが合成の下でも保持されることを示しており、これにより、容易に準備可能な入力状態のクラスに対して、標準的なユニバーサルCVゲートセット（回転、立方位相、および制御位相ゲートを含む）を効率的にシミュレートすることが可能になる。
• 高い成功確率： ポストセレクションの成功確率は少なくとも $(1-\varepsilon)^2$ であることが保証されており、誤差が減少するにつれて1に近づく。
• 応用：
  • CVダイナミクス： 本スキームは、立方位相ゲート ($e^{i\hat{X}^3}$)、回転 ($e^{i(\hat{X}^2+\hat{P}^2)}$)、およびダブルウェル進化 ($e^{i(\hat{X}^4-\hat{X}^2)}$) について数値的に検証され、有限のリソース収束性が示された。
  • 線形偏微分方程式（PDE）の解法： 本手法は、正則化された双調和方程式（線形楕円型PDE）の解決に適用された。2次元の例（$D=2$）では、偏極分解により、自明な寄与を除去した後で38個の非ゼロ単項式位相因子が得られた。すべての因子は可換であるため、同じ線形量子化に関連する項は単一の多項式位相ブロックにまとめられ、その結果、28個の多項式位相ブロックを用いた数値実装が行われた。著者らは、フーリエ・スペクトル参照値に対して忠実度 $F_{\text{sol}} = 0.9998$ で解をシミュレートしており、これはネイティブな4次または6次の振動子非線形性を必要としない。

意義
本論文は、量子ビット・振動子のラビ制御が、効率的で解析的にコンパイル可能であり、時間計算量の観点から準最適であるCV量子情報処理のプリミティブであることを確立している。線形相互作用を、精度に対する対数的なリソーススケーリングで正確な非線形変換へとコンパイルする体系的な方法を提供することで、本研究はユニバーサルなCV量子計算を実現するための主要なボトルネックに対処している。結果は、ハイブリッド・プラットフォームが、困難なエンジニアリングを要するネイティブな非線形性を必要とせずに、複雑なCVダイナミクスを効果的にシミュレートし、アルゴリズム（線形方程式ソルバーなど）を実行できることを示唆している。著者らは、上界と下界の間の残差 $\log^{o(1)}(1/\varepsilon)$ のギャップを埋めることが、今後の研究課題であると述べている。