EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

Similar documents
Math 333 Problem Set #2 Solution 14 February 2003

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Vehicle Arrival Models : Headway

CHAPTER 12 DIRECT CURRENT CIRCUITS

20. Applications of the Genetic-Drift Model

04. Kinetics of a second order reaction

Comparison between the Discrete and Continuous Time Models

ψ(t) = V x (0)V x (t)

NEWTON S SECOND LAW OF MOTION

Physics for Scientists & Engineers 2

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

72 Calculus and Structures

KINEMATICS IN ONE DIMENSION

Final Spring 2007

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Chapter 7 Response of First-order RL and RC Circuits

8. Basic RL and RC Circuits

Supplementary material

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

10. State Space Methods

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

AP Chemistry--Chapter 12: Chemical Kinetics

IV. Transport Phenomena. Lecture 19: Transient Diffusion. 1 Response to a Current Step. MIT Student. 1.1 Sand s Time

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

INDEX. Transient analysis 1 Initial Conditions 1

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

MOMENTUM CONSERVATION LAW

Phys1112: DC and RC circuits

Solutions to Assignment 1

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

Predator - Prey Model Trajectories and the nonlinear conservation law

5.1 - Logarithms and Their Properties

Chapter 10 INDUCTANCE Recommended Problems:

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

Sub Module 2.6. Measurement of transient temperature

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

Lecture 13 RC/RL Circuits, Time Dependent Op Amp Circuits

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

= ( ) ) or a system of differential equations with continuous parametrization (T = R

8.022 (E&M) Lecture 16

4.5 Constant Acceleration

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Homework 2 Solutions

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

Reliability of Technical Systems

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Linear Response Theory: The connection between QFT and experiments

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Chapter 13 Homework Answers

Chapter 8 The Complete Response of RL and RC Circuits

Solutionbank Edexcel AS and A Level Modular Mathematics

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

Physics 1402: Lecture 22 Today s Agenda

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

Chapter 2. First Order Scalar Equations

ENV 6015 Solution to Mixing Problem Set

EE100 Lab 3 Experiment Guide: RC Circuits

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Advanced Organic Chemistry

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

) were both constant and we brought them from under the integral.

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

( ) is the stretch factor, and x the

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

Section 2.2 Charge and Current 2.6 b) The current direction is designated as the direction of the movement of positive charges.

13.3 Term structure models

Math 10B: Mock Mid II. April 13, 2016

ln y t 2 t c where c is an arbitrary real constant

Effects of Coordinate Curvature on Integration

Some Basic Information about M-S-D Systems

Diffusion & Viscosity: Navier-Stokes Equation

2.4 Cuk converter example

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

IB Physics Kinematics Worksheet

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Hall effect. Formulae :- 1) Hall coefficient RH = cm / Coulumb. 2) Magnetic induction BY 2

6. Stochastic calculus with jump processes

4.1 - Logarithms and Their Properties

Class Meeting # 10: Introduction to the Wave Equation

CSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering

Lecture 14. The term Brownian motion derives its name from the botanist Robert Brown whom, in 1828, made careful

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections

Direct Current Circuits. February 19, 2014 Physics for Scientists & Engineers 2, Chapter 26 1

MEI Mechanics 1 General motion. Section 1: Using calculus

Differential Equations

Transcription:

EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he source and he drain erminals of a p-mos are grounded and a negaive gae bias (w. r.. source and drain) is applied. I is aribued o he breaking of Si- bonds a he Si-SiO inerface and hereby creaing inerfacial raps (). ere, he applied ve gae bias acs as a sress. When he above menioned bias configuraion is mainained, he densiy of ( ) follows a emporal variaion as shown qualiaively in Fig.. I is essenially an power-law variaion. One needs o have a quaniaive perspecive of such variaion o esimae and improve he device lifeime. When he sress (he gae bias) is removed he bonds sar `healing and he rap densiy reduces. The overall variaion is shown qualiaively in Fig.. ln Threshold for proper circui operaion n eed o be on his curve o increase lifeime Lifeime of circui Fig. : The variaion of rap densiy, ln wih ime (solid line). G V G S D V G 3 n A + ξ + Fig.. Trap densiy variaion under sress and wihou sress.

9. Reacion-Diffusion Model for BTI: In his lecure we will deduce he variaion of rap densiy ( ) under sress. The dynamics of he variaion can be expressed by he reacion diffusion equaions: d () k f ( () ) kr ( ) ( x, ) () d d d ( x, ) d ( ) ( ) ( x, ) D + x, µ E x δ () d dx d where, is he oal Si- bond densiy a he inerface, is he densiy of free ydrogen, D and µ are he diffusion consan and he mobiliy, respecively, of he same. E is he elecric field and δ is he widh of he small region near he inerface where he bond-breaking reacion akes place. In all he analysis in his lecure, x denoes he Si-SiO inerface and posiive x-direcion direcs owards he oxide as shown in Fig 3. ( x) Region where he reacion akes place ydrogen generaion due o bond breaking Free ydrogen curren due o drif and diffusion x x δ x Si SiO Fig. 3: Coninuiy of free hydrogen The firs erm on he righ hand side of equaion () denoes he forward reacion (i.e. breaking of bonds o increase he rap densiy) has been se proporional o he densiy of he unbroken bond. The second erm denoes he reverse reacion (i.e. resoraion of he bond) has been se proporional o he rap densiy (broken bond densiy) and he available ydrogen a he inerface. Equaion () is essenially he coninuiy equaion for he free ydrogen a he inerface. The firs erm on he lef hand side of his eq. () denoes he generaion rae of free ydrogen due o he bond-breaking and he erm inside he parenhesis denoes he curren of ydrogen a he inerface (Fig. 3). The difference of hese wo erms is proporional o he rae of increase of ydrogen concenraion a he inerface, which is zero a seady sae.

Anoher convenien form of eq. () is is inegral form: f () () ) dx where, he upper limi of he inegraion is a funcion of ime and will paramerically depend on D and µ for he case of diffusion and drif respecively. This equaion merely saes ha he oal amoun of free ydrogen inside he oxide is equal o he rap densiy. The exac funcional form of f () will become clearer as we proceed furher ino he discussion. We will now consider he soluion of he reacion diffusion equaion under differen condiion. 9.3 Soluion of Reacion-Diffusion Equaion in he Sress Phase 9.3. Righ afer he sress has been applied: We can neglec in comparison o in eq. () and also neglec he second erm on he righ hand side of eq. () as much of he bonds remain inac and no much free ydrogen is produced. Wih hese approximaions we ge from eq. (): () k f (4) which is an exponenial variaion wih an exponen of. 9.3. Afer some ime has elapsed beyond he applicaion of sress: ow, is no ha small o be discarded in he way as we did before as considerable amoun of bonds, by now, are broken and here is significan amoun of free ydrogen a he inerface. I is sill small compared o as in he previous case. Bu, unlike he previous case we canno neglec he reverse reacion anymore. Moreover, he ne rae of increase of rap densiy becomes very small as he reverse reacion sars. We ge, wih all hese assumpions, from eq. (): k f () ( x, ) (5) k r We will now consider hree cases: (i) ydrogen flow due o diffusion only: Diffusion is he dominan process when he free ydrogen is no ionized. Assuming ha here is no any kind of recombinaion of ydrogen inside he oxide we ge a linear profile of free ydrogen densiy inside he oxide as shown in Fig. 4. A any ime he profile exends o a disance of D. This is he funcional form of f () in eq. (3). Subsiuing his value of f () in eq. () and considering a linear profile as shown in Fig. 4 we ge: D Dividing eq. (6) by eq. (5) and hen solving for () we ge: (3) ) dx D ( x, ) (6)

/ /4 /4 ( ) kf D A kr which is again an exponenial variaion wih a fracional exponen of /4. (7) ) ( x, ) D x Fig. 4: Free ydrogen densiy profile under diffusion (ii) ydrogen flow due o drif only: Drif becomes he dominan process when he ydrogen is ionized (proon). In his case he ydrogen densiy profile looks fla all he way upo he edge of he prifle a x µ E (Fig. 5). So, he funcional form of f () is µ E. Subsiuing i in eq. (3) and considering a fla profile as shown in Fig. 5 we ge: µ E ) dx µ E Dividing eq. (8) by eq. (5) and solving for () we ge: ( x, ) (8) k µ E k r which is an exponenial variaion wih a fracional exponen of /. f () / (9) ) ( x, ) µ E Fig. 5: Free ydrogen densiy profile under drif x

(iii) eural diffusion (wih +! ) Since he molecules are charge neural diffusion is he dominan process. So we ge: D ) dx D (, ) () x The reason for muliplying he inegral by is ha here are wo ydrogen aoms in a x x, ydrogen molecule. We also need a relaion beween (, ) and ( ) k { ( x ) } k ) (), Finally, from eqns. (5), () and () we ge: / 6 ~ : The summary of cases (A), (B) and (C) is shown in Fig. 6 where we see ha he slope of he variaion becomes fracional afer iniial phase. The exac value of he slope depends on he underlying mechanism of ydrogen flow inside oxide and he mechanism of reacion. () ln n Drif of + n 4 n 6 Diffusion of Diffusion of n Righ afer he applicaion of sress ln Fig. 6: Exponenial variaion of rap densiy wih ime as per he reacion-diffusion model 9.3.3 Sauraion of rap densiy a he end of reacion: The variaion of rap densiy wih ime a he end of he reacion, when all he bonds are broken, becomes sauraed. From eq. (), we can wrie for neural ydrogen: d () d ( x, ) ( x, ) D D (3) d dx D

A he end of he reacion when approaches we canno neglec rap inside he firs erm of righ hand side of eq. (). We are sill neglecing he ne rae of increase of rap densiy. Wih all hese consideraions we ge from eq. (): k f ( ) (, ) (4) () x kr Subsiuing eq. (4) ino eq. (3) and solving for we ge: + / ln k (5) which can be approximaed by (easily verified by second order expansion in Taylor series) 4 η e (6) which shows ha as,, i.e. he rap densiy smoohly sauraes o he value of Si- bond densiy. The oal variaion of rap densiy is shown in Fig. 7. ln Drif of + Diffusion of Sauraion Diffusion of Righ afer he applicaion of sress Fig. 7: Toal variaion of rap densiy as per he reacion-diffusion model ln 9.4 Conclusions: Today we discussed one feaure of he emporal degradaion of BTI, namely he origin of power-exponen in he sress phase. We will consider he oher wo aspecs of emporal degradaion relaxaion and frequency-independence in he nex class.