文献解读 | 时空调制超表面：用非均匀泵浦光瞬态写下波前

论文简介

• 论文：Spatiotemporal Metasurface for Ultrafast All-Optical Wavefront Shaping
• 作者与机构：Mert Akturk、Giulia Crotti、Andrea Schirato、Katsuya Tanaka、Isabelle Staude、Thomas Pertsch、Giulio Cerullo、Giuseppe Della Valle、Margherita Maiuri；米兰理工大学（Politecnico di Milano）、莱斯大学（Rice University）、耶拿弗里德里希·席勒大学（Friedrich Schiller University Jena）
• 原文链接：arXiv:2608.00308：https://arxiv.org/abs/2608.00308 （2026-07-31 投稿，physics.optics，许可 CC BY 4.0）
• 一句话贡献：用一束空间上不均匀的飞秒泵浦光，在一片本来各单元完全等效的纳米线阵列上临时写出相位梯度，让探测光束在不到 1 ps 内散焦，并在大约 600 fs 后基本复原。

常规超透镜把相位剖面在加工阶段就刻进几何结构，出厂之后这枚器件能做什么已经定死。这篇工作追问的是另一件事：超表面的功能能不能由另一束光临场决定，而且决定得足够快？

动态超表面的老问题：改的是参数，不是相位分布

超表面（metasurface）是把亚波长尺度的谐振单元按周期或准周期排布成的准二维结构，几百纳米厚就能让光在振幅、相位、偏振上出现设计好的空间分布。它最成功的应用之一是波前调控：让每个单元提供不同相位延迟，整片结构合起来就等价于一个透镜或一块任意相位屏。

麻烦在于”出厂即固定”。围绕动态超表面，研究者试过很多手段：热调、机械拉伸、电学栅控、相变材料。它们能把调制速度拉开好几个数量级，但大多数方案改的是预先设计好的那个功能的参数——透射强度、谐振峰位、焦距大小——而相位剖面本身还是加工时就写好的。换句话说，外部刺激帮你把透镜调焦，却没法凭空换一种透镜。

全光方案瞄准的是速度上限。用超短泵浦脉冲激励材料里的”热载流子”（吸收光子后尚未弛豫到带边的非平衡电子-空穴对），折射率能在飞秒到皮秒的时间尺度上变化，对应的重构速率可以进入百 GHz 甚至更高的量级。

本文换了一条更直接的路线：不去调参数，而是用一束空间强度不均匀的泵浦光，在一片原本均匀周期性排列的超表面上直接”画”出一个介电常数梯度。这个梯度再被共振结构翻译成相位梯度，探测光束经过时就看到了一个此前不存在的柱面透镜。

先要找到只动相位的波长

材料被光激发后，复介电常数一变，透射光的振幅和相位通常一起变。想拿这种变化当透镜用，就必须躲开那些”振幅跟着乱动”的波长。论文的设计策略和静态超透镜的 lookup table 思路如出一辙：先做光谱扫描，看清”一份材料扰动会被结构翻译成多少相位变化、多少振幅变化”，再挑一个相位响应大、振幅响应小的窗口来工作。

论文用的是氢化非晶硅（a-Si:H）纳米线阵列，做在二氧化硅基底上，周期 P = 400 nm、线宽 w = 180 nm、线高 h = 220 nm。未激发时透射谱上有一条约 670 nm 的共振，表现为透射强度下陷、相位随之拐折。作者先用一个虚构的均匀扰动 Δε̄ = −1 + i0.5 代表强激发态（对应几 mJ/cm² 量级的泵浦通量），算出扰动前后的差分透射 ΔT 和差分相位 Δφ，两者都随波长剧烈变化，量级和符号都在变。

在共振红翼上的 691 nm，两个量刚好走到一个特殊组合：Δφ 逼近负的极值，而 ΔT 几乎归零——论文借”等吸收点”（isosbestic point）的说法描述这个位置。材料扰动在这里几乎全部转成相位变化。Δφ 取负值意味着相位延迟变小，被激发区域的透过相当于经过了一段更”薄”的光程。

