Estimating applied potentials in cold atom lattice simulators

技術要約：冷却原子格子シミュレータにおける印加ポテンシャルの推定

問題提起
光学格子中の冷却原子は、量子ガス顕微鏡（QGM）によるサイト分解能での読み出しが可能なハバードモデルを実現できる、多才な量子シミュレーション・プラットフォームとして機能する。重要な目的は、量子シミュレーションや情報処理への応用を可能にするため、任意のサイト依存的な光学ポテンシャルを実装することである。しかし、光学格子の回折限界（格子間隔が通常、回折限界を下回る）のため、これらのポテンシャルの精密な適用と較正は依然として困難である。

実際には、目標とするポテンシャル $V$ を適用しようとする実験者が、実際にはわずかに異なるポテンシャル $W$ を実装してしまうことがある。この差 $V - W$ は、変動するノイズではなく系統的なバイアスを構成する。この不一致は、時間の経過とともにシミュレーション誤差を蓄積させ、高精度な量子シミュレーションを妨げる。ハミルトニアンをダイナミクスから学習することは解決策となるが、既存の提案の多くは、複雑な状態準備、複数の局所基底での測定、あるいは多体ダイナミクスの古典的シミュレーションを必要とする。これらは、測定が通常、固定された粒子数基底（粒子の位置）に限定されている現在の極低温原子セットアップにおいては、実装が困難である。

手法
著者らは、粒子数読み出しのみを用いて、印加されたポテンシャル（およびホッピング率）を推定するための、シンプルかつ効率的な実験プロトコルを提案している。このプロトコルは、以下の2つの主要な物理的および計算的要素に基づいている：

1. 相互作用の抑制： このプロトコルは、特定の原子種がフェッシュバッハ共鳴を介して相互作用をオフにできる能力を利用し、系を非相互作用（自由フェルミオン）モデルへと減少させる。これにより、系の時間発展を古典的に効率よく計算することが可能になる。
2. データ収集： 実験では、既知の初期状態（具体的には電荷密度波）を準備し、未知のハミルトニアンの下でシステムを様々な時間 $t_j$ だけ進化させ、QGMを用いてサイト占有率のスナップショットを収集する。

著者らは、時間分解された局所密度 $D_i(t) = \langle c^\dagger_i(t)c_i(t) \rangle$ からポテンシャルを再構成するために、2つのアプローチを提示している：

• 厳密な解析的再構成： この手法は、相関行列の $t=0$ における二階微分を用いてハミルトニアンのパラメータを推論するために、測定データに多項式を適合させる。理論的には妥当であるが、数値解析によれば、この手法は統計的ノイズに対して非常に敏感であり、正確な再構成には膨大な数のスナップショット（$R \sim 10^7$）を必要とする。
• 効率的なヒューリスティック・プロトコル： ノイズへの敏感さを克服するため、著者らは勾配ベースの最適化を介してコスト関数の最小化を提案している。コスト関数 $K(v; C(0))$ は、試行ポテンシャル $v$ に基づく予測占有数と測定データの間の二乗差を測定する。このアプローチは、高次の微分計算を回避し、標準的な勾配降下アルゴリズム（例：L-BFGS）を利用する。

主な貢献と結果
本論文は、ヒューリスック・プロトコルが堅牢、スケーラブル、かつ実験的に実現可能であることを示している。数値実験からの主な知見は以下の通りである：

• 最適パラメータ： 再構成誤差は、最大進化時間 $t_{max} \approx 5$ （ホッピング時間の単位）および少数のサンプリング時間点 $S \gtrsim 5$ で飽和する。これらのパラメータをさらに増やしても、改善は無視できる程度である。
• スケーラビリティ： 平均再構成誤差は $\epsilon_{MRE} \sim 1/\sqrt{M}$ （$M$ は総スナップショット数）でスケールし、システムサイズ $N$ にはほとんど依存しない。最大100サイトのシステムにおいて、約 $3 \times 10^4$ 枚のスナップショットで5%以下の平均再構成誤差を達成できる。
• 誤差に対する堅牢性： プロトコルは、状態準備エラーやホッピング率の不確実性に対して堅牢である。ホッピング率が一致しない場合でも、オンサイトポテンシャルとホッピング率の両方について同時最適化を行うことで、結果として生じるバイアスを効果的に軽減できる。
• 初期化に対する堅牢性： 最適化は、初期の推測値が真のポテンシャルと大きく異なっていても（初期化の不確かさ $\lambda \lesssim 1$）、真のポテンシャルへと収束する。大きな偏差がある場合は複数の初期推測が必要になる場合があるが、手法が失敗することはない。

意義と主張
著者らは、このプロトコルが、ハミルトニアン較正のためのスケーラブルなツールを提供することで、任意のポテンシャルを用いた精密な量子シミュレーションへの道を開くと主張している。本手法の特筆すべき点は以下の通りである：

• 実験的実現可能性： 追加的な実験的変更（局所制御など）を必要とせず、標準的なQGM読み出しとグローバルなクエンチを利用する。
• 古典的効率性： 非相互作用のダイナミクスを利用することで、後処理の計算を効率化し、「次元の呪い」を回避している。
• 低いサンプル複雑性： プロトコルは、管理可能な数の実験ショットで高い精度を達成しており、多体局在やランダム化プロトコルのベンチマークに適している。

本論文は、ガウス型自由フェルミオンに対する真の逆写像を数値的に学習する、ロバストでアナログ互換性のあるシャドウ・トモグラフィ・スキームとして、この研究を位置づけている。現在の実装は非相互作用系と静的なポテンシャルに焦点を当てているが、著者らは、相互作用系や時間依存ポテンシャルへの拡張を今後の課題として挙げている。また、ランダム化測定に触発された「シャドウ較正」手法の有用な可能性についても示唆しているが、現在の研究ではこれらの拡張の実装には至っていない。