VirtualLab专题实验教程-3.二维分束超表面光栅

---

前言

今天想给大家看的是二维分束超表面光栅，结合了讯技光电官方给出的技术文档以及lighttrans公司前些天刚推出的研讨会。限试用版可操作。如果你看过专题实验教程的1，你难免会好奇仿真出了一个纳米柱的相位调制性能后，超表面的设计到后面该怎么做，本篇就可以解答你的疑惑。另外如果你看过2，可能也会对接下来要讲的内容更加熟悉。
研讨会的录屏链接如下link

一、二维分束超表面光栅是什么？

笔者浅要分析一下为什么说激光雷达会需要用到这个二维分束超表面光栅。

1、多线激光雷达及其当前存在的问题

我们可能经常会在市场上听到多线激光雷达的名头。它其实是指同时发射及接收多束激光的激光旋转测距雷达，市场上目前有4线、8线、16 线、32 线、64 线和128线之分，多线可以识别物体的高度信息并获取周围环境的3D扫描图。相比单线激光雷达自然是有一个多维的优势。单线扫一圈只会扫一行点，多线扫一圈会扫好几行（对应多线激光雷达的线束）点。
但是目前大部分多线激光雷达是存在弊端的：（在《基于阵列式激光器和探测器的多线激光雷达》这篇专利 link中写道）目前的多线激光雷达，无论Velodyne的16/32/64线，还是IBE0的4/8线激光雷达， 都是1个激光探测器对应1个激光发射器，有多少线就使用多少个激光器和探测器的组合。 每一线，也即每一个激光器和探测器的组合，均需要单独调整、单独固定。但如此装调难度是很大的，每一个激光器和探测器的对应需要单独调整对准，装调难度大、耗时长、生产效率低；并且因为每个激光器都需要驱动电路和机械固定也会导致整个激光器的体积变大。
专利中采取了下图这种解决措施，所述阵列式激光器1从整个晶片上切割而来包括至少二个激光器11，同理，阵列式探测器5从整个晶片上切割而来，包括至少二个探测器51，探测器51的数量与激光器11的数量对应，如此实现了激光器阵列与探测器阵列的晶片一体化。

如果这时你能猜到我们的二维分束超颖光栅要用在何处，你可能就会发现，最终实现的原理图竟有可能和上图非常相似。

2、将超颖光栅用到发射端

是的，我们应该要用在发射端上。不是多个单独的激光器，也不是专利中的集成在一片镜片上的多个激光器，而是一束光分成多束，从而实现一个等同于有多个激光器的效果。
将发射光分为x束，同时照射目标上x点。从x点反射回的信号投射到对应的x元探测上，通过处理可得到x个像素数的距离信息和强度信息。通过扫描器扫描，直到探测到所要求的像素数时，成像显示。
虽然说是这么说了，但是这项研究难度较大，对扫描器的要求较高，信号的探测也较为困难。那这样看来，在扫描器和信号探测上优势最大的反倒是单线激光雷达啦，它的技术会更成熟，但毕竟超表面这块也才刚发展起来，还是很有研究价值的。

待续。

二、2021试用版仿真

1、制作相位分布图

首先需要打开一个分束器的设计文档，进入衍射工具箱，这一步的最终目的是得出最终想要的相位分布图。

设计为940nm的工作波长（红外夜视灯采用的）

改成角谱模式。

希望目标的分束效果如下。将入射光分到3×3共9个级次上，各级次角度之差为30°。


这里取消勾选对相位的一个量化要求（不考虑是几个层级的量化）。

设置像素单元，即为纳米柱结构的单元。（400nm×400nm）

分束器设置结束，弹出IFTA设计窗口，交由傅里叶迭代法自行设计相位分布。

算好后查看结果，改成相位模式查看。
（每一次测算后的相位分布图的结果都不尽相同），最终我们选用如下的相位分布图作为后期纳米柱布局的依据。

2、超表面纳米柱直径仿真

接下来是构建出超表面的环节，我们进入光栅工具箱，选用超颖光栅的光路图搭建。

光源的波长设置为940nm。（大家应该很熟练了不放图了）
基板的厚度为0，选用折射率为1.5的材料。


使用后面的堆栈，导入超光栅作为堆栈的基础设置。

编辑该超光栅，这里先只设置一个单元（我们即将讨论的是改变这个单元纳米柱的直径，会有什么样的相位延迟，类似上一期）

