Construction of glass waveguide refractive index profiles by the effective-index finite-difference method

A numerical method is applied to construct the refractive index profiles of optical waveguides from the measured effective indices (EI). The method is based on choosing a proper analytical function for the refractive index profile and searching its unknown parameters using the simplex search algorithm. Simultaneously, the finite-difference method (FDM) is used to solve the semi-vectorial Helmholtz equation for the guided modes effective indices. The method is applied successfully to two particular Ag⁺---Na⁺ ion-exchanged glass slab waveguides. The results are as accurate as those obtained from from commonly used IWKB-based method. The EI-FDM in principle can be applied to both slab and channel waveguides and does not require that the index profiles are monotonically decreasing, like most of IWKB-based methods. The relation between the induced refractive index and silver concentration profile, measured by SIMS, is found to be almost linear.     

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.     

基于非协调边界元法的声辐射模态分析

利用非协调单元离散声学Helmholtz边界积分方程,采用极坐标变换法消除积分奇异性,通过CHIEF方法加Lagrange乘子法处理特征频率处解的不唯一性。在此基础上,应用非协调单元推导结构的声辐射功率和声辐射效率的表达式。以脉动球和辐射立方体为例,计算结构的声辐射功率、辐射效率、辐射模态、辐射模态效率等物理量,并与协调单元的计算结果做比较,取得较好的一致性。     

Attenuation in thin film optical waveguides due to roughness-induced mode coupling

Marcuse's mode-coupling theory for symmetrical planar optical waveguides was extended to cover the asymmetrical case. Computer programs were written to study the particular case of coupling between the zeroth order guided transverse electric (TE) mode and the radiation modes as a function of interface roughness.It is shown that a simple expression proposed by Tien gives results which correspond to the most pessimistic results obtained by the mode-coupling theory. R.f. sputtered thin film planar optical waveguides were prepared. The zeroth order TE guided mode having a vacuum wavelength of 0.6328 μm was excited in the thin film by means of a prism coupler. The attenuation of this mode was inferred from the reduction in intensity of the radiation modes which was measured with a photomultiplier. A Talystep was used to obtain the film thickness and also profiles of the roughness of the film-substrate and film-superstrate interfaces. The roughness profiles were computer-processed to determine the r.m.s. amplitude and the correlation length of the interface roughness. The refractive index was determined from the film thickness and the synchronous coupling angle using the reduced SYNCANG method. These data were used to compute the attenuation of the zeroth order guided TE mode from the mode-coupling theory. The correlation coefficient between the theoretical and experimental results was 0.96.     

Resonant scattering and mode coupling in two-dimensional textured planar waveguides

A heuristic formalism is developed for efficiently determining the specular reflectivity spectrum of two-dimensionally textured planar waveguides. The formalism is based on a Green's function approach wherein the electric fields are assumed to vary little over the thickness of the textured part of the waveguide. Its accuracy, when the thickness of the textured region is much smaller than the wavelength of relevant radiation, is verified by comparison with a much less efficient, exact finite difference solution of Maxwell's equations. In addition to its numerical efficiency, the formalism provides an intuitive explanation of Fano-like features evident in the specular reflectivity spectrum when the incident radiation is phase matched to excite leaky electromagnetic modes attached to the waveguide. By associating various Fourier components of the scattered field with bare slab modes, the dispersion, unique polarization properties, and lifetimes of these Fano-like features are explained in terms of photonic eigenmodes that reveal the renormalization of the slab modes due to interaction with the two-dimensional grating. An application of the formalism, in the analysis of polarization-insensitive notch filters, is also discussed.     

On the spectrum of the electric field integral equation and the convergence of the moment method

Existing convergence estimates for numerical scattering methods based on boundary integral equations are asymptotic in the limit of vanishing discretization length, and break down as the electrical size of the problem grows. In order to analyse the efficiency and accuracy of numerical methods for the large scattering problems of interest in computational electromagnetics, we study the spectrum of the electric field integral equation (EFIE) for an infinite, conducting strip for both the TM (weakly singular kernel) and TE polarizations (hypersingular kernel). Due to the self‐coupling of surface wave modes, the condition number of the discretized integral equation increases as the square root of the electrical size of the strip for both polarizations. From the spectrum of the EFIE, the solution error introduced by discretization of the integral equation can also be estimated. Away from the edge singularities of the solution, the error is second order in the discretization length for low‐order bases with exact integration of matrix elements, and is first order if an approximate quadrature rule is employed. Comparison with numerical results demonstrates the validity of these condition number and solution error estimates. The spectral theory offers insights into the behaviour of numerical methods commonly observed in computational electromagnetics. Copyright © 2001 John Wiley & Sons, Ltd.     

