本方案旨在为初学者提供处理、拟合和解读瞬态吸收光谱的入门指导。本方案的重点是数据集的准备,以及利用单波长动力学和全局寿命分析进行拟合。文中还讨论了瞬态吸收数据及其拟合过程中相关的挑战。
方法文章
本方案旨在为初学者提供处理、拟合和解读瞬态吸收光谱的入门指导。本方案的重点是数据集的准备,以及利用单波长动力学和全局寿命分析进行拟合。文中还讨论了瞬态吸收数据及其拟合过程中相关的挑战。
瞬态吸收(TA)光谱是一种强大的时间分辨光谱技术,通过监测系统吸收光谱的变化来追踪激发态过程的演化。早期的TA技术仅限于专业实验室使用,但随着商业化即用型系统的不断发展,该技术已逐渐普及至全球各地的研究团队。现代TA系统能够以高能量和高时间分辨率产生富含光物理信息的大量数据。然而,由于存在多种激发态特征和仪器伪影,TA光谱的处理、拟合与解释往往具有挑战性。在采集、处理和拟合TA数据时,必须仔细考虑诸多因素,以降低对何种模型或拟合参数组最能准确描述数据的不确定性。数据预处理与拟合的目标是在保留可用于分析的数据的同时,尽可能消除这些外部干扰因素。本方法为初学者提供了一套TA数据处理与准备的标准化流程,并简要介绍了若干选定的拟合方法与模型,特别是单波长拟合和全局寿命分析。文中还对数据预处理过程中常见的一些挑战及其应对方法进行了评述,随后讨论了这些简单拟合方法所面临的挑战与局限性。
瞬态吸收(TA)光谱是一种时间分辨光谱技术,通过监测样品在光脉冲激发后其吸收光谱随时间依赖的变化,来研究光激发物种的演化过程。由于TA是一种吸收型技术,因此能够识别并追踪那些既发生辐射跃迁(即通常会发射光子的状态)又发生非辐射跃迁(如通常不产生荧光、经历内转换、系间窜越或参与光反应的状态)的能态所产生的光谱信号1,2。根据激发光源和探测方法的具体配置,TA技术可覆盖从飞秒到微秒以上的时间尺度,光谱范围从紫外(UV)延伸至远红外(IR),使其成为一种多功能的光谱分析工具。在过去几十年中,TA光谱仪的商业化进程取得了显著进展,使得越来越多的实验室和科研机构能够使用这一强大技术2。
现代瞬态吸收(TA)系统能够以高能量和高时间分辨率产生大量数据集。这些数据集通常表现为一个二维矩阵,其数值为相对于激发脉冲的时间延迟和波长的透射率或吸光度差值。该数据集可被视作二维热图或三维拓扑图。随着研究人员力求在拟合过程中包含完整的数据集以最准确地描述其研究体系,这些数据的解读变得愈加复杂3。
尽管瞬态吸收光谱(TA)可覆盖广泛的波长范围和时间尺度,本实验方案聚焦于其中一种最容易实现的形式4:由飞秒脉冲激光驱动的紫外-可见光区宽带光谱。此类仪器的示意图见图15,6。实验首先从激光器中取出一个脉冲,并将其分为两束。其中一束称为“泵浦光”,用于激发样品,通常使用光学参量放大器(OPA)等装置将泵浦光脉冲转换为所需的激发波长5,7。另一束称为“探测光”,进入一个机械延迟装置,通过调节光脉冲传播的距离来改变泵浦光与探测光之间的延迟时间。随后,单波长的探测光脉冲通过蓝宝石或氟化钙(CaF2)晶体转换为白光连续谱8。白光脉冲穿过样品后,使用电荷耦合器件(CCD)相机等宽带探测器测量其光谱。通过测量有无泵浦光时白光脉冲光谱的变化,可获得泵浦光诱导的样品吸收光谱变化ΔA(T)。感兴趣的读者可参考这篇有用的综述9以了解有关探测过程的更多细节。
在所有形式的瞬态吸收(TA)光谱中,ΔA(t) 谱图是通过计算探针光的基态吸收 Aprobe 与在两束脉冲之间给定时间延迟 t 时激发态吸收 Apump+probe 之间的差值得到的2,5,9,10。
(1)
需要注意的是,Aprobe 等同于样品的稳态吸收光谱,且与时间无关;实验的时间分辨率来源于泵浦光与探测光之间延迟时间在 Apump+probe(t) 中的体现。这些数据的模拟结果如图2A所示。
与稳态吸收光谱不同,TA光谱由于采用了公式1中的差分方法,可能同时呈现正向和负向特征。正向特征源于泵浦光脉冲产生的新吸收物种,可代表激发态生色团、三重态、几何结构重排、溶剂化效应或激发态光产物3. 在讨论部分将介绍识别这些特征并将其归属到化学物种的一般原则。负特征可能源于基态漂白(GSB)或受激发射(SE)图2B。GSB 是由于泵浦脉冲吸收后基态粒子数减少所致。被激发至激发态的分子不再与其基态在同一区域吸收光;因此, 较少 探针脉冲的吸收导致式(1)中的差值在该区域可能为负值。基态漂白(GSB)的特征是其光谱形状与基态吸收相同,但符号相反。受激态发射(SE)信号来源于激发态物种在探针脉冲作用下产生的受激发射。3这些物种的发射导致 更多光 到达检测器,这相当于具有 吸收减少 在这些波长处,受激发射(SE)信号的光谱形状与该物种的自发发射光谱相似,但符号为负,且具有不同的频率权重。10.
除了有关激发态物种的信息外,瞬态吸收(TA)光谱还可能包含多种伪影和无关特征,这些因素可能会扭曲潜在的动力学过程,并妨碍对吸收带的准确归属11。若在数据预处理和分析过程中未能正确处理这些伪影,则可能导致对数据应用不恰当的光物理模型,从而得出误导性的结论11。因此,本实验方案的第一部分重点介绍在采集TA数据后如何对其进行正确处理。本节的目标是为刚接触TA技术的研究人员提供一套指导原则,帮助他们建立对数据严格准备与处理的直观理解与认识。
数据集处理完成后,可使用多种复杂程度和严谨程度不同的工具和模型对光谱进行拟合与解释10。本方案第二部分的目标是帮助读者掌握如何对数据应用单波长拟合和全局分析,并提供指导,说明在何种情况下这些模型适用于描述其数据。目前已有易于获取的商业软件可用于飞秒瞬态吸收(TA)数据的准备与处理,例如 Ultrafast Systems 公司开发的 Surface Xplorer12,13(可免费下载和使用,参见材料表)。此外,学术研究人员也发布了其他免费替代方案,例如 Glotaran14。Glotaran 是一款免费软件程序,专为时间分辨光谱学和显微镜数据的全局分析与目标分析而开发,其作为 R 语言程序包 TIMP14 的图形用户界面(GUI)。此外,用户也可使用 Python 等编程语言自行编写代码实现数据分析。这些拟合软件和编程解决方案各自具备独特优势,均为该领域的重要贡献。鉴于本研究的目的,我们仅能介绍其中一种软件用于本活动的可视化部分。对各类拟合软件的深入讨论已超出本文范围。
本文提供了一套逐步操作流程,用于(1)处理瞬态吸收(TA)数据,(2)利用单波长动力学和全局分析对TA数据进行拟合,以及(3)提取数据并将其拟合至其他模型。文中包含一组代表性TA数据,供读者练习使用(补充文件1和补充文件2)。该数据为在乙醇中浓度为165 µM的1,4-双(5-苯基噁唑-2-基)苯(POPOP)样品在330 nm激发下测得的信号,采集时间范围为−5 ps至5.5 ns。此外,在相同实验条件下还采集了仅含乙醇、不含样品的“空白”样品数据,时间范围为−5 ps至5 ps,该数据用于拟合前的数据预处理(步骤1)。光谱数据使用超快瞬态吸收光谱仪采集。样品置于光程为2 mm的比色皿中,并保持持续搅拌。所述数据处理与拟合流程基于Surface Xplorer软件,该软件可拟合*.ufs格式的数据,本文中简称其为“拟合程序”。已有程序可用于将其他格式的数据集转换为*.ufs文件15。尽管本方案的具体细节针对Surface Xplorer软件,但后续步骤可推广至任何商业或自主开发的软件包。此外,数据处理的结果也可导出并使用其他软件进行拟合。支持信息文件(补充文件3)提供了有关数据拟合的进一步建议。
1. 拟合前的数据准备
2. 进行拟合
图。
表面。左下图块将显示拟合生成的衰减关联差分光谱(DADS)。最后,右下图块用于设置拟合参数,包括使用的指数函数数量以及是否使用无限函数。
图中无明显特征或仅有极少特征,则可接受该拟合。通过更改寿命数量和/或勾选/取消“使用偏移(Ainf)”按钮,可轻松尝试多个拟合。在检查多种拟合参数变化后,应接受其中最佳的拟合结果。
图或主动力学轨迹拟合的任何信息。主动力学轨迹的保存方法见步骤 3.3。有关如何从原始数据窗口保存数据以供后续显示和绘图的信息,请参见步骤 3.0。3. 从拟合软件中提取原始数据和拟合结果用于绘图
注意:可通过单波长拟合或全局分析生成的原始数据或拟合结果可导出为 csv 文件,该文件可在多种其他程序中打开。
按照上述步骤制备并分析了溶于乙醇的1,4-双[2-(5-苯基恶唑基)]苯(POPOP)样品。使用超快瞬态吸收光谱仪进行测量,如图1所示,液体样品置于2 mm比色皿中,采用可调式比色皿支架和磁力搅拌器以确保溶液均匀混合。样品在常温常压条件下测量,未对温度或气氛进行额外控制。340 nm至680 nm的光学窗口通过氟化钙晶体产生。在−5 ps至约5500 ps范围内采集了250个时间点,共进行三次扫描并取平均值得到最终数据集,如图3所示。POPOP数据的处理按照实验方案所述进行。补充图16展示了 chirp 校正不理想的示例。对POPOP进行了单波长动力学拟合,选择632 nm作为目标波长。此外,还按照实验方案对POPOP进行了全局分析。
在632 nm处对POPOP进行单波长动力学拟合,得到两个寿命。允许这两个寿命自由变化,未作进一步调整。获得的最终参数如下:t0 = −0.1176 ps,IRF = 0.436 ps,A1 = 0.0956,t1 = 1.614 ps,A2 = 0.0646,t2 = 522.2 ps(图4)。这些结果与后续进行的全局分析以及文献报道的POPOP发射寿命值(τ = 1.35 ns)17高度一致。单波长拟合中寿命组分过少的一个示例见补充图16,并附有说明和讨论。
对POPOP进行全局分析拟合时,在执行SVD过程中选择了15个主成分(PCs),拟合后选取了两个寿命参数,且未固定任何参数。拟合得到的最终参数如下:t0 = −0.1586 ps,tp(IRF)= 0.4408 ps,t1 = 1459 ps,t2 = 267.5 ps。衰减关联差分光谱如图5所示。结果与在632 nm处进行单点动力学拟合所得结果以及POPOP的寿命值17高度一致。两个全局分析次优结果的示例见补充图16,并在文中进行了展示和讨论。

