Finite size scaling of bitstring probability distributions for Rydberg arrays

技術要約：リドベリ・アレイにおけるビット文字列確率分布の有限サイズスケーリング

問題提起
粒子、原子核、および物性物理学の量子シミュレーションにおいて、統計的な精度で観測量を推定するには、膨大な数の測定（ショット）が必要である。シミュレーションされるシステムサイズ（$N_q$、量子ビット数）が増大するにつれ、可能なビット文字列状態の数は指数関数的に増加し、個々の状態の確率（$p_{\{n\}}$）は減少するとともに、低確率領域において密度が高くなる。ここで、個々の確率の小ささは、その膨大な数によって補完され得るのか、そして、より大きなシステムに対して所望の統計的精度を得るために必要な最小ショット数はいくらなのか、という極めて重要な問いが生じる。先行研究は、大規模な分布を通じてシステムを特徴付けることが構造的な振る舞いを明らかにできることを示唆しているが、リドベリ原子アレイにおけるビット文字列確率の明示的なスケーリング則は、まだ十分に定量化されていない。

手法
著者らは、ラダー（梯子型）構成（垂直方向と水平方向の距離比 $\rho=2$ の2列の行）を持つリドベリ原子アレイを調査している。このシステムは、ラビ周波数（$\Omega$）、デチューニング（$\Delta$）、およびファンデルワールス相互作用（$V_{ij}$）を含むハミルトニアンによって記述される。

確率分布を分析するために、著者らは密度行列繰り込み群（DMRG）法を用いた行列積状態（MPS）を利用している。彼らは、4から22段（最大44原子）までのシステムに対して、ビット文字列の確率（$p_{\{n\}}$）の辞書を生成した。分析の核心は、以下の累積確率分布である：
$$\Sigma(p_\Lambda, N_q) = \sum_{\{n\}: p_{\{n\}} \leq p_\Lambda} p_{\{n\}}$$
この関数は、閾値 $p_\Lambda$ 以下の確率を持つすべての状態の総確率を表している。著者らは、異なるシステムサイズ（$N_q$）およびショット総数（$N_{sh} = 10^9$）にわたる $\Sigma$ のスケーリング挙動を分析している。

主な貢献と結果

1. 有限サイズスケーリングと曲線崩壊（Curve Collapse）：
   著者らは、異なるシステムサイズの累積確率分布が、「繰り返される振る舞い」を示すことを実証した。特に $0.2 < \Sigma < 0.8$ の範囲において顕著である。特性的な確率 $p^*$（分布が固定の $\Sigma$ 値に達する点）を特定し、独立変数を $(p_\Lambda/p^*)^B$ （ここで $B$ は線形領域の傾き）として再スケーリングすることで、様々な $N_q$ の分布が単一の普遍的な曲線へと崩壊（collapse）することを示した。

2. フェルミ・ディラックへの類似性（領域 II）：
   分布の中央領域（領域 II）において、崩壊した曲線はフェルミ・ディラック統計分布に酷似している。著者らは、$\Sigma(p_\Lambda)$ が以下のように近似できることを示している：
   $$\Sigma \simeq \frac{1}{1 + e^{\beta(\epsilon - \mu)}}$$
   ここで、$\epsilon = -\ln(p_\Lambda)$ および $\mu = -\ln(p^*)$ である。この類似性により、「変曲点」（$\Sigma = 1/2$ となる点、これを $p^*_{1/2}$ と表記）を特徴付けることが可能になる。$p^*_{1/2}$ のシステムサイズに対するスケーリングは指数関数的であることが判明しており、これは $N_q$ が増加するにつれて、確率質量が膨大な数の状態にわたってより一様に分布していくことを示している。

3. 最大確率のスケーリング（領域 I）：
   高確率状態である領域 I において、著者らは最大単一状態確率（$p_{\max}$）を追跡している。彼らは、$p_{\max}$ がシステムサイズに対して指数関数的にスケールすること（$p_{\max} \sim 2^{-N_q}$）を確認した。これにより、与えられたショット総数に対してアクセス可能なシステムサイズに根本的な限界があることが確立された。すなわち、最小分解可能確率は $1/N_{sh}$ によって制限される。

4. 低確率および単一出現状態（領域 III および IV）：

   • 領域 III: より低い確率においては、対数－対数プロット上で分布はべき乗則の形式 $\Sigma(p_\Lambda) \simeq A p_\Lambda^B$ に従う。著者らは、ここでの振る舞いが有限のショット総数（$10^9$）に強く影響されており、それが $p \approx 10^{-9}$ でハードカットオフを生じさせていることを指摘している。
   • 領域 IV: この領域は、単一回のみ出現する状態の「テイル（裾）」を表す。著者らは、正確に1回だけ観測された状態の数を定量化し、その数がシステムサイズとともに増大するものの、全ヒルベルト空間の次元（$2^{N_q}$）よりも遅い速度で増大することを示した。

5. ショット要求量の推定：
   $p^*_{1/2}$ および $p_{\max}$ のスケーリングを分析することで、著者らは将来のシミュレーションに必要なショット要求量を推定している。

   • $p_{\max}$ のスケーリングに基づくと、約 120 量子ビット のシステムが $10^9$ ショットでシミュレート可能である。
   • $p^*_{1/2}$ のスケーリング（累積確率が 0.5 に達する点）に基づくと、特定の精度を得るための限界は、同じショット総量に対して約 60 量子ビット 付近であると推定される。

意義と主張
本論文は、より小さな計算可能なシステムによる累積確率分布を分析することで、大規模な量子システムに必要なショット要求量を予測する手法を提供すると主張している。主要な意義は、スケーリングの $p^*_{1/2}$ および $p_{\max}$ を通じて定量化された、システムサイズ増加に伴う指数関数的なコストの特定にある。

しかし、著者らは結論において控えめな姿勢をとっており、コストは指数関数的ではあるものの、システムサイズの増加率のわずかな割合で増加することを指摘している。したがって、これはデバイスの使用を「ハードキャップ（絶対的な制限）」するものではなく、むしろ「リソースの知的な利用」を促すものである。本研究は、提案された崩壊（collapse）とフェルミ・ディラックの類似性を通じて、大規模なシステムの振る舞いがより小さなシステムの分布によって記述できることを確立しており、統計誤差の推定および、ますます大規模化するリドベリ・アレイのシミュレーションの実現可能性を判断するためのツールを提供している。