Nonequilibrium Photocarrier and Phonon Dynamics from First Principles: a Unified Treatment of Carrier-Carrier, Carrier-Phonon, and Phonon-Phonon Scattering

技术摘要：基于第一性原理的非平衡光载流子与声子动力学

问题陈述
超快泵浦-探测实验的最新进展使得在飞秒至皮秒时间尺度上高精度监测半导体中的光载流子和声子动力学成为可能。这些实验揭示了复杂的现象，如非热熔 Melt、相干声子动力学以及光诱导相变。然而，现有的理论框架缺乏一种统一、具有预测性且计算高效的方法来描述这些耦合的非平衡动力学。以往的方法（包括非平衡格林函数 (NEGF)、半导体布洛赫方程和半经典玻尔兹曼方法）往往忽略了三种关键散射机制中的一个或多个：载流子-载流子、载流子-声子以及声子-声子散射。此外，许多方法未能显式处理光-物质耦合，忽略了随时间变化的准粒子重整化，或者依赖于唯象参数而非从头算输入。由于缺乏一个能够将所有三个散射通道置于平等地位进行处理，同时兼顾相干原子运动和随时间变化的屏蔽效应的框架，限制了定量解释实验数据的能力。

方法论
作者开发了一个基于第一性原理的多体框架，用于描述半导体中光载流子和声子的实时动力学。该方法建立在以下支柱之上：

1. 哈密顿量与基组： 该方法始于波恩-奥本海默近似下的电子与原子核的全量子哈密顿量。电子动力学在最大局域化 Wannier 函数 (MLWF) 基组下进行处理，该基组源自密度泛函理论 (DFT) 和密度泛函微扰理论 (DFPT) 计算。这种选择确保了规范一致的散射积分，并通过直接插值方案实现了稠密的动量采样。
2. 运动方程 (EOMs)： 时间演化由电子占据数 ($f_{kn}$)、微观极化 ($p_{knm}$) 和声子占据数 ($n_{q\nu}$) 的量子动力学方程控制。这些 EOMs 源自密度矩阵形式，并包含：
   • 光-物质耦合： 通过偶极近似中的最小耦合替换，显式引入随时间变化的电场。
   • 准粒子重整化： 来自自能实部的电子本征值、拉比频率以及声子频率的随时间变化的修正。这包括 Hartree、COHSEX（库仑空穴加屏蔽交换）、Fan-Migdal 以及 Debye-Waller 项。
   • 碰撞积分： 利用广义 Kadanoff-Baym 假设 (GKBA) 和马尔可夫近似，从自能虚部推导出的散射项。其中包括：
     • 载流子-载流子： 屏蔽库仑相互作用（二阶 Born 近似）。
     • 载流子-声子： 采用静态屏蔽电子-声子顶角的 Fan-Migdal 近似。
     • 声子-声子： 源自三阶力常数的非谐性散射（Fermi 金律）。
3. 相干原子运动： 该框架包含了位移激发机制，即光激发载流子产生对离子的非平衡作用力。这些作用力使用线性响应 DFPT 方法计算，该方法考虑了由演化的电子分布引起的电子-声子顶点屏蔽，超越了轨道固定的标准 Ehrenfest 动力学。
4. 计算实现： 该方法在开源套件 epiq 中实现，并与 Quantum ESPRESSO 和 Yambo 对接。它利用自适应步长（Runge-Kutta Dormand-Prince）和并行化（MPI）技术，以处理在稠密布里渊区内评估碰撞积分的计算成本。

核心贡献
本文提出了一个统一的理论与计算框架，实现了以下目标：

• 统一的散射处理： 这是首个在实时动力学中同时对载流子-载流子、载流子-声子和声子-声子散射进行平等处理的从头算框架。
• 动态屏蔽与重整化： 它显式捕捉了由演化的光载流子分布和屏蔽效应驱动的随时间变化的准粒子能量重整化（带隙收缩）和声子频率偏移。
• 规范一致的插值： 通过利用具有直接插值方案的 Wannier 表示，该方法能够实现精确散射所需的超稠密动量采样，而不会产生直接从头算在每个时间步进行计算的巨大成本。
• 相干声子动力学： 它提供了相干原子运动和阻尼的第一性原理计算，其作用力源自屏蔽后的电子-声子耦合，而非固定的轨道。
• 可扩展性： 该实现具有可扩展性和计算效率，允许在现实的流强（fluence）和温度条件下，与实验泵浦-探测条件进行无参数对比。

结果
该框架应用于单层 MoS$_2$ 和六方氮化硼 (h-BN)：

• MoS$_2$（超快动力学）：

  • 载流子弛豫： 研究表明，包含载流子-载流子散射对于实现真实的超快载流子平衡（数十飞秒）至关重要。忽略它会导致弛豫时间比实际长一个数量级。
  • 热声子与热化： 载流子-声子散射最初会填充特定的“热”声子模式（如 $A_{1g}$）。然而，忽略声子-声子散射会导致错误的长时间晶格热化，并将相干 $A_{1g}$ 声子的阻尼时间高估两倍。
  • 重整化： 该方法捕捉到了显著的光诱导带隙减小（从 ~2.6 eV 降至 ~1.8 eV），这与时间分辨 ARPES 数据一致。
  • 相干运动： 计算出的位移作用力和 $A_{1g}$ 模式阻尼时间与实验瞬态吸收数据吻合良好，验证了屏蔽作用力处理的有效性。

• h-BN（准平衡性质）：

  • 屏蔽与能带结构： 在载流子热化后的准平衡态中，该框架量化了随流强增加的屏蔽增强效应。这导致了单粒子带隙的显著降低。
  • 激子熔解： 该方法展示了随着光载流子密度的增加，光学吸收光谱中激子特征（1s 和 2s 峰）的熔解过程，其中 1s 激子在 ~$3.8 \times 10^{13}$ cm$^{-2}$ 时仍然存在。
  • 声子重整化： 由于激发载流子导致的计算声子频率偏移与约束 DFPT (cDFPT) 基准高度一致，验证了该自能方法在预测光激发材料晶格响应方面的有效性。

意义与主张
作者声称，这项工作建立了一个通用的、具有定量预测能力的计算平台，用于解释半导体中的超快实验。通过将所有三种散射机制和随时间变化的重整化置于平等地位进行处理，该框架解决了简化模型（例如忽略声子-声子散射或载流子-载流子相互作用）中存在的差异。无需调节参数即可直接与实验观测值（反射率、吸收率、声子寿命）进行比较，这代表了光-物质相互作用理论描述迈出的重要一步。该论文将此方法定位为指导未来实验和研究非平衡现象（包括通过耦合约束 DFT 以研究更长时间尺度的光诱导结构相变）的工具。作者明确指出，其结果强调了简化描述在定性失效的机制，凸显了其统一处理对于准确建模的必要性。