4、强度自相关法（一）：从干涉仪到脉冲宽度估算

超短脉冲（飞秒至皮秒量级）的时间宽度远小于任何电子探测器的响应时间（通常为纳秒量级），因此无法直接用示波器测量。强度自相关法（Intensity Autocorrelation）是解决这一难题的最经典、最实用的手段之一。它通过将脉冲与其自身在时间上延迟的副本在非线性晶体中相互作用，产生一个与脉冲宽度相关的信号，从而间接“读出”脉冲的时间尺度。本节将系统讲解强度自相关的基本原理、迈克尔逊干涉仪结构、二阶自相关信号的产生过程，以及如何从自相关曲线估算脉冲宽度。

上图展示了强度自相关测量的完整信号链路：待测脉冲经分束镜分为两束，一束经固定镜反射，另一束经移动镜反射产生时间延迟 τ，两束光在非线性晶体（如BBO）中重合产生倍频信号，该信号的强度随延迟 τ 变化，最终由探测器记录并绘制出自相关曲线。

4.1 核心原理：为什么需要“自相关”？

超短脉冲的电场表达式可写为：

E(t) = A(t) · exp(iω₀t + iφ(t))

其中 A(t) 是脉冲包络（振幅），ω₀ 是中心角频率，φ(t) 是相位项。强度自相关法利用脉冲与自身延迟副本的非线性相互作用，将时间信息转换为空间位移信息（移动镜位置），再转换为可测量的光强信号。其数学基础是二阶强度自相关函数：

G₂(τ) = ∫ I(t) · I(t-τ) dt

式中 I(t) = |E(t)|² 为脉冲强度包络，τ 为两束光之间的时间延迟。G₂(τ) 描述了脉冲与自身延迟副本的重叠程度：当 τ=0 时，两脉冲完全重叠，信号最大；当 |τ| 大于脉冲宽度时，重叠消失，信号趋于零（对于单脉冲）或趋于背景值（对于脉冲序列）。

关键在于，G₂(τ) 的半高全宽（FWHM）与脉冲本身的宽度存在确定的比例关系。对于高斯脉冲，自相关宽度 Δτ_ac 与脉冲宽度 Δτ_p 的关系为：

Δτ_p = Δτ_ac / √2 ≈ 0.707 × Δτ_ac

对于双曲正割脉冲（sech²），比例系数为 0.648。因此，只要测量出自相关曲线的宽度，乘以相应的去卷积因子，即可得到脉冲宽度。这一原理避免了直接测量飞秒时间尺度，将问题转化为测量微米级的位移（1飞秒光程差约0.3微米），大大降低了测量难度。

4.2 迈克尔逊干涉仪结构：搭建与光路对准

强度自相关仪的核心是迈克尔逊干涉仪结构，其基本组成包括：

• 分束镜（BS）：将入射脉冲分为等强度的两束，通常使用50:50分束比，带宽需覆盖脉冲光谱。
• 固定反射镜（M1）：反射其中一束光，构成参考臂。
• 移动反射镜（M2）：安装在精密位移台（如步进电机或压电陶瓷）上，构成延迟臂，通过改变光程产生时间延迟 τ。
• 非线性晶体：常用BBO（β-硼酸钡）或KDP晶体，实现倍频（二次谐波产生，SHG）。
• 聚焦/准直透镜：保证两束光在晶体中空间重合，且满足相位匹配条件。
• 滤光片与探测器：滤除基频光，只让倍频信号进入光电倍增管（PMT）或硅光电二极管。

光路对准是实验成功的关键。具体操作流程如下：

1. 共路调节：先用连续激光（如He-Ne激光器）代替飞秒脉冲，调节M1和M2，使两束反射光在远场（如3米外的光屏）完全重合，确保两臂平行。
2. 晶体位置优化：将BBO晶体置于两束光的交点处，通过三维调节架微调晶体角度，使倍频信号最强。对于I类相位匹配，两束光偏振方向相同；对于II类匹配，偏振方向正交。
3. 延迟扫描：驱动位移台以恒定速度扫描M2，同时记录探测器输出。扫描范围应覆盖脉冲宽度的5-10倍，以确保完整记录自相关曲线。
4. 信号优化：观察示波器上的自相关信号，微调晶体角度和聚焦透镜位置，使信号幅度最大且波形对称。

4.3 二阶自相关信号产生过程：从光程差到倍频信号

二阶自相关信号的产生涉及非线性光学中的和频过程。当两束频率相同的脉冲在非线性晶体中时空重合时，会产生频率为 2ω₀ 的倍频信号。其强度正比于两束入射光强度的乘积：

I_SHG(τ) ∝ I₁(t) · I₂(t-τ)

探测器（响应时间远大于脉冲宽度）记录的是时间积分信号，即：

