基于严格耦合波理论的多层介质膜光栅衍射特性分析

基于严格耦合波理论建立了多层介质膜光栅的衍射机理模型,给出了TE波自准直条件下多层介质膜光栅衍射效率的表达式.以-1级衍射效率为评价函数,分析了表面浮雕结构分别为HfO2和SiO2材料的介质膜光栅获得衍射效率优于96%的结构参数.数值计算表明,顶层材料为HfO2的介质膜光栅具有更宽的结构选择范围.最后分析了介质膜光栅的制备容差和允许的入射角度范围.

High efficiency and broad bandwidth grating coupler between nanophotonic waveguide and fibre

A high efficiency and broad bandwidth grating coupler between a silicon-on-insulator (SOI) nanophotonic waveguide and fibre is designed and fabricated. Coupling efficiencies of 46% and 25% at a wavelength of 1.55 mu m are achieved by simulation and experiment, respectively. An optical 3 dB bandwidth of 45 nm from 1530 nm to 1575 nm is also obtained in experiment. Numerical calculation shows that a tolerance to fabrication error of 10 nm in etch depth is achievable. The measurement results indicate that the alignment error of +/-2 mu m results in less than 1 dB additional coupling loss.

Direct detection of the relative strength of Rashba and Dresselhaus spin-orbit interaction: Utilizing the SU(2) symmetry

We propose a simple method to detect the relative strength of Rashba and Dresselhaus spin-orbit interactions in quantum wells (QWs) without relying on the directional-dependent physical quantities. This method utilizes the two different critical gate voltages that leading to the remarkable signals of SU(2) symmetry, which happens to reflect the intrinsic-structure-inversion asymmetry of the QW. We support our proposal by the numerical calculation of in-plane relaxation times based on the self-consistent eight-band Kane model. We find that the two different critical gate voltages leading to the maximum spin-relaxation times [one effect of the SU(2) symmetry] can simply determine the ratio of the coefficients of Rashba and Dresselhaus terms. Our proposal can also be generalized to extract the relative strengths of the spin-orbit interactions in quantum-wire and quantum-dot structures.

Effect of defects on optical phonon Raman spectra in SiC nanorods

Fabricated one-dimensional (1D) materials often have abundant structural defects. Experimental observation and numerical calculation indicate that the broken translation symmetry due to structural defects may play a more important role than the quantum confinement effect in the Raman features of optical phonons in polar semiconductor quantum wires such as SiC nanorods, (C) 1999 Elsevier Science Ltd. All rights reserved.

Intrasubband and intersubband transitions in lightly and heavily doped GaAs/AlGa1-xAs multiple quantum wells

We use a polarizer to investigate quantum-well infrared absorption, and report experimental results as follows. The intrasubband transition was observed in GaAs/AlxGa1-xAs multiple quantum wells (MQWs) when the incident infrared radiation (IR) is polarized parallel to the MQW plane. According to the selection rule, an intrasubband transition is forbidden. Up to now, most studies have only observed the intersubband transition between two states with opposite parity. However, our experiment shows not only the intersubband transitions, but also the intrasubband transitions. In our study, we also found that for light doping in the well (4x10(18) cm(-3)), the intrasubband transition occurs only in the lowest subband, while for the heavy doping (8x10(18) cm(-3)), such a transition occurs not only in the lowest subband, but also in the first excited one, because of the electron subband filling. Further experimental results show a linear dependence of the intrasubband transition frequency on the root of the well doping density. These data are in good agreement with our numerical results. Thus we strongly suggest that such a transition can be attributed to plasma oscillation. Conversely, when the incident IR is polarized perpendicular to the MQW plane, intersubband-transition-induced signals appear, while the intrasubband-transition-induced spectra disappear for both light and heavy well dopings. A depolarization blueshift was also taken into account to evaluate the intersubband transition spectra at different well dopings. Furthermore, we performed a deep-level transient spectroscopy (DLTS) measurement to determine the subband energies at different well dopings. A good agreement between DLTS, infrared absorption, and numerical calculation was obtained. In our experiment, two important phenomena are noteworthy: (1) The polarized absorbance is one order of magnitude higher than the unpolarized spectra. This puzzling result is well explained in detail. (2) When the IR, polarized perpendicular to the well plane, normally irradiates the 45 degrees-beveled edge of the samples, we only observed intersubband transition spectra. However, the intrasubband transition signals caused by the in-plane electric-field component are significantly absent. The reason is that such in-plane electric-field components can cancel each other out everywhere during the light propagating in the samples. The spectral widths of bound-to-bound and bound-to-continuum transitions were also discussed, and quantitatively compared to the relaxation time tau, which is deduced from the electron mobility. The relaxation times deduced from spectral widths of bound-to-bound and bound-to-continuum transitions are also discussed, and quantitatively compared to the relaxation time deduced from electron mobility. [S0163-1829(98)01912-2].

