文献解读 | 双相位超表面算子：一块被动超表面完成微分、边缘检测与复振幅全息

一句话结论：复振幅拆成两套纯相位，靠偏振复用干涉还原，一块 1.8 mm 见方的 TiO₂ 超表面即可完成微分、边缘检测与体全息。

论文简介

这篇论文是 Linzhi Yu、Haobijam Johnson Singh、Jesse Pietila、Hümeyra Çağlayan 等人发表在 Light: Science & Applications（2026 年，卷 15）上的工作《Double-phase metasurface operators for all-optical image processing》，DOI 为 10.1038/s41377-025-02153-w。作者主要来自芬兰坦佩雷大学，合作单位包括荷兰埃因霍温理工大学。论文报告了一种可在可见光波段工作的“元算子（meta-operator）”：把双相位编码与偏振复用结合，让单块被动超表面承担多种图像处理运算，并顺带做复振幅体全息。

说它新在哪儿：过去一块超表面通常只对应一种运算。这项工作用一套编码方法，把多个复振幅滤波器装进同一块纳米结构，并在 4f 系统里把微分、互相关、角点与边缘增强的实测图像拍了出来。

图 1：超表面全光图像处理与复振幅全息的示意，涵盖边缘检测、目标检测与体全息生成。

全光图像处理卡在哪

图像处理的算力需求一直在涨。电子电路做图像变换时，光电信号要先完成模数转换，像素再被依次搬运、逐级处理；转换本身有延迟，芯片的热耗散又随着晶体管尺寸逼近极限而越来越难压。于是人们重新盯上光学：光本身是并行的，一张图像可以整体投射到一个光学元件上，无需逐像素扫描，也不需要模数转换。

想法不新。早在几十年前，用傅里叶光学做图像滤波就已经写进了教科书。真正卡住它的是体积。传统方案要透镜、滤波器、反射镜一堆分立元件，系统占一张光学平台，机动性几乎为零。

超表面提供了另一条路。它用亚波长尺度的“超原子”阵列调光，厚度不到波长量级，却能在亚波长分辨率上同时操控相位、振幅和偏振；望远镜、光谱仪、干涉仪都已经有超表面版本。把它放到傅里叶面当滤波器，理论上很自然。

问题出在“一个器件只干一件事”。论文梳理了已有工作：大多数超表面图像处理器要么只实现单一功能，要么依赖庞大的外围光路，要么工作在长波（太赫兹、红外）。可见光波段用得最广，与人眼视觉、高分辨率成像和显微系统都兼容，却偏偏少有人用单块超表面做出多个算子。这项工作要回答的，就是能不能在可见光、在单层结构里，把多种图像变换统一起来。

双相位编码：把复振幅拆成两套纯相位

要在傅里叶面实现任意滤波器，得同时控制光场的振幅和相位，也就是复振幅调制。超表面的麻烦在于：单个纳米柱能高效调的是相位，振幅往往被固定在接近 1。

论文采用的解法叫双相位编码（double-phase encoding），思路来自 1978 年的一篇全息文献。一个振幅为 \(A(x,y)\)、相位为 \(\theta(x,y)\) 的复振幅 \(U(x,y)\)，可以写成两个等幅纯相位分量的叠加：

直觉上，两个分量长度相同、方向关于 \(\theta\) 对称地偏开 \(\pm\psi\)。它们相加时，横向分量互相抵消，纵向分量相加，合矢量长度由夹角决定。振幅变小，就相当于把两个分量掰得更开一些；相位旋转，整对分量跟着一起转。于是“复振幅”这件事被翻译成了两套纯相位分布 \(\theta+\psi\) 与 \(\theta-\psi\)。

翻译完成后，剩下的任务是让超表面分别、独立地施加这两套相位。论文的做法是偏振复用：把两套相位分配给两个正交偏振通道，再让它们干涉还原出目标复振幅。

图 2：双相位编码超算子的工作原理与超原子设计。(a) 基于偏振复用的双相位编码图像处理算子示意；(b) 复振幅调制策略的概念；(c) TiO₂ 矩形纳米柱超原子；(d) x 偏振入射下透射率与相移随宽、长的变化；(e) y 偏振入射下的对应结果；(f) 从 LCP 到 RCP 的偏振交叉转换效率随宽、长的变化；(g) 长度 235 nm、宽度 115 nm 时交叉转换效率与几何相位随取向角（0°–180°）的变化。

超原子设计：TiO₂ 纳米柱如何同时管两种偏振

承载这两套相位的，是玻璃基底上一排排矩形的 TiO₂ 纳米柱。论文给出固定的几何骨架：柱高 \(H=600\) nm、周期 \(P=450\) nm，宽 \(W\) 与长 \(L\) 在 100–400 nm 之间可调。矩形截面带来结构双折射，x 偏振与 y 偏振看到不同的等效折射率，因此同一个柱子能给两个正交偏振施加两套独立的相移。

