Neutrino thermalization via randomization on a quantum processor

技術要約：量子プロセッサにおけるランダム化によるニュートリノの熱化

問題提起
超新星爆発の内核のような天体物理学的環境におけるニュートリノ・フレーバーの動的進化は、コヒーレントな前方散乱（$\nu-\nu$ 相互作用）に由来するランダムな結合を持つ、全対全（all-to-all）のスピン・ハミルトニアンによって支配されている。「高速」フレーバー振動の極限において、この系は熱化することが期待されており、その挙動は固有状態熱化仮説（ETH）によって予測され、小規模な古典シミュレーションによって確認されている。しかし、この2局所（two-local）ハミルトニアンの動力学を古典コンピュータで正確にシミュレートすることは困難である。なぜなら、エンタングルメントの急速な増大により、制御可能な精度を得るためのトロッター化（Trotterization）に必要な回路深さが、少なくともシステムサイズに対して線形に増加するためである。この制限により、熱化時間のスケーリングを決定的に判断できるほど大きな系を研究することが妨げられている。さらに、現在のノイズのある中規模量子（NISQ）デバイスは、数十個の量子ビットを超えるシステムに対して、厳密なトロッター化シミュレーションを実行できるほどの深さを備えていない。

手法
著者らは、ニュートリノガスの非局所的な動力学を模倣するために、ランダム化量子シミュレーション戦略を提案している。これは、システムサイズに依存しない回路深さで熱的挙動を再現することを目的としている。

1. 物理モデル： 系は、ニュートリノのランダムな速度ベクトルに依存する結合を持つ、全対全のハイゼンベルク型ハミルトニアンによってモデル化されている。著者らは、置換不変な観測量、特にフレーバー分散を、熱化のプロキシ（代理指標）として焦点を当てている。
2. ランダム化チャネル（プロトコルB）： $O(N)$ 層の進化ユニタリを必要とする完全なトロッター化された進化を実装する代わりに、著者らはランダム量子回路を構築する。このプロトコルでは、進化の各層における結合定数が、すべての可能な結合の集合から一様にランダムに選択される。層の数 $L$ は、システムサイズ $N$ から切り離されている。
3. 時間再スケーリング： ランダム回路の離散的な層を物理的な時間 $t$ にマッピングするために、再スケーリング係数 $\alpha$ が導出される。この係数は、初期のダイナミクスの曲率を保持するように、ランダム化された分散と厳密な分散の短時間展開を一致させることで決定される。
4. エラー抑制： IBMの超伝導量子ハードウェアを用いた実験では、ノイズ再正規化アプローチを採用している。ノイズをデポラリジング・チャネルとしてモデル化し、ダイナミクスに対して不変な初期状態を用いてノイズパラメータを推定することで、著者らはポストプロセッシングにおけるノイズを含む期待値を補正している。
5. 古典的ベンチマーク： ランダム化アプローチは、以下の手法と比較してベンチマークされている：
   • 厳密なトロッター化： 小規模な系（$N \le 22$）に対して。
   • 行列積状態（MPS）： エンタングルメント・エントロピーが時間とともに線形に増大するため、MPSがこの問題には不適切であることを示している（ボンド次元が指数関数的にスケールする必要があるため）。
   • 位相空間近似（PSA）： より大きな系で使用される半古典的な手法であり、経験的なベンチマークとして用いられる。
   • パウリ伝播（PP）： 注目領域において、必要な精度を得るには不十分であることが判明した。

主な結果

• 熱化のスケーリング： 量子ハードウェア上で最大 $N=112$ 量子ビットの系をシミュレートし、より大きなスケールに対してPSAを使用することで、著者らは熱化時間 $\tau_{th}$（フレーバー分散が閾値を下回るまでの時間と定義）が、システムサイズの平方根（$\sqrt{N}$）に近似してスケールすることを見出した。
• 手法の検証： ランダム化プロトコル（プロトコルB）は、典型的なランダム・ハミルトニアンのダイナミクスを中間時間において正常に再現しており、平均的なハミルトニアンのダイナミクスに収束してしまうプロトコルAよりも優れた性能を示した。
• ハードウェアの性能： エラー抑制を用いたIBM量子プロセッサ（ibm_fez および ibm_aachen）の結果は、PSAの結果および $\sqrt{N}$ のスケーリング則と整合している。エラー抑制がない場合、結果は熱化時間を系統的に過小評価する。
• 古典的な限界： 本研究は、線形なエンタングルメント増大のため、テンソルネットワーク手法（MPS）がこの特定の問題には適していないことを確認した。一方で、半古典的なPSAは依然として堅牢な経験的ツールである。

意義と主張
本論文は、厳密な古典シミュレーションが困難な領域において、近未来の量子デバイスが「経験的な検証器」として機能し得ることを示すものであると主張している。具体的には：

• 半古典的手法の検証： ランダム化量子シミュレーション、PSA、および $\sqrt{N}$ スケーリングの一致は、ニュートリノ・フレーバー熱化の半古典的な取り扱いの妥当性を強力に支持する証拠となる。
• ランダム回路の新しい応用： 本研究は、ランダム量子回路（RQC）の新しい応用を強調している。それは単に量子優位性を示すためだけでなく、古典的にシミュレートが困難な複雑な多体系のダイナミクスに対して物理的な洞察を提供するために使用されている。
• 量子ユーティリティ： 結果は「量子ユーティリティ」のパラダイムに適合しており、そこではNISQデバイスが、経験的な古典的手法をテスト、ベンチマーク、および洗練するために使用される。これにより、将来的にはより大きなシステムや、より複雑なシナリオ（例：3フレーバーの場合や前方散乱ではない場合）の探索が可能になる。

著者らは、ランダム化チャネルに関する厳密な解析的保証は現在の研究範囲を超えていること、および本手法が現在は特定の置換不変な観測量に限定されていることを認め、控えめな姿勢を保っている。しかし、本研究は、他の方法ではアクセスできない多体系におけるスケーリング則を調査するために、量子ハードウェアを使用するための経路を確立している。