二次谐波产生原理：倍频晶体与相位匹配

拆开一支普通的绿色激光笔，里面通常并不是一颗发出绿光的芯片。它先用激光二极管泵浦一块掺钕晶体，产生 1064 nm 的红外激光；这束红外光再穿过一小片非线性晶体，出来时波长变成了 532 nm——恰好是一半，于是我们看到了绿色。

红外光的波长是绿光的两倍，所以这件事在能量上是说得通的：两个 1064 nm 的光子，正好凑成一个 532 nm 光子的能量。但奇怪的地方在于，光子并不会像两块橡皮泥那样被”捏”在一起。光穿过普通玻璃、水或空气时，只会按折射、反射的规则改变方向，波长分毫不动。那么，是什么材料、什么条件，让介质把两个基频光子换成了一个倍频光子？

这就是二次谐波产生（second harmonic generation，SHG，中文常直接称为”倍频”）。它是非线性光学里最基础、也最容易在实验室里重复的现象之一。本文要回答三个问题：介质的什么性质允许它发生？为什么效率高低几乎由”相位匹配”一个条件决定？工程上又是怎么把这个条件做出来的？

非线性极化：强光下材料的响应不再只是线性

光进入介质时，电场会推动材料中的电子偏离平衡位置，形成电偶极矩。把单位体积内的偶极矩加起来，得到极化强度 $P$。在光强不高、介质响应温和的情况下，极化强度与电场强度 $E$ 成正比：

$$P=\varepsilon_0\chi^{(1)}E$$

$\chi^{(1)}$ 是一阶（线性）极化率，它决定了我们熟悉的折射率和吸收。这一步成立时，介质对光的响应是线性的：频率为 $\omega$ 的入射光进入，出来的还是频率为 $\omega$ 的光。

当电场足够强时，电子被推离平衡位置的距离已经不能算”微小”，上面的正比关系不再够用，极化强度需要按电场展开成级数：

$$P=\varepsilon_0\left(\chi^{(1)}E+\chi^{(2)}E^2+\chi^{(3)}E^3+\cdots\right)$$

其中 $\chi^{(2)}$ 是二阶非线性极化率，$\chi^{(3)}$ 是三阶非线性极化率。激光之所以让非线性光学成为可能，正是因为它能把电场强度推到足够高；反过来，普通光源的强度下，高次项小到可以忽略。

不过，光有强光还不够。二阶项 $\chi^{(2)}E^2$ 只在缺少反演对称中心的介质里才非零。所谓反演对称，是指把晶体中的每个原子都沿中心镜像翻转后，结构看起来完全一样。各向同性的气体、液体和玻璃通常满足这种对称性，二阶非线性在电偶极近似下被对称性直接”关掉”；因此倍频必须在特定的非线性晶体或人工结构里完成，常见的材料是铌酸锂（LiNbO₃）、磷酸钛氧钾（KTP）、三硼酸锂（LBO）和偏硼酸钡（BBO）。

一个容易被忽略的细节是：这些晶体的非线性系数并不是各向同性的一个数，而是张量，不同方向和不同偏振组合对应不同分量。后面讲准相位匹配时我们会看到，选对张量分量可以显著改变倍频效率。

从 $E^2$ 到 $2\omega$：二次谐波项从哪里来

倍频的数学来源其实只用到一步三角恒等式。设入射光是频率为 $\omega$ 的单色波 $E(t)=E_0\cos\omega t$，把它代入二阶极化项：

$$\chi^{(2)}E^2=\chi^{(2)}E_0^2\cos^2\omega t=\frac{\chi^{(2)}E_0^2}{2}\left(1+\cos 2\omega t\right)$$

结果里出现了两个分量：一个不随时间振荡的直流项，以及一个频率为 $2\omega$ 的分量。也就是说，”两个光子合成一个光子”并不需要光子之间发生碰撞——介质的非线性响应本身就把 $2\omega$ 的频率成分写进了极化强度，这个分量再辐射出去，就是我们看到的倍频光。

从量子图像看，这个过程可以画成两个入射光子通过一个虚能级同时被湮灭，同时产生一个频率为 $2\omega$ 的光子。虚能级并不是真实的能量本征态，中间态存在的时间极短，所以二次谐波产生不需要满足共振条件，但需要足够高的光强来提供足够的”同时性”。这也解释了为什么同样的晶体在连续弱光下几乎看不到倍频，而在调 Q 或锁模脉冲激光的峰值功率下就相当明显。

