文献解读 | 超表面集成 VCSEL 的偏振调控：降低各向异性把光学伊辛机精度提高一倍（arXiv 2026）

论文简介

题录：Wenjie Chen, Zifeng Yuan, Hong-Lin Lin, Luo Qi, Jiaru Chu, Aaron Danner, Yuhang Chen, Metasurface-integrated VCSEL designed for polarization control in optical Ising machines, arXiv:2609.21331（2026-09-18 投稿，physics.optics）。原文链接：https://arxiv.org/abs/2609.21331

作者机构：中国科学技术大学精密机械与精密仪器系、安徽省高等院校精密科学仪器重点实验室，以及新加坡国立大学电子与计算机工程系。

一句话贡献：作者把一层矩形纳米柱超表面刻蚀进 VCSEL（Vertical-Cavity Surface-Emitting Laser，垂直腔面发射激光器）的顶部反射镜，用几何尺寸去补偿两个正交偏振之间的反射率差，从而削弱 VCSEL 与生俱来的偏振偏好；在数值模型中，这种器件在 100 组随机耦合矩阵上的 3 位伊辛机求解精度由 30% 提升到 62%。

这篇文章的有意思之处，不在于”又做了一块超表面”，而在于它跨了两层：一面是纳米结构的光学设计，另一面是组合优化问题的物理求解。要理解它为什么值得做，得先弄清楚光学伊辛机到底把”比特”存在哪里。

光学伊辛机为什么需要一对简并的偏振态

伊辛模型原本描述磁性材料：一个由若干自旋组成的晶格，每个自旋只取 $+1$ 或 $-1$ 两个值。当这些自旋通过相互作用耦合在一起，系统会自发向能量最低的组态演化，这个最低能量就是哈密顿量

$$H=-\sum_{i<j}^{M}J_{i,j}\sigma_i\sigma_j+\sum_{i=1}^{M}h_i\sigma_i$$

的最小值。其中 $\sigma_i$ 是第 $i$ 个自旋的取值，$J_{i,j}$ 描述自旋 $i$ 与 $j$ 之间的耦合强度与符号，$h_i\sigma_i$ 是环境对单个自旋的外场项（Zeeman 项）。$M$ 个自旋代表 $M$ 个二值变量。

数学上的关键在于：二次无约束二值优化（Quadratic Unconstrained Binary Optimization，QUBO）问题的求解，可以等价地写成寻找伊辛模型基态。于是，”用物理系统自然演化的方式找到基态”就成了求解组合优化问题的一条捷径。伊辛机的实现路线大致有三类：电路式器件、光子器件，以及利用量子力学原理（量子退火）的器件。

光子伊辛机的思路是让光的某个二值物理量直接充当自旋变量，候选量包括偏振、幅度和相位。目前已有多种平台被尝试过：注入锁定激光器、带并行光反馈的 VCSEL 阵列、光参量振荡器、多芯光纤、微环谐振器和马赫—曾德尔干涉仪等。它们的共同点是利用光学系统本身的并行性和低延迟，让系统”自己”滚到低能态，而不是靠数字电路一步步迭代。

用激光偏振编码的自旋，有一个天然的好处：两个正交偏振态是同一个激光腔里的两种模式，它们共享增益介质、共享腔体，物理上最接近理想的”简并二能级”。当两台激光互相注入（即把各自的输出送入对方腔内）时，注入光会牵引对方的偏振；如果注入足够强，被注入激光就会”翻转”到注入偏振上。耦合系数 $J_{i,j}$ 的符号和大小，可以通过衰减器和半波片（Half-Wave Plate，HWP）在光路上设定，例如 $J_{i,j}>0$ 时两台激光倾向于同向偏振、$J_{i,j}<0$ 时倾向于正交偏振。一群这样互相注入的激光，会在模式竞争中自发收敛到一个低能组态——这就是”光在算”，而不是”电路在算”。

代价也随之而来：这套机制假设两个偏振态尽可能对称。如果器件本身强烈偏好某一个偏振，外来的注入光就很难把它的偏振掰过去，系统会在到达正确解之前就”卡”在偏好态上。

VCSEL 的”天生偏好”从哪里来

VCSEL 的结构像一个”垂直的三明治”：上下各一组 DBR（Distributed Bragg Reflector，分布式布拉格反射镜），中间夹着有源区。DBR 由折射率交替的薄膜对堆叠而成，每层厚度设计为四分之一波长，使反射光同相叠加形成高反射。本文的器件是 GaAs 基结构，标称激射波长 850 nm，由 16 对顶部 DBR 和 40 对底部 DBR 构成，每对 DBR 包含铝组分不同、折射率分别为 3.06 和 3.49 的两层 AlGaAs。

