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

Sintesi Tecnica: Fotogenerazione e Firme di Fononi Coerenti nella Spettroscopia Fotoemissione Risolta nel Tempo

Problematica
La dinamica reticolare coerente, o fononi, eccitati da impulsi di luce offre un meccanismo per modulare le proprietà elettroniche e dei materiali, con potenziali applicazioni nel controllo quantistico e nei sensori. Sebbene gli esperimenti di pump-probe ottici (ad esempio, la riflettività transitoria) possano rilevare i fononi coerenti, spesso mancano della risoluzione necessaria per determinare le intensità di accoppiamento elettrone-fonone (e-ph) risolte per stato. La spettroscopia fotoemissione risolta nel tempo e nell'angolo (TR-ARPES) fornisce l'accesso alla struttura elettronica nello spazio dei momenti, permettendo teoricamente l'estrazione dell'accoppiamento e-ph risolto per stato tramite l'analisi di Fourier dei segnali ARPES. Tuttavia, un quadro teorico completo basato sul primo principio che affronti simultaneamente la generazione dei fononi coerenti e la loro successiva modulazione degli spettri TR-ARPES — includendo le interazioni elettrone-elettrone, elettrone-buco ed e-ph — è rimasto una sfida significativa. I metodi esistenti, come i calcoli statici con spostamenti imposti o i modelli fenomenologici, non riescono a catturare l'interazione tra la dinamica elettronica e quella fononica o i meccanismi specifici di generazione (Impulsive Stimulated Raman Scattering [ISRS] vs. Displacive Excitation of Coherent Phonons [DECP]).

Metodologia
Gli autori presentano un approccio ab initio adiabatico GW tempo-dipendente (TD-aGW) esteso per includere gli accoppiamenti e-ph. Il metodo simula un setup pump-probe per il MoS₂ monostrato, un materiale con significativi effetti eccitonici.

• Hamiltoniana e Dinamica: Il sistema è modellato propagando la densità di matrice del singolo elettrone $\rho_{nm}^{\mathbf{k}}(t)$ nella base delle bande di Bloch e il valore di aspettazione dello spostamento fononico $\langle Q_\nu(t) \rangle$. L'Hamiltoniana include il termine di quasiparticella GW all'equilibrio, le variazioni dell'autoelemento di Hartree e GW (approssimazione statica COHSEX), l'accoppiamento con il campo di pump esterno e il termine di accoppiamento e-ph nell'approssimazione classica del fonone.
• Equazioni del Moto: La dinamica elettronica è governata dall'equazione di Schrödinger tempo-dipendente con l'autoelemento GW. Lo spostamento fononico segue l'equazione del moto per un oscillatore armonico guidato, dove la forza motrice è determinata dagli elementi di matrice di accoppiamento e-ph e dalla densità di matrice tempo-dipendente.
• Parametri di Simulazione: Il lampo di pump è modellato come una funzione seno al quadrato a semiperiodo troncata (durata di 10 fs) con frequenze variabili (da 0,6 eV a 1,7 eV) e intensità variabili. Lo spettro TR-ARPES è simulato calcolando la funzione spettrale elettronica $A(\mathbf{k}, \omega, t)$ a intervalli di sonda di 10 fs. La dissipazione è inizialmente trascurata per isolare la dinamica intrinseca, sebbene effetti di dephasing finito siano brevemente esplorati nel Materiale Supplementare.

Risultati Chiave