进入pillar medium下方的小铅笔进行编辑，我们的纳米柱高度为465nm，纳米柱材料的折射率为3.8。

调整纳米柱的分布，移除原有设置中默认的两项，添加新的一个纳米柱，添加成现在200nm直径，在单元正中心的纳米柱。


这里设置的是当前超表面的大小，因为只有一个纳米柱即为纳米住单元的大小，400nm*400nm。

设置完成后看看上方缩略图显示的是否正确。

其实这个时候去进行参数仿真已经可以得到比较不错的平滑曲线，但为了让曲线更平滑，我们微调一下算法。这里调整FMM的设置。

我们把结构分解改为手动设置，点与点之间的距离设置成1nm（更精细）。

到此，超颖光栅设置完成。光栅级次分析器的设置方法同上次试验: link。
我们按照上次超表面纳米柱的案例的思路，对它进行一个参数仿真，也即改变纳米柱直径，观察相位和透过率的曲线变化。（单元大小为400nm×400nm，我们从1nm取到399nm，步长取1nm）

Next至该界面后，只观察光栅级次分析器的结果。

仿真结果得到两张图，分别对应0级透过效率与瑞利系数（手册中解释瑞利系数为由傅里叶模态法计算出来的对应于各个光栅级的复振幅）

按照上次超表面教程的内容把两幅图合并起来。

调一下图例的位置，我们大致选用这个直径范围内的纳米柱（覆盖了2$\pi$的相位，同时具有比较高的透过率）

接下来需要把相位和纳米柱直径一一对应起来。（其实有编程模块，可以编写模组来查询，但我们没有，就采取人工查询法了）。之前得到的相位图的相位分布数值如下表（单位/rad）。

再回到提取的相位图，在刚刚选定的纳米柱直径范围内取点，去找到能符合各单元的纳米柱直径大小，记录25个单元在各自位置上的纳米柱的直径取值。

去索引纳米柱直径后得到下表（单位/nm）。

3、超颖光栅的正式设计

我们回到超颖光栅的设置界面，对它进行编辑，把刚才记录好的各个直径的纳米柱都放置到平面上。
我们重建一个超颖光栅的metagrating optical setup，重新设计。
同样先设计光源，波长为940nm。
超颖光栅元件的设置还是选择跟2中几乎相同的模式。

不同之处在于，这次纳米柱的分布要做详细地调整了。
在编辑窗口，依然进入第一个平面下pillar medium的小铅笔进行编辑，pillar distribution中添加25个纳米柱单元。这样添加出来的纳米柱即为按照相位需求添加的。

此外需要修改光栅周期，以放下这么多纳米柱为准，也即2um*2um。

傅里叶模态法的设置也需要再做修改，倏逝波的级次设为37。

结构分解上采用和纳米柱仿真相同的设定。

可以从上图右下角的decomposition view这里查看纳米柱分布情况。

光栅级次分析器这里更改一下设置，现在需要看到真正的级次的一个衍射的效果，而不是获得透过率等的信息。把级次采集再次勾选上，我们要观察的3×3级次，因此单个级次设置这里注意范围的设置。

才外还需要进入编辑窗口修改材料折射率为1.5（和基板的折射率相同）。也即平面波到超颖光栅之间没有经过其他材料，认为其是紧贴着超颖光栅的光源。

超颖光栅往后的材料我们就设置成空气就好，air in homogeneous medium。

然后选用光栅级次分析器作为仿真引擎，得到结果。

也可将观察平面选为笛卡尔角度。

结果变成下图。

在下方的结果框中可以查看各个级次的衍射效率。

附

之前有朋友说这样仿真出来的结果和案例文件的结果不一致，研究后发现在于编辑窗材料设置的问题，请大家务必注意，已经更正在正文里了。另外还有朋友遇到说显示的级次不只3*3，可能原因如下（我们问了原厂高级光学工程师Hartwig Crailsheim）。
这里order selection strategy这里没有选成above efficiency threshold，遇到这个问题的朋友不妨再试试看呢。

总结

这也是超表面的重要应用之一了（因为在表面上实现了相位跃变，根据广义的斯涅耳定律，这使得光束原本的传播方向发生了改变，因此实现了分束），超表面的神奇世界值得下一篇再开一个小小的探索。欢迎大家尝试后互相交流！