Stable and accurate method for modal analysis of multilayer waveguides using a graph approach

A method for calculating the propagation constants of allowed guided and leaky modes in multilayer planar waveguides is presented. We develop a two-way graph model to describe the tangential fields propagating in the waveguides. According to the special structure of the graph model, it is convenient to employ a topology scheme to derive analytical and closed-form dispersion equations for TE and TM modes. By comparing the dispersion equations formulated by series-expansion methods, approximation methods, and transfer-matrix methods, we find that the use of these equations for finding the eigenmodes has some benefits. First, this method can be easily employed to solve eigenmodes accurately in numerical computation without using series truncation. Second, the dispersion equations are exact. Moreover, all the eigenmodes can be determined according to the formulas without losing roots or causing numerical instability even for a waveguide with thick layers.     

Guiding mode expansion of a TE and TM transverse-mode integral equation for dielectric slab waveguides with an abrupt termination.

We propose a rigorous transverse-mode integral equation formulation for analyzing TE and TM electromagnetic radiation fields on the facet of dielectric slab waveguides with an abrupt termination in free space. Both exact waveguide guiding modes and discretized radiation modes are included in the kernels of the integral equation. To reduce the size of the matrix that approximates the exact integral equation, we expand the unknown field at the junctions in terms of guiding modes of a selected waveguide with sufficiently large normalized frequency and core thickness. By direct matrix inversion, we obtain numerical solutions of the scattered fields at the junctions. Our method can be used to study the field distribution as well as the energy reflection and transmission coefficients of dielectric waveguides with multiple step discontinuities.     

Reproduction in free space of fields confined by 3-dimensional optical waveguides

Abstract

It is shown by using the Fresnel-Kirchhoff diffraction theory and the method of images that a scalar field confined by a 3-dimensional optical waveguide can be generated in free space by a suitable light source. In the method the boundaries of a waveguide are replaced by virtual sources. This allows one to present the wave guiding as the interference and diffraction of multiple light beams produced in free space by the guide equivalent source (Fresnel waveguide). Thus, the scalar optics of a 3-dimensional waveguide is presented as the free-space beam optics. The method is illustrated by construction of the Fresnel sources for the triangular, rectangular and hexagonal waveguides. The numerical examples demonstrate the diffraction-free and self-imaging propagation in the free-space of the eigenmodes of the Fresnel rectangular-waveguide.     

Vectorial modal analysis of dielectric waveguides based on a coupled transverse-mode integral equation. II. Numerical analysis

We present numerical implementation and verification of a rigorous full-vector, integral-equation formulation suitable for analyzing modal characteristics of complex, two-dimensional (2D) rectangular-like dielectric waveguides. By dividing the waveguide into vertical slices, a system of integral equations we call vector-coupled transverse-mode integral equations (VCTMIE) is derived. The entire electromagnetic mode fields are completely determined by one-dimensional unknown field functions on the slice interfaces. To further reduce numerical computation, we expand these functions in terms of the guiding modes of a slab waveguide with a large normalized frequency. Through orthogonal projection the resulting nonlinear eigenvalue and eigenvector matrix formulation enables us to obtain the effective mode index with 10(-7) precision and to compute with high resolution the 2D vectorial mode field solutions of an open dielectric waveguide. We show stable and speedy convergence of our method as well as techniques to overcome the Gibbs phenomenon in the reconstruction of the transverse fields.     

Extended scattering-matrix method for efficient full parallel implementation of rigorous coupled-wave analysis

An extended and refined scattering-matrix method is proposed for the efficient full parallel implementation of rigorous coupled-wave analysis of multilayer structures. The total electromagnetic field distribution in the rigorous coupled-wave analysis is represented by the linear combination of the eigenmodes with their own coupling coefficients. In the proposed scheme, a refined recursion relation of the coupling coefficients of the eigenmodes is defined for complete parallel computation of the electromagnetic field distributions within multilayer structures.     

