Expanding RCWA capabilities with advanced S-matrix algorithms

当光在真空空间中传播时，其行为是可预测的；但一旦遇到薄膜堆叠或具有微小结构的表面，规则就会变得复杂。这就是层状光子结构领域，它是防反射涂层、镜面以及光纤网络内部复杂组件等现代光学器件的核心。当光线击中这些层时，它会在每个边界处发生反射、干涉和散射，从而产生一个难以预测的能量交织网。为了设计更好的器件，工程师需要一种能够精确计算光线穿过这些堆叠层时如何运动的方法。几十年来，科学家们一直依靠一种被称为散射矩阵法（scattering-matrix method）的强大数学工具来解决这个问题。这种方法之所以备受推崇，是因为即使在处理极厚层或光线迅速衰减的情况时，它也能保持稳定和精确，而传统的计算方法在这些情况下往往会失效。然而，尽管该工具在处理标准问题时表现良好，但其在探索复杂设计和理解光流深层物理机制方面的全部潜力尚未被完全释放。

来自法国波尔多大学和光学研究所（Institut d'Optique）的一个研究小组现在改进了这一工具，将其转化为一种更灵活、更强大的设计与发现工具。他们引入了一种定义散射矩阵的新方式，将重点从两层之间的连接转向了单层本身的属性。在传统方法中，针对特定层的计算都与紧邻其上方和下方的材料相关联，这意味着如果工程师想在不同的堆叠结构中测试同一层，就必须从头开始重新计算。新方法将每一层视为一个独立的实体，具有其自身的内在行为，这种行为由其处理穿过它的光线的方式来描述，而与其邻居无关。层与层之间的界面则被单独处理。这种分离使得层的属性可以被计算一次后，在不同的组合中进行无限次复用，就像一块可以在任何墙体中使用而无需为每个新设计重新塑形的砖块一样。

这种视角的转变不仅提高了计算速度，还揭示了这些结构内部隐藏的光路径。通过使用这种重新构建的方法，研究人员现在可以直接计算外部光波如何与周期性结构中自然存在的内部波（即布洛赫模态，Bloch modes）进行耦合。这些模态是光在图案化材料中传播的特定方式，类似于吉他弦只能以特定的频率振动。新方法提供了一组广义系数，它们就像一张地图，精确地展示了多少能量从入射光束转移到这些内部模态中，以及它们在边界处是如何相互作用的。为了证明这一点，团队分析了一个旨在实现宽波段高透明度的硅光栅。以往的解释侧重于单个纳米线产生的共振，而新的分析表明，透明度实际上是由仅有的两种特定内部光模态的激发和传播驱动的。一种模态主导了长波长下的行为，而第二种具有极低有效速度的异常模态则负责了短波长下的抗反射效应。这种细致程度让设计者能够看清究竟是哪些内部通道在承载能量，从而将模拟从一个只输出数值的“黑箱”，转变为观察器件物理现象的透明窗口。

这种重新构建带来的实际效益同样显著，特别是在设计优化的速度方面。由于新方法将层的厚度视为解析公式中的一个简单变量，它计算数百种不同厚度的光学响应所需的时间，仅相当于以前计算单一厚度所需的时间。研究人员通过模拟一个由二氧化硅间隔层分隔的两个硅光栅结构来测试这一点，并在其中改变了上层光栅的 400 种不同厚度值。整个扫描过程仅用了完成单次计算的一小部分时间，而对于需要为每个新参数进行完整运行的传统模拟方法来说，这是不可能实现的。这种能力还为处理厚材料中光相干性丢失的问题提供了一个巧妙的解决方案，这在太阳能电池和其他现实应用中是一个常见问题。研究团队并没有构建复杂的理论来处理这种随机性，而是展示了只需对许多厚度值的结果进行平均，就能得到准确的非相干行为图像，从而有效地平滑了由干涉引起的快速波动。

此外，该新框架允许同时分析从许多不同角度入射到结构上的光。传统的软件通常只计算单束光的计算结果，需要针对每个感兴趣的角度分别运行。然而，更新后的方法可以在一步之内计算出所有可能的入射和衍射方向的散射系数。这对于分析利用人工周期性建模的大型无序表面特别有用，因为单次计算即可揭示整个系统在各种照明条件下的行为。所有这些功能都已集成到一个名为 RETICOLO 的免费软件软件包中，作者已将其向科学界开放。通过提供用于生成其结果的详细代码和脚本，研究人员提供了一个实用的工具包，使其他科学家能够进行这些先进的模态分析，加速其设计流程，并对光如何在驱动我们光学技术之复杂层状结构中运动获得更深层的物理理解。