Angular output of hollow, metal-lined, waveguide Raman sensors

Hollow, metal-lined waveguides used as gas sensors based on spontaneous Raman scattering are capable of large angular collection. The collection of light from a large solid angle implies the collection of a large number of waveguide modes. An accurate estimation of the propagation losses for these modes is required to predict the total collected Raman power. We report a theory/experimental comparison of the Raman power collected as a function of the solid angle and waveguide length. New theoretical observations are compared with previous theory appropriate only for low-order modes. A cutback experiment is demonstrated to verify the validity of either theory. The angular distribution of Raman light is measured using aluminum and silver-lined waveguides of varying lengths.     

Characterization of supported cylinder-planar germanium waveguide sensors with synchrotron infrared radiation

Cylinder-planar Ge waveguides are being developed as evanescent-wave sensors for chemical microanalysis. The only non-planar surface is a cylinder section having a 300-mm radius of curvature. This confers a symmetric taper, allowing for direct coupling into and out of the waveguide's 1-mm(2) end faces while obtaining multiple reflections at the central <30-microm-thick sensing region. Ray-optic calculations indicate that the propagation angle at the central minimum has a strong nonlinear dependence on both angle and vertical position of the input ray. This results in rather inefficient coupling of input light into the off-axis modes that are most useful for evanescent-wave absorption spectroscopy. Mode-specific performance of the cylinder-planar waveguides has also been investigated experimentally. As compared to a blackbody source, the much greater brightness of synchrotron-generated infrared (IR) radiation allows a similar total energy throughput, but restricted to a smaller fraction of the allowed waveguide modes. However, such angle-selective excitation results in a strong oscillatory interference pattern in the transmission spectra. These spectral oscillations are the principal technical limitation on using synchrotron radiation to measure evanescent-wave absorption spectra with the thin waveguides.     

Modal behavior in single-mode slab waveguides consisting of inhomogeneous cores

We present an analytical method for TE and TM modes in weakly guiding inhomogeneous single-mode slab waveguides. Based on our results, the modal behavior or propagation constants depend on index profiles of the waveguides. It is important to know how the modal behavior depends on the index profile in single-mode waveguides, because it determines wave-front characteristics of propagating modes.     

Swift and heavy ion implanted chalcogenide laser glass waveguides and their different refractive index distributions

Planar waveguides have been fabricated in Nd- or Ho-doped gallium lanthanum sulfide laser glasses by 60?MeV Ar or 20?MeV N ion implantation. The refractive index profiles were reconstructed based on the results of prism coupling. The Ar implanted waveguides exhibit an approximate steplike distribution, while the N implanted ones show a "well + barrier" type. This difference can be attributed to the much lower dose of Ar ions. After annealing, the N implanted waveguides can support two modes at 1539?nm and have low propagation loss, which makes them candidates for novel waveguide lasers.     

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.     

Characterization of a compact silicon-based slot-to-strip waveguide crossing

A compact crossing scheme for a silicon-based slot-to-strip (or vice versa) waveguide is proposed and analyzed by using a finite-difference time-domain method, where a strip–multimode waveguide (SMW) crossing is used at the center and two logarithmically tapered slot-to-strip mode converters are incorporated into the ports with slot waveguides. For the input ports with slot waveguides and output ports with strip waveguides, the guided modes are efficiently transformed through the mode converter, and then enter into the SMW, where the fields converge at the center of the intersection due to the self-imaging effect. Hence, the size of the input beam is much smaller than the width of the SMW at the crossing center, leading to significant reductions of both crosstalk and radiation loss. The numerical results show that a slot-to-strip waveguide crossing operating at a wavelength of 1.55?μm can be achieved with insertion loss, crosstalk, and reflection of 0.134/0.182, ?36.18/?38.6, and ?35.8/?42.02?dB for input ports with slot/strip waveguides, respectively.     

Annealing behavior of barriers in ion-implanted LiNbO3 and LiTaO3 planar waveguides

