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

Resumen Técnico: Dinámica de Fotoportadores y Fonones Fuera del Equilibrio desde Primeros Principios

Planteamiento del Problema
Los avances recientes en experimentos de bomba-sonda ultrafast han permitido el monitoreo de la dinámica de fotoportadores y fonones en semiconductores con alta precisión en escalas de tiempo de femtosegundos a picosegundos. Estos experimentos revelan fenómenos complejos como la fusión no térmica, la dinámica de fonones coherentes y las transiciones de fase inducidas por la luz. Sin embargo, los marcos teóricos existentes carecen de un enfoque unificado, predictivo y computacionalmente eficiente para describir estas dinámicas acopladas fuera del equilibrio. Los métodos previos, incluyendo las Funciones de Green fuera del Equilibrio (NEGF), las ecuaciones de Bloch de semiconductor y los enfoques de Boltzmann semiclásicos, a menudo omiten uno o más de los tres mecanismos de dispersión críticos: la dispersión portador-portador, portador-fonón y fonón-fonón. Además, muchos enfoques no logran tratar el acoplamiento luz-materia de forma explícita, ignoran las renormalizaciones de cuasipartículas dependientes del tiempo o dependen de parámetros fenomenológicos en lugar de entradas ab initio. La ausencia de un marco que trate simultáneamente los tres canales de dispersión en igualdad de condiciones, mientras contabiliza el movimiento atómico coherente y el apantallamiento dependiente del tiempo, limita la capacidad de interpretar cuantitativamente los datos experimentales.

Metodología
Los autores desarrollan un marco de muchos cuerpos basado en primeros principios para describir la dinámica en tiempo real de fotoportadores y fonones en semiconductores. La metodología se construye sobre los siguientes pilares:

1. Hamiltoniano y Base: El enfoque parte del Hamiltoniano cuántico completo de electrones y núcleos bajo la aproximación de Born-Oppenheimer. La dinámica electrónica se trata en una base de funciones de Wannier localizadas máximamente (MLWF), derivada de cálculos de la Teoría del Funcional de la Densidad (DFT) y la Teoría de la Perturbación de la Funcional de la Densidad (DFPT). Esta elección asegura integrales de dispersión consistentes con la gauge y permite un muestreo denso de momento mediante esquemas de interpolación directa.
2. Ecuaciones de Movimiento (EOMs): La evolución temporal está gobernada por ecuaciones cinéticas cuánticas para las ocupaciones electrónicas ($f_{kn}$), polarizaciones microscópicas ($p_{knm}$) y ocupaciones de fonones ($n_{q\nu}$). Estas EOMs se derivan del formalismo de la matriz de densidad, incorporando:
   • Acoplamiento Luz-Materia: Inclusión explícita de un campo eléctrico dependiente del tiempo mediante la sustitución de acoplamiento mínimo en la aproximación de dipolo.
   • Renormalización de Cuasipartículas: Correcciones dependientes del tiempo a los autovalores electrónicos y frecuencias de Rabi, y frecuencias de fonones, derivadas de la parte real de la autoenergía. Esto incluye términos de Hartree, COHSEX (agujero de Coulomb más intercambio apantallado), Fan-Migdal y Debye-Waller.
   • Integrales de Colisión: Términos de dispersión derivados de la parte imaginaria de la autoenergía utilizando el Ansatz de Kadanoff-Baym Generalizado (GKBA) y la aproximación de Markov. Estos incluyen:
     • Portador-Portador: Interacción de Coulomb apantallada (aproximación de Born de segundo orden).
     • Portador-Fonón: Aproximación Fan-Migdal con vértices electrón-fonón estáticamente apantallados.
     • Fonón-Fonón: Dispersión anharmónica derivada de constantes de fuerza de tercer orden (regla de oro de Fermi).
3. Movimiento Atómico Coherente: El marco incluye un mecanismo de excitación displaciva donde los portadores fotoexcitados generan fuerzas fuera del equilibrio sobre los iones. Estas fuerzas se calculan utilizando un enfoque de DFPT de respuesta lineal que tiene en cuenta el apantallamiento del vértice electrón-fonón por la distribución electrónica en evolución, yendo más allá de la dinámica de Ehrenfest estándar donde los orbitales son fijos.
4. Implementación Computacional: El método se implementa en la suite de código abierto epiq, con interfaz con Quantum ESPRESSO y Yambo. Utiliza pasos de tiempo adaptativos (Runge-Kutta Dormand-Prince) y paralelización (MPI) para manejar el costo computacional de evaluar las integrales de colisión sobre zonas de Brillouin densas.

Contribuciones Clave
El artículo presenta un marco teórico y computacional unificado que logra lo siguiente:

