Non-adiabatic Effects Induced by Strong Light-Matter Coupling in Cavity QED

技術要約：キャビティQEDにおける強結合光・物質相互作用に起因する非断熱効果

問題提起
光学キャビティ内の量子化された電磁場と分子との相互作用は、外部駆動場を用いずに分子特性を修飾するメカニズムとして注目されている。量子電磁力学ハートリー・フォック（QED-HF）やQED結合クラスター理論といった、分子キャビティ量子電磁力学（QED）のためのab initio電子構造手法において、大きな進展が見られる一方で、これらの枠組みは通常、ボルン・オッペンハイマー（BO）近似に依存している。この近似では、核の動力学は間接的に扱われ、非断熱効果はしばしば無視される。具体的には、電子波動関数に対する核運動エネルギー演算子の期待値を考慮する対角ボルン・オッペンハイマー補正（DBOC）は、キャビティの設定下ではほとんど注目されてこなかった。キャビティ場は電子波動関数を非摂動的に変化させ得るため、DBOC自体も変調を受けることが予想される。しかし、キャビティ誘起のDBOC変調に関する系統的なab initio調査は欠如していた。

手法
著者らは、原子および分子が光学キャビティ内に埋め込まれた場合のDBOCを評価する系統的な研究を提示している。理論的枠組みは、標準的なQED-HFおよび強結合SC-QED-HF参照状態を用いる、長さゲージにおけるパウリ・フィールツ・ハミルトニアンに基づいている。

1. 参照状態: 本研究では、2つの平均場アプローチを利用する。
   • QED-HF: 光子場がグローバルな変位パラメータによって記述されるコヒーレント状態のアンザッツを使用する。
   • SC-QED-HF: 電荷を持つ系の並進対称性を回復し、強結合領域での精度を向上させるために、軌道依存のコヒーレント状態変位を導入する。
2. 相関の取り扱い: 電子－電子相関を考慮するために、著者らは単一および二重励起を含む量子電磁力学配置間相互作用（QED-CI）の枠組み、すなわちQED-CI(2,0)を採用している。このスキームでは、電子－光子相関は平均場レベル（偏極軌道を介して）で扱われるが、電子自由度と光子自由度の間の明示的な多体もつれは省略される。
3. DBOCの評価: DBOCは数値的な差分法を用いて評価される。複雑な結合摂動方程式を解く代わりに、著者らは核配置が変位した幾何構造（$\Psi_e(R_A \pm \delta R)$）における波動関数のオーバーラップを計算することで、核運動エネルギー演算子に必要とされる二階微分を近似している。この手法は、キャビティQEDの枠組み内におけるあらゆる電子構造手法に一般化可能である。
4. 対象系: 解析は以下の多様な系を対象としている。
   • 原子: ヘリウム（He）、水素陰イオン（H$^-$）、ベリリウム（Be）。
   • 二原子分子: H$_2$、LiH、HF。
   • 多原子分子: アンモニア（NH$_3$）、ホルムアルデヒド（CH$_2$O）。
     計算は、aug-cc-pV5Z基底関数を用いたOpenMSパッケージとPySCFのインターフェースを用いて行われた。

主な結果

• 原子系: He、H$^-$、およびBeにおいて、DBOCの導入は顕著な効果をもたらす。補正の大きさは実験分解能のオーダーに達する。H$^-$については、有限質量効果の摂動的取り扱い（QED-CI(2,0)における標準）が、強結合強度（$\lambda \ge 0.05$ a.u.）において非摂動的ベンチマークからの偏差を示すことから、光一電子系はこの領域において真に非摂動的な記述を必要とする可能性を示唆している。対照的に、HeとBeの結果は既存の非摂動的データとよく一致している。
• 二原子分子: キャビティ真空はポテンシャルエネルギー面とDBOCプロファイルを再形成し、数cm$^{-1}$程度の解離エネルギーのシフトを引き起こす。
  • H$_2$およびLiH: 分子軸に平行に偏極した電場に対して、標準的なQED-HF参照状態は、原点依存性により解離限界において非物理的な挙動を示す。これは、SC-QED-HF参照状態を使用することで解決され、正しい漸近挙動と並進対称性が回復される。H$_2$における最大DBOCシフトは約0.9 cm$^{-1}$（$\lambda = 0.1$ a.u.時）であり、結合が伸びた状態で発生する。
  • HF: 重いフッ素原子核のため効果はより小さい（約0.25 cm$^{-1}$）が、現代の高分解能分光法の分解能範囲内に留まっている。
• 多原子分子:
  • NH$_3$: 傘型反転モードに対するDBOCシフトは顕著であり、$\lambda = 0.1$ a.u.において0.32 cm$^{-1}$に達する。このシフトは反転遷移周波数の数％から数十％に相当し、キャビティ誘起のDBOC変調に対するトンネル動力学の敏感さを浮き彫りにしている。
  • CH$_2$O: C=O伸縮モードは、系統的かつ非ゼロのDBOCシフト（$\lambda = 0.1$ a.u.で0.018 cm$^{-1}$）を示し、結合強度とともに増大する。

意義と主張
本論文は、有限核質量効果はキャビティQEDにおける核動力学の不可欠な構成要素であり、単なる背景摂動として扱うべきではないと主張している。本研究の主要な意義は以下の通りである。

1. 系統的な定量化: 様々な原子および分子系にわたる、キャビティ誘起DBOC変調に関する初の系統的なab initio調査を提供したこと。
2. 分光学的関連性: 解離エネルギーや振動特性におけるキャビティ誘起シフト（数cm$^{-1}$のオーダー）が、現代の高分解能実験の分解能内に収まることを示したこと。
3. 手法の堅牢性: 平行電場構成における強結合領域での非物理的な漸近挙動を回避するために、SC-QED-HF参照状態が必要であることを確立したこと。
4. 精密解析: 強結合光・物質系における高精度分光のためには、有限質量効果を、自己エネルギーや真空偏極シフトといった他の量子電磁力学補正と同等の扱いで扱う必要があると論じている。

著者らは、彼らのQED-CI(2,0)フレームワークが、光学キャビティにおける高精度な分子分光のための厳密な基礎を提供するものであると結論付けているが、現実的な実験条件下での定量的予測のためには、完全に相関した電子－ボソン波動関数（例：キャビティFCI）への拡張が今後の課題であるとも述べている。