2025-11-13T00:22:10.836390

Expectation value estimation with parametrized quantum circuits

Wu, Kong, Yan et al.
Estimating properties of quantum states, such as fidelities, molecular energies, and correlation functions, is a fundamental task in quantum information science. Due to the limitation of practical quantum devices, including limited circuit depth and connectivity, estimating even linear properties encounters high sample complexity. To address this inefficiency, we propose a framework that optimizes sample complexity for estimating the expectation value of any observable using a shallow parameterized quantum circuit. Within this framework, we introduce two decomposition algorithms, a tensor network approach and a greedy projection approach that decompose the target observable into a linear combination of multiple observables, each of which can be diagonalized with the shallow circuit. Using this decomposition, we then apply an importance sampling algorithm to estimate the expectation value of the target observable. We numerically demonstrate the performance of our algorithm by estimating the expectation values of some specific Hamiltonians and inner product of a Slater determinant with a pure state, highlighting advantages compared to some conventional methods. Additionally, we derive the fundamental lower bound for the sample complexity required to estimate a target observable using a given shallow quantum circuit, thereby enhancing our understanding of the capabilities of shallow circuits in quantum learning tasks.
academic

パラメータ化量子回路を用いた期待値推定

基本情報

  • 論文ID: 2407.19499
  • タイトル: Expectation value estimation with parametrized quantum circuits
  • 著者: Bujiao Wu, Lingyu Kong, Yuxuan Yan, Fuchuan Wei, Zhenhuan Liu
  • 分類: quant-ph（量子物理学）
  • 発表時期: 2024年7月（arXiv プレプリント、v2版は2025年10月16日に更新）
  • 論文リンク: https://arxiv.org/abs/2407.19499

要旨

量子状態の性質推定（保真度、分子エネルギー、相関関数など）は量子情報科学における基礎的課題である。実際の量子デバイスの制限（有限な回路深度と接続性）により、線形性質の推定であっても高いサンプル複雑度に直面する。この非効率性を解決するため、本論文は浅層パラメータ化量子回路を用いて任意の可観測量の期待値推定のサンプル複雑度を最適化するフレームワークを提案する。このフレームワーク内で、2つの分解アルゴリズムが導入される：テンソルネットワーク法と貪欲投影法である。これらは目標可観測量を複数の可観測量の線形結合に分解し、各項は浅層回路で対角化可能である。この分解に基づいて、重要度サンプリングアルゴリズムを適用して目標可観測量の期待値を推定する。

研究背景と動機

問題定義

量子状態の線形性質推定Tr(ρH)は量子情報科学の中核的課題である。ここでρは量子状態、Hは可観測量である。このような問題は以下の領域に広く存在する：

  1. 量子化学：分子基底状態エネルギーの計算
  2. 多体物理学：相関関数の測定
  3. 量子情報：状態保真度の評価

既存手法の限界

  1. 古典的シャドウ(Classical Shadow)プロトコル：
    • ローカルCSプロトコルのk-局所可観測量に対するサンプル複雑度はO(4^k)
    • グローバルCSプロトコルはO(1)の複雑度を達成できるが、対数深度の回路が必要
    • いずれも「測定不依存」であり、目標可観測量の事前情報を活用していない
  2. Pauli分解法：
    • Clifford回路の実装に限定される
    • 分解はPauli可観測量に限定される
    • 深い回路または高いサンプル複雑度が必要な場合がある

研究動機

既存手法は近期量子デバイス上で以下の課題に直面している：

  • 回路深度の制限
  • 接続性の制約
  • ノイズの影響
  • 可観測量情報の不十分な活用

主要な貢献

  1. 統一フレームワーク：パラメータ化量子回路を用いた線形性質推定の汎用フレームワークを提案し、既存のPauli分解プロトコルを統一した
  2. 2つの分解アルゴリズム：
    • 貪欲投影分解(GPD)：一般的なハミルトニアンに適用可能
    • テンソルネットワーク分解(TND)：コンパクトなテンソルネットワーク表現を持つハミルトニアンに適用可能
  3. 理論的下界：与えられた浅層量子回路を用いて目標可観測量を推定するのに必要なサンプル複雑度の基本的下界を導出した
  4. 数値検証：疎密ハミルトニアンとSlater行列式内積推定において、アルゴリズムの優位性を検証した

方法の詳細

タスク定義

与えられるもの：

  • 未知の量子状態ρ
  • 目標可観測量H
  • L層深度のパラメータ化量子回路U_L(θ)
  • 精度要件εと成功確率1-δ

目標：Tr(ρH)を推定し、サンプル複雑度を最小化する

全体的フレームワーク

