Stress State Determination in Nanoelectronic Silicon Devices Coupling COMSOL Multiphysics and a Recursive Dynamical CBED Pattern Simulation

Similar documents
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

In this section is given an overview of the common elasticity models.

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Numerical Heat and Mass Transfer

General displacement arch-cantilever element method for stress analysis of arch dam

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD

Computer Based Porosity Design by Multi Phase Topology Optimization

Indeterminate pin-jointed frames (trusses)

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

Constitutive Modelling of Superplastic AA-5083

One Dimensional Axial Deformations

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

One-sided finite-difference approximations suitable for use with Richardson extrapolation

CHAPTER 9 CONCLUSIONS

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Second Order Analysis

Numerical Modeling of Woven Carbon Composite Failure

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

Report on Image warping

Tensor Smooth Length for SPH Modelling of High Speed Impact

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

Finite Element Modelling of truss/cable structures

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Chapter - 2. Distribution System Power Flow Analysis

Towards the modeling of microgalvanic coupling in aluminum alloys : the choice of boundary conditions

Lecture Note 3. Eshelby s Inclusion II

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Aging model for a 40 V Nch MOS, based on an innovative approach F. Alagi, R. Stella, E. Viganò

Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM

Module 3: Element Properties Lecture 1: Natural Coordinates

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

Structure and Drive Paul A. Jensen Copyright July 20, 2003

MECHANICS OF MATERIALS

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

GENERAL THEORY FOR THE MODELLING OF LAMINATED MULTI-LAYER PLATES AND ITS APPLICATION FOR THE ANALYSIS OF SANDWICH PANELS WITH PLY DROPS

Computer Simulations of Parallel-to-Series Conversion in Solid State Frame Transfer Image Sensors. J. Bisschop

Thermo-Calc Software. Modelling Multicomponent Precipitation Kinetics with CALPHAD-Based Tools. EUROMAT2013, September 8-13, 2013 Sevilla, Spain

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

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

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

Scroll Generation with Inductorless Chua s Circuit and Wien Bridge Oscillator

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

A Hybrid Variational Iteration Method for Blasius Equation

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

LAB 4: Modulus of elasticity

Supplementary information. Coulomb oscillations in a gate-controlled few-layer graphene quantum dot

DUE: WEDS FEB 21ST 2018

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

CHAPTER II THEORETICAL BACKGROUND

Uncertainty in measurements of power and energy on power networks

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

FUZZY FINITE ELEMENT METHOD

ANALYSIS OF CONTACT PROBLEM USING IMPROVED FAST MULTIPOLE BEM WITH VARIABLE ELEMENTS LENGTH THEORY

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

Chapter 13: Multiple Regression

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

A Fast Computer Aided Design Method for Filters

Grid Generation around a Cylinder by Complex Potential Functions

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

IC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

Chapter 9: Statistical Inference and the Relationship between Two Variables

Analysis of the Magnetomotive Force of a Three-Phase Winding with Concentrated Coils and Different Symmetry Features

Monte Carlo method II

Electron-Impact Double Ionization of the H 2

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

Department of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification

Analytical and Numerical Analysis of Free Bulge Tube Hydroforming

Supplemental document

DYNAMIC BEHAVIOR OF PILE GROUP CONSIDERING SOIL-PILE-CAP INTERACTION

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

CHAPTER-5 INFORMATION MEASURE OF FUZZY MATRIX AND FUZZY BINARY RELATION

Formulas for the Determinant

Chapter 11: Simple Linear Regression and Correlation

Inductance Calculation for Conductors of Arbitrary Shape

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

Handout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION

Three-dimensional eddy current analysis by the boundary element method using vector potential

Finite Wings Steady, incompressible flow

Multidimensional parallelepiped model a new type of non-probabilistic convex model for structural uncertainty analysis

Interconnect Modeling

THE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER.

Transcription:

