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

L3 · 第 6 章

怎样把保存的复信号变成二维谱

二维谱的每个像素都可以追溯到一组带相位的时间信号。先用本振测出实部和虚部,再沿两个相干时间做 Fourier 变换,最后把重聚与不重聚贡献按同一相位约定相加。采样窗、步长和脉冲频谱会决定这张图能看清什么。

探测器怎样量到一个场

探测器直接响应总光强。把弱信号 EsE_s 与强本振 ELE_L 叠加,扣去本振单独的强度后,留下主要交叉项。它对弱信号振幅是一阶,因此弱场也容易测量。

一个本振相位只给出信号的一个投影。将参考相位再转四分之一圈,可以读出另一个正交投影。保存两者,才组成复时间信号。单独保存光强的绝对值无法还原任意信号相位。

∣EL+Es∣2−∣EL∣2=2Re⁡(EL∗Es)+∣Es∣2|E_L+E_s|^2-|E_L|^2=2\operatorname{Re}(E_L^*E_s)+|E_s|^2
强本振、弱信号近似下,最后一项通常比交叉项小。实验仍需处理背景、相位漂移与探测器响应。

对应阅读:Tokmakoff · Park

先固定两根轴的正方向

横轴取激发波数,纵轴取探测波数;两根频率轴都向正值增加。以 RR∝e+iωat1−iωbt3R_{\mathrm R}\propto e^{+i\omega_at_1-i\omega_bt_3} 为重聚约定,以 RNR∝e−iωat1−iωbt3R_{\mathrm{NR}}\propto e^{-i\omega_at_1-i\omega_bt_3} 为不重聚约定。要让同一跃迁在两张图中落到同样的正频率,第一根变换轴必须采用相反符号。

SR(ω1,ω3)=∫0∞ ⁣∫0∞RRe−iω1t1+iω3t3 dt1dt3S_{\mathrm R}(\omega_1,\omega_3)=\int_0^\infty\!\int_0^\infty R_{\mathrm R}e^{-i\omega_1t_1+i\omega_3t_3}\,dt_1dt_3
重聚变换。实际记录有有限窗口,积分替换为带已知权重的求和。
SNR(ω1,ω3)=∫0∞ ⁣∫0∞RNRe+iω1t1+iω3t3 dt1dt3S_{\mathrm{NR}}(\omega_1,\omega_3)=\int_0^\infty\!\int_0^\infty R_{\mathrm{NR}}e^{+i\omega_1t_1+i\omega_3t_3}\,dt_1dt_3
不重聚变换。在脉冲顺序、公共相位和归一化一致后,取两者之和的实部得到所选吸收型谱。

对应阅读:Tokmakoff

完整算例:一步走多远,总共看多久

假设旋转参考系中的频率偏移每隔 Δt=0.05 ps\Delta t=0.05\,\mathrm{ps} 采样一次。可区分的偏移范围受 Nyquist 条件限制,半范围约为 333.6 cm−1333.6\,\mathrm{cm}^{-1}。若观测窗只有 3.2 ps3.2\,\mathrm{ps},频率网格特征尺度约为 10.4 cm−110.4\,\mathrm{cm}^{-1};真正分辨本领还受窗函数与衰减影响。

这里采的是相对载频的慢变化复包络。若直接采原始中红外光场,同样的时间步长可能连光学载波都采不准。旋转参考系必须写进记录,才能知道频谱坐标该加回哪个载频。

∣δν~∣<12cΔt,δν~grid≈1cTwindow|\delta\widetilde\nu|<\frac{1}{2c\Delta t},\qquad \delta\widetilde\nu_{\mathrm{grid}}\approx\frac{1}{cT_{\mathrm{window}}}
第一式控制混叠,第二式给出有限窗的频率网格尺度。增加补零点数会让图更平滑,不延长实际观测窗。

有限脉冲与窗函数都会改变峰

脉冲频谱像一个已知的频率选择器:带宽不够时,偏离载频的跃迁较弱。频谱边缘的弱峰不能直接与中央峰按原始峰高比较。有限采样窗又会切断仍未衰减的信号,产生旁瓣。

给记录乘平滑的窗函数可以压低旁瓣,代价是峰变宽。有噪声时,加长记录也未必总有帮助:后面可能主要在累计噪声。要比较两组等待时间,须用同一套窗与定相方法,并保留处理参数。

对应阅读:Park

保存哪些量才能重算

一份能复算的记录应包含:两个时间坐标、两类路径各自的实部与虚部、等待时间、载频、单位、脉冲权重、窗口、随机种子和模型版本。重建图要从这份记录计算;改显示色标不应改变原始数据。

还可以把时间域直接求和与另一份独立积分比较,检查峰的位置、符号与幅度。若改了采样步长,峰突然移到了别处,先检查混叠与 Fourier 符号,再讨论分子是否改变了。

动手检验预测

从导出记录重看正负峰

  1. 保存一份复信号与谱图,记下采样窗和载频约定。
  2. 保持固定采样网格,改变相干寿命或噪声,预测哪些峰最先变模糊。
  3. 把峰高变化与峰位置变化分开解释。

比较相同处理下的等待时间系列

  1. 保持带宽、窗口和噪声设置不变,浏览不同等待时间。
  2. 导出整个系列,检查各时刻是否共享相同的坐标和处理参数。

练习与解释性答案

  1. 补零让频率网格密了四倍,能把原来重叠的两条窄线分开吗?

    展开解释

    补零只插值已有有限窗的 Fourier 结果,不增加独立的时间信息。可靠地区分更近的频率通常需要更长的有效相干记录,也要有足够信噪比。

  2. 若只对重聚和不重聚变换使用同样的第一轴符号,会发生什么?

    展开解释

    同一正激发频率在其中一张图上会反号。直接相加可能把对应峰放到不同位置,破坏原本应当组合成吸收型线形的贡献。

来源与计算范围

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

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

  2. Park, Kwak and Fayer, Ultrafast 2D-IR vibrational echo spectroscopy (2007)

    Laser Physics Letters 4, 704—718,DOI 10.1002/lapl.200710046。实验综述,第 2—3 节及图 3 用于脉冲、本振和吸收谱定相。文中的另一种切片方向使用逆斜率;课程先固定坐标,再定义斜率。

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

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