Application of fast Pade approximation in simulating photonic crystal nanocavities by FDTD technology

The exact calculation of mode quality factor Q is a key problem in the design of high-Q photonic crystal nanocavity. On the basis of further investigation on conventional Pade approximation, FDM and DFT, Pade approximation with Baker's algorithm is enhanced through introducing multiple frequency search and parabola interpolation. Though Pade approximation is a nonlinear signal processing method and only short time sequence is needed, we find the different length of sequence requirements for 2D and 3D FDTD, which is very important to obtain convergent and accurate results. By using the modified Pade approximation method and 3D FDTD, the 2D slab photonic crystal nanocavity is analyzed and high-Q multimode can be solved quickly instead of large range high-resolution scanning. Monitor position has also been investigated. These results are very helpful to the design of photonic crystal nanocavity devices. (C) 2008 Elsevier B.V. All rights reserved.

Calculation of light delay for coupled microrings by FDTD technique and Pade approximation

The Pade approximation with Baker's algorithm is compared with the least-squares Prony method and the generalized pencil-of-functions (GPOF) method for calculating mode frequencies and mode Q factors for coupled optical microdisks by FDTD technique. Comparisons of intensity spectra and the corresponding mode frequencies and Q factors show that the Pade approximation can yield more stable results than the Prony and the GPOF methods, especially the intensity spectrum. The results of the Prony method and the GPOF method are greatly influenced by the selected number of resonant modes, which need to be optimized during the data processing, in addition to the length of the time response signal. Furthermore, the Pade approximation is applied to calculate light delay for embedded microring resonators from complex transmission spectra obtained by the Pade approximation from a FDTD output. The Prony and the GPOF methods cannot be applied to calculate the transmission spectra, because the transmission signal obtained by the FDTD simulation cannot be expressed as a sum of damped complex exponentials. (C) 2009 Optical Society of America

Calculation of propagation loss in photonic crystal waveguides by FDTD technique and Pade approximation

The propagation losses in single-line defect waveguides in a two-dimensional (2D) square-lattice photonic crystal (PC) consisted of infinite dielectric rods and a triangular-lattice photonic crystal slab with air holes are studied by finite-difference time-domain (FDTD) technique and a Pade approximation. The decaying constant beta of the fundamental guided mode is calculated from the mode frequency, the quality factor (Q-factor) and the group velocity v(g) as beta = omega/(2Qv(g)). In the 2D square-lattice photonic crystal waveguide (PCW), the decaying rate ranged from 10(3) to 10(-4) cm(-1) can be reliably obtained from 8 x 10(3)-item FDTD output with the FDTD computing time of 0.386 ps. And at most 1 ps is required for the mode with the Q-factor of 4 x 10(11) and the decaying rate of 10(-7) cm(-1). In the triangular-lattice photonic crystal slab, a 10(4)-item FDTD output is required to obtain a reliable spectrum with the Q-factor of 2.5 x 10(8) and the decaying rate of 0.05 cm(-1). (c) 2004 Elsevier B.V. All rights reserved.

Computation of resonant frequencies and quality factors of cavities by FDTD technique and Pade approximation

The finite-difference time domain (FDTD) technique and the Pade approximation with Baker's algorithm are used to calculate the mode frequencies and quality factors of cavities. Comparing with the fast Fourier transformation/Pade method, we find that the Fade approximation and the Baker's algorithm can obtain exact resonant frequencies and quality factors based on a much shorter time record of the FDTD output.

Resonant frequencies and quality factors for optical equilateral triangle resonators calculated by FDTD technique and the Pade approximation

