The Sample Complexity of Fidelity Estimation to a Known Rank-$r$ Reference State Is $\widetilde{\Theta}(r^2/\varepsilon^2)$

技術要約：既知のランク $r$ の参照状態に対する忠実度推定のサンプル複雑性

問題設定
本論文は、未知の量子状態 $\rho$ と既知のランク $r$ の参照状態 $\sigma$ との間のルート・ウルマン忠実度 $F(\rho, \sigma) = \text{tr}\sqrt{\sqrt{\sigma}\rho\sqrt{\sigma}}$ を推定するためのサンプル複雑性 $S(r, \varepsilon)$ を決定するという根本的な問いに取り組んでいる。先行研究では、$\Omega(r/\varepsilon^2)$ の下界と $O(r^2/\varepsilon^2)$ の上界が確立されていたが、因子 $r$ のギャップが未解決のまま残されていた。著者らは、忠実度のサンプル複雑性が（対数因子を除いて）本質的にランク $r$ に対して二次的であるのか、それとも線形の低界がタイトなのかを調査している。

手法
証明には、下界を確立するために「二重事前分布法（two-prior method）」を用いている。著者らは、以下の性質を持つ、量子状態に対する2つの異なる事前分布 $\Pi_0$ および $\Pi_1$ を構築する：

1. 参照状態 $\sigma$ に対する忠実度の値が一定のギャップで分離されている（すなわち、$\rho \sim \Pi_0$ に対して $F(\rho, \sigma) \le L$ かつ $\rho \sim \Pi_1$ に対して $F(\rho, \sigma) \ge H$ であり、$H - L > 2\varepsilon$）。
2. これらの事前分布に対応する $n$ コピーの平均状態が、主張された下界に達するまでの $n$ に対して、トレース距離において区別不能である（具体的には、$d_{\text{tr}}(\mathbb{E}_{\Pi_0}[\rho^{\otimes n}], \mathbb{E}_{\Pi_1}[\rho^{\otimes n}]) = o(1)$）。

この構成は、以下の3つの主要な技術的要素に基づいている：

• 厳密なスペクトルモーメント一致（Exact Spectral Moment Matching）: 著者らは、$r$ 次元上の2つの確率分布（スペクトル） $a^{(0)}$ と $a^{(1)}$ を、次数 $K \approx \Theta(\log r / \log \log r)$ まで同一のモーメントを共有するように構築する。一方の分布は疎（サポートが $\alpha r$ 個の要素）であり、もう一方は「良好な条件（well-conditioned）」（質量が $\beta r$ 個の要素に分散）である。これは、重心キャンセル技術（barycentric cancellation technique）に続く、等重性を保証するための定量的パケット補正（quantitative packet correction）を通じて達成される。
• サイズバイアス付き二重相関ウィシャートモデル（Size-Biased Doubly Correlated Wishart Model）: 著者らは、二重相関ウィシャートモデル（$G = U D^{1/2} Z D^{1/2} V^*$）を用いて行列 $G$ に対する事前分布を定義する。ここで、$\|G\|_F^{2m}$ に比例する放射状サイズバイアス（radial size bias）を導入する。このバイアスは、$m$ コピーの平均状態における正規化因子を相殺し、結果として得られる状態のシューア・ラベル分布が、シューア多項式 $s_\lambda(a)^2$ に比例するようにする。
• コーシー恒等式と長周期推定（Cauchy Identity and Long-Cycle Estimates）: コーシー恒等式とフロベニウス公式を用いて、2つの事前分布のヘリンジャー親和性と、重み付きランダム置換モデルとの関係を導出する。低次のモーメントの一致により、置換における短周期（short cycles）からの寄与が相殺される。状態の区別不能性は、この置換モデルにおける「長周期（long cycles）」の確率をバウンドすることに帰着され、これはコピー数 $m$ が二次的な閾値を下回る場合に無視できることが示される。

可換な状態から一般的な非可換なケースへと拡張するために、著者らは**直和埋め込み（direct-sum embedding）**を利用している。基本となる状態をより大きなシステム $T \oplus B$ に埋め込み、既知の参照状態を $\sigma = \frac{1}{r} P_T$ とする。$\sigma$ と非可換でありながら忠実度のギャップを維持する状態 $\omega_{G,q}$ を構築することで、可換性の仮定なしに下界が保持されることを保証する。最後に、**二項間引き（binomial thinning）**を用いて、コピー数を基本となる区別不能性の領域から、目的の $1/\varepsilon^2$ 依存性へとスケールさせる。

主な貢献と結果

1. 忠実度推定のギャップの解決: 本論文は、既知のランク $r$ の参照状態に対して、サンプル複雑性が $S(r, \varepsilon) = \tilde{\Theta}(r^2/\varepsilon^2)$ であることを証明している。具体的には、$\Omega\left(\frac{r^2}{\varepsilon^2 (\log r)^C}\right)$ の下界を確立し、Wangによる先行研究に残された因子 $r$ のギャップを閉じた。これは、固定された参照状態 $\sigma = \frac{1}{r} P_T$（$2r$ 次元のシステム上）および $\sigma$ と可換でない困難な状態族に対しても成立する。
2. スペクトル推定における近二次障壁: 同様の枠組みを適用することで、量子スペクトル（または正規化された平方根トレース $\text{Tr}(\rho) = \frac{1}{\sqrt{r}}\text{tr}\sqrt{\rho}$）を一定の精度で推定するための、$\tilde{\Omega}(r^2)$ という近二次的な下界を導出している。この結果は、Pelecanos–Spilecki–Tang–Wright による $O(r^2 (\frac{\log \log r}{\log r})^2)$ という最近の上界と、対数因子を除いて一致しており、この領域におけるサンプル複雑性の多項式次数を事実上決定している。
3. 両側ランク判定（Two-Sided Rank Testing）: これらの手法により、ランク $\le \alpha r$ の状態と、ランク $\alpha r$ の集合からの距離が大きい状態とを判別するために、$\tilde{\Omega}(r^2)$ コピーを必要とするという下界が得られる。

意義
本論文は、参照状態のランクに対する忠実度推定のサンプル複雑性の依存関係に関する未解決問題を解決したと主張している。忠実度推定およびスペクトル推定の両方において、複雑性が（対数因子を除いて）本質的に $r$ に対して二次的であることを示すことで、「近二次的な障壁（near-quadratic barrier）」を確立した。本研究は、サイズバイアス付きシューア測度を用いた厳密なモーメント一致と、重み付き置換モデルにおける長周期推定の組み合わせといった新しいツールを導入しており、これらは他の非多項式的なスペクトル関数を分析する際にも応用可能である。これらの結果は、無制限の集団測定（collective measurements）の下でも成立しており、この二次的なスケーリングが制限された測定戦略によるアーティファクトではなく、問題の本質的な限界であることを裏付けている。