结合GW/BSE实时波函数传播与开放量子系统的二维电子光谱理论框架

《Journal of Chemical Theory and Computation》：Time-Dependent Open-Quantum Approach to Two-Dimensional Electronic Spectroscopy within a GW/BSE Active Space

本文综述了（review）一种将GW/Bet-Salpeter方程（GW/BSE）电子结构、含时薛定谔方程（TDSE）/随机薛定谔方程（SSE）实时波函数传播、以及相位循环（phase-cycling）技术相结合的全新第一性原理计算框架，用于模拟真实激光脉冲激发条件下的二维电子光谱（2DES）。该方法在GW/BSE计算精度基础上，通过显式处理任意形状的超快激光脉冲，并引入开放量子系统理论描述退相干效应，实现了对分子体系（如苯、叶绿素b、苯酚二聚体）非线性光学响应的精确模拟，揭示了激发态吸收、受激发射及相干动力学等关键特征，为深入理解光合作用等光物理过程提供了强大的理论工具。

　　

本文致力于发展一个高精度、计算可行的理论计算框架，用以模拟真实实验条件下的二维电子光谱（2D Electronic Spectroscopy, 2DES）。2DES是一种强大的非线性光学技术，通过分析样品与三个超快激光脉冲序列相互作用后产生的非线性信号，可以得到一个关联激发频率与探测频率的二维谱图。这个谱图蕴含着丰富的量子相干和能量转移路径信息，在光合作用等光物理过程的研究中显示出巨大价值。然而，谱图中复杂的信号需要先进的理论工具进行精确解读。传统理论方法，如基于密度矩阵的微扰论，或基于线性组合的相位匹配（phase-matching）方案，通常采用脉冲极限（impulsive limit）近似，忽略了实际激光脉冲的形状和重叠效应。为了更贴近实验，本文提出了一种整合了多种前沿理论的全新方法。

本工作的核心理论框架（计算流程图可参考）主要包含四个紧密衔接的部分：GW/BSE电子活性空间构建、实时波函数传播、相位循环技术提取信号、以及二维傅里叶变换获得光谱。

首先，体系的电子结构和激发态信息由GW/Bet-Salpeter方程（GW/BSE） 方法提供。该方法建立在密度泛函理论（DFT）基础上，通过包含自能项G*W来改进电子结构，并求解BSE方程来精确获得包含激子效应的中性激发能E_λ和跃迁偶极矩?λ′|μ?^?|λ?。这些结果构成了后续实时动力学模拟的“电子活性空间”，即体系的基态|0?和一组激发态{|λ?}的完备集合，确保了模拟的精度。

其次，实时波函数传播是模拟光与物质相互作用动力学的核心。体系的含时波函数|Ψ(t)?在GW/BSE构建的活性空间{|λ?}中展开。当体系被视为封闭系统时，其演化遵循含时薛定谔方程（TDSE），其中的哈密顿量包含了外场相互作用项。入射光被建模为三个具有高斯包络形状的激光脉冲的叠加。通过求解TDSE，可以得到波函数随时间演化的系数C_λ(t)。为了描述环境（如溶剂、晶格振动）对体系的影响，本文还引入了开放量子系统的描述，即采用随机薛定谔方程（SSE） 来传播波函数。SSE中包含一个非厄米的耗散项和一个基于维纳过程的涨落项。本文重点考虑了纯电子退相位效应，通过设定退相位时间T₂（或退相位速率γ）来模拟。开放系统的密度矩阵需要通过平均大量独立的SSE随机轨迹来获得。

在得到实时演化的波函数后，体系的总极化强度可以通过计算跃迁偶极矩的期望值得到。然而，这个总极化包含了所有阶的响应。为了从中提取出与2DES实验相对应的三阶非线性信号，本文采用了相位循环（phase-cycling） 技术。具体而言，对三个入射脉冲设置12组不同的相位组合(φ₁, φ₂, φ₃)，分别进行波函数传播并计算总极化。然后，通过特定的线性组合（系数见文中Table 1），可以从这12个相位依赖的极化中分离出沿特定相位匹配方向（如-k₁+k₂+k₃和k₁-k₂+k₃）的信号，分别对应重聚（rephasing） 和非重聚（nonrephasing） 信号。

最后，将得到的时域信号S_R(τ, T, t′)和S_NR(τ, T, t′)对第一个延迟时间τ和第三个延迟时间t′进行二维傅里叶变换，即可得到频率域的重聚和非重聚2DES谱S_R(ω₁, T, ω₃)和S_NR(ω₁, T, ω₃)。通常，实验上报告的是二者的实部和，即吸收谱，其具有更好的分辨率。

本文将该方法应用于三个分子体系以验证其有效性。对苯分子的计算清晰地展示了2DES谱中基态漂白（ground-state bleach）、受激发射（stimulated emission）和激发态吸收（excited-state absorption） 的特征信号峰，其结果与已有理论和实验数据吻合。对叶绿素b的卟啉核心的计算，旨在模拟与光合作用相关的关键生色团的光物理过程，计算得到的谱图也成功揭示了其激发态动力学特征。此外，本文还研究了苯-苯酚二聚体模型，并通过结合GW/BSE动力学与SSE，模拟了环境诱导的电子退相干效应，分析了相干性随种群时间T的演化，为研究复杂生物体系中量子相干性的寿命和机制提供了理论洞察。

综上所述，本文提出的方法将第一性原理计算精度（GW/BSE）、对真实激光脉冲的显式处理、以及对开放量子系统效应的描述无缝集成到一个统一的框架中。这种方法不仅能够高精度地模拟和解释实验2DES谱，还能在原子层面揭示光激发后的超快动力学细节，包括量子相干和退相干过程，为从理论上探索光合作用、光催化、有机光电材料等领域中的复杂光物理和光化学现象提供了强有力的新工具。