文献解读 | 超表面全光加速度计：偏振编码位移读出与惯性加速度测量（Science Advances 2026）

论文简介

本文解读的工作是《All-optical accelerometer enabled by metasurface》，作者为 Haofeng Zang、Yichen Wang、Tianchen Zhang、Yonghua Lu 与 Pei Wang（中国科学技术大学），2026 年 9 月 4 日发表于 Science Advances 第 12 卷，原文链接为 https://doi.org/10.1126/sciadv.aef8922 ，开放全文见 https://pmc.ncbi.nlm.nih.gov/articles/PMC13544244/ 。

一句话概括这篇论文的贡献：它把一块”空间上偏振取向连续旋转”的超表面当作位移标尺，用光纤双程回射加检偏的方式直接读出位移，再把这块超表面装到塑料悬臂梁上，做成了位移噪声底 340 pm/√Hz（1 Hz）、加速度噪声 900 ng/√Hz（1 Hz）、动态范围 85.6 dB 的惯性加速度计。

要理解这项工作为什么值得单独讨论，需要先分清两件事：超表面能做”位移传感器”已经不是新闻，把位移传感器的分辨率原样搬成”加速度分辨率”却不容易。加速度要先经机械结构转成位移，这个机械环节有带宽、有阻尼、有模式耦合，它的传递函数会限定加速度噪声、量程和频率响应。论文的看点是这条从加速度到光信号的完整链条被压缩在一个光纤端的探头里，而且每一步都有实测数字。

位移怎样变成偏振旋转

超表面编码的是”随位置变化的偏振取向”。论文的样品由非晶硅（α-Si）椭圆纳米柱排成阵列，每个柱子相当于一个局域的波片，柱长轴相对 x 轴的取向角 $\theta$ 沿 x 方向线性变化，

$$\theta(x) = \frac{\pi x}{\Lambda},$$

其中 $\Lambda$ 是结构的空间周期。线性偏振光入射后可以分解为等幅的右旋圆偏振（RCP，right-handed circularly polarized）与左旋圆偏振（LCP，left-handed circularly polarized）分量，经过这种取向渐变的纳米柱后，两个分量获得的几何相位（Pancharatnam–Berry 相位）符号相反，输出重新合成时表现为线偏振方向的整体旋转。

单次通过按教科书里半波片的结论就是把线偏振方向转 $2\theta$，那么论文强调的双程结构为什么不是互相抵消？关键在于回射。光束被平面镜反射后按坐标变换 $x’=x$、$y’=-y$、$z’=-z$ 原路返回，第二次遇到的是”有效取向角为 $-\theta$”的同一根纳米柱，而镜面反射本身又让圆偏振的手性翻转——教材里有一句很好用的判据：法向反射时右旋圆偏振会变成左旋圆偏振。两者共同作用的结果是旋转量叠加而不是抵消，对线偏振输入 $|E_{\mathrm{in}}\rangle=[1,0]^{\mathsf T}$，论文给出的输出为

$$|E_{\mathrm{out}}\rangle = \begin{bmatrix}\cos 4\theta \ \sin 4\theta\end{bmatrix} = \begin{bmatrix}\cos\frac{4\pi x}{\Lambda} \ \sin\frac{4\pi x}{\Lambda}\end{bmatrix}.$$

这个式子把”位移 → 偏振方向”的映射写得很干净：横向移动 $x$，输出线偏振方向就线性旋转 $\frac{4\pi x}{\Lambda}$。后面接一个固定方向的检偏器，按马吕斯定律，透射光强随夹角按 $\cos^2$ 变化，光强就变成了位移的电学可读信号。

图 1：全光超表面加速度计的结构与原理。（A）概念布局；（B）器件细节，由偏振编码超表面与两组用于回射的四分之一波片（Q）和反射镜（M）组成；（C）偏振编码超表面由硅椭圆纳米柱阵列构成，柱长轴取向角 $\theta$ 沿 x 方向以周期 $\Lambda$ 线性变化；（D）光正向与反向通过同一纳米柱时有效取向反转，旋转角由 $\theta$ 变为 $-\theta$。

