文献解读 | 超材料模拟计算的起点：把数学运算写进一块平板
这篇论文解读的是哪项工作
本文解读的是 2014 年 1 月发表在 Science 上的论文：Performing Mathematical Operations with Metamaterials（DOI 10.1126/science.1242818），作者是宾夕法尼亚大学的 Alexandre Silva、Nader Engheta，德克萨斯大学奥斯汀分校的 Francesco Monticone、Andrea Alù，以及意大利桑尼奥大学的 Giuseppe Castaldi、Vincenzo Galdi [S1]。同期 Science 还配发了 Ari Sihvola 撰写的评论文章 [S2]。
一句话概括：这篇论文提出了”超材料模拟计算”（metamaterial analog computing）的概念——设计一种超材料块，让入射光场的空间轮廓在穿过它时自动完成指定的数学运算，比如空间微分、积分或卷积。 光不用先变成电信号、也不用经过数字芯片，运算在光传播的瞬间完成 [S1]。这是后来几乎所有”超表面光计算”工作的思想起点。
为什么是”模拟计算”：光天生适合做连续量的运算
先厘清一个概念：模拟计算（analog computing）和数字计算（digital computing）的区别。数字计算把信息编码成 0 和 1，通过门电路一步步算；模拟计算直接对连续物理量——这里是光场的振幅和相位——做运算。光场的空间轮廓本身就是一个连续函数 $f(x)$，对它求导、积分正是数学上最”模拟”的一类运算 [S4]。
模拟光计算的吸引力在于两点：其一，运算发生在光穿过器件的瞬间，没有模数转换延迟，也不消耗额外能量；其二，它天然适合图像处理这类”前端计算”任务——摄像头采集到的本来就是连续光场，与其转换成数字再算，不如在光学域直接处理 [S4]。论文的核心贡献，就是把”数学运算”翻译成”光场调控”问题，再交给超材料去执行 [S1]。
核心思想：把数学运算写进传递函数
理解这篇论文，关键是傅里叶光学的一条性质：空间微分在频域就是乘以波矢。
对输入光场 $f(x)$ 做傅里叶变换得到 $F(k_x)$，其中 $k_x$ 是横向空间频率（由入射角决定）。n 阶空间微分的傅里叶变换恰好是 $(ik_x)^n F(k_x)$ [S4]：
$$\frac{\partial^n f(x)}{\partial x^n} = \mathcal{F}^{-1}\left{(ik_x)^n \mathcal{F}{f(x)}\right}$$
也就是说，只要做一个”运算块”，让它在空间频域的传递函数 $H(k_x)$ 恰好等于 $(ik_x)^n$，光场穿过它之后，输出就是输入的第 n 阶导数 [S4]。
图 1 是论文对这套思想的概念示意：一个设计好的超材料块，让任意波信号在穿过它时执行期望的数学运算。
图 1：论文原图——超材料运算块对穿过它的波场执行数学运算的概念示意。来源：论文图 1（Silva 等，Science 343(6167), 2014，仅作文献解读引用）。
于是，”怎么设计会算数的材料”就变成了一个纯粹的电磁设计问题：怎么让一块平板的透射系数随入射角（即横向波矢 $k_x$）按目标函数变化 [S1][S4]。
论文提出的两种实现架构
Silva 等人给出了两条具体路线 [S1]。
路线一：傅里叶域处理器（GRIN 波导 + 超屏）
用渐变折射率（GRIN）波导充当透镜：第一段 GRIN 波导把输入场的空间轮廓变换到频域；中间放一块亚波长结构”超屏”（metascreen），它的透射系数随横向位置变化，等价于在频域乘上目标传递函数 $H(k_x)$；第二段 GRIN 波导再做一次逆变换，输出端就得到运算结果 [S1][S3]。这条路线像把传统的 4f 光学系统压扁成平面波导，能按像素级设计任意复杂的算子。
图 2 是路线一的完整设计：论文用渐变折射率波导夹一块亚波长超屏构成运算器，并给出了一阶/二阶微分、积分与卷积的仿真验证。
图 2：论文原图——路线一（超屏方案）：GRIN 波导夹 metascreen 的示意图，以及一阶/二阶微分、积分与卷积的仿真结果。来源：论文图 2（Silva 等，Science 343(6167), 2014，仅作文献解读引用）。
路线二：格林函数方案（多层平板）
更极致：直接用一块多层平板。平板横向均匀（两个方向材料不变）、纵向分层，通过设计每一层的介电常数和厚度，让平板整体的透射系数随入射角（横向波矢 $k_x$）按目标 $H(k_x)$ 变化 [S1][S7]。光场穿过平板时，每个平面波分量被乘上对应的系数，输出天然就是运算结果——不需要任何透镜。这正是后来”非局域超表面”的直接前身：运算靠器件整体的角度响应完成，而不是靠逐点调制。
论文还指出，这类设计可以用”体积等效参数”（volumetric constitutive parameters）描述超屏的等效行为，为工程实现给出了具体路径 [S1]。
图 3 是路线二：一块横向均匀、纵向分层的多层平板，透射响应随入射角按目标传递函数变化，光穿过平板即完成运算。
图 3：论文原图——路线二（格林函数方案）：多层平板示意与导数/积分运算仿真。来源：论文图 3（Silva 等，Science 343(6167), 2014，仅作文献解读引用）。
表 1 对比了这两种架构。
表 1：论文提出的两种实现架构对比。来源：论文正文与后续综述（S1、S4、S7）。
| 对比项 | 路线一：傅里叶域处理器 | 路线二：格林函数方案 |
|---|---|---|
| 结构 | GRIN 波导 + 超屏 | 多层平板（横向均匀、纵向分层） |
| 运算位置 | 频域平面（先变换 → 滤波 → 逆变换） | 光穿过平板的传播过程中 |
| 系统长度 | 较长，需两级变换 | 无需透镜，最紧凑 |
| 算子灵活度 | 高，可按像素设计任意 H(kx) | 受结构对称性限制（偶函数类易实现） |
| 后续发展 | 空间光调制器/4f 系统思路 | 非局域超表面（基于角响应的运算器件） |
论文验证了什么：一整套”光子微积分”
这是一篇仿真研究：作者用全波数值仿真验证了概念 [S1][S2]，演示了三类运算：
- 二阶空间微分：设计传递函数 $H(k_x) \propto -k_x^2$，输入轮廓通过后输出其二阶导数——这正是图像边缘检测所需的拉普拉斯类运算；
- 积分：在频域按积分核设计，并对零频附近的奇点做截断与归一化处理，仿真输出与解析积分结果吻合；
- 卷积：设计卷积核为宽 16 µm（约 5.33 个波长，对应中红外工作波长）的矩形函数，演示了空间卷积运算 [S1]。
细心的读者会发现，这些运算的传递函数在频域都是偶函数（如 $-k_x^2$、积分核的截断形式）。为什么强调这一点？因为一阶微分对应的 $H(k_x) = ik_x$ 是奇函数——从负频到正频需要一个 π 相位跳变，这要求结构本身破缺镜像对称，或改用斜入射，在当时无法实现 [S4]。这个限制后来才被突破：先是利用斜入射和布儒斯特效应实现一阶微分 [S7]，再是 2018–2019 年的非局域介质超表面在近正入射下实验完成一阶/二阶空间微分与边缘检测 [S5]。
图 2、图 3 的仿真结果正是这三类运算的演示：输入轮廓穿过运算块后，输出与解析结果吻合 [S1]。
意义与局限
意义：这篇论文把”计算”从电子学带进了电磁与光学设计。它确立的范式——把目标算子翻译成传递函数、再用材料结构去实现——被后来几乎所有超表面模拟计算工作继承。无论傅里叶域滤波还是非局域超表面，都可以看作”传递函数工程”的不同实现 [S4]；而衍射神经网络则把这一思想推广为多层可训练的结构。论文提出的格林函数方案，直接催生了非局域超表面边缘检测这一整个方向 [S4][S5]。
局限（对论文本身的评价）：
- 只做了数值仿真，没有实物实验——2014 年时，把一块平板按指定角度响应精确造出来还很困难 [S1][S2]；
- 一阶微分（奇函数算子）在对称结构下无法实现，论文只能覆盖偶函数类运算 [S4]；
- 格林函数方案的多层平板参数依赖数值优化，早期优化结果偏离常规可实现范围，工程价值受限 [S7]；
- 作为概念验证，论文聚焦于”能不能算”，没有系统讨论带宽、效率与加工容差 [S4]。
这些局限不是论文的失败，而是它给领域留下的”待办清单”：实验验证、对称性破缺、参数可制造化、性能评估——后来十几年的文献几乎是一一对应补全的。
参考资料
- [S1] Silva, A., Monticone, F., Castaldi, G., Galdi, V., Alù, A., Engheta, N. “Performing Mathematical Operations with Metamaterials.” Science 343(6167): 160–163 (2014). DOI 页面
- [S2] Sihvola, A. “Enabling Optical Analog Computing with Metamaterials.” Science 343(6167): 144–145 (2014). DOI 页面
- [S3] Zangeneh-Nejad, F., Sounas, D. L., Alù, A., Fleury, R. “Analogue computing with metamaterials.” Nature Reviews Materials 6(3): 207–225 (2021). Nature 页面
- [S4] Liu, T., Qiu, J., Shi, X., Liu, Q., Xiao, S. “Flat optics for analog computing: from fundamental mechanisms to advanced meta-processors.” arXiv:2604.16849 (2026). arXiv 页面
- [S5] Cordaro, A., Kwon, H., Sounas, D., Koenderink, A. F., Alù, A., Polman, A. “High-Index Dielectric Metasurfaces Performing Mathematical Operations.” Nano Letters 19(12): 8418–8423 (2019). DOI 页面
- [S6] Zhou, H., Zhao, C., He, C., Huang, L., Man, T., Wan, Y. “Optical computing metasurfaces: applications and advances.” Nanophotonics 13(4): 419–441 (2024). DOI 页面
- [S7] Zangeneh-Nejad, F., Khavasi, A., Rejaei, B. “Analog optical computing by half-wavelength slabs.” Optics Communications 407: 338–343 (2018). DOI 页面