Full-vectorial finite-difference beam propagation method based on the modified alternating direction implicit method

A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on $H$ fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance 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.     

Full-vectorial analysis of optical waveguides by the finite difference method based on polynomial interpolation

Based on the polynomial interpolation, a new finite difference (FD) method in solving the full-vectorial guided-modes for step-index optical waveguides is proposed. The discontinuities of the normal components of the electric field across abrupt dielectric interfaces are considered in the absence of the limitations of scalar and semivectorial approximation, and the present FD scheme can be applied to both uniform and non-uniform mesh grids. The modal propagation constants and field distributions for buried rectangular waveguides and optical rib waveguides are presented. The hybrid nature of the vectorial modes is demonstrated and the singular behaviours of the minor field components in the corners are observed. Moreover, solutions are in good agreement with those published early, which tests the validity of the present approach.     

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.     

全矢量有限差分法分析任意截面波导模式

针对不规则端口模式加载的需求, 改进了全矢量有限差分法, 推导了边界条件, 建立了模式求解器, 解决了不规则端口任意模式加载的难题. 该求解器对计算资源要求较低, 可以求解任意形状的波导模式. 计算了不同形状波导的模式特性, 得到了与解析解和商用软件结果相一致的计算结果. 关键词： 全矢量有限差分法 波导 模式     

Beam propagation analysis using higher-order full-vectorial finite-difference method

We develop a beam propagation method based on full-vectorial higher-order finite-difference formulation to model optical waveguides. Simulation of a photonic crystal fiber coupler is used as assessment of different schemes and convergence behaviors. Results including error convergence and field calculation show that they are improved by using higher-order scheme.     

基于插值的全矢量有限差分法求解阶跃光波导本征模式

给出了一种适合求解阶跃折射率分布光波导本征模式的全矢量有限差分法.该方法充分考虑了电场分量在折射率分布跃变处的不连续性, 在拉格朗日插值的基础上得到了一种精度较高的差分格式.为了验证其有效性,对一种典型结构、折射率阶跃分布的波导作了计算,并与其它方法的结果作了比较.最后对不理想刻蚀情况做了模拟计算.     

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.     

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

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

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.     

Enhanced transmission through periodic arrays of subwavelength holes: the role of localized waveguide resonances

By using the rigid full-vectorial three-dimensional finite-difference time-domain method, we show that the enhanced transmission through a metallic film with a periodic array of subwavelength holes results from two different resonances: (i) localized waveguide resonances where each air hole can be considered as a section of metallic waveguide with both ends open to free space, forming a low-quality-factor resonator, and (ii) well-recognized surface plasmon resonances due to the periodicity. These two different resonances can be characterized from electromagnetic band structures in the structured metal film. In addition, we show that the shape effect in the enhanced transmission through the Au film with subwavelength holes is attributed to the localized waveguide resonance.     

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

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

1D Bessel-like beam generation by highly directive transmission through a sub-wavelength slit embedded in periodic metallic grooves

A new method to generate a quasi-one dimensional (1D) Bessel-like beam whose non-diffracting length extends to macroscopic scale is proposed by utilizing highly directive transmission through a sub-wavelength single slit embedded in periodic metallic grooves. We employed finite-difference time-domain methods for full-vectorial diffraction and spectral analysis. In optimal conditions, unique quasi-1D Bessel-like beams were generated in transmission for the incident p-polarized plane wave such that highly asymmetric far-field distribution was achieved; Bessel-like beam along the slit axis and a flat-top super-Gaussian beam along the other perpendicular axis. Detailed parametric studies for the proposed structure are reported in terms of the operating spectral range and a general window of optimal conditions.     

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

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

Waveguide mode solver based on Neumann-to-Dirichlet operators and boundary integral equations

For optical waveguides with high index-contrast and sharp corners, existing full-vectorial mode solvers including those based on boundary integral equations typically have only second or third order of accuracy. In this paper, a new full-vectorial waveguide mode solver is developed based on a new formulation of boundary integral equations and the so-called Neumann-to-Dirichlet operators for sub-domains of constant refractive index. The method uses the normal derivatives of the two transverse magnetic field components as the basic unknown functions, and it offers higher order of accuracy where the order depends on a parameter used in a graded mesh for handling the corners. The method relies on a standard Nyström method for discretizing integral operators and it does not require analytic properties of the electromagnetic field (which are singular) at the corners.     

Enhancing surface plasmon polariton propagation by two-layer dielectric-loaded waveguides on silver surface

We report a study of a two-layer dielectric-loaded surface plasmon polariton waveguide (TDLSPPW) which consisted of two dielectric layers (high-index/low-index) on a silver film. The discontinuity of the electric field at the interfaces resulted in a concentrated field in the low-index region. It efficiently reduced the propagation loss of the surface plasmon polariton mode. The mode fields and corresponding complex propagation constants were calculated by a vector finite-difference method. The propagation properties were measured by a modified near-field optical microscope. It is confirmed that the propagation length of the proposed TDLSPPW was about 1.6 times longer than conventional single-layer SPP waveguides. In addition, a 90^° waveguide turn with 3 μm radius showed that the bending loss was smaller than 2 dB.     

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.     

Optical and microwave analysis of mushroom-type waveguides for traveling wave electroabsorption modulators based on asymmetric intra-step-barrier coupled double strained quantum wells by full-vectorial method

The finite difference method is exploited for a full-vectorial analysis of mushroom-type waveguides for traveling wave electroabsorption modulators (TWEAM) based on asymmetric intra-step-barrier coupled double strained quantum wells (AICD-SQW). In this analysis, the discontinuities of the normal components of the electric field across abrupt dielectric interfaces which are known as the limitations of scalar and semivectorial approximation methods are considered. The optical field distributions in mushroom-type TWEAM based on AICD-SQW and conventional ridge-type TWEAM of the same active region for 1.55 μm operation are presented. The important parameters in the high-frequency TWEAM design such as optical effective index which defines optical velocity and transverse mode confinement factor are calculated. Then, the transmission line microwave properties (microwave index, microwave loss, and characteristic impedance) of TWEAMs are obtained. The modulation response of mushroom-type TWEAM is calculated using circuit model by considering interaction between microwave and optical fields in waveguide and compared with conventional ridge-type TWEAM. It is found that increasing the width of p-cladding layer with the same i-layer to reduce the resistance in p-i-n mushroom-type waveguide of TWEAM based on AICD-SQW can improve the microwave propagation loss and thus the high-speed electro-optical response.