这里有一个容易被略过的细节：$4\theta$ 意味着光强变化的周期比结构周期短得多。取论文的 $\Lambda=5.84$ μm，输出偏振方向每 $\Lambda/2 \approx 2.92$ μm 就旋转一整圈，单路检偏后光强周期只有 $\Lambda/4 \approx 1.46$ μm。也就是说，直接读一路光强的话，位移每走 1.46 μm 信号就重复一次，量程会非常小。论文用两条探测光路解决这个问题：两路信号相位差固定为 $\pi/2$，构成一对正交信号（下文会看到它们的李萨如图是一个圆），于是位移可以跨周期连续解算。分辨率和量程这对矛盾，正是靠这一步解耦的。

器件与读出：一块超表面加两条回射光路

器件尺度的参数决定了它的波长和工艺窗口。设计用 1550 nm（通信波段）作为工作波长，非晶硅在该波长的折射率为 3.288（椭偏仪实测），纳米柱高度 740 nm，晶格周期 730 nm；为了使单根柱子尽量接近理想半波片，作者扫参数库找到了椭圆长短轴直径 $D_o=590$ nm、$D_e=260$ nm 的组合。整块超表面的超胞尺寸为 5.84 μm × 730 nm，实验中使用的样品横向口径为 800 μm。

表 1：论文中位移读出模块的关键设计参数（数据取自论文正文与 Methods）。

参数	数值	说明
工作波长	1550 nm	通信波段，与光纤链路匹配
柱材料与折射率	非晶硅（α-Si），$n=3.288$	椭偏仪实测
纳米柱高度 $H$	740 nm	固定高度扫直径
晶格周期 $P$	730 nm	仿真域与周期边界条件同尺寸
半波片条件	$D_o=590$ nm，$D_e=260$ nm	长轴/短轴直径
超胞尺寸	5.84 μm × 730 nm	对应结构周期 $\Lambda=5.84$ μm
超表面口径	800 μm	实验样品横向尺寸

实验链路是”全光纤进出”的：1550 nm 线偏振激光耦合进单模光纤，经偏振控制器与光纤环行器送到探头；准直后照射超表面，回光沿原路返回，再由环行器输出，经 50:50 非偏振分束器送到两路检偏与光电二极管（PD，photodiode）。超表面装在一台压电台（Coremorrow P11.XYZ100S）上，可以纳米级精度横向移动，用来提供标准振动；数据以 10 kHz 采样。

图 2：全光超表面加速度计的表征。（A）基于光纤的实验装置，激光波长为 1550 nm，PC 为偏振态控制器，OC 为光纤耦合器，Q 为四分之一波片，M 为反射镜，A 为检偏器，PD 为光电二极管；（B）超表面的 SEM 图像，显示取向空间渐变的硅纳米柱阵列，周期 5.84 μm，标尺 2 μm；（C）归一化光电二极管电压随超表面横向位移呈余弦变化；（D）两路探测电压构成的李萨如图，两路相位差 $\pi/2$；（E）固定位移下以 10 kHz 采样记录 5 min 的时间迹；（F）测得噪声的幅度谱密度。

图 2C 的余弦曲线与图 2D 的圆是这套读出方案的”体检报告”：前者说明偏振旋转确实被转成了干净的余弦强度信号，后者说明两路正交输出可以互补地确定位移落在余弦的哪一段。图 2E、2F 则是分辨率证据——在固定位置上连续采集，幅度谱密度（ASD，amplitude spectral density）在 1 Hz 处最低约 340 pm/√Hz。注意单位：这是”每根号赫兹”的谱密度，不是单次测量的峰峰值误差；它意味着在 1 Hz 附近、1 Hz 带宽内可分辨约 0.34 nm 的位移，而信号在 2.5 nm 振幅下就已经被测到，与这个量级自洽。

从 45 微米到 2.5 纳米：振动测量

位移分辨率说清楚了，接下来看它能不能跟得上运动。论文用压电台给出标准正弦振动，先试了大振幅：5 Hz、45 μm 振幅。此时位移在一个信号周期内跨越了十几个光学周期（单路信号周期约 1.46 μm），图 3A 上半部分的电压迹线看上去已经”糊”成一团快速振荡，肉眼无法判断波形；但用双路正交信号解算出的位移曲线（图 3A 下半部分）仍是一条干净的正弦，对应的频谱（图 3B）只在 5 Hz 处出现单峰。

另一端是小振幅：5 Hz、2.5 nm 振幅的正弦同样被重建出来（图 3C、3D）。把 2.5 nm 与单路 1.46 μm 的信号周期放在一起看，动态结构的跨度接近三个数量级，这也解释了为什么必须做跨周期解调——只靠单路余弦信号，2.5 nm 这样的小位移虽然能分辨，但无法确定它落在哪个周期上。