The mode wavelength and quality factor (Q-factor) for resonant modes in optical equilateral triangle resonators (ETR's) are calculated by the finite-difference time-domain (FDTD) technique and the Pade approximation, For an ETR with the side length of 3 mu m and the refractive index of 3.2, we get the mode wavelength interval of about 70 nm and the Q-factor of the fundamental mode over 10(3), The results show that the ETR is suitable to realize single-mode operation, and that the radiation loss in the corner regions of ETR is rather low, In addition, the numerical results of the mode wavelength agree very well with our analytical formula.

Calculation of Light Delay by FDTD Technique and Pade Approximation

The time delay for light transmission in a coupled microring waveguide structure is calculated from the phase shift of the transmission coefficient obtained by Pade approximation with Baker's algorithm from FDTD Output. The results show that the Pade approximation is a powerful tool for saving time in FDTD simulation.

Applicatioo of Pade Approximation in Simulating Photonic Crystals

To save finite-difference time-domain(FDTD) computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Padé approximation with Baker's algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000,the intensity spectrum obtained by the Padé approximation from a 28-item sequence output is more exact than that obtained by fast Fourier transformation from a 220-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method,and then the band diagrams are obatined. In addition,mode frequencies and Q-factors are calculated for photonic crystal microcavity.

High order unconditionally stable difference schemes for the Riesz space-fractional telegraph equation

In this paper, a class of unconditionally stable difference schemes based on the Pad´e approximation is presented for the Riesz space-fractional telegraph equation. Firstly, we introduce a new variable to transform the original dfferential equation to an equivalent differential equation system. Then, we apply a second order fractional central difference scheme to discretise the Riesz space-fractional operator. Finally, we use (1, 1), (2, 2) and (3, 3) Pad´e approximations to give a fully discrete difference scheme for the resulting linear system of ordinary differential equations. Matrix analysis is used to show the unconditional stability of the proposed algorithms. Two examples with known exact solutions are chosen to assess the proposed difference schemes. Numerical results demonstrate that these schemes provide accurate and efficient methods for solving a space-fractional hyperbolic equation.

Expansion of high pressure gas into air - A more realistic blast wave model

In this paper, we consider a more realistic model of a spherical blast wave of moderate strength. An arbitrary number of terms for the series solution in each of the regions behind the main shock - the expansion region, the nearly uniform region outside the main expansion and the region between the contact surface and the main shock, have been generated and matched across the boundaries. We then study the convergence of the solution by using Pade approximation. It constitutes a genuine analytic solution for a moderately strong explosion, which, however, does not involve a secondary shock. The pressure distribution behind the shock however shows some significant changes in the location of the tail of the rarefaction and the interface, in comparison to the planar problem. The theory developed for the spherical blasts is also extended to cylindrical blasts. The results are compared with the numerical solution.

Analysis of mode characteristics for equilateral triangle semiconductor microlasers with imperfect boundaries

The influence of imperfect boundaries on the mode quality factor is investigated for equilateral-triangle-resonator (ETR) semiconductor microlasers by the finite difference time domain technique and the Pade approximation with Baker's algorithm. For 2-D ETR with a refractive index of 3.2 and side length of 5 mum, the confined modes can still have a quality factor of about 1000 as small triangles with side length of 1 mum are cut from the vertices of the ETR. For a deformed 5 mum ETR with round vertices and curve sides, the simulated mode quality factors are comparable to the measured results.

Mode quality factor based on far-field emission for square resonators

The quality factors of modes in square resonators are calculated based on the far-field emission of the analytical field distribution. The obtained quality factors are in reasonable agreement with those calculated by the finite-difference time-domain (FDTD) technique and Pade approximation method. The emission power in the square diagonal directions for whispering-gallery-like modes in square resonators is zero due to the interference cancellation caused by the odd field distributions relative to the diagonal mirror planes, so they have larger quality factors than the modes with even field distribution.

Comparison of free spectral range and quality factor for two-dimensional square and circular microcavities

Free spectral range of whispering-gallery (WG)-like modes in a two-dimensional (2D) square microcavity is found to be twice that in a 2D circular microcavity. The quality factor of the WG-like mode with the low mode number in a 2D square microcavity, calculated by the finite-difference time-domain (FDTD) technique and the Pade approximation method, is found to exceed that of the WG mode in 2D circular microcavity with the same cavity dimension and close mode wavelength.

Modes in square resonators

Modes in square resonators are analyzed and classified according to the irreducible representations of the point group C-4v. If the mode numbers p and q that denote the number of wave nodes in the directions of two orthogonal square sides are unequal and have the same even-odd characteristics, the corresponding double modes are accidentally degenerate and can be combined into two new distributions with definite parities relative to the square diagonal mirror planes. The distributions with odd parities belong to the whispering-gallery-like modes in square resonators. The mode frequencies and quality factors are also calculated by the finite-difference time-domain technique and Pade approximation method. The numerically calculated mode frequencies agree with the theoretical ones very well and the whispering-gallery-like modes have quality factors much higher than other modes, including their accidentally degenerate counterparts in square resonators.

Mode characteristics for equilateral triangle optical resonators

Mode characteristics of equilateral triangle resonators (ETRs) are analyzed based on the symmetry operation of the point group C-3v. The results show that doubly degenerate eigenstates can be reduced to the A(1) and A(2) representations of C-3v, if the longitudinal mode number is a multiple of 6; otherwise, they form the E irreducible representation Of C-3v. And the one-period length for the mode light ray is half of the perimeter of the ETR. Mode Q-factors are calculated by the finite-difference time-domain (FDTD) technique and compared with those calculated from far-field emission based on the analytical near-field pattern for TE and TM modes. The results show that the far-field emission based on the analytical field distribution can be used to estimate the mode Q-factor, especially for TM modes. FDTD numerical results also show that Q-factor of TE modes reaches maximum value as the longitudinal mode number is a multiple of 7. In addition, photoluminescence spectra and measured Q-factors are presented for fabricated ETR with side lengths of 20 and 30 mu m, and the mode wavelength intervals are compared with the analytical results.

Analysis of mode quality factors and mode reflectivities for nanowire cavity by FDTD technique

The mode frequency and the quality factor of nanowire cavities are calculated from the intensity spectrum obtained by the finite-difference time-domain (FDTD) technique and the Pade approximation. In a free-standing nanowire cavity with dielectric constant epsilon = 6.0 and a length of 5 mu m, quality factors of 130, 159, and 151 are obtained for the HE11 modes with a wavelength around 375 nm, at cavity radius of 60, 75, and 90 nm, respectively. The corresponding quality factors reduce to 78, 94, and 86 for a nanowire cavity standing on a sapphire substrate with a refractive index of 1.8. The mode quality factors are also calculated for the TE01 and TM01 modes, and the mode reflectivities are calculated from the mode quality factors.