Observation of the third subband population in modulation-doped InGaAs/InAlAs heterostructure

The population of the third (n = 3) two-dimensional electron subband of InGaAs/InAlAs modulation-doped structures has been observed by means of Fourier transform photoluminescence (PL). Three well resolved PL peaks centred at 0.737, 0.908, and 0.980eV are observed, which are attributed to the transitions from the lowest three electron subbands to the n = 1 heavy-hole subband. The subband separations clearly exhibiting the features of the stepped quantum well with triangle and square potentials are consistent with numerical calculation. Thanks to the presence of Fermi cutoff, the population ratio of these three subbands can be estimated. Temperature- and excitation-dependent luminescences are also analyzed.

The third subband population in modulation-doped InGaAs/InAlAs heterostructures

We have observed the population of the third (n=3) two-dimensional electron subband of InGaAs/ InAlAs modulation-doped structures with very dense sheet carrier density by means of Fourier transform photoluminescence (PL). Three well-resolved PL peaks centered at 0.737, 0.908, and 0.980 eV are observed, which are attributed to the recombination transitions from the lowest three electron subbands to the n=1 heavy-hole subband. The subband separations clearly exhibit the features of the stepped quantum well with triangle and square potential, consistent with numerical calculation. Thanks to the presence of the Fermi cutoff, the population ratio of these three subbands can be estimated. Temperature and excitation intensity dependence of the quantum well luminescence intensity is also analyzed. (C) 1997 American Institute of Physics.

Mode field diameter and nonlinear properties of air-core nanowires

Exact solutions of Maxwell's equations describing the lightwave through 3-layer-structured cylindrical waveguide are obtained and the mode field diameter and nonlinear coefficient of air-core nanowires (ACNWs) are numerically calculated. The simulation results show that ACNWs offer some unique optical properties, such as tight field confining ability and extremely high nonlinearity. At a certain wavelength and air core radius, we optimize the waveguide design to maximize the nonlinear coefficient and minimize the mode field diameter. Our results show that the ACNWs may be powerful potential tools for novel micro-photonic devices in the near future.

Study on tapered crossed subwavelength gratings by Fourier modal method

Fourier modal method incorporating staircase approximation is used to study tapered crossed subwavelength gratings in this paper. Three intuitive formulations of eigenvalue functions originating from the prototype are presented, and their convergences are compared through numerical calculation. One of them is found to be suitable in modeling the diffraction efficiency of the circular tapered crossed subwavelength gratings without high absorption, and staircase approximation is further proven valid for non-highly-absorption tapered gratings. This approach is used to simulate the "moth-eye" antireflection surface on silicon, and the numerical result agrees well with the experimental one.

Simulation of Banded Spherulite Patterns by Coupled Logistic Map Lattice Model in Three Dimensions

Banded spherulite patterns are simulated in three dimensions by means of a Coupled Logistic map lattice model. The patterns obtained by numerical calculation are consistent with those in experiments. The simulation results also indicate that the hand spacing is decreased with the increase of parameter mu in the Logistic map and increased with the increase of the coupling parameter e for cube lattices, and increased with the increase of the thickness of the lattice for polymer film, which is quite similar to the results in some experiments. Spiral pattern in three dimensions is also shown in this paper, which helps us understand the form of banded spherulite in polymers.

Investigation of the gel effect in reactive extrusion processes for free radical polymerization

The gel effect in the reactive extrusion process for free radical polymerization in a closely intermeshing co-rotating twin screw extruder was investigated. First the reaction kinetic model was constructed mainly on the basis of entanglement theory. Next, numerical calculation expressions for the initiator and monomer concentrations, monomer conversion, average molecular weight and apparent viscosity were deduced. Finally, the evolution of the above variables were shown and discussed for the example of butyl methacrylate. The simulated results of the monomer conversion are in good agreement with experimental results.

层析静校正技术研究和应用