用能量守恒可以做一个定量检查。由 $hc=1239.84\ \mathrm{eV\cdot nm}$，1064 nm 光子能量为 $1239.84/1064\approx1.165\ \mathrm{eV}$，532 nm 光子能量为 $1239.84/532\approx2.331\ \mathrm{eV}$，恰好是前者的两倍：$2\times1.165=2.330$，与 2.331 只差千分之一，差值来自 1064 与 532 nm 的取值精度。倍频把波长减半、光子能量翻倍，能量账目是平的。

图 1：单个偶极子（上排）、两个平行偶极子（中排）与两个反平行偶极子（下排）产生的二次谐波信号对比。反平行排列时，各偶极子贡献的二次谐波相互抵消；这也对应了反演对称介质中 $\chi^{(2)}=0$ 的结论。

图 1 给出了”对称性决定信号强弱”的直观图像：二次谐波是从材料中被光驱动的大量偶极子辐射出来的，只有当这些贡献在观察方向上同相叠加时，才能得到可观的信号。反平行排列让相邻偶极子的辐射相消，整体信号被抵消——这正是反演对称介质不做倍频的微观原因。这也提醒我们：倍频效率不只取决于光强和材料，还取决于材料内部结构的排布方式。

相位匹配：两个频率必须以同样的速度传播

倍频光不是在某一个截面上突然产生的，而是在光线穿过晶体的整条路径上一点一点产生。晶体里每个位置都受到基频光驱动，向外辐射一份倍频光；观察者最终在出射面接收到的总幅度，是这些来自不同深度的贡献之和。

问题就出在求和上。基频光在介质中以波矢 $k(\omega)=n(\omega)\omega/c$ 传播，频率翻倍之后，新产生的那份倍频光对应的”驱动相位”每单位长度增加 $2k(\omega)$；而已有的倍频光波自身以 $k(2\omega)$ 传播。两者每走一段距离，就会积累相位差：

$$\Delta k=k(2\omega)-2k(\omega)=\frac{2\omega}{c}\left[n(2\omega)-n(\omega)\right]$$

只有当 $\Delta k=0$，也就是

$$n(2\omega)=n(\omega)$$

时，不同位置产生的贡献才会一直保持同相，越加越多。这个条件就叫相位匹配。

麻烦的是，材料的折射率通常随频率升高而增大（正常色散），所以 $n(2\omega)$ 往往略大于 $n(\omega)$，相位匹配并不是天然成立的。一旦 $\Delta k\neq0$，刚产生的倍频光走出一段距离后相位就反转，开始把能量”还”回基频光，功率在两条波之间来回振荡；每经过一个相干长度

$$L_c=\frac{\pi}{\Delta k}=\frac{\lambda}{4\left[n(2\omega)-n(\omega)\right]}$$

方向就翻转一次。结果是：晶体再长，倍频输出也不会持续增长。换句话说，晶体长度的意义不在于”越长吸收的光越多”，而在于能否在整个长度上保持同相叠加。

图 2：器件结构与工作原理。(a) 集成非线性光子器件示意，用梯度超表面（天线阵列）在 LiNbO₃ 波导中实现免相位匹配的二次谐波产生；(b) 器件截面示意；(c) 制备器件的扫描电子显微镜（SEM）照片，左图为不同长度的硅纳米棒天线阵列，右图为局部放大，比例尺均为 1 μm；(d) 免相位匹配的工作原理：泵浦先耦合到倍频的基波导模式，再由梯度超表面把倍频光耦合到高阶波导模式，形成单向功率转移。

图 2 把相位失配的后果画得很直接：在未经结构处理的波导里，倍频信号沿传播方向来回振荡，相干累积长度只有约 2 μm 量级；加上梯度超表面提供单向有效波矢后，倍频功率随传播距离单调上升，并且能在 1580–1650 nm 的泵浦范围内保持高效。这也说明相位匹配并非只能靠材料本身解决，还可以用人工结构去”改写”动量条件。

数值代入：KDP 的相位匹配角与相干长度

既然折射率决定了能不能匹配，一个自然的思路就是人为制造两个频率折射率相等的条件。单轴晶体正好提供了这种可能：同一块晶体里，光的偏振方向不同，感受到的折射率不同——寻常光感受 $n_o$，非常光感受 $n_e$，而且非常光的折射率随传播方向与光轴的夹角 $\theta$ 变化：