Possibility of excitation of azimuthal surface waves in a magnetized plasma by annular ion beams

Excitation of the azimuthal surface eigenmodes with the extraordinary polarization for the ion plasma component is shown to be possible in cylindrical waveguides with metal walls, filled partially by a cold magnetoactive neon plasma. The interaction of the oscillations with the flows of alpha particles rotating around the plasma column in a narrow layer that separates the plasma from the waveguide walls is studied. If the external magnetic field is strong enough, the resonance interaction of the beam with the waves can be realized for three minimum values of the azimuthal mode number when the wave frequency normalized by neon ion cyclotron frequency is close to the factor of five (m _Ne/m _He ≈ 5).     

基于双POD模型的空间相关三维随机风场数值模拟

POD法提供了一种高效、准确的风荷载模拟方法,通过对风场的功率谱密度矩阵进行Schur分解,得到一系列的特征值和特征向量,选取主要的几阶特征模态进行计算就可以得到比较精确的结果。该文讨论了具有空间相关性三维随机风场的数学模型,利用双POD模型和蒙特卡罗模拟法,详细描述了空间相关三维风场的数值模拟方法。通过大跨越输电塔三维风场的数值模拟研究表明,模拟的顺风向、横风向、竖直向风速的功率谱密度函数与理论值较为符合,并且具有较好的随机性。证实了该文提出的方法是一种高效、准确的结构三维风场模拟方法,并可应用于大跨空间结构、高层建筑及大跨度桥梁等结构之中。     

A procedure to obtain the complete electromagnetic scattering properties at discontinuities between integrated optical waveguides

Abstract

A method to obtain the complete electromagnetic scattering properties of discontinuities between arbitrary integrated optical waveguides is presented. The method involves a new generalized scattering matrix concept, together with the generalized telegraphist equations formulism and the modal matching technique. Radiation losses, as well as reflection and transmission coefficients between proper modes, can be obtained. Single and multiple discontinuities in planar and channel optical waveguides have been analysed. Numerical results of complex scattering coefficients are given. The possibilities of the method for analysing waveguide photonic crystals, as well as optical devices in waveguide periodic waveguide structures, are demonstrated.     

One-dimensional spatial dark soliton-induced channel waveguides in lithium niobate crystal

The anisotropic dependence of the formation of one-dimensional (1-D) spatial dark solitons on the orientation of intensity gradients in lithium niobate crystal is numerically specified. Based on this, we propose an approach to fabricate channel waveguides by employing 1-D spatial dark solitons. By exposure of two 1-D dark solitons with different orientations, channel waveguides can be created. The structures of the channel waveguides can be tuned by adjustment of the widths of the solitons and/or the angles between the two exposures. A square channel waveguide is experimentally demonstrated in an iron-doped lithium niobate crystal by exposure of two orthogonal 1-D dark solitons in sequence.     

Lossy multilayer channel optical waveguides analyzed by the transmission line matrix method

We demonstrate that the three-dimensional vectorial transmission line matrix (TLM) method is applicable to the analysis of lossy multilayer optical waveguiding structures. Any lossy multilayer waveguide geometry, including sharp discontinuities in the transverse plane, can be treated taking into account the coupling between all optical field components. The complex propagation constants (propagation constants and the attenuation coefficients) for the fundamental TE-like and TM-like modes can be determined. These parameters of the fundamental TM-like mode of a typical lossy multilayer rib dielectric waveguide are obtained as functions of free-space wavelength. Calculation of the electric-field pattern is also included. Numerical comparisons with the argument principle method (for the case of lossy multilayer slab waveguides) and the spectral-index technique (for the case of lossy multilayer rib waveguides) are also included, and it is shown that the application of the TLM method to lossy multilayer optical waveguides is accurate.     

Matrix method for two-dimensional waveguide mode solution