图1:本实验方案中描述的飞秒宽带瞬态吸收仪器的示意图。 请点击此处查看此图的放大版本。

图2:TAS数据的模拟及对TAS信号的贡献。(A)高斯曲线模拟的基态吸收光谱(蓝色虚线)和激发态光谱(红色虚线,相对于基态光谱发生红移)。这两条光谱的差值光谱(紫色实线)即为飞秒瞬态吸收(TA)实验中观测到的差分光谱。为便于说明,此处夸大了基态与激发态光谱之间的差异。(B)POPOP在330 nm激发后1.04 ps时具有代表性的TA差分光谱。虚线表示POPOP的基态吸收和稳态发射。高亮区域显示了该数据中常见的TA特征,包括基态漂白(GSB)、受激发射(SE)和激发态吸收(ESA)。请点击此处查看该图的放大版本。

图3:经步骤1“数据准备”处理后得到的POPOP数据。 数据以校正后的热图和代表性光谱形式展示。这些结果展示了在完成校正后数据应有的形态,表明此时已可对数据集进行拟合。 请点击此处查看该图的放大版本。

图 4:在执行步骤 2.3(单波长拟合)后,对 POPOP 在 632 nm 处进行单波长拟合的结果。 该图显示(上部)以表格形式列出的寿命值,(中部)实验数据(蓝色点)及对数据的拟合结果(红色拟合线),以及(下部)残差图。请注意,在“当前拟合”部分中,振幅(A)显示为来自数据的 ΔA 值,代表该特定寿命组分在 t₀ 时刻的贡献。然而,当当前拟合被保存并出现在“拟合系数”表格中时,默认设置为显示归一化振幅。此设置可通过取消勾选“拟合系数”标签旁标注为“归一化”的复选框进行更改。请点击此处查看该图的放大版本。