问题就藏在这套”垂直结构”里。VCSEL 的谐振腔在横向上并不是各向同性的：增益、腔长和层厚的微小差异都会在 X、Y 两个正交偏振之间造成差别。论文在自旋翻转模型（spin-flip model）里用两个参数来描述这种各向异性：腔二色性 $\gamma_a$（描述两偏振方向的损耗/增益差，也称线性二色性）与腔双折射 $\gamma_p$（描述两偏振方向的相位差，也称线性双折射）。

各向异性最直接的后果，是改变两个偏振的阈值增益。论文用有限差分时域法（Finite-Difference Time-Domain，FDTD）计算了不同顶部镜结构在 850 nm 附近的反射率：在纯 DBR 构成的顶部镜中，Y 偏振的反射率高于 X 偏振，$R_y>R_x$。反射率高意味着透出腔外的光少、腔内光场更容易积累，阈值增益更低——于是激光自发地先在一个偏振上起振，本文的情形就是 Y 偏振。

而在偏振编码的计算里，这种”偏好”直接把正确答案排除掉了：当伊辛问题的解要求某个比特取另一个偏振时，器件会抵抗翻转。论文的量化结论是，各向同性 VCSEL 在注入比 $E_{inj}/E_o$ 超过 0.0025 时就能锁定到注入偏振，而各向异性 VCSEL 需要把注入比提高到 0.07，约为前者的 28 倍。注入功率在系统中是稀缺资源——每台激光的输出还要分给多个通道去实现互注入，因此”更容易被锁定”直接等于”同等规模下算得更省”。

于是设计目标就明确了：不是把 VCSEL 做成完全没有各向异性（真实器件做不到），而是造出一个可以补偿 $R_x$、$R_y$ 差异的反射镜结构，把各向异性压低到不影响计算的程度。

超表面怎样改写顶部分布式布拉格反射镜

本文的做法是在顶部 GaAs 层里”刻”出一层超表面：一层矩形纳米柱阵列，等效光程设计为四分之一波长，直接嵌在顶部镜堆内部。这里的角色分工很清楚——DBR 提供高反射的基础，超表面负责在两个正交偏振之间制造一个可控、方向合适的反射率差，把纯 DBR 的 $R_y>R_x$ 补回来。

图 1：超表面集成 VCSEL 的结构示意，左上角插图为超表面一个周期的放大视图。刻蚀在顶部 GaAs 层中的矩形纳米柱阵列嵌在顶部镜堆内，下方依次是有源区与底部 DBR。

器件的设计自由度就落在纳米柱的几何尺寸上。论文沿用一条简单的缩放策略：沿 Y 方向的周期设为 $P_y=S\times200$ nm，纳米柱宽度设为 $W=S\times100$ nm，而沿 X 方向的周期固定为 $P_x=1000$ nm、纳米柱长度固定为 $L=800$ nm；$S$ 是统一的缩放因子，出射方向沿 Z 轴。改变 $S$ 就等于同时改变两件事——柱子的宽窄，以及它们在 Y 方向上的疏密。

表 1：超表面单元的设计定义与缩放因子。单元格数值取自论文正文与图注，$S$ 为统一缩放因子。

参数	符号	取值	说明
X 方向周期	$P_x$	1000 nm（固定）	不随 $S$ 变化
Y 方向周期	$P_y$	$S\times200$ nm	随 $S$ 缩放
纳米柱长度	$L$	800 nm（固定）	沿 X 方向
纳米柱宽度	$W$	$S\times100$ nm	随 $S$ 缩放
等效光程	—	$\lambda/4$	集成在顶部镜堆内的设计目标
设计波长	$\lambda$	850 nm	标称激射波长
双反射极值	—	$S=0.8$ 处最大	随 $S$ 增大到 1.25 逐步减小，之后再次上升

FDTD 计算给出了这套几何设计的效果。在纯 DBR 顶部镜中，Y 偏振的反射率高于 X 偏振；把超表面集成到各向同性 DBR 层之上后，关系被反转：X 偏振的反射率反超 Y 偏振。更关键的是，这种由超表面引入的双反射（bireflectance）可以在 850 nm 处通过缩放因子连续调谐——$S=0.8$ 时补偿量最大，随着 $S$ 增大到 1.25，补偿逐步减弱，之后又再次上升。也就是说，纳米柱的宽窄与间距提供了一个连续的”偏振天平”旋钮。

