文献解读 | 超表面直接充当 BB84 测量器件：单孔径交织线偏振与圆偏振四个探测端口（arXiv 2026）

论文简介

题录：Yan He, Adetunmise C. Dada, Secure Quantum Metasurfaces for Direct Quantum Key Distribution Measurements, arXiv:2609.22946（2026-09-19 提交，physics.optics / quant-ph）。原文链接：https://arxiv.org/abs/2609.22946

作者机构：英国格拉斯哥大学物理与天文学院（School of Physics and Astronomy, University of Glasgow）。

一句话贡献：作者把量子密钥分发（Quantum Key Distribution，QKD）中最常用的 BB84 协议所需的两套偏振基——线偏振基 H/V 与圆偏振基 R/L——以小于半个波长的周期交织在同一块介质超表面孔径里，让这块超表面直接充当接收机的测量器件，把完整测量映射到四个探测端口，在 780 nm、1550 nm、2000 nm 与 10.6 μm 四个波长上分别给出 97.96%–99.08% 的条件概率保真度。

这篇文章值得读的地方，不在于”又做了一个超表面器件”，而在于它把两件事接了起来：一边是器件级的光学量（透过率、焦点分布、探测器区域怎么划），另一边是协议级的密钥率（每个入射光子最终能留下多少秘密比特）。论文的结论之一正是这两者不能互相替代——只看光斑干不干净会选错设计。

为什么 BB84 接收机可以只用一块超表面

BB84 协议（Bennett–Brassard 1984）的思路很朴素：发送方 Alice 每次从四个偏振态里挑一个发出去——水平（H）与竖直（V）线偏振构成一套基，右旋（R）与左旋（L）圆偏振构成另一套共轭基；接收方 Bob 每次随机选一套基去测。两人事后公开比对基的选择，只保留基选对的那部分比特，这个过程叫筛选（sifting）。基选错的那些比特必须丢掉，因为对不兼容的基测量，结果本身是随机的，没有信息价值。

问题出在”怎么测”。H/V 基可以用偏振分束器直接分开，R/L 基则要先经过一片四分之一波片把圆偏振转成线偏振，再进偏振分束器——或者反过来。也就是说，一台常规 QKD 接收机里至少要有分束器、波片、偏振分束器、光纤耦合与探测器这几段级联的体光学（bulk optics），体积和装调复杂度都摆在那里。想把接收机做小，最自然的想法是：能不能让一块平面光学器件同时完成两套基的投影？

超表面正好适合做这件事。它由亚波长尺度的纳米结构阵列构成，每个”元原子”可以独立控制透射光的相位与幅度，因此能同时完成偏振分解与聚焦。论文的做法是把这件事推到极致：不做”先分偏振、再各自测量”的两段式结构，而是让整块孔径本身就输出四路聚焦光斑，四个光斑分别落在四个探测器上，直接给出 BB84 的四个测量结果。论文把这种器件叫”直接测量”的接收机——输出的已经是一个四符号的协议结果，不需要再反演斯托克斯参量去推断”这一下按的是哪个键”。

从测量算符的角度看，理想的无源 BB84 接收机对应四个算符 $\hat E_H,\hat E_V,\hat E_R,\hat E_L$，每个都是对应偏振态的投影算符乘以 $1/2$：

$$
\hat{E}_H=\tfrac12\left|H\right\rangle\left\langle H\right|,\quad
\hat{E}_V=\tfrac12\left|V\right\rangle\left\langle V\right|,\quad
\hat{E}_R=\tfrac12\left|R\right\rangle\left\langle R\right|,\quad
\hat{E}_L=\tfrac12\left|L\right\rangle\left\langle L\right|.
$$

这里的 $1/2$ 是基的选择权重，不是额外的插入损耗——四个算符之和仍是单位算符。把它写成”输入态 × 探测结果”的条件概率矩阵，理想情形就是下面这个四阶矩阵：对角线上的匹配端口给出 $1/2$，不兼容基的两个端口各分走 $1/4$。

$$
C^{\rm ideal}=
\begin{pmatrix}
1/2 & 0 & 1/4 & 1/4\
0 & 1/2 & 1/4 & 1/4\
1/4 & 1/4 & 1/2 & 0\
1/4 & 1/4 & 0 & 1/2
\end{pmatrix}.
$$

后面所有”这块超表面做得好不好”的讨论，本质上都是在问：实测或仿真的条件概率矩阵离这个目标矩阵有多远，以及为此损失了多少光。

交织孔径：把两套基的相位库放进同一个棋盘

