文献解读 | 逆向设计手性超表面：把螺旋二色性从 115% 仿真做到 62% 实测

论文简介

这篇论文是 Inverse-Designed Chiral Metasurfaces for Enhanced Helical Dichroism，2026 年 9 月 23 日上传到 arXiv（编号 2609.28000，原文链接 https://arxiv.org/abs/2609.28000 ），全文 17 页、4 张图。作者是 Chia-Chun Pan、Yu-Hsi Chen、Munseong Bae、Haejun Chung、Dan Smith、Daniel Fan、Sejeong Kim 与 Ranjith R Unnithan，单位包括墨尔本大学电子电气工程系、汉阳大学电子工程系与成均馆大学电子与计算机工程系。

一句话概括这项工作：作者把拉盖尔-高斯（Laguerre-Gaussian, LG）光束的空间形态直接写进逆向设计的目标函数，让拓扑优化自己去找结构，最终得到一片 6 μm 见方的氮化硅（Si₃N₄）手性超表面，在 ℓ = ±2 两种轨道角动量（orbital angular momentum, OAM）模式下给出约 115% 的模拟螺旋二色性（helical dichroism, HD），实物测量在 800 nm 附近达到 52.6%，峰值 61.7%。

需要说明的是，这是预印本，尚未经过同行评审。文中所有数值都来自作者自算与自测，目前还没有看到独立的重复实验。

螺旋二色性：把探针从自旋换成轨道

先分清楚两类”手性光学”信号，否则很容易把 HD 和常见的圆二色性（circular dichroism, CD）看成同一件事。

圆二色性是大家更熟悉的那一类：手性分子对左旋和右旋圆偏振光的吸收不同，用吸收差或透射差就能读出分子的手性。圆偏振光携带的是自旋角动量（spin angular momentum, SAM），通俗地说，电矢量在传播过程中原地打转。LG 光束提供的则是另一套自由度：光束横截面上的相位绕光轴转一圈，累积 $2\pi\ell$ 的相位，整数 $\ell$ 就是拓扑荷，正负号对应相位旋向相反的两种涡旋光。它携带的是轨道角动量，和圆偏振的自旋角动量彼此独立。

HD 测量的就是这两种涡旋光被同一个结构”区别对待”的程度。论文采用的计量方式是

$$ \mathrm{HD} = 2\times\frac{R_{\ell=-2}-R_{\ell=+2}}{R_{\ell=-2}+R_{\ell=+2}}\times 100% $$

其中 $R$ 是两种入射 OAM 模式各自的模态反射率，由反射监控面上坡印廷矢量对入射功率归一化后积分得到。写成对比度的形式，是为了让数值不依赖绝对反射率：只要正负两个模式的反射差别足够大，HD 就能取到很高的值，甚至超过 100%。这一点后面还会回到。

这里有一个容易忽略的前提。如果结构本身是几何对称的，那么它对 $\ell = -2$ 和 $\ell = +2$ 的响应几乎是镜像等价的，分子上的差自然接近零。所以高 HD 必须同时具备两件事：一片真正具备手性的结构，以及一束把相位结构”带进来”的涡旋光。论文专门做了对照实验来钉住这个前提，用一个几何对称的图形做非手性参照，测出来两个 OAM 态的反射谱几乎重合。

表 1：圆二色性与螺旋二色性的对比。

对比维度	圆二色性（CD）	螺旋二色性（HD）
探测手段	左/右旋圆偏振	正/负拓扑荷涡旋光
携带的角动量	自旋角动量（SAM）	轨道角动量（OAM）
对结构的要求	结构具备手性	结构具备手性，且光束相位结构与结构尺度匹配
本文指标	论文未涉及	模拟约 115%，实测峰值 61.7%

在全介质超表面里，CD 这一路已经相当成熟：二维手性超表面能够在较宽的波段上同时实现高 CD。HD 的路线起步更晚，但可以引入更多可调变量，比如涡旋光的拓扑荷、光束束腰与结构的相对尺寸，也更容易在设计阶段把这些空间参数一起优化。

设计：让梯度在 6 μm 见方的区域里找手性

设计问题本身很直白：在给定材料与版图约束下，让 $\ell=-2$ 的反射尽量高、$\ell=+2$ 的反射尽量低。麻烦在于，能满足这个要求的几何形状实在太多，人靠直觉画出来的螺旋、L 形或者缺口结构只是其中一小部分。