图 2：超表面集成顶部 DBR 的性能，绿色与橙色曲线分别表示 X、Y 偏振入射下的反射率。(a)(b) 分别为各向异性 DBR 层与超表面集成各向同性 DBR 层的双反射谱；(c) 为 850 nm 处不同缩放因子 $S$ 下超表面集成各向同性 DBR 层的双反射；(d)–(f) 为 850 nm 附近超表面集成各向异性 DBR 层的双反射谱，从左到右 $S$ 依次为 0.8、1.0、1.1。

至于双反射的物理来源，论文用近场电场分布给出了一个直观解释：在 X 偏振入射下，能量主要沿纳米柱的短边聚集；而在 Y 偏振入射下，相邻纳米柱之间的耦合增强，更多能量被局域在柱子附近、反射回腔内的比例下降。同一块结构对两个偏振呈现出不同的场分布，因此呈现出不同的反射率——这是形状各向异性在近场层面留下的痕迹，而不是某种材料吸收差异。值得一提的是，超表面不仅影响阈值增益（反射率差），论文的补充材料还显示它同样可以调节 VCSEL 的双折射，意味着 $\gamma_a$ 与 $\gamma_p$ 两个方向的各向异性都能被结构设计触及。

一块真的做出来的超表面

论文的实验部分相当克制，也需要读者注意它的边界：作者在 GaAs 衬底上制备了 $S=1$ 的超表面样片，用扫描电子显微镜（Scanning Electron Microscope，SEM）检查形貌，并测量了 900–1000 nm 的透射谱。SEM 显示整体结构与设计接近，实测纳米柱长 770 nm、宽 140 nm、高 61.5 nm，沿 X、Y 方向的周期分别为 1000 nm 与 210 nm。

两个数字值得停一下。设计值（$S=1$）为 $P_y=200$ nm、$W=100$ nm，而实测宽度为 140 nm、Y 方向周期为 210 nm，说明刻蚀后的柱子比设计更”胖”一些、间距略大一些；论文对此的表述是形貌与设计相符，没有给出这一偏差对双反射影响的定量分析。另外，实测透射谱覆盖的波段是 900–1000 nm，而器件设计波长是 850 nm——这为”结构确实具备偏振相关响应”提供了实验证据，但并不等同于在 850 nm 处实测到了设计所需的双反射补偿量。

更需要明确写在明面上的，是这项工作的验证层级：超表面与 VCSEL 的一体化制造尚未实现，论文明确把集成制造列为后续工作；此后所有关于注入锁定与伊辛机精度的结论，都来自自旋翻转模型与 FDTD 的数值计算，而不是对整机器件的测量。换句话说，这是一篇”器件设计 + 系统级数值验证”的工作：实验部分证明了关键单元可制造、具备偏振选择性，系统收益则由模型给出。

各向异性越低，偏振越容易被注入锁定

数值部分的第一步，是看单台激光在外部注入下的偏振演化。模型中的各向异性由 $\gamma_a$ 与 $\gamma_p$ 描述：各向同性情形取 0；各向异性情形取 $\gamma_a=\pm0.67$ ns$^{-1}$、$\gamma_p=\pm192$ ns$^{-1}$。外部注入在 50 ns 时刻开启，各向同性器件使用的注入比为 $E_{inj}/E_o=0.01$，各向异性器件为 0.1。

图 5：外部注入下 VCSEL 激射偏振的时间演化。(a)(b) 为各向同性 VCSEL 在 X 偏振与 Y 偏振注入下的时间响应；(c)(d) 为各向异性 VCSEL 在 X 偏振与 Y 偏振注入下的时间响应，其中 $\gamma_a$ 分别取 $0.67$ ns$^{-1}$ 与 $-0.67$ ns$^{-1}$，$\gamma_p$ 分别取 $192$ ns$^{-1}$ 与 $-192$ ns$^{-1}$。

结果与直觉一致，但给出了具体的门槛：各向同性 VCSEL 对两种偏振的注入响应完全相同，只要注入比超过 0.0025，偏振就会对齐到注入光；各向异性 VCSEL 则需要 0.07 才能完成翻转，是前者的 28 倍。$\gamma_p$ 的符号在这里决定了自由运转时”偏好”哪个偏振，而它的绝对值越大，器件维持初始偏振、抵抗外部翻转的能力就越强。补充材料中周期性注入的结果进一步印证了这一点：$|\gamma_p|$ 从 192 ns$^{-1}$ 依次增大到 215、500 ns$^{-1}$ 时，从”完全被锁定”逐步过渡到”周期性地在两个偏振之间切换”，最终几乎无法被锁定。

