Transport Coefficients For A GaAs Hydro dynamic Model Extracted From Inhomogeneous Monte Carlo Calculations

Similar documents
Thermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering

Event Shape Update. T. Doyle S. Hanlon I. Skillicorn. A. Everett A. Savin. Event Shapes, A. Everett, U. Wisconsin ZEUS Meeting, October 15,

Energy in Closed Systems

Multistage Median Ranked Set Sampling for Estimating the Population Median

Set of square-integrable function 2 L : function space F

19 The Born-Oppenheimer Approximation

P 365. r r r )...(1 365

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS

Generating Functions, Weighted and Non-Weighted Sums for Powers of Second-Order Recurrence Sequences

A Study about One-Dimensional Steady State. Heat Transfer in Cylindrical and. Spherical Coordinates

EE 5337 Computational Electromagnetics (CEM)

Physics Exam II Chapters 25-29

Stellar Astrophysics. dt dr. GM r. The current model for treating convection in stellar interiors is called mixing length theory:

Consequences of Long Term Transients in Large Area High Density Plasma Processing: A 3-Dimensional Computational Investigation*

Khintchine-Type Inequalities and Their Applications in Optimization

24-2: Electric Potential Energy. 24-1: What is physics

A Method of Reliability Target Setting for Electric Power Distribution Systems Using Data Envelopment Analysis

A. Thicknesses and Densities

2-DIMENSIONAL MODELING OF PULSED PLASMAS WITH AND WITHOUT SUBSTRATE BIAS USING MODERATE PARALLELISM*

Scalars and Vectors Scalar

Some Approximate Analytical Steady-State Solutions for Cylindrical Fin

MHD Oscillatory Flow in a Porous Plate

Rigid Bodies: Equivalent Systems of Forces

PHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle

NANOCAD Framework for Simulation of Quantum Effects in Nanoscale MOSFET Devices

Distinct 8-QAM+ Perfect Arrays Fanxin Zeng 1, a, Zhenyu Zhang 2,1, b, Linjie Qian 1, c

ALL QUESTIONS ARE WORTH 20 POINTS. WORK OUT FIVE PROBLEMS.

A Brief Guide to Recognizing and Coping With Failures of the Classical Regression Assumptions

8 Baire Category Theorem and Uniform Boundedness

Tian Zheng Department of Statistics Columbia University

V. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. "Flux": = da i. "Force": = -Â g a ik k = X i. Â J i X i (7.

Chapter Fifiteen. Surfaces Revisited

Bayesian Assessment of Availabilities and Unavailabilities of Multistate Monotone Systems

Thermoelastic Problem of a Long Annular Multilayered Cylinder

Large scale magnetic field generation by accelerated particles in galactic medium

(8) Gain Stage and Simple Output Stage

CFD Investigations of Spatial Arc Kinetic Influence on Fuel Burning- Out in the Tornado Combustor

UNIT10 PLANE OF REGRESSION

3.1 Electrostatic Potential Energy and Potential Difference

Test 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o?

Density Functional Theory I

Mechanics Physics 151

LASER ABLATION ICP-MS: DATA REDUCTION

Part V: Velocity and Acceleration Analysis of Mechanisms

If there are k binding constraints at x then re-label these constraints so that they are the first k constraints.

Molecular Dynamic Simulations of Nickel Nanowires at Various Temperatures

Engineering Mechanics. Force resultants, Torques, Scalar Products, Equivalent Force systems

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

PHYS Week 5. Reading Journals today from tables. WebAssign due Wed nite

The Greatest Deviation Correlation Coefficient and its Geometrical Interpretation

Impact of Polarimetric Dimensionality of Forest Parameter Estimation by Means of Polarimetric SAR interferometry

Analysis of the chemical equilibrium of combustion at constant volume

On Maneuvering Target Tracking with Online Observed Colored Glint Noise Parameter Estimation

Exact Simplification of Support Vector Solutions

Electrostatic Potential from Transmembrane Currents

E For K > 0. s s s s Fall Physical Chemistry (II) by M. Lim. singlet. triplet

Analytical and Numerical Solutions for a Rotating Annular Disk of Variable Thickness

Correspondence Analysis & Related Methods

Simulation of Surface Chemical Reactions in a Monolith Channel for Hydrogen Production

PHY126 Summer Session I, 2008

Physics Exam 3

4 Recursive Linear Predictor

A NOTE ON ELASTICITY ESTIMATION OF CENSORED DEMAND

Physics 1501 Lecture 19

Dynamic Performance, System Identification and Sensitivity Analysis of the Ladder Tracks. Ontario, Canada

3. A Review of Some Existing AW (BT, CT) Algorithms

Asymptotic Solutions of the Kinetic Boltzmann Equation and Multicomponent Non-Equilibrium Gas Dynamics

Physics 207 Lecture 16

PO with Modified Surface-normal Vectors for RCS calculation of Scatterers with Edges and Wedges

Rogerio Fernandes Brito Member, ABCM

Physics of the Earth and Planetary Interiors

Chapter 13 - Universal Gravitation

APPLICATIONS OF SEMIGENERALIZED -CLOSED SETS

Impact of donor-acceptor morphology on the charge carrier generation in organic photovoltaic devices

Rotational Kinematics. Rigid Object about a Fixed Axis Western HS AP Physics 1

Rotating Variable-Thickness Inhomogeneous Cylinders: Part II Viscoelastic Solutions and Applications

gravity r2,1 r2 r1 by m 2,1

Contact, information, consultations

Evaluation of Various Types of Wall Boundary Conditions for the Boltzmann Equation

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

9/12/2013. Microelectronics Circuit Analysis and Design. Modes of Operation. Cross Section of Integrated Circuit npn Transistor

ON THE FRESNEL SINE INTEGRAL AND THE CONVOLUTION

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Integral Vector Operations and Related Theorems Applications in Mechanics and E&M

GENERALIZATION OF AN IDENTITY INVOLVING THE GENERALIZED FIBONACCI NUMBERS AND ITS APPLICATIONS

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

CS649 Sensor Networks IP Track Lecture 3: Target/Source Localization in Sensor Networks

Solving the Dirac Equation: Using Fourier Transform

Summer Workshop on the Reaction Theory Exercise sheet 8. Classwork

Implementation in the ANSYS Finite Element Code of the Electric Vector Potential T-Ω,Ω Formulation

Slide 1. Quantum Mechanics: the Practice

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Rotating Disk Electrode -a hydrodynamic method

Machine Learning. Spectral Clustering. Lecture 23, April 14, Reading: Eric Xing 1

Groupoid and Topological Quotient Group

THE TIME-DEPENDENT CLOSE-COUPLING METHOD FOR ELECTRON-IMPACT DIFFERENTIAL IONIZATION CROSS SECTIONS FOR ATOMS AND MOLECULES

ONE-POINT CODES USING PLACES OF HIGHER DEGREE

Modeling of Electron Transport in Thin Films with Quantum and Scattering Effects

2/24/2014. The point mass. Impulse for a single collision The impulse of a force is a vector. The Center of Mass. System of particles

A Queuing Model for an Automated Workstation Receiving Jobs from an Automated Workstation

Transcription:

Tanspot Coeffcents Fo A GaAs Hydo dynamc Model Extacted Fom Inhomogeneous Monte Calo Calculatons MeKe Ieong and Tngwe Tang Depatment of Electcal and Compute Engneeng Unvesty of Massachusetts, Amhest MA 0003, USA Abstact The valdty of enegydependent elaxaton tmes often used n a GaAs hydodynamc model has been caefully examned usng the selfconsstent Monte Calo smulaton We have found that those tanspot coeffcents assocated wth the ntevalley tansfe fom the lowe to the uppe valley ae not snglevalued functons of the aveaged electon enegy n the valley If, nstead, the valley populaton ato o the aveage enegy weghted by the valley populaton ato s used a substantal mpovement n accuacy can be acheved I Intoducton Conventonally, macoscopc (aveaged) elaxaton tmes such as, p and TW appeang n the hydodynamc (HD) tanspot equatons ae assumed enegydependent and detemned by pefomng Monte Calo (MC) calculatons unde steady state and homogeneous feld condtons These expessons ae often extended to the case of nhomogeneous felds wthout any justfcaton Sandbon et al [], usng the MC smulaton, found that both the enegy and momentum elaxaton tmes n an equvalent snglevalley model unde the tansent condton dffe vey much fom the steady state and homogeneous feld values Yamada [2] also obseved the dscepancy n the elaxaton tmes between the homogeneous and nhomogeneous feld condtons He suggested that the elaxaton tmes should depend not only on the enegy but also on the valley populaton In ths study, nstead of a snglevalley model, a theevalley HD model fo GaAs has been developed The tanspot coeffcents appeang n the HD model fo homogeneous and nhomogeneos feld condtons ae evaluated by a sngle patcle MC smulaton pogam and a multpatcle selfconsstent MC smulaton pogam, espectvely A one dmensonal N + N N + GaAs ballstc dode was used as a test devce Ths appoach allows us to goously examne the valdty of the enegy dependence of each elaxaton tme as well as povdes us valuable nfomaton fo a moe appopate descpton of the elaxaton tmes II Moments of the Boltzmann Tanspot Equaton The HD tanspot equatons can be obtaned by takng vaous moments of the Boltzmann tanspot equaton (BTE) [3] Extendng the wok of [4], [5] to a multvalley system, we obtan the followng steadystate consevaton equatons fo the tth valley vm)[ + + JS], a) 'nj 'nk T"njt "^nk 65

