Coherent and Incoherent Time-Resolved Diffuse Scattering Intensity from First Principles

技術要約：第一原理に基づくコヒーレントおよびインコヒーレントな時間分解拡散散乱強度

問題提起
X線自由電子レーザー（XFEL）や超高速電子回折（UED）による時間分解拡散散乱は、光励起された固体の非平衡格子ダイナミクスへのアクセスを提供する。実験信号には、2つの明確な時間的変調が見られる。一つは、格子の加熱に伴いピコ秒オーダーで上昇する遅いインコヒーレントな背景（バックグラウンド）であり、もう一つは、フォノン周波数の2倍の周期を持つ速いコヒーレントな振動成分である。既存の理論的解釈の多くは、これらの成分を個別に扱う簡略化されたモデルに依存している。インコヒーレントな背景は、通常、マルチ温度モデル（例：アレンの二温度モデル）によって記述される。一方、コヒーレントな振動は、瞬時な周波数シフトを仮定する「突然のクエンチ（sudden-quench）」近似を用いてモデル化されることが多い。これら両方の寄与を等価に扱い、特に非平衡キャリア分布からの電子－フォノン散乱と光誘起周波数シフトの相互作用を考慮した、統一的な第一原理フレームワークが欠如している。

手法
著者らは、密度行列形式を用いた、時間依存する一次拡散強度の第一原理的な定式化を提示している。このアプローチは、プローブパルスのコヒーレンス時間が電子的な時間スケールよりも長く、イオンの運動の時間スケールよりも短いことを前提としている。

1. 理論的枠組み:

   • 散乱断面積は原子変位の項に展開される。時間依存する一次拡散強度は、変位－変位相関関数 $\langle u_{q\nu} u_{-q\nu'} \rangle_t$ に比例すると特定される。
   • 系は、時間依存する非平衡ボルン・オッペンハイマー・ポテンシャルに関連する動力学行列を、瞬時に対角化するフォノン正規モードによって記述される。
   • 変位－変位相関関数の運動方程式（EOM）が導出される。この階層は、4点相関関数のクラスター展開に続くマルコフ近似によって閉じられる。
   • この形式は、インコヒーレントな寄与（フォノン占有数 $n_{q\nu}$ によって支配される）と、コヒーレントな寄与（異常フォノン相関 $\langle \hat{a}_{q\nu} \hat{a}_{-q\nu} \rangle$ によって支配される）を区別する。

2. 主要な方程式:

   • 異常相関関数のEOM（式11）は、コヒーレントなダイナミクスが2つのソースによって駆動されることを明らかにしている。
     1. 光誘起されたフォノン周波数シフトに由来する、「スクイージング駆動（squeezing drive）」項 $\dot{\omega}_{q\nu}/\omega_{q\nu}$。
     2. 非平衡キャリア分布からの散乱に由来する電子－フォノンソース項。これは、格子不安定性が存在しない場合でも異常相関を生成し得る。
   • 全拡散強度（式12）は、インコヒーレント項（$2n_{q\nu} + 1$）とコヒーレント項（$2\text{Re}\langle \hat{a}_{q\nu} \hat{a}_{-q\nu} \rangle$）の和として表される。

3. 実装:

   • このスキームは、電子－フォノン行列要素の逆空間積分を効率的に評価するためにワニエ補間を用いて、EPIq コードに実装されている。
   • 初期状態および励起状態のフォノンスペクトルと変形ポテンシャルを得るために、2つの独立した第一原理計算（基底状態および光励起状態）が行われる。
   • 電子およびフォノンの占有数の時間発展は、スイッチング関数と時間依存格子温度を用いて、フェノメノロジー的にモデル化されている（電子占有数は約200 fs以内に準平衡状態へと遷移する）。

主な結果
本フレームワークは、光励起されたゲルマニウム（Ge）に適用され、フェムト秒拡散散乱測定（文献7）と比較された。

• 実験的特徴の再現: 計算された拡散強度は、以下の特徴を再現することに成功した。
  • 振動周波数: これは光励起されたフォノン枝の $2\omega_{q\nu}$ に対応する。
  • 振動振幅の波長依存性。
  • 上昇するインコヒーレントな背景: これは、エネルギーがキャリアから格子へと転移する際のフォノン占有数の増加から自然に生じる。この特徴は、標準的な突然のクエンチモデルでは欠落している。
• フォノン軟化: 計算された横波音響フォノン枝は、共有結合の弱体化と一致して、励起キャリア密度の増加に伴う単調な軟化を示す。軟化の大きさは波長依存性があり、ブリルアンゾーン境界付近で最大となる。
• 定量的合致: 励起キャリア密度を $n_{exc} = 0.0075$（単位胞あたり）、電子温度を 1000 K と仮定することで、モデルは実験のタイムトレース（欠陥や非調和性を完全に捉えきれないことを考慮するためのフェノメノロジー的な減衰項によってモデル化されたコヒーレント振動の減衰を含む）と良好な一致を示した。

意義および主張
本論文は、コヒーレント信号とインコヒーレント信号の扱いを統一することにより、時間分解拡散散乱を解釈するための体系的な第一原理的ルートを提供することを主張している。

• 統一: コヒーレントおよびインコヒーレントな寄与が、同一の密度行列形式から自然に生じることを示し、突然のクエンチの図式を限界ケースとして回収しつつ、電子－フォノン散乱効果へと拡張している。
• メカニズムの特定: 非自明なコヒーレントダイナミクスが、光誘起されたフォノン周波数シフトと、非平衡キャリア分布からの電子－フォノン散ブルの両方によって誘起されることを明示的に特定している。
• 限界と今後の展望: 著者らは、現在の実装が（スイッチング関数や共通の格子温度を用いることで）キャリアおよびフォノン集団の時間発展に対してフェノメノロジー的な取り扱いに依存していることを控えめに述べている。彼らは、フェノメノロジー的な近似を取り除き、背景の成長や非調和減衰の記述を改善するために、本スキームをEPIqに既に実装されている第一原理的なキャリアおよびフォノンダイナミクス（衝突積分）と統合することが次のステップであると特定している。また、「モード混合なし（no-mode-mixing）」近似はGeにおいては妥当であるが、対称性が低下する構造転移（例：VO$_2$）を伴う系では失敗することが予想されるため、より一般的なモード混合方程式が必要になると指摘している。