Analysis of surface relief diffraction gratings made of anisotropic material

A modal method for the analysis of surface relief gratings made with anisotropic material is presented. The structure is decomposed into a series of cascaded discontinuities between planar waveguides with stratified anisotropic dielectric. The basic problem is formulated by an integral equation which is solved numerically by the method of moments. The mode functions of the periodic region are assumed as basis functions to represent the unknown field on the junctions. Each junction is viewed as a waveguide junction problem and has been characterized by the generalized scattering matrix (GSM). The diffraction efficiencies of the grating are determined by combining the various GSM. In this way, the analysis method is stable and can be applied also to deep gratings.     

波长尺寸表面浮雕光栅的严格理论计算和分析

本文采用麦克斯韦(Maxwell)方程微分积分方程对波长尺寸的梯形反射光栅在TE模式(横向电场)偏振状态入射下,进行严格理论的计算和分析.表面电场计算表明,光栅刻槽斜边边界对于金属材料金和铝光栅表面电场分布具有明显效应,对于入射波长为λ=1.0μm附近的红外区域,金材料具有强烈的金属特性,此时表面电场的斜边边界效应更加明显.表面电场及其远场衍射效率伴随刻槽深度变化,而振荡特性进一步阐明了此效应对近,远场的作用.该方程适用于任意形状光栅或单一散射体对任意入射方向,任意偏振态的衍射和散射问题.     

Canonical Quantization of Electromagnetic Field in the Presence of Nonlinear Anisotropic Magnetodielectric Medium with Spatial-Temporal Dispersion

Modeling a nonlinear anisotropic magnetodielectric medium with spatial-temporal dispersion by two continuum collections of three dimensional harmonic oscillators, a fully canonical quantization of the electromagnetic field is demonstrated in the presence of such a medium. Some coupling tensors of various ranks are introduced that couple the magnetodielectric medium with the electromagnetic field. The polarization and magnetization fields of the medium are defined in terms of the coupling tensors and the oscillators modeling the medium. The electric and magnetic susceptibility tensors of the medium are obtained in terms of the coupling tensors. It is shown that the electric field satisfy an integral equation in frequency domain. The integral equation is solved by an iteration method and the electric field is found up to an arbitrary accuracy.     

Application of the beam propagation method to diffused transmission gratings

In this paper results for the diffraction characteristics of phase gratings calculated by the beam propagation method and by a rigorous integral equation method are presented, with the latter results serving as a reference. The configurations under consideration have refractive index changes which are realistic for ion exchange processes in glass. Thus, multiple reflections have only little influence which makes the beam propagation method applicable to the problem. Typically, the diffraction efficiencies computed by both methods differ by less than 2%. In contrast to the integral equation method, the beam propagation method works with a definable lateral field profile, enabling the modelling of a grating within a complete device design including the near field pattern. This revised version was published online in November 2006 with corrections to the Cover Date.     

各向异性介质三维电磁响应模拟的Ho-GEBA算法

本文基于积分方程法研究并建立了一种模拟横向同性介质中任意各向异性异常 体三维电磁响应的高阶广义扩展Born近似(Ho-GEBA)算法. 首先利用逐次迭代技术给出积分方程的广义级数展开解, 为保证其收敛性, 引入一种各向异性条件下满足压缩映射的迭代算子. 然后利用异常体区域分解技术, 并结合扩展Born近似原理, 得到各向异性介质三维电磁响应的Ho-GEBA解. 为提高效率, 计算过程中采用并矢Green函数的解析表达式. 最后通过数值计算实例对比验证了本文算法的有效性. 关键词： 高阶广义扩展Born近似 积分方程 电磁模拟 解析Green函数     

Anisotropic Bragg diffraction of finite-sized volume holographic grating in photorefractive crystals

Anisotropic diffraction of uniform plane wave by finite-sized volume holographic grating in photorefractive crystals is considered. It is found that the anisotropic diffraction can take place when some special conditions are satisfied. The diffracted image is obtained in experiment for the anisotropic Bragg diffraction in Fe:LiNbO₃ crystals. A coupled wave analysis is presented to study the properties of anisotropic diffraction. An analytical integral solution for the amplitudes of the diffracted beams is submitted. A trade off between high diffraction efficiency and the deterioration of reconstruction fidelity is analyzed. Numerical evaluations also show that the finite-sized anisotropic volume grating exhibits strong angular and wavelength selectivity. All the results are useful for the optimizing design of VHOE based on finite-sized volume grating structures.     

