Search NASA⌕ Search

SEARCH · Search NASA

Results for “drift-diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Improved millimeter-wave mixer performance analysis at cryogenic temperatures

The results of a 600-GHz mixer performance analysis using an improved model for computing the Schottky diode capacitance-voltage (C-V) relationship are presented. The computed C-V data for a realistic submillimeter-waver mixer diode are given as a function of physical temperature and compared to the standard analytic expression based on a solution of Poisson's equation. Both C-V relationships are used to predict the performance of an ideally terminated 600 GHz mixer operating at 300, 140, 80, and 30 K. It is shown that the drift-diffusion model more accurately describes the mixer performance when the physical temperature is reduced below roughly 100 K.

Siegel, Peter H.↗

Thermionic emission current in a single barrier varactor

From I-V measurements on Single Barrier Varactors (SBV) at different temperatures we concluded that thermionic emission across the barrier of the actual device is mainly due to transport through the X band. The same structure was also modeled with a one-dimensional drift-diffusion model, including a 'boundary condition' for thermionic emission across the heterojunction interface. By including thermionic field emission through the top of the triangular barrier of a biased diode and the effect of a non-abrupt interface at the heterojunction, we obtained good agreement between the modeled and measured I-V characteristics.

Hjelmgren, Hans↗

A Parallel Multilevel Spectral Element Scheme

A parallel multilevel strategy is developed using spectral (p) finite elements. Hierarchic bases are particularly well suited since the element matrices and vectors are nested and the multilevel projections easiliy performed. Since the basis degree is used to specify the multigrid level, an EBE strategy is natural br the multilevel technique. Results are presented for two candidate nonlinear elliptic transport problems: the augmented drift-diffusion equations of semiconductor device modeling and the stream function-vorticity equations of incompressible fluid dynamics.

Davis, M. B.↗

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗

Comparison of Non-Parabolic Hydrodynamic Simulations for Semiconductor Devices

Parabolic drift-diffusion simulators are common engineering level design tools for semiconductor devices. Hydrodynamic simulators, based on the parabolic band approximation, are becoming more prevalent as device dimensions shrink and energy transport effects begin to dominate device characteristic. However, band structure effects present in state-of-the-art devices necessitate relaxing the parabolic band approximation. This paper presents simulations of ballistic diodes, a benchmark device, of Si and GaAs using two different non-parabolic hydrodynamic formulations. The first formulation uses the Kane dispersion relationship in the derivation of the conservation equations. The second model uses a power law dispersion relation {(hk)(exp 2)/2m = xW(exp Y)}. Current-voltage relations show that for the ballistic diodes considered. the non-parabolic formulations predict less current than the parabolic case. Explanations of this will be provided by examination of velocity and energy profiles. At low bias, the simulations based on the Kane formulation predict greater current flow than the power law formulation. As the bias is increased this trend changes and the power law predicts greater current than the Kane formulation. It will be shown that the non-parabolicity and energy range of the hydrodynamic model based on the Kane dispersion relation are limited due to the binomial approximation which was utilized in the derivation.

Smith, A. W.↗

Theoretical Study of Response of InGaAs Metal-Semiconductor-Metal Photodetectors

We present a theoretical investigation of the response of metal-semiconductor-metal (MSM) photodetectors made of InGaAs lattice-matched to InP using a two-dimensional drift-diffusion model with a thermionic-field emission boundary condition for the heterojunctions. The effect of including a top InAlAs layer to increase the effective barrier height of the metal fingers on the InGaAs active layer is thoroughly examined and found to limit the collection of the photocurrent signal due to the electron and hole barriers that it forms with InGaAs. Due to the thickness and height of the InAlAs barrier layer in existing designs, the tunneling current obtained from the model is found to be negligibly small to significantly affect the output signal current. In an attempt to obtain a better response, different design structures including one where a quasi-Schottky contact is utilized are studied and their speed of response, breakdown voltage, and dark current are compared to that of the usual InGaAs device.

Salem, Ali F.↗

Efficient Multi-Dimensional Simulation of Quantum Confinement Effects in Advanced MOS Devices

We investigate the density-gradient (DG) transport model for efficient multi-dimensional simulation of quantum confinement effects in advanced MOS devices. The formulation of the DG model is described as a quantum correction to the classical drift-diffusion model. Quantum confinement effects are shown to be significant in sub-100nm MOSFETs. In thin-oxide MOS capacitors, quantum effects may reduce gate capacitance by 25% or more. As a result, the inclusion or quantum effects in simulations dramatically improves the match between C-V simulations and measurements for oxide thickness down to 2 nm. Significant quantum corrections also occur in the I-V characteristics of short-channel (30 to 100 nm) n-MOSFETs, with current drive reduced by up to 70%. This effect is shown to result from reduced inversion charge due to quantum confinement of electrons in the channel. Also, subthreshold slope is degraded by 15 to 20 mV/decade with the inclusion of quantum effects via the density-gradient model, and short channel effects (in particular, drain-induced barrier lowering) are noticeably increased.

Biegel, Bryan A.↗

Hierarchical Approach to 'Atomistic' 3-D MOSFET Simulation

