文献解读 | 超表面慢光纳米涂层:飞秒激光脉冲的色散控制与压缩
一句话结论:610 nm 厚的硅纳米柱涂层用慢光给出 −71 fs² 异常色散,把 48.3 fs 脉冲压到 37.5 fs。
论文简介
这篇论文是 M. Ossiander、Y.-W. Huang、W. T. Chen、Z. Wang、X. Yin、Y. A. Ibrahim、M. Schultze 与 F. Capasso 发表在 Nature Communications 第 12 卷、文章号 6518(2021 年 11 月 11 日)上的 Slow light nanocoatings for ultrashort pulse compression,DOI 为 10.1038/s41467-021-26920-6。作者主要来自哈佛大学 John A. Paulson 工程与应用科学学院,合作单位包括台湾阳明交通大学、格拉茨技术大学与滑铁卢大学。
论文演示了一片 610 nm 厚的非晶硅纳米柱涂层:它在 800 nm 附近的全工作范围内给出 −71 fs² 的异常群延迟色散,覆盖约 80 nm 带宽,把 48.3 fs 的钛宝石脉冲压到 37.5 fs。整片元件面积 3.1 mm²,直接插在光路里,不需要光栅或透镜先把光谱在空间上分开。它的贡献落在色散控制这一环:把通常由棱镜对、光栅对或啁啾镜承担的补偿,压缩成一片可镀在现有光学元件上的超薄涂层。
这篇论文发表于 2021 年,按“优先近 3 年”的常规口径偏旧。把它作为主读,是因为它较早把透射式慢光色散从理论方案做成可实测元件,并给出单层结构的色散上限;色散强度、带宽、加工容差和时域验证在这里形成完整闭环,同方向后来向多层结构扩展的工作目前仍以预印本为主。对激光链路来说,它代表的是增益介质之外、脉冲在到达实验之前必须经过的色散控制环节。
一片涂层要补多少色散
先看一个具体的量级。论文的数值算例设定是:一个 10 fs 的高斯脉冲穿过 1.6 mm 熔石英后,会被展宽到 19 fs。原因不是玻璃吸收光,而是不同频率的光在玻璃里速度不同。可见光和近红外波段,蓝光成分的群延迟通常大于红光成分,脉冲前后沿的相对位置被拉开,时域包络变长。
把这件事写成公式,光谱相位记为 $\varphi(\omega)$,群延迟和群延迟色散分别是
$$ \mathrm{GD}=\frac{d\varphi}{d\omega},\qquad \mathrm{GDD}=\frac{d\mathrm{GD}}{d\omega}=\frac{d^2\varphi}{d\omega^2}. $$
GDD 的单位是 fs²。常规透明材料在 800 nm 附近是正常色散,GDD 为正;要让已经展宽的脉冲重新变短,涂层需要提供符号相反、量级相当的负 GDD。对一个带宽为 $\Delta\omega$ 的脉冲,群延迟在带宽两端的变化量近似为
$$ \Delta\tau \approx \mathrm{GDD}\cdot \Delta\omega. $$
熔石英在 800 nm 的 GDD 约为 +36.2 fs²/mm。2 mm 厚的窗口给出 +72.4 fs²。若带宽取 80 nm,中心波长 800 nm,则
$$ \Delta\omega \approx \frac{2\pi c\,\Delta\lambda}{\lambda_0^2} = \frac{2\pi \times 3.0\times10^8 \times 80\times10^{-9}}{(800\times10^{-9})^2} \approx 2.36\times10^{14}\ \mathrm{rad/s} = 0.236\ \mathrm{rad/fs}. $$
代入上式,2 mm 熔石英在 80 nm 带宽内造成的群延迟变化约为 $72.4\times0.236\approx17.1$ fs。要把它抵消,涂层需要约 −71 fs² 的 GDD。这个数字就是后面器件的设计目标:在 80 nm 带宽内提供接近恒定的负 GDD,让不同频率的成分重新对齐;单纯把色散数字做大没有意义。

