A Modified full-vectorial finite-difference beam propagation method based on H-fields for optical waveguides with step-index profiles

A modified full-vectorial finite-difference beam propagation method based on H-fields in solving the guided-modes for optical waveguides with step-index profiles is described. The propagation is split into two substeps. In the first substep, the field propagates in the absence of the cross-coupling terms, and then they are evaluated and double used in the second substep. An improved six-point finite-difference scheme is constructed to approximate the cross-coupling terms along the transverse directions. By using the imaginary-distance procedure, the field patterns and the normalized propagation constants of the fundamental modes for a buried rectangular waveguide and rib waveguide are presented, and the hybrid nature of the full-vectorial guided-modes is demonstrated. Solutions are in good agreement with the benchmark results from film mode matching method, which tests the validity and utility of the present method.     

Full-vectorial analysis of bending waveguides using finite difference method based on H-fields in cylindrical coordinate systems

A full-vectorial (FV) analysis of optical dielectric waveguide bends by using finite difference (FD) method in terms of magnetic field components is developed in a local cylindrical coordinate system. The perfectly matched layer absorbing boundary conditions via the complex coordinate stretching technique are incorporated into the FV wave equations for effectively demonstrating the leaky nature of waveguide bends, and a six-point FD scheme is constructed to approximate the cross-coupling terms for improving the convergent behavior. The leaky modes of a typical rib waveguide bend are calculated and the complex propagation constants and the field patterns for TE- and TM-like modes are obtained. Solutions are good agreement with those from the film mode matching method, which shows the validity and utility of the established method.     

Full-vectorial mode solver for anisotropic optical waveguides using multidomain spectral collocation method

A full-vectorial mode solver in terms of the transverse magnetic field components for optical waveguides with transverse anisotropy is described by using the multidomain spectral collocation method based on Chebyshev polynomials. The waveguide cross section surrounded by the perfectly matched layers is divided into suitable number of homogeneous rectangles, and then connected with by imposing the continuities of the longitudinal field components at the dielectric interfaces shared by the adjacent rectangles, resulting in a generalized matrix eigenvalue problem. To validate the established method, results of an anisotropic square waveguide and a magnetooptic rib waveguide are presented and compared with those from the full-vectorial finite difference method, full-vectorial beam propagation method, and the experimental data.     

Computation of higher-order waveguide modes by imaginary-distance beam propagation method

Abstact Yevick and Hermansson presented an efficient numerical method for calculating fundamental modes of optical waveguides. We extend their technique to higher-order eigenmodes. Using a finite-difference beam propagation method, we obtain propagation constants and field profiles for the three lowest-order TE modes in an asymmetric rib waveguide.     

Modified surface plasmonic waveguide formed by nanometric parallel lines

In this paper, two kinds of modified surface plasmonic waveguides formed by nanometric parallel lines are proposed. The finite-difference frequency-domain method is used to study propagation properties of the fundamental mode supported by these surface plasmonic waveguide structures. Results show that the transverse magnetic field of the fundamental mode is mainly distributed in the face to face region formed by two rods. With the same geometrical parameters and the same working wavelength of 632.8~nm, in the case of rods with a triangular cross-section, the degree of localization of field is strong, i.e. the mode area is small, but the fraction of the modal power in the metal increases, so the effective index increases and the propagation length of the mode decreases. With the same geometrical parameters, relative to the case of a working wavelength of 632.8~nm, when working wavelength is large, the mode area of transverse magnetic field distribution is large, i.e. the degree of localization of field is weak, and the interaction of field and silver is weak too, then the effective index decreases, so the propagation length increases. The rounded radii of rods have a great influence on the performance of the surface plasmonic waveguides with rounded triangular cross-sections, but have little influence on the performance of surface plasmonic waveguides with rounded square cross-sections. Since the distribution of transverse magnetic field, effective index, propagation length and the mode area can be adjusted by the geometrical parameters, this kind of modified surface plasmonic waveguide can be applied to the field of photonic device integration and sensors.     

Highly efficient noniterative full-vectorial beam propagation method based on the alternating-direction implicit scheme

A novel alternating-direction implicit method is developed for the full-vectorial wave equation. The cross-coupling terms are split off from the full-vectorial propagation operator, and their inversions are performed separately by a noniterative method at each step. The present method shows good performance in efficiency by comparison with the iterative method. It can be proved that the algorithm is of second-order accuracy along the propagating direction. To examine the accuracy of the present algorithm, it is compared with previously presented formulas, and the results show that it maintains high accuracy.     

一种水平变化波导中声传播问题的耦合模态法