We present a hierarchical approach to the 'atomistic' simulation of aggressively scaled sub-0.1 micron MOSFET's. These devices are so small that their characteristics depend on the precise location of dopant atoms within them, not just on their average density. A full-scale three-dimensional drift-diffusion atomistic simulation approach is first described and used to verify more economical, but restricted, options. To reduce processor time and memory requirements at high drain voltage, we have developed a self-consistent option based on a solution of the current continuity equation restricted to a thin slab of the channel. This is coupled to the solution of the Poisson equation in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison to the full self-consistent solution. At low drain voltage, a single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. In this case, the current is calculated by solving the current continuity equation in a drift approximation only, also in a thin slab containing the MOSFET channel. The regions of applicability for the different components of this hierarchical approach are illustrated in example simulations covering the random dopant-induced threshold voltage fluctuations, threshold voltage lowering, threshold voltage asymmetry, and drain current fluctuations.

Asenov, Asen↗

Electro-Thermal Simulation Studies of SiC Junction Diodes Containing Screw Dislocations Under High Reverse-Bias Operation

The objective of this work was to conduct a modeling study of SiC P-N junction diodes operating under high reverse biased conditions. Analytical models and numerical simulation capabilities were to be developed for self-consistent electro-thermal analysis of the diode current-voltage (I-V) characteristics. Data from GRC indicate that screw dislocations are unavoidable in large area SiC devices, and lead to changes in the SiC diode electrical response characteristics under high field conditions. For example, device instability and failures linked to internal current filamentation have been observed. The physical origin of these processes is not well understood, and quantitative projections of the electrical behavior under high field and temperature conditions are lacking. Thermal calculations for SiC devices have not been reported in the literature either. So estimates or projections of peak device temperatures and power limitations do not exist. This numerical study and simulation analysis was aimed at resolving some of the above issues. The following tasks were successfully accomplished: (1) Development of physically based models using one- and two-dimensional drift-diffusion theory for the transport behavior and I-V characteristics; (2) One- and two-dimensional heat flow to account for internal device heating. This led to calculations of the internal temperature profiles, which in turn, were used to update the electrical transport parameters for a self-consistent analysis. The temperature profiles and the peak values were thus obtainable for a given device operating condition; (3) Inclusion of traps assumed to model the presence of internal screw dislocations running along the longitudinal direction; (4) Predictions of the operating characteristics with and without heating as a function of applied bias with and without traps. Both one and two-dimensional cases were implemented; (5) Assessment of device stability based on the operating characteristics. The presence of internal non-uniformities, particularly filamentary structures, was probed and demonstrated; (6) Cause and physical origins of filamentary behavior and unstable I-V characteristics were made transparent; (7) It was demonstrated that diodes containing defects would be more prone to thermal breakdown associated with the temperature dependent decrease in the thermal conductivity; and (8) Finally, negative differential resistance (S-shaped NDR) which can potential lead to device instability and filamentary behavior was shown to occur for diodes containing a line of defects such as could be associated with a screw dislocation line.

Joshi, R. P.↗

A Self-Consistent Plasma-Sheath Model for the Inductively Coupled Plasma Reactor

Accurate determination of ion flux on a wafer requires a self-consistent, multidimensional modeling of plasma reactor that adequately resolves the sheath region adjoining the wafer. This level of modeling is difficult to achieve since non-collisional sheath lengths are usually 3-4 orders of magnitude smaller than the reactor scale. Also, the drift-diffusion equations used for ion transport becomes invalid in the sheath since the ion frictional force is no longer in equilibrium with drift and diffusion forces. The alternative is to use a full momentum equation for each ionic species. In this work we will present results from a self-consistent reactor scale-sheath scale model for 2D inductively coupled plasmas. The goal of this study is to improve the modeling capabilities and assess the importance of additional physics in determining important reactor performance features, such as the ion flux uniformity, coil frequency and configuration effects, etc. Effect of numerical dissipation on the solution quality will also be discussed.

Bose, Deepak↗

Quantum Mechanical Study of Nanoscale MOSFET

The steady state characteristics of MOSFETS that are of practical Interest are the drive current, off-current, dope of drain current versus drain voltage, and threshold voltage. In this section, we show that quantum mechanical simulations yield significantly different results from drift-diffusion based methods. These differences arise because of the following quantum mechanical features: (I) polysilicon gate depletion in a manner opposite to the classical case (II) dependence of the resonant levels in the channel on the gate voltage, (III) tunneling of charge across the gate oxide and from source to drain, (IV) quasi-ballistic flow of electrons. Conclusions dI/dV versus V does not increase in a manner commensurate with the increase in number of subbands. - The increase in dI/dV with bias is much smaller then the increase in the number of subbands - a consequence of bragg reflection. Our calculations show an increase in transmission with length of contact, as seen in experiments. It is desirable for molecular electronics applications to have a small contact area, yet large coupling. In this case, the circumferential dependence of the nanotube wave function dictates: - Transmission in armchair tubes saturates around unity - Transmission in zigzag tubes saturates at two.

Svizhenko, Alexei↗

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗

Simulation of Ultra-Small Electronic Devices: The Classical-Quantum Transition Region