S(τ) = ∫ I₁(t) · I₂(t-τ) dt = G₂(τ)

这一过程的关键物理图像是：只有当两脉冲在时间上重叠时，晶体中才会产生显著的倍频信号。通过移动M2改变延迟 τ，相当于“扫描”两脉冲的重叠程度，从而获得完整的自相关曲线。

值得注意的是，强度自相关（非干涉自相关）要求两束光的偏振方向正交（II类相位匹配），以避免基频光与倍频光之间的干涉条纹。如果两束光偏振平行（I类匹配），探测器将记录到带有干涉条纹的信号，这属于干涉自相关（Interferometric Autocorrelation），可用于提取脉冲的相位信息，但分析更为复杂。本节聚焦于强度自相关，因此采用正交偏振配置。

4.4 从自相关曲线估算脉冲宽度：数据处理与去卷积

获得自相关曲线后，脉冲宽度的估算分为以下步骤：

1. 数据预处理：扣除背景噪声（暗电流和散射光），对曲线进行平滑处理（如移动平均或Savitzky-Golay滤波）。
2. 确定半高全宽：找到曲线的最大值 I_max，计算 I_max/2 对应的两个横坐标位置，其差值即为自相关宽度 Δτ_ac。
3. 去卷积因子选择：根据脉冲形状选择相应的去卷积因子。常用脉冲形状及其因子如下：

# 示例：假设测得 Δτ_ac = 150 fs，脉冲形状为 sech²
Δτ_p = 150 fs / 1.543 ≈ 97.2 fs

需要注意的是，去卷积因子的选择依赖于对脉冲形状的先验假设。对于飞秒激光器，通常假设为 sech² 形状（锁模激光器典型输出）；对于高斯脉冲（如经滤波后的脉冲），则使用 √2 因子。如果无法确定脉冲形状，可以同时给出两种假设下的估算值，并注明假设条件。

4.5 案例拆解：测量钛宝石激光器输出的脉冲宽度

背景与目标：某实验室拥有一台钛宝石锁模激光器，标称中心波长 800 nm，重复频率 80 MHz，预期脉冲宽度约 100 fs。需要精确测量其脉冲宽度以验证激光器性能。

实验条件：

• 待测脉冲平均功率：500 mW，单脉冲能量约 6.25 nJ。
• BBO晶体：厚度 50 μm，I类相位匹配，切割角 29.2°。
• 位移台：步进电机驱动，最小步长 0.5 μm（对应 3.3 fs 时间延迟）。
• 探测器：硅光电二极管，响应时间 1 ns。

操作过程：

1. 将激光器输出导入自相关仪，调节分束镜使两束光功率相等（各约 250 mW）。
2. 用He-Ne激光器进行共路调节，确保两束光在3米外完全重合。
3. 将BBO晶体置于焦点处，微调角度使倍频信号最大。观察到绿色倍频光（400 nm）后，用滤光片滤除基频光。
4. 设置位移台扫描范围 ±2 ps，步长 5 fs，扫描速度 10 fs/s。记录探测器输出信号。
5. 重复扫描3次，取平均以降低噪声。

结果与分析：测得自相关曲线半高全宽 Δτ_ac = 155 fs。假设脉冲为 sech² 形状，则脉冲宽度 Δτ_p = 155 / 1.543 ≈ 100.5 fs。与标称值 100 fs 吻合良好。若假设为高斯形状，则 Δτ_p = 155 / 1.414 ≈ 109.6 fs，偏差约 9%。这表明脉冲形状更接近 sech²。

风险与注意事项：实验中观察到信号基线不为零（约为峰值的5%），这可能是由于晶体中非相位匹配的倍频过程或散射光所致。处理时需扣除基线，否则会导致自相关宽度偏大。此外，位移台步长过大（>10 fs）会降低曲线分辨率，导致宽度测量偏大。

4.6 总结与复盘

强度自相关法作为超短脉冲测量的入门技术，具有结构简单、成本较低、操作直观等优点。其核心要点可归纳为：

• 原理本质：利用非线性倍频过程，将时间延迟转换为可测量的光强变化。
• 结构关键：迈克尔逊干涉仪的两臂光程差控制精度直接影响测量分辨率。
• 信号解读：自相关曲线的宽度需乘以去卷积因子才能得到真实脉冲宽度，因子选择依赖脉冲形状。
• 局限认知：强度自相关不提供相位信息，无法判断脉冲是否带有啁啾；对于复杂脉冲形状，测量结果可能产生误导。

在实际应用中，建议先使用强度自相关快速评估脉冲宽度，若需进一步优化脉冲（如压缩或整形），再引入FROG等更高级的测量手段。掌握本节内容后，学习者应能独立搭建强度自相关仪，并准确估算飞秒脉冲的宽度。