跳到主要内容
Sandbox Physics
从振动走向二维光谱

L4 · 第 8 章

读懂中心线斜率与频率相关

打开 Fenn 与 Fayer 2011 年发表于 Journal of Chemical Physics 的论文,135 卷 074502。先看图 1 的切片方向,再看图 2 的时间曲线,最后读图 3 为什么把混合物问题变难。我们复算式 (12) 的单组分 Gaussian 关系,并用自定合成参数检验它;水的实验数据不在这次复现范围内。

图 1 先找横轴、纵轴与切片

论文图 1 位于印刷页 074502-2。横轴是激发频率 ωτ\omega_\tau,纵轴是探测频率 ωm\omega_m,图中数值用倒厘米表示。两幅图对应等待时间 0.2 ps0.2\,\mathrm{ps} 与 2 ps2\,\mathrm{ps}。图注将它们标为计算的二维谱,采用水的已知频率相关函数。

实线标出沿探测轴取切片的方向,虚线连接各切片的最大值。晚时的峰更接近圆形,虚线更平。沿着正确方向切片,才得到文中使用的中心线斜率;另一类文献会交换切片方向再取逆斜率。

对应阅读:Fenn and Fayer

图 2 的起点为什么低于一

图 2 在下一页,把各等待时间的中心线斜率连成衰减曲线。短时斜率已经明显低于一,文中把这与较大的均匀展宽贡献联系起来。最快的涨落会进入运动窄化部分,不能把图上第一个点机械地补到一,再按单条指数读出所有涨落强度。

接着读式 (1) 与式 (4):作者把频率相关分成快贡献与若干指数,再通过双重时间积分联系到线形函数。中心线曲线提供等待时间信息,线性吸收谱还帮助确定快贡献。两个观测各有作用,少一个会改变可推断的参数。

对应阅读:Fenn and Fayer

完整复算:式 (12) 的单组分极限

式 (12) 给出特定单组分 Gaussian 近似下的关系:沿探测轴找最大值所得的中心线斜率等于 C(T)/C(0)C(T)/C(0)。第七章已经从条件 Gaussian 最大值推导过。现在选标准差 18 cm−118\,\mathrm{cm}^{-1}、相关系数 0.60.6,横轴偏移取 10 cm−110\,\mathrm{cm}^{-1}。预测纵向最大值在 6 cm−16\,\mathrm{cm}^{-1}。

独立数值检查不调用响应函数中的斜率公式:直接在二维 Gaussian 密度上寻找每个切片最大值,再对最大值连线求斜率。随后用另一份频率积分检查有限 Gaussian 频谱权重如何改变协方差。这两种计算检验不同的步骤,原始合成记录、随机种子和版本则允许重建整个信号系列。

再加一个独立的纵向 Gaussian 展宽,标准差为 6 cm−16\,\mathrm{cm}^{-1}。若它只增加探测轴方差,固定初始频率时的条件均值斜率仍由原有横轴方差决定;若两根轴都加同样的独立展宽,斜率才变为 0.6×324/(324+36)=0.540.6\times324/(324+36)=0.54。这说明校准必须先讲清哪个过程作用在哪根轴上。

CLS=Cov⁡(x,y)Var⁡(x)=C(T)C(0)\mathrm{CLS}=\frac{\operatorname{Cov}(x,y)}{\operatorname{Var}(x)}=\frac{C(T)}{C(0)}
最后一个等号采用理想等边际、单组分 Gaussian 频率对。实际窗函数、有限脉冲与其他路径可能改变观测关系。

对应阅读:Fenn and Fayer

图 3 把一个峰拆成两个环境

图 3 位于印刷页 074502-5。黑线是反胶束中 HOD 的 OD 伸缩吸收,蓝线与红线分别表示核区与界面组分的谱。式 (14) 用权重 a1=0.56a_1=0.56 组合两个吸收谱。这里的 0.56 是这张图的分解系数,不能搬给另一份混合物当普适参数。

重要变化是:横轴上不同位置混入两个组分的比例不同。某个切片的最大值还会随两个组分的寿命改变。混合中心线因此一般不再是两个相关函数的简单加权平均;它还要纳入随频率与时间变化的组分比例。