In this paper, we show that the transfer matrix theory of multilayer optics can be used to solve the modes of any two-dimensional (2D) waveguide for their effective indices and field distributions. A 2D waveguide, even composed of numerous layers, is essentially a multilayer stack and the transmission through the stack can be analysed using the transfer matrix theory. The result is a transfer matrix with four complex value elements, namely A, B, C and D. The effective index of a guided mode satisfies two conditions: (1) evanescent waves exist simultaneously in the first (cladding) layer and last (substrate) layer, and (2) the complex element D vanishes. For a given mode, the field distribution in the waveguide is the result of a ‘folded’ plane wave. In each layer, there is only propagation and absorption; at each boundary, only reflection and refraction occur, which can be calculated according to the Fresnel equations. As examples, we show that this method can be used to solve modes supported by the multilayer step-index dielectric waveguide, slot waveguide, gradient-index waveguide and various plasmonic waveguides. The results indicate the transfer matrix method is effective for 2D waveguide mode solution in general.     

Ion-exchanged glass waveguides with low birefringence for a broad range of waveguide widths

Optical communications networks require integrated photonic components with negligible polarization dependence, which typically means that the waveguides must feature very low birefringence. Recent studies have shown that waveguides with low birefringence can be obtained, e.g., by use of silica-on-silicon waveguides or buried ion-exchanged glass waveguides. However, many integrated photonic circuits consist of waveguides with varying widths. Therefore low birefringence is consequently required for waveguides having different widths. This is a difficult task for most waveguide fabrication technologies. We present experimental results on waveguide birefringence for buried silver-sodium ion-exchanged glass waveguides. We show that the waveguide birefringence of the order of 10(-6) for waveguide mask opening widths ranging from 2 to 10 microm can be obtained by postprocessing the sample through annealing at an elevated temperature. The measured values are in agreement with the values calculated with our modeling software for ion-exchanged glass waveguides. This unique feature of ion-exchanged waveguides may be of significant importance in a wide variety of integrated photonic circuits requiring polarization-independent operation.     

Reflectivity of buried slab waveguides

The scattering properties of an abruptly ended buried slab waveguide for both TE and TM modes are examined by an improved iteration technique that is based on the integral equation method with "accelerating" parameters. The waveguide is considered a symmetrical slab, for which the weakly guiding conditions are invalid, and it is embedded in a different dielectric material. The tangential electric field distribution on the terminal plane, the reflection coefficient of the first TE and TM guided modes, and the far-field radiation pattern are computed. Numerical results are presented for several ended waveguides, while special attention is given to the far-field radiation pattern rotation and the terminal field distributions.     

Bend loss effects in diffused, buried waveguides

Bend loss effects can be a significant concern in the design and performance of diffused, buried waveguide devices. Since diffused, buried waveguides typically do not have analytical mode solutions, the bend mode must be expressed as an expansion of straight waveguide modes. For the case of buried ion-exchanged waveguides, the bend loss is affected by bend radius, the duration of the ion exchange and burial processes, as well as the size of the mask opening used to create the waveguides and applied field during burial. The bend loss effects for each of these variables are explored under typical fabrication conditions.     

Scattering of SH-waves by an interface cavity

Summary. The scattering of the SH-wave and dynamic stress concentrations near an arbitrary cavity situated at the planar interface separating two different elastic media are investigated. The total wave field can be obtained by superposition of the free field and the scattered field. The free field is composed of the incident, reflected and refracted waves. The scattered wave fields in adjacent media are expressed respectively, and the method of wave functions expansion is applied to obtain the solutions for these fields. The scattered wave functions can be expanded into Hankel-Fourier series with unknown coefficients. In solving for the unknown coefficients according to the boundary conditions for the total wave field at the interface and at the cavity wall, the non-orthogonality makes the system of equations for the unknown coefficients infinite and coupling each other. Another key point is to extend each scattered wave field from its own half-plane domain into the full plane domain by a certain way keeping the total wave field unchanged for the non-orthogonal Fourier integrals around the cavity. Finally, the scattering of the SH wave by an interface ellipse with different ratios between long and short axis is considered, and the distributions of dynamic stress concentration factors at the cavity wall are presented.     