Concern is increasing about how quantum effects will impact electronic device operation as down-scaling continues along the SIA Roadmap through 2010. This document describes part of a new semiconductor device modeling (SDM) program at NAS to investigate these concerns by utilizing advanced NAS and third-party numerical computation software to rapidly implement and investigate electronic device models including quantum effects. This SDM project will investigate quantum effects in devices in the classical-quantum transition region, especially sub-0.1 mm MOSFETs. Specific tasks planned for this project include the use of quantum corrections to the classical drift-diffusion and hydrodynamic models of electron transport, arid the use of nominally quantum models including significant scattering.

Biegel, Bryan A.↗

Multi-Dimensional Quantum Tunneling and Transport Using the Density-Gradient Model

We show that quantum effects are likely to significantly degrade the performance of MOSFETs (metal oxide semiconductor field effect transistor) as these devices are scaled below 100 nm channel length and 2 nm oxide thickness over the next decade. A general and computationally efficient electronic device model including quantum effects would allow us to monitor and mitigate these effects. Full quantum models are too expensive in multi-dimensions. Using a general but efficient PDE solver called PROPHET, we implemented the density-gradient (DG) quantum correction to the industry-dominant classical drift-diffusion (DD) model. The DG model efficiently includes quantum carrier profile smoothing and tunneling in multi-dimensions and for any electronic device structure. We show that the DG model reduces DD model error from as much as 50% down to a few percent in comparison to thin oxide MOS capacitance measurements. We also show the first DG simulations of gate oxide tunneling and transverse current flow in ultra-scaled MOSFETs. The advantages of rapid model implementation using the PDE solver approach will be demonstrated, as well as the applicability of the DG model to any electronic device structure.

Biegel, Bryan A.↗

Efficient Multi-Dimensional Simulation of Quantum Confinement Effects in Advanced MOS Devices

We investigate the density-gradient (DG) transport model for efficient multi-dimensional simulation of quantum confinement effects in advanced MOS devices. The formulation of the DG model is described as a quantum correction ot the classical drift-diffusion model. Quantum confinement effects are shown to be significant in sub-100nm MOSFETs. In thin-oxide MOS capacitors, quantum effects may reduce gate capacitance by 25% or more. As a result, the inclusion of quantum effects may reduce gate capacitance by 25% or more. As a result, the inclusion of quantum effects in simulations dramatically improves the match between C-V simulations and measurements for oxide thickness down to 2 nm. Significant quantum corrections also occur in the I-V characteristics of short-channel (30 to 100 nm) n-MOSFETs, with current drive reduced by up to 70%. This effect is shown to result from reduced inversion charge due to quantum confinement of electrons in the channel. Also, subthreshold slope is degraded by 15 to 20 mV/decade with the inclusion of quantum effects via the density-gradient model, and short channel effects (in particular, drain-induced barrier lowering) are noticeably increased.

Biegel, Bryan A.↗

Hypersonic Magneto-Fluid-Dynamic Compression in Cylindrical Inlet

Hypersonic magneto-fluid-dynamic interaction has been successfully performed as a virtual leading-edge strake and a virtual cowl of a cylindrical inlet. In a side-by-side experimental and computational study, the magnitude of the induced compression was found to be depended on configuration and electrode placement. To better understand the interacting phenomenon the present investigation is focused on a direct current discharge at the leading edge of a cylindrical inlet for which validating experimental data is available. The present computational result is obtained by solving the magneto-fluid-dynamics equations at the low magnetic Reynolds number limit and using a nonequilibrium weakly ionized gas model based on the drift-diffusion theory. The numerical simulation provides a detailed description of the intriguing physics. After validation with experimental measurements, the computed results further quantify the effectiveness of a magnet-fluid-dynamic compression for a hypersonic cylindrical inlet. At a minuscule power input to a direct current surface discharge of 8.14 watts per square centimeter of electrode area produces an additional compression of 6.7 percent for a constant cross-section cylindrical inlet.

Shang, Joseph S.↗

A theoretical analysis of steady-state charge collection in simple diodes under high-injection conditions

A previous rigorous mathematical analysis of drift-diffusion equations was used to investigate collected charge in a simple reverse-biased p-n junction diode exposed to an ionization source that liberates carriers (electron-hole pairs) in a quasi-neutral region within the diode. Each of two simple models was found to agree with the more rigorous analysis when carrier liberation is sufficiently intense. One is the sensitive volume (SV) model, and the other was called “ambipolar diffusion with a cutoff” (ADC). The earlier rigorous analysis was worked out in detail only for a localized source, i.e., a point source of carrier liberation, so it was able to validate the applicability of each simple model only for that case. The present paper treats an arbitrary spatial distribution of carrier generation and concludes that the ADC model remains valid for this more general case, but the SV model does

Edmonds, Larry D.↗

Effect of Plasma Sheaths on Earth Re-entry MHD Processes

The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. Results obtained indicate that the cathode sheath can exert significant resistance to current flow and thus reduce the MHD-generated force by 40% or more. The effect of the sheath on the MHD force is found to be more pronounced at lower applied magnetic field.

Earth Re-entry↗