《Methods and Protocols》：A Reproducible Two-Stage Protocol for Global Fitting of Time-Resolved Fourier Transform Infrared Difference Spectra
编辑推荐：
摘要专业翻译
时间分辨光谱的全局分析（global analysis）在所有光谱通道上拟合共享的动力学参数，但非线性拟合可能依赖于起始值和停止设置。本研究描述了一种可复现的两阶段协议UFiτ（UFit），用于时间分辨傅里叶变换红外（FTIR）差谱。第一阶段（S
摘要专业翻译
时间分辨光谱的全局分析（global analysis）在所有光谱通道上拟合共享的动力学参数，但非线性拟合可能依赖于起始值和停止设置。本研究描述了一种可复现的两阶段协议UFiτ（UFit），用于时间分辨傅里叶变换红外（FTIR）差谱。第一阶段（Stage 1）拟合峰-峰变化最大的M个通道，以估计共享寿命和拉伸指数；这些估计值用于初始化第二阶段（Stage 2），即对全矩阵进行可变投影（variable projection）拟合。研究人员利用6940个参数已知的独立生成矩阵，将该top-M初始化方法与固定通用起始值的直接拟合、奇异值分解（singular value decomposition，SVD）辅助预拟合、均匀及随机通道子集以及十次随机起点最优法进行了比较。各方法起始值误差相差可达约二十倍，但在大多数条件下最终恢复结果相似。在一项预先声明的近邻寿命测试中，top-M和SVD初始化分别在200个矩阵中的118个和113个中恢复了全部组分（p = 0.27）；在探索性比较中，top-M在十二种条件中的两种条件下恢复的矩阵数较少，且排序可能遗漏一个弱晚期组分。终点一致性主要取决于停止设置：采用缩放残差和更严格的容差后，直接拟合、连续拟合、SVD辅助拟合和top-M拟合对一个实测光系统I记录均达到共同的数值终点屏幕；而均匀间隔和随机子集预拟合则以略低的目标函数值终止于一个不同的、未确定的慢组分。研究人员提供了该协议、组分解析检查清单以及基准代码和结果。
论文解读
一、研究背景与问题
时间分辨光谱（time-resolved spectroscopy，TRS）通过跟踪扰动后光谱可观测量的时间演化来解析化学和物理过程，广泛应用于光解、混合、温度跳变、电化学阶跃和场切换等实验场景。当多个过程同时或依次发生时，单个光谱通道的动力学迹线通常是多个衰减和振幅贡献的叠加。全局分析（global analysis）通过在所有光谱通道上拟合共享动力学模型，同时允许各通道具有独立的振幅，为解析组分提供了一种规范方法。拉伸指数弛豫函数可追溯至1854年Rudolf Kohlrausch对莱顿瓶放电的描述，后被Williams和Watts推广至介电弛豫和聚合物物理领域，即Kohlrausch–Williams–Watts（KWW）函数。在物理上，指数函数是线性弛豫方程的自然解，而拉伸指数则源于微观速率或弛豫时间的分布，反映了异质环境、弥散输运或分布式俘获/复合路径。
尽管全局拟合概念上简单，但在常规使用中存在诸多困难，包括：优化随光谱通道数增加而变慢或不稳定；收敛到局部极小值或近重合/未解析组分；结果依赖起始值；部分软件的实际限制促使光谱子采样和重复试拟合。这些问题在共享参数为非线性（如拉伸指数、分布动力学或耦合速率模型）时尤为突出。选取具有清晰时间演化的动力学迹线并以其拟合时间常数作为全局分析的起点是既有实践，可变投影（variable projection）方法也已在常用工具中实现。然而，该实践缺乏完整规格和系统性验证。本研究的目的在于将这一实践的一个版本完全具体化并加以测试：UFiτ（UFit）按峰-峰时间变化对光谱通道排序，拟合排名前M的通道以估计共享寿命和拉伸指数（第一阶段），再以这些估计作为全矩阵可变投影拟合的初值（第二阶段），并记录所有设置、边界和停止准则。研究进一步通过配对基准测试，评估排序通道初始化是否改变拟合结果。
二、研究内容与主要结论
研究人员开发了UFiτ两阶段初始化协议，并对其进行了系统基准验证。第一阶段根据各通道在拟合窗口内的峰-峰振幅计算“动态性”得分并排序，选取前M个通道（默认M=50）进行预拟合，获得共享非线性动力学参数的初始估计。第二阶段以该估计初始化全矩阵的可变投影拟合，对每个试验参数集通过条件加权线性最小二乘重新计算所有通道的振幅。研究强调，初始化不能为数据添加信息、不能补偿不充分的动力学模型，也不能使超过观测窗口的寿命变得可识别；它只能影响优化器到达哪个局部解以及所需计算量。
通过6940个独立生成的合成矩阵，研究人员比较了top-M初始化与直接固定通用起点拟合、SVD辅助预拟合、均匀/随机通道子集以及十次随机起点最优法。结果表明，起始值误差在不同方法间相差可达约二十倍，但最终恢复结果在大多数条件下相似。在预先声明的近邻寿命测试中（2400个新矩阵），top-M与SVD的恢复率无显著差异（118对113，p = 0.27），均匀和随机子集表现相当。探索性比较中，top-M在特定条件下恢复更少，且排序可能遗漏弱晚期组分。停止设置是决定初始化方法是否达成一致终点的关键因素：以原始ΔA单位应用默认容差时，不同初始化的拟合在不同点停止；采用mOD缩放残差和更严格容差后，直接、连续、SVD辅助和top-M拟合在实测数据上达到共同终点。研究还发现，当光谱密度增加而M固定时，top-M选择的通道可能完全遗漏某一组分，导致初始化质量下降。
三、关键技术方法
本研究采用的主要技术方法包括：其一，可变投影算法，用于可分离非线性最小二乘问题，使优化器仅作用于共享非线性动力学参数，振幅通过条件线性最小二乘获得；其二，SVD辅助预拟合，对前K个时间奇异向量（按奇异值缩放）从通用向量拟合；其三，top-M通道排序预拟合，按峰-峰变化选取高动态性通道；其四，拉丁超立方采样多起点优化（十次随机起点取最优）；其五，合成数据基准测试框架，独立指定寿命、拉伸指数和光谱，生成带高斯噪声或一阶自回归结构化噪声的模拟矩阵。实测数据为77 K下Bruker Vertex 80 FTIR光谱仪采集的时间分辨步进扫描FTIR差谱，样品为menB缺失突变体光系统I（PSI）蛋白复合物，掺入叶绿醌（PhQ），532 nm纳秒激光激发，时间分辨率6 μs，覆盖3.5 ms窗口。
四、研究结果
4.1 实测数据集。 研究人员分析了代表性的TRSS FTIR差谱数据集（331个通道，1800–1120 cm?1，0–1200 μs）。早期时间平均谱呈现尖锐结构特征，后期振幅较小且演化较慢。某些光谱区间的吸收变化携带不成比例的大振幅时间变化，这正是第一阶段排序的直观依据。
4.2 起始值不同但最终结果大多一致。 固定通用起点的最差组分寿命误差在分离寿命场景为0.44、近邻寿命场景为7.3，而SVD辅助和top-M起始值的中位误差分别为0.049–0.053和0.39–0.40。然而，全拟合后的最终误差几乎不可区分。在分离寿命场景中，直接、SVD、top-M和十次随机起点策略分别恢复了100个矩阵中的91、91、87和85个；近邻寿命场景中均为16–21个。在强早期/弱晚期组分场景中，top-M起始值明显更差（中位误差0.86对SVD的0.25），最终误差也更高（0.38对0.23），与所选通道遗漏该组分一致。
4.3 排序未提供一致优势。 预先声明的近邻寿命复制研究中，top-M和SVD分别恢复了200个矩阵中的118和113个（p = 0.27）。均匀子集恢复125个，随机子集恢复117个。没有任何支持top-M优于均匀或随机选择的比较在Holm校正后仍显著，而一项比较反而支持随机选择（缩放寿命、替代光谱、0.5%噪声条件下，top-M恢复99个而随机恢复132个，Holm校正p = 4.6 × 10??）。
4.4 停止设置决定初始化方法是否一致。 采用已发表设置时，不同起始值的拟合在不同点停止，寿命差异最高达约9%。采用推荐设置后，直接、连续、SVD和top-M拟合在每次重复中均达到相同终点；均匀和随机预拟合则以略低的目标函数值终止，主要差异在未确定的慢组分。top-M拟合耗时0.190 s，比直接拟合（0.233 s）少18.5%，SVD辅助拟合为0.119 s。
4.5 menB+PhQ数据集的全局分析。 应用协议于图1数据（K=3，独立β）得到τ? = 12.19 μs（β? = 0.952）、τ? = 336.8 μs（β? = 0.821）、τ? = 2.47 ms（β? = 1.000）。12.19 μs组分暂归因于激光闪光后的样品加热或温度平衡；336.8 μs组分与P700?A??自由基对重组一致；A3为唯象慢贡献，其寿命超过拟合窗口，标记为未确定（n.d.）。
4.6 具有加热样谱形的重现性快组分。 A1的DAS与单独测量的PSI温度诱导差谱高度相关（归一化后相关系数0.87）。在八个存档TRSS数据集中均发现具有相似谱形的快组分，快寿命范围为约9至21 μs。八个归一化谱的平均标准误在每个波数处低于0.02，主负带达?0.30，表明这些异质存档曲线具有共享的谱形特征。
4.7 检查清单应用于典型记录。 将检查清单应用于menB+PhQ记录：M=16、50和100给出相同寿命（差异≤0.02 μs）；窗口延长至3378 μs时A1稳定（11.86 μs），A2从337延长至384–386 μs，A3移至120,000 μs上界；缩短至600 μs时拟合塌缩为两个近重合的约6 μs组分；排除第一个时间点后快寿命变为6.62 μs（排除两个点为8.24 μs），A1谱形保持（相关系数0.95–0.96）。残差的时间滞后一阶相关为?0.08（1200 μs窗口），相邻通道间相关为0.83–0.89，与测量噪声结构一致。
五、讨论与结论
讨论部分指出，初始化显著改变了起始值（近邻寿命时中位最差组分寿命误差相差约二十倍），但很少改变最终结果。在共享可变投影目标、共同边界和充分停止准则下，直接、SVD辅助和top-M初始化通常收敛到相同解；恢复不佳时对所有策略均不佳。两个有实际影响的发现均涉及协议而非起始值：原始ΔA单位的停止准则导致不同初始化在不同点停止，缩放残差和收紧容差后各方法达到共同终点；峰-峰排序可能遗漏弱或光谱分离的组分，尤其在M固定而光谱密度增加时。因此，本研究的贡献在于一个完全规格化的协议而非更优的优化器。与pyglotaran相比，UFiτ更窄，仅拟合单矩阵的衰减关联模型，无仪器响应卷积或目标分析；其价值在于初始化程序确定且可记录，而非更准确。协议原则上可推广至其他时间分辨光谱，但受两个条件限制：无仪器响应卷积，以及排序假设所选通道代表所有组分。结论部分总结：UFiτ是一个完全规格化的两阶段全局拟合协议，基准测试表明预拟合通常产生比固定通用起点更好的起始值，但最终恢复与直接、SVD辅助和未排序子集初始化相似；最清晰的实际因素是停止准则和排序可能遗漏组分的风险；应用于PSI步进扫描FTIR数据时，协议识别出一个具有稳定加热样谱形的快组分、展示了八个存档记录中的共享快谱形，并识别了慢组分不可确定的情形。
