激子驱动的Floquet工程第一性原理研究

《npj Quantum Materials》：Exciton-driven Floquet engineering from first principles

　　

Floquet工程利用周期性驱动在非平衡态下裁剪电子结构和材料性质。激子为这种调控提供了一条独特途径：其电子-空穴动力学可作为内部驱动，使Floquet工程超越直接激光缀饰。对该区域的微观描述需要同时获取激子动力学及其动量分辨光电子发射特征，这在真实材料中仍具挑战性。本研究开发了一个高效的第一性原理框架（framework），将实时含时密度泛函理论（real-time TDDFT）与修饰的长程修正交换关联核（MLRC xc kernel）、Floquet能带分析以及时间分辨和角分辨光电子能谱（TrARPES）模拟相结合。以单层六方氮化硼（hBN）作为基准体系，研究人员捕获了一个带隙内的激子共振（subgap excitonic resonance），并证明由此产生的振荡激子场即使在泵浦与探测脉冲无时间重叠的情况下，也能产生清晰的Floquet边带（Floquet sidebands）和能带重整化（band renormalization）。该研究方法为预测量子材料中由集体激发（collective excitations）驱动的非平衡Floquet现象建立了一条实用且计算高效的途径。

研究背景方面，在超快时间尺度上调控量子材料的电子结构是非平衡材料科学的核心目标之一。Floquet工程通过相干周期性运动创造平衡态下不存在的有效电子结构，传统上采用时间周期性驱动场缀饰Bloch电子。时间分辨和角分辨光电子能谱（TrARPES）已成为直接分辨瞬态Floquet能带结构的强大工具，已在多种量子材料中得到应用。从驱动源和缀饰自由度的角度看，驱动源不一定是外部激光场，也可以来自固体中的集体玻色激发，如相干声子和激子。实验上已在单层WS₂中报道了激子驱动的Floquet工程，基于动力学投影算符方法的模型研究也表明共振泵浦可在泵浦脉冲后维持激子序参量的相干振荡。然而，描述激子驱动的Floquet光谱仍面临重大理论挑战，涉及两个相互交织的任务：激子本身的描述（通常通过GW-BSE或含修饰交换关联核的TDDFT实现）以及超快光电子谱的模拟（可通过非平衡格林函数方法或实时波函数方法结合光电子提取方案实现）。这些方法在计算成本和实际适用性方面存在差异，因此亟需一种兼顾第一性原理精度和计算效率的策略。本研究采用OCTOPUS代码中实现的修饰长程修正（MLRC）交换关联核，在实时TDDFT框架内描述超快激子动力学，以单层hBN为基准体系，计算光吸收谱并识别带隙内激子共振，进而模拟激子场诱导的TrARPES强度图和Floquet能带结构。

在技术方法方面，研究人员主要采用以下关键技术：实时含时密度泛函理论（real-time TDDFT）在实空间网格和速度规范下求解含时Kohn-Sham方程；修饰长程修正（MLRC）交换关联矢势方法描述长程电子-空穴相互作用以捕获束缚激子；时间分辨和角分辨光电子能谱（TrARPES）模拟采用半周期系统含时表面通量方法（t-SURFFP）从光电子通量直接获取动量分辨光电子信号；Floquet能带分析通过实时Floquet分析框架从时间传播的Kohn-Sham态重构单周期演化算符并对角化获取准能和Floquet态。计算使用开源的OCTOPUS代码完成，单层hBN采用12×12×1的k点网格，DFT带隙为4.19 eV。