Lv(n U l )F=^^ q A (2) n pa fjj Pk lv (^)«!# = Il^!»5 5 + a]& + =B, (3) n TW Ttyj T~Wk ^ T e j Tl T e k («) (Wj + &) F = ^ ^ *, (4) n /*«M«J M«fe whee F; = (t>),7 = ( & ), W; = fa), St = (ve) and R = (e^k) In the conventonal HD tanspot model, the tanspot coeffcents ( n, f,tw,etc) ae usually assumed to depend on the aveage enegy,w{) Ths appoach gnoes the dependence of the tanspot coeffcents on the shape of the dstbuton functon III SelfConsstent Monte Calo Smulaton To examne the accuacy of enegydependent tanspot coeffcents, we begn wth a gous soluton of the steadystate BTE by the MC method The smulaton pogam uses an analytcal multvalley, nonpaabolc band The followng types of scatteng ae taken nto account acoustc phonon scatteng, optcal phonon scatteng, pola optcal phonon scatteng, onzed mputy scatteng, equvalent and nonequvalent ntevalley scatteng The vaous scatteng paametes ae smla to those used n [6] In ths wok, a onedmensonal N + N N + GaAs stuctue wth a 05 fm Negon was examned The dopng denstes of the thee layes wee N d = X 0 7 cm' 3, x 0 6 cm" 3 and x 0 7 cm~ 3, espectvely The appled bas was 20 volts Fg dsplays the dopng densty and electc feld pofles wthn the devce as obtaned fom the selfconsstent MC (SCMC) smulaton Fg 2 shows the Y valley velocty and enegy pofles I Results and Dscussons At each poston wthn the devce we evaluated the tanspot coeffcents and the aveage enegy n each valley by the SCMC pogam We found that these coeffcents ae geneally a functon of the local aveage enegy n the valley except fo those due to the ntevalley tansfe fom the lowe to the uppe valley (e T > L, * X, and L» X) Fgs3,5 and 7 espectvely dsplay the Y * L ntevalley tansfe coeffcents, n x,, and TWTL vesus the aveage T valley enegy, W, fo both the homogeneous and the nhomogeneous feld calculaton The "hysteess" loops clealy ndcate that none of them can be descbed as a snglevalued functon of W Fo these cases, the enegydependent tanspot coeffcents whch wee obtaned fom the homogeneous feld calculaton always undeestmate the actual one n the nceasng feld egon and oveestmate the same n the deceasng feld egon If, nstead, W s weghted by the valley populaton ato (e, ^W) then the hystess loop fo ntl can be sgnfcantly educed (see Fg 4) Ths s motvated by the fact that the valley populaton ato moe o less eflects the facton of electon populaton whch have suffcent enegy to tansfe fom the lowe to the uppe valley We also found that usng the valley populaton ato alone (e, = ) the hystess a seen loop can be consdeably educed fo /pl and TWTL * n Fgs6 and 8,espectvely 66