1. Firme TR-ARPES dei Fononi Coerenti: Le simulazioni del MoS₂ monostrato rivelano che i fononi coerenti inducono modulazioni periodiche nell'intensità e nell'energia TR-ARPES. Nello specifico, lo spettro di differenza $\Delta A(\mathbf{k}, \omega, t)$ mostra spostamenti temporali nelle energie delle quasiparticelle. L'analisi di Fourier conferma che queste modulazioni corrispondono alla frequenza del modo fononico $A'_1$.
2. Estrazione dell'Accoppiamento e-ph Risolto per Stato: Lo studio dimostra una relazione lineare tra l'ampiezza dello spostamento spettrale ($\Delta A_{amp}$) e l'ampiezza del fonone coerente ($Q_{amp}$). Fondamentalmente, la pendenza di questa dipendenza lineare è direttamente proporzionale all'elemento di matrice di accoppiamento e-ph risolto per stato $g_{nn'}^{\nu}(\mathbf{k}, \mathbf{q}=0)$. Ciò suggerisce che gli esperimenti TR-ARPES possono determinare quantitativamente le intensità di accoppiamento e-ph per stati specifici.
3. Identificazione dei Meccanismi di Generazione: Le simulazioni identificano firme sia del meccanismo ISRS che del DECP, a seconda della frequenza di pump:
   • DECP (Eccitazione Displaciva): A frequenze di pump vicine al gap ottico (ad esempio, 1,7 eV, risonante con l'eccitone A), il sistema esibisce uno spostamento fononico medio finito ($Q_{avg} \neq 0$) e una forza motrice media non nulla derivante dalla densità di portatori fotoeccitati. Questo si allinea con il meccanismo DECP, in cui i portatori eccitati esercitano una forza che sposta la posizione di equilibrio.
   • ISRS (Scattering Raman Stimolato Impulsivo): A frequenze di pump più basse (ad esempio, 0,6 eV, sotto il gap), lo spostamento e la forza medi sono nulli, e la fase dell'oscillazione subisce uno shift di $\pi/2$ rispetto al caso DECP. Questo comportamento è coerente con il meccanismo ISRS, in cui la forza motrice è impulsiva e proporzionale all'intensità del pump senza lasciare una densità di portatori residua.
   • Transizione: La transizione tra questi regimi non è abrupta; frequenze intermedie (ad esempio, 1,2 eV) mostrano caratteristiche miste.
4. Regole di Selezione e Scattering Raman: Gli autori chiariscono la relazione tra la generazione di fononi coerenti e lo scattering Raman. Entrambi i processi condividono una regola di selezione derivata dallo stesso prodotto di matrice dei dipoli di transizione ottica e degli elementi di accoppiamento e-ph. Di conseguenza, i modi fononici che sono attivi al Raman sono generalmente generati. Tuttavia, la dipendenza in frequenza della forza motrice (generazione di fononi coerenti) differisce dall'intensità dello scattering Raman, in particolare al di sotto del gap ottico e sopra l'energia dell'eccitone C, indicando strutture di risonanza distinte nonostante la condivisione della stessa regola di selezione.

Significatività e Rivendicazioni
L'articolo sostiene di fornire la prima simulazione ab initio tempo-dipendente che integra completamente le interazioni elettrone-buco (tramite GW-BSE) e gli accoppiamenti e-ph per modellare sia la generazione di fononi coerenti che il loro impatto su TR-ARPES.

• Analisi Quantitativa: Il lavoro stabilisce un metodo per estrarre le intensità di accoppiamento e-ph risolte per stato direttamente dai dati TR-ARPES analizzando le modulazioni dei fononi coerenti.
• Chiarezza Meccanicistica: Il lavoro risolve il dibattito sui meccanismi di generazione dimostrando che ISRS e DECP derivano dallo stesso termine di driving nella risposta della densità di matrice, ma si manifestano diversamente a seconda del regime di frequenza del pump (trasparente vs. opaco).
• Approfondimento sulle Regole di Selezione: Lo studio chiarisce che, sebbene la generazione di fononi coerenti e l'intensità dello scattering Raman condividano le regole di selezione, i loro comportamenti spettrali non sono identici, offrendo una comprensione più sfumata delle interazioni luce-materia.
• Utilità Futura: Gli autori affermano che questo approccio supporta l'analisi della dinamica dei fononi coerenti in TR-ARPES e consente l'"ingegneria quantitativa della struttura a bande attraverso i fononi coerenti", pur rimanendo modesti riguardo a specifiche proposte sperimentali, notando che le applicazioni nel mondo reale devono tenere conto di dephasing ed effetti di depopolamento che complicano l'identificazione dei meccanismi di generazione nel regime dei cento femtosecondi.