$$\frac{1}{n_e^2(\theta)}=\frac{\cos^2\theta}{n_o^2}+\frac{\sin^2\theta}{n_e^2}$$

让基频光走寻常光、倍频光走非常光，再转动晶体改变 $\theta$，就可能让两者的折射率相等。教材给出的 KDP 晶体数据是 $n_o(\omega)=1.505$、$n_e(\omega)=1.465$、$n_o(2\omega)=1.534$、$n_e(2\omega)=1.487$。相位匹配要求 $n_e(2\omega,\theta)=n_o(\omega)=1.505$，代入上式：

$$\frac{1}{1.505^2}=\frac{\cos^2\theta}{1.534^2}+\frac{\sin^2\theta}{1.487^2}$$

三个平方倒数分别是 $1/1.505^2=0.4415$、$1/1.534^2=0.4250$、$1/1.487^2=0.4523$，整理后

$$\sin^2\theta=\frac{0.4415-0.4250}{0.4523-0.4250}\approx0.606$$

得到 $\theta\approx51.1^\circ$。也就是说，把 KDP 晶体的光轴转到与被倍频红外光成约 51° 的位置，两个偏振方向上的折射率就大致对齐了。实际器件还会用温度微调，把角度固定在接近 90° 的位置（非临界相位匹配），以减小走离效应。

反过来也可以估算”不做相位匹配会怎样”。这组数据里 $n(2\omega)-n(\omega)=1.534-1.505=0.029$，代入相干长度公式：红宝石激光（694.3 nm）对应的相干长度

$$L_c=\frac{694.3\ \mathrm{nm}}{4\times0.029}\approx6.0\ \mathrm{\mu m}$$

若把波长换成 1064 nm，则 $L_c\approx9.2\ \mathrm{\mu m}$。两者都只有微米量级，而晶体往往是毫米到厘米长——这意味着没有相位匹配时，绝大部分晶体长度都在做”无用功”，甚至把已经产生的倍频光又转换回去。这也解释了为什么倍频效率对相位匹配如此敏感，以及为什么器件上常常要精确控制晶体的角度和温度。

准相位匹配：每走过一个相干长度就把符号翻过来

双折射相位匹配虽然经典，但限制不少：需要晶体双折射足够大，可能出现空间走离，能够匹配的波长范围也偏窄。

准相位匹配（quasi-phase matching，QPM）换了一个思路：既然相位失配无法消除，那就让它”每一段都重新开始”。在每经过一个相干长度 $L_c$ 之后，把介质的非线性系数符号反转一次，本来会相消的那部分贡献就变成相加，倍频幅度于是可以继续增长。实现方式通常是周期性极化——用外加电场周期性反转铁电晶体（如铌酸锂）的自发极化方向，让非线性系数沿传播方向交替取正负，极化周期约为 $2L_c$，常见器件的周期在数十微米量级。

准相位匹配的代价是有效非线性系数要乘上 $2/\pi$（效率按系数平方下降）。但它换来两个好处：一是三个波可以取相同的偏振方向，从而使用张量中最大的分量——以 LiNbO₃ 为例，双折射相位匹配常用 $d_{31}=4.35\ \mathrm{pm/V}$，而准相位匹配使用 $d_{33}=27\ \mathrm{pm/V}$，折算 $2/\pi$ 之后仍相当于约 17 pm/V；二是传播方向可沿晶轴，避免空间走离，接受角也更大。对泵浦功率有限的场景，这笔交换通常很划算，这也是周期极化铌酸锂（PPLN）成为绿光、蓝光倍频和光学参量振荡器主流器件的原因。

表 1：三条实现二阶非线性频率转换的技术路线对比。

路线	核心做法	关键条件	主要优势	适用边界
双折射相位匹配	让基频光与倍频光取不同偏振，转动晶体角度或调节温度	晶体双折射足够大，满足 $n(2\omega)=n(\omega)$	结构简单，块状晶体即可实现，不需要额外工艺	可用波长受双折射大小限制；可能出现空间走离；相位匹配带宽窄
准相位匹配（PPLN 等）	周期性反转非线性系数符号，周期约为 $2L_c$	极化周期与波长、温度、色散匹配	可用最大张量分量 $d_{33}$，允许同偏振相互作用，无空间走离	需要周期性极化工艺；周期对温度敏感；带宽仍然有限
梯度超表面等结构方案	用人工结构提供额外单向有效波矢，抑制逆过程	结构波矢与波导模式匹配，制造精度足够	绕开传统相位匹配条件，可在较宽泵浦波长范围内工作	器件加工复杂，效率、损耗与功率承受能力仍在优化