The esult fo T C TL s vey smla to TWTL an d that fo f t L s smla to fl Once the hystess loop s educed, these tanspot coeffcents can now be modelled empcally as snglevalued functons of the valley populaton ato o the enegy weghted by the valley populaton ato Conclusons A SCMC smulaton pogam was used to examne the conventonal assumpton of enegydependent tanspot coeffcents n a multvalley system We found that the tanspot coeffcents elated to the ntevalley tansfe fom the lowe to the uppe valley (e, + L, T X and L + X) ae not a snglevalued functon of the aveage enegy n the valley A substantal mpovement n the accuacy can be acheved f the valley populaton ato o the aveage enegy weghted by the valley populaton ato s used Acknowledgement The authos acknowledge Pof D Navon fo caefully eadng the manuscpt Ths wok was suppoted n pat by NSF Gant ECS900358 Refeences [] PA Sandbon, A Rao and P A Blakey, IEEE, Tans Electon Devces, ol36, No7, pp244253, Jul, 989 [2] Yoshno Yamada, IEICE Tans, ole73, No2, pp255259, Feb, 990 [3] Blotekjag, K, IEEE Tans Electon Devces,ol ED7, Nol, pp3847, Jan, 970 [4] SC Lee and Tw Tang, Sold State Elec, ol35, P56569, Ap, 992 [5] Tw Tang,S Ramaswamy and J Nam,to appea n IEEE Tans Electon Devces [6] MA Lttlejohn, JR Hause, TH Glsson, J Appl Phys, ol48, pp45874590, Nov, 977 67

