[論文レビュー] Evolution of non-thermal phonon and electron populations in photo-excited germanium on picosecond timescales

F. Murphy‐Armando, Éamonn Murray+4 more|arXiv (Cornell University)|Nov 27, 2019

Advanced Electron Microscopy Techniques and Applications被引用数 4

ひとこと要約

本研究は、LCLSにおける時間分解型X線散乱散乱実験と第一原理計算を組み合わせ、光励起状態のゲルマニウムにおける非平衡フォノンおよび電子のダイナミクスをピコ秒スケールでモデル化する。非平衡フォノンの主な発生源は、ブリユアンゾーンのL点近くにおけるΔおよびL伝導帯谷間のインターバレー電子-フォノン散乱であり、調整パラメータなしに実験データと一致するシミュレーションが得られた。

ABSTRACT

We investigate from first-principles theory and experiment the generation of phonons on picosecond timescales and the relaxation of carriers in multiple conduction band valleys of photo-excited Ge by inter-valley electron-phonon scattering. We provide a full description of the phonon and electron relaxation dynamics without adjustable parameters. Simulations of the time-evolution of phonon populations, based on first-principles band structure and electron-phonon and phonon-phonon matrix elements, are compared with data from time-resolved x-ray diffuse scattering experiments, performed at the LCLS x-ray free-electron laser facility, which measures the diffuse scattering intensity following photo-excitation by a 50 fs near-infrared optical pulse. Comparing calculations and measurements show that the intensity of the non-thermal x-ray diffuse scattering signal, that is observed to grow substantially near the L-point of the Brillouin zone over 3-5 ps, is due to phonons generated by scattering of carriers between the $Δ$ and L valleys. Non-thermal phonon populations throughout the Brillouin zone are observed and simulated from first principles without adjustable parameters for times up to 10 ps. With inclusion of phonon decay through 3-phonon processes, the simulations also account for other non-thermal features observed in the x-ray diffuse scattering intensity, which are due to anharmonic phonon-phonon scattering of the phonons initially generated by electron-phonon scattering.

研究の動機と目的

• 時間分解型X線散乱を用いて観測された光励起ゲルマニウムにおける非平衡フォノン分布の微視的起源を解明すること。
• ピコ秒スケールで複数の伝導帯谷（Γ、Δ、L）を介した電子およびフォノンの緩和ダイナミクスを解明すること。
• エネルギー保存および運動量保存を正確に反映する、調整パラメータなしの第一原理的シミュレーションフレームワークを構築すること。
• ブリユアンゾーンのL点近くで観測された散乱強度の3–5 psにわたる徐々な増加を説明すること。

提案手法

• 密度汎関数理論から導出した第一原理的バンド構造および電子-フォノン結合行列要素を用いた。
• 特にΔおよびL谷の周辺のような、小さな局所的領域における散乱過程を解像するため、適応的kポイントグリッドを採用した。
• 電子およびフォノン分布の連立レート方程式を定式化し、谷間および谷内電子-フォノン散乱および3フォノン非調和的崩壊を組み込んだ。
• エネルギー保存に配慮した重みを用いた四面体ベースの統合スキームを導入し、運動量空間全域における散乱核およびフォノン分布の計算を実施した。
• LCLS実験データとの直接比較を可能にするために、シミュレートされたフォノン分布を用いて時間依存のX線散乱構造因子を計算した。
• 散乱過程におけるエネルギー保存を、エネルギー保存表面に近い四面体頂点に重みを付けて厳密に保持した。

実験結果

リサーチクエスチョン

• RQ1 光励起GeにおいてブリユアンゾーンのL点近くに観測された非平衡フォノン分布の微視的起源は何か？
• RQ2 なぜL点近くの散乱強度は光励起直後ではなく、3–5 psにわたって徐々に増加するのか？
• RQ3 ΔおよびL伝導帯谷間のインターバレー電子-フォノン散乱プロセスは、フォノン生成および緩和にどのように寄与するか？
• RQ4 第一原理的シミュレーションは、調整パラメータなしに非平衡フォノンおよび電子分布の時間発展をどれほど正確に再現できるか？
• RQ5 3フォノン非調和的プロセスは、観測された非平衡フォノン特徴（散乱強度の広がりや非対称性など）をどのように形成するか？

主な発見

• L点近くの非平衡フォノン分布は、主にΔおよびL伝導帯谷間のインターバレー電子-フォノン散乱によって生成される。
• L点近くの散乱強度の増加は、行列要素が弱く、散乱の開始が遅れるという緩やかなプロセスに起因し、電子分布が最初にΔおよびL谷に蓄積されるため、3–5 psにわたって進行する。
• 第一原理パラメータに基づくシミュレーションは、調整パラメータなしに実験的X線散乱強度と極めて良好に一致した。
• 3フォノン非調和的フォノン-フォノン散乱は、散乱強度の広がりや非対称性といった、追加の非平衡的特徴を形成する。
• 適応的kポイントグリッドおよび四面体ベースのエネルギー保存統合スキームにより、非平衡ダイナミクスに不可欠な局所的散乱過程の正確なモデリングが可能になった。
• 四面体頂点の重み付き平均化スキームにより、散乱過程におけるエネルギー保存が厳密に維持され、シミュレーションの物理的整合性が保証された。