-
UID:317649
-
- 注册时间2020-06-19
- 最后登录2024-05-09
- 在线时间1260小时
-
-
访问TA的空间加好友用道具
|
SITE ZHANG,1,2,*CHRISTIAN HELLMANN,2 AND FRANK WYROWSKI1 )!'Fa_$ e A*G
)CG
1Applied Computational Optics Group, Institute of Applied Physics, Friedrich Schiller University Jena, Max-Wien-Platz 1, 07743 Jena, Germany eo>/ 2Wyrowski Photonics UG, Kahlaische Straße 4, 07745 Jena, Germany fui4@ *Corresponding author: site.zhang@uni‑jena.de <WZ1- Oqd"0Qt- Received 8 March 2017; revised 29 April 2017; accepted 3 May 2017; posted 3 May 2017 (Doc. ID 290298); published 19 May 2017 |De!ti &oon'q5; Received 8 March 2017; revised 29 April 2017; accepted 3 May 2017; posted 3 May 2017 (Doc. ID 290298); published 19 May 2017 IOSuaLH^ LvSP #$f 通过使用平面波谱分析以及S矩阵方法,我们研究了一般电磁场经过由各向同性介质或者各向异性介质构成的光学层状结构的传播。我们还开发了一种基于快速傅里叶变换技术的算法,具有数值高效的采样规则。通过将此算法与其它系统建模技术相结合,我们展示了一些仿真案例,如经过一个各向同性法珀标准具的光场传输以及具有任意方位和光轴方向的单轴晶体平板的光场传输。 4 *Bp BM|-GErE OCIS codes: (260.0260) Physical optics; (260.2110) Electromagnetic optics; (260.1440) Birefringence; (230.4170) Multilayers. NFrNm'v N$N;Sw https://doi.org/10.1364/AO.56.004566 NvQY7C ,k4z; 1.引言 )?!vJb" #]5&mKi 由平行平面构成的光学层在光学中广泛应用。层状结构可以用作许多不同情况的模型,像平板和标准具。基于这个事实,光与层状结构相互作用的主题一直引起大家的注意并且对此已经进行了大量的研究。 04I6-}6 &>&dhdTQ 在这类研究中,大多数观点都侧重于平面波,然而仅仅少数的研究使用了平面波谱方法(SPW)来考虑一般的电磁场。例如,参考文献[1-6]中研究了各向同性-各向同性的界面上,高斯光束的反射率和透射率;在参考文献[7-11]中研究了各向同性层或者平板的情况;参考文献[12-22]讨论了各向同性-各向异性界面的情况,在参考文献[23-26]中则讨论了各向异性层或者平板的情况。 [
c ~LY4: ~C}(\8g 上面所提到的许多研究都用于特定的研究主题,像[1,3,5]中研究了高斯光束全内反射的横向偏移,并且他们常常关注于具体的配置。因此,将这些方法推广到更一般的情况的可能性受到了限制。 54`bE$:+ #)W8. 在这篇文章中,我们从一个更一般的观点来考虑此问题。光学层几乎不会单独使用;相反,他们常常是一个光学系统的一部分并且和其他的元件一起使用,如图1中所示。基于此事实,我们遵循场追迹的概念[27],并使用不同的场追迹算子组合[28-32],如图1中所示,以对一个包含了层介质元件的系统进行物理光学模拟。考虑到模拟是对整个系统而不是单个元件,仿真层结构必须与系统的前后部分相连接。这要求我们传播步骤(图1中的P)进行适当的考虑,将前一个元件的输出连接到当前元件的输入,并将当前元件的输出传递到下一个元件。一般情况下,这样的传输步骤会出现在平行或者非平行平面之间。在参考文献[28,29]中已经提到了平行平面间几种有效的传输方法,在参考文献[33]中则可以找到对非平行平面间传输的一个详细的讨论。在这篇文章中,我们不会研究传输步骤,但会关注层状结构的元件算子C。 Kj~>&WU \e`6=Q% 此外,从数值计算的观点出发,为了执行一个连续且有效的系统模拟,要求元件算子C 7xLo4 se<i5JsSV 正确地处理采样场数据并和其他的算子以一种统一的格式传递场数据; c%q}"Y0oh U@o2gjGN 优化数值计算的效率。 nu\ -_y~rx
> 考虑到上述两个标准,我们开发了一种具有自动数值采样规则的SPW方法。与之前一些利用积分方法对空间和角谱相关的傅里叶变换进评估的研究相比(如参考文献[23]中的二维中点规则和参考文献[12-14,20,25]中的Stamnes–Spjelkavik–Pedersen方法[34]),我们使用了快速傅里叶变换(FFT)技术,此技术在大部分数值软件包中容易访问并且效率高。再加上在角谱域中经过深入考虑的数值采样规则,我们的方法具有一般适用性,对层元件和入射场没有任何限制。因此,此算法可以直接包含在一个物理光学系统模拟之中。 I
ybl;u d[s;a. 图1.结合使用不同的场追踪算子来模拟光学系统: C是元件算子,P是相邻元件之间的传输算子。 <7vI h0 2.理论 ki[;ZmQqY y8<lp+ 如图2所示,层状结构分别由两个位于和的平行平面构成。和的区域充满了复折射率为和的均匀各向同性介质。参考文献[27]中表明使用横向分量Ex和Ey已足够表征均匀各向同性介质中电磁场了。因此,我们可以使用以下表达式来描述此问题: NYSj^k;^(z ^BF@j4*~ 其中,分别在平面和处定义输入和输出横向电场矢量,(两者位于界面的数学位置,但总是认为在均匀介质的一侧),由下式给出 ^
I{R[O'8 M!xm1-,[ 其中 。方程(1)中的元件算子是一个2x2的矩阵形式, OhSt6&+ 4f@havFIJ gE-w]/1zD5 图2.层状结构分别由两个位于 和 的平行平面构成。 和 的区域由均匀各向同性介质填充,其折射率分别是 和 。输出场和输出场在层表面进行定义,但总是在相应的各向同性介质的一侧。 }Q{
=:X9 在这章节,我们的目标是找到C的精确的形式,以连接层介质元件的输入和输出场。为了研究与层结构的相互作用,我们对输入横向场分量进行了一个傅里叶变换,并获得了 adRNrt*! "Cb.cO$i; 其中, F表示二维傅里叶变换, @&~OB/7B: q;1VF;<"vH DuTlYXM2^ M^|"be~{' 。逆傅里叶变换定义如下 +Op%,,Db *3w/`R<\ 方程(6)中的积分可以解释为将分解为具有不同横向波矢分量κ的平面波。因此,在我们的情况下,每个输入平面波都可以单独处理——我们首先计算每个输入平面波的输出,然后进行求和从而获得输出场。 NPB ,q& Th 此外,根据边界条件对电磁场施加的连续性要求,可以显示出一个给定的输入平面波在与层结构相互作用的过程中其横向波矢分量κ必定保持不变。同样可以显示出,通过叠加原理的有效性,不同的κ之间没有耦合。因此,对于输出角谱,我们可以写下 |1_$!
p F8|5_214' 其中 vOvxQS}dBp ^Ff~j&L@{ )rs|=M=Xk !#@4xeBPo Y\z^\k 公式(8)中分别是透射和反射系数矩阵。为了计算T(κ) 或者R(κ),我们选择使用数值稳定S矩阵方法。为了计算S矩阵,首先必须确定每个各向异性层的平面波。基于文献[35]中Berreman的4x4矩阵公式,Landry和Maldonado开发并展示了一种数值友好的形式,见参考文献[23]。我们采用了他们的方法,对于每一层,求解了参考文献[23]中由方程(28)所描述的特征系统的特征值和特征向量。 =`7#^7Q9 不同于[23,25]中直接使用本征解来构建一个转换矩阵,另外,我们还需要根据他们的传输方向整理出平面波,这是为计算S矩阵所做的一个必需的准备。为此,我们遵循[36]中4.3部分由Li所提出的标准。 C*W.9 然后可以应用递归S矩阵公式。我们在这篇文中不再重复给出已发展成熟的S矩阵方法,读者可以参考文献[37]中的方程(5)-(8)以获得更多的信息。在我们的情况中,由于没有反向传输输入场,我们仅对正向透射或者反向反射感兴趣,因此这篇文章中的矩阵系数T(κ) 和R(κ)对应于[37]中方程(5c)或者(5d)中的子矩阵T_uu 或者R_du。 #VQGN2bK. 一旦获得了矩阵系数,通过方程(7)即可获得输出角谱。对输出角谱进行一个逆傅里叶变换,我们获得了输出横向场矢量 [T}%q"< ;qT!fuN; JR?
)SGB ]AFM Y<mB 通过联合方程(4),(7),(9),我们可以写出从输入场到输出场整个计算流程,如果 s60:0 > 我们以透射的情况作为例子,则 ;EE&~&*w O5G<O(,\ K-"HcHuF ^ RcIE ( 因此,方程(1)中元件算子C的精确形式如下 ])$."g rHC+nou <WHs
S2At$47v 通过使用系数矩阵R(κ)代替 T(κ),可以获得反射情况的表达式。 V:0uy> H/^TXqQ8 3.算法 rr07\; 按照方程(10)的顺序,我们可以应用一种数值算法以计算场经过分层介质元件的传输。让我们从横向输入场矢量和,以图3(a)中的均匀网格进行采样。这种网格定义为,其中和作为索引数,δx和δy是x方向和y方向的采样距离。初始采样参数应该受到合适的控制以使他们符合先前算子的奈奎斯特-香农采样定理。方程(10)中的F为联系空间域和频率域的算子,可以使用不同的数值方法实现,像广泛使用的快速傅里叶变换(FFT)技术,以及包含了能够进一步提高数值效率的半解析傅里叶变换[38]和啁啾z变换[39-41]的更高级的方法。在此篇文章中,我们使用了FFT技术,并以此获得输入角谱。但我们的算法不受限于该技术。 +}`O^#<qLX | `)V^e_ n&3iv^ 图3.在角谱域定义均匀采样网格。(a),初始网格定义,5x5采样点,和为采样距离;(b),使用5x5采样点,沿垂直方向生成一个测试网格,和=0.5为采样距离;(c)中使用9x5个采样点定义了沿水平方向的一个测试网格,=0.5和为采样距离。((b)和(c)中的实心点是在(a)的初始网格中出现的共同的采样点位置,而空心圆环与初始网格中不一样的采样点。 VtzyB 9GtVI^] (8@hF#N1 f)AW !/ [)J49 >DL-Q\U 其中角谱域中的均匀网格定义为,其中和为Kx和Ky的采样距离。利用FFT计算的结果,在两个域中的采样点是一样的,因此我们有 。 (Q h7bfd 接下来,将输入角谱乘以系数矩阵,然后我们可以获得输出角谱场,例如,对于透射 cmr6,3_ "zeJ4f WDH[kJ &,Zz 角谱的采样通过傅里叶变换关系自动确保。然而,的采样没有必要进行确保,因为他们包含了系数和输入角谱之间的点乘。这个问题将在第四部分的案例中进行清楚地说明。 x^='pEt{ 一般来说,采样距离,和带宽必须进行适当调整以确保正确的采样。在我们的情况中,方程(13)中的运算没有改变谱范围。因此,我们仅仅只需要找到合适的采样距离即可。此问题并没有解析解,因为系数矩阵T或者R是通过递归S矩阵方法数值上获得的。确保采样的唯一方法是实行数值测试。 :,}:c%-^" FkxhEat8 让我们在角谱域定义一个测试采样网格,如如3(b)或者3(c)所示。测试网格定义为,其中为指数;为了不改变角谱范围,要求。与初始网格相比,它也需要进行细化,这意味着。
'QekQ]; 一方面,我们以一种严谨的方式来计算测试网格上的输出横向角谱分量。由于输入横向角谱分量的合适采样,根据需要可以进行插值,例如,在测试网格上,我们获得 u\1Wkxj iu6WGmR =-Q Z#TgFQ3u 其中“INTRPL”代表的是数值插值运算。为了在一个更精细的网格上严格地获得输入场,如方程(14)中,我们总是采用基于FFT的Sinc插值方法,把它们代入到方程(13)中,可以严格地获得输出横向角谱量 @k:f}-t Qq FfR# q/w6sQx$ ko+M,kjwR 另一方面,通过对的插值,我们也可以获得测试网格上的输出横向角谱量,并将插值结果表示为。我们将此过程描述为 8O.:3%D~
t #/B~G.+( tN'- qdm E/L?D 值得强调的是,在方程(16)中无需再使用严格的Sinc插值。因为这些结果仅仅用于评估一下情况。在这篇文章中,我们使用三次插值。上面的插值没必要给出正确的输出值,因此,在方程(16)的左边,我们对那些量使用上标“intrpl”。下一步,我们定义 CK!pH{n+ A^r
[_dyZ /nMqEHCyg R0YC:rAt 作为插值结果和严格结果之间的相对偏差。只要在初始网格上的采样满足奈奎斯特-香农采样定理,基于他们的插值结果就不会对严格的插值方法表现出一个大的偏差。在这种情况下,根据所选的插值方法,两种结果的偏差应该在数值误差级次,由σ_0表示。使用上述标准,通过一步步的减小采样间距来测试场采样,直至σ<σ_0。 \")YKN=W e/HX,sf_g 以上的测试过程对应着一个循环,在每个测试循环中都会执行方程(14)到(17)的计算。为了在每个测试循环中充分地使用测试数据,如方程(15)中的结果,我们总是使用图3(b)和3(c)中的测试网格。这样规则的测试网格定义为 k+Ma_H` C1Pt3 s!o<Pd yJK y/sWy1P7 对每个循环内,需要对先前的采样距离减半。使用这个测试网格具有如下好处:当计算时,仅仅需要计算那些在空心圆圈位置处的值,而实心点位置处的值在之前的测试循环中已经计算过了,并且让这些值简单地进行传输即可。在这种方式下,这些基于测试目的所进行严格计算的值不会被丢掉,而会作为下一次循环中的起始点以重新使用。因此,在一个完整的计算中所执行的是严格的计算,从而验证合适的采样,并且这些计算的值能够有效地用于构建最终输出场。我们将上面的方法总结在算法1中。 9J;H.:WH 算法1:通过合适的采样控制来进行输出角谱计算的数值流程 MlC-Aad( 1)在初始网格严格地计算 ++{+
#s6 2)初始化相对误差值σ=+∞ _7j/[ 3)判断σ>σ_0, ⊳如果为真,则采样不合适 zk 'e6 4)沿Kx(或者Ky)二等分采样距离 kDpZnXP 5)应用新的采样距离定义测试网格 ry'^1~, 6)在测试网格中对输入角谱插值,根据方程(14)获得 IrwF
B 7)根据方程(15),在测试网格上严格地计算输出角谱 y1"^S 8)根据方程(16),通过插值获得 >A'!T'"~ 9)根据方程(17),计算误差σ *}iT6OJ 10)如果σ>σ_(0 ) ,则⊳将目前的输出场设置为下一次循环的起始点 D16;6K'{ 11)返回 Ko]h r +zdq+<9X 为了有效地处理非对称情况,例如,光束在x方向和y方向有不同的发散角,算法1中的测试需要按顺序沿两个方向进行。开始方向的选择是任意的,在我们的情况中,我们是沿y方向开始测试的。 p2NB~t7Z Y)j,(9 4.示例 8s5ru) $E~Lu$| 在VirtualLab Fusion[42]软件中,我们将第三节中提出的算法实现在“可编程元件”的编程界面中。这个元件可以与VirtualLab Fusion中其它的物理光学仿真技术进行联合仿真。接下来,我们展示了四个案例:前两个主要关注元件本身并以一种严格的数值方式检查此算法;后两个案例中,元件将用于光学系统中,例如,此算法与其它仿真技术一起使用。 HG3>RcB CD<u@l,1 在进入实际的案例之前,设置方程(17)中的迭代终止标准σ_0很重要。对于在此文中所使用的三次插值,我们预先检查了它在一般情况下的表现并在我们的数值环境中找到了一个0.01的基准值。 ?e_}X3{ wPO@f~[Ji A.各向同性标准具 e/ % ; 36Z`.E>~L 第一个案例模拟了一个线性偏振高斯光束经过一个标准具的传输。通过这个案例,我们将清楚的指出第三节中所说的采样问题并描述了算法1的工作原理。标准具由熔融石英制成,两侧有多层薄膜,如图4所示。关于标准具的光学参数和结构的更多信息,请见表1。输入场为波长633nm,x方向线偏光的高斯光束。其在元件的输入平面定义为E_⊥^in,且束腰半径为(2um,2um)。在经过傅里叶变换后,我们获得其角谱,同样具有高斯轮廓,如图5所示。按照方程(10)中的操作算子序列,输入角谱将乘以透射或者反射率系数。我们仍以透射作为例子,并且对于线性偏振输入场,我们使用t_xx和t_yx乘以E ̃_x^in,以获得输出角谱分量。 eNC5' Z (_n8$3T75 \7tvNa,C 图4.由熔融石英制成,两侧有多层膜的标准具。其结构和光学参数如表1中所示。 (PrPH/$ tAu|8aL 表1 标准具的结构和光学参数 !RUo:b+ Dp6"I!L<| #/v_h6$ FivaCNA -I-&<+7v 图5 .(a)输入高斯场分量的振幅;(b)对应的角谱分量。由于输入光场为沿x方向的线性偏振光,因此仅显示Ex分量。 Sl<-)a: &fy8,} 如第三节中所指出,乘积的采样不能自动得到保证,此案例中将显示该现象。标准具由于其频率选择功能(频谱或角频率)而得到广泛的使用。在我们的案例中,角频率选择可以解释为系数txx和tyx以一种方式调制输入角谱,以使特定的角频率加强而其它的减弱。这种调制可以出现在一种非常精细的频率水平上。因此,需要使用更精细的采样以在输出角谱中解析这样一个精细的调制。为了获得需要的采样间距,我们遵循算法1,图6中显示了部分结果。 ; o@`l$O /sC[5G% 通过图6和表2,我们根据算法1中的步骤描述了工作流程,如下: nT.2jk+ }]GK@nn7 第一步:从图6(a)中所示的输入角谱开始,计算各个系数并乘上元件矩阵以生成初始化的输出角谱; |2&mvjk@H 第二步:初始化相对偏差σ=+∞; z`:^e1vG
第三步:开始测试循环; _ktSTzH0 第四&五步:将采样距离沿κ_x或者κ_y方向减半,以定义测试网格,对应的采样点如表2中所示; 7<Js'\Z 第六步:对输入角谱在测试网格上插值; Mc3h
R0 第七步:在测试网格上严格的计算输出角谱,在此例中,对应图6(e); HY| SLk/E 第八步:执行插值以获得,在此例中对应图6(b)-6(d); ! tPK"k 第九步:比较严格仿真和插值结果,并计算相对误差; IguG03:.N 第十步:对于较差的插值结果,如6(b)和6(c),表2中的0-3行,其结果是σ>σ_0,严格的结果将会传递给下一个循环并用于输入;否则,返回当前的结果。 &+{xR79+& gV44PI6h _,?H rL9 图6.算法1中不同步骤时沿κ_x方向的一维提取结果:(a)输入角谱振幅,(b)-(d),在测试循环中的插值角谱振幅以及(e)在最后的循环中严格地计算输出角谱。所有的子图中的值都缩放到相同的范围内。 0m!ZJH e v;qL?_:=c 从表2中我们也可以看到测试首先是沿y方向,之后再沿x方向执行的,如第三节的最后所提到的。从第0轮到第4轮测试,采样距离δκx并没有改变,因此采样点Mx^,保持不变;在第四轮测试后σ<σ0,沿y方向的测试终止,意味着场数据已经可以从先前一轮的结果中恢复。因此,沿y方向的采样点数是705,在第四轮中额外的704个数据仅仅是用于测试目的,对最终的输出并没有贡献。然后以一种类似的方式沿x方向开始测试并在第11轮终止,同样,场数据可以从先前一轮的结果中恢复。因此,最终输出的采样点数固定在2817x705。在表2的测试轮中,包含的总的采样点数是5633x705+45x704,数据45x704来自于沿y方向的最后测试轮。再次强调一下,为了在角谱域进行合适的采样,有必要采用如此大的采样点数。并且,除了沿x方向或者y方向最后一轮测试,用于测试目的而所有严格计算地值都用于构建最终输出场。 9/KQAc* cWy0N 表2 Etalon模拟中每个测试轮次中的采样参数和误差 4wD^?S!p K:50?r_-6 3K:Xxkk 文章中的模拟,使用的是一台Intel Core i7-4910MQ处理器,2.9GHz,32Gb物理内存的电脑。表2中所显示的示例所需的总的计算时间是118s。注意计算时间最多的部分是花在S矩阵计算上,而最后步骤的傅里叶逆变换仅花了约0.5s。 C)^\?DH '0M0F'R 从表2中也可以看出沿x方向和沿y方向的测试轮次数也不一样,即沿两个方向所需的采样间距不同。因此,我们在算法中更倾向沿每个方向分别执行采样测试。 $`{q = ghQsS|)p. 通过这个算法,我们获得了标准具的透过场,如图7所示。除了图7中透射场的尺寸远大于图5中的输入场,非零输出场分量也值得注意,尽管其振幅要远小于。Ex和Ey之间出现的偏振串扰是由于系数矩阵T或者R的非对角形式造成的。从的振幅分布中,还可以看出精细的同心圆,这是由于标准具多层结构中发生的多重反射造成的。 m/B9)JzY J;Xz'0 I%*Zj,> 图7.以任意单位给出的标准具输出场分量的振幅,(a):|E_x^out |,(b):|E_y^out | kV%y%l(6 a#GqJ?nY B.单轴晶体板 s"F,=]HQ!G 54_m{&hb 我们的方法对各向同性介质层和各向异性介质层都同样适用。在这个案例中,我们演示了聚光干涉仪的原理,其可以用于精确的测量单轴晶体的光轴倾斜角度[43]。当使用会聚单色光照射正交偏振器间的晶体时,即可产生聚光干涉。图8中展示了聚光干涉仪的简要原理。在我们的模拟中,光源是一个Ex线性偏振,波长为633nm,NA=0.25的球面波,在晶体板前20mm。光源将输入场传递到晶体板的前表面,后表面的输出场将被分析。在这个子节中,我们没有讨论晶体板元件外部区域的传播步骤。 2-2LmxLG *<5lx[:4/x H$iMP.AK 图8.使用聚光干涉仪测量光轴的角度。偏振片P1沿x方向,偏振片P2沿y方向,光轴表示为“o.a”。 J@{Bv% BU\NBvX$ 我们使用聚光干涉仪,测试了一个6mm厚的向列液晶板,其光轴与z轴成θ角。其中n_e=1.7,n_0=1.5[43],设θ=0°,2°以及5°,我们观察了晶体板后y方向的输出场,如图9所示。当θ=0°,晶体的光轴垂直于晶体板表面,同心环的原点位于中心,如图9(a)所示;当θ≠0°时,如此例中的图9(b)和 `fVA.% 图9(c),同心环的原点相应的发生了横向偏移。我们在模拟中获得的值与参考文献[43]中分析的结果一致。 +mPB?5 >'g60 R[ iv`-)UsE 图9.针对不同的θ值,在晶体板后表面的输出场分布(振幅平方)。使用沿y方向的偏振片P2分析输出场,因此仅显示了E_y分量。红色的十字标出了同心圆环的原点,为(0,0),(210um,0)以及(525um,0),分别对应子图(a),(b)和(c)。 o"CqVRR 1@@]h!>k: 为了获得图9中的透射场,算法1再次被用于确定采样距离。为了完整性,表3中给出了θ=5°时此例中采样距离的数值测试的细节。 md_s2d 在表3所显示的测试轮次中,包含了3969x1985+63x1984个采样点,对于每个采样点,计算了各向异性S矩阵,耗时416s。 n!orM5=:O >%A=b}VS 表3.光轴为5°的单轴晶体板在每个测试轮次中的采样距离和偏差 T U6EE `b$I)UUm $jL.TraV7 C.倾斜的各向异性晶体板系统 7cQw?C yE/I)GOQjs 如第一节所述,我们开发的算法是系统建模的一部分。在这个例子中,我们模拟了一个包含方解石晶体板的光学系统以研究偏振转换以及涡旋的生成[44],如图10所示。因为是针对整个系统进行建模,因此在处理过程使用了几种不同的场追迹算子[27]:偏振片使用琼斯矩阵算子,透镜被当作是理想透镜,到倾斜方解石板表面的光场传输使用参考文献[33]中的技术,方解石板中的光场传输则使用了此文中的方法。 _\]UA?0 根据参考文献[44],给出了方解石晶体折射率,其中 j+v)I= ]qQB+]WN #f@53Pxb 9{&x-ugM ~$Yuxo 其中λ单位为微米。 "q8'tN>< |