Photogeneration and signatures of coherent phonons in time-resolved photoemission spectroscopy: First-principles time-dependent adiabatic GW approach

Resumen Técnico: Fotogeneración y Firmas de Fonones Coherentes en Espectroscopía de Fotoemisión Resuelta en el Tiempo

Planteamiento del Problema
La dinámica de red coherente, o fonones, excitada por pulsos de luz ofrece un mecanismo para modular las propiedades electrónicas y materiales, con aplicaciones potenciales en el control cuántico y sensores. Aunque los experimentos ópticos de bomba-sonda (por ejemplo, reflectividad transitoria) pueden detectar fonones coherentes, a menudo carecen de la resolución para determinar las fuerzas de acoplamiento electrón-fonón (e-ph) resueltas por estado. La espectroscopía de fotoemisión resuelta en el tiempo y en el ángulo (TR-ARPES) proporciona acceso a la estructura electrónica en el espacio de momentos, lo que teóricamente permite la extracción del acoplamiento e-ph resuelto por estado mediante el análisis de Fourier de las señales de ARPES. Sin embargo, un marco teórico integral de primeros principios que aborde simultáneamente la generación de fonones coherentes y su posterior modulación de los espectros de TR-ARPES —incluyendo las interacciones electrón-electrón, electrón-hueco y e-ph— ha representado un desafío significativo. Los métodos existentes, como los cálculos estáticos con desplazamientos impuestos o los modelos fenomenológicos, no logran capturar la interacción entre la dinámica electrónica y de fonones o los mecanismos de generación específicos (Dispersión Raman Estimulada Impulsiva [ISRS] frente a Excitación Displaciva de Fonones Coherentes [DECP]).

Metodología
Los autores presentan un enfoque de GW adiabático dependiente del tiempo (TD-aGW) ab initio extendido para incluir acoplamientos e-ph. El método simula un montaje de bomba-sonda para MoS₂ monocapa, un material con efectos excitónicos significativos.

• Hamiltono y Dinámica: El sistema se modela propagando la densidad de matriz de un solo electrón $\rho_{nm}^{\mathbf{k}}(t)$ en la base de bandas de Bloch y el valor de la esperanza del desplazamiento de fonón $\langle Q_\nu(t) \rangle$. El Hamiltoniano incluye el término de cuasipartícula GW en equilibrio, cambios en la autoenergía de Hartree y GW (aproximación COHSEX estática), el acoplamiento al campo de bomba externo y el término de acoplamiento e-ph en la aproximación de fonón clásico.
• Ecuaciones de Movimiento: La dinámica electrónica está gobernada por la ecuación de Schrödinger dependiente del tiempo con la autoenergía GW. El desplazamiento del fonón sigue la ecuación de movimiento para un oscilador armónico impulsado, donde la fuerza impulsora está determinada por los elementos de matriz de acoplamiento e-ph y la densidad de matriz dependiente del tiempo.
• Parámetros de Simulación: El pulso de bomba se modela como una función seno al cuadrado de medio periodo truncada (duración de 10 fs) con frecuencias variables (0.6 eV a 1.7 eV) e intensidades. El espectro de TR-ARPES se imita computando la función espectral electrónica $A(\mathbf{k}, \omega, t)$ en intervalos de sonda de 10 fs. Inicialmente se desprecia la disipación para aislar la dinámica intrínseca, aunque en el Material Suplementario se exploran brevemente los efectos de desfasaje finito.

Resultados Clave