参数是这样定的。工作波长 532 nm，周期 450 nm 小于波长，阵列处在亚波长区，正常入射时不会激起除零级以外的高阶衍射。TiO₂ 由椭偏测得折射率 2.47、消光为零，与玻璃基底 1.53 形成足够高的折射率差，600 nm 的柱高约为 1.13 个波长，相位累积量做得开。仿真用 COMSOL 的射频模块扫宽和长，自由空间传播用角谱法，整套设计在 MATLAB 里完成。

线性偏振复用下，把纳米柱取向角设为 0，x 与 y 两个通道就能分别拿到 \(\theta+\psi\) 与 \(\theta-\psi\)，干涉后还原复振幅。

圆偏振复用走的是另一条物理路线。此时纳米柱充当半波片，尺寸固定、靠取向角调相位。当长短轴相移差满足 \(\phi_y-\phi_x=\pi\) 时，一束线偏振光被分解成 RCP 和 LCP 两个圆偏振分量，它们拿到的相位延迟大小相等、符号相反，这就是庞加莱球上的几何相位（Pancharatnam–Berry phase，PB 相位）。论文扫描后指出，长 235 nm、宽 115 nm、高 600 nm 的柱子交叉偏振转换效率高，且几何相位随取向角 0–180° 平滑变化。哪些算子适合线性、哪些适合圆偏振，取决于算子在空间反演下的对称性：二阶微分这类算子是偶对称的，交给 PB 相位更顺手。

一阶微分与互相关：把算子搬进傅里叶面

微分在图像里对应的就是边缘。数学上，空间域的一阶微分在傅里叶域只差一个因子：

这意味着求导可以整体搬进频域完成：先做一次傅里叶变换，把每个空间频率乘上一个复系数，再逆变换回来。透镜本身就是傅里叶变换器，问题只剩下“在傅里叶面放一个什么样的复滤波器”。论文的元算子就在这个位置，用双相位编码把它做成一块超表面。

实验用 4f 系统验证。输入图像先经过起偏器，得到等量的 x、y 偏振分量；超表面放在 4f 系统的傅里叶面；出射端再加一个检偏器，让两个偏振分量干涉，把复振幅调制转换成可拍的强度。测试图案是一维条带和二维轮辐。经过一维微分算子后，条带的横向边缘被显著提亮，截面强度曲线在边界处出现尖峰；二维算子则让轮辐图案在各个方向上都显边。

互相关是同一思路的延伸。识别某个模板，等价于把图像与翻转后的模板做卷积，在频域就是乘以模板频谱的复共轭：

论文为字母 T、A、U 各设计了一块元算子。把含“TAU”三个字母的图像送进去，对应算子只在目标字母的位置留下一个亮斑，其余位置被压低。相关运算在光速下并行完成，不需要数字后处理。

图 3：基于线性偏振复用的一阶微分元算子。(a) 一维一阶微分算子的光学显微与 SEM（俯视与倾斜）图像；(b) 二维一阶微分算子的对应图像；(c) 一维算子在 x、y 偏振照明下的相位分布；(d)(e) 条带图案处理前后的图像；(f) 二维算子在两种偏振下的相位分布；(g)(h) 轮辐图案处理前后的图像；(i)(j) 分别对应 (d)(e) 与 (g)(h) 的截面强度曲线。标尺：光学图像 200 μm，SEM 图像 1 μm，实验图像 200 μm。

二阶微分：换圆偏振，做顶点检测与拉普拉斯

一阶微分强调梯度，二阶微分则对曲率和更精细的结构变化敏感。论文选了两种算子：顶点检测与拉普拉斯。

顶点检测用的是混合二阶导数 \(\partial^2/\partial x\partial y\)。它只对 x、y 两个方向同时有强度变化的成分加权，因此对线条的交点和角点特别敏感，一条直边经过它会被抑制，拐角反而亮起来。拉普拉斯算子 \(\nabla^2=\partial_x^2+\partial_y^2\) 是各向同性的高通滤波器，它不挑方向，把所有高频细节一起抬起来，效果是锐化边缘、增强对比。

这两个算子都放在圆偏振复用下实现。原因在上一节的对称性讨论里：\(\partial_x\partial_y\) 与 \(\nabla^2\) 在空间反演下是偶运算，正好匹配 PB 相位对 RCP、LCP 施加的相反相移。实验里，十字图案经过顶点检测算子后四个角被明显提亮；汉字“YU”经过拉普拉斯算子后笔画边缘被锐化，截面上出现陡峭的高频起伏。论文还提到，圆偏振复用路线对波长变化有较强的容忍度，这在实际系统中是个不坏的附带好处。

复振幅全息：同一套编码延伸到体三维重建

图像处理之外，论文把同一套偏振复用编码用到了全息上。全息要重建的是完整的复振幅光场，振幅和相位缺一不可，这恰好是双相位编码擅长的事。

演示器件重建的是一个螺旋点阵，点在 0–800 μm 的深度范围内以 8 μm 为间隔排布。0 到 800 之间有 100 个间隔，对应 101 个深度面；每个面上的强度分布由数值反向传播确定，再经双相位编码写成 x、y 两个偏振通道的相位分布。器件由 1000×1000 个超原子组成，物理尺寸 450 μm × 450 μm。模拟与实验拍到的不同深度强度分布吻合，说明这套编码不仅能做二维滤波，也能把复振幅信息投到三维空间里。

