FDTD叠层结构光学计算概述

      人们通常可以使用分析方法来确定多层堆叠的光学响应，而不是直接模拟麦克斯韦方程组。Lunerical Stack光学计算器提供了一系列脚本函数（stackdipole，stackpurcell）来进行这些计算。此外，FDTD中为STACK仿真模型为脚本命令提供了GUI（graphical user interface，图形用户界面）。

      Stackrt和stackfield脚本命令对于分析多层结构也十分有用。这些命令包含在所有的Lumerical的产品中。

1. 用于平面波照明的多层堆栈（Stack for Plane Wave Illumination）

Stackrt脚本命令将会使用分析转移矩阵方法，来返回一个平面波通过任意多层堆栈结构的反射率和透过率。该函数返回S和P偏振光的透过的和反射的能量的比例（Ts, Tp, Rs, Rp），和 复数反射和透射系数（ts, tp, rs,rp）。Stackfield脚本命令返回堆栈结构的场信息（Es, Ep, Hs, Hp）。

       对于stackrt和stackfield脚本命令，用户必须制定光线的入射角和频率范围，还有每层结构的厚度和折射率。为了对比采用分析转移矩阵方法和一维FDTD模拟获得的结果，可以参考4层堆栈结构应用示例（https://optics.ansys.com/hc/en-us/articles/360042089573）。

(1) 复数系数

        复数系数定义为反射或透射后输出波的振幅除以入射波的振幅。界面处平面波的菲涅尔分析将3D空间的矢量问题缩减为1D，其中所有重要信息包含在复数系数中。有关坐标的更彻底的处理可以参考Source Rotations in 3D FDTD (https://optics.ansys.com/hc/en-us/articles/1500002383802)。维度上的缩减是一个非常有用的简化；然而，它依赖于一些常规，这些常规在解释结果时很重要。应该理解两个重要的概念：入射面和偏振态。

（2）入射面

   为了充分指定结果，我们需要选择要使用的坐标系。尽管坐标系是任意的，但利用平面波是横电磁波这个事实是有帮助的，这意味着波矢k，电场E和磁场B是正交的。显然选择沿着主轴排列的坐标系将简化大大简化分析。历史上曾考虑过不同的选择，但大多数现代分析都使用了入射平面的有用概念。入射平面定义为包含波矢量k （波传播方向）和表面法线n的平面。这样便可以进行两个直观的线性偏振态S-pol (TE,横电波：电场方向垂直于入射面，磁场方向平行于入射面)和P-pol（TM,横磁波：磁场方向垂直于入射面，电场方向平行于入射面），在分析时它们需要单独处理。这些基本状态分别定义为具有入射面外的电场E或者磁场B（即TE和TM的定义）。一旦确定了这些，第三个分量将根据Maxwell方程得出，并且在入射平面内。这些基础信息可以在下图中看到。

 （3）基本约定（Basic Convertions )

在反射或折射后，我们可以看到电场E，磁场B和波矢k将经过一个简单的旋转。旋转的角度可以由Snell定律轻易获得，这些角度将会由计算器反馈。大部分著者同意在上图中的定义；然而，应该被强调的是在选择反射P偏振光的基准存在些变动。这个基本可以归结为根据反射来选择E和B的正向方向，即箭头指向那个方向。我们使用在上图ii中的定义，但是许多重要的光学教科书比如Hecht, Born,和Wolf使用了一个基准电场E朝下和磁场B射出纸面（为圆点，与图ii中的x相反）。这个效应可以视为额外增加了一个π相位因子，或者对反射系数rp乘以一个-1的因子，这有时会引起混淆。

（4）传播相量 Propagation Phasor

复数系数中包含来自于界面的相位信息以及界面间的附加相量传播。当考虑多个反射时，比如在微腔中，重要的是要考虑累加相位。应当记住，当你设定第一层或者最后一层的厚度时，将会在距离第一个和最后一个界面一定距离处观察到相位的变化。通常，就像在椭圆偏振测量中一样，人们对复系数的比率感兴趣，该比率不会随着距离界面的距离而变化。如果你想要将单个界面的振幅系数和一个参考对比，那么需要将第一层设置非常薄来避免额外的相位移。

偶极子发生的叠层结构（STACK for Dipole Emission）

        对于许多薄膜应用比如OLEDs和其它电致发光器件，量子效率和提取效率均受到层间的界面干涉影响。Stackdipole（https://optics.ansys.com/hc/en-us/articles/360034926493）和stackpurcell（https://optics.ansys.com/hc/en-us/articles/360034926513）脚本命令非常适合研究这种围墙效应，因为它根据偶极源在堆栈中的位置来解析计算偶极源的发射特性。

       Stackdipole函数返回亮度(cd/m^2)和辐射度(W/steradian/m^2)，二者随着发射角而变化，此外还会返回相应的三刺激值X，Y，Z。CIE1931色函数用来计算三刺激值X，Y，Z。为了确定堆叠几何结构如何分别影响量子效率和提取效率，应使用stackpurcell函数。

(1) 厚层薄膜中的相干效应（Coherence Effects in Thick Layers）

当研究微腔中量子效率和提取效率时，必须对多个反射间的相干作用进行建模；然而，当一层厚度在相干长度级别时，这些相干效应应该禁用。可以使用incoherent_propagation选项来完成这个设置。在许多发射器件中，活性区域会被封装在一个较厚的玻璃保护层中。为了模拟谱线，和依赖于角度的耦合出玻璃的提取效率，应该在该区域（较厚的玻璃层中）采用非相干传播。

参考：

1. K. A. Neyts, "Simulation of light emission from thin-film microcavities", J. Opt. Soc. Am. A, Vol. 15, No. 4 (1998)
2. W. Lukosz and R. E. Kunz, "Light emission by magnetic and electric dipoles close to a plane interface", J. Opt. Soc. Am. 67, 1607-1619 (1977)
3. CIE Proceedings (1932), 1931. Cambridge: Cambridge University Press.

相关话题

See Video: Optimize Your OLED/LED in One Second for examples on how to use the Stack optical solver to simulate OLEDs and LEDs.

See stackdipole for an example of a multilayer stack dipole emission calculation.

For mode analysis of an OLED layer structure, see OLED slab mode analysis.
See OLED application examples for more examples.