Analysis of Halo Implanted MOSFETs

Similar documents
Using Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems

Direct Computation of Generator Internal Dynamic States from Terminal Measurements

Including the Consumer Function. 1.0 Constant demand as consumer utility function

Lecture 12: FET AC Properties

Torque OBJECTIVE INTRODUCTION APPARATUS THEORY

V DD. M 1 M 2 V i2. V o2 R 1 R 2 C C

Capacitance and Dielectrics

A simple model for the small-strain behaviour of soils

Applying Axiomatic Design Theory to the Multi-objective Optimization of Disk Brake *

EE 330 Lecture 15. Devices in Semiconductor Processes. Diodes Capacitors MOSFETs

An Improvement on Extended Impedance Method towards Efficient Steady-State Analysis of High-Frequency Class-E Resonant Inverters

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

Math 1B, lecture 8: Integration by parts

EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

SUPPLEMENTARY INFORMATION

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

1 CHAPTER 7 PROJECTILES. 7.1 No Air Resistance

DAMAGE LOCALIZATION IN OUTPUT-ONLY SYSTEMS: A FLEXIBILITY BASED APPROACH

Two Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices

Gravity wave propagation in a nonisothermal atmosphere with height varying background wind

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science

Spoof Plasmons in Corrugated Metamaterials

SYNCHRONOUS SEQUENTIAL CIRCUITS

5. Network Analysis. 5.1 Introduction

EE 435. Lecture 10: Current Mirror Op Amps

Thermal runaway during blocking

Comparison of Efflux Time between Cylindrical and Conical Tanks Through an Exit Pipe

ECEG 351 Electronics II Spring 2017

arxiv:physics/ v4 [physics.class-ph] 9 Jul 1999

Revisiting the Charge Concept in HBT/BJT Models

SiC-based Power Converters for High Temperature Applications

Lab 4: Frequency Response of CG and CD Amplifiers.

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Least-Squares Regression on Sparse Spaces

arxiv: v1 [math.co] 3 Apr 2019

EECE488: Analog CMOS Integrated Circuit Design. Introduction and Background

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009

UNDERSTANDING INTEGRATION

Section 7.1: Integration by Parts

ANALYZE In all three cases (a) (c), the reading on the scale is. w = mg = (11.0 kg) (9.8 m/s 2 ) = 108 N.

A Parametric Device Study for SiC Power Electronics

Introduction to MOSFET operation

Impact of sidewall spacer on gate leakage behavior of nano-scale MOSFETs

Where A is the plate area and d is the plate separation.

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

EE 330 Lecture 14. Devices in Semiconductor Processes. Diodes Capacitors MOSFETs

28.1 Parametric Yield Estimation Considering Leakage Variability

Chapter 6: Energy-Momentum Tensors

EECE488: Analog CMOS Integrated Circuit Design. Set 2: Background

Homework 7 Due 18 November at 6:00 pm

Analysis and Design of Analog Integrated Circuits Lecture 7. Differential Amplifiers

Lecture contents. Metal-semiconductor contact

2Algebraic ONLINE PAGE PROOFS. foundations

Implicit Differentiation

Capacity Analysis of MIMO Systems with Unknown Channel State Information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

BENFORD S LAW AND HOSMER-LEMESHOW TEST

Why Bernstein Polynomials Are Better: Fuzzy-Inspired Justification

Mathematical Analysis of Efficiencies in Hydraulic Pumps for Automatic Transmissions

Mixed-Signal IC Design Notes set 1: Quick Summary of Device Models

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue

Stopping-Set Enumerator Approximations for Finite-Length Protograph LDPC Codes

Analytical model of Threshold Voltage and Sub-threshold Slope of SOI and SON MOSFETs: A comparative study

Unified kinetic model of dopant segregation during vapor-phase growth

5 Shallow water Q-G theory.

Polytech Montpellier MEA4 M2 EEA Systèmes Microélectroniques. Analog IC Design

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Quantum mechanical approaches to the virial

J.L. Kirtley Jr. September 4, 2010