针对介质参数及海底边界水平变化波导中的声传播问题,本文基于多模态导纳法提出一种能量守恒且便于数值稳定求解的耦合模态方法.将声压表示为一组正交完备的本地本征函数之和,对声压满足的Helmholtz方程在本地本征函数上作投影,推导出关于声压模态系数的二阶耦合模态方程组.耦合矩阵直观描述水平变化因素对模态耦合的贡献.为避免直接求解二阶耦合模态方程组可能遇到的数值发散问题,将其重构为两个耦合的一阶演化方程组,引入导纳矩阵并使用Magnus数值积分方法获得稳定的声场解.利用该耦合模态方法数值计算水平变化波导中的声场,并与COMSOL参考解比较,结果表明该耦合模态理论能够精确求解水平变化波导中的点源及分布源传播问题.     

Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method

Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer(CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain(ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor(PEC) waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.     

Solutions of the full- and quasi-vector wave equations based on H-field for optical waveguides by using mapped Galerkin method

The mapped Galerkin method in solving the full-vector and quasi-vector wave equations in terms of transverse magnetic fields (H-formulation) for optical waveguides with step-index profiles is described. By transforming the whole x-y space onto a unit square and using two-dimensional Fourier series expansion, the modal distributions and propagation constants for optical waveguides are obtained in the absence of boundary truncation. Results for step-index circular fiber, buried rectangular waveguide, and optical rib waveguide are presented. Solutions are good agreed with exact solutions and numerical results by using vector nonlinear iterative method, Fourier operator transform method, and vector beam propagation method.     

Semi-vectorial analysis of a compact wavelength demultiplexer based on the tapered multimode interference coupler

Based on a parabolically tapered multimode interference (MMI) coupler with a deep-etched SiO₂/SiON rib waveguide, a compact wavelength demultiplexer operating at 1.30 and 1.55\mum wavelengths is proposed and analysed by using three-dimensional semi-vectorial finite-difference beam propagation method (3D-SV-FD-BPM). The results show that a MMI section of 330.0\mum in length, which is only 76% length of a straight MMI coupler, is achieved with the contrasts of 42.3 and 39.2dB in quasi-TE mode, and 38.4 and 37.8dB in quasi-TM mode at wavelengths 1.30 and 1.55\mum, respectively, and the insertion losses below 0.2dB at both wavelengths and in both polarization states. The alternating direction implicit algorithm with the Crank--Nicholson scheme is applied to the discretization of the 3D-SV-FD-BPM formulation along the longitudinal direction. Moreover, a modified FD scheme is constructed to approximate the resulting equations along the transverse directions, in which the discontinuities of the derivatives of magnetic field components H_y and H_x along the vertical and horizontal interfaces, respectively, are involved.     

用基于磁场的改进型有限差分法分析光波导的全矢量本征模

提出一种改进的有限差分法,用以求解全矢量磁场波方程,分析光波导承载的全矢量本征模.离散交叉项时,采用与波导结构无关的六点差分格式,考虑磁场分量的导数在芯包分界处的不连续性,比传统四点差分格式有更高的计算精度.分析阶跃型光纤、矩形和脊形光波导的全矢量本征模问题,给出全矢量基模的磁场分布及其归一化传播常数,揭示全矢量本征模的混合特性,所得结果与解析法、高精度模横向谐振法的结果吻合,验证了方法的收敛性与计算精度.     

基于广角FD-BPM的PML边界处理方法

本文提出了基于三阶Padé近似广角有限差分光束传输法(FD-BPM)的简单有效的完美匹配层(PML)边界处理方法.分别给出了该方法应用于斜波导和多模干涉(MMI)数值模拟的结果,并与用完全透明边界条件处理的模拟计算结果进行了比较.本文最后对基于广角FD-BPM的PML边界处理方法进行了优化.     

Propagation of soliton waves in nonlinear dielectric slab waveguides of parabolic index of refraction

The pulse propagation in a non-linear slab waveguide of parabolic index of refraction is treated by using differential equation techniques. A graded index dielectric slab waveguide free of material dispersion with a cubic order non-linearity is considered. The electromagnetic wave inside the waveguide is described in terms of a non-linear equation. Slowly varying envelope function representation is employed to develop a non-linear partial differential equation for the unknown envelope function of the electric field. An averaging method over the transverse direction is applied to reduce the unknown envelope function non-linear differential equation into a form resembling the well known non-linear Schrödinger differential equation. This equation is solved by applying the Inverse Scattering Method. The N-soliton solution is developed and presented explicitly for the practical case of the single mode dielectric slab waveguide. Numerical results presenting single and double soliton propagation are also given.     

基于道格拉斯格式的宽角有限差分光束传播法分析大角度交叉波导

交叉波导是集成光学中的一种重要的器件单元，针对大角度交叉波导的结构及折射率分布，运用基于广义道格拉斯算子的宽角多步有限差分光束传播法模拟大角度交叉波导的传输特性，数值计算出波导中的光场的分布。与Crank-Nicholson格式的有限差分光束传播法相比，该算法在几乎不增加计算时间和计算机内存资源的情况下，结合了广义道格拉斯格式的有限差分光束传播法截断误差小[只有o(△x)^4]和宽角多步法所具有较高精度的优点，是光波导器件理想的数值模拟工具。在此基础上计算分析了大角度交叉波导的输出分支波导的能量与交叉角度的关系，以及交叉角度对波导间能量耦合的影响。这为实际的设计制作光器件的工作提供了极好的基础。     

Numerical calculations of modes of optical waveguides with two-dimensional refractive index profiles by a field correction method

Abstact We describe a numerical method devised to calculate the lower-order modes of optical waveguides with two-dimensional transverse index profiles. This method is based on finite-difference approximation of the scalar wave equation and the minimization of the error in the field shape of a mode. The field correction method can be applied to waveguides with large index steps and generates the mode's scalar propagation constant and its field. To test the field correction method, it is applied to a three-layer optical fibre and the results are compared with exact solutions. The error in the calculations of the modal effective indices is found to be of the order of 10^-6. Example mode calculations are given for fused couplers and rib waveguides.     

Modelling of self-aligned total internal reflection waveguide mirrors: an interlaboratory comparison

Results of modelling of light propagation in 45° self-aligned total internal reflection rib waveguide mirrors on InP substrate are compared. Six laboratories participated in the comparison with the following six modelling methods: the standard fast-Fouriertransform beam propagation method (BPM), the standard finite-difference (FD) BPM using the Crank-Nicholson scheme (two laboratories), the FD-BPM with the correction for the slowly varying envelope approximation, the method of lines, the eigenmode expansion and propagation method, and a simple method based on the field overlap. All the laboratories used the effective-index method to reduce the three-dimensional problem to two dimensions. The differences among the results obtained by different methods are briefly discussed and qualitatively compared to measured values.     

修正的频域有限差分法在二维金属光子晶体分析中的应用

与介质光子晶体相比,金属光子晶体的带隙特性在毫米波和亚毫米波波段有着重要的应用价值.基于Yee网格的频域有限差分法推导得出的本征模方程,求解后能方便而又可靠地得出介质光子晶体的带隙图和场分布.但由于金属与介质的本质差异,该方法不能直接应用于金属光子晶体.文中引入了金属表面边界条件,推导了二维金属周期结构的光子带隙本征模方程.通过数值计算,得出了不同晶格结构（正方/三角格子）下两种模式（TE/TM）的全禁带特性,并与介质周期结构的禁带特性进行对比,分析了金属周期结构在模式选择和器件集成方面的优点.     

修正的频域有限差分法在二维金属光子晶体分析中的应用

与介质光子晶体相比,金属光子晶体的带隙特性在毫米波和亚毫米波波段有着重要的应用价值.基于Yee网格的频域有限差分法推导得出的本征模方程,求解后能方便而又可靠地得出介质光子晶体的带隙图和场分布.但由于金属与介质的本质差异,该方法不能直接应用于金属光子晶体.文中引入了金属表面边界条件,推导了二维金属周期结构的光子带隙本征模方程.通过数值计算,得出了不同晶格结构（正方/三角格子）下两种模式（TE/TM）的全禁带特性,并与介质周期结构的禁带特性进行对比,分析了金属周期结构在模式选择和器件集成方面的优点. 关键词： 金属光子晶体 频域有限差分法 全禁带     

基于有限差分束传播法的Mach-Zehnder型波导电光调制器的模拟设计

针对掺钛铌酸锂Mach-Zehnder型波导电光调制器的结构及折射率分布,运用有限差分束传播方法(FD BPM)进行模拟计算,取得了二维情况下的最佳设计尺寸.在此基础上,把所得的设计结果与原有的用快速傅里叶束传播方法(FFT BPM)所得的结果进行比较和分析.结果显示,运用有限差分束传播方法不仅取得了预期的结果,而且,这种方法比原有的快速傅里叶束传播方法更快捷和方便,是器件设计的理想和现实的数值模拟工具.     

一种模拟倾斜折射率界面光波导的新方法

针对传统离散格式下，束传播方法(BPM)模拟倾斜折射率界面出现的问题，提出了一种简明且易于编程实现的改进方案. 在横向上，通过坐标系变换和插值处理，以新颖的7点差分格式代替传统的5点差分方式；在纵向上，以四阶显式Runge-Kutta方法(RKBPM)代替二阶Crank-Nicholson算法(CNBPM)，避免了求解不规则矩阵方程，从而使计算效率显著提高. 关键词： 束传播方法 Runge-Kutta方法 光波导 脊形波导