First-principles electron-phonon scattering in real-time TDDFT

Resumen Técnico: Dispersión Electrón-Fonón de Primeros Principios en Tiempo Real mediante TDDFT

Planteamiento del Problema
La teoría del funcional de la densidad dependiente del tiempo en tiempo real (rt-TDDFT) es una herramienta establecida de primeros principios para simular la dinámica coherente de electrones en sólidos excitados por láser. Sin embargo, su formulación estándar depende de una evolución unitaria de los orbitales de Kohn-Sham, lo que describe un sistema cuántico cerrado. En consecuencia, la rt-TDDFT convencional no puede capturar procesos irreversibles como la relajación de portadores, la termalización y la decoherencia derivadas del acoplamiento de los electrones excitados con su entorno. En materiales reales, las interacciones electrón-fonón (e-ph) son el mecanismo dominante para estos procesos disipativos, pero integrar estas interacciones en un marco de primeros principios en tiempo real, manteniendo el acceso a observables espectroscópicos, ha seguido siendo un desafío central. Los enfoques existentes suelen depender de modelos fenomenológicos o ecuaciones de Boltzmann semiclásicas que solo rastrean las poblaciones de portadores, fallando al describir la degradación de las coherencias cuánticas, esenciales para la dinámica ultrarrápida.

Metodología
Los autores desarrollan un marco de rt-TDDFT disipativo que trata el subsistema electrónico como un sistema cuántico abierto débilmente acoplado a un baño de fonones térmicos. La metodología integra tres componentes clave:

1. Formulación de la Matriz de Densidad de Sistema Abierto: Partiendo de los elementos de matriz e-ph de primeros principios, los autores derivan un integral de colisión para la matriz de densidad de un cuerpo reducida dentro de la aproximación de Born–Markov. La ecuación de movimiento para la matriz de densidad reducida $\hat{\rho}^{(1)}$ combina la evolución coherente de Liouville–von Neumann impulsada por el Hamiltoniano de Kohn-Sham dependiente del tiempo ($\hat{H}_{KS}$) con un término disipativo $I_{e-ph}$:
   $$\frac{\partial \hat{\rho}^{(1)}(t)}{\partial t} = -i [\hat{H}_{KS}(t), \hat{\rho}^{(1)}(t)] + I_{e-ph}[\hat{\rho}^{(1)}(t)]$$
   El integral de colisión $I_{e-ph}$ se construye para describir simultáneamente la redistribución de población asistida por fonones (términos de ganancia-pérdida) y la degradación de las coherencias interbanda (amortiguamiento fuera de la diagonal), asegurando que la evolución sea completamente positiva y preservadora de la traza (CPTP).

2. Propagación Híbrida de Espacio Real/Espacio Activo: Para mantener la viabilidad computacional preservando las capacidades completas de la rt-TDDFT en espacio real, los autores emplean una representación híbrida.

   • Dinámica Coherente: La propagación unitaria de los orbitales de Kohn-Sham se realiza en el espacio real utilizando el código Octopus, capturando interacciones luz-materia no lineales y efectos de continuo.
   • Actualizaciones Disipativas: El integral de colisión se aplica en un "espacio activo" seleccionado definido por una base adiabática dependiente del tiempo (autovalores instantáneos de $\hat{H}_{KS}$). La matriz de densidad se proyecta sobre esta base, se actualiza mediante el integral de colisión y luego se mapea de nuevo a la representación de orbitales en el espacio real.
   • Construcción de la Base: Para evitar la pérdida de norma y ambigüedades de calibre, el espacio de Hilbert se particiona en un espacio activo ($P$) y un espacio residual ($Q$). Un procedimiento de fijación de calibre basado en el algoritmo de Procrustes ortogonal asegura la continuidad de los orbitales entre pasos de tiempo.