Genetic algorithm-based design method for multilevel anisotropic diffraction gratings

We developed a method for the design of multilevel anisotropic diffraction gratings based on a genetic algorithm. The method is used to design the multilevel anisotropic diffraction gratings based on input data that represent the output from the required grating. The validity of the proposed method was evaluated by designing a multilevel anisotropic diffraction grating using the outputs from an orthogonal circular polarization grating. The design results corresponded to the orthogonal circular polarization grating structures that were used to provide outputs to act as the input data for the process. Comparison with existing design methods shows that the proposed method can reduce the number of human processes that are required to design multilevel anisotropic diffraction gratings. Additionally, the method will be able to design complex structures without any requirement for subsequent examination by a human designer. The method can contribute to the development of optical elements by designing multilevel anisotropic diffraction gratings.     

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.     

Collinear diffraction of divergent optical beams

Two variants of anisotropic collinear diffraction of light on a sinusoidal diffraction grating are studied theoretically with forth and back light scattering relative to the direction of the incident beam. For both variants, two-dimensional transfer functions are calculated and their modification with optical wavelength and diffraction grating period is analyzed. The dependence of integral diffraction efficiency and filter passband on optical beam divergence is examined.     

Analysis of stress field in a single anisotropic diamond/square shaped inclusion using the volume integral equation method and the numerical equivalent inclusion method

A volume integral equation method (VIEM) is used to study elastostatic problems in an unbounded elastic solid containing a single diamond/square shaped inclusion subject to uniform tensile stress at infinity. The inclusion is assumed to be a long parallel diamond/square cylinder composed of isotropic or anisotropic elastic materials and perfectly bonded to the isotropic matrix. The solid is assumed to be under plane strain on the plane normal to the cylinder. A detailed analysis of the stress field at the interface between the isotropic matrix and the single isotropic/orthotropic diamond/square shaped inclusion is carried out. The effects of a single isotropic/orthotropic diamond/square shaped inclusion on the stress field at the interface between the matrix and the inclusion are investigated in detail. The accuracy of the volume integral equation method for the interfacial stress field is validated and compared by the numerical equivalent inclusion method (NEIM) and the finite element method (FEM) using ADINA. Through detailed analysis of plane elastostatic problems using the parallel volume integral equation method (PVIEM) in an unbounded isotropic matrix with multiple isotropic diamond shaped inclusions under uniform remote tensile loading, it is demonstrated that the volume integral equation method can also be applied to solve general two- and three-dimensional elastostatic problems involving multiple isotropic/anisotropic inclusions whose shape and number are arbitrary.     

Diffraction of cylindrical waves by an ideally conducting wedge in an anisotropic plasma

The problem of diffraction of cylindrical waves by an ideally conducting wedge in an anisotropic plasma is formulated and solved. The integral equations for the field are reduced to function equations, which are solved with the aid of a special function that is introduced. The properties of this function are studied. The general solution is represented as a double contour integral in the plane of a complex variable. The radiation field and surface waves for a number of special cases are analyzed: a source of cylindrical waves on an edge; at infinity; etc. Diffraction in a half-plane is studied separately.     

Determination of unsteady heat release distribution from acoustic pressure measurements: a reformulation of the inverse problem

An integral method is developed to solve the inverse problem of determining the oscillatory heat release distribution from the knowledge of the acoustic pressure field within a combustor. Unlike earlier approaches, in which the problem is formulated in terms of Fredholm integral equation, the inverse problem is reformulated in terms of Volterra integral equation. This reformulation, valid for low Mach numbers (M2 < 1), facilitates the recovery of heat release at all frequencies. The resulting Volterra integral equation is solved using both direct numerical method and implicit least-squares method. The results show that the implicit least-squares method is superior to the direct numerical method and yields accurate determination of heat release at all frequencies.     

Stress analysis of multiple anisotropic elliptical inclusions in composites

A volume integral equation method is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solid containing interacting multiple anisotropic elliptical inclusions subject to uniform remote tension or remote in-plane shear. This method is applied to two-dimensional problems involving long parallel elliptical cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central elliptical inclusion is carried out for square and hexagonal packing of anisotropic inclusions. The effects of the number of anisotropic inclusions and various inclusion volume fractions on the stress field at the interface between the isotropic matrix and the central elliptical cylindrical inclusion are investigated in detail. The stress field at the interface between the isotropic matrix and the central elliptical inclusion is also compared with that between the isotropic matrix and the central circular inclusion.     