1. Firmas de TR-ARPES de Fonones Coherentes: Las simulaciones de MoS₂ monocapa revelan que los fonones coherentes inducen modulaciones periódicas en la intensidad y la energía de TR-ARPES. Específicamente, el espectro de diferencia $\Delta A(\mathbf{k}, \omega, t)$ muestra desplazamientos dependientes del tiempo en las energías de las cuasipartículas. El análisis de Fourier confirma que estas modulaciones coinciden con la frecuencia del modo fonónico $A'_1$.
2. Extracción del Acoplamiento e-ph Resuelto por Estado: El estudio demuestra una relación lineal entre la amplitud del desplazamiento espectral ($\Delta A_{amp}$) y la amplitud del fonón coherente ($Q_{amp}$). Crucialmente, la pendiente de esta dependencia lineal es directamente proporcional al elemento de matriz de acoplamiento e-ph resuelto por estado $g_{nn'}^{\nu}(\mathbf{k}, \mathbf{q}=0)$. Esto sugiere que los experimentos de TR-ARPES pueden determinar cuantitativamente las fuerzas de acoplamiento e-ph para estados electrónicos específicos.
3. Identificación de Mecanismos de Generación: Las simulaciones identifican firmas de ambos mecanismos, ISRS y DECP, dependiendo de la frecuencia de la bomba:
   • DECP (Excitación Displaciva): A frecuencias de bomba cercanas a la brecha óptica (p. ej., 1.7 eV, resonante con el excitón A), el sistema exhibe un desplazamiento promedio de fonón finito ($Q_{avg} \neq 0$) y una fuerza impulsora promedio no nula derivada de la densidad de portadores fotoexcitados. Esto se alinea con el mecanismo DECP, donde los portadores excitados ejercen una fuerza que desplaza la posición de equilibrio.
   • ISRS (Dispersión Raman Estimulada Impulsiva): A frecuencias de bomba más bajas (p. ej., 0.6 eV, por debajo de la brecha), el desplazamiento promedio y la fuerza son cero, y la fase de la oscilación se desplaza $\pi/2$ respecto al caso DECP. Este comportamiento es consistente con el mecanismo ISRS, donde la fuerza impulsora es impulsiva y proporcional a la intensidad de la bomba sin dejar una densidad de portadores residual.
   • Transición: La transición entre estos regímenes no es abrupta; las frecuencias intermedias (p. ej., 1.2 eV) muestran características mixtas.
4. Reglas de Selección y Dispersión Raman: Los autores aclaran la relación entre la generación de fonones coherentes y la dispersión Raman. Ambos procesos comparten una regla de selección derivada del mismo producto de matrices de los dipolos de transición óptica y los elementos de acoplamiento e-ph. En consecuencia, los modos de fonones que son activos en Raman se generan generalmente. Sin embargo, la dependencia de la frecuencia de la fuerza impulsora (generación de fonones coherentes) difiere de la intensidad de la dispersión Raman, particularmente por debajo de la brecha óptica y por encima de la energía del excitón C, indicando estructuras de resonancia distintas a pesar de compartir la misma regla de selección.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera simulación dependiente del tiempo de primeros principios que integra plenamente las interacciones electrón-hueco (vía GW-BSE) y los acoplamientos e-ph para modelar tanto la generación de fonones coherentes como su impacto en TR-ARPES.

• Análisis Cuantitativo: El trabajo establece un método para extraer las fuerzas de acoplamiento e-ph resueltas por estado directamente de los datos de TR-ARPES mediante el análisis de las modulaciones de fonones coherentes.
• Claridad Mecanística: Resuelve el debate sobre los mecanismos de generación al demostrar que ISRS y DECP surgen del mismo término de respuesta de la densidad de matriz pero se manifiestan de manera diferente dependiendo del régimen de frecuencia de la bomba (transparente frente a opaco).
• Perspectiva de las Reglas de Selección: El estudio aclara que, si bien la generación de fonones coherentes y la intensidad de la dispersión Raman comparten reglas de selección, sus comportamientos espectrales no son idénticos, ofreciendo una comprensión más matizada de las interacciones luz-materia.
• Utilidad Futura: Los autores afirman que este enfoque respalda el análisis de la dinámica de fonones coherentes en TR-ARPES y permite la "ingeniería cuantitativa de la estructura de bandas mediante fonones coherentes", aunque mantienen la modestia respecto a propuestas experimentales específicas, señalando que las aplicaciones del mundo real deben tener en cuenta los efectos de desfasaje y depopulación que complican la identificación de los mecanismos de generación en el régimen de los cientos de femtosegundos.