First-principles Floquet analysis from real-time propagation

Technische Samenvatting: First-Principles Floquet-analyse vanuit Real-Time Propagatie

Probleemstelling
De praktische toepassing van de Floquet-theorie op realistische materialen binnen first-principles simulaties staat voor aanzienlijke computationele en interpretatieve hindernissen. Traditionele benaderingen vereisen de expliciete constructie van een vergrote Floquet-Hamiltoniaan, die elektronische en fotonische vrijheidsgraden combineert in een massieve Hilbertruimte. Deze constructie brengt aanzienlijke computationele overhead met zich mee en vereist een willekeurige afkaping van de Fourierreeks (het kiezen van een maximaal aantal Floquet-componenten, $n_{max}$). Bovendien leveren standaard first-principles methoden zoals Time-Dependent Density Functional Theory (TDDFT) van nature real-time gepropageerde golffuncties op in plaats van Floquet-eigenstanden. Daarnaast omvatten experimentele scenario's vaak eindige pulsen en transiënte regimes die geen strikte tijdsperiodiciteit vertonen, condities waaronder het formele Floquet-kader strikt genomen niet van toepassing is, maar waaronder "Floquet-achtige" elektronische structuren frequent worden waargenomen in pump-probe experimenten.

Methodologie
De auteurs stellen een real-time Floquet-analyseframework voor dat quasi-energieën en Floquet-toestanden direct extraheert uit gepropageerde golffuncties, waardoor de constructie van de vergrote Floquet-Hamiltoniaan wordt omzeild. De methodologie bestaat uit drie primaire componenten:

1. Reconstructie van de Eén-Periode Evolutie-operator:
   In plaats van een statische Floquet-Hamiltoniaan te diagonaliseren, reconstrueert de methode de één-periode tijdsevolutie-operator, $U(t+T, t)$, met behulp van overlap tussen geëvolueerde toestanden. Gegeven een set orthonormale toestanden $\{|\psi_i(t)\rangle\}$ (bijv. Kohn-Sham orbitalen) op tijdstip $t$, worden de matrixelementen van de evolutie-operator berekend als:
   $$U_{ij}(t+T, t) = \langle \psi_i(t) | \psi_j(t+T) \rangle$$
   Het diagonaliseren van deze matrix levert de Floquet-quasi-energieën ($E_\alpha = -\arg(\lambda_\alpha)/T$) en de periodieke delen van de Floquet-toestanden ($|\phi_\alpha(t)\rangle$) binnen de subspace van de gepropageerde toestanden op. Deze aanpak voegt een verwaarloosbare computationele kosten toe aan bestaande tijd-propagatieschema's.

2. Floquet-band Unfolding via Harmonische Decompositie:
   Om de ambiguïteit van de gereduceerde zone-representatie op te lossen (waarbij quasi-energieën gedefinieerd zijn modulo de drijf-frequentie $\Omega$), introduceren de auteurs een unfolding-procedure. Ze definiëren een gewicht $w_\alpha$ gebaseerd op de norm van de nulde-harmonische component ($n=0$) van de Floquet-toestand:
   $$w_\alpha = |u^{(0)}_\alpha|$$
   waarbij $|u^{(0)}_\alpha\rangle$ verkregen wordt door de periodieke Floquet-toestand $|\phi_\alpha(t)\rangle$ te projecteren op de nulde harmonische. Dit gewicht kwantificeert hoe sterk een Floquet-toestand het karakter van de onderliggende evenwichtband behoudt. Door het spectrum te ontplooien naar een uitgebreid zone-schema en de banden te wegen met $w_\alpha$, kunnen fysiek relevante takken (die een grote $w_\alpha$ hebben) worden onderscheiden van hogere-orde foton-replica's (die een kleine $w_\alpha$ hebben).

