Loading [MathJax]/jax/output/HTML-CSS/jax.js

Search

Article

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Modal analysis of planar waveguides by the immersed interface method

Xiao Yong Cao Yu-Sheng

Xiao Yong, Cao Yu-Sheng. Modal analysis of planar waveguides by the immersed interface method. Acta Phys. Sin., 2023, 72(14): 140201. doi: 10.7498/aps.72.20230595
Citation: Xiao Yong, Cao Yu-Sheng. Modal analysis of planar waveguides by the immersed interface method. Acta Phys. Sin., 2023, 72(14): 140201. doi: 10.7498/aps.72.20230595

Modal analysis of planar waveguides by the immersed interface method

Xiao Yong, Cao Yu-Sheng
Article Text (iFLYTEK Translation)
PDF
HTML
Get Citation
  • Modal analysis of optical waveguides is a basic task in the design of advanced waveguide devices and optical circuits. How to deal with the problem of electromagnetic heterogeneous interface and absorption boundary condition are two major difficulties in implementing efficient numerical analysis of optical waveguides. Existing high-order accurate finite-difference modal analysis methods do not take into consideration the absorption boundary problem, which, thus, makes it difficult to accurately simulate leakage and radiation modes. Based on the immersed interface method and perfectly matched layer absorption boundary condition, a finite-difference method with the second- and fourth-order accuracy is proposed in this work. By using this method, the single-interface plasmonic waveguide mode, planar symmetric waveguide mode, and one-dimensional photonic crystal waveguide mode are analyzed. Numerico-experimental results show that the convergence rate of the second- and fourth-order algorithm are consistent with the anticipated order of the guided mode, leakage mode and radiation mode. The second-order algorithm provides an ultimate accuracy of about 109 for the relative error of effective refractive index, when the normalized step size is 104. The fourth-order algorithm provides an ultimate accuracy of about 1010 for the relative error of effective refractive index, when the normalized step size is 103. Through the study of field distribution of guided mode and cladding mode in a one-dimensional photonic crystal waveguide, we show that the continuity of the field of transverse electric mode and its first derivative across the interface, and the continuity of the field of transverse magnetic mode and the discontinuity of its first derivative across interface, can be analyzed accurately. The method proposed in this paper can be used to calculate any mode for any refractive index profile, only by using the value of refractive index, independent of the specific functional representation of modal fields. The method provides a simple and efficient tool for implementing the modal analysis of step-index planar waveguides.
      PACS:
      Corresponding author: Cao Yu-Sheng, yscao@outlook.com
    • Funds: Project supported by the Natural Science Foundation of Ningxia, China (Grant No. 2020AAC03071) and the National Natural Science Foundation of China (Grant No. 11504190)

    电磁异质结构对于集成光学和光纤光学中多数光子器件功能的实现至关重要. 在先进波导器件和光波线路的设计中, 光波导模式的全矢量求解是一项基本任务[1-3]. 在各种数值方法中, 有限差分法(finite-difference method, FDM)易于实现而被普遍用于波导模式分析[4,5]. 有限差分法一般通过求解稀疏矩阵本征值问题可同时高效地确定多个模式的传播常数以及模场. 波导存在折射率跳变的界面, 导致模场及其导数一般是不连续的. 标准有限差分法基于微分方程的解及其导数的连续性, 直接应用于界面问题时, 会降低差分格式的精度. 针对界面问题的高精度有限差分求解, LeVeque和Li[6,7]提出了浸入界面方法(immersed interface method, IIM), 该方法可以系统地生成所需精度的差分格式. Horikis和Kath[8-10]将IIM应用于光波导模式分析, 分别发展了具有任意截面形状二维波导的二阶精度算法以及光纤的二阶和四阶精度算法. 在IIM提出之前, Stern[11]利用泰勒级数展开模场然后匹配界面边界条件的方法, 来改进阶跃折射率光波导模式分析的有限差分法. 模场展开和匹配边界条件也是IIM的两个基本步骤, IIM和Stern方法仅在差分格式的系数修正如何表示方面存在微小差异, IIM将系数修正表示为矩阵方程的解, Stern方法直接给出系数修正的解析表达式. 在Stern工作的基础上, 研究者经过二三十年的努力, 最终开发出了光波导的全矢量任意阶精度模式求解器[12-18]. 高阶格式提供更快的收敛速度, 极大降低计算资源需求, 提高计算精度. 然而在现有波导模式分析的高阶方法中, 均未考虑吸收边界条件, 需要设置较大的计算区域以令导模的倏逝场衰减至可用硬壁截断, 而漏模则难以精确模拟.

    本文针对阶跃折射率平面波导模式的全矢量求解, 采用IIM和完美匹配层(perfectly matched layer, PML)[19,20]吸收边界条件发展二阶和四阶精度的有限差分法, 该方法对导模、漏模和辐射模均能高效模拟. 本文的组织结构为: 第2节列出发展算法所需的解析方程; 第3节详细推导二阶精度有限差分格式, 四阶精度格式的推导与二阶精度情形类似, 主要结果罗列在附录中; 第4节通过多种典型平面波导中各类模式的模拟, 验证算法的有效性; 第5节给出结论.

    考虑如图1所示的多层平行平面系统, 每层填充各向同性介质. 平行平面系统的模式分为横电(transverse electric, TE)和横磁(transverse magnetic, TM)偏振, 非零场分量分别为(Ey,Hz,Hx)(Hy,Ez,Ex). 取时谐因子为eiωt, 并令/y=0, TE偏振满足的场方程为

    图 1 平面波导示意图及坐标系放置\r\nFig. 1. Depiction of a planar waveguide and layout of the coordinate system.
    图 1  平面波导示意图及坐标系放置
    Fig. 1.  Depiction of a planar waveguide and layout of the coordinate system.
    Eyz=iωμHx, (1)
    Eyx=iωμHz, (2)
    (2z2+2x2+ω2εμ)Ey=0. (3)

    边界条件为EyHz在界面处连续. TE和TM偏振有对偶关系, 二者满足的场方程和边界条件在变换

    EH,HE,εμ,με (4)

    下均互换. 因此, 可以先求解TE偏振, TM偏振的相应结果直接由对偶关系获得.

    令导波模式的传播因子为eiβz, 简记EyE, iη0HzHk0xˉx, 其中η0=μ0/ε0为真空波阻抗、k0=ωε0μ0为真空波数, (2)式和(3)式归一化为

    Eˉx=μrH, (5)
    2Eˉx2+n2E=n2effE, (6)

    其中n=εrμr为介质的折射率, εrμr分别为相对介电常数和相对磁导率, neff=β/k0为模式的有效折射率. 边界条件为EH在界面处连续, 导波模式的有效折射率及电场分布由(6)式描述的本征问题解出.

    值得指出的是, 经典电磁学中没有长度量纲的基本常数, 麦克斯韦方程满足缩放不变性[21]. 这种缩放不变性在本问题中表现为, (5)式和(6)式中的空间坐标x和电磁波长λ均可收缩入归一化坐标ˉx=k0x=2πx/λ. 在将空间坐标归一化后, 计算结果对于波导结构和波长满足相同比例关系的所有情形均适用.

    对于给定的平面波导结构, 选取包含所有界面的计算区域, 并用N个等间隔的格点将其离散, 电场和磁场位于相同的格点, 如图2所示, 其中竖线表示界面. 记归一化空间步长为h, 格点i的归一化坐标为ˉxi=ih. 对于一般的第i个格点, 用三点格式

    图 2 计算区域离散\r\nFig. 2. Discretization of the computational domain.
    图 2  计算区域离散
    Fig. 2.  Discretization of the computational domain.
    C1Ei1+C0Ei+C+1Ei+1=n2effEi (7)

    近似波动方程(6).

    若格点i及其邻近的两个格点i1i+1均在同一介质层内, 则i为常规格点, 对二阶导数2/ˉx2取标准的中心差分格式近似将获得二阶精度, 相应的差分方程(7)的系数为

    C1=1h2,C0=n22h2,C+1=1h2. (8)

    对于如图2所示的某一界面α两侧的格点jj+1, 三点格式跨越界面, 由于E的一阶导数不连续, 若依然采用(8)式作为差分方程(7)的系数, 则全局精度一般降为零阶. 格点jj+1称为非常规格点, 相应的系数C±1,0需要修正以获得全局的二阶精度.

    采用IIM获得非常规格点处的修正系数C±1,0, 需要E及其一阶和二阶导数在界面处的跳跃条件. 引入记号[f]=0表示函数f在界面处连续, 即f=f+, 其中“”和“+”分别表示界面左侧和右侧. 边界条件直接给出E的跳跃条件:

    [E]=0. (9)

    H的边界条件及(5)式导出E/ˉx的跳跃条件:

    [1μrEˉx]=0. (10)

    E的边界条件及(6)式导出2E/ˉx2的跳跃条件:

    [2Eˉx2+n2E]=0. (11)

    Ej1, EjEj+1关于ˉxα=αh泰勒展开得

    Ej1=E+(ˉxj1ˉxα)Eˉx+12(ˉxj1ˉxα)22Eˉx2+O(h3), (12)
    Ej=E+(ˉxjˉxα)Eˉx+12(ˉxjˉxα)22Eˉx2+O(h3), (13)
    Ej+1=E++(ˉxj+1ˉxα)E+ˉx+12(ˉxj+1ˉxα)22E+ˉx2+O(h3). (14)

    利用跳跃条件(9)—(11)式, 将(14)式改写为

    Ej+1=δ+(ˉxj+1)E+μ+(ˉxj+1ˉxα)Eˉx+12(ˉxj+1ˉxα)22Eˉx2+O(h3), (15)

    其中

    δ+(ˉx)=1+12(ˉxˉxα)2n+, (16)
    n+=(n)2(n+)2, (17)
    μ+=μ+rμr. (18)

    将波动方程(6)应用至界面左侧得

    2Eˉx2+(n)2E=n2effE. (19)

    在格点j, 差分方程(7)的局域截断误差估计为

    T=C1Ej1+C0Ej+C+1Ej+1n2effE+O(h). (20)

    将(12)式、(13)式、(15)式和(19)式代入(20)式, 并令(20)式等号右边O(h)前的部分为零, 得到

    [11δ+(ˉxj+1)ˉxj1ˉxαˉxjˉxα(ˉxj+1ˉxα)μ+(ˉxj1ˉxα)2(ˉxjˉxα)2(ˉxj+1ˉxα)2]×[C1C0C+1]=[(n)202], (21)

    从中求得的C±1,0使得格点j处的局域截断误差为T=O(h), 即局域一阶精度. 由于界面比计算区域低一维, 非常规格点具有局域一阶精度即可以保证算法的全局二阶精度[6]. 引入记号

    V(ˉx)=[1ˉxˉxα(ˉxˉxα)2],C=[C1C0C+1],b(n)=[n202], (22)

    将(21)式改写为紧凑形式

    [V(ˉxj1),V(ˉxj),diag([δ+(ˉxj+1),μ+,1])×V(ˉxj+1)]C=b(n). (23)

    类似地, 可以求得非常规格点j+1处的修正系数C±1,0满足的方程:

    [diag([δ+(ˉxj),μ+,1])V(ˉxj),V(ˉxj+1),V(ˉxj+2)]C=b(n+), (24)

    δ+μ+的定义式(16)—(18)式中所有上标“–”替换为“+”、“+” 替换为“–”得到δ+μ+.

    在上述非常规格点差分格式系数修正的推导中, 未限制电磁参数是否为实数或复数(介质可具有增益或损耗), 亦未限制其在界面两侧的跳跃量. 因此对于任意电磁参数剖面的平面波导, 本文描述的方法均适用.

    波导的第一层和最后一层为开放介质层, 需要做电磁波的吸收处理以降低计算区域截断带来的误差. PML是目前计算电磁学中较为普遍采用的吸收边界条件. Berenger[19]最初通过引入非物理分裂场的方式实现PML, 而利用麦克斯韦方程在复平面上的解析延拓性质, 通过复坐标拖拽实现PML的等价方式[20]在数学上更加优美且易于扩展至其他应用场景.

    使用PML吸收边界, 在吸收层内做复坐标拖拽ˉx˜x(ˉx), 波动方程(6)变为

    2E˜x2+n2E=n2effE, (25)

    其中˜xˉx的复函数, 具体函数形式由吸波需求决定, ˜x的实部用于吸收倏逝波, 虚部用于吸收行波. 可以推得

    2˜x2=1˜χ22ˉx2˜γ˜χ3ˉx, (26)

    其中˜χ=d˜x/dˉx, ˜γ=d˜χ/dˉx. 将(26)式代入(25)式并取标准的中心差分格式近似, 得PML中差分方程(7)的系数为

    C1=1˜χ2h(1h+˜γ2˜χ),C0=n22˜χ2h2,C+1=1˜χ2h(1h˜γ2˜χ). (27)

    电磁波被有效吸收后, 在第一个格点和最后一个格点外应用硬壁条件E=0, 最终得到N×N的三对角矩阵本征值问题, 利用稀疏矩阵技术可以高效地求得期望模式的有效折射率和电场分布.

    算法可以自然地扩展至高阶精度, 附录给出了四阶精度的相应公式. 对TE偏振的算法做对偶变换((4)式), 即得TM偏振的相应结果.

    通过典型平面波导, 如单界面等离子波导、对称波导和一维光子晶体波导的模式分析, 验证算法的有效性. 需要指出的是, 本节算例中涉及的坐标、厚度和步长均为相应量乘以k0=2π/λ归一化后的数值. 简单起见, 仅考虑模式在最外层主要表现为倏逝波的情形. 如图3所示, PML吸收边界可通过如下的实坐标拖拽实现:

    图 3 PML吸收边界条件设置示意图\r\nFig. 3. Schematic diagram of the setup of PML absorption boundary conditions.
    图 3  PML吸收边界条件设置示意图
    Fig. 3.  Schematic diagram of the setup of PML absorption boundary conditions.
    ˜x(ˉx)=ˉx+(ˉxeˉxsδ)(ˉxˉxsδ)m, (28)

    其中ˉxs为拖拽起点, 可取为最外层第一个常规格点所在位置; ˉxe描述拖拽程度, δ为吸收层厚度; m需不小于格式精度, 统一取为4. ˉxe由下式设定:

    ˉxe=ˉxα±logαPML(n2eff,PMLεrμr), (29)

    其中“±”的“+”和“–”分别对应左侧和右侧最外层, ˉxα为相应的物理界面位置, εrμr为层中电磁参数; αPML为坐标拖拽后最外侧场相对界面处场的比值, 统一取为108; neff,PML为待求模式有效折射率的估计值, 或者为使得所有期望模式有效吸收的数值, 例如取为有效折射率的最小估计值; 表示求复数的实部, 用于处理漏模和辐射模. 算例采用通信波长1.55 μm处典型电磁材料的参数, 如εAu=104.2+3.7i, εair=1, εSi=12.25, εSiO2=1.96, 材料的相对磁导率均取为1.

    金属和介质构成的单界面表面等离子波导中, TM模式的有效折射率可由解析式给出:

    neff=εdεmεd+εm. (30)

    对于金属为金、介质为空气的情形, 精确解为neff=1.00482710586+0.00017264861i. 该波导中金属层的折射率为nAu=0.18+10.21i, 其含有较大的虚部, 且与空气层的折射率差异极大. 三层平面对称波导中, TE模式的有效折射率由解析的特征方程

    (μr2μr1α1α2μr1μr2α2α1)sin(α1d)=2cos(α1d) (31)

    定出, 其中α1=n21n2eff, α2=n2effn22, 12分别代表芯层和包层; d为归一化芯层厚度. 对于芯层为硅, 包层为空气, d=1的情形, 前两个模式有效折射率的精确解为neff,1=2.92535519956neff,2=1.05265908179.

    图4给出有效折射率的相对误差

    图 4 收敛曲线, 单界面波导($ \text{TM}_\text{s} $)和对称波导($ \text{TE}_\text{1, 2} $) 模式\r\nFig. 4. Convergence curves: Modes of single interface waveguide ($ \text{TM}_\text{s} $) and symmetric waveguide ($ \text{TE}_\text{1, 2} $).
    图 4  收敛曲线, 单界面波导(TMs)和对称波导(TE1, 2) 模式
    Fig. 4.  Convergence curves: Modes of single interface waveguide (TMs) and symmetric waveguide (TE1, 2).
    RE=|neff,numericalneff,analytical||neff,analytical| (32)

    与计算步长h的关系, 其中所有材料层的归一化厚度均为1, 单界面波导的neff,PML=1.004, 对称波导的neff,PML=1.05. 作为对比, 图中的虚线为理想的二阶和四阶收敛曲线. 可以看出, 对于三个模式, 二阶和四阶算法的收敛速度均与期望的阶数相符. 计算误差来源包含舍入误差和截断误差. 舍入误差仅与计算机硬件(本文计算在CPU主频3.40 GHz 和双精度下执行)和计算操作次数有关, 步长越小则操作次数越多, 误差累积越大. 对于稳定的有限差分算法, 截断误差随步长减小而减小. 在大步长和小步长区间, 计算误差分别由截断误差和舍入误差主导. 两种误差相等的步长即对应最优步长, 计算误差降为极限值[22]. 二阶算法约在步长h=104时提供约为RE=109的极限精度; 四阶算法约在步长h=103时提供约为RE=1010 的极限精度; 进一步减小步长将使舍入误差超过模型截断误差. 对于给定的步长, 四阶方法相对二阶方法在数量级上降低误差, 而计算耗时至多仅增加约(53)/3 67%. 例如, h=9.3750×104时, TE1模式在二阶方法下RE=6.7717×107, 计算耗时0.181 s; 在四阶方法下RE=4.1112×1011, 计算耗时0.244 s, 计算耗时增加约35%.

    常规有限差分法在处理界面问题时通常使用参数平滑或者平均方法. 对于椭圆微分方程界面问题, 不连续系数的调和平均方法比平滑方法更为精确, 然而使用调和平均的有限差分法仅对满足特定条件的一维椭圆界面问题提供二阶精度[23]. 图5给出电磁参数未做平均处理和施加调和平均处理两种条件下, 对称波导中TE1模式和单界面波导中TMs模式的收敛曲线. TE1模式的收敛速度在两种处理下无论对于二阶还是四阶格式均降为一阶, 未做平均处理时收敛曲线抖动剧烈, 调和平均处理后收敛曲线趋于平稳, 但误差增加约1个数量级. 对于TMs模式, 两种处理条件下二阶和四阶格式的精度均降为零阶, 且误差水平基本一致, 其中调和平均处理下的收敛曲线未在图5中给出.

    图 5 收敛曲线, 电磁参数未做处理(-non)和施加调和平均处理(-avg)\r\nFig. 5. Convergence curves: non-averaging (-non) and harmonic averaging (-avg) for electromagnetic parameters.
    图 5  收敛曲线, 电磁参数未做处理(-non)和施加调和平均处理(-avg)
    Fig. 5.  Convergence curves: non-averaging (-non) and harmonic averaging (-avg) for electromagnetic parameters.

    图5结果可解释如下. 对于TE模式, 三点式((7)式)仅在(9)—(11)式给出的电场跳跃条件全部满足的情形下提供二阶精度收敛. 在介质非磁时, 电磁参数未做处理使电场及其一阶导数的跳跃条件满足, 而二阶导数跳跃条件不满足, 因而给出一阶收敛速度. 此时, 五点式((A1)式)亦仅给出一阶精度. 电磁参数调和平均处理后, 尽管收敛曲线的平稳性得到改善, 但由于电场一阶导数的跳跃条件不再满足, 致使误差上升. 对于TM模式, 磁场的一阶和二阶导数的跳跃条件在电磁参数未做处理下均不满足, 三点式((7)式)和五点式((A1)式)只能给出零阶精度. 本文利用的IIM方法精确匹配所需的跳跃条件, 对TE和TM模式都能提供期望的收敛阶.

    考虑更为复杂, 解析方法难于处理的一维光子晶体波导. 取波导堆叠为A(LH)MC(HL)MA, 其中低折射率层L为二氧化硅, 高折射率层H为硅, 芯层C的介电常数εC=4, 周围介质A为空气. 算例取M=10, 归一化层宽wC=5, wL, H由1/4波长堆叠条件给出:

    wL=π2εCεL,wH=π2εHεC. (33)

    定义当前步长h下有效折射率的相对误差为

    RE=|neff,hneff,h1||neff,h1|, (34)

    其中h1为小于h的临近步长. 计算取wA=1, neff, PML=1.5, 有效折射率的相对误差与步长的关系如图6 所示. 算法对包含23层介质的一维光子晶体波导同样有效, 四阶格式同样约在步长h=103时给出约为RE=1010的精度.

    图 6 收敛曲线: 一维光子晶体波导模式($ \text{TE}_\text{c} $为包层模式)\r\nFig. 6. Convergence curves: Modes of one-dimensional photonic crystal waveguide ($ \text{TE}_\text{c} $ is cladding mode).
    图 6  收敛曲线: 一维光子晶体波导模式(TEc为包层模式)
    Fig. 6.  Convergence curves: Modes of one-dimensional photonic crystal waveguide (TEc is cladding mode).

    为了清楚地展示模场分布, 令M=3. 图7给出了3个模式在对称轴左侧的场分布, 其中取了二阶算法, 步长h=1.688×102, 其他计算参数不变. 由图可见, 界面处TE模式的场及其一阶导数的连续, TM模式的场的连续及其一阶导数的不连续, 均可以被精确解析.

    图 7 模场分布: 一维光子晶体波导, 竖直虚线标示界面位置\r\nFig. 7. Distribution of modal fields: One dimensional photonic crystal waveguide, with dashed lines indicating the positions of interfaces.
    图 7  模场分布: 一维光子晶体波导, 竖直虚线标示界面位置
    Fig. 7.  Distribution of modal fields: One dimensional photonic crystal waveguide, with dashed lines indicating the positions of interfaces.

    针对平面波导的高效模式分析问题, 本文基于浸入界面方法和完美匹配层吸收边界条件, 提出一种具有二阶和四阶精度的全矢量计算方法. 本方法只需利用取样点上折射率的值而不依赖于模场的特定函数展开形式, 对于任意折射率剖面下的任意模式(导模、漏模和辐射模)均能简单而高效求解. 算法不仅可进一步向高阶扩展, 且其技术框架亦可应用于其他波动系统的数值模拟中界面问题和吸收边界条件的处理.

    本附录罗列TE偏振模式的四阶精度公式. 对于一般的第i个格点, 用五点格式

    C2Ei2+C1Ei1+C0Ei+C+1Ei+1+C+2Ei+2=n2effEi

    近似波动方程(6).

    电场E的三阶和四阶导数的跳跃条件为

    [1μr(3Eˉx3+n2Eˉx)]=0,
    [4Eˉx4+2n22Eˉx2+n4E]=0.

    引入中间函数

    ϕ+(ˉx)=δ+(ˉx)+(ˉxˉxα)424n+,
    κ+(ˉx)=μ+θ+(ˉx),
    θ+(ˉx)=13[2+δ+(ˉx)],

    将(22)式扩展为

    V(ˉx)=[1ˉxˉxα(ˉxˉxα)2(ˉxˉxα)3(ˉxˉxα)4],C=[C2C1C0C+1C+2],b(ˉx,n)=[n2(ˉxˉxα)n22+(ˉxˉxα)2n26(ˉxˉxα)12(ˉxˉxα)2],

    非常规格点j处, 浸没界面修正的差分方程系数满足

    [V(ˉxj2),V(ˉxj1),V(ˉxj),diag([ϕ+(ˉxj+1),κ+(ˉxj+1),θ+(ˉxj+1),μ+,1])V(ˉxj+1),  diag([ϕ+(ˉxj+2),κ+(ˉxj+2),θ+(ˉxj+2),μ+,1])V(ˉxj+2)]C=b(ˉxj,n).

    非常规格点j1处为

    [V(ˉxj3),V(ˉxj2),V(ˉxj1),V(ˉxj),  diag([ϕ+(ˉxj+1),κ+(ˉxj+1),θ+(ˉxj+1),μ+,1])V(ˉxj+1)]C=b(ˉxj1,n).

    非常规格点j+1处为

    [diag([ϕ+(ˉxj1),κ+(ˉxj1),θ+(ˉxj1),μ+,1])V(ˉxj1),  diag([ϕ+(ˉxj),κ+(ˉxj),θ+(ˉxj),μ+,1])V(ˉxj),V(ˉxj+1),V(ˉxj+2),V(ˉxj+3)]C=b(ˉxj+1,n+).

    非常规格点j+2处为

    [diag([ϕ+(ˉxj),κ+(ˉxj),θ+(ˉxj),μ+,1])V(ˉxj),  V(ˉxj+1),V(ˉxj+2),V(ˉxj+3),V(ˉxj+4)]C=b(ˉxj+2,n+).

    利用(A7)式定义的C将常规格点处的系数紧凑地表示为

    CT=[1,16,30,16,1]12h2+[0,0,n2,0,0],

    PML中格点的系数为

    CT=[1,16,30,16,1]12h2˜χ2[1,8,0,8,1]˜γ12h˜χ3+[0,0,n2,0,0].

    [1]

    陈宪锋, 沈小明, 蒋美萍, 金铱 2008 物理学报 57 3578Google Scholar

    Chen X F, Shen X M, Jiang M P, Jin Y 2008 Acta Phys. Sin. 57 3578Google Scholar

    [2]

    陈云天, 王经纬, 陈伟锦, 徐竞 2020 物理学报 69 154206Google Scholar

    Chen Y T, Wang J W, Chen W J, Xu J 2020 Acta Phys. Sin. 69 154206Google Scholar

    [3]

    李慧慧, 薛文瑞, 李宁, 杜易达, 李昌勇 2022 物理学报 71 108101Google Scholar

    Li H H, Xue W R, Li N, Du Y D, Li C Y 2022 Acta Phys. Sin. 71 108101Google Scholar

    [4]

    乔海亮, 王玥, 陈再高, 张殿辉 2013 物理学报 62 070204Google Scholar

    Qiao H L, Wang Y, Chen Z G, Zhang D H 2013 Acta Phys. Sin. 62 070204Google Scholar

    [5]

    Wu D, Ohnishi R, Uemura R, Yamaguchi R, Ohnuki S 2018 IEEE Photonics Technol. Lett. 30 1024Google Scholar

    [6]

    LeVeque R J, Li Z L 1994 SIAM J. Numer. Anal. 31 1019Google Scholar

    [7]

    LeVeque R J, Li Z L 1995 SIAM J. Numer. Anal. 32 1704Google Scholar

    [8]

    Horikis T P, Kath W L 2006 Opt. Lett. 31 3417Google Scholar

    [9]

    Horikis T P 2011 Int. J. Mod. Phys. C 22 687Google Scholar

    [10]

    Horikis T P 2013 Appl. Math. Model. 37 5080Google Scholar

    [11]

    Stern M S 1988 IEE Proc. J. (Optoelectron.) 135 56Google Scholar

    [12]

    Vassallo C 1992 IEE Proc. J. (Optoelectron.) 139 137Google Scholar

    [13]

    Yamauchi J, Sekiguchi M, Uchiyama O, Shibayama J, Nakano H 1997 IEEE Photonics Technol. Lett. 9 961Google Scholar

    [14]

    Yamauchi J, Takahashi G, Nakano H 1998 IEEE Photonics Technol. Lett. 10 1127Google Scholar

    [15]

    Chiang Y C, Chiou Y P, Chang H C 2002 J. Lightwave Technol. 20 1609Google Scholar

    [16]

    Sujecki S 2010 Opt. Lett. 35 4115Google Scholar

    [17]

    Chiou Y P, Du C H 2011 J. Lightwave Technol. 29 3445Google Scholar

    [18]

    Cao Y 2022 Opt. Express 30 4680Google Scholar

    [19]

    Berenger J P 1994 J. Comput. Phys. 114 185Google Scholar

    [20]

    Chew W C, Weedon W H 1994 Microw. Opt. Technol. Lett. 7 599Google Scholar

    [21]

    Joannopoulos J D, Johnson S G, Winn J N, Meade R D 2008 Photonic Crystals: Molding the Flow of Light (Princeton: Princeton University Press) p20

    [22]

    Press W H, Teukolsky S A, Vetterling W T, Flannery B P 2007 Numerical Recipes: The Art of Scientific Computing (Cambridge: Camebridge University Press) p10

    [23]

    Li Z L, Ito K 2006 The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Philadephia: SIAM) p10

  • 图 1  平面波导示意图及坐标系放置

    Figure 1.  Depiction of a planar waveguide and layout of the coordinate system.

    图 2  计算区域离散

    Figure 2.  Discretization of the computational domain.

    图 3  PML吸收边界条件设置示意图

    Figure 3.  Schematic diagram of the setup of PML absorption boundary conditions.

    图 4  收敛曲线, 单界面波导(TMs)和对称波导(TE1, 2) 模式

    Figure 4.  Convergence curves: Modes of single interface waveguide (TMs) and symmetric waveguide (TE1, 2).

    图 5  收敛曲线, 电磁参数未做处理(-non)和施加调和平均处理(-avg)

    Figure 5.  Convergence curves: non-averaging (-non) and harmonic averaging (-avg) for electromagnetic parameters.

    图 6  收敛曲线: 一维光子晶体波导模式(TEc为包层模式)

    Figure 6.  Convergence curves: Modes of one-dimensional photonic crystal waveguide (TEc is cladding mode).

    图 7  模场分布: 一维光子晶体波导, 竖直虚线标示界面位置

    Figure 7.  Distribution of modal fields: One dimensional photonic crystal waveguide, with dashed lines indicating the positions of interfaces.

  • [1]

    陈宪锋, 沈小明, 蒋美萍, 金铱 2008 物理学报 57 3578Google Scholar

    Chen X F, Shen X M, Jiang M P, Jin Y 2008 Acta Phys. Sin. 57 3578Google Scholar

    [2]

    陈云天, 王经纬, 陈伟锦, 徐竞 2020 物理学报 69 154206Google Scholar

    Chen Y T, Wang J W, Chen W J, Xu J 2020 Acta Phys. Sin. 69 154206Google Scholar

    [3]

    李慧慧, 薛文瑞, 李宁, 杜易达, 李昌勇 2022 物理学报 71 108101Google Scholar

    Li H H, Xue W R, Li N, Du Y D, Li C Y 2022 Acta Phys. Sin. 71 108101Google Scholar

    [4]

    乔海亮, 王玥, 陈再高, 张殿辉 2013 物理学报 62 070204Google Scholar

    Qiao H L, Wang Y, Chen Z G, Zhang D H 2013 Acta Phys. Sin. 62 070204Google Scholar

    [5]

    Wu D, Ohnishi R, Uemura R, Yamaguchi R, Ohnuki S 2018 IEEE Photonics Technol. Lett. 30 1024Google Scholar

    [6]

    LeVeque R J, Li Z L 1994 SIAM J. Numer. Anal. 31 1019Google Scholar

    [7]

    LeVeque R J, Li Z L 1995 SIAM J. Numer. Anal. 32 1704Google Scholar

    [8]

    Horikis T P, Kath W L 2006 Opt. Lett. 31 3417Google Scholar

    [9]

    Horikis T P 2011 Int. J. Mod. Phys. C 22 687Google Scholar

    [10]

    Horikis T P 2013 Appl. Math. Model. 37 5080Google Scholar

    [11]

    Stern M S 1988 IEE Proc. J. (Optoelectron.) 135 56Google Scholar

    [12]

    Vassallo C 1992 IEE Proc. J. (Optoelectron.) 139 137Google Scholar

    [13]

    Yamauchi J, Sekiguchi M, Uchiyama O, Shibayama J, Nakano H 1997 IEEE Photonics Technol. Lett. 9 961Google Scholar

    [14]

    Yamauchi J, Takahashi G, Nakano H 1998 IEEE Photonics Technol. Lett. 10 1127Google Scholar

    [15]

    Chiang Y C, Chiou Y P, Chang H C 2002 J. Lightwave Technol. 20 1609Google Scholar

    [16]

    Sujecki S 2010 Opt. Lett. 35 4115Google Scholar

    [17]

    Chiou Y P, Du C H 2011 J. Lightwave Technol. 29 3445Google Scholar

    [18]

    Cao Y 2022 Opt. Express 30 4680Google Scholar

    [19]

    Berenger J P 1994 J. Comput. Phys. 114 185Google Scholar

    [20]

    Chew W C, Weedon W H 1994 Microw. Opt. Technol. Lett. 7 599Google Scholar

    [21]

    Joannopoulos J D, Johnson S G, Winn J N, Meade R D 2008 Photonic Crystals: Molding the Flow of Light (Princeton: Princeton University Press) p20

    [22]

    Press W H, Teukolsky S A, Vetterling W T, Flannery B P 2007 Numerical Recipes: The Art of Scientific Computing (Cambridge: Camebridge University Press) p10

    [23]

    Li Z L, Ito K 2006 The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Philadephia: SIAM) p10

  • [1] Feng Nai-Xing, Wang Huan, Zhu Zi-Xian, Dong Chun-Zhi, Li Hong-Yang, Zhang Yu-Xian, Yang Li-Xia, Huang Zhi-Xiang. Intelligent algorithm of extreme gradient boosting based perfectly matched monolayer and its efficient absorption on airborne transient electromagnetics problems. Acta Physica Sinica, 2024, 73(6): 065201. doi: 10.7498/aps.73.20231724
    [2] Gao Dong-Bao, Zhu Ji-Lin, Zhang Sai, Zhou He-Feng, Zeng Xin-Wu. Rayleigh-Bloch mode based monolayer bend waveguide. Acta Physica Sinica, 2021, 70(3): 034301. doi: 10.7498/aps.70.20201270
    [3] Chen Yun-Tian, Wang Jing-Wei, Chen Wei-Jin, Xu Jing. Reciprocal waveguide coupled mode theory. Acta Physica Sinica, 2020, 69(15): 154206. doi: 10.7498/aps.69.20200194
    [4] Tao Si-Cen, Chen Huan-Yang. Absorption characteristics of perfect absorber, electromagnetic “black hole” and inner perfectly matched layer. Acta Physica Sinica, 2020, 69(15): 154201. doi: 10.7498/aps.69.20200110
    [5] Li Xiao-Man, Zhang Ming-Hui, Zhang Hai-Gang, Piao Sheng-Chun, Liu Ya-Qin, Zhou Jian-Bo. A passive range method of broadband impulse source based on matched-mode processing. Acta Physica Sinica, 2017, 66(9): 094302. doi: 10.7498/aps.66.094302
    [6] Liu Jing, Xu Wei-Jiang, Hu Wen-Xiang. Evolution of the ultrasonic resonance modes in a three-layer structure with change of material and interface adhesion properties. Acta Physica Sinica, 2016, 65(7): 074301. doi: 10.7498/aps.65.074301
    [7] Chen Yu, Liu Long, Huang Zhong, Tu Lin-Lin, Zhan Peng. Great enhancement of transversal magneto-optical Kerr effect for magnetic dielectric film embedded by one-dimensional metallic grating. Acta Physica Sinica, 2016, 65(14): 147302. doi: 10.7498/aps.65.147302
    [8] Yang Qing, Cao Shu-Yang, Liu Shi-Yi. Numerical simulation of flow around two elongated rectangles in tandem arrangement using an immersed boundary method. Acta Physica Sinica, 2014, 63(21): 214702. doi: 10.7498/aps.63.214702
    [9] Zhang Tong-Wei, Yang Kun-De. A virtual time reversal method for passive source localization in a range-dependent waveguide. Acta Physica Sinica, 2014, 63(21): 214303. doi: 10.7498/aps.63.214303
    [10] Shen Wen-Yuan, Wang Hu, Geng Zhi-Hui, Du Chao-Hai, Liu Pu-Kun. Study of a W-band TE62 mode generator by the waveguide mode transformation. Acta Physica Sinica, 2013, 62(23): 238403. doi: 10.7498/aps.62.238403
    [11] Lai Ying-Xin, Yang Lei, Zhang Shi-Chang. Mode-matching analytic method of a coaxial Bragg structure corrugated with rectangular ripples and its experimental verification. Acta Physica Sinica, 2013, 62(20): 208402. doi: 10.7498/aps.62.208402
    [12] Zhang Rui, Yang Ya-Ping. Analysis of modes in  pairing of single-negative material. Acta Physica Sinica, 2010, 59(4): 2451-2456. doi: 10.7498/aps.59.2451
    [13] Chen Xian-Feng, Shen Xiao-Ming, Jiang Mei-Ping, Jin Yi. The mode characteristics of a symmetrical slab waveguide with single-negative-material cladding. Acta Physica Sinica, 2008, 57(6): 3578-3582. doi: 10.7498/aps.57.3578
    [14] Lin Min, Mao Qian-Min, Zheng Yong-Jun, Li Dong-Sheng. Frequency matching method for stochastic resonance control. Acta Physica Sinica, 2007, 56(9): 5021-5025. doi: 10.7498/aps.56.5021
    [15] Cai Xin-Lun, Huang De-Xiu, Zhang Xin-Liang. Application of full vectorial mode matching method to the calculation of eigenmodes of three-dimensional bent waveguide. Acta Physica Sinica, 2007, 56(4): 2268-2274. doi: 10.7498/aps.56.2268
    [16] Wang Jian, Lei Nai-Guang, Yu Chong-Xiu. Analysis of confinement loss in microstructured optical fibers with elliptical air holes. Acta Physica Sinica, 2007, 56(2): 946-951. doi: 10.7498/aps.56.946
    [17] Yu He-Jun, Xia Jin-Song, Yu Jin-Zhong. A highly efficient beam propagation method for modeling step-index waveguides with tilt interface. Acta Physica Sinica, 2006, 55(3): 1023-1028. doi: 10.7498/aps.55.1023
    [18] Xue Ting, Yu Jian, Yang Tian-Xin, Ni Wen-Jun, Li Shi-Chen. . Acta Physica Sinica, 2002, 51(3): 565-572. doi: 10.7498/aps.51.565
    [19] Guo Qi, Shi Zhi-Wei. . Acta Physica Sinica, 2002, 51(8): 1716-1723. doi: 10.7498/aps.51.1716
    [20] АНТИВОЛНОВОДНОЕ РАСПРОСТРАНЕНИЕ В НЕОДНОРОДНОМ СЛОЕ. Acta Physica Sinica, 1961, 17(4): 180-190. doi: 10.7498/aps.17.180
Metrics
  • Abstract views:  4518
  • PDF Downloads:  105
Publishing process
  • Received Date:  14 April 2023
  • Accepted Date:  11 May 2023
  • Available Online:  18 May 2023
  • Published Online:  20 July 2023

/

返回文章
返回