图 1:纳米涂层压缩器的工作原理与设计。(a) 透明材料色散对飞秒脉冲的影响,以及涂层的两种用法:镀在薄衬底上压缩已经展宽的脉冲,或镀在厚光学元件上补偿其群延迟色散;白色为短脉冲,彩色为被材料色散展宽的脉冲,浅蓝为正常色散材料,紫色为涂层。(b) 针对 800 nm 优化的设计。(c) 垂直于传播方向的二维截面中主要模式的色散:HE11-like(蓝线)与 HE12-like(红线),并给出真空(灰线)与硅(紫线)的色散;HE12-like 模在 0.36 PHz(825 nm)处 $k_z=0$,接近截止时群速度下降,形成异常群速度色散(红箭头)。(d) 自由空间到两种模式的功率耦合比例;在 0.36 PHz 处耦合由 HE12-like 转向 HE11-like,从而拓宽工作带宽;绿色为被反射的功率。(e) 800 nm 处两种模式的横向电场分布,HE11-like 没有拐点,HE12-like 有一个拐点。
图 1a 把两种用法画在一起:涂层可以单独镀在薄衬底上,作为压缩器;也可以镀在分束器、波片、窗口这类厚元件上,抵消元件本身的正常色散。两种用法对涂层的要求相同,都是让红光成分比蓝光成分多走一点等效光程。区别只在于被补偿的色散来自哪里。
慢光从哪里来
器件的结构很简单:0.5 mm 厚熔石英衬底上立着一层 610 nm 厚的非晶硅纳米柱,柱子按 475 nm 周期的方阵排列,直径在 147–170 nm 之间。这个周期比光在空气和衬底中的最短波长都小:760 nm 的光在熔石英里对应约 524 nm,而阵列周期只有 475 nm,所以除零级以外所有衍射级都是倏逝波,光束的横向轮廓不会被阵列改写。涂层改变的是相位和群延迟,不把光散到别的方向。
要理解它为什么能给出负 GDD,得看这个二维阵列支持的两个模式。论文把它们按介质波导的混合模式命名为 HE11-like 与 HE12-like。HE11-like 模主要被限制在硅柱里;HE12-like 模有一部分能量在柱外的自由空间中传播,因此在真空光线以上仍能有实的传播常数。关键在 825 nm 附近:HE12-like 模的 $k_z$ 降到零,进入截止。根据互易性,光在正向和反向传播时速度相同,所以 $k_z\to0^{+}$ 时群速度必须趋于零:
$$ v_g=\frac{d\omega}{dk_z}\to0^{+}. $$
这就是慢光。群速度在截止点附近迅速下降,群延迟随频率剧烈变化,对应的群速度色散(group velocity dispersion,GVD)为
$$ \mathrm{GVD}=\frac{d^2k_z}{d\omega^2} =\frac{d(v_g^{-1})}{d\omega}. $$
在图 1c 里,HE12-like 色散曲线从高频侧离开真空光线,弯向 $k_z=0$ 的慢光区,这段弯曲就是异常 GVD 的来源。
如果光只耦合到 HE12-like 一个模式,群延迟会在截止点发散,器件就不可能宽带工作。实际结构在截止点附近把耦合迅速转给 HE11-like 模,如图 1d 所示。两个模式的混合把相位曲率摊平,得到一段近似恒定的异常 GDD。论文认为,关键在模式耦合的转移,不在单一谐振的相位变化;这也是它能在 80 nm 带宽内保持群延迟近似线性的原因。
图 1b 给出最终设计的截面:柱高 610 nm,周期 475 nm,柱径在 150 nm 量级。图 1c 是模式色散,图 1d 是自由空间到两个模式的功率耦合比例,图 1e 是 800 nm 处的横向电场分布。把这三张分图连起来看,逻辑是:入射光同时耦合到两个模式,HE12-like 模在接近截止时群速度下降,提供异常色散;截止点之后耦合转向 HE11-like 模,避免群延迟发散。设计要调的,就是让这个转移发生在工作波段内,并且转移得足够平缓。
设计、实测与误差
设计分两步。先用严格耦合波分析在柱径、柱高、周期的大参数空间里找候选,再换到时域有限差分(finite-difference time-domain,FDTD)细调,目标同时包括大的负 GDD、低高阶色散和高透射。FDTD 给出的设计窗口是 760–840 nm:全范围内平均 GDD 为 −71 fs²,群延迟偏离线性小于 2 fs;如果只看 780–825 nm,平均 GDD 为 −82 fs²,偏离小于 0.5 fs。仿真还预测全范围透射接近 1,入射强度超过 50 GW/cm² 时双光子吸收仍低于 1%。这几个数字放在一起,说明设计不是用低透射换大色散。
样品用电子束光刻加电感耦合等离子体反应离子刻蚀做成,非晶硅层厚 610 nm,衬底是 0.5 mm 熔石英,演示面积 3.1 mm²。论文特别提到,所用材料与特征尺寸可以由商用深紫外浸没光刻量产;这意味着它具备走晶圆级涂层工艺路线的条件。
实测用白光迈克尔逊干涉仪。钨灯的白光经过准直和起偏后分成两臂,一臂两次通过样品,另一臂由延迟线扫描;两臂的光谱干涉条纹由光栅光谱仪记录,条纹的相位随频率的变化就给出样品的群延迟。三片不同柱径的样品覆盖了不同的中心波长,结果列在表 1 中。
表 1:三片压缩器的实测参数。GDD 为对应波段的平均值,透射率为相对熔石英衬底的平均值;工作波段由柱径调谐。
| 柱径 (nm) | 全工作范围 (nm) | 全范围平均 GDD (fs²) | 线性范围 (nm) | 线性窗口 GDD (fs²) | 平均透射 |
|---|---|---|---|---|---|
| 147 ± 7 | 705–775 | −71 ± 1 | 735–770 | −120 ± 2 | 81% |
| 158 ± 7 | 740–810 | −61 ± 2 | 775–805 | −128 ± 6 | 79% |
| 170 ± 6 | 770–850 | −64 ± 2 | 800–840 | −127 ± 4 | 80% |
表 1 里最需要分清的是“全范围”和“线性窗口”两列。全范围平均 GDD 在 −61 到 −71 fs² 之间,对应 70–80 nm 的可用带宽;线性窗口更窄,但 GDD 更大,达到 −120 到 −128 fs²。窄窗口来自一片被抑制的法布里—珀罗谐振:它在透射谱上留下不到 1 nm 宽的凹陷,压窄了高频侧的线性范围,却也把该窗口内的异常 GDD 推高。低频侧另有一个宽带透射凹陷,把平均透射拉到约 80%。把 4 nm 的柱径加工容差代入仿真,可以复现这些透射特征;论文据此判断,改进加工后透射还有接近 1 的空间。