フレームワークは古典段階と量子段階の2つに分かれている：

古典段階：目標可観測量を以下のように分解する H≈∑k=1KUL(θ(k))†ΛkUL(θ(k))H \approx \sum_{k=1}^K U_L(\theta^{(k)})^\dagger \Lambda_k U_L(\theta^{(k)}) ここでΛ_kは実対角行列である

量子段階：重要度サンプリングを用いて期待値を推定する

  • 確率p_k ∝ ||Λ_k||_2に従って項kをサンプリング
  • U_L(θ^{(k)})を実行し、計算基底で測定
  • 中央値平均法を適用して最終推定値を得る

貪欲投影分解(GPD)アルゴリズム

核心的考え方：最適近似項U_L(θ)†ΛU_L(θ)を反復的に見つける

アルゴリズムの流れ：

  1. H^{(0)} = H、k = 0で初期化
  2. ||H^{(k)}||_2 ≥ εの場合：
    • 最適化問題を解く：θ^{(k)} = argmin_θ ||U_L(θ)H^{(k)}U_L†(θ) - diagU_L(θ)H^{(k)}U_L†(θ)||_F
    • Λ_k = diagU_L(θ^{(k)})H^{(k)}U_L†(θ^{(k)})を設定
    • H^{(k+1)} = H^{(k)} - U_L†(θ^{(k)})Λ_k U_L(θ^{(k)})を更新
    • k = k + 1

複雑度分析：古典処理時間はO(poly(n)·2^{ωn})である。ここでω ≈ 2.37は行列乗算指数である

テンソルネットワーク分解(TND)アルゴリズム

適用場面：目標ハミルトニアンが効率的な行列積演算子(MPO)表現を持つ場合

最適化目標：損失関数を最小化する L=∣∣H−∑kUL(θk)†\ΛkUL(θk)∣∣F2L = ||H - \sum_k U_L(\theta_k)^\dagger \Λ_k U_L(\theta_k)||_F^2

主要技術：

  • U_L(θ_k)を深度Lのユニタリテンソルネットワークとして表現
  • Λ_kをMPO形式で表現
  • テンソルネットワーク縮約を用いて損失関数を計算
  • 勾配降下法でパラメータ{θ^{(k)}, Λ_k}を最適化

サンプル複雑度分析

上界：アルゴリズム1は T = O(||Λ||_1^2 log(1/δ)/ε_2^2) 個のサンプルを必要とする。ここで||Λ||_1はすべての||Λ_k||_2の和である

下界：パラメータ化電路U_L(θ)を用いる任意の単一コピー適応戦略は以下を必要とする T=Ω(Tr(H02)2ε2δ(H0)4n)T = Ω\left(\frac{\text{Tr}(H_0^2)^2}{\varepsilon^2 \delta(H_0) 4^n}\right) ここでH_0はHのトレースレス部分であり、δ(H_0)はH_0の到達可能状態集合上の最大期待値の二乗である

実験設定

実験シナリオ

  1. 疎ハミルトニアン基底状態エネルギー推定：8量子ビットシステム、64個の非ゼロ要素
  2. 密ハミルトニアン期待値推定：4量子ビットランダムHermitian行列
  3. Slater行列式内積推定：3量子ビットシステム、τ-Slater行列式と純粋状態の内積

比較手法

  • 古典的シャドウプロトコル：グローバルCSとローカルCS
  • Pauli分解法：Derandomized、C-LBCS、SG、Adaptive、OGMなど
  • 専門的手法：フェルミオン古典的シャドウ(FCS)

実装の詳細

  • パラメータ化ゲート：iSWAPゲート + 2つの任意単一量子ビットゲートのテンソル積
  • GPDアルゴリズム：L=4層、K=20または80個の分解項
  • TNDアルゴリズム：L=1層、K=3個の分解項

実験結果

主要な結果

疎ハミルトニアン（8量子ビット）：

  • 25848個のサンプルの下で、GPDの誤差は0.030であり、最良の比較手法OGMの0.097を大幅に上回る
  • サンプル数の増加に伴い、GPDは常に最低の誤差を維持する

密ハミルトニアン（4量子ビット）：

  • 25848個のサンプルの下で、GPDの誤差は0.046であり、最良の比較手法OGMの0.053を上回る
  • より少ないサンプル数の場合、優位性がより顕著である

Slater行列式内積（3量子ビット）：

  • GPDはすべてのサンプル数で最低の誤差を達成
  • 25848サンプル時の誤差は0.009であり、最良の比較手法は0.012

収束性分析