设计流程从单元库开始。作者在周期边界下扫描矩形介质纳米条的长度与宽度，记录线偏振输入的透射幅度与相位，把相位折叠到 $0$–$2\pi$ 区间，每个相位库只保留 8 个离散相位级。线偏振子晶格直接用一对尺寸分别逼近指定的 $x$、$y$ 相位；圆偏振子晶格则先筛出相位延迟约等于 $\pi$ 的纳米条，再靠它们的传播相位与面内旋转角分别提供右旋、左旋输入所需的平均相位与差分相位。

有了两套库，接下来给四路输出各分配一个聚焦相位。第 $m$ 个焦点的目标相位取理想无像差形式

$$
\phi_m(x,y)=-\frac{2\pi}{\lambda}\left[\sqrt{(x-x_m)^2+(y-y_m)^2+f^2}-f\right],
$$

其中 $(x_m,y_m)$ 是第 $m$ 个焦点的横向位置，$f$ 是焦距。这个表达式而不是常见的抛物近似，是因为大数值孔径下高阶相位项不能丢。相位量化到可用的库之后，线偏振单元与圆偏振单元被安排在同一个棋盘孔径的互补宇称格点上。

关键的尺寸关系在这里：四个选定设计的周期与波长之比 $p/\lambda$ 分别是 0.487、0.484、0.450 与 0.443，都小于 1/2；而任何一套宇称子晶格自身的重复距离是 $\sqrt{2}p$，对应 $\sqrt{2}\times380/780=0.689$、0.684、0.636 与 0.627，仍然小于一个波长。也就是说，两套基在空间上的分工尺度比它们各自作为周期阵列的尺度还要小。

图 1：BB84 测量的焦平面重构。(a) 交织的线偏振子晶格与圆偏振子晶格。(b) 选定 1550 nm 器件的 H 输入远场，叠加经标定固定的四个区域积分边界。(c) 由积分后的探测器功率得到的条件响应，概率保真度 99.08%。(d) 从单元扫描到 BB84 指标的分析流程：每个端口中心由匹配态光斑的「椭圆高斯 + 常数基线」局部拟合得到；同一器件使用同一个正方形，取保留每个拟合核心至少 99%、扣背景能量至少 95%、拟合基线贡献不超过 10% 的最小尺寸；随后对原始强度积分、把四个功率按行归一化，并用四个输入态的完整响应拟合出 H、V、R、L 的物理端口分配。

图 1(d) 里的标定步骤值得多说一句，因为它决定了后面所有数字是否可比。模拟得到的焦平面是一张 $200\times200$ 的强度图，四个光斑要划成四个积分区域。如果每个端口各划一个大小不同的框，响应矩阵就会被”哪个框大哪个数大”污染。作者的做法是：先用四块 $41\times41$ 像素的窗口做粗初始化，再用椭圆高斯加常数基线拟合每个匹配态光斑的真实中心；然后强制同一器件的四个端口使用同一个正方形边长，取满足三条判据的最小尺寸——至少保留 99% 的拟合核心、至少保留宽参考框内 95% 的扣背景能量、常数基线占比不超过 10%。最终边长是 23、39、29、21 像素（对应 780 nm、1550 nm、2000 nm、10.6 μm），最大拟合基线占比不超过 4.9%。最后还要遍历 24 种”物理端口 → 逻辑结果”的映射，取响应矩阵均方根误差最小的那一种。

四个波长：仿真远场是否真的分成四个焦点

论文给出 780 nm 到 10.6 μm 四个独立设计的全波（FDTD，Finite-Difference Time-Domain，时域有限差分）仿真结果。选定几何参数如下表，材料都是硅（Si），只有 780 nm 另有一个二氧化钛（TiO$_2$）候选设计用于对比。

表 1：四个选定波长对应的设计参数，$f$ 为焦距，”焦点偏移”给出四个焦点相对中心的横向位置。

波长	材料	单元周期	结构高度	焦距 $f$	焦点偏移
780 nm	Si	380 nm	450 nm	20 μm	±4.5 μm
1550 nm	Si	750 nm	1000 nm	100 μm	±10 μm
2000 nm	Si	900 nm	1100 nm	100 μm	±12.5 μm
10.6 μm	Si	4700 nm	10 000 nm	150 μm	±60 μm

在匹配基输入下，每一个设计的焦平面都出现”一个端口占优、同基正交端口被抑制”的图样，同时功率按协议要求在不兼容基的两个端口之间分配。四者的差异也很清楚：780 nm 的同基泄漏最大；1550 nm 的场最干净；2000 nm 与 10.6 μm 保持了很低的同基泄漏，但两个不兼容基端口之间的分配更不均匀。