图 2:压缩涂层的 FDTD 模拟结果。(a) 透射相位(红线)以及全工作范围(浅绿)与线性工作范围(深绿)的抛物线拟合,拟合标注的曲率就是群延迟色散;向下弯曲表示异常色散。(b) 模拟的群延迟(红线);黑线为线性窗口的线性拟合,全工作范围内群延迟偏离线性小于 2 fs;紫色箭头标出被抑制的法布里—珀罗谐振留下的窄带残余。(c) 镀膜衬底表面的透射(红)、反射(蓝)与线性吸收(绿);入射强度超过 50 GW/cm² 时双光子吸收仍低于 1%,灰线为未镀膜表面的透射。
图 2 把“相位曲率就是 GDD”画得很直接。图 2a 中红线的向下弯曲越明显,二阶色散越大;两条拟合线分别对应全范围和线性窗口。图 2b 是相位的导数,也就是群延迟,斜率就是 GDD;黑线越接近直线,高阶色散越小。图 2c 则给出透射、反射和吸收的三条曲线,说明损耗主要来自哪里。读这三张图时要注意,仿真的“接近 1 透射”是设计值,不是实测值;实测的平均透射在表 1 中是 79%–81%。

图 3:不同工作波长压缩器的白光干涉测量。(a) 柱径 158 nm 样品的扫描电子显微镜图。(b) 白光迈克尔逊干涉仪:一路光两次通过样品,另一路被延迟,两臂的光谱干涉由光栅光谱仪记录。(c) 柱径 170 nm 样品的原始数据:干涉谱(绿)、参考臂挡住时透过样品的光谱(蓝)、用熔石英衬底替代样品时的光谱(灰)与灯谱(红)。(d)–(f) 白光干涉测得的群延迟(蓝)与计入 4 nm 柱径容差的 FDTD(红);黑线为线性窗口的拟合,标注其 GDD;浅绿为全工作范围,深绿为线性范围;e 中红色区域是涂层给出正 GDD 的波段,灰线为熔石英衬底的群延迟。(g)–(i) 相对熔石英衬底的实测(蓝)与模拟(红)透射率。
图 3 提供了三条相互独立的证据。图 3a 的扫描电子显微镜(scanning electron microscope,SEM)图给出柱阵的实际形貌,柱径 158 nm 的样品与设计量级一致;图 3b、3c 说明群延迟是从干涉条纹的相位差反演出来的,不是从仿真曲线抄来的;图 3d–3f 则把实测(蓝)与计入 4 nm 容差的仿真(红)叠在一起,三片样品的色散量级和斜率都对得上。图 3g–3i 的透射曲线解释了为什么平均透射只有 80%:低频侧有一个宽凹陷,不是单点坏掉。把图 2 到图 4 连起来看,设计、结构、光谱、时域四条证据链比单独一个漂亮数字更有说服力。
把 48 fs 压到 37.5 fs
色散曲线只是频域证据,真正要回答的是脉冲有没有变短。论文把柱径 162±6 nm 的样品插进一台锁模钛宝石振荡器的光路,用二次谐波频率分辨光学门(second-harmonic frequency-resolved optical gating,SH-FROG)测量时域结构。
SH-FROG 的做法是:把待测脉冲分成两束,一束相对另一束延迟,再让它们在 β-硼酸钡(beta-barium borate,BBO)晶体里非共线重叠。只有当两束脉冲同时到达晶体时,才会产生二次谐波;把不同延迟下的二次谐波谱记录下来,就得到一张二维谱图。谱图里同时包含脉冲的强度包络和相位信息,再用迭代叠层(ptychographic)反演算法把脉冲还原出来。论文的反演误差 $G'<0.05$,说明谱图与还原结果吻合得不错。
这里有一个测量上的细节。二次谐波过程混合的是两个相同脉冲复制品中的光子,所以谱图对延迟是对称的,$S(\tau)=S(-\tau)$;仅凭谱图,GDD 的正负号会有歧义。论文的处理办法是让入射脉冲本身的 GDD 量值大于涂层提供的 GDD,这样补偿过程中不会发生符号翻转,反演结果就不会落在错误的分支上。论文在正文里交代了这个条件。
时域结果是:只经过熔石英衬底时,脉冲半高全宽为 48.3±1.2 fs;加上涂层后变为 37.5±1.3 fs,缩短约 11 fs。按这两个数算,缩短比例是 $(48.3-37.5)/48.3\approx22\%$。论文正文用“四分之一”描述这一步;按给出的数字算是 22%,本文保留原始数字。压缩没有把脉冲压到变换极限;论文指出激光带宽超出了压缩器的线性工作范围,并另外给出一个带宽匹配的数值外推。
数值外推的设定是:一个 10 fs 的高斯脉冲先被 1.6 mm 熔石英展宽到 19 fs,再让测得的涂层色散作用上去。结果是去掉超过 65% 的展宽,脉冲回到 13 fs 以下。把 65% 代进去,19 fs 与 10 fs 之间相差 9 fs,去掉 5.85 fs 后剩 13.15 fs;论文说“低于 13 fs”,说明实际补偿比例略高于 65%,与它“超过 65%”的表述一致。这个外推不是实验直接测到的 10 fs 脉冲,正文把它和实测的 48.3→37.5 fs 分开写。

