Variational Polaron Theory for Ground States of Strongly Coupled Light-Matter and Electron-Phonon Systems

技術要約：強結合光・物質系および電子・フォノン系の基底状態のための変分ポラロン理論

問題提起
物質の自由度（電子、スピン、または振動）とボゾンモード（光子またはフォノン）との間の強結合は、仮想的なボゾン励起によって「ドレスされた」基底状態を生成する。結合強度が特性的なボゾン周波数や物質励起エネルギーと同等、あるいはそれを超える極強結合（ultrastrong-coupling）および深強結合（deep-strong-coupling）領域では、標準的な理論的手法は破綻する。ボゾンのドレス効果がもはや微小ではないため、摂動論的な取り扱いは崩壊し、一方でボゾン・ヒルベルト空間の裸の状態による打ち切り（truncation）は非効率的または制御不能になる。逆に、従来の強結合ポラロン描像は、中間領域において定量的な正確性に欠けることが多い。既存の変換フレームワークの手法は、制限されたボゾン参照系（変換された真空やコヒーレント状態など）に依存することが多く、残差ボゾン・ヒルベルト空間を明示的にスキャンしないため、残留励起を見逃す可能性がある。さらに、多くの既存手法は、無限結合極限における変換ハミルトニアンの漸近的なデカップリング（非結合化）を、変分最適化の指針として明示的に利用していない。

手法
著者らは、状態依存型ポラロン変換に基づく非摂動的な変分基底状態フレームワークを導入している。本手法の核となる構成要素は以下の通りである：

1. 状態依存型変換： ユニタリ変換 $\hat{U}(\lambda)$ は、特定の物質状態への射影演算子 $\hat{P}_\mu$ に依存する変位振幅 $\lambda_{\alpha,\mu}$ を用いて定義される。この変換は、支配的な線形物質・ボゾン結合を基底の中に吸収する。
2. 漸近的デカップリング： 著者らは、無限結合極限において、最適化された変位パラメータが主要な対角線型線形結合を相殺することを証明している。同時に、非対角的な物質遷移は、変位の差に伴って指数関数的に減衰する変位振動子のオーバーラップによって抑制される。その結果、変換されたハミルトニアンは、ボゾンのみの項とブロック対角な物質項からなる、漸近的にデカップリングされたものとなる。
3. 変分積状態アンザッツ： 零次近似の基底状態は、変換されたフレーム内での物質セクターとボゾンセクターの間の積状態（$|\Phi_0\rangle = |\Phi_M\rangle \otimes |\Phi_B\rangle$）として構成される。決定的なのは、ボゾン成分 $|\Phi_B\rangle$ が変換された真空や独立したモードに限定されないことである。これは、非真空励起やモード相関を捉えるために、保持されたマルチモード・フォック空間上で展開される。
4. 二次の摂動補正： 有限の結合において存続する残留的な物質・ボゾンもつれを考慮するために、二次のエネルギー補正が適用される。これには、平均場零次近似ハミルトニアンを定義し、残留結合演算子 $\delta \hat{H}$ に基づいて補正を計算するプロセスが含まれる。
5. 最適化： 変分エネルギー $E_{var}$ を最小化することにより、変位パラメータ、物質係数、およびボゾン係数を同時に最適化する。

主な貢献と結果
本論文は、このフレームワークを、ディッケ・モデル（集団的な光・物質結合）とホルシュタイン・モデル（電子・フォノン結合）という二つの典型的なモデルに対してベンチマークを行っている。

• ディッケ・モデルのベンチマーク：

  • 単一エミッタ（$N=1$）の場合、零次の変分結果による相対エネルギー誤差の最大値は約 $1.2\%$ である。二次の補正を含めることで、これは $0.2\%$ 未満に減少し、フィデリティは $0.999$ を超える。
  • 本手法はシステムサイズに対して良好なスケーリングを示す。$N=100$ の場合、相対エネルギー誤差は約 $10^{-4}$ まで低下する。
  • 本フレームワークは超放射量子相転移を正常に再現しており、基底状態エネルギーの二階微分における鋭い特徴や、臨界結合における光子凝縮の開始を捉えている。
  • 結果は、「実用的な開始」としての有効デカップリング領域を定量化しており、誤差が弱結合（$g \to 0$）と強結合（$g \to \infty$）の両方の極限で消失し、中間領域でのみピークに達することを示している。

• ホルシュタイン・モデルのベンチマーク：

  • 二サイト・シングルエレクトロン・モデルに適用した場合、本手法は変分レベルで最大約 $1.4\%$ のエネルギー誤差を達成し、これは二次の摂動補正後に約 $0.5\%$ まで減少する。
  • 本研究は、並進対称性の役割を強調している。対称性の制約を緩和すると、補正前の変分エネルギーは低下するが、対称性の破れによりフィデリティが急激に低下する。対称性制約のある分岐の方が、補正前の正確な基底状態に近い状態を維持しており、これは、残留相関が顕著な場合、必ずしも「より低い補正前エネルギー」が「より優れた物理的状態」を意味するわけではないことを示唆している。

意義と主張
本論文は、このドレスされた基引数フレームワークが、弱結合、中間結合、および強結合の全領域において、強結合系をモデリングするためのコンパクトかつ定量的に信頼できる経路を提供すると主張している。その意義は、以下の具体的な側面に存在する：

• 補間： 弱結合（裸の状態）と強結合（ポラロン）の記述の間を正常に補間し、固定されたポラロン変換が最も信頼性を欠く中間領域においても精度を維持する。
• 物理的基礎： 漸近的デカップリングの特性は、単なる形式的な極限としてではなく、積状態アンザッツの物理的正当化、およびコンパクトな記述が十分であるかを判断するための診断ツールとして機能する。
• モジュール性： 零次レベルで物質セクターとボゾンセクターを因数分解する（一方で相関のあるマルチモード・ボゾン波関数を保持する）ことにより、物質セクターに対して問題固有のソルバー（例：構成間相互作用、テンソルネットワーク、または電子構造計算法）を使用できる一方で、ボゾンによるドレス効果を変分的に扱うことが可能になる。
• 非摂動的な正確性： 変分最適化と二次の摂動補正を組み合わせることで、裸の状態による打ち切りや制限された変換真空参照に頼ることなく、基底状態の非摂動的なモデリングを可能にする。

著者らは、このアプローチが、マルチモード分子ポラロンの基底状態、構造化された電子・フォノン環境、および強結合した分子や材料系の第一原理モデルにとって特に有用であると結論付けている。