这里的物理图像可以概括成一场”牵引力对回正力”的竞争：注入光提供把偏振掰向目标方向的牵引，而腔二色性与双折射提供把偏振拉回偏好方向的回正力矩。降低各向异性，相当于削弱回正力矩，于是同样强度的注入就能主导偏振状态——这正是”超表面补偿反射率差”最终要换来的东西。

三位伊辛机的精度：30% 与 62%

数值验证的最后一层，是把三台这样的激光接成 3 位伊辛机：每台 VCSEL 的偏振态代表一个自旋，$\sigma=+1$ 对应 X 偏振、$\sigma=-1$ 对应 Y 偏振；三台之间通过光路互注入，衰减器与半波片设定耦合矩阵 $J$ 的强度与符号，外场项 $h_i$ 取 0，互注入从 50 ns 开始，所有组合的互注入耦合系数统一取 3.55 ns$^{-1}$，频率失谐取零。

论文先给出一组具体的耦合矩阵用于演示，它的两个正确解是 $[+1,-1,-1]$ 与 $[-1,+1,+1]$（这两个组态对应同一个低能态的两个整体翻转版本）。随后作者换了三档各向异性参数，观察系统如何收敛。

图 6：VCSEL 各向异性对 3 位光子伊辛机计算精度的影响。(a) 光学装置示意图；(b)–(d) 在随机相互作用矩阵 $J$ 下，由不同各向异性 VCSEL 构成的 3 位伊辛系统的时间演化。

表 2：三档各向异性参数下 3 位伊辛机的数值表现。精度为 100 组随机对称耦合矩阵（对角为零、非对角元在 $[-1,1]$ 内均匀分布）上的统计结果。

器件档位	腔二色性 $\gamma_a$	腔双折射 $\gamma_p$	时间演化表现	计算精度
低各向异性	$0.67$ ns$^{-1}$	$192$ ns$^{-1}$	容易演化到正确解	62%
中等各向异性	$0.67$ ns$^{-1}$	$500$ ns$^{-1}$	出现起伏但仍能收敛	精度下降
高各向异性	$6.7$ ns$^{-1}$	$1920$ ns$^{-1}$	倾向于锁定在偏好偏振	30%

结论很直接：随着各向异性增大，系统从”轻松收敛”变成”有起伏但仍收敛”，再变成”锁在偏好偏振、算不出来”。在 100 组随机耦合矩阵上的统计结果是精度从 62% 降到 30%，也就是低各向异性器件的精度提高了一倍以上（$62/30\approx2.07$）。

论文还强调了两点并不显眼、但对规模化重要的影响。其一，降低各向异性同时降低了完成偏振翻转所需的注入功率。由于每台激光的输出都要经过分束送到多个通道才能实现互注入，注入功率是随比特数急剧消耗的资源，”同等耦合强度下更省功率”意味着更容易向更大规模的系统扩展。其二，当伊辛问题存在多个等价解时，器件初始的各向异性会成为解空间里的”偏向力”，把系统推向某一个特定解；各向异性越大，这种偏向越明显。也就是说，各向异性不只是影响”算得准不准”，还影响”在多个正确答案中选哪一个”。

把数字代进去核对一遍

读这类”设计 + 仿真”的论文时，最实用的检查是把几个关键数字重新算一遍，看看量级是否自洽。下面三处核对都能用文中的参数独立完成。

第一处：四分之一波长的尺寸自洽。 超表面的等效光程被设计为 $\lambda/4$，设计波长 $\lambda=850$ nm。若把实测的纳米柱高度 $h=61.5$ nm 当作等效的四分之一波长厚度，反推等效折射率：

$$n_\mathrm{eff}=\frac{\lambda}{4h}=\frac{850\ \mathrm{nm}}{4\times61.5\ \mathrm{nm}}\approx3.46$$

这个值与论文给出的 DBR 材料折射率区间（AlGaAs 的 3.06 与 3.49）相符，说明”61.5 nm 的高度”与”850 nm 下的四分之一波长设计”在量级上自洽，而不是两组彼此矛盾的参数。