Analysis of photonic crystal waveguide gratings with coupled-mode theory and a finite-element method

Both rod and air-hole types of photonic crystal waveguide gratings are proposed and their coupling coefficients and transmission characteristics are effectively investigated by using a simple coupled-mode theory combined with a finite-element method. The results obtained are compared with the results obtained by using a full numerical simulation method. A new definition for unperturbed waveguides is introduced to obtain accurate coupling coefficients. It is shown that, by using a pi-phase-shifted waveguide structure in the case of an air-hole type of photonic crystal waveguide grating, the coupling coefficient is strongly enhanced. The accuracy of the method is discussed through numerical examples of high-index-contrast waveguide gratings.     

Direct integration of nanoscale conventional and slot waveguides

Slot waveguides can provide high optical confinement in a nanoscale low-index layer. While a conventional waveguide has a Gaussian-like Eigenmode profile, the Eigenmode profile of a slot waveguide is quite non-Gaussian type, due to the large discontinuity of refractive indices and thus the transverse electric field component between the high and low index layers of a slot waveguide. Although the field profiles of the two types of waveguides seem different, here we show that direct integration of conventional and slot waveguides yields efficient coupling of light into and out of slot waveguides using the rigorous finite-difference time domain method. The proposed direct coupling method has comparable performance to recently proposed taper based coupling methods, while having advantages in easier integration with conventional waveguide optics and higher integration density. We also show that coupling efficiency is not sensitive to the symmetricity of the slot waveguide, resulting in good manufacturing tolerance. The proposed direct coupler may have a number of applications in lightwave interconnects, sensing and data storage.     

Design of an ultrashort Si-nanowaveguide-based multimode interference coupler of arbitrary shape

A design procedure for an arbitrarily tapered multimode interference (MMI) coupler based on Si nanowaveguides is presented. First a series of the effective indices of the zeroth and first eigenmodes in multimode waveguides are obtained as the core width increases by using a full-vectorial finite-difference method. Two polynomial functions are used to fit the two relations between the effective indices and the core width. The phase difference Delta phi between the zeroth and first eigenmodes can then be easily calculated when the light goes through any given arbitrarily tapered MMI section. By making the phase difference Delta phi equal to a certain value Delta phi 0 required for an N-fold self-imaging, the length of a MMI coupler is determined. With the present design procedure, an ultrashort 2 x 2 parabolic MMI coupler is designed as an example. The size of the designed ultrashort MMI section is only approximately 1.4 microm x 4.7 microm.     

主支式排烟道中主烟道横截面尺寸的确定

高层住宅主支式排烟道的主烟道中烟气流动可看作是一维定常垂直摩擦管流。通过计算管流流动相与颗粒相的物理参数沿管程分布情况,确定了颗粒相流动情况,进而判断烟气排放情况。这受到排烟系统高度、烟道管径,及材料组成等因素的影响,通过预设排烟系统高度并保证良好排放情况,求得主烟道管径应满足的条件。这使得布设排烟道时可以找到一个较好的横截面尺寸,既满足烟气排放条件,又节约修建材料。     

Transmission performance of one waveguide and double micro-ring resonator using 3×3 optical fiber coupler

Abstract

This paper investigates theoretically the transmission characteristics of one waveguide and double micro-ring resonator using 3 × 3 optical fibre coupler. Our analytical solution of transmittance is suitable for either linearly distributed coupler or circularly symmetric distributed coupler. The all-optical analogue to electromagnetic inducted transparency spectrum of one waveguide and double micro-ring resonators can be controlled by changing the coupling strength between waveguide and micro-rings, the absorption coefficient around micro-rings, as well as the asymmetric coupling coefficients between non-adjacent waveguides. The curves show that the transitions of transmission spectra sensitively depend on asymmetric coupling coefficients.