论文给出的调制分辨率下限是 450 nm，也就是一个超原子的周期。更细的结构要靠继续优化超原子设计去争取。

图 4：基于线性偏振复用的复振幅体全息。(a) 螺旋图案复振幅全息器的光学显微与 SEM 图像；(b) x、y 偏振照明下的相位分布；(c) (d) 空间不同深度处重建体全息的模拟与实验强度分布。标尺：光学图像 50 μm，SEM 图像 1 μm，全息图像 50 μm。

实验实现与关键参数

器件制造走的是标准纳米加工路线：玻璃基底上先沉积 600 nm 厚的 TiO₂ 薄膜，甩电子束光刻胶、电子束曝光显影出二维掩模图形，再沉积铬做硬掩模，用反应离子刻蚀把图形转到 TiO₂ 里，最后去掉铬层。图像处理实验用 4f 系统，元算子放在傅里叶面，入射端和出射端各有一片偏振片；全息实验则用准直激光照射，经物镜和管镜把空间中的重建强度成像到传感器上。

表 1：论文报告的关键结构与工作参数。

这里可以做一次量纲与数量级的自查，确认参数彼此自洽。图像处理阵列的物理边长是 \(4000\times450\ \text{nm}=1.8\times10^{6}\ \text{nm}=1.8\) mm；全息阵列是 \(1000\times450\ \text{nm}=4.5\times10^{5}\ \text{nm}=450\) μm，都与论文给的毫米/微米尺寸对得上。全息深度范围 \(800\ \mu\text{m}\div8\ \mu\text{m}=100\) 个间隔，含两端共 101 个深度面。周期 \(P=450\) nm 小于工作波长 \(\lambda=532\) nm，比值 \(P/\lambda\approx0.85<1\)，阵列处于亚波长区，这也解释了为什么可以只考虑零级透射。

意义与展望

这项工作的价值，在于把“一个超表面只做一件事”推进到“一个设计框架承载多种复振幅算子”。微分、互相关、角点检测、各向同性锐化，本质上都是频域里的线性滤波，只要能把目标复传递函数写成两套纯相位，就能用同一层纳米柱实现。对于机器视觉、无损检测、显微成像这类希望把运算搬到光路里、减轻后端数字负担的场景，这条路子有吸引力；论文也点到半导体光刻缺陷检测、定量相位成像和表面计量这些具体方向。

从原理验证走向实际应用，还有几处需要推进。当前每块器件对应固定算子，想要可编程得引入液晶、相变材料或电光调谐；工作点集中在 532 nm 附近，宽带成像下的色散需要单独处理；论文的效率与分辨率细节放在补充材料，数据也未公开，第三方独立复现的完整链条尚未建立。不过被动、单层、可见光这三个特征本身已经把系统体积压得很低，剩下的更多是工程化问题，而非原理障碍。

常见误解

误以为这块超表面像一块光学 GPU，能执行任意计算。 它做的是线性滤波：把输入图像的频谱乘以一个固定的复传递函数。算子内容在设计和加工时就固化进相位分布，它不存储数据、不迭代，也不会自己“学习”。（S1）

误以为双相位编码需要两层器件。 两套相位由同一层 TiO₂ 纳米柱的两个正交偏振通道分别承担，出射端用检偏器让它们干涉，复振幅才被还原出来。整块器件是一层平板。（S1）

误以为超表面做图像处理可以随心切换功能。 论文演示的每块器件对应一个固定算子，换算子就要换一块器件。要做到同一块器件可重构，还得引入液晶、相变材料或电调机制，这是论文指出的后续方向。（S1）

误以为“全息”只是平面的另一个名字。 这里重建的是复振幅光场在空间中的三维分布：螺旋点阵分布在 0–800 μm 的深度范围内，需要同时控制振幅和相位，才能让不同深度的点各自对焦。（S1）

参考资料

1. Yu, L., Singh, H. J., Pietila, J., et al. Double-phase metasurface operators for all-optical image processing. Light: Science & Applications 15 (2026). https://doi.org/10.1038/s41377-025-02153-w

2. Meta-operators: all optical and wireless image processing via metasurfaces. Light: Science & Applications (2026). https://doi.org/10.1038/s41377-026-02318-1

3. Optical computing metasurfaces: applications and advances. Nanophotonics (2024). https://doi.org/10.1515/nanoph-2023-0871

4. Hsueh, C. K., Sawchuk, A. A. Computer-generated double-phase holograms. Applied Optics 17, 3874–3883 (1978). https://doi.org/10.1364/AO.17.003874

5. Double-phase metasurface operators for all-optical image processing, Europe PMC 开放全文（PMC12929636）. https://europepmc.org/article/PMC/PMC12929636

6. 论文原图来源页（Springer Nature，CC BY 4.0）. https://link.springer.com/article/10.1038/s41377-025-02153-w