Calculation of the band structure of semiconductor quantum wells using scattering matrices

The transfer-matrix method widely used in the calculation of the band structure of semiconductor quantum wells is found to have limitations due to its intrinsic numerical instability. It is pointed out that the numerical instability arises from free-propagating transfer matrices. A new scattering-matrix method is developed for the multiple-band Kane model within the envelope-function approximation. Compared with the transfer-matrix method, the proposed algorithm is found to be more efficient and stable. A four-band Kane model is used to check the validity of the method and the results are found to be in good agreement with earlier calculations.

On the soft wall guiding potentials in realistic quantum waveguides

A transfer matrix method is presented for the study of electron conduction in a quantum waveguide with soft wall lateral confinement. By transforming the two-dimensional Schrodinger equation into a set of second order ordinary differential equations, the total transfer matrix is obtained and the scattering probability amplitudes are calculated. The proposed method is applied to the evaluation of the electron transmission in two types of cavity structure with finite-height square-well confinement. The results obtained by our method, which are found to be in excellent agreement with those from another transfer matrix method, suggest that the infinite square-well potential is a good approximation to finite-height square-well confinement for electrons propagating in the ground transverse mode, but softening of the walls has an obvious effect on the electron transmission and mode-mixing for propagating in the excited transverse mode. (C) 1996 American Institute of Physics.

Theoretical Study on the Butt-coupling of Semiconductor Laser and Fiber Bragg Grating

The butt-coupling between a semiconductor laser diode and a fiber Bragg grating external cavity acts a key roll on the laser characteristics. The scatter matrix method considering the butt-coupling efficiency is used to analyze the butt-coupling between them. It is found that the butt-coupling distance and coupling efficiency determine the laser characteristics. For strong feedback, the single lasing wavelength changes in the reflection bandwidth of the effective reflectivity ( approximately the Bragg region of the fiber Bragg grating) as the distances change. For weak feedback condition, some different results are obtained. The SMSRs in the two conditions are presented and analyzed. These results can provide important design guidance of device parameters for the practical fabrication.

3μm GaInNAs/GasAs Quantum Well Resonant Cavity Enhanced Photodetector

A 1.3μm GaInNAs resonant cavity enhanced (RCE) photodetector (PD) has been grown by molecular beam epitaxy (MBE) monolithically on (100) GaAs substrate using a home-made ion-removed dc-plasma cell as nitrogen source. A transfer matrix method was used to optimize the device structure. The absorption region is composed of three GaInNAs quantum wells separated by GaAs layers. Devices were isolated by etching 130μm-diameter mesas and filling polyamide into grooves. The maximal quantum efficiency of the device is about 12% at 1.293μm. Full width at half maximum (FWHM) is 5.8nm and 3dB bandwidth is 304MHz. Dark current is 2 * 10~(-11) A at zero bias voltage. Further improvement of the performance of the RCE PD can be obtained by optimizing of the structure design and MBE growth conditions.

Evaluation of the Uncertainties Caused by the Forward Scattering in Turbidity Measurement of the Coagulation Rate

The forward scattering light (FSL) received by the detector can cause uncertainties in turbidity measurement of the coagulation rate of colloidal dispersion, and this effect becomes more significant for large particles. In this study, the effect of FSL is investigated on the basis of calculations using the T-matrix method, an exact technique for the computation of nonspherical scattering. The theoretical formulation and relevant numerical implementation for predicting the contribution of FSL in the turbidity measurement is presented. To quantitatively estimate the degree of the influence of FSL, an influence ratio comparing the contribution of FSL to the pure transmitted light in the turbidity measurement is introduced. The influence ratios evaluated under various parametric conditions and the relevant analyses provide a guideline for properly choosing particle size, measuring wavelength to minimize the effect of FSL in turbidity measurement of coagulation rate.

Theoretical calculation of photoionization cross sections of B-like ions: N2+, O3+ and F4+

There can be found some notable discrepancies with regard to the resonance structures when R-matrix calculations from the Opacity Project and other sources are compared with recent absolute experimental measurements of Bizau et al [Astron. Astrophts. 439 387 (2005)] for B-like ions N2+, O3+ and F4+. We performed close-coupling calculations based on the R-matrix formalism for the photoionizations of ions mentioned above both for the ground states and first excited states in the near threshold regions. The present results are compared with experimental ones given by Bizau et al and earlier theoretical ones. Excellent agreement is obtained between our theoretical results and the experimental photoionization cross sections. The present calculations show a significant improvement over the previous theoretical results.