图5:完成步骤2.4拟合后对POPOP进行全局分析拟合的结果。 请点击此处查看该图的放大版本。

图6:在实验TA数据中观察到的受激拉曼散射示例。 在本数据集(非本教程中展示的POPOP数据)中,样品在550 nm处被激发(由虚线表示)。拉曼散射出现在时间零点附近,通常出现在泵浦激发波长的蓝移侧(反斯托克斯拉曼散射)和红移侧(斯托克斯拉曼散射)。受激拉曼散射寿命较短,通常仅约~200 fs,因为它发生在探针光束与泵浦光束同时与样品相互作用时,从而激发拉曼过程。然而,该特征无法避免,应通过裁剪予以去除。请点击此处查看该图的放大版本。
补充图1:数据首次加载时的主菜单。通常情况下,在此拟合软件中,当鼠标悬停在窗口上时,窗口的右下角会出现选项框。这些选项框可用于更改光标与窗口的交互方式,例如移动选区、缩放或平移选定的窗口。此外,还可以调整窗口的x轴和y轴刻度的精度,并可在线性与对数显示模式之间切换。坐标轴也可进行锁定或解锁操作。当坐标轴被锁定时,其缩放级别或数值范围将保持不变;解锁后,坐标轴范围将自动调整以包含全部数据。其他选项还包括调整数值显示格式以及(如果存在)网格线的颜色。请点击此处下载该文件。
补充图2:执行散射光减除的表面菜单。 请点击此处下载该文件。
补充图3:减去散射光屏幕——设置光谱取平均值。 请点击此处下载该文件。
补充图4:设置裁剪范围。 请点击此处下载该文件。
补充图5:裁剪数据集。 裁剪后的数据将被永久删除。请点击此处下载该文件。
补充图6:拟合溶剂响应窗口及溶剂响应拟合曲线(红线)与实验数据(蓝点)的拟合示例。 请点击此处下载该文件。
补充图7:使用“空白”样本数据的啁啾校正窗口。 请点击此处下载该文件。
补充图8:显示与数据集一起的啁啾校正。 请点击此处下载此文件。
补充图9:单动力学拟合窗口。此外,顶部的滑块允许用户选择用于执行拟合的波长。选择合适的波长需要了解该系统的光谱特性,以确定感兴趣过程发生的区域。例如电荷转移产物、三重态形成或光产物形成,其光谱特征已知与特定波长相关联。可对这些特定波长进行拟合,以获得相应特定事件的寿命。此外,也可选择多个波长用于验证选定的全局拟合模型。请点击此处下载该文件。
补充图10:设置单动力学拟合参数。 如有需要,可手动固定或调整各个拟合参数,以优化拟合效果并最小化残差的标准偏差。注意:可通过单击参数框中的数值,然后使用滑块调节或直接输入数值来轻松调整参数。在修改数值时,显示的拟合曲线将实时更新。当获得满意的拟合结果后,可通过右键单击拟合界面导出拟合结果,支持将数据导出至剪贴板以便粘贴到目标程序中,或以图片形式导出以便快速查看。请点击此处下载该文件。
补充图11:如果一个或多个参数是已知的,则可以对单次动力学拟合参数的值进行约束和固定。 请点击此处下载该文件。
补充图12:单值分解窗口,显示主成分,展示当添加足够多的主成分时,轨迹和各成分的形态。 请点击此处下载该文件。
补充图13:全局拟合分析生成前的全局拟合窗口。 请点击此处下载该文件。
补充图14:显示全局分析拟合结果的全局拟合窗口。 请点击此处下载该文件。
补充图15:用于保存和导出文件的文件菜单。 请点击此处下载该文件。
补充图16: chirp 校正和拟合不佳的示例。(A)展示了一个较差的单指数动力学拟合。A图中所示的结构化残差通常表明需要引入额外的寿命分量才能拟合数据。请注意,残差在较短的时间点低于零线,而在较长寿命区域则高于零线。(B)展示了使用预览按钮进行的正确 chirp 校正,其中特征线被拉直,数据中无任何弯曲。(C)展示了一种错误的 chirp 校正,其中光谱的蓝色区域存在明显弯曲,表明 chirp 函数在该区域过度校正。(D)展示了一种较差的全局寿命分析拟合,其中过度拟合(参数过多)导致出现“对称”的差分吸收动力学谱(DADS),这些谱图在x轴上下互为镜像,寿命范围相似,实际上相互抵消。当观察到此类特征时,不应采用该拟合结果。(E)展示了一种较差的全局拟合,其中参数过多导致出现一个振幅极大但寿命极短的组分。图(E)中的问题也可能在零点附近伪影未得到适当校正时出现,此时拟合过程会过度强调一个无实际物理意义的极短寿命,以最小化残差。请点击此处下载该文件。
补充图17:包含拉曼散射的瞬态吸收光谱(TAS)数据示例。 在时间零点附近存在散射,且与泵浦激发波长重合。散射表现为一系列尖锐的峰,泵浦激发波长蓝移处有一个非常强烈的正峰,红移处有一个负峰。该特征无法合理避免,应从数据中裁剪去除,以避免干扰拟合结果。请点击此处下载该文件。
补充文件 1:包含本教程数据集的文件(POPOP data_POPOP-inEtOH.ufs)。 请点击此处下载该文件。
补充文件 2:本教程所用的空白数据集文件(POPOP data_BLANK.ufs)。 请点击此处下载该文件。
补充文件3:包含有关溶剂响应拟合、散射校正和表面扣除的附加说明的支持信息文件。 请点击此处下载该文件。
数据准备的一般注意事项
时间分辨吸收(TA)数据的拟合乍看之下似乎相对简单直接,人们可能期望对于给定的数据集会得到一个明确且正确的“答案”。然而,正如本实验方案中所强调的,数据采集、数据准备和数据分析过程中存在许多需要仔细考虑的因素,这些因素可能导致难以确定哪一个模型或哪一组拟合参数最能准确描述数据。数据准备与拟合的目标是在保留可用于分析的数据的同时,尽可能减少这些外部因素的影响。对于初学者而言,这项任务可能看似艰巨,因为需要考虑的内容较多。为了建立对拟合过程的直观理解,建议初学者尝试以略微不同的方式多次从头开始准备同一组数据,以检验数据准备步骤对最佳拟合结果的影响程度。此外,两位不同的研究人员可以分别准备并拟合同一组数据,然后比较结果。这一过程在最初几次可能较为耗时,但这样做将帮助初学者逐步建立起一致地处理未来样品数据的能力。如同任何技能一样,数据准备与拟合需要时间来掌握,建议初学者在探索和学习该过程时保持耐心和严谨。本研究中使用的数据集已提供,以便初学者能够跟随教程进行实际拟合,并直接将自己的结果与教程中的结果进行比较。
数据可能包含 背景特征 在所有时间延迟下均存在的补充图2 和 补充图3)例如泵浦光束的散射以及样品的自发发射。必须去除这些干扰因素,才能分离出目标物的瞬态吸收信号11通过选择、平均并去除若干负时间延迟光谱的贡献来实现此类特征的消除。在选择背景光谱时,必须确保不包含任何可能属于目标过程本身的特征,以免误删。在瞬态吸收(TA)数据中,还可能观察到由溶剂引起的相关背景信号,例如杂质或溶剂本身的吸收。当溶剂产生信号时,需从样品数据中减去一份“空白”数据集,该数据集应在与样品完全相同的实验条件下仅使用溶剂单独测得。有关此操作的具体细节包含在 补充文件 3.
另一个需要仔细考虑的因素是啁啾校正。当探测脉冲传播到样品时,由于转向镜的不完美或穿过色散光学元件(如透镜或滤光片),脉冲会变宽,从而产生啁啾效应。最终结果是探测脉冲中能量较低的光子(即探测光谱的红端)比能量较高的光子(即探测光谱的蓝端)更早到达样品。这导致瞬态吸收(TA)光谱的“零时间点”在数飞秒或皮秒范围内被展宽18,在原始数据集中表现为从蓝光波长区域开始并向红光方向逐渐趋于平缓的独特曲线(补充图7)。啁啾效应在超快TA所能达到的较短时间尺度上尤为明显。这种波长依赖的零时间点偏移可根据本方案所述方法进行校正,但该过程的实施可能较为复杂且具有主观性。使用“空白”样品或溶剂克尔响应的测量结果,可以降低手动选择用于拟合多项式以校正啁啾的点所带来的主观性。啁啾校正的目标是消除零时间点处明显的“曲线”特征。可能需要多次尝试拟合啁啾,才能获得最佳的啁啾校正数据。可通过应用不同的啁啾校正进行多次拟合,以理解啁啾校正对短寿命TA组分数值的影响。
在“零时刻”出现的伪影
在瞬态吸收(TA)数据中,接近“零时刻”时可观察到多种伪影,包括瑞利散射、受激拉曼散射和交叉相位调制。泵浦光束的瑞利散射是一种弹性散射,不伴随能量变化,该特征出现在与泵浦脉冲相同的波长处。受激拉曼散射可能伴随泵浦散射信号同时出现19。拉曼散射源于泵浦光子的非弹性散射,会在入射泵浦光能量的高能端(反斯托克斯峰)和低能端(斯托克斯峰)产生峰信号。在TA数据中,由于样品同时受到泵浦光束和探测光束的照射,当探测光束与泵浦光束同时与样品相互作用时,会激发拉曼过程,因此受激拉曼散射出现在零时刻附近,并在前几百飞秒内的光谱中产生额外的峰信号(图6,在高亮区域中较深蓝色光谱中可见,以及补充图17)。交叉相位调制则源于溶剂折射率在强脉冲电场作用下的调制。
受激拉曼散射与交叉相位调制的区别在于,拉曼峰出现在对应于溶剂振动模式的特定频率处。由于这是一种拉曼过程,因此在激发光两侧均可观察到斯托克斯线和反斯托克斯线。氯化溶剂(如二氯甲烷)由于氯原子具有较大的极化率,会表现出非常显著的拉曼谱带。交叉相位调制的光谱特征对特定溶剂具有唯一性,但相比拉曼散射特征而言更难以预测。
根据所测样品的动力学特性,瑞利散射、拉曼散射和交叉相位调制可能与瞬态吸收(TA)数据的早期特征重叠,且难以从数据中完全去除。原则上,这些特征可在纯溶剂测量中观察到,并从数据中扣除;数据分析程序可能包含用于拟合这些特征的函数,但在实际操作中往往较为困难。当去除这些伪影会严重影响样品数据时,更合适的做法可能是裁剪掉时间零点附近的受影响光谱,以消除伪影。这样做虽会不可避免地丢失时间零点前后约300飞秒(fs)的数据,但可提高后续拟合的可靠性。在分析多个相同或不同样品的数据集过程中,初学者将逐渐掌握在扣除背景信号与裁剪初始100–200飞秒数据之间取得平衡的实践经验。
通用裁剪 可能是必要的,用于处理信噪比较低的光谱部分。探针光束在某些区域的不稳定性、探针光强度过低、样品浓度过高(从而阻挡了大部分入射探针光)、泵浦光强度不足以及样品的吸收截面,都是导致信噪比低下的常见原因,这些因素会使数据拟合变得困难。在这种情况下,通过在光学窗口两侧裁剪数据集,以达到所需的信噪比水平,有助于改善拟合过程。
当数据集经过充分裁剪以去除质量较差的部分、完成啁啾校正,并对背景光谱进行平均和扣除后,即可用于分析。该步骤应确保数据仅包含与所关注的光物理和光化学过程最相关的部分。显然,此过程在一定程度上具有主观性。数据预处理的目标是平衡去除伪影以避免其干扰拟合结果,同时避免过度去除而损害数据集的完整性,从而影响对数据的正确解读。找到这一平衡点需要时间和经验的积累,以逐步建立区分伪影与真实信号的直觉判断能力。在不同日期多次对同一组数据进行拟合(或重新拟合),或由两位研究人员分别拟合同一组数据,有助于减少人为误差以及数据预处理和分析过程中的主观性。
拟合与解释的一般注意事项
在对原始瞬态吸收(TA)光谱进行处理后,必须对其进行解释和建模,以提取有关体系中存在物种及其动力学过程的信息。该过程可分为三个步骤:初步光谱解释、定量建模/拟合,以及将光谱解释结果归因于模型/拟合结果。
初始光谱解析:在光谱解析步骤中,目标是将瞬态吸收(TA)光谱中出现的特征归属于系统在光物理或光化学演化过程中所涉及的电子态。首先,应识别出各种不同的态。在本研究中,态指的是参与系统光物理或光化学演化过程的独特电子态。一个态可由一条特定的势能曲线(PEC)表示,具有代表其吸收光谱的一组特征峰。发生在单一态内部的变化称为过程。光物理过程在TA光谱中可能表现为峰位移动或光谱宽度的变化。过程的关键特征是该态的布居数保持不变(即过程发生在给定的PEC内部);发生变化的是态内部的能量分布。若某一态的布居数发生变化,则称之为跃迁。在跃迁过程中,系统演化至另一条PEC(即电子态)。跃迁可能包括内转换(IC)、系间窜越(ISC)、电荷转移、能量转移、新产物的生成,或返回基态。下文将讨论对态、过程和跃迁进行指认的指导原则。
指认能态
该过程的第一步是将光谱特征指认为特定的化学物种或状态。瞬态吸收(TA)中的 S1 态应表现出与利用时间分辨发射光谱测得的荧光寿命相匹配的寿命。若三重态的存在受到怀疑,可通过氧气淬灭其寿命来加以验证。如果在光物理演化过程中怀疑存在自由基阴离子或阳离子,可进行光谱电化学或化学氧化/还原反应以生成相应的自由基物种,并获取该物种的吸收光谱,进而与TA谱带形状进行比对。可通过电子顺磁共振(ESR)光谱实验验证自由基的存在。美国化学会无机化学分会主办的一场优秀教学讲座对TA技术及在指认光谱特征时需考虑的因素进行了全面概述20。在将谱带指认为相应物种后,解释TA光谱的下一步是定性描述体系中发生的动力学过程。此步骤至关重要,因为它可帮助研究人员判断哪些模型适用于描述其体系,并为其拟合参数提供一个比较的基准。
单态内的变化
振动冷却、几何构型重排或溶剂化是极快的过程(亚皮秒到数十皮秒),可通过瞬态吸收光谱(TA)进行观测。振动冷却在数皮秒时间尺度上表现为TA光谱的快速蓝移21,22,23。几何构型重排可在数十皮秒的时间尺度上发生。在常规偶极液体中,溶剂化动力学表现为光谱在数皮秒内的红移和谱带变窄;但在高黏度溶剂(如甘油、聚乙二醇(PEG)、离子液体和深共晶溶剂)中,溶剂化动力学可能持续多个纳秒24,25,26。
状态种群的变化
反应的特征表现为谱带强度的变化,其中强度降低对应于该化学物种浓度的降低,反之,强度增加则对应浓度升高。vice versa 在某些情况下,反应物和产物物种在光谱中均可观察到;而在其他情况下,产物状态存在时间过短或红移过大而无法观测。通常,通过光谱中等吸光点(isosbestic point)的存在可以观察到状态到状态的转变。
定量建模/拟合:必须对数据进行模型拟合,以提取有关系统动力学的定量信息。如引言中先前所述,可选用的模型种类繁多。本方案重点介绍两种最常用的方法:单波长拟合和全局分析。单波长方法涉及将光谱中单个波长的时间轨迹拟合为某种函数形式,通常为指数函数之和:
(2)
其中,ΔA(t) 表示在选定波长下的 TA 信号,n 为指数组分的数量,ai 是第 i 个指数组分的幅值,其时间常数为 τi。可逐步增加组分数量,直至拟合结果能够再现实验数据。任何拟合过程的目标都是使用足够数量的寿命组分以良好地再现数据,同时避免因组分过多而导致过拟合。因此,常采用加权拟合优度参数,例如
,以帮助判断拟合结果是否在实验不确定度范围内符合数据5。
当衰减过程被充分拟合后,模型参数可用于表征系统的动力学行为。随后可提取并解析所得的时间常数。然而,由于瞬态吸收(TA)光谱中存在大量重叠特征,光谱中的单一波长可能包含多种不同物种的动力学信息,这些物种的光谱特征相互重叠,导致从单波长拟合中提取的时间常数可能代表多个同时发生过程的复合结果。此外,任何谱带形状和位置的变化也会影响从单波长拟合中获得的振幅和时间常数。在某些情况下,可通过一种称为“谱带形状分析”(bandshape analysis)的拟合方法来规避上述问题,该方法需确定或假设系统中每种吸光物种的TA谱带的函数形式。这些谱带形状随后由随时间变化的振幅加权,并叠加以重现观测到的光谱。该方法常用于时间分辨荧光光谱的分析,但由于TA谱带形状更为复杂且组分重叠严重,此方法仅在少数简单情况下可行,具体细节见其他文献10。
单波长拟合的另一个缺点是,它无法内在地利用现代瞬态吸收(TA)实验所提供的宽光谱范围。原则上,可以系统地对光谱中的每个单独波长进行拟合,但这种分析方法繁琐、耗时且计算成本高昂。为应对这一挑战,可采用一种称为“全局分析”的方法,同时对整组TA光谱进行拟合并共享一组动力学参数4。全局分析以及一种密切相关的方法——目标分析——是成功且被广泛使用的方法,但它们也各自具有一系列独特的缺点和局限性。与任何模型一样,必须清楚理解构建该模型所依据的假设及其所带来的局限性。
在全局分析中,瞬态吸收光谱(TA光谱)用一个 m × n 矩阵表示,其中 m 代表每条光谱中测量的波长数量,n 代表采集的时间点数量。该矩阵随后被假设可分解为另外两个矩阵的乘积:
(3)
其中 C(t) 是一个 n×k 矩阵,S(λ) 是一个 m×k 矩阵。k 表示用于重构光谱的不同光谱组分的数量。每个组分代表一种具有独特光谱特征和动力学行为的吸光物质。S(λ) 矩阵表示这 k 个组分的瞬态吸收(TA)光谱,C(t) 表示它们随时间变化的浓度。在全局分析最简单且最常见的实现方式中,假设每个组分具有单指数动力学行为(即公式 2 中 i = 1),且每个组分被赋予独立的时间常数。综上所述,完整的 TA 光谱可表示为 k 个光谱组分的叠加,每个组分具有其特有的吸收光谱和单指数衰减行为。
当对瞬态吸收(TA)光谱进行拟合时,用户需先估计所需组分的数量(即k的取值),并初步猜测这些物种单指数衰减对应的时间常数。拟合程序随后生成Cguess(t),并求解方程3以得到Sfit(t)。接着,按照方程3将Sfit(λ)与Cguess(t)相乘,得到拟合后的光谱ΔA(λ,t)fit。最后,通过最小化残差ΔA(λ,t)exp − ΔA(λ,t)fit,返回最优的Sfit(λ)和时间常数。全局分析方法的相对简洁性在于,它使用少量时间常数和固定的光谱组分来表征整套光谱数据,因此成为解析TA光谱中复杂谱带形状和动力学行为的一种有效且成功的方法。然而,必须谨慎判断全局分析是否适用于当前研究体系。
全局分析中的一个关键假设(如公式3所示)是动力学中波长部分与时间部分完全可分离,这一特性称为“双线性”。该假设要求各组分的谱带形状不随时间变化(即它们具有固定的光谱形状,不随时间发生改变或位移)。实验过程中唯一变化的是各组分的相对布居数,由C(t)表示。在较长的时间尺度上(约~1 ns),这一假设通常成立,因此可放心使用全局分析。然而,在飞秒瞬态吸收(TA)可探测的超快时间尺度上,激发态过程(如振动冷却和溶剂化动力学)会导致物种光谱特征随时间变化,从而破坏双线性。这并不意味着全局分析无法拟合数据集——事实上,只要使用足够多的组分数目,总能得到令人满意的拟合结果。问题在于如何解释各组分光谱,并将时间常数归因于特定的激发态过程,因为此时各组分可能已不再对应于明确的吸光物种。因此,在无法假设双线性成立的情况下应用全局分析时,必须始终谨慎对待。
将光谱解释指派给模型或拟合:一旦获得拟合结果,必须将光谱解释映射到拟合所得的寿命参数上。拟合得到的寿命应对应于在光谱初始解释中识别出的各种过程和反应。然而,光谱的初步分析结果与模型拟合所得的寿命数量可能无法立即一一对应。在这种(常见!)情况下,拟合者需要回溯并重新评估最初的光谱解释。也许在初始分析中遗漏了诸如振动冷却或其他过程,但在建模和拟合过程中被识别出来;或者可能存在两组不同的拟合参数均能较好地重现数据,此时初始的光谱解释可指导选择哪一组拟合参数更为合理。在此最后阶段,拟合者必须在光谱解释与拟合结果之间反复比对,以获得对体系中物种及其动力学行为的合理光物理指认。此外,也可尝试其他包含顺序拟合模型的拟合程序(例如目标分析),以补充本文所介绍的全局分析方法及拟合软件所提供的拟合结果4。
总之,本方案讨论了瞬态吸收数据的准备与拟合方法。其目的是强调该过程中相关的挑战,并就如何以实际可行的方式避免或减轻这些挑战进行说明。与技术领域中大多数数据拟合一样,瞬态吸收数据的拟合可能具有难度,有时还带有主观性。因此,充分了解数据本身的特性与局限性、数据预处理过程,以及用于建模和赋予数据意义的数学工具,至关重要。科研人员必须以批判性的眼光对待数据与建模工作。
可以尝试降低拟合结果的主观性。例如,可以从不同的起始点、在不同日期对数据进行预处理和拟合,以确保得到一致的拟合结果。还可以比较不同日期、不同样品制备条件下获得的数据。多位研究人员可对同一组数据进行拟合,并相互比较结果。随着时间的积累,研究人员能够基于自身实验装置和实验参数的具体情况,逐渐形成对所获数据的直观理解,从而对其拟合结果更有信心。
关于瞬态吸收(TA)数据拟合以及本文所讨论模型的细节,有许多内容需要学习。我们强烈推荐几篇优秀的综述文章,这些文章深入探讨了这一主题4,10,27。本实验方案旨在为初学者提供一个进入数据分析与拟合过程的入门途径,激发人们对深入理解该过程的兴趣。
作者无任何利益冲突需要披露。
本研究工作得益于美国国家科学基金会(NSF)重大研究仪器项目的支持,该项目建立了用于瞬态吸收的多用户激光设施(项目编号:CHE-1428633)。 本材料基于美国国家科学基金会资助项目(项目编号:CHE-2313290)支持下开展的研究工作。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 200% 优级纯乙醇(EtOH) | Decon Laboratories Inc | CAS 64-17-5 | 用于样品制备的溶剂 |
| Helios 瞬态吸收光谱仪 | Ultrafast systems | https://ultrafast.systems/products/spectrometers-accessories/helios/ | 瞬态吸收光谱仪 |
| POPOP 1,4-双[2-(5-苯基恶唑基)]苯 | 东京化成工业株式会社(Tokyo Chemical Industry) | CAS 1806-34-4 | 示例中所用样品 |
| Surface Xplorer | Ultrafast systems | https://ultrafast.systems/products/spectrometers-accessories/surface-xplorer/ | 拟合程序 |
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可