对应阅读:Fenn and Fayer

寿命差异会让混合权重随时间漂移

式 (15) 进一步写出频率与等待时间相关的权重。用 S1(x),S2(x)S_1(x),S_2(x) 表示某一激发频率处两个组分的初始谱贡献,以 L1,L2L_1,L_2 表示它们的布居寿命。这里把寿命改用字母 LL,避免和第一段相干时间混淆。

做一个新的教学算例:固定某个频率,两组分初始贡献分别为 0.56 和 0.44,寿命分别为 1 ps1\,\mathrm{ps} 与 2 ps2\,\mathrm{ps}。等待 1 ps1\,\mathrm{ps} 后,第一组分在剩余信号中的比例约为 0.436。即便组分各自的环境完全不变,混合信号的权重仍在变。

f1(x,T)=S1(x)e−T/L1S1(x)e−T/L1+S2(x)e−T/L2f_1(x,T)=\frac{S_1(x)e^{-T/L_1}}{S_1(x)e^{-T/L_1}+S_2(x)e^{-T/L_2}}
按论文式 (15) 改写符号。上面的寿命为自定算例值,未声称是论文所测寿命。

对应阅读:Fenn and Fayer

把限定复现与分子推断分开

本路线的复现对象是式 (12) 的解析 Gaussian 极限与其独立数值检验。交互实验另外生成有有限带宽、有限时间窗和噪声的复信号,再比较两类相关模型的留出预测。它没有数字化论文中的水谱图,也没有重新估计反胶束中两组分的频率相关。

连续液体涨落在相干时间内也会改变频率。当前频率对模型把这段变化近似冻结;当相关时间接近或短于相干观测时间时,它只是一种有明确定义的现象模型。要做真实溶剂的定量推断,应恢复完整随机频率过程与响应卷积。

对应阅读:Tokmakoff

下一篇论文可以沿这条顺序读

先记下版本、图号、轴与单位,确认每张图是测量、计算还是理想参考。随后找到从原始信号到图的处理步骤,检查相位、窗口与切片方向。最后再列出哪些参数由哪些观测决定,是否有独立数据检查预测。

如果两种微观模型在现有带宽与噪声下给出同样的谱,下一步应设计能把它们分开的观测。换偏振、拉长等待窗口或加入独立结构信息,比在同一张图上继续增加参数更有说服力。

动手检验预测

对照式 (12) 与有噪声的推断

  1. 先以低噪声、宽带宽查看中心线随等待时间的变化。
  2. 读出拟合模型的留出误差,比较哪些晚时记录更能分辨静态部分与慢衰减。
  3. 增大噪声后,解释偶尔越出物理范围的估计为什么应作为诊断保留,不能直接裁回范围再声称拟合良好。

练习与解释性答案

  1. 混合样品的中心线变平,能否直接解释成每个组分都失去了同样多的频率记忆?

    展开解释

    不能。组分重叠、频率相关的混合比例与不同布居寿命都能改变中心线。需要组分吸收、寿命和适合混合体系的响应分析。

  2. 复算式 (12) 为什么还不能称作复现论文的全部结果?

    展开解释

    它只验证了作者分析中的一个明确极限。论文还研究两组分条件、具体谱线和反算算法,需要额外参数、实验数据与实现。本路线逐项标明采用与未采用的部分。

来源与计算范围

  1. Tokmakoff, Two-dimensional spectroscopy chapter (June 2011)

    同一份存档讲义中,页码重新从 1 开始的二维光谱部分:页 4—5、式 (1)—(3) 连接响应、Fourier 变换与异频探测;页 10—11 讨论耦合振动的成对峰。不转载原图。

  2. Fenn and Fayer, Extracting 2D IR frequency-frequency correlation functions from two component systems (2011)

    Journal of Chemical Physics 135, 074502,出版版本,DOI 10.1063/1.3625278。核读图 1—3 与式 (1)—(4)、(12)、(14)—(15)。限定数值基准采用式 (12) 的单组分 Gaussian 斜率关系,参数为自定教学值,不拟合论文中的水实验。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策