Static correction is one of the indispensable steps in the conventional onshore seismic data processing, particularly in the western part of China; it is theoretically and practically significant to resolve the issue of static correction. Conventional refraction static correction is put forward under the assumption that layered medium is horizontal and evenly distributed. The complicated nature of the near surface from western part of China is far from the assumption. Therefore, the essential way to resolve the static correction problem from the complex area is to develop a new theory. In this paper, a high-precision non-linear first arrival tomography is applied to solve the problem, it moved beyond the conventional refraction algorithm based on the layered medium and can be used to modeling the complex near surface. Some of the new and creative work done is as follows: One. In the process of first arrival tomographic image modeling, a fast high-order step algorithm is used to calculate the travel time for first arrival and ray path and various factors concerning the fast step ray tracing algorithm is analyzed. Then the second-order and third-order differential format is applied to the step algorithm which greatly increased the calculation precision of the ray tracing and there is no constraint to the velocity distribution from the complex areas. This method has very strong adaptability and it can meet the needs of great velocity variation from the complicated areas. Based on the numerical calculation, a fast high-order step is a fast, non-conditional and stable high-precision tomographic modeling algorithm. Two, in the tomographic inversion, due to the uneven fold coverage and insufficient information, the inversion result is unstable and less reliable. In the paper, wavelet transform is applied to the tomographic inversion which has achieved a good result. Based on the result of the inversion from the real data, wavelet tomographic inversion has increased the reliability and stability of the inversion. Three. Apply the constrained high-precision wavelet tomographic image to the static correction processing from the complex area. During tomographic imaging, by using uphole survey, refraction shooting or other weathering layer method, weathering layer can be identified before the image. Because the group interval for the shot first arrival is relatively big, there is a lack of precision for the near surface inversion. In this paper, an inversion method of the layer constraint and well constraint is put forward, which can be used to compensate the shallow velocity of the inversion for the shot first arrival and increase the precision of the tomographic inversion. Key words: Tomography ，Fast marching method，Wavelet transform, Static corrections, First break

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

温度、压力耦合背景下多相流体运移的数学模型研究

The Mathematical modeling of multiphase fluid flow is an important aspect of basin simulation, and also is a topic of geological frontier. Based on coupling relation of temperature, pressure and fluid flow, this dissertation discusses the modeling which conform to geological regularities of fluid migration. The modeling that is multi-field and multiphase includes heat transport equation, pressure evolvement equation, solution transport equation and fluid transport equation. The finite element method is effective numerical calculation methods. Author applies it to solve modeling and implements the finite element program, and the modeling is applied to Ying-Qiong Basin. The channels of fluid vertical migration are fault, fracture and other high penetrability area. In this thesis, parallel fracture model and columnar channel model have been discussed, and a characteristic time content and a characteristic space content been obtained to illustrate the influences of stratigraphic and hydrodynamic factors on the process. The elliptoid fracture model is established and its approximately solution in theory is gotten. Three kinds of modeling are applied to analyze the transient variation process of fluid pressure in the connected permeable formations. The elliptoid fracture model is the most similar geology model comparing with the other fracture models so the research on this fracture model can enhance the understanding to fluid pressure. In the non-hydrodynamic condition, because of the difference between water density and nature gas density, nature gas can migrate upon by float force. A one-dimension mathematical model of nature gas migration by float force is established and also applied to analyze the change in the saturation of gas. In the process of gas migration its saturation is non-continuous. Fluid flow is an important factor which influences the distribution of the temperature-field, the change of temperature can influence fluid property (including density, viscidity, and solubility),a nd the temperature field has coupling relations to the fluid pressure field. In this dissertation one-dimension and two-dimension thermal convection modeling is developed and also applied to analyze convective and conductive heat transfer. Author has established one-dimension and two-dimension mathematical modeling in which fluid is a mixture of water and nature gas based on the coupling relation between temperature and pressure, discussed mixture fluid convection heat transfer in different gas saturation, and analyzed overpressure form mechanism. Based on geothermal abnormity and pore pressure distribution in Dongfong 1-1, Yinggehai Basin, South China Sea, one-dimension mathematical modeling of coupling temperature and pressure is established. The modeling simulates the process that fluid migrates from deep to shallow and overpressure forms in shallow. When overpressure is so large that fractures appear and overpressure is released. As deep fluid flow to shallow, the high geothermal then forms in shallow. Based on the geological characteristics in Ya13-1, two-dimension mathematical modeling of coupling temperature and pressure is established. Fluid vertically flows in fault and then laterally migrates in reservoir. The modeling simulates the geothermal abnormity and pore pressure distribution in reservoir.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.