We formed planar waveguides in LiNbO(3) and LiTaO(3) crystals by megaelectron volt He-ion implantation. The dark modes of both waveguides are measured and their refractive-index profiles are described according to the parameterized index profile reconstruction method. The extraordinary indices of both ion-implanted waveguides exhibit quite different profiles. We compare the thermal stability of barriers in ion-implanted LiNbO(3) and LiTaO(3) waveguides by annealing at different temperatures. The results show that the barrier in a LiTaO(3) planar waveguide has higher thermal stability than that in a LiNbO(3) waveguide. The experiments also show that annealing at a temperature higher than 400 degrees C results in recrystallization of the barrier.     

基于有限点振速分布设计误差传感策略

误差传感策略是实现主动结构声学控制的关键的一环。在中低频率时，控制振动结构前几阶声辐射模态可以有效控制总声功率，基于声辐射模态进行板结构主动声学控制的关键是如何获得前几阶声辐射模态伴随系数。在振动平板上测量少数点振动速度分布，利用声辐射模态性质，通过求解欠定方程，得到所需要的前N阶声辐射模态伴随系数。利用得到的前N阶声辐射模态伴随系数作为控制器的输入，形成基于声辐射模态的主动控制策略和相应的误差传感策略。以固定支撑板为例，从主动控制效果分析得到的结果与理论值一致，说明利用上述误差传感策略得到前几阶声辐射模态伴随系数是可行的。     

Full vectorial simulation of multilayer anisotropic waveguides with an accurate and automated finite-element program

An efficient, accurate, and automated vectorial finite-element software package (named WAVEGIDE), which is implemented within a PDE/Protran problem-solving environment, has been extended to general multilayer anisotropic waveguides. With our system, through an interactive question-and-answer session, the problem can be simply defined with high-level PDE/Protran commands. The problem can then be solved easily and quickly by the main processor within this intelligent environment. In particular, in our system the eigenvalue of waveguide problems may be either a propagation constant (β) or an operated light frequency (F). Furthermore, the cutoff frequencies of propagation modes in waveguides can be calculated. As an application of this approach, numerical results for both scalar and hybrid modes in multilayer anisotropic waveguides are presented and are also compared with results obtained with the domain-integral method. These results clearly illustrate the unique flexibility, accuracy, and the ease of use f the WAVEGIDE program.     

Greatly reduced radiation loss in planar waveguides with two-dimensional conducting interfaces

A new strategy for radiation loss reduction in curved slab waveguides is presented. The proposed strategy is based on the proper modification of the boundary conditions at the core-to-cladding interface, whereupon extremely thin conductive nanolayers with non-zero surface conductance are imposed. The obtained numerical results show a noticeable decrease in the overall loss level.     

Theoretical analysis of phase-matched second-harmonic generation and optical parametric oscillation in birefringent semiconductor waveguides

We analyze the phase-matching conditions for second-harmonic generation (SHG) and optical parametric oscillation (OPO) in birefringent nonlinear semiconductor waveguides and apply these results to the model system of ZnGeP(2) on a GaP substrate. The analyses and numerical results show that phase matching can be achieved for OPO and SHG for reasonable guide thicknesses throughout much of the infrared, indicating significant potential applications for nonlinear birefringent waveguides. For the fundamental mode of a relatively thick guide the region of phase matching and the phase-matching angles are similar to those in bulk material. However, the waveguide has the added flexibility that phase-matched coupling can occur between the various modes of the guide. For example, the phase-matching region for SHG can be considerably extended by coupling the pump into the guide in the fundamental, m =, mode and phase matching to the m = mode of the second harmonic. Significantly, the results indicate, among other things, that ZnGeP(2) waveguides with harmonic output in the m = mode can be used for efficient SHG from input radiation in the 9.6-10.6-mum region where bulk efficiencies in this wavelength range are too small to be useful.     

Synthesis of symmetric and asymmetric planar optical waveguides

A novel and accurate refractive index profile synthesis method for planar optical waveguides is presented and demonstrated using the transmitted near-electric-field-data. This method is based on the inverse transmission-line (TL) technique. From Maxwell's equations, a TL equivalent circuit (electric T-circuit) for the refractive index profile of a planar optical waveguide is derived. The authors demonstrate how to use this model to carry out the inverse problem and synthesise the exact refractive index profile numerically from near-field-data. The TL method can reconstruct arbitrary refractive index profiles for planar optical waveguides that support singlemode or multi- modes. The cases of both symmetric and asymmetric arbitrary refractive index profile planar waveguides are discussed. The accuracy of the reconstructed waveguides is examined numerically.