3. Tasas de Transición de Gran Escala (Coarse-Grained): Dado que la dispersión e-ph de primeros principios requiere mallas de momento ultra finas para resolver energías de fonones de la escala de meV (prohibitivo para la propagación en tiempo real), los autores implementan un esquema de gran escala. Las tasas de transición microscópicas calculadas en una malla densa (usando el código EPW) se promedian sobre parches de momento para definir tasas efectivas para la malla más gruesa utilizada en la simulación de rt-TDDFT.

Resultados Clave
El marco fue validado mediante dos aplicaciones de referencia:

• Relajación de Portadores Ultrarrápida en Silicio Volumétrico (Bulk):

  • Al simular un pulso láser de 2.5 eV, el método capturó la redistribución de los portadores fotoexcitados.
  • Relajación de Energía vs. Momento: Los resultados revelaron una separación de escalas temporales: la redistribución de momento (dispersión fuera de los valles) ocurrió en una escala de tiempo de 8–14 fs, mientras que la relajación de energía (enfriamiento) fue más lenta (40–106 fs). Esto respalda la formación de un conjunto de electrones calientes en cuasi-equilibrio antes del enfriamiento macroscópico.
  • Acuerdo Cuantitativo: Los tiempos de vida de energía extraídos (40 fs a 2.6 eV, 106 fs a 1.6 eV) mostraron un acuerdo cuantitativo con las mediciones experimentales, validando el enfoque disipativo.

• Dinámica de Valles en Monocapa de WS$_2$:

  • El estudio modeló la dispersión entre valles tras la excitación con luz circularmente polarizada (CP) y linealmente polarizada (LP).
  • Despolarización de Valle: Bajo luz CP, la polarización de valle inicial (valle K) decayó a medida que los portadores se dispersaban hacia los valles K', Q y Q'. El tiempo de dispersión entre valles extraído fue de 20 fs (17 fs para LP), consistente con valores experimentales recientes (16 fs).
  • Conexión con tr-ARPES: Al acoplar la dinámica disipativa con el marco de fotoemisión tSURFF, los autores generaron espectros de tr-ARPES totalmente ab initio. Estas simulaciones reprodujeron con éxito la excitación selectiva de valle y la subsiguiente redistribución, vinculando las tasas de dispersión microscópicas con señales de fotoemisión medibles.

Significancia y Reivindicaciones
El artículo afirma proporcionar una ruta práctica de primeros principios para simular la dinámica cuántica disipativa en materiales cristalinos realistas sin recurrir a parámetros fenomenológicos. Sus principales contribuciones son:

1. Descripción Unificada de Relajación y Decoherencia: A diferencia de los enfoques de Boltzmann semiclásicos que solo rastrean poblaciones, este marco retiene la matriz de densidad de partícula única completa, permitiendo la descripción simultánea de la relajación de población y la degradación de la coherencia impulsada por los mismos canales de dispersión microscópica.
2. Integración con Espectroscopía: Al incrustar la actualización disipativa dentro del marco de rt-TDDFT en espacio real del código Octopus, el método preserva el acceso a todo el conjunto de observables, incluyendo las respuestas ópticas no lineales y las firmas espectroscópicas resueltas en el tiempo (ej. tr-ARPES).
3. Escalabilidad: La formulación de espacio activo limita el costo computacional del integral de colisión al subconjunto relevante de estados, haciendo que el enfoque sea viable para materiales realistas manteniendo la conexión rigurosa con la teoría de e-ph de primeros principios.

Los autores posicionan este trabajo como una realización controlada de un marco general de rt-TDDFT de sistema abierto, capaz de abordar fenómenos ultrarrápidos donde la excitación coherente y la dispersión disipativa compiten, como en la generación de armónicos altos y la espectroscopía de bomba-sonda. Señalan que, aunque la implementación actual asume un baño de fonones térmico y aproximaciones de Born–Markov, la estrategia es extensible a otros canales ambientales (ej. fotones, magnones) y regímenes que involucran efectos no markovianos fuertes.