A New Conjugate Gradient Method. with Exact Line Search

Quantum Search on the Spatial Grid

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

One Dimensional Convection: Interpolation Models for CFD

This is a repository copy of The effects of gas-phase and in-depth radiation absorption on ignition and steady burning rate of PMMA.

Wire antenna model of the vertical grounding electrode

arxiv: v1 [hep-lat] 19 Nov 2013

v( t) g 2 v 0 sin θ ( ) ( ) g t ( ) = 0

Lyapunov Functions. V. J. Venkataramanan and Xiaojun Lin. Center for Wireless Systems and Applications. School of Electrical and Computer Engineering,

A Course in Machine Learning

PHYSICS BASED CHARGE AND DRAIN CURRENT MODEL FOR AlGaN/GaN HEMT DEVICES

First- and Second Order Phase Transitions in the Holstein- Hubbard Model

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Short Intro to Coordinate Transformation

EE 435. Lecture 18. Two-Stage Op Amp with LHP Zero Loop Gain - Breaking the Loop

Lecture 5 Processing microarray data

Design of Chevron Gusset Plates

Fluid Kinematics. by Dr. Nor Azlina binti Alias Faculty of Civil and Earth Resources Engineering

Linear-Time Encodable and Decodable Error-Correcting Codes MIT. ipped by the channel. codewords of the code are the individual images under

Quantum Mechanics in Three Dimensions

Equivalent rocking systems: Fundamental rocking parameters

Outline. Calculus for the Life Sciences II. Introduction. Tides Introduction. Lecture Notes Differentiation of Trigonometric Functions

A Multigrid-like Technique for Power Grid Analysis

SUPPLEMENTARY INFORMATION

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

Admin BACKPROPAGATION. Neural network. Neural network 11/3/16. Assignment 7. Assignment 8 Goals today. David Kauchak CS158 Fall 2016

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

Thermodynamic Limits of Solar Cells with Non-ideal Optical Response

Transcription:

Analysis of alo Implante MOSFETs olin McAnrew an Patrick G Drennan Freescale Semiconuctor, Tempe, AZ, olinmcanrew@freescalecom ABSTAT MOSFETs with heavily ope reions at one or both ens of the channel exhibit qualitative ifferences in electrical behavior compare to evices with laterally uniform channel opin These ifferences inclue a istinct peakiness in the transconuctance near threshol, asymmetries in capacitances, an a surprisin ecrease in the statistical variation of the peak ain factor as channel lenth ecreases istorically, accurate moelin of such evices is best one with sectional MOSFET moels ere we present an analytic moel of the behavior of the current an transconuctance of a (unilaterally or bilaterally halo implante MOSFET an show that it preicts the ecrease in variation of ain factor with channel lenth Keywors: MOSFET moelin; SPIE moelin; statistical moelin; halo implant INTODUTION ow voltae MOSFETs historically have uniform lateral opin alo (pocket implants [] were introuce to improve performance an suppress some short-channel effects, an are now stanar Specific heavily ope reions were introuce at the ens of the channel in GMOS technoloies [] to boost performance, an these coul be one selectively at source an/or rain ens, leain to so-calle unilateral an bilateral evices, havin one or both ens ope heavily, respectively These laterally non-uniformly ope evices exhibit qualitatively ifferent behavior compare to laterally uniformly ope evices First, as the ate voltae increases from below to above threshol, the lihtly ope central reion turns on at a lower ate voltae than the more heavily ope en reion(s, an therefore when the en reions start to turn on the channel lenth is effectively the lenth of these heavily ope reion(s As the ate voltae further increases an the level of inversion becomes more uniform in the evice, the channel lenth becomes more the lenth of both the heavily an lihtly ope reions The (D electrical effective channel lenth eff is therefore shorter near threshol than for well above threshol, an this ives a istinct peake m s b behavior near threshol (we assume 0 here A very interestin analysis of halo implante MOSFETs was presente in [3] This work provie a physical threshol voltae moel for bilaterally ope evices, an unequivocally showe that a 3-section MOS moel is neee to capture fine etails of halo implante MOS behavior; merely moelin threshol voltae for a -section moel is insufficient Althouh the peake m behavior is not explicitly note in [3] it is apparent in Fis 0 an in [3] But a moel explainin this behavior an its statistical variation was not provie source s 0 N T ate N Fi MOSFET compose of 3 separately ope reions For a unilaterally ope evice (assumin the heavy opin is on the source sie of the evice, as the evice turns on the rain en inverts before the source en, hence there is a istinct peak in near threshol, an the evice is hihly asymmetric an s, even at 0 Aain, this behavior cannot be capture by a -section MOSFET moel, especially with the emphasis of the past ecae in improvin the symmetry of MOSFET moels For MOSFETs with constant lateral opin, the statistical variation in the peak transconuctance m µ 0oxW chanes with channel lenth; for short evices it is reater than for lon evices ariation in the oxie thickness T ox an the low fiel mobility µ 0 contribute equally to the variations in m for both lon an short evices, but variations in eff affect a short evice more than a lon evice, so the overall variation in m increases as lenth ecreases This behavior is not observe in evices with nonuniform lateral opin In unilateral GMOS evices the measure variation in m is less for short evices than for lon evices We present an analysis that explains this counterintuitive behavior The basic moel unerlyin the analysis is simple, an is not intene to provie a hihly accurate escription of all aspects of evice behavior But it is amenable to analysis that eluciates the unerlyin cause of the unexpecte ecrease in m as channel lenth ecreases Note that DMOS transistors also have laterally nonuniform opin, but this is because their boy reion is forme by out-iffusion an not because of a separate heavy opin reion at the source or rain The analysis here is not irectly applicable to DMOS evices, but oes capture some of the qualitative behavior they exhibit, because a unilaterally ope GMOS evice is qualitatively similar to an DMOS evice NUMEIA BIATEA AO MODE onsier a MOSFET semente into 3 reions, see Fi Ajacent to the source is a reion of lenth an threshol voltae T, in the mile is a reion of lenth an N T rain 594 NSTI-Nanotech 007, wwwnstior, ISBN 4006844 ol 3, 007

threshol voltae, an ajacent to the rain is a reion of lenth an threshol voltae T The source an bulk are assume to be at zero olts an voltaes an are applie to the ate an rain For operation in stron inversion non-saturation, at small, the currents (normalize to µ 0 ox W in each of the 3 sections, which must be equal, are approximately T 05 ( i + B 05 ( i T + B 05( (3 i ( where an are the voltaes between the left (source an center sections, an between the center an riht (rain sections, respectively, an B an B are the boy effect factors for the central an riht reion, respectively The terms involvin the last two parameters are require because the effective backate bias for the center an riht sections (takin the source as the reference is not zero an in ( throuh (3 vary with ate bias Iteratively solvin for these, as a function of ate bias an overall ate lenth + +, then computin the current an m (an normalizin with respect to the m at just above threshol, ives the results shown in Fi one to specifically eluciate the effects of the non-uniform lateral opin an to not have them confoune with other physical phenomena As observe in experimental ata, the m exhibits a peake characteristic (aain, this is not from series resistance or vertical fiel mobility eraation, which increase the peakiness, the latter uniformly over lenth an the former more for short than for lon evices; it is purely from the intrinsic nature of the laterally non-uniform opin for a laterally uniformly ope evice m woul be inepenent of bias in Fi Note that the peakiness is non-monotonic in lenth Asymptotically as 0 the evice appears, for N N, to be laterally uniformly ope, an as the relative effect of the en reions becomes vanishinly small, therefore the m m, max curve shoul be flat for both of these cases Fi 3 shows the peakiness (the ratio of the peak m to its value for lare as a function of (reciprocal channel lenth Fi 3 m peakiness vs reciprocal channel lenth Startin from short evices (rihtmost abscissa values in Fi 3, as the channel lenth increases (movin to the left, the peakiness increases, as is observe in practice For very lare lenths the peakiness aain ecreases, as is expecte from the qualitative analysis above Fi Normalize m vs for various lenths For the simulations use to enerate Fi, 008µ m (the simulations here were base on the technoloy of [], which is a 04µm MOS technoloy with a minimum eff of 05µm from the GMOS structure, T T 04 an 0 Note that for the simulations series resistance is inore an mobility reuction ue to both vertical an lateral fiels is not inclue, so surface rouhness scatterin an velocity saturation are inore; this is 3 ANAYTIA BIATEA AO MODE Instea of numerically solvin ( throuh (3, analytically ifferentiatin with respect to, an roppin neliible terms,, an are small, ives (4 (5 + T m, + m, NSTI-Nanotech 007, wwwnstior, ISBN 4006844 ol 3, 007 595

(6 T m, (etaile calculations that keep the omitte terms for now an eliminate them later are extremely teious an lea to the same results iven below Because the current throuh each section must be the same, the chane in each current with must the same, so equatin (4 an (6 ives (7 Equatin (5 an (6 an usin (7 ives (8 T T + T + T T + T T Substitutin (8 into (4 an rearranin, ives (aain normalize to µ 0 ox W (9 T + T( T + ( T( m ( T + T( T + T( Note the inherent symmetry an consistency of (9; settin any of reion lenth to zero (notin that if 0 then 0, if 0 then, an if 0 then reuces the expression to an obviously similar form; as oes settin the parameters of one reion to those of an ajacent reion Assumin a symmetric halo implant, an T T, then for small, equatin ( an (, an notin that at low for a symmetric evice, (0 + T This inicates that moelin the evice as the series connection of separate reion channel resistances, as in [3], is reasonable owever, our explicit solution for the internal noe voltaes verifies this approach alculation of the transconuctance ives ( m + ( + This is simpler than the form (9 yet almost as accurate owever, it is still not intuitive how this varies with eometry an statistical fluctuations Aain for left an riht en reions with the same physical characteristics (lenth, opin, threshol voltae, for a lon evice from (9 m (3 0,, + Because >, this approximation preicts that m increases as increases, which is qualitatively ifferent from the behavior in Fi This is because it oes not take into account the better approximation (0 A tihter upper boun is (4 m + For a short evice (5 m 0 5, 05, + (As 0, because the central reion has a lower opin an hence lower threshol than the en reions, it is turne on more than the ens, an has a lower resistance; thus approaches (which is clearly the expecte result for The current in the transistor is then, normalize to µ 0 ox W, ( I + Fi 4 Normalize m from analytic approximations 596 NSTI-Nanotech 007, wwwnstior, ISBN 4006844 ol 3, 007

Fi 4 compares the simple approximations (4 an (5 with values calculate from numerical simulations base on iterative solution of ( throuh (3 (results from ( are not inclue; they are almost ientical to those from the numerical simulations; results from (5 are also not inclue for the 0µm lon evice as this approximation is only applicable to short evices The lon channel simple approximation (4 is clearly somewhat inaccurate owever, the areement of the very simple analytic moel (5 with the more etaile numerical moel from the iterative solution of ( throuh (3 is quite reasonable 4 ANAYSIS Because the halo reion is more heavily ope than the center of the channel >, an the ratio of the halo to central reion ate overrives in the enominator of (5 is less than The effect of this ratio on the overall evice transconuctance ecreases as ecreases, hence the smaller the the shorter the effective channel lenth in the enominator of (5, an the hiher the m This is the root cause of the peakiness of the m characteristic, an this simple an intuitive theoretical analysis arees both with the qualitative escription in the Introuction an with experimental observation that m is more peake for short than for lon evices More interestinly, consier what happens when,,, or in (5 vary Increasin any of these parameters woul normally, as a first response, be expecte to lea to a reuction in current, an hence m This is seen to be the case for the parameters,, an, but is unexpectely not true if increases Increasin causes to ecrease (at a fixe ate voltae, which reuces the overall enominator in (5 an therefore increases m An alternative viewpoint is that if increases compare to, the central reion is turne on harer when the halo reions invert, this makes the effective channel lenth appear shorter, an hence enhances m Numerical simulations confirm this counter-intuitive behavior Statistical variations in the ain factors come from statistical variations in T ox an µ 0, which (inorin the uncorrelate, or mismatch, component are the same for short an lon evices, an from statistical variations in,,, an For illustrative an comparison purposes we inore the effect of the variations in T ox an µ 0 an consier only the variations in the lenth an threshol parameters Propaation of variance (Po analysis [4] of (5 therefore ives (6 δm m e + 4 δ + δ + where + e δ δ Applyin this Po analysis to the 5µm an 04µm channel lenth evices, for which Fi 4 shows that the approximation (5 is reasonable, ives the results shown in Fi 5 an Fi 6 The contribution of the variations in each of the parameters, at the - level, is shown in these fiures The counter-intuitive conclusion from this analysis arees with experimental observations: the statistical variation in the peak m for halo evices is reater for loner evices than for shorter evices Fi 5 Statistical variations in m for the 5µm evice Fi 6 Statistical variations in m for the 04µm evice Most strikinly, it is the ecrease of the influence of variation in as lenth ecreases on the ain factor that stans out Analytically, the 4 th term in the riht-han sie of (6 shoul ecrease as ecreases, an this is observe in numerical simulations an ata owever, overall contributions in practice are more ifficult to compartmentalize The empirical conclusion from Fi 5 an Fi 6 is that statistical variations in peak m ecrease as channel lenth ecreases, for lenths in a certain rane This conclusion is supporte by experimental ata, which show that the peak m variation for GMOS type evices ecreases as ecreases NSTI-Nanotech 007, wwwnstior, ISBN 4006844 ol 3, 007 597

5 UNIATEA AO ANAYSIS For a evice with an implant only at the source en of the channel, which in some ways is analoous to an DMOS evice, (7 T + T m ( + T an for lon evices 0 an for short evices, hence (3 an (5 still apply, with substitute for ; qualitatively the behavior is ientical The analyses here have been tarete to bilateral halo implante evices, but the oriinal analysis we evelope was base on unilateral evices, to try to unerstan physically why the observe statistical variation in the peak m was smaller for shorter evices than for moerate lenth evices 6 ONUSIONS We have provie an analysis of the rain current of halo implante MOSFETs, at low s, base on equality of the currents flowin in ifferent reions of the evice This enable expressions for the ain factor m s to be erive The analysis eluciates the reason for the qualitative chane in the shape of m s plots, an also for the observe ecrease in statistical variation as channel lenth ecreases Our analyses support the conclusion of [3], that there are qualitative features of the characteristics of halo implante evices that cannot be moele usin a sinle MOSFET moel; a multi-section moel (two or three sections, epenin on the structure of the evice is neee owever, rather than analyzin a evice from the channel resistance perspective, as note in [3], we foun it avantaeous to consier equality of currents in ifferent reions of the evice This le to the physical unerstanin of why the peak m variation for short evices is unexpectely less than for moerate lenth evices, an what physical parameter variations unerlie this behavior EFEENES [] Y Okumura et al, A Novel Source-to-Drain Nonuniformly Dope hannel (NUD MOSFET for ih urrent Drivability an Threshol oltae ontrollability, IEEE IEDM Tech Di, pp 39-394, Dec 990 [] F K hai et al, A ost-effective 05µm eff BiMOS Technoloy Featurin Grae-hannel MOS (GMOS an a Quasi-Self-Aline (QSA NPN for F Wireless Applications, Proc IEEE BTM, pp 0-3, Sep 000 [3] ios et al, A Three-Transistor Threshol oltae Moel for alo Processes, IEEE IEDM Tech Di, pp 3-6, Dec 00 [4] McAnrew an P G Drennan, Unifie Statistical Moelin for ircuit Simulation, Proc MSM-WM, pp 75-78, 00 598 NSTI-Nanotech 007, wwwnstior, ISBN 4006844 ol 3, 007