Photogeneration and signatures of coherent phonons in time-resolved photoemission spectroscopy: First-principles time-dependent adiabatic GW approach

技術要約：時間分解光電子分光におけるフォトジェネレーションとコヒーレントフォノンのシグネチャー

問題提起
光パルスによって励起されるコヒーレントな格子ダイナミクス（フォノン）は、電子状態や材料特性を変調するメカニズムを提供し、量子制御やセンサーへの応用が期待されている。光ポンプ・プローブ実験（例：過渡反射率）はコヒーレントフォノンを検出できるが、状態分解された電子－フォノン（e-ph）結合強度の決定に必要な解像度が不足していることが多い。時間・角度分解光電子分光（TR-ARPES）は、運動量空間における電子構造へのアクセスを可能にし、理論的には、ARPES信号のフーリエ解析を通じて状態分解されたe-ph結合を抽出することができる。しかし、コヒーレントフォノンの生成と、それに続くTR-ARPESスペクトルの変調を同時に扱い、かつ電子－電子、電子－正孔、およびe-ph相互作用を含む包括的な第一原理的理論フレームワークは、依然として大きな課題であった。静的な変位を課した計算や現象論的なモデルといった既存の手法は、電子とフォノンのダイナミクスの相互作用や、特定の生成メカニズム（インパルシブ・スティミュレイテッド・ラマン散乱［ISRS］対、変位励起コヒーレントフォノン［DECP］）を捉えることができない。

手法
著者らは、e-ph結合を含むように拡張された、アブイニシオ時間依存的断熱的GW（TD-aGW）アプローチを提示している。この手法は、顕著な励起子効果を持つ単層MoS₂をモデルとしたポンプ・プローブ設定をシミュレートする。

• ハミルトニアンとダイナミクス： 系は、ブロッホバンド基底における単一電子密度行列 $\rho_{nm}^{\mathbf{k}}(t)$ と、フォノン変位の期待値 $\langle Q_\nu(t) \rangle$ を伝播させることでモデル化される。ハミルトニアンには、平衡状態のGW準粒子項、ハートリーおよびGW自己エネルギーの変化（静的なCOHSEX近似）、外部ポンプ場への結合、および古典的フォノン近似におけるe-ph結合項が含まれる。
• 運動方程式： 電子のダイナミクスは、GW自己エネルギーを用いた時間依存シュレディンガー方程式によって支配される。フォノン変位は、駆動力を決定するe-ph結合行列要素と時間依存密度行列によって駆動される、強制調和振動子の運動方程式に従う。
• シミュレーションパラメータ： ポンプパルスは、異なる周波数（0.6 eVから1.7 eV）および強度を持つ、切断された半周期のサイン二乗関数（持続時間10 fs）としてモデル化される。TR-ARPESスペクトルは、10 fsの間隔で電子スペクトル関数 $A(\mathbf{k}, \omega, t)$ を計算することによって模倣される。本質的なダイナミクスを孤立させるために、初期段階では散逸は無視されているが、有限のデフェージング効果についても補足資料において簡潔に検討されている。

主な結果

1. コヒーレントフォノンのTR-ARPESシグネチャー： 単層MoS₂のシミュレーションにより、コヒーレントフォノンがTR-ARPESの強度とエネルギーに周期的な変調を引き起こすことが明らかになった。具体的には、差スペクトル $\Delta A(\mathbf{k}, \omega, t)$ は、時間依存の準粒子エネルギーのシフトを示す。フーリエ解析により、これらの変調が$A'_1$フォノンモードの周波数と一致することが確認された。
2. 状態分解e-ph結合の抽出： 本研究は、スペクトルシフトの振幅（$\Delta A_{amp}$）とコヒーレントフォノンの振幅（$Q_{amp}$）の間に線形関係があることを示している。極めて重要なことに、この線形依存性の傾きは、状態分解されたe-ph結合行列要素 $g_{nn'}^{\nu}(\mathbf{k}, \mathbf{q}=0)$ に直接比例する。これは、TR-ARPES実験によって、特定の電子状態に対するe-ph結合強度を定量的に決定できることを示唆している。
3. 生成メカニズムの特定： シミュレーションは、ポンプ周波数に依存して、ISRSとDECPの両方のシグネチャーを特定している。
   • DECP（変位励起）： 光学ギャップ付近（例：A励起子に共鳴する1.7 eV）のポンプ周波数では、系は有限の平均フォノン変位（$Q_{avg} \neq 0$）を示し、光励起キャリア密度に由来する非ゼロの平均駆動力を示す。これは、励起されたキャリアが平衡位置をシフトさせる力を及ぼすDECPメカニズムと一致する。
   • ISRS（インパルシブ・スティミュレイテッド・ラマン散乱）： 低いポンプ周波数（例：ギャップ以下の0.6 eV）では、平均変位と力はゼロであり、振動位相はDECPの場合に対して$\pi/2$位相シフトする。この挙動は、駆動力がインパルシブであり、残留キャリア密度を残すことなくポンプ強度に比例するISRSメカニズムと一致している。
   • 遷移： これらの領域間の遷移は急激ではなく、中間的な周波数（例：1.2 eV）では混合した特性を示す。
4. 選択則とラマン散乱： 著者らは、コヒーレントフォノン生成とラマン散乱の関係を明確にしている。両方のプロセスは、光学遷移双極子とe-ph結合要素の同一の行列積から導かれる選択則を共有している。したがって、ラマン活性なフォノンモードは一般に生成される。しかし、駆動力の周波数依存性（コヒーレントフォノン生成）は、ラマン散乱強度とは異なり、特に光学ギャップ以下およびC励起子エネルギー以上において、共通の選択則を持ちながらも異なる共鳴構造を示す。

意義と主張
本論文は、コヒーレントフォノン生成とそのTR-AR明示への影響をモデル化するために、電子－正孔相互作用（GW-BSE経由）とe-ph結合を完全に統合した、最初の第一原理的時間依存シミュレーションを提供すると主張している。

• 定量的解析： 本研究は、コヒーレントフォノンの変調を解析することにより、TR-ARPESデータから状態分解されたe-ph結合強度を直接抽出する手法を確立した。
• メカニズムの明確化： 本研究は、ISRSとDECPが密度行列の応答における同一の駆動項から生じるものの、ポンプ周波数（透明領域か不透明領域か）に応じて異なる現れ方をすることを示すことで、生成メカニズムに関する議論を解決した。
• 選択則の洞察： コヒーレントフォノン生成とラマン散乱強度は選択則を共有しているものの、そのスペクトル挙動は同一ではないことを示し、材料における光－物質相互作用に関するより微細な理解を提供した。
• 将来の有用性： 著者らは、このアプローチがTR-ARPESにおけるコヒーレントフォノンダイナミクスの解析をサポートし、「コヒーレントフォノンによるバンド構造の定量的エンジニアリング」を可能にすると述べている。ただし、実世界の応用においては、生成メカニズムの特定を困難にするデフェージングや脱占有効果を考慮する必要があるとして、具体的な実験提案については控えめな表現にとどめている。