Kibble--Zurek Scaling in the Dicke Model at Mesoscopic Scales

技術要約：メソスコピックスケールにおけるディッケモデルのキブル・ズurek・スケーリング

問題提起
ディッケモデルは、集団的な光・物質相互作用および超放射相転移の典型的な枠組みとして機能する。その基礎となるハミルトニアンは単純であるが、実験的にアクセス可能なメソスコピック領域（有限の粒子数 $N$）において、臨界指数を定量的に抽出することは依然として大きな課題である。二つの主要な障害が漸近的なスケーリング挙動を不明瞭にしている：

1. 相関時間の緩やかな発散： 全結合系においては、臨界相関時間が $N$ に対して緩やかに成長するため、有限サイズ補正が支配的となり、動的指数の推定値にバイアスを与える。
2. 散逸によるクロスオーバー： 回避不可能な光子の損失は、有効な摂動として作用し、システムを異なる散逸普遍クラスへと駆動する。メソスコピックなサイズにおいては、コヒーレント（閉じた）極限と散逸的（開いた）極関数の間のクロスオーバーは幅広く、見かけ上の指数はどちらの漸近的な固定点とも一致しないことが多い。
   さらに、既存の全結合モデルに対する動的スケーリング則は漸近的な大-$N$極限に対して定式化されており、強い有限サイズ・ドリフトや、サイズ・散逸・ランプ速度が絡み合ったクロスオーバー効果が解析を複雑にする有限系における、制御された抽出手法の欠如が課題となっている。

手法
著者らは、閉じたディッケモデルと開いたディッケモデルの両方に対して、静的、スペクトル的、および駆動（ランプ）ダイナミクスを等価に扱う、統一されたメソスコピック・スケーリング・フレームワークを開発した。

1. 大-$N$場理論： 著者らは、臨界指数を分析するために統一されたケルディッシュ場理論を利用している。主要な $1/N$ 補正を明示的に保持することで、以下の異なる固定点を特定した：

   • 閉じた場合 ($\kappa=0$): システムは全結合イジング普遍クラスに属し、指数は $z=1/3$ および $\nu=3/2$ である。
   • 開いた場合 ($\kappa \neq 0$): システムは Model-A 普遍クラスへと流れ、指数は $z=1/2$ および $\nu=2$ である。
   • この理論は、キブル・ズurek（KZ）スケーリング関係式 $1/\mu = 1/\nu + z$ を確立している（ここで $\mu$ はランプ速度を特徴付ける）。

2. 統一スケーリング・プロトコル： メソスコピックなデータから指数を抽出するために、著者らは主要な無関係な補正（指数 $\omega$）を明示的に考慮した二段階のプロトコルを提案している：

   • ステージ I ( $\omega$ の固定): 無次元の観測量（例：バインダー比 $U$）を臨界点において分析する。有限-$N$ ドリフトを $N^{-\omega}$ のべき級数の切断系列としてフィッティングし、無関係な指数を決定する。
   • ステージ II ( $z, \nu, \mu$ の抽出): 固定された $\omega$ を用いて、以下の項目から残りの指数を独立して抽出する：
     • $z$: 臨界ギャップ（開いたシステムではリウヴィリアン・ギャップ）のスケーリング。
     • $\nu$: 調律パラメータに関する無次元観測量の微分のサイズ依存性。
     • $\mu$: ランプダイナミクスにおける主要な速度補正（例：$\partial_s U$ または $\partial^2_s U$）のサイズ依存性。

主な結果
数値シミュレーションと大-$N$展開が行われ、フレームワークの検証がなされた：

• 閉じたディッケモデル:

  • 無関係な指数は $\omega \approx 0.372$ と決定された。
  • 抽出された指数は、理論的な大-$N$値である $z \approx 0.33$（理論値 $1/3$）、$\nu \approx 1.52$（理論値 $3/2$）、$\mu \approx 0.99$（理論値 $1$）を回収した。
  • KZ関係式 $1/\mu = 1/\nu + z$ が検証された。
  • ランプダイナミクスは、時間反転対称性の欠如により、二次的な速度依存性を示した（線形補正の不在）。

• 開いたディッケモデル:

  • 散逸率 $\kappa \in [0.1, 3]$ にわたる結合フィットにより、$\omega \approx 0.363$ が得られた。
  • $\kappa=1$ において、厳密対角化は $z \approx 0.48$ を与えた。しかし、一般的な $\kappa$ に対して、厳密対角化でアクセス可能な小さなサイズ（$N \le 24$）では、モード交差や再配置が発生し、漸近的なスケーリングを不明瞭にした。
  • これを解決するため、$N \in [10^5, 10^7]$ に対して数値的な大-$N$展開が行われ、抽出された指数が異なる $\kappa$ に対して理論値 $z=1/2$ に近づくことが確認された。
  • 指数 $\nu$ は、$\kappa=1$ で $\nu \approx 2.02$ として抽出され、異なる $\kappa$ に対して期待値の周囲で変動した。
  • 指数 $\mu$ は、$\kappa=1$ で $\approx 0.96$、異なる $\kappa$ に対して $\approx 1.02(5)$ であった。
  • KZ関係式は開いたケースでも成立した。
  • ランプダイナミクスは、時間反転対称性の破れにより、線形速度補正を示した。

• クロスオーバー・ダイナミクス:

  • 本研究は、有限サイズ（$N$）、散逸（$\kappa$）、およびランプ速度（$s$）の競合を明らかにした。
  • クロスオーバー・スケーリング変数 $\kappa N^{1/3}$ が特定され、「閉じた系に近い」振る舞いを示す小規模・弱散逸系と、「開いた系に近い」振る舞いを区別した。
  • KZタイムスケール（$t_{KZ} \sim s^{-1/3}$）と散逸タイムスケール（$t_\kappa \sim \kappa^{-1}$）の相互作用が、生のランプデータが複雑に見える理由を説明しており、これらを分離するために統一的なスケーリング処理が必要であることを示した。

意義および主張
本論文は、閉じた系および開いた系の両方における、長距離相互作用を持つメソスコピック系における静的および動的な臨界スケーリングを解決するための統一的フレームワークを確立したと主張している。その主な貢献は以下の通りである：

1. メソスコピックな曖昧さの解消: 主要な無関係な補正を取り入れることで、標準的な有限サイズ・スケーリングが失敗したり、バイアスのかかった推定値を与えたりする可能性があるメソスコピック・データから、正しい熱力学的指数を回収するプロトコルを提示した。
2. KZスケーリングの検証: 空間的なドメイン形成を伴わない完全結合モデルにおいて、キブル・ズurek・スケーリング（$1/\mu = 1/\nu + z$）の初の制御された検証を行い、コヒーレントおよび散逸の両方の領域においてこの関係を確認した。
3. 広範な適用可能性: 著者らは、このフレームワークがキャビティQEDに限定されるものではなく、トラップされたイオンやSachdev–Ye–Kitaev型モデルなど、制御された大-$N$理論が利用できず、直接的なデータ抽出が不可欠な他のメソスコピック系にも広く適用可能であると述べている。

本研究は、熱力学的極限は原理的には明確であるが、実験的なアクセスはしばしばメソスコピックなスケールに制限されることを強調しており、本プロトコルはその領域において普遍クラスを探索するための実用的な経路を提供するものである。