• Tratamiento Unificado de la Dispersión: Es el primer marco ab initio que trata simultáneamente la dispersión portador-portador, portador-fonón y fonón-fonón en igualdad de condiciones en la dinámica en tiempo real.
• Apantallamiento Dinámico y Renormalización: Captura explícitamente las renormalizaciones de energía de cuasipartículas dependientes del tiempo (reducción de la brecha de banda o bandgap) y los desplazamientos de frecuencia de fonones impulsados por la evolución de la distribución de fotoportadores y el apantallamiento.
• Interpolación Consistente con la Gauge: Al utilizar una representación de Wannier con un esquema de interpolación directa, el método permite un muestreo de momento ultra-denso necesario para integrales de dispersión precisas sin el costo prohibitivo de los cálculos ab initio directos en cada paso de tiempo.
• Dinámica de Fonones Coherentes: Proporciona un cálculo de primeros principios del movimiento atómico coherente y su amortiguamiento, derivando fuerzas del acoplamiento electrón-fonón apantallado en lugar de orbitales fijos.
• Escalabilidad: La implementación es escalable y computacionalmente eficiente, permitiendo comparaciones sin parámetros ajustables con condiciones experimentales de bomba-sonda a fluencias y temperaturas realistas.

Resultados
El marco se aplica a la monocapa de MoS$_2$ y nitruro de boro hexagonal (h-BN):

• MoS$_2$ (Dinámica Ultrafast):

  • Relajación de Portadores: Se demuestra que la inclusión de la dispersión portador-portador es esencial para una equilibración de portadores ultrafast realista (decenas de femtosegundos). Omitirla conduce a tiempos de relajación un orden de magnitud más largos.
  • Fonones Calientes y Termalización: La dispersión portador-fonón puebla inicialmente modos de fonones "calientes" específicos (por ejemplo, $A_{1g}$). Sin embargo, omitir la dispersión fonón-fonón resulta en una termalización de la red incorrecta a largo plazo y sobreestima el tiempo de amortiguamiento del fonón coherente $A_{1g}$ por un factor de dos.
  • Renormalización: El método captura una reducción significativa de la brecha de banda inducida por la fotoexcitación (de ~2.6 eV a ~1.8 eV) debido a los cambios en el apantallamiento, consistente con los datos de ARPES resueltos en el tiempo.
  • Movimiento Coherente: Los desplazamientos de fuerza calculados y los tiempos de amortiguamiento del modo $A_{1g}$ concuerdan bien con los datos experimentales de absorción transitoria, validando el tratamiento de las fuerzas apantalladas.

• h-BN (Propiedades de Cuasi-Equilibrio):

  • Apantallamiento y Estructura de Bandas: En el estado de cuasi-equilibrio tras la termalización de portadores, el marco cuantifica el aumento del apantallamiento dependiente de la fluencia. Esto conduce a una reducción significativa de la brecha de banda de partícula única.
  • Fusión de Excitones: El método demuestra la fusión de las características excitónicas (picos 1s y 2s) en el espectro de absorción óptica a medida que aumenta la densidad de fotoportadores, con el excitón 1s persistiendo hasta ~$3.8 \times 10^{13}$ cm$^{-2}$.
  • Renormalización de Fonones: Los desplazamientos de frecuencia de fonones calculados debido a los portadores excitados concuerdan estrechamente con los puntos de referencia de DFPT restringido (cDFPT), validando el enfoque de autoenergía para predecir las respuestas de la red en materiales fotoexcitados.

Significancia y Reivindicaciones
Los autores afirman que este trabajo establece una plataforma computacional versátil y cuantitativamente predictiva para interpretar experimentos ultrafast en semiconductores. Al tratar los tres mecanismos de dispersión y las renormalizaciones dependientes del tiempo en igualdad de condiciones, el marco resuelve las discrepancias encontradas en modelos simplificados (por ejemplo, al omitir la dispersión fonón-fonón o las interacciones portador-portador). La capacidad de comparar directamente con observables experimentales (reflectividad, absorción, vidas medias de fonones) sin parámetros ajustables representa un paso significativo hacia adelante en la descripción teórica de las interacciones luz-materia. El artículo posiciona este enfoque como una herramienta para guiar futuros experimentos e investigar fenómenos fuera del equilibrio, incluyendo las transiciones de fase estructural inducidas por la luz, mediante el acoplamiento con DFT restringido para escalas de tiempo más largas. Los autores declaran explícamente que sus resultados resaltan los regímenes donde las descripciones simplificadas fallan cualitativamente, enfatizando la necesidad de su tratamiento unificado para un modelado preciso.