论文的做法是把结构”溶掉”，用伴随法（adjoint method）拓扑优化来搜索。设计区域里每个位置的介电常数在空气（$n=1$）与 Si₃N₄（$n=2.024$）之间连续取值，初始给一个带手性的种子形状。每轮迭代要做两次电磁仿真：一次是入射涡旋光激发的正向场，一次是由目标函数反向传播得到的伴随场，两者的场分布重叠部分给出梯度，再由 Adam 优化器更新材料分布。随着迭代推进，材料分布被逐步二值化，最终收敛为可以直接送去光刻的二元版图。整个过程用开源 FDTD 软件 MEEP 完成，工作波长 800 nm，空间分辨率 20 像素/μm。

具体参数值得记下来，因为它们决定了这件器件的尺寸感：计算域 9.2×9.2×6.6 μm³，四周是 0.8 μm 厚的完美匹配层；真正可自由变化的有效设计区只有 6×6×1.2 μm³，结构高度固定 1.2 μm，下面接半无限厚、折射率 1.45 的二氧化硅衬底。为了保证刻得出来，优化过程中强制最小特征尺寸和最小间隙都不小于 0.2 μm。这几个数字放在一起看，就知道所谓”超表面”在这里其实是平面工艺里一块 6 μm 见方、1.2 μm 高的小凸台。

图 1：逆向设计手性超表面的演化过程与其手性光学响应。(a) 优化轨迹：结构从连续材料分布逐步收敛为最终二值版图（上排），下方为对应的反射率与螺旋二色性（HD）表现，背景阴影表示施加的二值化程度，最终结构在 800 nm 处取得约 115% 的 HD。(b) 优化结构在轨道角动量模式 ℓ = −2（左）与 ℓ = +2（右）照射下的近场电场强度与相位分布。

图 1(a) 里的曲线值得多看两眼。随着二值化程度不断提高，HD 一度出现回落：把连续材料分布强行推向”非黑即白”的过程会扰动散射相位，目标函数跟着抖动，这是拓扑优化里常见的一幕。优化器最终仍然收敛到全二值版图，并在这个版图上给出约 115% 的模拟 HD。图 1(b) 给出的近场强度与相位分布解释了原因：同样的结构，$\ell=-2$ 入射时反射场被强烈激发，$\ell=+2$ 时则被压制，两个模式在结构内部的场分布走的是完全不同的路。

顺带把 115% 这个数理解清楚，它并不代表效率。按定义，设 $R_{\ell=+2}=1$（只看比值），$R_{\ell=-2}=R$，令 $\mathrm{HD}=1.15$ 有 $2(R-1)/(R+1)=1.15$，解出 $R\approx 3.7$。也就是说，模拟结果对应的是 $\ell=-2$ 的反射率约为 $\ell=+2$ 的 3.7 倍，两个模式的反射差别很大，但单看某个模式的绝对反射效率，仍然由场的归一化结果决定，论文并未宣称超过 100% 的效率。

作者选择 $\ell=\pm 2$ 也不是拍脑袋。他们对不同 OAM 阶数做过一轮扫描，$\ell=\pm 2$ 给出的 HD 最高，于是把自由度和优化目标都锁定在这个模式上。此前文献里，基于伴随法的逆向设计在 $|\ell|=3$、800 nm 条件下报告过约 107% 的模拟 HD，参数化共振几何的变形六角硅结构报告过约 108%，这篇文章的 115% 大体处在同一量级，真正往前推了一步的是后面的实物验证。

制备与测量：EBL、SLM 与 1.13 μm 的束腰

仿真里的连续灰度分布没法直接刻出来，必须落到工艺上。样品做在石英衬底上，用电子束光刻（electron-beam lithography, EBL）定义图形，蒸发一层镍作为刻蚀掩模，再用电感耦合等离子体反应离子刻蚀把图案转移到 Si₃N₄ 层，最后用 Piranha 溶液去掉残余镍。为了便于逐个测量，结构做成了阵列，相邻单元间距保持在 20 μm 以上，把近场耦合的影响压下去，同时也留出加工冗余。

