First-principles Floquet analysis from real-time propagation

Sintesi Tecnica: Analisi di Floquet dai primi principi tramite propagazione in tempo reale

Problematica
L'applicazione pratica della teoria di Floquet a materiali realistici all'interno di simulazioni basate sui primi principi affronta significativi ostacoli computazionali e interpretativi. Gli approcci tradizionali richiedono la costruzione esplicita di un Hamiltoniano di Floquet ampliato, che combina i gradi di libertà elettronici e fotonici in uno spazio di Hilbert massiccio. Questa costruzione comporta un overhead computazionale sostanziale e necessita di una troncatura arbitraria della serie di Fourier (scelta di un numero massimo di componenti di Floquet, $n_{max}$). Inoltre, i metodi standard basati sui primi principi, come la Teoria del Funzionale della Densità Dipendente dal Tempo (TDDFT), produiscono naturalmente funzioni d'onda propagate nel tempo anziché autostati di Floquet. Inoltre, gli scenari sperimentali spesso coinvolgono impulsi finiti e regimi transitori che mancano di una stretta periodicità temporale, condizioni in cui il formalismo formale di Floquet non è strettamente applicabile, sebbene siano frequentemente osservate strutture elettroniche "tipo Floquet" negli esperimenti di pump-probe.

Metodologia
Gli autori propongono un framework di analisi di Floquet in tempo reale che estrae quasi-energie e stati di Floquet direttamente dalle funzioni d'onda propagate, evitando la costruzione dell'ampliato Hamiltoniano di Floquet. La metodologia consiste in tre componenti primarie:

1. Ricostruzione dell'operatore di evoluzione di un periodo:
   Invece di diagonalizzare un Hamiltoniano di Floquet statico, il metodo ricostruisce l'operatore di evoluzione temporale di un periodo, $U(t+T, t)$, utilizzando le sovrapposizioni tra stati evoluti nel tempo. Dato un insieme di stati ortonormali $\{|\psi_i(t)\rangle\}$ (ad esempio, orbitali di Kohn-Sham) al tempo $t$, gli elementi di matrice dell'operatore di evoluzione sono calcolati come:
   $$U_{ij}(t+T, t) = \langle \psi_i(t) | \psi_j(t+T) \rangle$$
   La diagonalizzazione di questa matrice fornisce le quasi-energie di Floquet ($E_\alpha = -\arg(\lambda_\alpha)/T$) e le parti periodiche degli stati di Floquet ($|\phi_\alpha(t)\rangle$) all'interno del sottospazio degli stati propagati. Questo approccio aggiunge un costo computazionale trascurabile agli esistenti schemi di propagazione temporale.

2. Unfolding delle bande di Floquet tramite decomposizione armonica:
   Per risolvere l'ambiguità della rappresentazione della zona ridotta (dove le quasi-energie sono definite modulo la frequenza di guida $\Omega$), gli autori introducono una procedura di unfolding (sviluppo). Definiscono un peso $w_\alpha$ basato sulla norma della componente dell'armonica zero ($n=0$) dello stato di Floquet:
   $$w_\alpha = |u^{(0)}_\alpha|$$
   dove $|u^{(0)}_\alpha\rangle$ è ottenuto proiettando lo stato periodico di Floquet $|\phi_\alpha(t)\rangle$ sull'armonica zero. Questo peso quantifica quanto fortemente uno stato di Floquet mantenga il carattere della banda sottostante all'equilibrio. L'unfolding dello spettro in uno schema di zona estesa, pesando le bande con $w_\alpha$, permette di distinguere i rami fisicamente rilevanti (quelli con grande $w_\alpha$) dalle repliche fotoniche di ordine superiore (che hanno un piccolo $w_\alpha$).

