Time-resolved ARPES in pumped excitonic systems: Floquet physics induced by excitonic fields

技術要約：励起されたエキシトニック系における時間分解ARPES

問題提起
本論文は、相互作用する電子系、特にポンプ励起されたエキシトニック系において、時間分解角分解光電子分光（TR-ARPES）を用いて超高速ダイナミクスを理論的にモデル化するという課題に取り組んでいる。実験手法はフェムト秒およびサブフェムト秒のダイナミクスを捉えるまでに進歩しているが、既存の理論的アプローチには限界がある。時間依存密度汎関数理論（TD-DFT）は計算コストが高く、物理的な洞察に欠けることが多い。一方で、標準的なモデル・ハミルトニアンによるアプローチは、過度な単純化に陥る傾向がある。さらに、固体におけるフローケ物理（光誘起によるバンド混成）の実験的実現には大きな障壁が存在するが、近年の提案では、ボゾン励起、特にエキシトン（励起子）がフローケ的な現象を駆動できることが示唆されている。著者らは、これらの系におけるTR-ARPES信号をシミュレートするための、効率的かつ正確な理論的枠組みを開発することを目指している。具体的には、外部ポンプレーザーによって誘起されるサイドバンドと、内部のエキシトニック場によって誘起されるサイドバンドを区別することを目的としている。

手法
著者らは、現実的なマルチバンド系における超高速現象の研究のために以前開発された、演算子ベースの形式論である**動的射影演算子アプローチ（DPOA）**を採用している。この手法は以下のステップで構成される：

1. ハミルトニアンの定式化： 系は、価電子帯と伝導帯を持つ二次元二バンド半導体としてモデル化される。ハミルトニアンには、非相互作用項（$H_0$）とクーロン電子間相互作用項（$H_V$）が含まれる。
2. ハートリー＝フォック（HF）近似： クーロン相互作用は、HFデカップリングを用いた平均場レベルで扱われる。極めて重要な点として、イオン背景の除去により実効的なHF相互作用が引力的となり、エキシトニック束縛状態の形成が可能となる。HF自己エネルギー（$\Xi_{HF}$）は、一粒子密度行列（SPDM）から導出される。
3. 時間発展： 系は、時間依存ベクトルポテンシャルによって記述されるポンプパルスにさらされる。完全な時間依存ハミルトニアンには、非相互作用バンド、ポンプ場（最小結合およびパイエルス置換を介して）、および瞬時的なHFポテンシャルが含まれる。生成・消滅演算子の発展は、発展行列 $P(k,t)$ の運動方程式を用い、四次ルンゲ＝クッタ法を用いて解かれる。
4. TR-ARPESの計算： TR-ARPES信号は、非平衡グリーン関数を用いて計算される。著者らは、電子信号（$I_e$）と正孔信号（$I_h$）の両方を定義している。計算コストを削減するため、二重積分の直接計算を回避する補助行列 $Q(k, \omega, t_{pr})$ を利用しており、これにより信号の数値評価が大幅に高速化される。
5. システムパラメータ： 本研究では、特定のホッピングパラメータと、3D埋め込み格子から導出された長距離クーロン相互作用ポテンシャルを持つ典型的な2Dモデルを使用している。長距離相互作用と局所的（運動量に依存しない）相互作用の限界の両方が分析されている。

主要な貢献と結果

• 平衡相図： 著者らは、クーロン相互作用の強さ（$V$）および遮蔽因子（$\kappa$）の関数として平衡相図をマッピングした。その結果、$V$が増加するか$\kappa$が減少するにつれて、半導体相（エキシトニック秩序パラメータ $\tilde{\Delta}_\Gamma = 0$）からエキシトニック絶縁体相（$\tilde{\Delta}_\Gamma \neq 0$）への転移を特定した。
• エキシトン周波数： 半導体相において、著者らは弱いインパルス・ポンプを適用することで平衡エキシトン周波数（$\hbar\omega_{exc}$）を決定した。$\hbar\omega_{exc}$は、相互作用の強さが増大するか、より長距離になるにつれて減少することが判明した。
• エキシトン場誘起フローケ・サイドバンド：
  • 半導体相： ポンプ周波数がエキシトンモードと共鳴する場合、ポンプパルスが停止した後もエキシトニック秩序パラメータのコヒーレントな振動が持続する。これらの振動は、ポンプ後のTR-AR配向スペクトルにおいて明確なフローケ・サイドバンドを生成する。これらのサイドバンドは、ポンプレーザー場自体によって引き起こされるものとは明確に異なる。
  • 強度依存性： ポンプ強度が弱い場合、エキシトン周波数は平衡値に近い状態を維持する。中間および強い強度の場合、ポンプは電子分布を変化させ、電子と正孔の間の実効的な引力を変化させる。これにより、エキシトン周波数のシフトと、価電子帯における「メキシカンハット」型の分散形状の出現がもたらされる。
  • 運動量依存性： エキシトン場誘起サイドバンドは$\Gamma$点付近で最も強く現れる。著者らは、これが単に相互作用の運動量依存性に起因するのではなく、主にエキシトン周波数が$\Gamma$における局所的なバンドギャップに近いこと（小さなデチューニング条件）に起ころうとしていることを明らかにしている。
  • バンド共鳴サイドバンド： 著者らは、真のエキシトン場サイドバンドと「バンド共鳴誘起サイドバンド」を区別している。後者は、ポンプ周波数が特定の運動量におけるバンドギャップと共鳴した際に発生し、ポンプ周波数で振動する残留コヒーレンスを残すものである。
• エキシトニック絶縁体相： 平衡エキシトニック絶縁体相においては、ポンプパルスはエキシトニック秩序を部分的に抑制（エキシトンの融解）するが、コヒーレントに振動する成分を残す。これにより、価電子帯と伝導帯の両方に平行なサイドバンドが生じ、各バンドは非相互作用時の位置へとシフトする。
• 局所相互作用限界： 本研究は、観察された現象が局所的なクーロン相互作用（運動量に依存しない）の限界においても持続することを確認しており、サイドバンドの運動量依存性が相互作用の範囲ではなく、デチューニング条件によって駆動されていることを補強している。

意義と主張
本論文は、DPOAフレームワークが、相互作用する電子系における超高速現象をシミュレートするための効率的かつ正確な手法であり、微視的な観測量へのアクセスを維持しつつ、ab initio 時間依存GW法に対する計算コストの低い代替手段を提供すると主張している。

本研究の主な意義は、以下の点にある：

1. 実験的観察の裏付け： 得られた結果は、単層遷移金属ダイカルコゲナイドにおけるエキシトン駆動のフローケ的サイドバンドの観測に関する最近の実験結果（特に文献[82]）を支持している。
2. 明確なシグネチャーの特定： 本研究は、外部ポンプレーザー場に由来するサイドバンドと、内部のエキシトニック場によって誘起されるサイドバンドとの間の明確な理論的区別を確立しており、エキシトニック・コヒーレンスの「直接的な指紋」を提供している。
3. 物理的描像の検証： エキシトニック束縛状態の形成とその非平衡下でのダイナミクスの進化を成功裏にモデル化することで、本研究はエキシトン場誘起フローケ物理の物理的描像を検証している。
4. 解釈のためのツールの提供： このフレームワークは、ポンプされた半導体における進行中および将来のTR-ARPES実験を解釈するための基礎となり、非平衡相におけるクーロン相互作用の役割を理解するための基盤となる。

著者らは新しい相を発見したと主張しているのではなく、既存の実験データを説明し、これらの複雑な相互作用系をシミュレートするための堅牢な理論的記述を提供したのである。