图 1：光致动态金属透镜的设计。(a) 非均匀泵浦照明及其对探测波前影响示意图；下插图为超表面单元的 xz 截面，由二氧化硅基底上的 a-Si:H 纳米线构成，几何参数 P = 400 nm、w = 180 nm、h = 220 nm；上插图为泵浦激发示意。(b) 未扰动条件（黑线）与施加均匀介电常数扰动 Δε̄ = −1 + i0.5（蓝线）时超表面的模拟透射强度 T（实线，上）与相位 φ（虚线，下）；竖线对应后续分析的波长 λ₁ = 691 nm。(c) 扰动与未扰动条件之间的差分透射强度 ΔT（红色实线，左轴）与差分相位 Δφ（绿色虚线，右轴）。(d) 空间非均匀介电常数扰动 Δε(x)（黄色为实部、紫色为虚部）随超表面平面内 x 坐标的变化。(e) λ₁ 下距超表面 2 µm 处透射探测场的相位 φ 空间分布：未扰动（黑线）保持平坦，施加 (d) 所示 Δε 后（蓝线）呈高斯形状。

图 1(c) 里那条曲线值得多看两眼：同一个 Δε̄ 在不同的探测波长下，被结构翻译成的 ΔT 与 Δφ 可以差出很多倍，还能反号。超表面在这里扮演了一个”波长选择性转换器”的角色，这也是后面所有实验得以成立的前提。

让光斑自己撑开

波长选定之后，接下来把空间维度加进去。实验里泵浦光被柱透镜整形成强像散的椭圆光斑：沿周期方向的强度半高全宽只有 45 µm，沿纳米线纵向则有 4 mm，于是激发在一条方向上是陡峭的高斯型剖面，在另一方向上几乎均匀。

仿真把这份非均匀激发写成沿 x 方向的高斯型介电常数扰动 Δε(x)，半高全宽同样取 45 µm。有限口径的计算包含 300 µm 宽、超过 730 个元胞，得到的透射相位剖面不再平坦：中心处相位延迟减小最多，向两翼逐渐恢复，整体呈上凹的高斯形状。中心”更薄”、边缘”更厚”，这片结构在实验中就等价于一枚发散柱面透镜，透射光束因此散焦。

45 µm 的泵浦光斑跨越了约 110 个周期单元，相位梯度是缓变的。这个尺度关系很重要——扰动是铺在一片连续等效的介质上，而不是只砸中几个孤立单元，光束看到的因此是一枚透镜，而不是一层粗糙的散射体。

实验用 400 nm、约 100 fs 的泵浦（光子能量约 3.1 eV，远高于硅的带隙），通量 3.5 mJ/cm²；探测光是一束 70 µm 的高斯光束，偏振平行于纳米线，信号用中心 690 nm、带宽 10 nm 的滤光片提取。成像结果显示，泵浦到达之前（t = −1.4 ps）探测光斑在探测面上的强度半高全宽是 0.28 mm；激发后 0.1 ps，光斑沿 x 方向撑到约 0.44 mm。

图 2：实验装置与强度分布图。(a) 超快成像装置示意图；插图为样品顶视 SEM 图像（水平 xy 面）。(b–c) CMOS 相机在不同泵浦-探测延迟下测得的光强分布：−1.4 ps（b，即未扰动系统）与 0.1 ps（c）。轮廓线标出对应光斑的强度半高全宽，估计为 FWHM ≈ 0.28 mm（b）与 FWHM ≈ 0.44 mm（c）。

