对PT对称波导的模拟重现及探讨浅析.docx

 PAGE \* MERGEFORMAT 26 PT对称波导的电磁仿真方法的探讨 应用物理 摘要：通过精心设计，光学波导亦可展现出PT对称性，然而要对其进行分析则需要一定的理论基础。本文拟尝试讨论传输矩阵法，时域有限差分法（FDTD）等数值电磁仿真针对该类问题分析的可行性。基于传输矩阵法，我们使用matlab编程试图对特殊设计的光波导问题进行模拟，结果表明该类问题应用传输矩阵法的条件不匹配，效果不理想，宜使用FDTD算法分析，而进一步的探讨则应建立在量子理论基础上。 关键词：PT对称；多层平板波导；传输矩阵法；FDTD；matlab 目录 引文 3 1.PT对称性. 3 2.电磁基础理论及传输矩阵法 4 2.1电磁基础理论 4 2.2传输矩阵法 5 3.时域有限差分法（FDTD） 9 3.1差分近似 9 3.2数值色散关系 10 3.3数值相速和数值群速 11 3.4数值稳定性 12 3.5激励源 13 4.多层平板波导和传输矩阵法模拟 13 5.模拟结果和讨论分析 17 6.PT对称的应用和展望 21 7.附录 21 8.参考文献 25 引文 众所周知，物理学规律常常表现出一定的对称性，这种对称性在数学上则表现为拉格朗日函数或哈密顿量对于某一变换的不变性。在诸多的对称性之中，1998年由华盛顿大学Bender教授创建的PT对称量子理论尤其受到关注。这里的P和T分别指代空间反演变换和时间反演变换，PT对称则是指某一系统在空间和时间联合反演变换下表现出不变性。PT对称理论主要应用于量子力学和粒子物理中，然而，在Bender教授和Boettcher教授所著论文[1]的基础上，近年来也有不少文章开始针对光学系统可能表现出的PT对称性进行探讨。其中，有一些文章提出了使用复折射率或边界条件的光学波导来赋予系统PT对称性的建议。本文即是以其中一篇文献[2]为基础，针对复现该文献所述内容使用何种电磁仿真方法进行分析，并对多层平板波导展现出的PT对称性进行了部分的探讨。 PT对称性 PT对称性，即在对某一系统的哈密顿量或拉格朗日函数做时空反演时保持不变的性质。这里首先要提一下空间反演变换和时间反演变换。 空间反演变换，顾名思义，是将某一系统内所有的坐标和坐标相关项改变符号的变换。最简单的是x变为-x，动量p变为-p。这是在经典力学下对空间反演变换的描述；在量子力学背景下，空间反演变换则是指对波函数的自变量r改变符号，即变为其镜像。在这样的情况下，若镜像变换下其波函数Ψ同未做变换时一致，则认为具有空间反演对称性。 同空间反演变换类似，时间反演变换是指将某一系统内所有的时间量改变符号的变换，即t变为-t，动量p变为-p。这是在经典力学下对时间反演变换的描述；在量子力学背景下，时间反演变换则是指对含时波函数的变量t改变符号。同样的，若改变后和改变前的波函数一致，则认为具有时间反演对称性。 按照量子力学定义，P变换（空间反演）和T变换（时间反演）可以分别用如下形式描述[3]： Pfx,t=f-x,t （seq 公式1.1） Tfx,t=fx,-t （1.2） ?f∈L2Fn×R （1.3） 很明显可以看出，两个变换都是连续的。空间反演变换P是线性的，而时间反演变换T则是共轭线性的。同样，据量子力学定义，在P，T是复希尔伯特空间K中的算子时，有： P2=T2=I，PT=TP （1.4） 由上可知，当系统满足PT对称性时，其哈密顿量H应可与PT联合算符交换位置。换而言之，有下式： PTH＝HTP=H （1.5） 根据量子力学，系统的哈密顿量可以写成如下形式： H=p22+Vx （1.6） 对其做PT联合变换，有： PTH=p22+V*-x （1.7） 按照PT对称性的要求，则有： p22+V*-x=p22+Vx （1.8） 换而言之，一个满足PT对称性的系统要求： V*-x=Vx （1.9） 电磁基础理论及传输矩阵法 2.1电磁基础理论 在探讨系统是否满足PT对称性之前，首先来回顾一下电磁基础理论。经典电磁理论是以麦克斯韦方程组为基础建立的，该方程组高度凝练了电磁场的行为和本质，在此列出[4]： ?×E=-?B?t （2.1.1） ?×H=jf+?D?t （2.1.2） ??D=ρf （2.1.3） ??B=0 （2.1.4） 这里使用了非稳恒态下的麦克斯韦方程组是由于要构造的波导结构中，存在复折射率区域。换而言之，电磁场在其中传播时存在着损耗。在这种情况下，可将介电常数化为复介电常数来保持原有平面波解形式不变，同样折射率取为复数。则在此时，平面波传播中折射率的公式 n=μ? （2.1.5） 依然保持原有形式。 在之后所要探讨的波导结构中，存在着多层周期性排布的介质区域。对于这种多层结构，使用传输矩阵法