First-principles Floquet analysis from real-time propagation

技术摘要：基于实时演化的第一性原理 Floquet 分析

问题陈述
将 Floquet 理论应用于实际材料面临着显著的计算与解释障碍，这在第一性原理模拟中尤为突出。传统方法需要显式构建一个扩大的 Floquet 哈密顿量，该哈密顿量将电子和光子自由度结合进一个巨大的希尔伯特空间。这种构建过程会产生巨大的计算开销，并且必须对傅里叶级数进行人为截断（选择最大 Floquet 分量数 $n_{max}$）。此外，标准的 first-principles 方法（如含时密度泛函理论 TDDFT）自然输出的是实时传播的波函数，而非 Floquet 本征态。此外，实验场景通常涉及有限脉冲和瞬态机制，这些情况并不具备严格的时间周期性，在这些条件下，形式上的 Floquet 框架并不严格适用，但在泵浦-探测实验中经常能观察到“类 Floquet”的电子结构。

方法论
作者提出了一种实时 Floquet 分析框架，直接从传播的波函数中提取准能量和 Floquet 态，从而绕过了构建扩大的 Floquet 哈密顿量的步骤。该方法由三个主要部分组成：

1. 单周期演化算符的重构：
   不同于对静态 Floquet 哈密顿量进行对角化，该方法通过演化态之间的重叠来重构单周期时间演化算符 $U(t+T, t)$。给定一组在时刻 $t$ 的正交态 $\{|\psi_i(t)\rangle\}$（例如 Kohn-Sham 轨道），演化算符的矩阵元计算如下：
   $$U_{ij}(t+T, t) = \langle \psi_i(t) | \psi_j(t+T) \rangle$$
   对该矩阵进行对角化即可得到 Floquet 准能量（$E_\alpha = -\arg(\lambda_\alpha)/T$）以及在传播态子空间内的周期性部分 Floquet 态（$|\phi_\alpha(t)\rangle$）。这种方法对现有的时间传播方案几乎不增加计算成本。

2. 通过谐波分解进行的 Floquet 能带展开：
   为了解决约化区表示中的歧义（其中准能量是模掉驱动频率 $\Omega$ 定义的），作者引入了一种展开程序。他们定义了一个权重 $w_\alpha$，该权重基于 Floquet 态零阶谐波分量的范数：
   $$w_\alpha = |u^{(0)}_\alpha|$$
   其中 $|u^{(0)}_\alpha\rangle$ 是通过将周期性 Floquet 态 $|\phi_\alpha(t)\rangle$ 投影到零阶谐波上获得的。该权重量化了 Floquet 态保留底层平衡态能带特征的强度。通过将能谱展开到扩展区方案中，并利用 $w_\alpha$ 对能带进行加权，可以区分出具有物理意义的分支（具有较大 $w_\alpha$ 的分支）与高阶光子副本（具有较小 $w_\alpha$ 的分支）。

3. 侧带的对称性分析：
   重构的实空间 Floquet 波函数允许直接分析单个光诱导侧带的对称性性质。通过检查波函数的空间分布，可以确定特定 Floquet 侧带（$n=0, \pm 1, \dots$）的对称性特征（例如滑移镜面对称性），并验证光学选择定则。

关键结果
该框架通过对四种不同材料系统的 TDDFT 第一性原理模拟得到了验证：

• 石墨烯 (2D)： 在圆偏振光 (CPL) 下，该方法成功识别了狄拉克点处的拓扑能隙开启。研究表明，展开程序是必不可少的；如果没有它，第一 Floquet 布里渊区 (FBZ) 中的能谱会呈现为密集且难以理解的折叠能带集合。加权展开清晰地揭示了杂化能隙和拓扑能隙。
• 单层黑磷 (2D)： 该方法展示了分辨光诱导侧带对称性特征的能力。对于扶手椅方向极化的驱动，分析确认了 $n=0$ 价带和导带保留了其平衡态的滑移镜面对称性（分别为奇和偶），而 $n=\pm 1$ 侧带则表现出相反的对称性。这为观察到的能带排斥提供了直接的对称性允许杂化的解释。
• 单层六方氮化硼 (2D)： 该方法被应用于有限脉冲 CPL，证明了其在严格时间周期性之外的实用性。通过在脉冲内的不同时间窗口（第二、第五和第八个周期）内进行分析，作者捕捉到了瞬态 Floquet 特征。结果显示，谷依赖的能带重整化（对于左旋 CPL，在 $-K$ 谷处更强）随脉冲包络的出现而出现，随之消失，验证了该方法在瞬态非平衡机制下的有效性。
• 体相硫化锡 (3D)： 该方法被扩展到完全周期的 3 维体相半导体。模拟揭示了极化相关的能带重整化现象：扶手椅方向极化的驱动影响 $\Gamma X$ 谷，而锯齿方向极化的驱动影响 $\Gamma Y$ 谷，这与之前的基于对称性的分析一致。

意义与主张
本文声称提供了一个实用的、计算高效的框架，用于将实时模拟与 Floquet 可观测量联系起来。其主要贡献在于：

1. 计算效率： 它避免了构建和对角化扩大的 Floquet 哈密顿量所带来的高昂成本，仅利用了在标准时间传播过程中已经计算出的状态重叠。
2. 可解释性： 基于谐波权重的展开方案解决了折叠光子副本的歧义，使得识别物理能带边缘和杂化能隙变得简单直观。
3. 通用性： 该方法既适用于二维单层材料，也适用于三维体相材料，并可扩展到不存在严格时间周期性的有限脉冲机制。
4. 分析深度： 它能够从实空间波函数中直接提取单个光诱导侧带的对称性属性和选择定则。

作者得出结论，Floquet 理论不仅是构建周期性驱动稳态的工具，也是研究现实材料中瞬态准粒子动力学和非平衡电子结构的强大分析框架。