Study of unbound nuclei N-11 resonance energy level

The excitation functions of elastic scattering proton which were measured with inverse kinematics of elastic resonance scattering reactions in GANIL and MSU have been fitted by the multi-energy level R-matrix theory. The final result shows that the new energy levels order for nucleus N-11 should be 1/2(+), 1/2(-), 5/2(+), 3/2(+), 3/2(-), 5/2(+), 7/2(-), which is consistent with the experimental results of Be-11 (the mirror nucleus of N-11) and the theoretical calculation of N-11 with GCM theory.

Effects of rubber content and temperature on dynamic fracture toughness of ABS materials

The effects of rubber content and temperature on dynamic fracture toughness of ABS materials have been investigated based on the J-integral and crack opening displacement (COD, delta) concepts by an instrumented Charpy impact test. A multiple specimens R-curve method and stop block technique are used. It is shown that the materials exhibit a different toughness behavior, depending on rubber content and temperature. The resistance against stable crack initiation (J(0.2) or delta(0.2)) increases with increasing rubber content. However, J(0.2) first increased with increasing temperature until reaching the maximum value; after that, it decreases with further increasing the temperature. (C) 2000 John Wiley & Sons, Inc.

The preparation and bioactivity research of agaro-oligosaccharides

Agaro-oligosaccharides were hydrolytically obtained from agar using hydrochloric acid, citric acid, and cationic exchange resin (solid acid). The FT-IR and NMR data showed that the hydrolysate has the structure of agaro-oligomers. Orthogonal matrix method was applied to optimize the preparation conditions based on alpha-naphthylamine end-labeled HPLC analysis method. The optimal way for oligosaccharides with different degree of polymerization (DP) was achieved by using solid acid degradation, which could give high yield and avoid solution neutralization process. Agaro-oligosaccharides with high purity were consequently obtained by activated carbon column isolation. Furthermore, the antioxidant and alpha-glucosidase inhibitory activity of three fractions were also investigated. The result indicated that 8% ethanol-eluted fraction showed highest activity against alpha-glucosidase with IC50 of 8.84 mg/mL, while 25% ethanol-eluted fraction possessed excellent antioxidant ability.

计算层状介质中电磁场的层矩阵法

As an important branch of electrical prospecting method, the artificial source frequency domain electromagnetism method has received more and more attention. But when conducts the fundamental research, people often isolated study some concrete method, so the research results of one method are very difficult to apply to another method directly. This article will possess the artificial source frequency domain EM method to an 1D model simply. It is stratified medium model, with an electric or magnetic source in or outside of it. Then take the horizontal electric dipole source as an example to introduce how to computing the EM field in stratified medium. Because layer matrix is the key of establishing equations, so we call it the layer-matrix method. The key of layer-matrix method is establishing equations by using layer matrixes in wavenumber(kx, ky, z) domain, then obtains the electromagnetic field value of wavenumber domain. After Fourier transform, we can get electromagnetic field of any position in spatial domain. The layer matrix technique theoretically can calculate electromagnetic field of any position for any source, is suitable for many kinds of electromagnetic method. After introduction of the layer matrix method, this article has done some CSAMT, MCSEM and Wireless Electro-Magnetic Method (WEM) modeling with layer matrix method separately. In CSAMT modeling, we get electromagnetic field dissemination characteristics considering wave number of the air, and obtain three-dimensional distribution characteristics of the electromagnetic field. In MCSEM modeling, we get electromagnetic field dissemination characteristics with and without considering the airwave, and obtain three-dimensional distribution characteristics of electromagnetic field. In WEM modeling, we get electromagnetic field’s difference between considering the ionosphere and not considering it, and recognize the ionosphere’s influence of electromagnetic field. With the layer matrix technique, we have got some new understandings of EM dissemination rules of different situations. All analysis results indicate that the layer-matrix technique is credible and effective, and are worthy of further thorough research and development.

中强地震震源的破裂过程