3. Symmetrie-analyse van Sidebands:
   De gereconstrueerde real-space Floquet-golffuncties maken directe analyse mogelijk van de symmetrie-eigenschappen van individuele licht-geïnduceerde sidebands. Door de ruimtelijke distributie van de golffuncties te onderzoeken, kan men de symmetrie-karakter (bijv. glide-mirror symmetrie) van specifieke Floquet sidebands ($n=0, \pm 1, \dots$) bepalen en optische selectieregels verifiëren.

Belangrijkste Resultaten
Het framework werd gevalideerd via first-principles TDDFT-simulaties op vier verschillende materiaalsystemen:

• Grafeen (2D): Onder circulair gepolariseerd licht (CPL) identificeerde de methode succesvol de topologische gap-opening bij het Dirac-punt. De unfolding-procedure bleek essentieel; zonder deze verscheen het spectrum in de eerste Floquet Brillouin zone (FBZ) als een dichte, onbegrijpelijke verzameling gevouwen banden. De gewogen unfolding onthulde duidelijk de hybridisatie-gaps en de topologische gap.
• Monolaag Zwart Fosfor (2D): De methode demonstreerde het vermogen om de symmetrie-karakter van licht-geïnduceerde sidebands te resolveren. Voor armstoel-gepolariseerd drijven bevestigde de analyse dat de $n=0$ valentie- en conductiebanden hun evenwicht glide-mirror symmetrieën behielden (respectievelijk oneven en even), terwijl de $n=\pm 1$ sidebands een omgekeerde symmetrie vertoonden. Dit maakte een directe verklaring mogelijk voor geobserveerde band-repulsie via symmetrie-toegestane hybridisatie.
• Monolaag Hexagonaal Boriumnitride (2D): De benadering werd toegepast op een eindige CPL-puls, wat het nut ervan buiten strikte tijdsperiodiciteit aantoont. Door het systeem te analyseren in verschillende tijdvensters binnen de puls (tweede, vijfde en achtste cyclus), legden de auteurs transiënte Floquet-kenmerken vast. De resultaten toonden aan dat valley-afhankelijke band-renormalisatie (sterker bij de $-K$ valley voor links-CPL) opkwam en vervaagde in overeenstemming met de puls-envelope, wat de methode valideert voor transiënte niet-evenwicht regimes.
• Bulk Tin Sulfide (3D): De methode werd uitgebreid naar een volledig periodiek 3D bulk halfgeleider. De simulaties onthulden polarisatie-afhankelijke band-renormalisatie, waarbij armstoel-gepolariseerd drijven de $\Gamma X$ valley beïnvloedde en zigzag-gepolariseerd drijven de $\Gamma Y$ valley beïnvloedde, consistent met eerdere symmetrie-gebaseerde analyses.

Betekenis en Claims
Het artikel claimt een praktisch en computationeel efficiënt framework te bieden voor het verbinden van real-time simulaties met Floquet-observabelen. De primaire bijdragen zijn:

1. Computationele Efficiëntie: Het vermijdt de dure constructie en diagonalisatie van de vergrote Floquet-Hamiltoniaan, door enkel gebruik te maken van de overlaps van toestanden die al berekend zijn tijdens de standaard tijd-propagatie.
2. Interpreteerbaarheid: Het unfolding-schema gebaseerd op harmonische gewichten lost de ambiguïteit van gevouwen foton-replica's op, waardoor de identificatie van fysieke bandranden en hybridisatie-gaps eenvoudig wordt.
3. Generaliseerbaarheid: De methode is toepasbaar op zowel 2D monolagen als 3D bulkmaterialen en strekt zich uit tot eindige-puls regimes waar strikte tijdsperiodiciteit afwezig is.
4. Analytische Diepgang: Het maakt directe extractie van symmetrie-eigenschappen en selectieregels voor individuele licht-geïnduceerde sidebands mogelijk vanuit real-space golffuncties.

De auteurs concluderen dat de Floquet-theorie niet alleen dient als een hulpmiddel voor het construeren van periodiek gedreven steady states, maar ook als een krachtig analyseframework voor transiënte quasiparticle dynamica en niet-evenwicht elektronische structuren in realistische materialen.