3. Analisi della simmetria delle sideband:
   Le funzioni d'onda di Floquet ricostruite nello spazio reale permettono l'analisi diretta delle proprietà di simmetria dei singoli sideband indotti dalla luce. Esaminando la distribuzione spaziale delle funzioni d'onda, è possibile determinare il carattere di simmetria (ad esempio, la simmetria di specchio di scorrimento o glide-mirror) di specifici sideband di Floquet ($n=0, \pm 1, \dots$) e verificare le regole di selezione ottica.

Risultati Chiave
Il framework è stato validato attraverso simulazioni TDDFT sui primi principi su quattro distinti sistemi di materiali:

• Grafene (2D): Sotto luce polarizzata circolarmente (CPL), il metodo ha identificato con successo l'apertura del gap topologico al punto di Dirac. La procedura di unfolding è stata dimostrata essenziale; senza di essa, lo spettro nella prima zona di Brillouin di Floquet (FBZ) appariva come un insieme denso e inintelligibile di bande ripiegate. L'unfolding pesato ha rivelato chiaramente i gap di ibridazione e il gap topologico.
• Fosforo Nero Monostrato (2D): Il metodo ha dimostrato la capacità di risolvere il carattere di simmetria dei sideband indotti dalla luce. Per la guida con polarizzazione armchair, l'analisi ha confermato che le bande di valenza e conduzione $n=0$ mantengono le rispettive simmetrie di specchio di scorrimento all'equilibrio (dispari e pari, rispettivamente), mentre i sideband $n=\pm 1$ mostrano simmetrie invertite. Ciò ha permesso una spiegazione diretta della repulsione tra bande osservata tramite l'ibridazione permessa dalla simmetria.
• Nitruro di Boro Esagonale Monostrato (2D): L'approccio è stato applicato a un impulso CPL finito, dimostrando la sua utilità oltre la stretta periodicità temporale. Analizzando il sistema in diversi intervalli temporali all'interno dell'impulso (secondo, quinto e ottavo ciclo), gli autori hanno catturato caratteristiche di Floquet transitorie. I risultati hanno mostrato che la rinormalizzazione della banda dipendente dalla valle (più forte nella valle $-K$ per la sinistra-CPL) emergeva e svaniva in accordo con l'inviluppo dell'impulso, validando il metodo per regimi non all'equilibrio transitori.
• Stagno Solfuro Bulk (3D): Il metodo è stato esteso a un semiconduttore bulk 3D completamente periodico. Le simulazioni hanno rivelato una rinormalizzazione della banda dipendente dalla polarizzazione, dove la guida armchair influenzava la valle $\Gamma X$ e la guida zigzag influenzava la valle $\Gamma Y$, in coerenza con precedenti analisi basate sulla simmetria.

Significatività e Rivendicazioni
L'articolo sostiene di fornire un framework pratico ed efficiente dal punto di vista computazionale per connettere le simulazioni in tempo reale con gli osservabili di Floquet. I suoi contributi primari sono:

1. Efficienza Computazionale: Evita la costosa costruzione e diagonalizzazione dell'ampliato Hamiltoniano di Floquet, utilizzando solo le sovrapposizioni degli stati già calcolati durante la normale propagazione temporale.
2. Interpretabilità: Lo schema di unfolding basato sui pesi armonici risolve l'ambiguità delle repliche fotoniche ripiegate, rendendo la identificazione dei bordi di banda fisici e dei gap di ibridazione immediata.
3. Generalità: Il metodo è applicabile sia a monostrati 2D che a materiali bulk 3D e si estende ai regimi di impulsi finiti dove la stretta periodicità temporale è assente.
4. Profondità Analitica: Permette l'estrazione diretta delle proprietà di simmetria e delle regole di selezione per i singoli sideband indotti dalla luce dalle funzioni d'onda nello spazio reale.

Gli autori concludono che la teoria di Floquet funge non solo da strumento per costruire stati stazionari guidati periodicamente, ma anche da potente framework di analisi per la dinamica dei quasi-particelle transitorie e le strutture elettroniche fuori dall'equilibrio in materiali realistici.