E 20 20 60 80 00 j f / / 0 02 04 06 08 2 IA 6 x(um) l 0'» Fgue Imputy dopng pofle and ts selfconsstent electc feld fo an abupt N + N N + stuctue wth an appled bas of 20 volts 0" 06 05 0A 03 02 0 ' I II 0 02 0A 06 08 2 IA A x(jun) * 0 * ST Fgue 2 Aveage electon enegy and velocty n the T valley, 0 0 H m " ';' ' ' ^ ' I 0 0 lllll lllll 0 bomot ldx u s^ 0 02 03 0A 03 ^ nthomo t^^ homo %n 2 3 (njn T )Wf(ev) l"l Fgue 3 Cae elaxaton tme, vesus the aveage T valley enegy,^ Fgue 4 Cae elaxaton tme, nl, vesus 68

00 0 'T~ " S " ' ^ "* ^ txllll '!! 4,» 'A '"x >, ^ N ">v N»»»»»* "^^ 0 nhomo fa 0 oa OJ W (ev) *" 03 0 nhomo ^fl homo fl 4 («d*) ~z Fgue 5 Cae moblty, fl, vesus the aveage T valley enegy, Wp Fgue 6 Cae moblty, f^l, vesus (j*) 0 p [ t ;" j IOOJO IOJO ' 3 jj IJO J* ^ K " ** 4 3 ~z W g IJO 0 nhon 0 02 03 W (cv) 0^ ] OA OJ 0 OJO nhomo "C^YL ~ homo t^yl 4 ' ^ ^ t, c,,,,, Fgue 8 Enegy elaxaton tme, TWTL, vesus Fgue 7 Enegy elaxaton tme, TWTL, vesus the 6J * aveage T valley enegy, Wp ^'' 69