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

Design of shallow acceptors in ZnO through compensated donor-acceptor complexes: A density functional calculation

The intrinsic large electronegativity of O 2p character of the valence-band maximum (VBM) of ZnO renders it extremely difficult to be doped p type. We show from density functional calculation that such VBM characteristic can be altered by compensated donor-acceptor pairs, thus improve the p-type dopability. By incorporating (Ti+C) or (Zr+C) into ZnO simultaneously, a fully occupied impurity band that has the C 2p character is created above the VBM of host ZnO. Subsequent doping by N in ZnO: (Ti+C) and ZnO: (Zr+C) lead to the acceptor ionization energies of 0.18 and 0.13 eV, respectively, which is about 200 meV lower than it is in pure ZnO.

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.

Calculation of energy levels in InGaAs/GaAs quantum dot array

The subbands of the ground state E-c1, the first excited state E-c2 and heavy hole state E-HH1 are calculated by solving the eigenvalues of effective-mass Hamiltonian H-0 which is derived from eight-band k . p theory and the calculations are performed at k(x) = k, = k = 0 for the three-dimensional array of InGaAs/GaAs quantum dots (QDs). With indium content in InGaAs QDs gradually increasing from 30% to 100%,the intersubband transition wavelength of E-c2 to E-c1, blue-shifts from 18.50 to 11.87 mu m,while the transition wavelength of E-c1, to E-HH1, red-shifts from 1. 04 to 1. 73 mu m. With the sizes of Ir-0.5 Ga-0.5 As and InAs QDs increasing from 1.0 to 5.0 nm, the intersubband transition from E-c1, to E-C2 transforms from bound-state-to-continuum-state to bound-state-to-bound-state, and the corresponding intersubband transition wavelengths red-shift from 8.12 pm (5.90 pm) to 53.47 mu m (31.87 pm), respectively, and the transition wavelengths of E-C1 to E-HH1 red-shift from 1. 13 mu m (1.60 mu m) to 1.27 mu m (2.01 mu m), respectively.

Calculation of the current noise spectrum in mesoscopic transport: A quantum master equation approach

Based on our recent work on quantum transport [X. Q. Li , Phys. Rev. B 71, 205304 (2005)], we show how an efficient calculation can be performed for the current noise spectrum. Compared to the classical rate equation or the quantum trajectory method, the proposed approach is capable of tackling both the many-body Coulomb interaction and quantum coherence on an equal footing. The practical applications are illustrated by transport through quantum dots. We find that this alternative approach is in a certain sense simpler and more straightforward than the well-known Landauer-Buttiker scattering matrix theory.

Application of plane wave method to the calculation of electronic states of nano-structures

The electronic states of nano-structures are studied in the framework of effective-mass envelope-function theory using the plane wave basis. The barrier width and the number of plane waves are proposed to be 2.5 times the effective Bohr radius and 15(n), respectively, for n-dimensional nano-structures (n = 1,2,3). Our proposals can be widely applied in the design of various nano-structure devices.

Calculation of the hourly and daily radiation incident on three step tracking planes

Mathematical formulas for estimating the hourly and daily radiation incident on planes of azimuth three step tracking and hour angle three step tracking have been derived in this paper. Based on the hourly solar radiation data of an average day in each month at Er-Lian-Hao-Te city, the hourly and monthly radiation received by planes of these two kinds of tracking have been calculated. The results show that in this district, one axis azimuth three step tracking and hour angle three step tracking could, respectively, obtain 66.5% and 63.3% higher radiation than that on the horizontal surface all year. Moreover, a two axis azimuth three step tracking plane could receive 72% more radiation than the horizontal surface. (C) 2002 Elsevier Science Ltd. All rights reserved.

The content calculation of hexagonal phase inclusions in cubic GaN films on GaAs(001) substrates grown by metalorganic chemical vapor deposition

Cubic GaN(c-GaN) films are grown on GaAs(001) substrates by metalorganic chemical vapor deposition (MOCVD). Two GaN samples were grown with different buffer layer, the deposition time of each was 1 and 3 min, respectively. 4-circle X-ray double crystal diffraction (XRDCD) was used to study the secondary crystallographic phases presented in the c-GaN films. The phase composition of the epilayers was determined by X-ray reciprocal space mapping. The intensities of the c-GaN(002) and h-GaN(10 (1) over bar 1) planes detected in the mapping were investigated by omega scans. The content of the hexagonal phase inclusions in the c-GaN films was calculated to about 1.6 and 7.9%, respectively. The thicker buffer layer is not preferable for growing high quality pure c-GaN films. (C) 2000 Elsevier Science S.A. All rights reserved.

Ground-state properties of an artificial molecule: a simple model calculation for the double quantum dot

The ground state of a double quantum-dot structure is studied by a simplified Anderson-type model. Numerical calculations reveal that the ground-state level of this artificial molecule increases with the increasing single particle level of the dot, and also increases with the decreasing transfer integrals. We show the staircase feature of the electron occupation and the properties of the ground-state eigenvector by varying the;single particle level of the dot.

Calculation of valence subband structures for strained GaInP/AlGaInP quantum wells without axial approximation

Usually in the calculation of valence subband structure for III-V direct bandgap material, axial approximation had been used in the Luttinger-Kohn model to simplify the computational efforts. In this letter, the valence subband structure for the GaInP/AlGaInP strained and lattice-matched quantum wells was calculated without axial approximation, on the basis of 6x6 Luttinger-Kohn Hamiltonian including strain and spin-orbit splitting effects. The numerical simulation results were presented with help of the finite-difference methods. The calculation results with/without axial approximation were compared and the effect of axial approximation on the valence subband structure was discussed in detail. The results indicated that there was a strong warping in the GaInP valence band, and axial approximation can lead to an error when k was not equal to zero, especially for compressively strained and lattice-matched GaInP/AlGaInP quantum wells.

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.

Fast calculation technique for scattering in T-matrix method

In scattering calculations using the T-matrix method, the calculation of the T-matrix involves multiplication and inversion of matrices. These two types of matrix operations are time-consuming, especially for the matrices with large size. Petrov et al. [D. Petrov, Y. Shkuratov, G. Videen, Opt. Lett. 32 (2007) 1168] proposed an optimized matrix inversion technique, which suggests the inversion of two matrices, each of which contains half the number of rows. This technique reduces time-consumption significantly. On the basis of this approach, we propose another fast calculation technique for scattering in the T-matrix method, which obtains the scattered fields through carrying out only the operations between matrices and the incident field coefficient. Numerical results show that this technique can decrease time-consumption by more than half that of the optimized matrix inversion technique by Petrov et al. (c) 2008 Elsevier B.V. All rights reserved.