这里可以做一次简单的量纲与比例核对。光斑半高全宽从 0.28 mm 变到 0.44 mm，相当于在同一次测量里放大了 0.44/0.28 ≈ 1.57 倍，也就是横向尺寸增加约 57%。按高斯光束的传播关系，光斑在探测面上的尺寸由波前曲率和发散角共同决定，相位剖面上多加一层上凹的柱面分量，直接后果就是远端光斑变宽，数字方向与图 2 中肉眼可见的轮廓变化一致。另一个可复核的换算来自时间尺度：约 600 fs 的恢复时间对应 1/600 fs ≈ 1.7 THz 的特征速率，把有效带宽放在百 GHz 到 THz 之间是合理的量级估计，而不是精确的调制带宽。

为什么必须是超表面

到这一步只剩一个关键疑问：既然原理是材料被光激发后折射率变了，那随便拿一片等厚度的硅薄膜是不是也能做出同样的效果？论文专门做了这个对照实验。同一套光路、同样的泵浦通量、同样厚度的未图形化 a-Si:H 薄膜，测到的差分信号是一条均匀的负值曲线：透射被整体压暗，但沿 x 方向没有任何结构，谈不上波前再分布。

对照结果说明，空间非均匀的光照本身只带来”哪里光强高、哪里吸收多”，要靠结构把这份扰动放大成相位差，才可能改变波前。超表面的共振在这里起到了放大器的作用：同一份介电常数扰动，落在共振附近时能换出远大于平整薄膜的相位改变。

为了把瞬态响应量化，论文定义了一个归一化的差分量：

$$\frac{\Delta I}{I_0}(x,t)=\frac{I(x,t)-I_0(x)}{I_0(0)}$$

其中 I(x,t) 是某延迟时刻沿 y 方向平均后的探测强度，I₀(x) 取激发前的参考剖面，分母 I₀(0) 是束中心的最大值。散焦在这一指标下有一个很好认的指纹：中心强度下降、两侧上升，剖面呈双峰。更重要的是它随时间演化的样子：双峰在泵浦与探测脉冲时间重合时最强，约 600 fs 内衰减殆尽；延迟拉到 1.6 ps 时信号翻成负凹陷，说明此时光致吸收造成的整体衰减盖过了相位再分配。

图 3：超快波前调控。(a–b) 探测波长 λ = 690 nm（实验）与 λ = 691 nm（仿真）下，不同泵浦-探测延迟时超表面的实测 (a) 与模拟 (b) 差分信号 ΔI/I₀ 随 CMOS 相机上空间坐标 x 的变化。(c–d) 光激发超表面的实验 (c) 与仿真 (d) ΔI/I₀ 二维图，横轴为相机上的 x 坐标、纵轴为泵浦-探测延迟。(e–f) 相同探测波长下对等厚度平整薄膜所做对照实验（e）与仿真（f）的对应结果。

图 3 的下半部分把两组数据放在一起对比的效果很直接：上面三格是超表面，实验与仿真的双峰结构和衰减节奏对得上；下面两格是平整薄膜，只剩一条均匀的暗带。

一套把载流子动力学接到电磁仿真上的模型

要让仿真复现这种时空图，光算电磁场不够，还得把材料内部那几百飞秒里发生的事写进去。论文的模型分两层。

第一层是材料侧。泵浦光在硅里通过单光子和双光子吸收产生非平衡载流子，随后载流子经缺陷辅助复合与双分子复合弛豫，每次复合把能量交给晶格，抬升晶格温度。对应的介电常数变化由三部分叠加：双光子吸收贡献一个纯虚部，自由载流子密度 N(t) 贡献 Drude 型项（实虚部都有），晶格温度 Θ(t) 通过热光效应贡献第三项。载流子密度与晶格温度的演化由双温模型（Two-Temperature Model, TTM）描述，参数取值包括陷阱辅助复合时间 30 ps、双分子复合系数 7×10⁻¹⁵ m³/s、晶格热容 1.66 J K⁻¹ cm⁻³，400 nm 处的线性吸收约 0.55，有效双光子吸收系数 5.2×10⁻¹⁰ m/W。