Dispersion characteristics of asymmetric coupled anisotropic dielectric waveguides using FDFD

The formulation is developed in the frequency domain and the finite difference method is used for the numerical solution of the scalar wave equation, written in terms of the transverse components of the magnetic field. As a result a conventional eigenvalue problem is obtained without the presence of spurious modes due to the implicit inclusion of the divergence of the magnetic field equal to zero. The formulation is developed to include biaxial anisotropic dielectrics with an index profile varying arbitrarily in the cross section of the waveguide under analysis. This formulation is then applied to the analysis of the influence on the dispersion characteristics of the dimensions of asymmetric coupled rectangular uniaxial anisotropic dielectric waveguides. As expected, the reduction of the height or the width of one of the rectangular dielectric waveguides causes the dispersion curves to move towards higher frequencies.     

Diffraction theory of an open resonator embedded in a uniaxial medium

A two-dimensional diffraction theory of an open strip resonator embedded in a umaxially anisotropic medium is presented. Applying the modified Kirchhoff approximation and a Green's function technique, the problem is reduced to the solution of a system of two coupled, homogencous integral equations for the unknown field distributions on the mirrors. The eigenvalues of this system of integral equations are related to the characteristic (complex) frequencies of the modes of the open resonator by an implicit relationship. The eigenvalues as well as the corresponding characteristic frequencies are calculated numerically. From the characteristic frequency of a mode, its resonance frequency, its diffraction loss and its quality factor follow immediatedly. For various resonator configurations the influence of the orientation of the optical axis of the uniaxial medium on the properties of the modes is studied in detail.     

Modelling of finite-size grating waveguide filters on the basis of a free-space diffraction approach

The field reflection coefficientR of a high strength and short length waveguide grating section is obtained by resorting to a free-space diffraction approach whereby the waveguide grating is excited by means of a virtual high index prism. The effective index of the guided mode, the coupling coefficient between contradirectionally propagating modes of the grating waveguide as well as the modal field are first obtained in the case of an infinite grating by interpreting the results given by an exact diffraction solution by means of a three-wave coupled mode formalism. These results are then introduced in the coupled wave analysis of the structure of finite extent to provide the field at the entrance of the waveguide grating section where the usual overlap integral with the incident field delivers the reflection coefficient. We neglect only the result of radiation mode propagation. This assumption can be made when the grating section is longer than the propagating length of all the leaky modes.     

Modelling of finite-size grating waveguide filters on the basis of a free-space diffraction approach

The field reflection coefficient R of a high strength and short length waveguide grating section is obtained by resorting to a free-space diffraction approach whereby the waveguide grating is excited by means of a virtual high index prism. The effective index of the guided mode, the coupling coefficient between contradirectionally propagating modes of the grating waveguide as well as the modal field are first obtained in the case of an infinite grating by interpreting the results given by an exact diffraction solution by means of a three-wave coupled mode formalism. These results are then introduced in the coupled wave analysis of the structure of finite extent to provide the field at the entrance of the waveguide grating section where the usual overlap integral with the incident field delivers the reflection coefficient. We neglect only the result of radiation mode propagation. This assumption can be made when the grating section is longer than the propagating length of all the leaky modes.     

Light diffraction effects from Kevlar aramid fibers

The light diffraction phenomenon from Kevlar 49 fibers exhibits unusual polarization effects and its angular and wavelength dependence cannot be described by the grating equation. This unexpected behavior is analyzed by modeling the fiber as an anisotropic grating of finite thickness in which the diffraction arises from a periodic change in the orientation of the molecular dipoles. Theoretical predictions on this basis explain all available experimental observations. These findings are discussed in regard to currently proposed structural models for the fiber.     

Directional properties of an array of sound receivers positioned in an impedance screen recess

Expressions for calculating the directional characteristics of an array of sound receivers positioned in a waveguide with impedance walls are obtained from the solution to the problem on the diffraction of a plane sound wave by the waveguide open end with impedance flanges. The waveguide can be of a finite length, and, in this case, it can be considered as an open cavity in an impedance screen. The solution of the integral equation for the sound pressure distribution over the opening area is reduced to the solution of an infinite system of algebraic equations for the coefficients of the field expansion in normal waveguide waves. Examples of calculated directional characteristics are presented for arrays with receivers positioned at different distances from the opening and for different values of the impedances of the waveguide walls and flanges.