The influence of 1 nm AlN interlayer on properties of the Al0.3Ga0.7N/AlN/GaN HEMT structure

The sheet carrier concentrations, conduction band profiles and amount of free carriers in the barriers have been determined by solving coupled Schrodinger and Poisson equation self-consistently for coherently grown Al0.3Ga0.7N/GaN and Al0.3Ga0.7N/AlN/GaN structures on thick GaN. The Al0.3Ga0.7N/GaN heterojunction structures with and without 1 nm AlN interlayer have been grown by MOCVD on sapphire substrate, the physical properties for these two structures have been investigated by various instruments such as Hall measurement and X-ray diffraction. By comparison of the theoretical and experimental results, we demonstrate that the sheet carrier concentration and the electrons mobility would be improved by the introduction of an AlN interlayer for Al0.3Ga0.7N/GaN structure. Mechanisms for the increasing of the sheet carrier concentration and the electrons mobility will be discussed in this paper. (C) 2007 Elsevier Ltd. All rights reserved.

Properties of AlyGa1-yN/AlxGa1-xN/AlN/GaN Double-Barrier High Electron Mobility Transistor Structure

Electrical properties of AlyGa1-yN/AlxGa1-xN/AlN/GaN structure are investigated by solving coupled Schrodinger and Poisson equation self-consistently. Our calculations show that the two-dimensional electron gas (2DEG) density will decrease with the thickness of the second barrier (AlyGa1-yN) once the AlN content of the second barrier is smaller than a critical value y(c), and will increase with the thickness of the second barrier (AlyGa1-yN) when the critical AlN content of the second barrier y(c) is exceeded. Our calculations also show that the critical AlN content of the second barrier y(c) will increase with the AlN content and the thickness of the first barrier layer (AlxGa1-xN).

Effect of spontaneous and piezoelectric polarization on intersubband transition in AlxGa1-xN-GaN quantum well

A self-consistent solution of conduction band profile and subband energies for AlxGa1-xN-GaN quantum well is presented by solving the Schrodinger and Poisson equations. A new method is introduced to deal with the accumulation of the immobile charges at the AlxGa1-xN-GaN interface caused by spontaneous and piezoelectric polarization in the process of solving the Poisson equation. The effect of spontaneous and piezoelectric polarization is taken into account in the calculation. It also includes the effect of exchange-correlation to the one electron potential on the Coulomb interaction. Our analysis is based on the one electron effective-mass approximation and charge conservation condition. Based on this model, the electron wave functions and the conduction band structure are derived. We calculate the intersubband transition wavelength lambda(21) for different Al molar fraction of barrier and thickness of well. The calculated result can fit to the experimental data well. The dependence of the absorption coefficient a on the well width and the doping density is also investigated theoretically. (C) 2004 American Vacuum Society.

Observations on subband electron properties in In0.65Ga0.35As/In0.52Al0.48As MM-HEMT with Si delta-doped on the barriers

Magneto-transport measurements have been carried out on a Si delta-doped In0.65Ga0.35As/In0.52Al0.48As metamorphic high-electron-mobility transistor with InP substrate in a temperature range between 1.5 and 60 K under magnetic field up to 13 T. We studied the Shubnikov-de Haas (SdH) effect and the Hall effect for the In0.65Ga0.35As/In0.52Al0.48As single quantum well occupied by two subbands and obtained the electron concentration and energy levels respectively. We solve the Schrodinger-Kohn-Sham equation in conjunction with the Poisson equation self-consistently and obtain the configuration of conduction band, the distribution of carriers concentration, the energy level of every subband and the Fermi energy. The calculational results are well consistent with the results of experiments. Both experimental and calculational results indicate that almost all of the delta-doped electrons transfer into the quantum well in the temperature range between 1.5 and 60 K.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Quantum Mechanical Study on Tunnelling and Ballistic Transport of Nanometer Si MOSFETs

Using self-consistent calculations of million-atom Schrodinger-Poisson equations, we investigate the I-V characteristics of tunnelling and ballistic transport of nanometer metal oxide semiconductor field effect transistors (MOSFET) based on a full 3-D quantum mechanical simulation under nonequilibtium condition. Atomistic empirical pseudopotentials are used to describe the device Hamiltonian and the underlying bulk band structure. We find that the ballistic transport dominates the I-V characteristics, whereas the effects of tunnelling cannot be neglected with the maximal value up to 0.8mA/mu m when the channel length of MOSFET scales down to 25 nm. The effects of tunnelling transport lower the threshold voltage V-t. The ballistic current based on fully 3-D quantum mechanical simulation is relatively large and has small on-off ratio compared with results derived from the calculation methods of Luo et al.