数値結果は以下を示している：

  1. 分解項数Kを固定した場合、Frobenius距離は回路深度Lの増加に伴い減少する
  2. 回路深度を固定した場合、Frobenius距離は分解項数Kの指数関数的に減少する

テンソルネットワーク法の性能

低ボンド次元ハミルトニアンの場合：

  • TND法は3個の分解項と1層の回路深度のみを使用
  • 18000ステップ時の誤差は0.050であり、従来の手法を上回る

関連研究

量子状態学習

  • 量子層析：量子状態の完全な再構成、複雑度は指数関数的に増加
  • シャドウ層析：状態の古典的記述を提供し、複数の性質推定をサポート

ランダム測定プロトコル

  • ローカル測定：単一量子ビットClifford群、局所可観測量に適切
  • グローバル測定：グローバルClifford群、深い回路が必要
  • 浅層回路：折衷案だが、可観測量情報をまだ十分に活用していない

Pauli分解法

  • 可観測量のPauli展開に基づく
  • Clifford回路と計算基底測定を通じて実装
  • 本論文のフレームワークはこれらの手法を統一した

結論と考察

主要な結論

  1. 提案されたフレームワークは既存の測定プロトコルを成功裏に統一し、一般的なパラメータ化回路に拡張した
  2. GPDおよびTNDアルゴリズムは多くのシナリオで既存手法を大幅に上回る
  3. 確立された理論的下界は、浅層回路が量子学習タスクにおける基本的な制限を明らかにする

限界

  1. GPDアルゴリズム：
    • 古典的最適化複雑度は依然として高い
    • 貪欲戦略は全体最適性を保証しない
    • 分解項数Kの理論的分析は困難
  2. TNDアルゴリズム：
    • 効率的なMPO表現を持つハミルトニアンにのみ適用可能
    • 追加のテンソルネットワーク最適化技術が必要
  3. 理論的下界：
    • 低ランク可観測量（保真度など）に対して十分に厳密でない可能性
    • 回路能力パラメータδ(H_0)の正確な推定に依存

今後の方向性

  1. アルゴリズムの最適化：
    • 機械学習ベースのより効率的な分解アルゴリズムの開発
    • 非貪欲な全体最適化戦略の探索
  2. 理論の完善：
    • より厳密なサンプル複雑度下界の確立
    • 分解項数Kと回路能力の関係の分析
  3. 応用の拡張：
    • 非線形性質推定への拡張
    • 量子メモリを備えたプロトコル設計
    • ハードウェア構成切り替え回数の削減

深い評価

利点

  1. 理論的貢献：
    • 複数の既存手法を統合する統一フレームワークを提供
    • 浅層回路の能力に関する理解を深める重要な理論的下界を確立
  2. 方法の革新性：
    • GPDアルゴリズムは一般的なハミルトニアンに適用可能で実用性が高い
    • TNDアルゴリズムは特定の構造に最適化され効率が高い
    • 可観測量の事前情報を十分に活用
  3. 実験の充実：
    • 複数の応用シナリオ（疎密ハミルトニアン、内積推定）をカバー
    • 複数の主流手法との比較で説得力のある結果
    • 収束性と平均性能の分析を提供

不足点

  1. スケーラビリティの問題：
    • GPDの古典的最適化複雑度は量子ビット数に対して指数関数的に増加
    • 大規模システムでの実用性は検証が必要
  2. 実験の限界：
    • 数値実験の規模は比較的小さい（最大8量子ビット）
    • 実際の量子デバイスでの検証が欠落
    • ノイズがアルゴリズム性能に与える影響を考慮していない
  3. 理論的ギャップ：
    • 上界と下界の間に大きな隔たりが存在
    • 分解項数Kの収束速度に対する厳密な理論的保証が欠落

影響力

  1. 学術的価値：
    • 量子状態学習に新しい理論的フレームワークを提供
    • 浅層量子回路の能力に関する理解を推進
    • 近期量子計算応用に実用的なツールを提供
  2. 実用的価値：
    • NISQ デバイスに適したアルゴリズム設計の理念
    • 量子化学と多体物理学での潜在的応用
    • 量子優位性検証のためのベンチマークツール

適用場面

  1. 量子化学：分子基底状態エネルギーと性質計算
  2. 量子シミュレーション：多体系の相関関数測定
  3. 量子機械学習：特徴マッピングとカーネル法
  4. 量子最適化：目的関数の期待値推定
  5. 量子誤り訂正：符号語保真度と誤り訂正性能評価

参考文献

本論文は66篇の関連文献を引用しており、量子状態学習、ランダム測定、古典的シャドウ、Pauli分解などの中核領域における重要な研究をカバーしており、研究に堅実な理論的基礎を提供している。