图 3：标准振动信号的测量结果。（A）5 Hz、45 μm 振幅正弦运动下，两路归一化电压信号（上）与由两路信号重建的振动波形（下）；（B）对应的振动波形傅里叶谱；（C）5 Hz、2.5 nm 振幅正弦运动下的两路电压信号（上）与重建波形（下）；（D）对应的傅里叶谱；（E）0.01–100 Hz 多种驱动频率下测得的标准振动信号频谱。

频率覆盖是这套读出的另一个看点。论文把驱动频率从 0.01 Hz 扫到 100 Hz，测得频率与设定值吻合；其中最低的 0.01 Hz（一个周期 100 s）连续监测了近 5 小时，用于考察低频段的长期稳定性。需要说清楚的是，这一轮实验的上限由压电台的驱动带宽决定，而不是偏振读出机制本身：光学读出没有需要维持的干涉工作点，因此不受”偏离线性区”这类约束。

论文还进一步验证了波形的普适性：让超表面按 1 Hz 的正弦、方波、三角波、锯齿波运动，重建出的位移曲线保持了尖角特征，说明读数对非谐、非光滑运动同样有效。由位移再微分得到速度与加速度时，论文做了一步必要的处理——二阶微分会把远离基频的高频电子噪声放大，因此对频谱施加了 100 Hz 以下的低通滤波。这是所有”由位移推加速度”方案共同面对的信号处理问题，论文的处理方式与滤波器带宽都写得很明确。

从位移到加速度：悬臂梁与倾角标定

要把位移计变成加速度计，需要一个把惯性力转成位移的机械环节。论文用 ABS（acrylonitrile butadiene styrene，丙烯腈-丁二烯-苯乙烯共聚物）塑料 3D 打印了一段 S 形截面的悬臂梁，超表面固定在质量块上，敏感轴与悬臂梁的 x 轴一致。有限元分析（FEA，finite element analysis）给出 x 向平动的一阶固有频率为 24.2 Hz，明显低于 y、z 向平动与各转动模态的频率，这样在低频段测量的主要是敏感方向的运动，轴间串扰被压低。

标定用的是倾角法：把器件装到精密倾斜台上，敏感轴与安装面平行，从 −0.2° 扫到 +0.2°，每步 0.01°。重力沿敏感轴的分量正比于倾角的正弦，因此 ±0.2° 对应 ±3.5 mg（$g\cdot\sin 0.2° = 9.81 \times 3.49\times 10^{-3} = 3.49$ mg 的加速度变化）。用这条标定曲线拟合出位移到加速度的灵敏度

$$k_s = 4.74 \times 10^{-4}\ \mathrm{m/g},$$

即 1 个重力加速度对应约 474 μm 的位移——这个数字后面还会用到，它把”位移分辨率”和”加速度分辨率”联系起来。

分辨率与动态范围：数字成立的条件

装到悬臂梁上之后，论文报告了三组关键指标：1 Hz 处加速度噪声 900 ng/√Hz（ng 为 nano-g，即 $10^{-9}g$）、10 Hz 处噪声底 1 μg/√Hz、实测动态范围 85.6 dB。这些数字内部是自洽的，可以当场验算。

动态范围的定义是最大可测加速度与噪声底之比的分贝值。取实测上限 0.019 g 与 10 Hz 处噪声底 $1\times10^{-6}$ g，

$$20\log_{10}\frac{0.019}{1\times10^{-6}} = 20\log_{10}(1.9\times10^{4}) = 85.58\ \mathrm{dB},$$

与论文给出的 85.6 dB 一致，量纲上分子分母同为重力加速度单位，比值无量纲。论文进一步估算：若把量程推到超表面口径允许的极限，最大加速度约 1.93 g，理论动态范围为 $20\log_{10}(1.93\times10^{6}) \approx 125.7$ dB——这个估算的前提是”最大可测位移取决于超表面横向口径”，即量程由器件尺寸而不是读出原理限定。