第二层是电磁侧。作者把超表面口径沿 x 离散采样，在每个采样位置单独求解一次 TTM，得到局域的 Δε(x,t)，再喂给与实验同口径的全波仿真；透射场提取出来后，用标量衍射公式沿实验成像光路传播到探测面，最终算出与相机记录的同一个物理量 ΔI/I₀。这样一来，从光子吸收、载流子弛豫、晶格升温，到相位剖面成形、成像系统把相位差转成强度差，整条链条都在同一个模型里闭合。实验与仿真的吻合程度，构成了对上述机制解释最有力的支持。

顺带一提，实验里相机被刻意从标称像面挪开了 15 mm。这并非操作失误：把相机放在离焦位置，可以让瞬态相位变化尽量多地转成强度对比，否则相位差在理想像面上几乎看不出来。

意义与展望

这项工作把相位剖面从加工阶段的固定输入，变成了实验中可以随时写入的变量。论文中已经验证的是散焦这一半，另一半也在同一套框架里：在共振蓝翼约 666 nm 处，ΔT 再次过零而 Δφ 出现正峰，仿真表明此时激发会产生会聚相位，把探测光在皮秒尺度上聚焦。于是一种别致的操作方式出现了，同一片超表面只靠改变探测波长，就能在发散透镜与会聚透镜之间切换。

把视线放远一些，这类时空调制超表面还有几处值得继续推进的方向。实验验证目前是一维周期阵列加像散光斑的简化构型，扩展到二维超表面后，泵浦光的空间分布可以写出的相位图案会丰富得多，理论上能临时拼装出更复杂的波前。材料侧也有余量：非晶硅的载流子寿命决定了约 600 fs 这一恢复窗口，换成复合更快的材料或提高泵浦强度，响应时间还能继续压缩。反过来，超快也意味着功能只在极短窗口内存在，同步精度、泵浦能量和口径尺寸之间需要仔细权衡，这正是从原理验证走向实际应用时最需要打磨的地方。

表 1：三类超表面调控路线的对比。

路线	相位剖面的来源	外部刺激改变的量	响应时间尺度
静态超透镜	加工时写入几何结构	无	不适用
参数可调超表面（热、机械、电学、相变）	加工时写入几何结构	已有功能的强度、谐振位置或焦距	热、机械、电学与相变方案跨越数个数量级
时空调制超表面（本文）	空间结构化的泵浦光临时写入	相位剖面本身（中心变薄，光束散焦）	亚皮秒，约 600 fs 恢复

后续若能把这套光写入相位的思路与可编程泵浦整形结合起来，一幅由光实时绘制的平面光学元件就不再只是概念。那时的超表面更像一块白板，写什么由当次实验决定，而白板的擦除速度，取决于材料里那几百飞秒的热载流子动力学。

参考资料

1. Akturk M., Crotti G., Schirato A., et al. Spatiotemporal Metasurface for Ultrafast All-Optical Wavefront Shaping. arXiv:2608.00308（2026）. https://arxiv.org/abs/2608.00308
2. Della Valle G., Hopkins B., Ganzer L., et al. Nonlinear anisotropic dielectric metasurfaces for ultrafast nanophotonics. arXiv:1706.00869（2017）；Nano Letters 17, 2334–2339. https://arxiv.org/abs/1706.00869
3. Iyer P. P., Karl N., Addamane S., et al. Sub picosecond steering of ultrafast incoherent emission from semiconductor metasurfaces. arXiv:2206.04645（2022）；Nature Photonics 17, 588–593. https://arxiv.org/abs/2206.04645
4. Sadafi M. M., Inampudi S., Mosallaei H. Floquet-Bloch Theory for Dispersive Time-varying Metasurfaces. arXiv:2609.27109（2026）. https://arxiv.org/abs/2609.27109
5. Lippitz M. Optics Lecture Notes, Chapter 2 “Gaussian Beams”（CC BY-SA 4.0）. https://github.com/lippitz/optics-lecture-notes