Excerpt from the Proceedngs of the COMSO Conference 2008 Hannover Stress State Determnaton n Nanoelectronc Slcon Devces Couplng COMSO Multphyscs and a Recursve Dynamcal CBED Pattern Smulaton A. Spessot* 1,2, S. Frabbon 1, A. Armglato 3 and R. Balbon 3 1 Numonyx, Advanced R&D, NVMTD-FTM, Va Olvett 2, 20041 Agrate Branza (Italy) 2 Natonal Research Center S3, CNR-INFM and Department of Physcs, Unversty of Modena e Reggo Emla, Va G. Camp 213/A, 41100 Modena, (Italy) 3 CNR-IMM Secton of Bologna, Va P. Gobett 101, 40129 Bologna, (Italy) *Correspondng author: alesso.spessot@numonyx.com Abstract: Straned technology s beng promoted as the best way to extend the performance of semconductor transstors. An nhomogeneous layer deposted on top of a slcon devce can nduce a strong modfcaton n the real slcon stran state, and consequently n ts electronc performance. Couplng the fnte elements analyss done by COMSO wth a recursve CBED and ACBED dynamcal smulaton, we are able to explan the observed dffracton pattern modfcaton, reconstructng the stran feld n the devce. Keywords: Straned slcon technology, devce, CBED, ACBED, TEM 1. Introducton Straned technology s beng promoted as the best way to extend the performance of semconductor transstors, and the downscalng of the slcon technology contnues the shrnkng of the devce geometry. hen the devces are straned by a non unform stressor source, the usual process devce smulators are no more able to descrbe the stran dstrbuton, and a more accurate descrpton of the resultng stran dstrbuton s requred. Ths s the case of almost all the recent mcroelectronc devces, where the actual dmensons are reduced down to 45 nm and even below for both memory cell elements and transstor devce [1]. To mprove ther performance, the straned slcon technology s currently used by the electronc ndustres [2]. The most used stressor sources are ether the lateral stran elements (.e.: SC or SGe source and dran contact) or the deposted-on-top stressor layer (.e.: tensle or compressve ntrde layer). An example of stressed devce s gven n Fg. 1, whch shows a sketch of the slcon substrate and the stressor source. The acronym CBED and ACBED refer to two TEM methods of electron dffracton and wll explaned n the followng. The correspondng electronc devce s vsble n Fg. 2 and. In ths knd of problem, a mere process devce smulaton s not Stressor sources Gate ACBED data drecton Slcon substrate CBED data drecton STI Source Slcon substrate Dran STI Fg. 1 A sketch of the analyzed sample, wth the stressor sources on top of the slcon substrate Fg. 2 A lateral and a top vew of the electronc devce; the man electronc elements are evdenced, as the geometrcal dmensons (wdth, length )

accurate enough to reproduce the stran feld, and an expermental method should be employed. 2. Expermental problems It has been recently proved that quanttatve stran mappng can be expermental obtaned by the Convergent Beam Electron Dffracton (CBED) technque of the Transmsson Electron Mcroscopy (TEM). Ths technque [3] s based on the analyss of dffracton lnes, called HOZ (Hgh Order aue Zone), whch are very senstve to the varaton n the crystal lattce parameters: a spatal resoluton n the nanometre range s feasble, [4], even n the case of hghly straned samples [5], lke the one presented n Fg. 1. Furthermore, n ths case the lnes are splt n two components, whch must be ftted to a computer smulated pattern. To perform the smulatons, a model of the nduced deformaton s needed. e set up a recursve approach, based on the refnement of some parameters nvolved n the dsplacement smulaton, untl the reconstructon of the observed sample s reached. An expermental vew of the stressng layer deposted on top of the substrate s gven by Fg. 3: four gran boundary nterfaces, generatng the dsplacement feld n the underlyng slcon, are vsble. In Fg. 4a) s presented an example of expermental pattern, wth the correspondng smulaton n Fg. 4b). An addtonal dffculty s gven by the requrement for the TEM/CBED analyss of a thnned sample, that should be transparent to the electron beam. After the thnnng, the stran state of the sample can be modfed by the elastc relaxaton wth respect to the orgnal bulk-lke case (whch s our real nterest). Therefore, to reconstruct the real sample confguraton before the thnnng, we need also addtonal nformaton from a lateral vew, coupled wth a recursve smulaton of ths case. The expermental pattern should be acqured wth a bgger feld of nvestgaton, whle a smaller resoluton s allowed. 3. Governng Equatons 3.1 Thn sample stran state The modelng of the thnned sample starts from the reproducton of the stressor source n the expermental slce. e modeled the dfferent stran state of the stressor layer onto the slcon substrate usng the Structural Mechancs Module. The deformaton nduced n the slcon substrate by the gran boundares (shown n Fg. 3) are well descrbed by a lattce msmatch [6]. e reproduced the morphology of the expermental sample observed n the TEM mage, then we used a statc analyss of a thermal relaxaton. e recursvely adjusted some selected parameters (thermal expanson coeffcents, geometrcal parameters) for the varous elements nvolved, n order to reproduce the stran state of the devces. By ths way, we obtan the nput for the Fg. 3 A strpe of gran boundary, whch are the nduced deformaton sources. The bamboo-lke dstrbuton s clearly vsble; four nterfaces between adjacent grans are shown by the arrows. Fg. 4 a) A detal of an expermental CBED pattern, ncludng splt dffracton lnes, and b) the correspondng smulaton, resultng from the recursve stran feld reconstructon.