研究结果方面，首先，激子共振的 characterization。通过比较MLRC和TD-ALDA方法计算的光吸收谱，TD-ALDA谱在约4.19 eV带隙附近呈现弱吸收起始，随后在5至7 eV之间出现宽泛的带间跃迁特征。相比之下，MLRC计算在约3.6 eV处产生一个强烈而窄的峰，远低于TD-ALDA吸收边，对应于约0.59 eV的激子结合能。这种从高能带间连续态到带隙内共振的显著振子强度重分布是强束缚激子的特征，表明电子-空穴相互作用定性改变了hBN的光学响应，TD-ALDA无法捕获。结果与先前GW-BSE计算定性一致。其次，激子驱动Floquet态的光电子指纹。研究人员进行了共振和非共振激发条件下的TrARPES模拟。共振情况下，泵浦光子能量为3.6 eV（与激子能量共振），持续时间为20 fs，泵浦结束后用40.8 eV的弱探测场（24 fs）探测。MLRC计算中泵浦诱导的交换关联矢势在泵浦脉冲结束后继续振荡，而TD-ALDA中诱导矢势迅速衰减。MLRC的TrARPES谱沿Γ?K?Γ路径显示出一个额外的价带复制品，与原始能带相隔约激子能量，由于泵浦和探测脉冲无时间重叠，该特征可归因于持续激子场产生的Floquet边带而非外部泵浦场。TD-ALDA谱仅在平衡能带结构附近显示可观的谱权重，带隙内无清晰特征。MLRC谱还在K点价带顶（VBM）及其Floquet复制品附近表现出显著的能带重整化，表明导带和价带复制品之间存在强杂化。非共振情况下（泵浦光子能量3.2 eV），VBM的Floquet复制品显著弱于共振情况，表明非共振泵浦产生的振荡激子极化弱得多。K点处TrARPES谱权重的定量比较显示，TD-ALDA谱由一个约?5.9 eV的窄峰主导，带隙区域基本无谱权重；非共振MLRC结果保持类似的主价带峰，仅在约?2.3 eV处出现弱卫星峰；共振激发则导致显著的谱权重重分布，主价带峰移至约?6.2 eV，同时在约?2.6 eV处出现强附加峰，能量间隔与特征激子驱动能量相当，支持将高能特征解释为激子场诱导的Floquet复制品。第三，激子缀饰的Floquet能带分析。Floquet能带结构显示共振泵浦-探测模拟中存在清晰的带隙内边带，而平衡DFT能带结构中不存在。在K点附近，VBM的零阶Floquet复制品位于约?6.2 eV，而带隙内约?2.6 eV处出现额外的Floquet分支，该分支紧密重现价带的动量色散并向上移动约3.6 eV，与激子能量一致，可确认为价带的一阶激子诱导Floquet复制品。该分支在K点附近相对较强的Floquet权重及其与平衡导带的邻近性表明激子位移的价带复制品与导带态之间存在可观耦合。由于Floquet分析在外部泵浦结束后的探测区间进行，带隙内复制品提供了持续激子极化作为内部周期性驱动场维持电子结构Floquet缀饰的直接证据。非共振MLRC计算中带隙内仍可识别微弱价带复制品，但Floquet权重明显降低，平衡能带及其复制品几乎不变；TD-ALDA计算中几乎仅保留零阶Floquet能带，高阶复制品被强烈抑制。

讨论部分总结如下：本研究从第一性原理出发研究了单层hBN中激子驱动的Floquet物理，在同一实时框架内处理激子动力学和Floquet响应，建立了内部激子场、Floquet边带和能带重整化之间的直接联系。通过比较共振和非共振MLRC计算与TD-ALDA参考计算，研究人员确定了所观察Floquet特征的微观起源。共振激发下MLRC中交换关联矢势在外部泵浦消失后继续振荡，提供内部周期性驱动，产生显著的带隙内Floquet边带、TrARPES谱权重重分布和通过与激子诱导复制品杂化实现的能带重整化。非共振泵浦仅产生弱残余激子极化，TD-ALDA参考计算中高阶Floquet复制品被强烈抑制。TrARPES特征与从同一传播中提取的Floquet能带之间的一致性提供了激子场驱动缀饰的一致微观图像。方法上，该框架通过结合实时TDDFT与MLRC交换关联场、t-SURFFP光电子模拟和OCTOPUS中的直接Floquet分析，从同一时间传播中一致地获取激子动力学、TrARPES谱和Floquet态，无需完整的非平衡格林函数处理。该MLRC方法存在重要局限性：由于长程电子-空穴相互作用通过仅由少量材料依赖参数表征的宏观交换关联矢势表示，该方法主要能够重现单个主导束缚激子共振而非完整的激子态流形；与长程修正核的一般情况类似，获得物理合理的激子结合能可能伴随激子振子强度的过估计，因此激子峰高度不应视为定量精确，Floquet边带的谱权重和能带重整化幅度应在定性层面解释。此外，对于激子响应强各向异性的材料，简单宏观MLRC参数化可能不足以捕获电子-空穴相互作用的完整方向依赖性。未来可结合杂化或屏蔽交换泛函与t-SURFFP和实时Floquet分析以提供更准确的描述，并纳入第一性原理声子介导退相干和有限激子寿命以更真实地描述激子极化的持续和衰减。研究结论翻译如下：该研究方法为获取激子驱动的Floquet现象及其动量分辨特征提供了一条高效的第一性原理途径，通过将实时TDDFT与MLRC交换关联场、t-SURFFP光电子发射模拟和OCTOPUS中的直接Floquet分析相结合，可从同一时间传播中一致地获得激子动力学、TrARPES谱和Floquet态，而无需完整的非平衡格林函数处理。这种计算效率、物理透明性和对实验相关可观测量的直接获取之间的平衡，使该方法适合将激子驱动的Floquet研究扩展到更复杂的材料。