第二处：各向异性参数的时间尺度。 $\gamma_a$ 与 $\gamma_p$ 的单位是 ns$^{-1}$（速率），取倒数就是时间尺度：$1/\gamma_a=1/0.67\ \mathrm{ns^{-1}}\approx1.5$ ns，$1/\gamma_p=1/192\ \mathrm{ns^{-1}}\approx5.2$ ps。前者给出二色性导致偏振弛豫的纳秒量级时间，与图 5、图 6 中”注入自 50 ns 开始、系统在数十纳秒内完成翻转或收敛”的观测窗口一致；后者说明双折射带来的相位演化比注入过程快三个数量级——在一个 50 ns 的观测窗口内，相位累积的圈数约为 $\gamma_p t/2\pi=192\times50/2\pi\approx1.5\times10^3$ 圈，因此模型观察到的是快速相位平均之后的偏振竞争结果。这也解释了为什么各向异性档位需要放大到 $\gamma_p=1920$ ns$^{-1}$ 才能明显拖垮计算精度。

第三处：阈值与精度的比例关系。 注入锁定阈值由 0.0025 提高到 0.07，比值为 $0.07/0.0025=28$；精度由 62% 降到 30%，比值为 $62/30\approx2.07$。两个数字描述的是不同层面的代价：前者是单台器件的物理门槛（需要多大的注入才能翻转偏振），后者是系统层面的统计结果（100 组随机问题里解对了多少次）。它们不能互相换算，但方向一致——各向异性越大，偏振越难被外部主导，系统越容易停在错误组态上。

意义与展望

这篇文章的价值，在于把”计算精度”这一系统级指标，与”纳米柱有多宽、间距多大”这一结构级参数连成了一条可设计的链条：几何尺寸 → 两个偏振的反射率差（双反射）→ 阈值增益差 → 腔二色性/双折射 → 注入锁定门槛 → 伊辛机求解精度。对做偏振编码计算的人来说，这条链条提供了一个具体的器件侧抓手：不需要换掉 VCSEL 工艺，而是在顶部镜堆里加一层可设计的光学结构。

它同时给出了一个有用的方法论示范：当计算系统的瓶颈落在”物理对称性”上时，超表面这类亚波长结构可以充当”精密调平”的工具，把几十纳米尺度上的不对称补偿掉。这样的思路并不限于伊辛机，凡是依赖简并态的物理计算方案（例如依赖正交模式、双稳状态的方案）都可能需要类似的补偿手段。

从原理验证走向实际应用，还有几条明确的推进方向。首先是集成制造：论文目前只制备并表征了超表面样片，把这层结构稳定地做进真实的 VCSEL 腔面、并保证与有源区对准，是下一步必须跨过的工程门槛。其次是实验闭环：注入锁定门槛、$\gamma_a$ 与 $\gamma_p$ 的实际取值、以及系统精度，都需要在整机器件上实测，才能确认数值模型里”62% 对 30%”的差距能有多大程度落到现实。第三是规模与一致性：本文的验证是 3 位、器件参数理想一致、频率失谐为零的小系统；当比特数增加、注入功率需要在通道间分配、器件之间存在参数离散、温度变化带来波长漂移（论文也提到更高注入功率会引起热致红移、从而增大主从激光之间的失谐）时，如何在整片阵列上维持低各向异性与高强度互注入，会成为决定可扩展性的关键问题。

参考资料

1. Chen W., Yuan Z., Lin H.-L., Qi L., Chu J., Danner A., Chen Y. Metasurface-integrated VCSEL designed for polarization control in optical Ising machines. arXiv:2609.21331, 2026. https://arxiv.org/abs/2609.21331
2. Aadhi A., Al-Kayed N., Austin T., St-Arnault C., Rotenberg N., Ballani H., Plant D. V., Ng E., Yamamoto Y., McMahon P. L., Shastri B. J. Photonic Ising machines toward and beyond a million spins. arXiv:2607.13446, 2026. https://arxiv.org/abs/2607.13446
3. Yuan Z., Zhang D., Lin H.-L., Danner A. Engineering Polarization Switching in VCSELs with Custom Aperture Shapes. arXiv:2608.16597, 2026. https://arxiv.org/abs/2608.16597
4. Xie Y.-Y., Ni P.-N., Wang Q.-H., Kan Q., Briere G., Chen P.-P., Zhao Z.-Z., Delga A., Ren H.-R., Chen H.-D., Xu C., Genevet P. Metasurface integrated Vertical Cavity Surface Emitting Lasers for programmable directional lasing emissions. arXiv:1906.12078, 2019. https://arxiv.org/abs/1906.12078
5. Yuan Z., Shan W., Chen T., Lin B., Danner A. Mesa Orientation Engineering for Polarization Locking in VCSELs. arXiv:2608.24075, 2026. https://arxiv.org/abs/2608.24075