calculaton of the recursve dynamcal CBED smulatons, n order to reproduce the expermental patterns n the 3D sample. e started assumng for the stressor source a purely radal shape (and consequently a radal deformaton). Then we ntroduced a more complex ellpsodal deformaton, adjustng the rato between the sem-axes. The resultng dsplacement feld, obtaned by COMSO, s shown n Fg. 5. 3.2 Entre sample stran state reconstructon After the reproducton of the stran state n the thnned sample, we can ntate the reconstructon of the deformaton observed n the thck case. For these case, we should use the expermental nformaton taken from a dfferent vew of the sample, acqured wth the ACBED (arge Angle CBED) technque [7]. By ths method, we can obtan nformaton about the stran state from a larger area wth respect to the CBED technque. Startng from the reconstructed deformaton model n the thn sample, we use the expermental data acqured by ACBED n the lateral vew to extend the deformaton model n the 3D space, agan usng a recursve method. 4. Theory The deformaton nduced nto the substrate by the stressor source can be computed by usng a pseudo thermal expanson, assumng a lattce msmatch between the gran boundary and the substrate. The nduced stran s computed by the followng relaton: ε = α T where α s the thermal expanson coeffcent and T s the temperature varaton. In our case, we use a value of T=1K, and then we recursvely modfed the value of the parameter α. The relaton between the resultng stran components ε j and the dsplacement s calculated accountng for the large deformaton by the Green-agrange stran equaton: ε = j 1 ( 2 R R j R Rk (1) + + ) j k where the R represents the dsplacement components along the x drectons, and the thrd term takes nto account the non-lnearty. After the stran calculaton, we compute the stress n the materal usng the consttutve relaton: σ = Dε (2) where D s the 6x6 elastcty matrx and σ represents the stress tensor. 5. Numercal Model j Ry +1.5A 200-5A nm Fg. 5 Dsplacement feld reconstructed n the thnned sample. The vertcal component s shown. Exp. Sm. Fg. 6 a) Expermental, lateral vew of the sample. Two bendng of the Bragg contour (connected to the deformaton nduced n the crstallographc plane by the stressor sources) are evdenced by the red arrows. b) Correspondng smulaton, obtaned by the ACBED smulaton after the recursve reconstructon n the thck sample.