图 2：制备所得手性超表面的形貌与光学表征。(a) 超表面阵列的扫描电子显微镜（SEM）图像，下方标注对应的阵列编号。(b) 单个手性结构的放大 SEM 图像。(c) 白光照明下由彩色 CMOS 相机拍摄的阵列光学显微图像，插图为 ℓ = −2 的 LG 光束聚焦在阵列旁裸衬底上的焦斑轮廓。(d) 聚焦 LG 光束与所选手性超表面对准的光学显微图像，插图为关闭白光后放大的光束—结构重叠区域。

图 2 里最该注意的是 (b) 的放大 SEM 图。它显示设计的手性几何确实被复现出来了，但侧壁有倾斜、边缘有粗糙、转角有圆化。这几处看起来不起眼的偏差，正是后面实验值低于模拟值的主要原因，我们放到结果一节再展开。

测量端要解决的是”如何可控地把 OAM 加载到光束上，并且知道光打在结构的什么位置”。作者用超连续激光经波长滤波器取出中心 800 nm、带宽 10 nm 的激发光，偏振由半波片加线偏振片固定。LG 光束不是靠螺旋相位板或叉形光栅做出来的，而是把计算好的 LG 相位与一个倾斜平面波的相位叠加成闪耀光栅全息图，加载到液晶硅基空间光调制器（Santec SLM-200）上；全息图引入 0.5° 的衍射角，把需要的一级衍射与零级分量在空间上分开，避免零级直接漏进探测器。

聚焦条件也需要算一笔账。系统用的是 50× 物镜（NA = 0.55，焦距 4.0 mm），但入射光束直径只有约 1.8 mm，比物镜的通光孔径小，物镜实际处于欠填充状态，有效数值孔径只有 $D/(2f)\approx 0.225$。对应的束腰按 $w_0=\lambda/(\pi\mathrm{NA}_{\mathrm{eff}})$ 估算，在 800 nm 处约为 1.13 μm。

这个数字和结构尺寸的关系很关键：手性结构本身约 6 μm 见方，而聚焦后的光束束腰只有 1 μm 量级，也就是说光束照亮的是结构的一部分，而不是整片均匀覆盖。论文把这一点当作模式匹配问题来处理，先把聚焦光斑对准阵列旁的裸衬底确认焦斑轮廓，再平移样品台把光束移到目标结构上，并在白光关闭的条件下用放大的重叠图像确认对准。反射信号经同一个物镜收集，耦合进多模光纤送到光谱仪，测量时关掉所有室内与卤素光源、扣除背景，并以金参考谱归一化。

实验结果：52.6% 与 61.7%

先说对照。作者挑了一个相邻 EBL 标识上几何对称的图形，在完全相同的条件下测量。两个 OAM 态的反射谱几乎重合，导出的 HD 谱在大部分测量范围内贴着零走，只在设计波长附近有小幅起伏。这一步没有新物理，但把”观测到的差异到底来自结构手性，还是来自光路漂移或背景噪声”这个疑问排除掉了，后面测到的高 HD 才有意义。

再看手性结构。把聚焦的 LG 光束对准单个手性超表面之后，$\ell=-2$ 与 $\ell=+2$ 的反射谱出现了清楚的分离，$\ell=-2$ 在 800 nm 附近反射更强。由此算出的 HD 在 800 nm 处为 52.6%，在 801.6 nm 达到 61.7% 的峰值，摘要中表述为”最高约 62%”。

图 3：螺旋二色性（HD）的实验测量结果。(a) 相邻 EBL 标识中一个几何对称结构（作为非手性参照）在 ℓ = −2（绿）与 ℓ = +2（黄）LG 光束照射下测得的反射光谱。(b) 由 (a) 得到的对应 HD 谱，非手性参照的手性光学响应几乎为零。(c) 同一条件下单个手性超表面测得的反射光谱。(d) 手性超表面相应的 HD 谱：在 800 nm 处为 52.6%，并在 801.6 nm 达到 61.7% 的最大值。

同样做一次代入，把 61.7% 翻译成反射率的比值：令 $R_{\ell=+2}=1$、$R_{\ell=-2}=R$，$2(R-1)/(R+1)=0.617$ 解得 $R\approx 1.9$，也就是 $\ell=-2$ 的反射约为 $\ell=+2$ 的 1.9 倍。相比模拟里 3.7 倍的差距，实测的对比度明显缩水，但方向与量级关系是对的：结构确实更”喜欢”其中一个涡旋态。

