First-principles Floquet analysis from real-time propagation

Resumen Técnico: Análisis de Floquet desde Primeros Principios a partir de la Propagación en Tiempo Real

Planteamiento del Problema
La aplicación práctica de la teoría de Floquet a materiales realistas dentro de simulaciones de primeros principios enfrenta obstáculos computacionales e interpretativos significativos. Los enfoques tradicionales requieren la construcción explícita de un Hamiltoniano de Floquet ampliado, que combina los grados de libertad electrónicos y fotónicos en un espacio de Hilbert masivo. Esta construcción incurre en una sobrecarga computacional sustancial y requiere la truncación arbitraria de la serie de Fourier (eligiendo un número máximo de componentes de Floquet, $n_{max}$). Además, los métodos estándar de primeros principios como la Teoría del Funcional de la Densidad Dependiente del Tiempo (TDDFT) producen naturalmente funciones de onda propagadas en el tiempo real en lugar de autoestados de Floquet. Asimismo, los escenarios experimentales suelen involucrar pulsos finitos y regímenes transitorios que carecen de una periodicidad temporal estricta, condiciones bajo las cuales el marco formal de Floquet no es estrictamente aplicable, pero donde se observan frecuentemente estructuras electrónicas "tipo Floquet" en experimentos de bomba-sonda (pump-probe).

Metodología
Los autores proponen un marco de análisis de Floquet en tiempo real que extrae cuasienergías y estados de Floquet directamente de las funciones de onda propagadas, evitando la construcción del Hamiltoniano de Floquet ampliado. La metodología consta de tres componentes principales:

1. Reconstrucción del Operador de Evolución de un Periodo:
   En lugar de diagonalizar un Hamiltoniano de Floquet estático, el método reconstruye el operador de evolución de un periodo, $U(t+T, t)$, utilizando los solapamientos entre estados evolucionados en el tiempo. Dado un conjunto de estados ortonormales $\{|\psi_i(t)\rangle\}$ (por ejemplo, orbitales de Kohn-Sham) en un tiempo $t$, los elementos de matriz del operador de evolución se calculan como:
   $$U_{ij}(t+T, t) = \langle \psi_i(t) | \psi_j(t+T) \rangle$$
   La diagonalización de esta matriz produce las cuasienergías de Floquet ($E_\alpha = -\arg(\lambda_\alpha)/T$) y las partes periódicas de los estados de Floquet ($|\phi_\alpha(t)\rangle$) dentro del subespacio de los estados propagados. Este enfoque añade un costo computacional insignificante a los esquemas de propagación temporal existentes.

2. Desplegado de Bandas de Floquet mediante Descomposición Armónica:
   Para resolver la ambigüedad de la representación de zona reducida (donde las cuasienergías se definen módulo la frecuencia de conducción $\Omega$), los autores introducen un procedimiento de desplegado (unfolding). Definen un peso $w_\alpha$ basado en la norma del componente de armónico cero ($n=0$) del estado de Floquet:
   $$w_\alpha = |u^{(0)}_\alpha|$$
   donde $|u^{(0)}_\alpha\rangle$ se obtiene proyectando el estado de Floquet periódico $|\phi_\alpha(t)\rangle$ sobre el armónico cero. Este peso cuantifica qué tan fuertemente un estado de Floquet retiene el carácter de la banda de equilibrio subyacente. Al desplegar el espectro en un esquema de zona extendida y ponderar las bandas con $w_\alpha$, se distinguen las ramas físicamente relevantes (aquellas con un $w_\alpha$ grande) de los réplicas de fotones de orden superior (que tienen un $w_\alpha$ pequeño).

3. Análisis de Simetría de las Bandas Laterales (Sidebands):
   Las funciones de onda de Floquet reconstruidas en el espacio real permiten el análisis directo de las propiedades de simetría de las bandas laterales inducidas por la luz individuales. Al examinar la distribución espacial de las funciones de onda, se puede determinar el carácter de simetría (por ejemplo, la simetría de espejo de deslizamiento o glide-mirror) de las bandas laterales de Floquet específicas ($n=0, \pm 1, \dots$) y verificar las reglas de selección óptica.

Resultados Clave
El marco fue validado mediante simulaciones TDDFT de primeros principios en cuatro sistemas de materiales distintos:

• Grafeno (2D): Bajo luz polarizada circularmente (CPL), el método identificó con éxito la apertura de la brecha topológica en el punto de Dirac. Se demostró que el procedimiento de desplegado es esencial; sin él, el espectro en la primera zona de Brillouin de Floquet (FBZ) aparecía como un conjunto denso e ininteligible de bandas plegadas. El desplegado ponderado reveló claramente las brechas de hibridación y la brecha topológica.
• Fósforo Negro Monocapa (2D): El método demostró la capacidad de resolver el carácter de simetría de las bandas laterales inducidas por la luz. Para la conducción polarizada en dirección armchair, el análisis confirmó que las bandas de valencia y conducción $n=0$ mantuvieron sus simetrías de espejo de deslizamiento de equilibrio (impares y pares, respectivamente), mientras que las bandas laterales $n=\pm 1$ exhibieron simetrías invertidas. Esto permitió una explicación directa de la repulsión de bandas observada mediante la hibridación permitida por simetría.
• Nitruro de Boro Hexagonal Monocapa (2D): El enfoque se aplicó a un pulso de CPL finito, demostrando su utilidad más allá de la periodicidad temporal estricta. Al analizar el sistema en diferentes ventanas de tiempo dentro del pulso (segundo, quinto y octavo ciclo), los autores capturaron características de Floquet transitorias. Los resultados mostraron que la renormalización de bandas dependiente del valle (más fuerte en el valle $-K$ para CPL a la izquierda) emergía y desaparecía de acuerdo con la envolvente del pulso, validando el método para regímenes no equilibrados transitorios.
• Sulfuro de Estaño Masivo (3D): El método se extendió a un semiconductor masivo (bulk) totalmente periódico en 3D. Las simulaciones revelaron una renormalización de bandas dependiente de la polarización, donde la conducción polarizada en armchair afectaba al valle $\Gamma X$ y la conducción en zigzag afectaba al valle $\Gamma Y$, consistente con análisis previos basados en simetría.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco práctico y computacionalmente eficiente para conectar las simulaciones en tiempo real con los observables de Floquet. Sus contribuciones primarias son:

1. Eficiencia Computacional: Evita la costosa construcción y diagonalización del Hamiltoniano de Floquet ampliado, utilizando únicamente los solapamientos de los estados ya computados durante la propagación temporal estándar.
2. Interpretabilidad: El esquema de desplegado basado en pesos armónicos resuelve la ambigüedad de las réplicas de fotones plegadas, haciendo que la identificación de los bordes de banda físicos y las brechas de hibridación sea sencilla.
3. Generalidad: El método es aplicable tanto a monocapas 2D como a materiales masivos 3D y se extiende a regímenes de pulsos finitos donde la periodicidad temporal estricta está ausente.
4. Profundidad Analítica: Permite la extracción directa de las propiedades de simetría y las reglas de selección para las bandas laterales de Floque individuales a partir de las funciones de onda en el espacio real.

Los autores concluyen que la teoría de Floquet sirve no solo como una herramienta para construir estados estacionarios conducidos periódicamente, sino también como un poderoso marco de análisis para la dinámica de cuasipartículas transitorias y las estructuras electrónicas fuera del equilibrio en materiales realistas.