Multiple valley couplings in nanometer Si metal-oxide-semiconductor field-effect transistors

We investigate the couplings between different energy band valleys in a metal-oxide-semiconductor field-effect transistor (MOSFET) device using self-consistent calculations of million-atom Schrodinger-Poisson equations. Atomistic empirical pseudopotentials are used to describe the device Hamiltonian and the underlying bulk band structure. The MOSFET device is under nonequilibrium condition with a source-drain bias up to 2 V and a gate potential close to the threshold potential. We find that all the intervalley couplings are small, with the coupling constants less than 3 meV. As a result, the system eigenstates derived from different bulk valleys can be calculated separately. This will significantly reduce the simulation time because the diagonalization of the Hamiltonian matrix scales as the third power of the total number of basis functions. (C) 2008 American Institute of Physics.

Spin states in InAs/AlSb/GaSb semiconductor quantum wells

We investigate theoretically the spin states in InAs/AlSb/GaSb broken-gap quantum wells by solving the Kane model and the Poisson equation self-consistently. The spin states in InAs/AlSb/GaSb quantum wells are quite different from those obtained by the single-band Rashba model due to the electron-hole hybridization. The Rashba spin splitting of the lowest conduction subband shows an oscillating behavior. The D'yakonov-Perel' spin-relaxation time shows several peaks with increasing the Fermi wave vector. By inserting an AlSb barrier between the InAs and GaSb layers, the hybridization can be greatly reduced. Consequently, the spin orientation, the spin splitting, and the D'yakonov-Perel' spin-relaxation time can be tuned significantly by changing the thickness of the AlSb barrier.

Dyakonov-Perel spin relaxation in InSb/AlxIn(1-x)Sb quantum wells

We investigate theoretically the Dyakonov-Perel spin relaxation time by solving the eight-band Kane model and Poisson equation self-consistently. Our results show distinct behavior with the single-band model due to the anomalous spin-orbit interactions in narrow band-gap semiconductors, and agree well with the experiment values reported in recent experiment [K. L. Litvinenko et al., New J. Phys. 8, 49 (2006)]. We find a strong resonant enhancement of the spin relaxation time appears for spin align along [1 (1) over bar0] at a certain electron density at 4 K. This resonant peak is smeared out with increasing the temperature.

Self-consistent analysis of double-delta-doped InAlAs/InGaAs/InP HEMTs

We have carried out a theoretical study of double-delta-doped InAlAs/InGaAs/InP high electron mobility transistor (HEMT) by means of the finite differential method. The electronic states in the quantum well of the HEMT are calculated self-consistently. Instead of boundary conditions, initial conditions are used to solve the Poisson equation. The concentration of two-dimensional electron gas (2DEG) and its distribution in the HEMT have been obtained. By changing the doping density of upper and lower impurity layers we find that the 2DEG concentration confined in the channel is greatly affected by these two doping layers. But the electrons depleted by the Schottky contact are hardly affected by the lower impurity layer. It is only related to the doping density of upper impurity layer. This means that we can deal with the doping concentrations of the two impurity layers and optimize them separately. Considering the sheet concentration and the mobility of the electrons in the channel, the optimized doping densities are found to be 5 x 10(12) and 3 x 10(12) cm(-2) for the upper and lower impurity layers, respectively, in the double-delta-doped InAlAs/InGaAs/InP HEMTs.

双δ掺杂In_(0.65)Ga_(0.35)As/In_(0.52)Al_(0.48)As赝型高迁移率晶体管材料子带电子特性研究

研究了基于InP基的In_(0.65)Ga_(0.35)As/In_(0.52)Al_(0.48)As赝型高迁移率晶体管材料中纵向磁电阻的Shubniko-de Haas (SdH)振荡效应和霍耳效应,通过对纵向磁电阻SdH振荡的快速傅里叶变换分析,获得了各子带电子的浓度,并因此求得了各子带能级相对于费米能级的位置.联立求解Schrodinger方程和Poisson方程,自洽计算了样品的导带形状、载流子浓度分布以及各子带能级和费米能级位置.理论计算和实验结果很好符合.实验和理论计算均表明,势垒层的掺杂电子几乎全部转移到了量子阱中,转移率在95%以上.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.