从激光笔到显微成像：倍频在工程中的用法

倍频最广为人知的用途就是把红外激光搬到可见光。最常见的 532 nm 绿光由 1064 nm 的 Nd:YAG 或 Nd:YVO₄ 激光倍频得到，倍频晶体可以放在谐振腔外，也可以直接放进腔内利用高腔内功率；小功率场合则使用倍频微片激光器或倍频激光二极管——绿色激光笔就是这条路线。继续往短波走，通过倍频与和频可以进入蓝光与紫外波段，此时材料透明范围和损伤阈值成为主要约束。

另一条路是把倍频当作成像手段。二次谐波显微成像的信号来自样品内部具有非中心对称性的结构排列，因此特别适合观察胶原蛋白等生物组织，并发展出偏振分辨、干涉型等多种变体。典型装置由激光源、功率控制单元、扫描系统、物镜和探测器组成，探测器通常采用光电倍增管（PMT），并与计算机同步完成扫描采集。

图 3：二次谐波显微成像的典型装置，包含激光源、功率控制单元、扫描系统、物镜与探测器；探测器与计算机相连，用于控制显微镜并把激光扫描与信号采集同步起来。

在集成光学里，倍频还有一种低功率的实现方式：把光限制在很小的模式面积内，让非线性相互作用在厘米级长度上持续进行。非线性波导在铌酸锂、钽酸锂等材料上已有成熟工艺，厘米长器件的倍频效率可以超过 100%/W，配合周期极化结构后，甚至能对连续光实现可观的转换效率。

适用边界与待推进方向

二次谐波产生的门槛写在它的名字里：必须找到缺少反演对称中心的介质，并且把光强推到足够高。低功率下倍频强度随泵浦强度近似平方增长，因此用连续弱光单次穿过块状晶体，效率通常很低；工程上要么使用调 Q、锁模脉冲激光提高峰值功率，要么把晶体放进谐振腔，或改用波导结构。相位匹配本身也带来边界：相位匹配带宽有限，空间走离与群速度失配会限制超短脉冲的倍频，往紫外走还要同时面对材料透明范围和损伤阈值。

正在推进的方向大致有三条。其一是把非线性作用限制在波导或微腔中，用更长的有效作用长度换取更低的泵浦功率；其二是通过周期极化、啁啾结构等设计，让相位匹配条件在更宽的波长范围或更高的温度稳定性下成立；其三是用梯度超表面等人工结构直接改写动量条件，绕开传统相位匹配，并已在薄膜铌酸锂波导中实现宽带倍频。对大多数实验来说，理解 $n(2\omega)=n(\omega)$ 这一个条件，就已经能解释倍频器件为什么需要精确的角度与温度控制，也能解释为什么效率可以从几乎为零跳到 50% 以上。

参考资料

1. RP Photonics Encyclopedia, “Frequency Doubling”. https://www.rp-photonics.com/frequency_doubling.html
2. RP Photonics Encyclopedia, “Quasi-phase Matching”. https://www.rp-photonics.com/quasi_phase_matching.html
3. RP Photonics Encyclopedia, “Periodic Poling”. https://www.rp-photonics.com/periodic_poling.html
4. RP Photonics Encyclopedia, “Green Lasers”. https://www.rp-photonics.com/green_lasers.html
5. RP Photonics Encyclopedia, “Nonlinear Crystal Materials”. https://www.rp-photonics.com/nonlinear_crystal_materials.html
6. J. Peatross & M. Ware, Physics of Light and Optics（BYU，作者免费分发），第 5–8 章复习题 R40：KDP 晶体双折射相位匹配算例.
7. A. Aghigh, S. Bancelin, M. Rivard, M. Pinsard, H. Ibrahim, F. Légaré, “Second harmonic generation microscopy: a powerful tool for bio-imaging”, Biophysical Reviews, 2023. https://doi.org/10.1007/s12551-022-01041-6
8. C. Wang, Z. Li, M.-H. Kim, X. Xiong, X.-F. Ren, G.-C. Guo, N. Yu, M. Lončar, “Metasurface-assisted phase-matching-free second harmonic generation in lithium niobate waveguides”, Nature Communications 8, 2098 (2017). https://doi.org/10.1038/s41467-017-02189-6