图 2：选定四端口设计的全波焦平面强度。行为分别在 780、1550、2000 与 10 600 nm 优化的独立超表面，列为 H、V、R、L 四种输入。为便于空间对比，各面板的强度在面板内单独归一化，绝对透过率与收集效率另行定量给出。矩形标出定量分析中使用的、经标定固定的探测器区域。

这里有一个容易被忽略的判读细节：单看”光斑是不是够亮、够干净”并不足以说明设计成功。协议要求的是一个完整的概率分布——四乘四共十六个输入-结果概率都要对。图 2 中同时表达了两件事：匹配端要占优（否则误码升高），不兼容基的两个端口又要接近平分（否则基之间的相关性出现偏置）。这两条约束会拉扯设计自由度，后面的密钥率正是在这个拉扯里找折中点。

响应矩阵与内禀 QBER

把四张远场强度图按标定好的区域积分，就得到条件概率矩阵 $C_{i,m}=P(m\mid i)$：$i$ 是送入超表面的偏振态，$m$ 是探测器结果 H、V、R、L。以最干净的 1550 nm 设计为例，H 与 V 的匹配端口概率分别是 0.531 与 0.504，而同基正交端口的泄漏只有 0.005 与 0.008；R 与 L 的匹配端口概率为 0.467 与 0.465，正交泄漏均为 0.010。也就是说，被”按错键”的匹配基事件被压到了 1% 以下。

器件本身引入的误码用”按基条件的量子比特误码率”（Quantum Bit Error Rate，QBER）来量化。设 $N_{HV}$、$N_{RL}$ 分别是 H/V 基与 R/L 基内被接受的概率总和，则

$$
Q_{HV}=\frac{C_{H,V}+C_{V,H}}{N_{HV}},\qquad
Q_{RL}=\frac{C_{R,L}+C_{L,R}}{N_{RL}},
$$

平均内禀 QBER 取两者的算术平均。四个设计的结果是：780 nm 为 3.96%，1550 nm 为 1.68%，2 μm 为 2.01%，10.6 μm 为 1.42%。论文把这一项称为”内禀”，是因为它只包含器件自身造成的错误，不含信道衰减、探测器暗计数与光纤传输等环节。

图 3：条件四结果响应矩阵。(a) 理想无源 BB84 响应。(b–e) 由选定全波强度分别在 780、1550、2000 与 10 600 nm 计算得到的响应，概率保真度分别为 0.980、0.991、0.983 与 0.989。行为输入态，列为探测器结果，顺序为 H、V、R、L，每行在四个被接受区域内加和为 1。对角元是匹配结果；同一基内的正交端口代表比特误码；共轭基的两个结果构成理想情形下的均匀随机分配。

衡量矩阵与目标差多远，论文用的是平方 Bhattacharyya 系数

$$
F_i=\left(\sum_m\sqrt{C_{i,m}C^{\rm ideal}_{i,m}}\right)^{2},
$$

对四个输入态取平均得到平均概率保真度 $\overline F$。注意这个指标只看”归一化之后的分布形状”，不看绝对功率——这也是下一节要单独处理透过率的原因。

从光学指标到密钥率

密钥率计算是全文最有工程味的部分。对等基概率的 BB84，一个无损无错的接收机每次入射光子最多贡献 $Y_{\rm ideal}=1/2$ 比特：一半光子因基匹配被保留，每个保留光子带来 1 比特。由此，如果收集效率是 $\eta$，无错时的密钥量就是 $\eta/2$。

有错时就要扣两项。纠错（error correction）要在公开信道上交流信息，这一部分泄漏必须从密钥长度里减去；保密放大（privacy amplification）再压缩密钥，以消除窃听者可能掌握的信息。论文采用的光学基准保密分数是

$$
r_{\rm sec}=\max\left{0, 1-f_{\rm EC}h_2(Q_{HV})-h_2(Q_{RL})\right},
$$

其中 $h_2$ 是二元熵函数，纠错低效因子 $f_{\rm EC}=1.16$ 表示实际纠错公开的信息比理论最小值多 16%；对基 $Q_{RL}$ 的熵扣除则用反基误码率作为相位误码的替代估计。设备级密钥量再乘上平均筛选探测概率 $\overline p_{\rm sift}$：

$$
Y_{\rm sec}=\overline p_{\rm sift}r_{\rm sec}.
$$

