《JACS Au》：Transient X-ray Absorption Signatures of Photodissociation Pathways and Conical Intersection Dynamics in Phenol
苯酚：理解超快非辐射动力学的典范模型
苯酚是理解生物芳香体系中超快非辐射动力学和光稳定性的典范模型。其基本光化学过程——紫外激发后的超快O-H键裂解——一直是大量实验和理论研究的焦点。作为核碱基的最简单结构类似物，苯酚展现出由羟基介导的丰富激发态势能面，该羟基耦合了所有主要轨道类别：芳香π/π系统、O-H σ/σ键、氧孤对电子(n)和里德伯轨道。这种全面的轨道相互作用使苯酚成为光化学的经典教科书模型。在紫外激发到明亮的ππ态后，强耦合且能量接近的πσ态成为O-H键裂解的主要通道。
低位的单重激发态，特别是最低的三个(S0–S2)，被发现对O-H键裂解机制至关重要。O-H距离(ROH)无疑是唯一的、能通过势能曲线(PECs)捕捉基本解离化学的主要反应坐标，尽管其他内坐标(如C-O键长和C-O-H角)也调制该过程。依赖于ROH的PECs中的光化学涉及两个关键的锥交叉(CoIns)，它们引导激发态布居：S2/S1CoIn (CI1)和S1/S0CoIn (CI2)。在初始通过CI1弛豫后，S1PEC上演化的波包在CI2处分叉，导致三条竞争路径：
- •
Path 1：在排斥的S1势能面上直接解离，生成H + PhO•(?)。
- •
Path 2：非绝热跃迁到S0，随后快速解离为H + PhO•(X?)。
- •
Path 3：一个重穿通道，波包在CI2附近暂时被困。
机理洞察的一个基石是动能释放(KER)谱的测量，其对于氢消除通道显示出一致的双峰分布，由分离良好的高动能和低动能成分标记。然而，这些成分的归属一直是不断演变的解释主题。解决所有竞争路径之间的精确分支，包括CI2处的分支(Paths 1–3)，以及从S1直接内转换后的统计基态衰变路径，需要一种能够以电子和原子特异性区分它们的探针。
瞬态X射线吸收光谱：通往原子分辨动力学的窗口
在超快技术中，泵浦-探测瞬态X射线吸收光谱(TXAS)在捕捉光诱导过程中的电子和核动力学方面尤其强大。然而，解释需要可靠的理论支持。当前的策略涉及沿分子动力学轨迹或沿静态势能面计算核心激发态。轨迹策略面临可扩展性挑战，而非绝热分子动力学(NAMD)方法虽然能直接获取时间依赖的密度矩阵，但需要大量的轨迹采样以进行可靠的系综平均，计算成本随传播时间显著增长。相比之下，静态方法更高效，沿反应坐标计算垂直核心激发。它已成功模拟了多种系统的TXAS。
本文采用静态“沿PESs的快照”方法模拟了光激发苯酚的氧K边TXAS，重点关注其O-H解离动力学。TXAS谱作为反应坐标ROH的函数进行计算，产生STXAS(ω, ROH)，这提供了沿解离路径的连续光谱演化图。这个理论光谱图可以与未来的实验光谱STXAS(ω, τ)直接关联，以构建潜在超快动力学的全面且分辨率高的图像。目标是提供所有可能反应路径的决定性光谱特征，并为未来的实验建立清晰的机理框架。为了解码光谱演化并将其与价态电子结构的变化联系起来，本文采用了自然跃迁轨道(NTO)分析。
价态势能曲线与非绝热路径
计算的S0–S2态沿ROH的态平均完全活性空间自洽场(CASSCF) PECs揭示了两个关键的锥交叉CI1(ROH= 1.10 ?) 和 CI2(ROH= 1.96 ?)，整体形貌与先前报道的MRCISD PECs一致。通过两个锥交叉的电子态特征的演化由NTO分析阐明。在Franck-Condon几何结构，S1和S2分别呈现ππ和π-里德伯特征。拉伸O-H键并通过CI1后，S1态获得主导的σ特征，标志着驱动O-H键断裂的排斥πσ*态的开始。通过CI2后，该态（现在变为S0）显示出主导的氧孤对电子(n)特征，这与PhO•自由基产物的预期电子结构一致。
O1s势能曲线
最低10个O1s激发态(e1–e10)的PECs通过态平均限制活性空间SCF(SA-RASSCF)方法模拟。核心激发态PECs与其价态对应物相比，表现出明显不同的演化模式。这表明核心能级提供了一个独立的、高度局域化的解离过程视图。虽然价态主导核动力学，但O1s核心激发态对氧原子周围的局域电子结构异常敏感，使其光谱位移成为O-H键断裂的直接特征。最低的O1s激发态e1表现出显著的解离特征，在长O-H距离处接近约532 eV。
早期S1动力学(min S0→ CI1)的TXAS特征
沿S1PEC的早期光动力学，跨越II → III → IV，在模拟的TXAS谱中捕获。在初始的II → III区域，出现三个不同的峰(A–C)，每个都表现出特征性的能量位移。峰B是三者中最强的，经历了显著的~2 eV红移(536 → 534 eV)，而较弱的峰A和C显示出~1 eV蓝移(533 → 534 eV 和 538 → 539 eV)。因此，当接近CI1(ROH= 1.10 ?)时，峰A和B在能量上变得更近，标志着在这个关键几何结构处的电子重组。NTO分析阐明了这些特征的电子起源。峰A对应两个潜在的O1s → π跃迁。对于峰B，跃迁特征随O-H伸长显著演化：该特征在较短ROH处源于O1s → 里德伯激发，随着键的延伸，它逐渐转变为主要的O1s → σ特征，伴随着在解离的O-H键上电子局域化的增加。峰C归属于O1s → π跃迁。
早期S1动力学(CI1→ CI2)的TXAS特征
在两个锥交叉之间的区域，跨越III → IV，识别出四个新峰(D–G)。峰D是一个以528–530 eV为中心的显著特征，发生红移并增强强度；它被指定为O1s → π跃迁。峰E出现在532–534 eV范围，对应O1s → π跃迁。峰F在535–537 eV处出现，表现出特征的π特征。峰G在537–540 eV处观察到，被指定为O1s → σ*跃迁；其能量和强度都沿反应坐标减小。当体系接近CI2时，峰F和G合并成一个单一的宽而弱的吸收特征。
TXAS模拟捕获了沿S1PEC的完整动力学演化，包括CI1附近。虽然实际的隧穿过程可能不精确地达到CI1几何结构，但计算的CI1附近的光谱变化可以作为隧穿是否发生（相对于O-H伸缩坐标）的指示器。
后CI2分支路径的TXAS特征
沿三条后CI2路径(Paths 1–3)的模拟O1s TXAS谱揭示了不同的光谱特征。
- •
沿Path 1，一个非常强的峰H(527–528 eV)源于S1→ e1跃迁，沿坐标有轻微红移。其特征被描述为随着ROH增加，从O1s → σ演变为主要的O1s → n特征，反映了氢原子后退留下局域在氧孤对轨道n上空穴的物理图景。大的振子强度源于氧局域跃迁的偶极允许性质。此外，在CI2附近出现一个弱的瞬态特征I，具有π特征，但沿此解离坐标迅速减弱。
- •
对于Path 2，谱图表现出三个持续的特征。一个突出的峰J在528–529 eV，主要归属于O1s → n跃迁，与相邻原子上的少量π特征混合。与Path 1中的峰H相比，J的强度大约减半，并且能量几乎恒定。此外，在531–534 eV有一个具有π特征的较弱峰K，以及在535–536 eV的一个更弱的峰L，主要是π特征，在短距离处有少量σ*混合，两者都保持几乎恒定的能量和强度。
- •
Path 3显示出最动态的光谱演化。一个强峰M经历了约7 eV的显著蓝移(从CI2处的~529 eV到min S0处的~536 eV)，伴随着强度降低。它被指定为具有局域在O-H键上主要σ特征的S0→ e1共振。此外，在更高能量区域，出现两个π介导的特征：峰N从541 eV红移，而较弱的峰O出现在min S0附近。
所有三条路径都在低能量区域显示出不同的、强的特征，这些是路径识别的清晰指纹。
直接S1→ S0内转换后S0*光衰变的TXAS特征
为了完全映射统计衰变的光谱特征，定义了Path 4，代表直接内转换到势能束缚侧的振动热基态S0*，随后是振动弛豫。因此，Paths 3和4代表了相同潜在统计过程的互补初始条件，区别在于分子内振动能量再分配(IVR)和最终解离是发生在预解离构型还是束缚构型。沿Path 4的模拟TXAS谱跟踪了从高度收缩的O-H距离(ROH= 0.6 ?)回到平衡几何结构(0.96 ?)的振动弛豫。识别出四个峰(P–S)，其中P和Q最强。最突出的光谱演化是535–537 eV处的一个宽特征分裂成两个不同的峰(P和Q)，随着ROH增加，它们略微向相反方向移动。棒状谱图揭示分裂源于S0→ e1(混合σ/里德伯)和S0→ e2(π)跃迁之间不断增长的能量分离，S0→ e1强度先减小后增大而S0→ e2强度几乎保持恒定，加剧了这一分裂。在更高能量处观察到两个额外的弱π*峰，R(~541 eV) 和 S(545–547 eV)。
光谱指纹的独特性和实验可分辨性
已识别光谱指纹的诊断能力得到了其显著能量分离的支持，这远超过了典型的展宽效应。峰M(Path 3)经历从529到536 eV的蓝移，导致与Paths 1和2的特征峰(峰H和J，都在528–529 eV附近)最终分离约7 eV。即使是更接近的峰H和J也相隔~1 eV，并且表现出明显不同的强度。这些分离远超过了振动展宽或核心激发非平衡核运动通常引起的亚电子伏特位移。考虑到现代X射线自由电子激光器(XFELs)在O K边常规实现优于0.5 eV的分辨率，这些显著的能量分离确保了光谱特征在时间分辨实验中将保持清晰可辨，不受动力学展宽影响。
额外核坐标的影响
目前呈现的计算集中在O-H解离坐标ROH上，所有其他内坐标固定。虽然这种受限方法有效地捕捉了基本的电子演化，并且对于构建平滑的核心激发态PECs是独特实用的，但进行了额外的计算来评估其他核自由度的影响。具体来说，价态PECs是沿着一个具有弛豫分子几何结构的、更具物理代表性的路径计算的。
关键的是，从弛豫坐标路径获得的价PECs保留了与受限ROH扫描相同的定性拓扑结构。完全优化的锥交叉仅比其未弛豫对应物稳定0.3–0.4 eV，证实了一维未弛豫坐标捕捉了解离的主导能量学。虽然未沿弛豫路径重复完整的核心激发映射，但在优化的锥交叉III'和IV'处计算了TXAS谱。虽然每对之间(III vs III'；IV vs IV')可辨别光谱差异，但整体光谱轮廓保持定性相似。例如，在ROH= 1.72 ?的冻结几何结构处计算的谱图在能量和线形上介于IV和IV'的谱图之间，证明了ROH是控制光谱演化的主要坐标。NTO分析证实了三个结构处相应的跃迁显示出非常相似的价态和核心激发。
总之，额外核坐标的弛豫并没有定性地改变整体PEC拓扑结构或光谱演化。本研究的主要目标是为紫外激发苯酚中的氢原子解离建立路径特异的XAS特征。这里达到的精度水平足以区分所提出的路径。因此得出结论，受限的ROH方法提供了一种有效且实用的手段来捕捉光解离过程的基本特征，并且不影响本研究的主要结论。
精度与效率的平衡
准确预测核心激发/电离能和X射线光谱即使在基态几何结构下也是量子化学中的一个重大挑战。当为光化学反应映射整个PESs时，这一困难加剧，这需要一种方法能够准确捕捉一系列核构型中的电子关联，同时保持计算可行性。对于轻元素(C, N, O)的K边光谱，广泛使用的ΔKohn-Sham(ΔKS) DFT方法提供了一个实用的折衷方案，通常在激发能方面实现0.5–1.0 eV的误差。虽然它可以扩展到Franck-Condon区域甚至略微拉伸的几何结构，成功再现PECs和振动分辨的XAS谱，但其单参考性质在描述完整光化学动力学中遇到的强关联电子态和复杂的非绝热耦合时变得不充分。
多组态方法是这个体系的自然选择。RASSCF方法为捕捉静态关联——光化学过程中控制PES拓扑结构和态交叉的主导效应——提供了一个稳健的框架。虽然它通常在绝对核心激发能方面产生几电子伏特的精度，但它可靠地再现了实验光谱轮廓，并已成功应用于映射PECs和解释光化学动力学。RASPT2方法还通过微扰论包含了动态关联。对于过渡金属L边XAS计算，动态关联通常是再现弱光谱特征所必需的。然而，在有机K边光谱中，光谱形状和非绝热路径主要由静态关联决定。已有证明，对于碳-氮体系，RASPT2和RASSCF产生相似的价态和C1s核心态PECs以及可比的高分辨率XAS谱。
方程运动耦合簇单双(EOM-CCSD)提供了改进的动态关联处理，通常被认为是核心能级光谱建模的参考，其应用已扩展到代表性的激发态结构。同时，MRCI已证明具有高精度(核心激发的平均绝对误差~0.25 eV)，但当前应用主要限于基态几何结构。需要进一步的方法学发展，以充分实现它们在光化学中构建PECs的潜力，特别是在管理计算需求和确保大小一致性方面。
考虑到精度、计算成本和对激发态轨迹适用性之间的权衡，本文采用的RASSCF方法为模拟光化学过程提供了一个平衡的折衷方案。该方法已被证明能再现结构类似系统(如尿嘧啶)的MRCI PECs，并且RASSCF模拟的尿嘧啶S1暗态特征(526–527 eV)在O1s TXAS中与在胸腺嘧啶XFEL实验中观察到的结果非常吻合。具体对于苯酚，RASSCF在基态几何结构产生的谱图与实验吻合良好。本文的RASSCF价态PECs与报道的MRCISD PECs基本相似。因此，认为RASSCF方法对于识别沿光解离坐标的路径特异光谱特征的当前目标既高效又足够准确。
结论
本文系统地模拟了苯酚在紫外激发到S1(ππ)态后沿O-H解离路径的O1s瞬态X射线吸收光谱(TXAS)。研究结果在三个方面取得了关键进展，解决了引言中概述的机理问题。首先，早期S1动力学(min S0→ CI1→ CI2)期间的TXAS演化揭示了CI1处的突然光谱变化，作为ππ→ πσ*跃迁的清晰特征，并能够追踪第一个锥交叉的通过。其次，三条后CI2路径——直接解离到H + PhO•(?) (Path 1)、通过S0快速解离到H + PhO•(X?) (Path 2)和S0重穿(Path 3)——每条都表现出独特的光谱指纹。这些特征允许在未来的TXAS实验中进行明确的机理区分。第三，本文识别了通过在势能曲线束缚侧直接内转换形成的振动热基态S0*(Path 4)的独特TXAS特征。与来自解离侧(Path 3)的特征一起，这些特征为探究难以捉摸的基态统计衰变通道提供了互补的实验手段。
展望
“沿PESs的快照”方法为表示统计平均提供了一个实用的解决方案，补充了非绝热分子动力学在光化学和基态热反应的玻恩-奥本海默分子动力学(BOMD)中捕捉随机事件的明确时间依赖性的优势。此策略对于具有主导一维反应坐标的系统特别有效，正如本文在苯酚O-H裂解的教科书案例中所彻底证明的。类似的1D理念适用于质子转移和其他可以定义主导集体坐标的场景，展示了该方法的广泛适用性。这个概念框架为解决更复杂的光化学问题奠定了基础。
从一维到多维光化学
一个自然的问题是，这种静态方法是否可以推广到具有两个或多个耦合反应坐标的系统。原则上，可以在预定义的网格上计算静态光谱图STXAS(ω, ?, θ)。一个看似直观的简化是定义一条单一、预定的光化学反应路径，并沿其计算光谱。然而，STXAS(ω, ?, θ)数据无法确定唯一的动力学序列。例如，在锥交叉处，动力学通过连续的路径进行，而不是单一的轨迹，使得任何预定路径都不足以表示系综。因此，如果没有额外的动力学输入，无法直接将静态图STXAS(ω, ?, θ)转换为时间分辨光谱STXAS(ω, τ)。
为了克服这一点，同时保持静态光谱计算的准确性，本文提出了一种结合非绝热分子动力学统计采样的关键快照协议。这种混合策略将构象采样与光谱计算解耦：1) 高效NAMD生成轨迹系综；2) 聚类识别主导区域并选择代表性的“关键快照”；3) 仅在这些快照处计算基准光谱，构建一个紧凑的库{Sk(ω)}。然后光谱被重构为STXAS(ω,τ) ≈ ∑kwk(τ)Sk(ω)，权重wk(τ)来自轨迹系综。这种混合策略保留了静态方法的效率，同时结合了系综统计。它克服了纯NAMD方法的主要限制，在纯NAMD方法中，耗时的X射线光谱计算需要在所有轨迹上进行。这个概念框架的实际验证仍是未来研究的主题。
机器学习加速的动力学和光谱计算
机器学习(ML)势为获得统计收敛的动力学参数(如分支比)提供了一条变革性途径。通过在量子化学参考上训练神经网络势，ML加速的NAMD可以以比直接进行量子力学计算低几个数量级的计算成本生成统计收敛的轨迹系综，使得以前无法进行的模拟变得可行。此外，生成光谱库的后续步骤同样可以通过ML加速。加速光谱模拟的方法正在成为一项补充策略。
扩展到复杂环境
该方法可以通过量子力学/分子力学(QM/MM)等嵌入方案扩展到异质环境。正如在基态热化学中所证明的，环境效应——包括来自溶剂、蛋白质口袋或材料基质的环境效应——可以通过具有适当静电嵌入的MM区域纳入。结合ML加速动力学，QM/MM嵌入为将X射线光谱应用于现实环境中的超快光化学过程提供了一条实用途径。因此，光谱图STXAS(ω, ξ)作为将实验STXAS(ω, t)数据转化为机理洞察的实用理论参考。