位移与加速度之间的换算也应当对得上。用上文的灵敏度 $k_s$，1 μg 对应的位移是 $4.74\times10^{-4} \times 10^{-6}\ \mathrm{m} = 4.74\times10^{-10}\ \mathrm{m} = 474$ pm；反过来，900 ng/√Hz 折算成位移约 $900\times10^{-9} \times 4.74\times10^{-4} = 4.27\times10^{-10}\ \mathrm{m} = 427$ pm/√Hz。它与直接测得的 340 pm/√Hz 位移噪声底处于同一量级而不是同一数值，这符合预期：加速度噪声里包含了悬臂梁的机械噪声与环境扰动，位移噪声则是在压电台上单独测得的，两者不是同一环节的读数。

图 4：惯性加速度测量结果。（A）惯性超表面加速度计中 3D 打印 ABS 悬臂梁的结构示意；（B）把测得超表面位移映射到加速度信号的标定曲线，倾角增量为 0.01°；（C）显示谐振峰的加速度计动态响应；（D）静止状态下记录的加速度噪声；（E）超表面加速度计的加速度幅度谱密度；（F）10 Hz 不同振幅激励下重建的加速度信号。

图 4C 至 4E 交代了这些指标的适用条件。响应曲线里的谐振峰位置对应悬臂梁一阶模态，而在远高于谐振的频率上，响应并没有按理论预期那样快速衰减，原因是安装后压电台内部机械结构引入了一个约 85 Hz 的寄生谐振，把高频段抬了起来。图 4D、4E 是静止状态下的加速度噪声与其谱密度，也就是 1 μg/√Hz 这个数字的直接来源。图 4F 是 10 Hz 不同振幅激励下的重建信号，跨度内线性度 $R^2 = 0.998$，说明在实测范围内响应是线性的。

论文对噪声来源的交代值得单独一提：作者明确说明当前噪声底主要反映环境扰动，而不是通过交叉相关等算法提取出的器件本征自噪声，实验也未做真空或温控隔离。因此 340 pm/√Hz 与 900 ng/√Hz 应理解为”本实验条件下”的系统指标；这既说明这套架构的潜力尚未被完全释放，也提示横向对比不同论文的加速度噪声数字时，要先看清测量环境与噪声分离方法是否可比。

意义与展望

这篇工作的价值不在于把某一项指标推到极限，而在于把”超表面位移计量”和”惯性传感”之间缺的那一环补上了：超表面提供亚纳米级、跨周期的位移标尺，悬臂梁提供加速度到位移的机械转换，光纤把两者连成一个可以远程放置的探头。整套方案不需要维持干涉工作点，也不依赖成像光路，读数来自偏振基的选择，这让量程与分辨率可以在同一器件上兼顾。

从原理验证走向实际应用，论文给出了几条明确的推进方向。机械侧可以更换蠕变与滞回更小的材料，并优化悬臂几何与阻尼设计；环境侧可以引入真空与温控封装，把现在占主导的环境噪声压下去，同时拓展低频响应；光学侧则可以借用矩阵超表面实现多维位移感知、用 Stokes 参数超表面直接在光纤端面完成偏振探测，从而做出矢量加速度计。由于超表面与硅 MEMS 工艺兼容，把这些环节集成到同一芯片上的路径也是清楚的。对地震监测、结构健康监测和导航这类既要求低噪声、又要求低功耗与小体积的场景来说，这条技术路线提供了一个值得持续跟踪的选择。

参考资料

1. Zang H., Wang Y., Zhang T., Lu Y., Wang P. All-optical accelerometer enabled by metasurface. Science Advances 12 (2026). https://doi.org/10.1126/sciadv.aef8922 （开放全文：https://pmc.ncbi.nlm.nih.gov/articles/PMC13544244/ ）
2. Reschovsky B. J., Long D. A., Zhou F., Bao Y., Allen R. A., LeBrun T. W., Gorman J. J. Intrinsically accurate sensing with an optomechanical accelerometer. Optics Express 30 (2022). https://doi.org/10.1364/oe.457499
3. Meng Y., Hu Q., Vatankhah E., Liu Z., Wang W., Hall N. A. An integrated monolithic photonic accelerometer with on-chip fringe-counting. Microsystems & Nanoengineering 12 (2026). https://doi.org/10.1038/s41378-026-01410-1
4. Peatross J., Ware M. Physics of Light and Optics（BYU，作者免费分发）：第 6 章「Polarization of Light」§6.4 线性偏振片与琼斯矩阵、§6.6 波片琼斯矩阵与法向反射的手性反转. https://optics.byu.edu/