这里可以做一次可以手算复核的数值代入。取 1550 nm 设计的 $Q_{HV}=0.0129$ 与 $Q_{RL}=0.0208$：$h_2(0.0129)\approx0.0995$，乘以 1.16 得 0.1154；$h_2(0.0208)\approx0.1459$。两者相减，$r_{\rm sec}\approx1-0.1154-0.1459=0.7387$，与论文给出的 0.73884 一致（差值来自熵函数的舍入）。再看筛选概率：该设计的四区域收集效率 $\eta^{\rm coll}=40.94%$，按 $\eta/2$ 估算应得 $0.4094/2=0.2047$，论文给出的 $\overline p_{\rm sift}=0.20472$，两者吻合，说明”基筛选取走一半”这个直觉在该模型下是成立的。两者相乘 $0.20472\times0.73884\approx0.1513$，正是论文报告的 $Y_{\rm sec}=0.15126$ 比特/入射光子。

表 2：四个选定设计的光学与协议指标。误码按基给出；$\overline F$ 为平均概率保真度，$\lvert T\rvert$ 为平均透过率幅度，$f^{\rm ROI}$ 为四区域收集分数（ROI，Region of Interest，感兴趣区域），$\eta^{\rm coll}=\lvert T\rvert f^{\rm ROI}$ 为收集效率，$r_{\rm sec}$ 为保密分数，$Y_{\rm sec}$ 为每次入射光子的密钥量；保密分数与密钥量为渐近模型下、纠错低效因子 1.16 的设备级数值，不含信道、光源、探测器与有限密钥效应。

波长	$Q_{HV}$	$Q_{RL}$	$\overline F$	$\lvert T\rvert$	$f^{\rm ROI}$	$\eta^{\rm coll}$	$r_{\rm sec}$	$Y_{\rm sec}$
780 nm	5.97%	1.96%	97.96%	37.55%	60.78%	22.82%	48.26%	0.0551
1550 nm	1.29%	2.08%	99.08%	53.63%	76.34%	40.94%	73.88%	0.1513
2000 nm	2.32%	1.70%	98.30%	48.01%	71.94%	34.54%	69.09%	0.1193
10.6 μm	1.64%	1.19%	98.85%	30.20%	57.71%	17.43%	76.68%	0.0668

图 4：选定设计的光学与协议指标。(a) 按基条件给出的量子比特误码率 QBER。(b) 条件概率保真度。(c) 平均透过率、四区域焦平面收集分数与收集效率。(d) 筛选后的保密分数，以及渐近模型下每次入射光子可获得的密钥比特数。

表 2 里藏着两个反直觉的例子。第一，10.6 μm 设计的保密分数最高（76.68%），但收集效率只有 17.43%，最终密钥量反而最低（0.0668）——误码漂亮并不等于出钥匙多，因为大部分光子根本没被探测到。第二，780 nm 处硅 P380 的收集效率低于二氧化钛 P360 候选（22.82% 对 31.50%），但它的响应分布更接近目标矩阵，预测密钥量反而更高（0.0551 对 0.0516）。论文据此强调：只拿概率保真度排序候选设计是不够的，必须把绝对透过率与”探测器区域怎么划”一起代入密钥率指标。

六态层析：把接收机当成待表征的测量算子

前面的密钥量都建立在四态光学基准上，即用同一套保密分数公式处理 H/V 与 R/L 两套基的结果。但真实接收机存在偏振相关损耗、两套基权重不平衡与错误端口泄漏，这些都可以用更严格的”测量算子”语言描述。

作者因此补测了 $\left|D\right\rangle=(\left|H\right\rangle+\left|V\right\rangle)/\sqrt2$ 与 $\left|A\right\rangle=(\left|H\right\rangle-\left|V\right\rangle)/\sqrt2$ 两个对角基输入，凑齐六个输入态，做接收机层析（detector tomography），把测量算子重构为”正算子值测度”（Positive Operator-Valued Measure，POVM），其中包括一个”无点击”算符 $\hat E_{\rm loss}^{\rm rec}=\hat I-\sum_j\hat E_j^{\rm rec}$，用来承接那些既没进匹配端口也没进正交端口、实际丢失掉的光子。重构获得的六态概率直接给出每个基的匹配点击概率 $p_b$ 与比特误码率 $e_b$，再配合虚滤波（virtual filter）与量子存储下的不确定性关系，得到相位误码率 $e_{{\rm ph},b}$ 的下界，于是每个基的保密分数变为

$$
r_b=\left[1-h_2\left(e_{{\rm ph},b}\right)-1.16h_2(e_b)\right]+ ,
\qquad
Y{\rm tom}=\tfrac12 p_Z r_Z+\tfrac12 p_Y r_Y,
$$

其中下标 $Z$ 指 H/V 基、$Y$ 指 R/L 基，$[x]_+=\max(0,x)$。这个式子不假设接收机的两个基正好各占一半权重，也允许两基错误率不同——这正是普通四态公式做不到的地方。