缩水的原因，论文归到几何差异上。SEM 图里的侧壁倾斜、边缘粗糙与转角圆化会改变局部的电磁场分布和共振条件，而逆向设计给出的版图是按理想矩形边界算出来的。这类误差在亚波长结构上会被放大：0.2 μm 量级的最小特征与间隙，本身就处在曝光与刻蚀误差的量级上。加上模式匹配对光度测量很敏感，束腰约 1.13 μm 的光要稳稳落在 6 μm 结构的目标位置上，对准偏差也会直接吃掉落差。

适用边界与待推进方向

关于这项结果能推多远，有几条边界需要交代。

它的验证发生在单波长附近。器件为 800 nm 设计，实验用带宽 10 nm 的光源验证了这个工作点，论文没有给出完整的宽带 HD 谱，因此把”高 HD”推广到宽带器件上还缺数据。

它的验证对象是单个结构，而不是阵列器件。测量到的信号来自一片 6 μm 见方的手性单元，阵列主要是为了测量冗余与定位方便。要变成一个能成像、能大面积读出的器件，需要处理阵列之间的一致性、串扰与相位关系。

论文没有报告绝对反射率与完整带宽的对应关系，所以”HD 高”与”总反射高”是两件事。HD 作为对比度指标，可以在总反射不高的情况下取到漂亮的数值；评价器件时要把两个指标放在一起看。

最后，这是预印本结果，还没有第三方复现。文中给出的 115% 模拟值与 62% 实测值都建立在作者自己的仿真设置与测量系统上，独立验证出现之前，把它们当作”该设计路线的可行性证据”是合适的，当作成熟器件的性能指标则偏早。

从这些边界往外看，真正待推进的方向也不难列：把工艺扰动纳入优化循环，让设计对侧壁倾角、边缘粗糙更宽容；把带宽写进目标函数，换取更宽的工作窗口；把光束束腰、对准误差这类系统参数与结构一起优化，减小模式失配带来的损失；以及用阵列化把单结构的响应扩展成可读出的图像信号。

意义与展望

螺旋二色性并不新，把 HD 幅度做大的努力也持续了多年，宏观器件层面甚至已经出现”giant HD”的说法。这篇文章的位置在于，它把携带 OAM 的光束的空间形态当成设计变量的一部分，交给逆向设计去处理，然后用 EBL 制备的 Si₃N₄ 结构在实验上验证了这条路走得通。这项工作把两件事拆开处理：几何形状交给梯度优化决定，手性光学响应的判据由涡旋光提供。

这类工作最终指向的应用方向是手性传感与手性分辨测量。传统 CD 测量依赖圆偏振，而 OAM 提供了一组正交且可连续调节的模式，理论上能给同一份样品更多维度的读数。当前 62% 的实测 HD 是在单一结构、单波长附近取得的，距离实用器件还有工程距离，但”设计结构 + 设计光场”的组合已经显出了它的空间。

参考资料

1. Pan C-C, Chen Y-H, Bae M, Chung H, Smith D, Fan D, Kim S, Unnithan R R. Inverse-Designed Chiral Metasurfaces for Enhanced Helical Dichroism. arXiv:2609.28000, 2026. https://arxiv.org/abs/2609.28000
2. Jin Y, Wang X, Xia Z, et al. Multidimensional helical dichroism from a chiral molecular nanoassembly. Nature Communications, 2026, 17. https://doi.org/10.1038/s41467-026-68540-y
3. Helfrich C, Frosz M, Tani F. Giant Helical Dichroism in Twisted Hollow-Core Photonic Crystal Fibers. ACS Photonics, 2025, 12. https://pmc.ncbi.nlm.nih.gov/articles/PMC11844252/
4. Wang R, Wang C, Sun T, Hu X, Wang C. Simultaneous broadband and high circular dichroism with two-dimensional all-dielectric chiral metasurface. Nanophotonics, 2023, 12. https://doi.org/10.1515/nanoph-2023-0407
5. OpenStax, Rice University. College Physics 2e, §27.8 Polarization. https://openstax.org/books/college-physics-2e/pages/27-8-polarization
6. Tang Y, Zhao S, Luo L, Tang T, Xu H, Li J. Spin-multiplexed terahertz image edge detection based on extrinsic chirality in all-dielectric metasurface. arXiv:2609.25147, 2026. https://arxiv.org/abs/2609.25147