图 4:飞秒激光脉冲的压缩。(a)(b) 仅通过熔石英衬底 (a) 与通过镀膜衬底 (b) 后的 SH-FROG 谱图,样品柱径 162 nm。(c)(d) 反演得到的时域强度(蓝线)与瞬时频率(红线),箭头标出半高全宽;标准偏差由 bootstrap 方法给出。(e) SH-FROG 反演的群延迟(蓝线)与同一片样品的白光干涉结果(红线),以及全工作范围的最小二乘拟合(黑线);灰线为熔石英衬底的 SH-FROG 群延迟。(f) SH-FROG 装置:脉冲经样品后分束、延迟,非共线聚焦到 BBO 晶体,二次谐波由光栅光谱仪记录。
图 4 的读法可以分三层。先看 a、b 两张谱图:补偿后谱图在延迟方向上的能量分布更集中,说明脉冲在时间上更短。再看 c、d 两条时域曲线:蓝线的宽度从 48.3 fs 收到 37.5 fs,红线(瞬时频率)的斜率也变平,说明残余啁啾减小。最后看 e:SH-FROG 反演的群延迟与白光干涉测得的群延迟在不确定度内重合,两条完全不同的测量路径给出同一结果,这比单独一条曲线更可信。
它在激光链路里替换了什么
把这片涂层放回飞秒激光链路,它的位置很清楚:在振荡器或放大器之后、实验样品之前。它不提供增益,不构成谐振腔,也不选择纵模;它处理的是脉冲在传输和经过光学元件时积累的相位。换句话说,它承担的是“色散补偿”这个功能块,不碰激光器的增益或反馈。
传统上,这个功能块由棱镜对、光栅对、啁啾镜或光子晶体光纤承担。棱镜对和光栅对利用角色散,需要让不同频率的光在空间上分开一段距离再合束,光路长、对准敏感。啁啾镜是反射式多层膜,能在相近带宽内给出更大的 GDD 总量,但膜系厚、设计复杂,而且反射式使用会限制光路布局。另一条路线是把超表面放进傅里叶变换(4f)装置:光栅先把光谱在空间上分开,透镜把各频率成分送到超表面的不同位置,超表面逐点改振幅和相位,再由第二个透镜和光栅合束。它能合成更复杂的脉冲形状,但装置仍然依赖角色散和较长的光路。
本篇涂层的优势在于透射式、单片刻蚀、直接插在光路里,不需要空间分光;论文的对比也落在“单位厚度提供的 GDD”上:在相近带宽内,10 μm 量级的啁啾镜能给出更大的色散总量,但按厚度折算,这片 610 nm 涂层的色散密度高得多。
它的边界同样清楚。第一,带宽是 70–80 nm,中心波长在 800 nm 附近,用不同柱径可以移到 705–850 nm 的不同窗口,但一片器件覆盖不了整个可见到近红外。第二,平均透射是 79%–81%,对高功率或弱信号实验,损耗需要计入链路预算。第三,色散在制成后固定:换柱径调中心波长,换柱高调色散量,但这篇工作没有可调电极,不能像空间光调制器那样逐点改相位。第四,演示面积 3.1 mm²,光束口径更大时要扩大加工面积或拼接。第五,入射角在 s 偏振、5° 以内时特性基本不变,p 偏振和更大角度没有给出同等结论。第六,设计仿真的群延迟在 80 nm 内偏离线性小于 2 fs;对更短的少周期脉冲,残余三阶色散会成为下一个要处理的项。
还有一个措辞上的边界。从应用定位看,它接近人们说的“片上色散控制”,但这篇器件是自由空间透射式涂层,镀在熔石英或既有光学元件上,不是把波导、分束器和探测器集成在同一块芯片上的集成光路。更准确的说法是“芯片级厚度的在线色散元件”:它的厚度和加工方式接近薄膜工艺,但光仍在自由空间传播。这个区分不影响它的工程价值,却决定了它不能直接替代波导中的色散工程。
从应用看,它最直接解决的是“现成光学元件的色散拖累”。一片 2 mm 熔石英窗口、一个分束器、一片波片,单独看都不厚,但在 10–50 fs 的脉冲链里,累积的 GDD 足以把脉冲拉长。把涂层镀到这些元件上,或者单独插一片在光路里,就能在不大改光路的前提下补回大部分色散。论文还指出,把工作波段移到长波侧,同一设计可以给出最高 +83 fs² 的正 GDD,用作紧凑的脉冲展宽器;这说明“慢光涂层”不是只能做压缩,同一片涂层的工作波长决定它是压缩还是展宽。
意义与展望
这篇论文把色散控制从“一桌棱镜和光栅”变成“一片 610 nm 厚的涂层”。它的色散由结构的几何参数决定,材料本征色散不是主角,所以换柱径、换柱高就能移动波长和调节色散量;加工又落在电子束光刻和深紫外光刻的工艺范围内,具备做大的可能。对飞秒激光系统来说,这意味着色散元件可以像镀膜一样贴在现有光学件上,而不必每次重新搭一套压缩器。
下一步的推进方向在论文里也写得很具体:单层慢光涂层的 GDD 与带宽受同一个慢光上限约束,80 nm 带宽下理论极限约 −264 fs²/μm,实际器件大约只能用到一半;要提高色散量,可以把多层结构叠起来,或用逆向设计、机器学习去同时裁剪色散曲线和高阶项。更长远的目标从压缩一个高斯脉冲转向合成更复杂的光谱相位,用于相干控制化学反应、量子体系或高次谐波等非线性过程。这些方向都还没有在这篇 2021 年的器件上实现,但它们的起点,是同一片能把红光和蓝光重新对齐的纳米柱阵列。
参考资料
- M. Ossiander, Y.-W. Huang, W. T. Chen, Z. Wang, X. Yin, Y. A. Ibrahim, M. Schultze, F. Capasso. Slow light nanocoatings for ultrashort pulse compression. Nature Communications 12, 6518 (2021). 论文原文
- Author Correction: Slow light nanocoatings for ultrashort pulse compression. Nature Communications 12, 7076 (2021). 更正说明
- S. Divitt, W. Zhu, C. Zhang, H. J. Lezec, A. K. Agrawal. Ultrafast optical pulse shaping using dielectric metasurfaces. Science 364, 890–894 (2019). Science 论文
- R. Paschotta. Group Delay Dispersion. RP Photonics Encyclopedia. RP Photonics 词条