结果显示，六个被重构的接收机模型都给出正的渐近密钥界限：780 nm 为 0.06901，1550 nm P750 为 0.15164，2000 nm 为 0.13599，10.6 μm 为 0.06130；另外两个候选中，1550 nm P700 为 0.15935（在重构口径下取到最高值），780 nm 二氧化钛 P360 为 0.04846。注意 P750 与 P700 的排序在这里发生了反转：光学基准下 P750 的 0.1513 高于 P700 的 0.1471，而在重构接收机口径下 P700 反超。这不是笔误，而是两套计算口径的差别——一套用统一保密分数，另一套用逐基比特误码、算子推出的相位误码与各自的验收权重。论文没有把两种结果调和平滑，而是并列呈现。

亚波长交织的空间访问约束

把两套基以小于半波长的周期交织，除了让器件更紧凑，论文还给出了一个安全层面的理由。

如果测量基占据彼此分开、宏观可分辨的区域，那么空间上就存在一种潜在的攻击路线：用基选择性的掩模或照明，只作用于其中一套基的分析区域，把它照亮、遮挡或衰减，从而独立操纵一套基的统计。但当线偏振与圆偏振功能以 $p<\lambda/2$ 交替铺满整个孔径时，普通的远场掩模无法只分辨并透过其中一套宇称子晶格；若真按单元尺度打孔去分离，需要的开口本身就在亚波长尺度，结果是入射场被强烈衰减或散射，而不是形成一条通往单一基的独立传播通道。换句话说，亚波长交织在紧凑性之外，还把”两套基必须一起被照明”变成了物理约束。

需要说明的是，这是一条基于远场分辨极限的定性推理，论文并没有在这个模型下量化出对应的攻击上限。同一节里作者也交代了模型的信任假设：理想单光子源、可信接收机、渐近集体攻击，且不含信道衰减、探测器效率、暗计数、多光子发射与有限密钥效应。论文没有建模纠缠光源或成对接收机，但指出同一局部测量原则上可以放进 BBM92 协议，把两块这样的接收机分别放在纠缠对的两臂上。

意义与展望

这项工作把”超表面承担偏振测量”往前推了一步：此前的相关研究多是把两套基同时投影到正交空间模式、缩短测量时间，器件仍旧是分析器；而这里输出的直接就是 BB84 的四符号协议结果，且把亚波长交织带来的空间访问约束一起纳入论证。工程上，它给出了从纳米条查找表、相位库、四焦点布局、全波传播、标定区积分到密钥率的完整链路，并证明这条链路可以在 780 nm 到 10.6 μm 这样跨度很大的波段上复用——换的是单元库与焦点几何，不变的是测量目标与验证方法。

从原理验证走向实际应用，还有几件事需要推进。本文的全部结论来自全波仿真与协议模型，尚未有器件实测；把模拟的”内禀 QBER”换成实验数据，需要处理制备公差、光谱漂移与入射模式变化对固定模式表征的影响。模型层面，信道衰减、探测器效率、暗计数与有限密钥效应都在现有假设之外，诱骗态估计与更贴近实际的窃听模型也需要补上。在集成方向上，四个波长的设计说明这套方法可以覆盖通信波段与中红外波段，若能把光源、探测器与这块孔径一起做进同一平台，紧凑接收机的体量优势才会真正体现出来。总体而言，论文提供的价值在于把”器件形状”与”协议密钥率”用一条可复核的数值链连在一起，让接收机的取舍有了可比较的标尺。

参考资料

1. Yan He, Adetunmise C. Dada, Secure Quantum Metasurfaces for Direct Quantum Key Distribution Measurements, arXiv:2609.22946（2026）. https://arxiv.org/abs/2609.22946
2. 论文 arXiv 源包（LaTeX 源与矢量图件）：https://arxiv.org/e-print/2609.22946
3. Nicolas Gisin, Grégoire Ribordy, Wolfgang Tittel, Hugo Zbinden, Quantum Cryptography, Reviews of Modern Physics 74, 145 (2002). https://arxiv.org/abs/quant-ph/0101098
4. Mohamed ElKabbash, Metasurface-Based Dual-Basis Polarization Beam Splitter for efficient entanglement witnessing, arXiv:2510.03931（2025）. https://arxiv.org/abs/2510.03931
5. Peter W. Shor, John Preskill, Simple Proof of Security of the BB84 Quantum Key Distribution Protocol, Physical Review Letters 85, 441 (2000). https://arxiv.org/abs/quant-ph/0003004
6. OpenStax, College Physics — 27.8 Polarization, Rice University. https://openstax.org/books/physics/pages/27-8-polarization