With the large developments of the seismic sources theory, computing technologies and survey instruments, we can model and rebuild the rupture process of earthquakes more realistically. On which earthquake sources' properties and tectonic activities law are realized more clearly. The researches in this domain have been done in this paper as follows. Based on the generalized ray method, expressions for displacement on the surface of a half-space due to an arbitrary oriented shear and tensile dislocation are also obtained. Kinematically, fault-normal motion is equivalent to tensile faulting. There is some evidence that such motion occurs in many earthquakes. The expressions for static displacements on the surface of a layered half-space due to static point moment tensor source are given in terms of the generalized reflection and transmission coefficient matrix method. The validity and precision of the new method is illustrated by comparing the consistency of our results with the analytical solution given by Okada's code employing same point source and homogenous half-space model. The computed vertical ground displacement using the moment tensor solution of the Lanchang_Gengma earthquake displays considerable difference with that of a double couple component .The effect of a soft layer at the top of the homogenous half-space on a shallow normal-faulting earthquake is also analyzed. Our results show that more seismic information would be obtained utilizing seismic moment tensor source and layered half-space model. The rupture process of 1999 Chi-Chi, Taiwan, earthquake investigated by using co-seismic surface displacement GPS observations and far field P-wave records. In according to the tectonic analysis and distributions of aftershock, we introduce a three-segment bending fault planes into our model. Both elastic half-space models and layered-earth models to invert the distribution of co-seismic slip along the Chi-Chi earthquake rupture. The results indicate that the shear slip model can not fit horizontal and vertical co-seismic displacements together, unless we add the fault-normal motion (tensile component) in inversions. And then, the Chi Chi earthquake rupture process was obtained by inversion using the seismograms and GPS observations. Fault normal motions determined by inversion, concentrate on the shallow northern bending fault from Fengyuan to Shuangji where the surface earthquake ruptures reveal more complexity and the developed flexural slip folding structures than the other portions of the rupture zone For understanding the perturbation of surface displacements caused by near-surface complex structures, We have taken a numeric test to synthesize and inverse the surface displacements for a pop-up structure that is composed of a main thrust and a back thrust. Our result indicates that the pop-up structure, the typical shallow complex rupture that occurred in the northern bending fault zone form Fengyuan to Shuangji, can be modeled better by a thrust fault added negative tensile component than by a simple thrust fault. We interpret the negative tensile distributions, that concentrate on the shallow northern bending fault from Fengyuan to Shuangji, as a the synthetic effect including the complexities of property and geometry of rupture. The earthquake rupture process also reveal the more spatial and temporal complexities form Fenyuan to SHuangji. According to the three-components teleseismic records, the S-wave velocity structure beneath the 59 teleseismic stations of Taiwan obtained by using the transform function method and the SA techniques. The integrated results, the 3D crustal structure of Taiwan reveal that the thickest part of crustal local in the western Central Range. This conclusion is consistent with the result form the Bouguer gravity anomaly. The orogenic evolution of Taiwan is young period, and the developing foot of Central Range dose not in static balancing. The crustal of Taiwan stays in the course of dynamic equilibrium. The rupture process of 2003)2,24,Jiashi, Xinjiang earthquake was estimated by the finite fault model using far field broadband P wave records of CDSN and IRIS. The results indicate that the earthquake focal is north dip trust fault including some left-lateral strike slip. The focal mechanism of this earthquake is different form that of earthquakes occurred in 1997 and 1998, but similar to that of 1996, Artux, Xinjiang earthquake. We interpreted that the earthquake caused trust fault due to the Tarim basin pushing northward and orogeny of Tianshan mountain. In the end, give a brief of future research subject: Building the Real Time Distribute System for rupture process of Large Earthquakes Based on Internet.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

Theory for determining alignment parameters of symmetric top molecule using (n+1) LIF

Expressions used for extracting the population and alignment parameters of a symmetric top molecule from (n + 1) laser-induced fluorescence (LIF) are derived by employing the tensor density matrix method. The molecular population and alignment are described by molecular state multipoles. The LIF intensity is a complex function of the initial molecular state multipoles, the dynamic factors, and the excitation-detection geometrical factors. The problem of how to extract the initial molecular state multipoles from (2 + 1) LIF, as an example, is discussed in detail. (C) 2000 American Institute of Physics. [S0021-9606(00)30744-9].

Nonlinear parameter retrieval from three- and four-wave mixing in metamaterials

We present a generalized nonlinear susceptibility retrieval method for metamaterials based on transfer matrices and valid in the nondepleted pump approximation. We construct a general formalism to describe the transfer matrix method for nonlinear media and apply it to the processes of three- and four-wave mixing. The accuracy of this approach is verified via finite element simulations. The method is then reversed to give a set of equations for retrieving the nonlinear susceptibility. Finally, we apply the proposed retrieval operation to a three-wave mixing transmission experiment performed on a varactor loaded split ring resonator metamaterial sample and find quantitative agreement with an analytical effective medium theory model. © 2010 The American Physical Society.