非局部超表面：偏振切换的相衬与明场双模式成像
这篇论文解读的是哪项工作
本文解读的是 2026 年 6 月发布在 arXiv 上的一项工作：Multifunctional Imaging with an Inverse-Designed Nonlocal Metasurface（arXiv:2606.21051），作者是墨尔本大学的 Lincoln Clark 等 9 位研究者（ARC Centre of Excellence for Transformative Meta-Optical Systems 等机构），v2 版本发布于 2026-06-24。
一句话概括：一片只有 300×300 µm 的介质超表面，通过旋转入射光的偏振，就能在”明场成像”和”相位成像（伪 3D 相衬）”两种模式之间切换——不需要移动任何光学元件，也不需要计算后处理。演示对象是未染色的细胞和组织切片。
引言：透明样本为什么”看不见”
明场显微镜看透明生物样本有一个根本困难：细胞和未染色的组织对光几乎不吸收，强度对比度极低。但样本内部并非”没有信息”——折射率与厚度的微小变化会让透射光带上空间变化的相位，只是普通的强度探测器对相位完全不敏感 [S1][S2]。
把相位变化转成可见的强度对比，是显微术里一个成熟的方向：Zernike 相衬显微镜（1953 年诺贝尔物理学奖）和微分干涉对比（DIC，又称 Nomarski）显微镜都是为此设计的 [S2][S3]。但它们通常意味着庞大的光学系统：相衬需要相位板，DIC 需要偏光棱镜组，而且切换成像模式时要手动把元件插入或移出光路。另一种思路是计算成像（如强度传输方程、数字全息），但往往需要多帧拍摄、机械移动或大量计算 [S1]。
超表面给了一个不同的答案：既然空间频率与入射角一一对应，那么让透射率随入射角变化，就等于在物平面直接对图像做空间滤波——光经过器件的一瞬间，运算已经完成 [S4][S5]。这类”非局域”超表面之前已经实现了边缘检测等单种运算，但大多只能做一件事 [S5]。这篇论文的增量在于：用一片器件同时承载两种运算，用偏振来切换，而且是在正入射下工作。
核心原理：让 OTF 替你做一阶导数
线性、空间不变的成像系统可以用光学传递函数（OTF）$H(k_x,k_y)$ 完整描述：输出场是输入场的空间傅里叶变换乘上 OTF 再逆变换回来 [S1]：
$$U(x,y)=\mathcal{F}^{-1}\left[H(k_x,k_y)\mathcal{F}[U_0(x,y)]\right]$$
其中 $k_x,k_y$ 是空间频率，与平面波的传播方向一一对应（$k_x=k_0\sin\theta\cos\varphi$）。所以”工程化 OTF”就等于”工程化角度相关的透射”——这正是非局域超表面擅长的事 [S1][S5]。
考虑一个特殊的标量 OTF：$H(k_x)=\alpha k_x+\beta$，其中 $\alpha,\beta$ 是实常数。当输入是纯相位对象 $U_0=A e^{i\phi(x,y)}$ 时，输出场和输出强度分别是 [S1]：
$$U(x,y)=U_0(x,y)\left(\beta+\alpha\frac{\partial\phi}{\partial x}\right)$$
$$I(x,y)=A^2\left[\alpha^2\left(\frac{\partial\phi}{\partial x}\right)^2+\beta^2+2\alpha\beta\frac{\partial\phi}{\partial x}\right]$$
这里的关键是交叉项 $2\alpha\beta(\partial\phi/\partial x)$：相位梯度被直接转成了强度变化，而且非零的 $\beta$ 保留了梯度的符号——向上坡亮、向下坡暗，于是产生类似 DIC 的”伪 3D”立体感 [S1]。当 OTF 是常数（$\alpha=0$，即 $H=\gamma$）时，上式退化为均匀强度 $I=\gamma^2A^2$，就是普通的明场成像。
论文把这两个算子同时做进了同一片超表面：
- x 偏振：OTF 为常数 $\gamma=0.8$，执行恒等运算 → 明场成像；
- y 偏振：OTF 为线性函数 $H=\alpha k_x/k_0+\beta$，对相位做一阶导数 → 相衬成像。
目标参数取 $\beta=0.5$、$\alpha=0.5/\mathrm{NA}$，在 $\pm3°$ 范围内（对应 NA = 0.05）优化 [S1]。图 1 是论文原图，展示了这套偏振切换的整体方案。
图 1：论文原图——偏振切换成像方案示意。转动偏振片即可在明场与相衬两种模式间切换。来源：论文图 1（arXiv:2606.21051v2，作者版权，仅作文献解读引用）。
为什么这种方案比传统显微镜值得关注？除了免移动元件，超表面还天然把光学系统压缩到近乎一张薄膜的厚度，并且光运算不消耗能量 [S4][S5]。表 1 给出本文超表面与常规 DIC 显微镜的对比。
表 1：本文超表面与常规 DIC 显微镜的对比。DIC 对比数据来自论文 S1 的实验配置（Olympus BX60，NA = 0.4）。
| 对比项 | 常规 DIC/相衬显微镜 | 本文超表面（演示） |
|---|---|---|
| 模式切换 | 需插入/移除光学元件 | 旋转偏振片即可 |
| 处理位置 | 系统级（物镜/聚光镜组合） | 物平面，紧贴样本 |
| 数值孔径 | DIC 参考配置 NA = 0.4 | NA = 0.05 |
| 对比方向 | DIC 沿 45° 对角线 | 沿 x 轴（导数方向） |
| 器件尺寸 | 完整显微镜系统 | 300×300 µm 单片 |
| 定量能力 | 以定性对比为主 | 可定量恢复相位梯度 |
非对称响应从哪来：衍射 + 全内反射
一个物理难题藏在前面：介质超表面通常在正入射下是左右对称的，而一阶导数 OTF 恰恰需要关于正入射不对称。之前不少介质方案靠倾斜照明来”作弊”，论文的做法是用衍射在正入射下自己制造不对称 [S1]。
设计单元在 x 方向约 850 nm、y 方向约 500 nm（亚波长周期，防止 y 方向产生衍射；论文正文把该数值打印为 µm，与 425×250 像素网格、2 nm 电子束步长和 300 µm 器件尺寸矛盾，按机理应为 nm）。当光正入射时，x 方向的周期结构把一部分能量耦合进 $\pm1$ 阶衍射，衍射角大于石英衬底的全内反射临界角，这些高阶模就在衬底里反复反射、被”困”住，充当波导把能量带走 [S1]。
结果就是：零级透射随入射角出现非对称变化——某个方向的空间频率被削掉更多，另一个方向保留更多，形成近似线性的 OTF。对于 x 偏振，这个衍射响应几乎消失，OTF 恢复平坦，两种模式由此分家。论文图 2 完整给出了这个机制的设计与仿真验证。
图 2：论文原图——器件设计与仿真。a 拓扑优化的最终设计；b 衍射级次示意（±1 级在衬底内全内反射、0 级出射）；c x/y 偏振 OTF 与设计目标对比（RCWA/FEM）；d 衍射效率随空间频率的变化。来源：论文图 2（arXiv:2606.21051v2，作者版权，仅作文献解读引用）。
设计方法：拓扑优化与鲁棒性
这么复杂的自由形状结构显然不是手画的。论文把单元划分成 425×250 个像素（约 2 nm 像素尺寸），让算法直接优化每个像素处的材料密度：0 代表空气，1 代表非晶硅 [S1]。
优化循环是标准的逆向设计套路 [S1]：
- 用严格耦合波分析（RCWA，TORCWA 实现）算出当前设计的角度相关 OTF；
- 损失函数取当前 OTF 与目标 OTF 的均方根误差（RMSE）；
- 用 PyTorch 的自动微分求每个像素密度的梯度，交给 ADAM 优化器更新；
- 通过投影函数逐步把灰度密度压向 0/1 二值化。
工程上更关键的一步是鲁棒性设计：从第 350 次迭代开始，每个优化步同时评估”膨胀”和”腐蚀”两个制造偏差变体（边界位移约 4 nm，与电子束光刻的工艺偏差同量级），损失函数取三个变体的加权平均 [S1]。这样优化出来的结构即使过刻或欠刻几纳米，性能也不会骤变——实测三种变体在 30 nm 波长范围内都保持单调的 y 偏振 OTF。
器件用 300 nm 非晶硅（PECVD 沉积）做在石英衬底上，电子束光刻加 Al₂O₃ 硬掩模、伪-Bosch 刻蚀加工。加工完成的器件用氦离子显微镜（HIM）检查，实测透射谱与 OTF 和设计预期基本一致（图 3）。设计波长是 900 nm（此时硅低损耗、硅探测器仍可用），但 900 nm 处实测 OTF 在高空间频率端偏离线性，所以成像实验在 886 nm 进行 [S1]。
图 3：论文原图——加工后的器件与实测响应。a、b 为氦离子显微镜（HIM）图；c 透射谱随归一化空间频率的变化；d 886 nm 下的实测 OTF。来源：论文图 3（arXiv:2606.21051v2，作者版权，仅作文献解读引用）。
实验结果：从相位靶到生物样本
论文用两种自制相位测试靶验证定量能力：辐条靶（spoke）和牛眼靶（bullseye），在 886 nm 下相位行程约 $\pi$。实测 OTF 标定为 [S1]：
$$t_x=0.954\pm0.006,\qquad t_y=(3.464\pm0.150)\frac{k_x}{k_0}+(0.706\pm0.005)$$
用 x 偏振图像当背景、y 偏振图像当信号，可以定量解出相位梯度（求解式 (6) 的二次方程）[S1]：
$$\frac{\partial\phi}{\partial x}=\frac{1}{\alpha}\left[\sqrt{\frac{\gamma^2 I_y-\beta^2 I_x}{I_x}+\beta^2}-\beta\right]$$
其中 $I_x,I_y$ 分别是两种偏振下的强度图像，$\alpha,\beta,\gamma$ 取实验标定值。该式可化简为 $(\sqrt{\gamma^2 I_y/I_x}-\beta)/\alpha$，即相位梯度直接由两种偏振图像的比值决定。两个测试靶的恢复结果与解析真值对比，RMSE 分别只有 0.0548 rad/µm（辐条靶）和 0.0617 rad/µm（牛眼靶）[S1]。图 4 给出了两个相位靶的成像与定量恢复结果。
图 4：论文原图——相位测试靶成像。a 辐条靶、b 牛眼靶；i 相位对象（AFM 转相位），ii 有无超表面的成像对比，iii 截面强度，iv 恢复的相位梯度。来源：论文图 4（arXiv:2606.21051v2，作者版权，仅作文献解读引用）。
更接近实际的是生物样本实验：未染色的 HeLa 细胞在明场下几乎看不见，切到 y 偏振后细胞位置和形态立刻清晰，左右两侧明暗分明；恢复出的相位梯度图与奥林巴斯 BX60 显微镜的 DIC 图像高度吻合（图 5）[S1]。图 5 中 b、c 分别是商用显微镜的明场与 DIC 图像，d 是超表面恢复的相位梯度，两者对细胞轮廓的呈现一致；DIC 沿 45° 对角线产生对比，而超表面沿 x 轴。对人黏液性卵巢腺癌切片（4 µm 薄切，相邻切片用 H&E 染色作对照），y 偏振还能分辨组织内的囊性间隙、个别细胞以及肿瘤细胞的环形排列——这些在明场下都难以辨认（图 6）。
图 5：论文原图——未染色 HeLa 细胞成像。a 有无超表面的对比；b、c 商用显微镜明场与 DIC；d 沿 x 方向恢复的相位梯度。来源：论文图 5（arXiv:2606.21051v2，作者版权，仅作文献解读引用）。
图 6：论文原图——人黏液性卵巢腺癌切片成像。a 有无超表面的对比；b 相邻切片 H&E 染色；c DIC；d 恢复的相位梯度。来源：论文图 6（arXiv:2606.21051v2，作者版权，仅作文献解读引用）。
论文也坦率地交代了噪声来源：恢复误差主要来自器件小缺陷的衍射伪影，生物样本成像时通过把超表面平移小偏移、取 10 张平均来软化边缘衍射影响（非局域响应在面内均匀，平均不会改变运算本身）[S1]。表 2 汇总了关键结果。
表 2：论文关键实验结果汇总。来源：论文 S1（正文与补充材料）。
| 实验项目 | 对象/条件 | 结果 |
|---|---|---|
| OTF 标定（886 nm） | x 偏振 | $t=0.954\pm0.006$（常数 → 明场） |
| OTF 标定（886 nm） | y 偏振 | $t=(3.464\pm0.150)k_x/k_0+0.706\pm0.005$（线性 → 相衬） |
| 相位梯度恢复 | 辐条靶（$\pi$ 相位） | RMSE = 0.0548 rad/µm |
| 相位梯度恢复 | 牛眼靶（$\pi$ 相位） | RMSE = 0.0617 rad/µm |
| 生物成像 | 未染色 HeLa 细胞 | y 偏振下形态清晰，与 DIC 吻合 |
| 生物成像 | 卵巢腺癌切片 | 囊性间隙、细胞环形排列可见 |
意义与局限
从论文本身的价值看，这项工作的贡献主要有四点。第一，它把”偏振可重构”第一次做进了非局域介质超表面成像器件：旋转偏振片就能切换明场与相衬，免除了显微镜里最容易磨损、最占空间的机械切换机构 [S1]。第二，它解决了介质超表面必须斜入射的历史包袱，正入射工作意味着可以像普通光学元件一样随意插入光路。第三，它展示了拓扑优化在成像器件上的落地能力——不是设计单个单元的相位，而是直接优化整个 OTF，并且把制造容差显式放进目标函数。第四，它不止做了定性展示，还给出了可复现的定量恢复流程（两种偏振图像的比值 + 标定系数），把”相衬成像”升级成了”模拟光学计算”。
局限也需要如实说。最直接的是分辨率：演示器件 NA 只有 0.05，远低于商用 DIC 物镜的 NA 0.4，细胞级细节可以看，亚细胞结构则超出能力 [S1]。其次，当前只计算 x 方向的一阶导数，y 方向信息被丢弃，要对任意样本做定量相位成像还需要双方向导数。第三，方案需要空间相干照明，普通 LED 光源要先经过滤波。第四，对比度与 NA 之间存在原理性权衡：OTF 斜率越大、相衬越强，可用的空间频率范围就越窄，参数需要针对样本尺寸定制。最后，器件缺陷引起的衍射噪声目前靠多帧平均缓解，更好的工艺或更聪明的鲁棒性模型才是根本出路 [S1]。
论文给出的后续方向也值得留意：宽带非局域滤波配合 LED 照明的低成本成像系统、贴近图像传感器的”计算相机”、以及按目标特征定制的 OTF（比如只响应特定细胞形态），这些都在作者自己的展望清单里 [S1]。对平台读者来说，这篇文章和 2026-08-05 的《超表面逆向设计：伴随方法如何让光子器件自动进化》正好构成一对：那篇讲的是”如何让算法自动设计结构”，这篇展示的是逆向设计的产物如何解决一个真实的生物成像问题。
参考资料
- [S1] Clark, L. et al. “Multifunctional Imaging with an Inverse-Designed Nonlocal Metasurface.” arXiv:2606.21051v2 (2026). arXiv 摘要页 · PDF
- [S2] 维基百科：”Phase-contrast microscopy.” Wikipedia
- [S3] 维基百科：”Differential interference contrast microscopy.” Wikipedia
- [S4] Liu, T. et al. “Flat optics for analog computing: from fundamental mechanisms to advanced meta-processors.” arXiv:2604.16849 (2026). arXiv 摘要页
- [S5] Kwon, H. et al. “Nonlocal metasurfaces for optical signal processing.” Phys. Rev. Lett. 121, 173004 (2018). APS 页面