Channel Ry 4A Stress_ 0 Channel 1500 nm Bulk -9A Fg. 7 The resultng dsplacement feld obtaned by the presented method, dsplayed n the drecton. The green arrows ndcate the stressor source on top of the slcon substrate. The current flow occurs on top of the selected regon. 500 nm -1GPa Fg. 8 The resultng dsplacement feld obtaned by the presented method, dsplayed n the drecton, representng the startng pont for the calculaton of the moblty varaton. The green arrows ndcate the stressor source on top of the slcon substrate. A peak of compressve stress up to 1GPa s reached. The dsplacement deformaton nduced nto the thn sample are calculated by a 3D geometry, usng two dfferent elements: the stressor source and the substrate. The substrate was smulated usng the standard slcon propertes, wth the assumpton of the correct ansotropy of the elastc stffness n the dfferent used drectons. The stressor source was modeled wth an ellpsodal shape, nserted on the hghest part of the substrate, as presented n the sketch of Fg. 1. Dfferent confguratons of the stran source were used, untl the observed deformaton nduced nto the substrate are ftted. e analyzed dfferent thnned samples, wth thckness n the range between 400 and 500 nm (calculated n the drecton, see Fg. 1). Then we are able to reconstruct the thck sample also n the drecton (Fg. 1). The bulk-lke was smulated wth a 3D model, usng the same thckness of the thn sample and a lateral dmenson n the mcrometer range. 6. Results and Dscusson e found that the problem of the stran feld generated by a non homogeneous stressor layer n a slcon devce can be descrbed by an approprate seres of stressors, wth dfferent thermal relaxaton parameters. Frst, we reconstruct the straned confguraton of the thn sample, by usng an approprate recursve smulaton code for the HOZ pattern smulaton [8]. The good agreement between the expermental data of Fg. 4 and the correspondng smulaton of Fg. 4 s evdent. The resultng dsplacement feld s gven n Fg. 5. Then, we can reconstruct the complete real sample stran state, usng the nformaton of the lateral vew to recursvely reconstruct the expermental ACBED pattern of the thck sample. The results of the teratve procedure s vsble n Fg. 6, whch shows the expermental data (Fg. 6a) and the correspondng smulaton (Fg. 6b). The resultng dsplacement feld, obtaned by the COMSO smulaton, s gven n Fg. 7. From these results, the stress confguraton close to the nterface can be calculated, for both the and drectons (ndcated n the sketch of Fg. 1). By the stress feld computaton n the regon of the channel, the carrer moblty can be evaluated, gvng fundamental nformaton about the electronc performance of the devces. [9]. 8. Conclusons The varous components of the stran feld n the actve area of a straned slcon devce can be obtaned, couplng the fnte elements analyss by COMSO wth a recursve CBED pattern dynamcal smulaton. Startng from ths stress state results, the carrer moblty varatons can be calculated.

8. References 1. ITRS Roadmap 2007, www.trs.net 2. S.E.Thompson, G. Sun, Y.S. Cho, T. Nshda Unaxal process nduced straned S: extendng the CMOS roadmap, IEEE Trans. El. Dev., 53,1010 (2006) 3. J.C.H. Spence, J.M. Zuo, Electron Mcrodffracton Plenum Press, New York (1992). 4. A. Armglato, R. Balbon, S. Frabbon, Improvng spatal resoluton of convergent beam electron dffracton stran mappng n slcon mcrostructures Appl. Phys. ett., 86, 63508 (2005) 5. A. Armglato, A. Spessot, R. Balbon, A. Benedett, G.P. Carnevale, S. Frabbon, G. Mastraccho and G. Pava, TEM/CBED nvestgaton of stran nduced by T-salcdes n shallow trench S solaton structures, J. of Appled Physcs, 99, 64504 (2006) 6. D. Jackob, Y. Androuss, T. Benabbas, P. Francos, A. efebvre, Surface relaxaton of straned semconductor heterostructures revealed by fnte element calculatons and transmsson electron mcroscopy, Phl Mag A, Vol. 78, 879 (1998) 7. J. P. Mornrol, arge-angle Convergent- Beam electron dffracton, SF, Pars. (2002) 8. A. Spessot, S. Frabbon, R. Balbon, A. Armglato, Method for determnaton of the dsplacement feld n patterned nanostructures by TEM/CBED analyss of splt hgh-order aue zone lne profles, J. of Mcroscopy, 223, 249 (2007) 9. P.Y. Yu and M. Cardona, Fundamentals of Semconductors, Sprnger, Berln(1996) 9. Acknowledgements One of the authors (A.S.) would lke to thank Paolo Fantn and Andrea Marmrol (of Numonyx Technology and Devce Modellng Group) for frutful dscussons.