Processes Using Advanced Physical and


 Angelina Greene
 1 years ago
 Views:
Transcription
1 156 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. ED32, NO. 2, FEBRUARY 1985 Simultion of C:llritil IC Fbrition Proesses Using Advned Physil nd Numeril Methods ERNER JUNGLING, PETER PICHLER, S : EGFRIED SELBERHERR, SENIOR MEMBER, IEEE, EDGAR GUERRERO, AND HANS. POTZL, MEMBER, IEEE Abstrt Critil steps of IC fbritionresimultedby ont  nd the eletrostti potentil, et. In the first step, the pplitwodimensionl omputer progrms using dvned physil modek. Our tion is minly sientifi nd the ode is used to study odes delwithnrbitrrynumberofphysilquntitiessuh! on dvned models. e investigte the vlidity rnge of entrtions of dopnts, vnies, interstitils nd lusters, the eletr:)s,tti potentil, nd so on. Furthermore, they esily permit the exhnge or urrent models nd estimte, s losely s possible, if vrition of the physilmodelsunderonsidertion. As typill:plli improvements of more omplited models re justified by tions phenomen of oupled diffusion in one nd two dimensions nd the dditionl mount of omputtion time nd memory dynmi rseni lustering re investigted. The differenes used Iy the resoures. The progrm hs been developed to quikly vry models of the ero spehrge pproximtionndthesolution (1 the physil models (e.g., vrying the number nd kind of ext Poisson eqution restudied by exmples of AsB diffusioiwith vrious doping onentrtions t different tempertures. A dynmi l~ster physil quntities, modifying some of the physil prmmodel developed for the simultion of thermlly nneled As implnt uions eters, et). Sine the min purpose of the progrm is is ompred to mesureddt of lser nneling experiments. A r,,hort outline of the mthemtil nd the numeril problems is given to!8bhow the mount of sophistition neessry for uptodte proess simul ion. INTRODUCTION I. HE RECENT dvnes in devie miniturition lve T lled for better understnding of physil proxmes nd their effet upon the dopnt distribution during dtwie nd iruit fbrition, espeilly in VLSI. To dte, fir mount of work hs been done in terms of experiment,: 1 ion ' nd in the development of dvned physil models /:.g., dynmi luster/preipittion models, enhnement of: Idiffusion by vnies, interstitils nd the eletril fkld, et.). This reserh requires effiient omputer progrnu; to evlute, to ompre, to verify, or to rejet the vrilms models. In ddition, urte simultions of proess!itrps re often desirble to ompute unknown prmeters IC] to lulte more preisely roughly estimted prmeters of physil models. The omplited struture of dvned physil mdels often exeeds the pbilities of ommonly used simultion progrms. The implementtion of models often req~il~res simplifitions whih diminish the expeted improver1:nts by the models. To ountert this sitution we hve developed 1;I:nerlpurpose simultion progrm whih trets n rbit :ilry number of physil quntities suhs onentrtion; of donors, eptors, vnies, interstitils nd lusters, lnd Mnusript reeived Otober 18, This work ws sponsored blr the Siemens Reserh Lbortories, Munih, Germny, nd the Fond ur Forderung der wissenshftlihen Forshung, Projet S 43/1. The uthors re with the Institut fur Allgemeine und Elektr,:~~ik, Abteilung fur Physiklishe Elektronik, 14ein, Austri. onentrted on the physil problem, the mesh neessry for the disretition is reted utomtilly nd dpted with respet to the speifi problem. This feture frees the user from mthemtil nd numeril problems. In Setion I1 we explin the pbilities of our ode nd disuss the differentil equtions whih n be solved by our progrms. In Setion I11 we use our ode to study some phenomen of oupled diffusion nd ompre the results of the qusineutrl pproximtion to those using the ext Poisson eqution. Setion IV ompres dynmi luster models for shllow highly doped rseni implnttions. A short outline of the numeril methods nd pbilities is given in Setion V. 11. SPECIFICATIONS OF THE PROGRAMS On speifying the pbilities of our ode, we hd to find ompromise between wide rnge of pplitions nd n eptble mount of omputtion time nd memory requirement. The speifitions to be desribed in the following re onerned with ion implnttion nd redistribution of dopnts by diffusion. Ion implnttion is one of the most pplied doping tehniques in IC fbrition, espeilly in VLSI. The simultion of implnttion must be bsi pbility of ny proess simultion progrm. Advned simultion progrms tend to use distribution funtions whih re speified by higher number of moments to desribe the impurity onentrtions, e.g., Person IV distributions. Sine our ode is lid out to hve minly sientifi pplitions, we hve implemented ll ommonly used distribution funtions. This llows for investigtions of differenes /85/2~156$ IEEE used by the ssumption of different impurity distribution
2 JUNGLING et i. : SIMULATION OF IC FABRICATION PROCESSES 157 funtions nd for omprison with ny published simultions, In the first step, we hve instlled the simple Gussin, the joined hlf Gussin, nd the Person IV distributions. Ion implnttion n be simulted in single lyer s well s in multilyer strutures. The lyers my onsist of rbitrry ordering of silion, silion nitride, nd silion dioxide. A summry of the distribution funtions nd their mthemtil bkground n be found in [ll]. Sometimes it is neessry to use experimentl dt, obtined by SIMS, n(), or spreding resistne mesurements s n initil solution for further proess simultion. It should be noted tht mesured dt often need mthemtil pretretment to yield suffiiently ontinuous distribution funtions for the simultion of following diffusion step. An exmple my be short time nneling step using C or pulsed lser. Monte Crlo methods hve been omitted t present beuse they led to very time onsuming omputtions. This level of implementtion permits the omputtion of nerly ll impurity distributions whih my be the initil solution for following diffusion step. Impurity redistribution during IC fbrition will be the min pplition of the first version of our ode desribed in this pper. The physil model on whih ll onsidertions re bsed is genexl ontinuity eqution (1) nd generl flux reltion (2). The number n of physil quntities in (1) nd (2) is rbitrry nd is only restrited by omputer resoures. The ordering nd the physil mening of the quntities must be speified by the user through few prmeter definitions. This feture offers the possibility to optimie the progrm for speifi problem (e.g., omitting unneessry quntities or modifying some physil prmeters). ith respet to fixed mount of memory, the user n hoose between dditionl physil quntities or dditionl grid points to obtin higher ury of the simultion. If the physil quntity is onentrtion, (1) nd (2) stte tht the onentrtion n be modified by diffusion flux, field flux, or genertion or reombintion proess. The first two effets re inluded in Fik s lw nd re ommonly used to desribe dopnt redistribution wheres the other effets beome essentil when time dependent trnsformtions between physil quntities our (e.g., vniesinterstitils, eletrillytiversenilustered rseni, et.). Certinly ny of these effets n be inluded or exluded if neessitted by the orresponding physil model. It n be esily seen tht (1) nd (2) lso inlude the Poisson eqution. The progrm solves system of n prtil differentil equtions (PDE) of the form n, n i]. =  pij. divj, + G, (1) j=1 J, =  Di.grdCi +/ p,.ci.grdt,l (2) where ij, pij, Di, pin nd Gi re funtions of spe, time nd ny of the quntities Ck nd + represents the eletrostti pohtil. The index i runs from 1 to n nd indites the ith eqution of the differentil eqution system. It is ler tht unusul funtions, whih would hnge the type of the differentil equtions, must be exluded. This onstrint, however, does not ffet the ommonly used modeling of the prmeters. This set of PDE s inludes the possibility to simulte simple diffusion step, Poisson eqution, oupled diffusion of dopnts withsmeor different hrge, dynmi lustering of dopnts, the effet of preipittion nd enhned diffusion by vnies or interstitils. The vribles Ck in (1) nd (2) represent not only onentrtion or potentil, but n lso represent the rdius; of preipittes or ny other physil quntity. By hnging the funtions, p, D, p, nd G we n simply vry physil models nd pply the min prt of our ode unhnged to very mny problems. The initil solution for the PDE s n be omputed from the modeling of ion implnttion, n be given by the bkground impurities of the wfer or n be omputed from mesured dt. It is worthwhile now to disuss the boundry onditions in more detil. In the first level of our ode we onstrin ourselves to rigid boundries, i.e., we exlude the simultion of oxidtions with moving boundries. The boundry onditions for eh of the differentil equtions re given by (3). j=l n j=l tij*c, n + qi+*j, + B, (3) where tij, qij re funtions of spe nd time nd B, is funtion of spe, time, nd ny of the quntities k. It is worth noting tht tij nd qij re independent from ij nd Pij This formultion inludes Dirihlet, Neumnn, nd mixed boundry ondition. From the physil point of view, these implementtions offer the possibility to simu lte n inert diffusion step (eg, Ji = ), predeposition step (C, = C,), n idel sink for vnies or interstitils (C, = C;q), onentrtion dependent soure of interstitils Jr = J,(k), et COUPLEDIFFUSION Field enhnement strongly ffets the migrtion of dopnts in diffusion steps. The effet domintes in the viinity of pn juntions, espeilly when the rrier onentrtions ner the juntion exeed the ilntrinsi number. Negleting hevy doping effets, the ext formultion of the field enhnement is given by (4)(6).  As t,  div ( DAs. grd CAS + pas Cis grd IJ ) (4)  div(d,.grdc,p,.c,.grdrc/) (5) t This orret implementtion is, for the ske of simpliity,
3 158 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. ED32, NO. 2, FEBRUARY 1985 seldomly used in simultion progrms. hen only eqully hrged dopnts re present the influene of the int:rnl field is modeled by the fieldenhnement ftor o' the diffusion oeffiient (7). h=l+&.($+l) .5 If pn juntions re present, the nlyti solution :' the qusineutrl pproximtion, (3), reples the Poi won eqution. Both methods redue the mount of omputtion tinll: t lest by ftor two sine the mount of omputtion time is proportionl to bout the squre of the number o ' the physil quntities. The use of the fieldenhnement, Ttor, however, dditionlly simplifies the struture 1 the ontinuity eqution, using n dditionl nonliner t8l:rm in the diffusivities. Considering simultion of rs :ni, boron, nd the eletrostti potentil we n redus: the mount of omputtion time nd memory by ftor 4/9 using (8) insted of (6). These fts hve supported the implementtion of the qusineutrl pproximtion into simultion progrms for engineering pplition where high performne nd smll mount of memory re powxful rguments. Some publitions, e.g., [1][3], hve delt with oupled diffusion nd the question if the qusinelltrl pproximtion my reple the ext Poisson equtio 3, or not. The results presented in these ppers re not sl.i;fying. In [3] the ext Poisson eqution ws repled b! the depletion ssumption, e.g., set n =, p =, nd t) ornputed by (6) in the viinity of the juntion. In the other domins the qusineutrl pproximtion hs been used to reple the Poisson eqution. The boundry between the two domins wsposed t the end of the speh rge region. In [2] n ext solution of the Poisson equtim is presented; however, the Computtions do not involve pn juntion. Corret implementtion of fieldenhned diffusiw is not the only point of interest. Sometimes it is limed :ht the influene of the eletrostti potentil on the hitlrge stte of vnies neessittes preise knowledge of tlme physil quntities or t lest n estimtion of the differenes used by vrious models. As it is often prop rled tht the diffusion of, e.g., boron, rseni, et. ours m.illy by hrged vnies, the eletrostti potentil n influene the enhnement of diffusion vi point defets, Sine reently developed models for diffusion hve heome more dvned nd urte, omprison betnjelen the results of the two models seems neessry. As the SI: e hrge uses the differene between the two models, 'we will without loss of generlity onentrte our invest!:tions to boron nd rseni pn juntions with firly ljgh onentrtions t the juntions I I I. I I I I I, I, I S.6.7 DEPTH [pml Fig. 1. Arseni nd boron profilefter ion implnttions. As: 4~1'~ m2 with 9 kev, LSS, Gussin; B: 4X 1'5 m2 with 8 kev, LSS, Gussin. The eletrostti potentil is omputed from the ext Poisson eqution t 95 C. The first exmple demonstrtes worst se estimtion. Ion implnttions, As 4X lo1' m2 with 9keV nd boron 4X 115 m2 with 8 kev, led to boron dopnt onentrtions whihexeed those ommonly used. Furthermore, we neglet ny lustering effets of boron nd rseni to obtin n extreme eletril field nd onsiderble spe hrges. The diffusivities re defined by (9) nd (1). They re modeled temperture dependent nd inlude enhned diffusion by hrged point defets. D, = exp m2 (  3;; ev ) S exp m2 S ( 3.46; qt) m2 (; ; 3 e~ ) DAs =.66  exp s exp S ( kt Fig. 1 shows the initil distribution of boron nd rseni nd the solution of the potentil due to (6) for temperture of 95 C. Fig. 2 revels the finl distribution nd the eletrostti potentil fter n inert diffusion of 3 min t 95OC. The trnsient behvior of the dopnts is worth being investigted more refully. Figs. 35 revel survey of the time dependent migrtion of the boron nd rseni dopnts in the viinity of the juntion. They indite tht the diffusion onsists of two periods, short fieldontrolled period (. 1 s) followed by n mbipolr diffusion. The strong fieldontrolled migrtion of both dopnts in the viinity of the juntion leds to steepening of the dopnt profiles. This uses the trnsition from the n domin to the pdomin to beome nrrower nd inreses, therefore, the eletril field. These two effets
4 JUNGLING et i. : SIMULATION OF IC FABRICATION PROCESSES 159.b.4 u 1" v) Ly U i IAs 1 i i b.7 DEPTH [pl > I J. 5 I Fig. 2. Arseni nd boron distribution nd eletrostti potentil fter n inert diffusion step of 3 min t 95 C. The initil distribution is given in Fig. 1. t ll.12 DEPTH [pml.1 m e Fig. 4. Arseni nd boron distributions fter 3s inert diffusion t 95OOC due to the models P nd Q Irr E I y r, u y m 4 u 3 5 VI y U 3! U y u s lb : m.1 r7 e 1. 5 I 1 t DEPTH [pml Fig. 3. Arseni nd boron distributions fter 1s inert diffusion t 95OoC due to the models P nd Q..1.9.ll.12 DEPTH [pl Fig. 5. Arseni nd boron distributions fter 1s inert diffusion t 95 C due tothe models P nd Q. enhne eh other nd led to stedily inresing flux of dopnts, espeilly of boron toms whih hve higher mo'bility. The proess will weken s soon s the dopnt migrtion willhveredued the brupt trnsient in the juntion nd will hve repled it by smooth juntion. The deresing eletril field wekens the fieldontrolled flux. Figs. 35 show the As nd B onentrtions, nd 4 for both models t 1, 3, nd loos, respetively. The potentil, due to the models, differs from the beginning of the diffusion step. The steeperslope of 4 omputed by model Q (qusineutrl pproximtion (4),(5),(8))leds to stronger eletril field. Therefore the fieldontrolled migrtion of the boron nd rseni due to model Q is more distint thin tht of model P (ext Poisson eqution (4)(6)). The boron profiles revel the d.ifferenes more lerly beuse the mobility is lrger nd the grdient is smller thn tht of the rseni profiles. Fig. 6 revels the trnsient dependene of the spe hrge nd onfirms our findings. Sine the spe hrge, using the min differene between the two models Q nd P, disppers fter 1 s, the differenes between the two models will ertinly not inrese. On the ontrry, Figs. 35 indite tht the smll differenes whih hve been built up in the first period re too smll to resist the following diffusion step. Fig. 7 supports the suggestion. Fig. 7 shows the time dependene of the quntities t depth of.99 pm. This depth oinides with the nside of the pn juntions lose to the juntion. Differenes of the dopnt onentrtions
5 16 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. ED32, NO. 2, FEBRUARY 1985 Fig. 6. Time dependent survey of the positive spe hrge dt;nsity whih is defined s the mximum of nd n p + N,,  ND 111 the viinity of the pn juntion for n inert diffusion step t 95 C. lnitil nd finl distribution re given by Figs. 1 nd 2, respetively. 9 2 Fig. 7. 1s 1 12 lo3 TIME [sel 1.3, Trnsient behvior of rseni, boron nd the eletrostti pt~l entil t depth of.99 pm for model P nd model Q. begin to build up fter 1 s, remin stble until 2 s, l.nd then slowly dispper. The dopnt migrtion during the following time peiod is minly hrteried by n mbipolr diffusion in the viinity of the pn juntion. This mens tht the dopmts migrte in suh wy tht neither signifint spe hrge nor signifint eletri field n be built up. The mbipolr diffusion seems ontrolled by the diffusivity of Ithe rseni toms whih leds to strongly retrded diffu:;~~on of the boron toms. Awy from the juntion the diffu4on of the mjority dopnts revels norml results wheres i.he migrtion of minority dopnts is strongly retrded. Th :i is used by the extreme vlues of the potentil whih nhnes mjority diffusivities but retrds the minority di.fpusivities. A omprison of the initil profile (Fig. 1) nd the finl profile (Fig. 2) revels the impression tht the tbwo dopnt distributions moved towrds eh other, Fkt whih is not unexpeted, sine the dopnts hve diffel nt hrge. It is worth noting tht lthough the models 1e.I to different dopnt onentrtions in the viinity of the jl ]ICtion, the juntion depth is not ffeted by these differen :es. It might be interesting if the smll differenes t the pn juntions whih our t the beginning of the diffusion!i tep re suffiiently lrge to use different eletril ltiirteristi of the devie. Our seond exmple shows more relisti proess. Two implnttions, As 4X lo1 m with 9keV nd boron 7 ~ 1 m ~ with ~ 4 kev, re followed by n 1 C nneling step for 4 min. In this exmple the boron onentrtion is fr too smll to influene the eletrostti potentil t the ndoped side. Sine the pek vlue of the boron onentrtion is of the sme order of mgnitude s the intrinsi number (6X1Ol8 m3 t 1 C), we n lmost neglet the influene of the boron onentrtion on the omputtion of the potentil wy from the pek vlue. The influene of the field n lerly be seen by vlley in the boron profile (Fig. 8) used by the field ontrolled flux. The eletril field pushes the boron toms towrds the surfe ginst the diffusion flux. Certinly the vlley in the boron profile moves s the steep grdient of the rseni profile shifts into the bulk. The spe hrge, Fig. 9, indites gin two periods. A short time domin hrteried by remrkble spe hrge nd eletril field nd following period dominted by norml diffusion of the rseni nd the retrded diffusion of the boron toms. In ontrry to tlie former exmple we nnot spek of mbipolr diffusion. The rseni migrtion is hrdly influened by the boron onentrtion sine rseni exeeds the intrinsi number nd the boron onentrtion by some orders of mgnitude. On the other hnd, the boron migrtion is ontrolled by the rseni in two wys. Firstly, the movement of the boron toms isfield ontrolled in the viinity of the pn juntion nd seondly strongly retrded in the regions of high rseni onentrtion. (The positive eletrostti potentil suppresses the seond term in the diffusion oeffiient of (9).) Here gin, the qusineutrl pproximtion omputes steeper slope of the eletrostti potentil, using stronger eletril field. Therefore, the vlley^' in the boron profile omputed by model Q is more distint thn tht of model P. The mximum differene between the two models n be expeted when the spe hrge rehes its mximum, whih ours t 3 s. Fig. 1 shows the time dependene t depth of.12 pm. The differenes used by the two models pper s soon s the juntion moves into the viinity of the depth under onsidertion nd disppers s soon s the juntion moves wy. Fig. 11 shows the pn juntion t time of 3 s. The stronger eletril field uses deeper vlley in the boron profile for model Q thn for model P. Awy from the pn juntion no differenes n be seen. The depth of the pn juntion is not dependent on the model beuse the steep grdient of the rseni ompenstes the different boron profiles. The two exmples justify the use of the qusineutrl pproximtion insted of the ext Poisson eqution when the durtion of the diffusion step is long enough to exeed the first, fieldontrolled period of dopnt migrtion. The first exmple revels tht the differenes diminish fter t lest 2 s t temperture of 95 C. e ould not observe ny signifint differenes for typil diffusion steps in IC fbrition t tempertures beyond 95 C. At lower tempertures one hs to be more reful. To estimte the temperture dependene of the oupled diffusion we hve simulted the dopnt migrtion of exmple one t
6 JUNGLING et /. : SIMULATION OF IC FABRICATION PROCESSES t B i ' Fig. 8. Timede endent boron rofile for n inert diffusion step of 4 min t 1 8 The initil hribution hs been obtined by ion implnttions. As: 4X lo1' m', 9 kev nd B: 7x1l3 m2, 4 kev. Fig. 9. Trnsient behvior of the positive spe hrge density whih is defined s the mximum of nd np+ N,  N, for the proess of Fig n E I D! m n U V  v) D! U '* 1 ' Depth .12pm TInE [sel.1 3 _I 4 5. ' lo2 '.l Fig. 1. Time dependene of rseni, boron nd the eletrostti potentil due to the models P nd Q for n inert diffusion t 1oOO"C t depth of.12 pm. temperture of 8 C. The diffusivities t 8" nd 95 C differ typilly by ftor of 1 exp( 3$ ") = { 6.92 X for T = 8 C 6.64~ for T = 95OOC. G? o m 4 VI D! equt ion As. B.3  P _I.2 4.I. .l.1.15 DEPTH [pml 1' ' '  Fig. 11. Arseni, boron, nd the eletrostti potentil due to the models P nd Q for n inert diffusion t 1 C: fter 3 s. For diffusionontrolled proess we would expet the 8 C proess to be retrded by ftor 1 ompred.to the 95 C proess. But, beuse the first period is fieldontrolled we hve to modify our retrdtion ftor. The reltion between the field flux nd the diffusion flux n _be hosen to estimte the time dependene in the.fi'r$ period. If we onfine our interest to the viiiiity of the'pn juntion nd ssume tht CAS = CB, DAs << D, d CAS: CB << 2. ni, some lgebr leds to JF: JD = ni. Sine ni evlutes to 4.4X 1l8 m3 nd d m3 t 95" nd 8"C, respetlvely, we, tnerefore, expet ll diffusion proesses to be retrded by ftor 1 wheres the field ontrolled proesses re only retrded by ftor of 1/3 = 33. Furthermore, the effets of fieldontrolled dopnt migrtion should be more distint. Our simultions onfirm our rough estimtions. For 8 C the spe hrge shows its mximum t 17 s in ontrst to the 5 s of Fig. 6 for the 95 C step. The dopnt distnitions in the viinity of the juntion t 3, 1, nd 3 s for the 8 C proess re similr to the distributions t 1, 3, nd 1 s for the 95 C step. The qulittive behvior of the two proessd differs little. The lst exmple onfirms tht even t low tempertures, the qusineutrl pproximtion desribes the dopnt migrtion qulittively well in the viinity of pn juntion nd exellently in the other domins. The :smll differenes t the juntion re probbly not suffiiently lrge to use onsiderble effets on the eletril behvior of the devie. As simple exmple of our twodimensionl ode we hve simulted oupled rseni boron diffusion. Equtions (4) nd (5) re the ontinuity equtions for boron nd rseni, respetively, nd (6) represents Poisson's eqution. The diffusion oeffiients re modeled s given by (7) nd (8). e simulte n inert nneling step t 1.O"C temperture. As initil stte we ssume homogeneously doped boron, substrte mp3) implnted with rseni mp2, 13 kev, LSS, Gussin) through msk (35 nm)
7 162 IEEE TRANSACTIONS ON ELECTRONDEVICES, VOL. ED32, NO. 2, FEBRUARY Fig. 12. Twodimensionl distribution of the eletrostti potentil id'ter n inert diffusion of 1 s t 1OOO"C. Fig. 14. Twodimensionl boron distribution fter n inert diffusion of 1 s t 1OOO"C. 3% $< t x i? 8% LI. T? Fig. 15. Twodimensionl boron distribution fter n inert diffusion of 5 s t 1OOO"C. Fig. 13. Twodimensionl distribution of the boron dopnts fte' n inert diffusion of 1 s of 1OOO"C. with n infinitely steep edge to field oxide (t tt, origixt of the oordinte xis in the following figures). The influene of the eletril field on the boron prlile is the min purpose of this simultion. Sine the rseni onentrtion exeeds the intrinsi number nd the bmn onentrtion by severl orders of mgnitude we do riot expet notible effets of the eletril field on the slmpe of the rseni profile. The diffusion of the rseni profil,: is, therefore, of minor interest nd we onentrte on 1,he disussion of the eletrostti potentil nd the boron profile. Fig. 12 shows the eletrostti potentil. t the beginrmg of the diffusion. The field term in the ontinuity equ:l~on for boron uses flux of doping toms in the diretion of positive grdients of the potentil. Fig. 13 indites tht in the beginning of the diffusion the boron profile forn:u; Fig. 16. Twodimensionl distribution of the eletrostti potentil fter n inert diffusion of 15 s t 1OOO"C. The potentil profile spreds with inresing diffusion time orresponding to the spreding rseni profile s indited in Figs. 12 nd 16. The mximum of the boron profile lose to the surfe is, therefore, redued (Fig. 14) nd vnishes finlly (Figs. 15 nd 17), wheres the mximum in the bulk inreses nd spreds. It should be mentioned tht the minimum of the boron profile remins t the surfe next to the msk edge (origin of the oordi Ushped mximum of inresed onentrtions t Ithe surfe nd in the bulk. The mximum ner to the surfe is smller thn the mximum in the bulk. This is used by the less ];~'I:onouned grdient of the eletrostti potentil t the nte system) used by the lterl field ontrolled flux of surfe. Sine the mbient ondition t the surfe (Jl;,ron boron. Affeted by the field in vertil nd lterl dire = ) inhibits flux of boron perpendiulr to the surl'ile tion, the boron profileforms distint pekvluewith the onentrtion of boron is reltively strongly deplete13 in inresing diffusion time. This extreme umultion of this region. At the edge of the Ushped mximum the boron profile is depleted. The minim of the boron prdile re used by the grdients of the potentil in lterl md vertil diretion t the origin of the oordinte system. boron toms n only be explined by twodimensionl simultion. The shpe of the potentil uses lterl nd vertil fluxes of boron whih umulte the dopnts in enter nd use the pek vlues. After 15 s diffusion
8 JUNGLING et l. : SIMULATION OF IC FABRICATION PROCESSES ' i Fig. 17. Twodimensionl distribution of the boron onentrtion fter n inert diffusion of 15 s t 1 C. time no signifint hnge in the qulittive behvior of diffusion will tke ple. An experimentl verifition of this firly simple exmple is, most unfortuntely, not possible, sine the resolution of even the most modern, e.g., SIMS equipment is too orse for tht purpose. However, we hve ertinly verified the numeril ury of our progrm with problems whose solutions re known DEPTH [pml Fig. 18. Mesured nd omputed dt of the totl nd eletrilly tive rseni onentrtion fter 2min nneling step t 1 C. IV. ARSENICLUSTER MODELS Arseni impurities in silion hve temperture dependent solubility, e.g., [4] or [5]. Arseni onentrtions below this solubility limit re eletrilly tive. Beyond this limit they form lusters nd/or preipitte nd beome eletrilly intive. Therefore, the eletrilly mesured profiles (e.g., spreding resistne mesurements) differ from tomisti mesurements (e.g., SIMS mesurements). The usge of shllow high dose As implnttions for the fbrition of MOS soure nd drin regions in VLSI devies neessittes dynmi luster models for proess simultion. A fir mount of work hs been done to evlute equilibrium solubilities but reltively few experiments hve been done reently to determine the trnsient behvior of lustering nd delustering. e hve investigted the experiments published in [SI nd [6]. In [5] mesurements nd simultions of typil MOS soure fbrition proesses re simulted: An implnttion of 2 X 1l6 ern?, rseni with n energy of 14 kev through 25nm oxide followed by n nneling step of 2 min t 1 C (+ further nneling step t 8 C for 6 min). The physil model used for the simultion is well desribed in [5], ((11)(14) nd (6)) llowing us to proeed with the simultion nd test the bilities of our ode. The only differene from t,he model of [5] is the use of the ext Poisson eqution insted of the qusineutrl pproximtion. The exhnge of the rseni diffusivity from (1) to (14) offers no problems nd uses only miniml modifitions of the ode.  'As t m.as+ + k. e + Cl('"k)+ (11)  div( DAs. grd CAS + pas* CAS. grd $J) + k,.c,, k,.c,",.nk (12) Fig. 19. Trnsient dependene of the eletrilly tive rseni 2 min nneling step t 1 C. m.= l t k,ec, + k,.c,",.nk during ,4.1 ev l+loo.n/ni DAs = 22.9 exp ( kt )*mr m2.s'. (14) In (13) nd (14) k, nd k hve to be orreted by ftor m to be onsistent with the formuls in [5]. The only prmeter to be mthedhsbeen the intrinsi number whih equls 7.2 X lls m3 t 1 C in our lultion. Fig. 18 shows the result fter the first nneling step whih oinides well with the mesured dt nd the published results in [5]. Figs. 19 nd 2show the trnsient dependene of the tive nd the lustered rseni. Clustering nd delustering re very fst proesses ompred to diffusion t 1 C. Fig. 19 revels tht the equilibrium between the tive nd the lustered rseni hs well been obtined fter 2 min. The luster onentrtion shows signifint mximum
9 164 IEEE TR4NSAJCTIONS ON ELECTRON DEVICES, VOL. ED32, NO. 2, FEBRUARY 1985 Fig. 2. Trnsient dependene of the lustered rseni onentrtion during 2min nneling step t 1 C. 1 OZ2 Fig. 22. Trnsient dependene of the eletrilly tive rseni 6min nneling step t 8 C. during lo1 ;1 119 DI t u g lol U 116. u o etetr. tive As totl A5, DEPTH Lpml Fig. 21. Totl nd eletrilly tive rseni onentrtions ftelr I:O nneling steps. fter 3 s of nneling nd n esily be explined. In Ithe first 3 s, lustering is the dominting proess nd re;ults in inresing luster onentrtions. Then the diffusioll of the tive rseni leds to deresing tive onentrt~ons whih re immeditely ompensted by delustering. :%is leds to derese of the luster onentrtion. The flt profiles in the tive rseni onentrlion results from the ft tht the hnge of CAS with tin].: is proportionl to the devitions from the equilibrium (onentrtion. The simultion of the following 8 C nneling r:lep for 6 min is shown in Figs. 21 nd 22. The results ointide gin well with the experimentl dt nd the res Jlts obtined by [5]. The desending slope of C,, in Fig. 22 indites tht the equilibrium between lustered nd tive rseni hs not yet been hieved. The use of stti insted of dynlni luster model would hve resulted in too smll live rseni onentrtion. One differene between the ode presented in [5] nd lour,ode is worth being disussed more refully. In [5] i is stted tht lustering effets re omputed individu~il y nd onentrtions re updted with smll time intervls. This indites tht (12)(14) re evluted in two steps ( diffusion step followed by lustering step). To the ontrry, our ode solves (12) to (14) simultneously. Clerly the differenes between the two methods beome negligible if the time nd spe disretition issuffiientlysmll. Nevertheless, we believe our method to be the better one from physil nd mthemtil point of view. The differenes will our when the grid width in spe nd time is inresed, in order to obtin fster nd less memory intensive progrm. e lsownt to pronoune tht dynmi phse trnsformtion, suh s the rseni lustering, hs not to be implemented into the ode but represents typil physil effet for whih our ode hs been developed. The model of Tsi is bsed on experiments using only therml tretment of the wfer to investigte the lustering of rseni. It seems interesting to ompre this model to experiments whih hve been published in [6]. These experiments use lser nneling fter the implnttion to be sure tht ll rseni rests on lttie sites nd is eletrilly tive. Further therml tretment by C C,lser nd onventionl furne is performed to investigte the trnsient behvior of the lustering. The experimentl results re shown in Fig. 23. The xxis shows the nneling time in logrithmi sle, the yxis the sheet rrier onentrtion per m2 surfe. e hve hosen the impurity distribution fter the lser nneling s the initil solution for our lultion nd hve simulted the nneling experiments tking into ount diffusion s well s lustering. The dshed lines show the simulted results using the model of [5] for 9, 8, nd 7 C. Unfortuntely mesurements nd simultion differ gretly. This my be used by the different nneling methods pplied to the wfers fter implnttions nd/or by the different therml post tretment. The different totl rseni onentrtion of the two experiments (. 2 X lo2 m3 in [5] nd 4X lo2 m3 in [6]) will use different vlues of the stble tive rseni onentrtion. This is strongly pronouned in the liner sle of the sheet rrier
10 JUNGLING et l. : SIMULATION OF IC FABNCATION PROCESSES 165 1s 9 C I ' 12 1' 1' 185 ANNEALING TIME Isel Fig. 23. Mesured dt of the sheet rrier onentrtion of [6] t four tempertures. Results due to the luster model of [5] nd n optiml fit to mesured dt. onentrtion in Fig.23. The differenesrevel tht the dynmi properties of rseni lustering re not yet fully understood in terms of developing model whih overs wide rnge of proesses nd n be instlled in generl proess simultion progrms. e hve used our ode to find vlues for k nd k, to gree with mesurements optimlly. Our omputtions revel tht the vlue of the minimum is minly ffeted by the rtio of k/k, nd tht the time dependene of the minimum n be ontrolled by modifying k, nd k, in the sme wy by ftor. It is worth noting tht the slope is nerly independent of k nd k, nd yields fr too steep lines for m = 3 nd k = 1. After severl ttempts m = 7 nd k = 1 yields the best results. The vlues of k, nd k,, = k,/k, used for the optiml fit of tlhe mesured dt in Fig. 23 re listed below. T ("C) kd (sl> keq (m21> X 5.5 X x 4. X X 1~ X The temperture dependene of the vlues of keq, k,, nd k n be desribed by Arrhenius lws. The 1 C proess ould not be optimied due to the shortge of experimentl dt. The poor greement t the beginning of the nneling step is most likely used by very simple initil solution. e hve ssumed tht ll of the rseni ws tivted by the lser nneling step. Our omputtions, using the mesured dt of Fig. 23, indite trend to lrger luster sies. Sine the vlues of r nerly equls tht of CAS, we ould hve lso set m = 6, k = 2; m = 5, k = 3 or m = 4, k = 4 without hnging the term C.nk nd, therefore, the dynmis of (12) or (13). V. SOME MATHEMATICAL AND NUMERICAL ASPECTS A high level of sophistition on the numeril side is neessry to produe ode whih is useful for the investi gtions of physil problems. Sine the ]min purpose of our ode is the investigtion of physil models, we try to free the user from s mny mthemtil problems s possible, e.g., the mesh in spe nd time is reted nd dpted utomtilly. The retion nd modifitions of the grid re minly influened by the ide to minimie the truntion errors. In the spe domin this implies optiml distribution of the gridpoints, in the time domin the steps hve to be suffiiently smll. The deisiions in both domins re ritil nd represent ompromise between ury nd the omputer pbilities. The distribution of the grid points is ontrolled by the following three riteri. ) Qusiuniformity: (only in the 1D ode) In [ll] it is shown tht the sptil truntion error of qusiuniform mesh redues proportionl to the squre of the grid distnes if the mximum rtio between two djent distnes minus unity is smll ompred to unity. b) The minimition of the seond derivtive of the physil quntities. ) The mximum rtio/differene riterion. The use of qusiuniform mesh ombines the dvntge of the smll disretition error of uniform mesh nd the possibility to umulte gridpoints t domins of physil interest. The high dynmi rnge of the diffusion problem in devie fbrition steps (dopnt onentrtions of interest re in the rnge from 1O1O m3 to lo2 ~m~) leds to problems in the design of sptil mesh. On the one hnd, dose onstny is required even t longtime diffusions, on the other hnd the dopnt distribution in the viinity of pn juntion is of speil interest for the eletril performne of the devie. Clultions reveled tht b) retes grid whih permits long time inert diffusions without ny loss of dopnts. Unfortuntely this riterion ignores the low onentrtion domim whih re very often the domins of eletril interest. Criterion )sets points in the viinity of steep grdients nd lose to the pn juntions (mu. differene riterion for the potentil) nd ompenstes the disdvntges of riterion b). A minimum mesh length hs been invoked to void too fine grids ner disontinuities (e.g., initil solution of predeposition). Sine the dopnt migrtion during diffusion step leds to strongly vrying profiles, rigid grid seems to be unsuitble for n dvneed simultion progrm. Our mesh is utomtilly modified during the diffusion step. Observtions during simultions revel tht devitions from the initil doses minly our during grid modifitions. If the grid ontins enough points to fulfillll riteri the vritions turn out to be sttistil deviti.ons nd do not ugment beyond the speified vlues. For the omputtion of the time step we use bkwrd differentition formuls (BDF) similr to those proposed in [12]. In the one dimensionl ode we use BDF of sixth order to obtin n optiml trnsient behvior of our omputtion. In the twodimensionl progrm limited omputer resoures onstrin us to use BDF of third order. Here gin the wide dynmi rnge of interest in the physil quntities neessittes modifitions of the ommonly proposed model. The originl error estimtion in
11 1.66 IEEE TRPAISAGTIQNS Oh1 ELECTRON DEVICES, VOL. ED32, NO. 2, FEARIJARY 1985 [12] turns out to be fr too restritive to yield reson,[hle differenes between the qusineutrl pproximtion nd time steps. the ext Poisson eqution re smll in the viinity of pn In the following tble the number of time steps,.[,he juntions nd negligible elsewhere. The differenes beome number of newton itertions, the number of newton it :r smller with inresing durtion of the diffusion step nd tions per time step nd the omputtion times re listeil For with inresing proess temperture. e n, therefore, the omputtions in Fig. 23. The sptil grid onsist:; of reommend the use of the qusineutrl pproximtion pproximtely 12 points, the prmeter is the u.ly whenever fst performne nd limited memory resoures for the trnsient integrtion. neessitte simplifitions of physil models. In our ode we will, nevertheless, use the Poisson eqution nd investigte osionlly the differenes between the two number of number of models t vrious pplitions. number of Newton itertions omputtio~l A seond pplition dels with the phenomenon of time steps itertions per time step time in se rseni lustering. The omprison with published results revels tht our progrm n urtely simulte uptodte models for dynmi rseni lustering without exeeding its speifitions. The strong devitions of time dependent mesurements revel tht the physil bkground of r seni lustering/preipittion is not yet understood nd._ neessittes further investigtions. The ommonly proposed luster models [5],[8], or,[9] re well estblished but The numbers in the tble indite tht there is no li yer do not led to onentrtion rndependent CAS s it ws reltion between the number of time steps nd the orllpupointed out in, e.g., [4]. On the ontrry [lo] proposes tht ttion time nd tht less urte trnsient integril tion preipittion is the physil effet responsible for the difdoes not signifintly derese the omputtion times.. The ferenes between eletrillytive nd the totl rseni lrger time steps led to less urte preditions for the onentrtion. This model would led to doping indepeninitil solution of the Newton itertion. This in turn nt::esdent tive rseni onentrtion. Furthermore, the insittes more Newton itertions steps nd ompensteu. the fluene of different nneling methods on the dynmis of intended redution of omputtion time by using v,ider rseni hs not been investigted stisftorily. The time time steps. dependent mesurements of [6] nd [7] revel ontrditory results for the dynmis of lustering fter lser nneling For the disretition of the spe opertors we USE: the using different tehniques. Our simultions of lser nnelmethods of finite differenes. Equtions identil to (:I>I to ing proess indite tht luster sies between 5 to 7 toms (3) re used in devie simultion. e hve, therelore, led to betltl greement between mesured nd omimplemented gret del of sophistited numeril methputed dt thn luster sie of 3 toms. ods developed for devie simultion into our Code. The formule for the disretition of (1)(3) hve been txl ren From the numeril point of view the first pplitions of our ode lerly revei tht nrgh level of numeril nd from [ll]. An exponentil fitting ftor for the disre: ::i.stion of the flux reltion (2), whih hs been suessli~lly mthemtil sophistition is bsolutely neessry for pplied in the devie simultion, hs lso been used. The sientifi proess simultion progrm. An optiml disreti Newton lineristion of the disretition formuls le(:.>; to tion of the prtil differentil equtions nd mtomti bloktridigonl (1D) or blokpentgonl (;!.D) lly generted grid in spe nd time gurntee miniml eqution system whih is solved by Gussin elimin,;ltion disretition errors nd fine resolutions of sptil nd (1D) nd n itertive solver (2D), respetively. An ir (:retrnsient domins of interest. ment dmping method developed by Deufelhrd [l.i i is implemented to void detrimentl overshooting effets dur ACKNOLEDGMENT ing the itertion proess. VI. SUMMARY AND CONCLUSIONS This pper presents sientifi generl purpose pro2;rm for proess simultions in one nd two dimensions. The number nd kind of physil quntities to be simultd is speified by few prmeter definitions. This llows l or wide rnge of pplitions nd simple modifitions o the physil models under onsidertion. First pplitions revel tht the qusineutrl pproximtion seems to be n exellent pproh to des1::i:ibe the field enhned diffusion in engineering progrms, The The uthors wouldlike to thnk the Interuniversit ie Rehenentrum for the generous mount of omputer resoures it provided. REFERENCES S. M. Hu, Diffusion in hevily doped semiondutors, J. Appl. Phys., vol. 43, no. 4, pp , R. Shrivstv nd A. Mrshk, Chrge neutrlity nd the internl eletri field produed by impurity diffusion, Solidstte Eletron., VO~. 23, pp , 198. P. J. Anthony, Altertion of Diffusion Profiles in semiondutors due to pn juntions, Solidstte Eletron., vol. 25, no. 1, pp
12 JGLING et [. : SIMULATION OF IC FABRICATION PROCESSES 161 A. Lietoil, J. F. Gibbons, nd T.. Sigmon, The solid solubility nd therml behvior of metstble onentrtions of As in Si, Appl. Phys. Lett., vol. 36, no. 9, pp , My 1, 198. M. Y. Tsi, F. F. Morehed, J. E.E. Bglin, nd A. E. Mihel, Shllow juntions by highdose As implnts in Si, J. Appl. Phys., vol: 51, no. 6, pp , June 198. J. Goetlih, P. H. Tsien, G. Henghuber, nd H. Ryssel, CO, lser nneling of ionimplnted silion, Springer Series in Elerrophysis, vol. 11, pp , P. H. Tsien, H. Ryssel, D. Rshenthler, nd I. Ruge, Nd:YAG lser nneling of rseniimplnted silion, J. Appl. Phys., vol. 52(4), Apr E. Guerrero, H. Ptl, R. Tielert, M. Grsserbuer, nd G. Stingeder, Generlied model for the lustering of As dopnts in Si, J. Eletrohem. SOC., vol. 139, no. 8, pp , Aug [9] R. Fir, GR. eber, Effet of omplex formtion on diffusion of rseni in silion, J. Appl. Phys., vol. 44, no. 1, pp , Jn [lo] D. Nobili, A. Arbelos, G. Celotti, nd S. Solmi., Preipittion s the phenomenon responsible for the eletrilly intive rseni in silion J. Eletrohem. So., vol. 13, no. 4, pp , [ll] S. Selberherr, Anlysis nd Simultion of Semiondutor Devies. New York: SpringerVerlg, [12] R. K. Bryton, F. G. Gustvson, nd G. D. Hhtel, A new effiient lgorithm for solving differentillgebri systems using impliit bkwrd differentition formuls, E ro. IEEE, vol. 6, pp ,1972. [13] P. Deufelhrd, A modified Newton method!.!or the solution of illonditioned systems of nonliner equtions with pplitions to multiple shooting, Numer. Mrh., vol. 22, pp. ;!89315, 1974.
Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1)
Green s Theorem Mth 3B isussion Session Week 8 Notes Februry 8 nd Mrh, 7 Very shortly fter you lerned how to integrte singlevrible funtions, you lerned the Fundmentl Theorem of lulus the wy most integrtion
More informationSolutions to Assignment 1
MTHE 237 Fll 2015 Solutions to Assignment 1 Problem 1 Find the order of the differentil eqution: t d3 y dt 3 +t2 y = os(t. Is the differentil eqution liner? Is the eqution homogeneous? b Repet the bove
More informationLecture 27: Diffusion of Ions: Part 2: coupled diffusion of cations and
Leture 7: iffusion of Ions: Prt : oupled diffusion of tions nd nions s desried y NernstPlnk Eqution Tody s topis Continue to understnd the fundmentl kinetis prmeters of diffusion of ions within n eletrilly
More informationA Mathematical Model for UnemploymentTaking an Action without Delay
Advnes in Dynmil Systems nd Applitions. ISSN 97353 Volume Number (7) pp. 8 Reserh Indi Publitions http://www.ripublition.om A Mthemtil Model for UnemploymentTking n Ation without Dely Gulbnu Pthn Diretorte
More informationLecture 1  Introduction and Basic Facts about PDEs
* 18.15  Introdution to PDEs, Fll 004 Prof. Gigliol Stffilni Leture 1  Introdution nd Bsi Fts bout PDEs The Content of the Course Definition of Prtil Differentil Eqution (PDE) Liner PDEs VVVVVVVVVVVVVVVVVVVV
More information6.5 Improper integrals
Eerpt from "Clulus" 3 AoPS In. www.rtofprolemsolving.om 6.5. IMPROPER INTEGRALS 6.5 Improper integrls As we ve seen, we use the definite integrl R f to ompute the re of the region under the grph of y =
More informationProject 6: Minigoals Towards Simplifying and Rewriting Expressions
MAT 51 Wldis Projet 6: Minigols Towrds Simplifying nd Rewriting Expressions The distriutive property nd like terms You hve proly lerned in previous lsses out dding like terms ut one prolem with the wy
More informationLine Integrals and Entire Functions
Line Integrls nd Entire Funtions Defining n Integrl for omplex Vlued Funtions In the following setions, our min gol is to show tht every entire funtion n be represented s n everywhere onvergent power series
More informationT b a(f) [f ] +. P b a(f) = Conclude that if f is in AC then it is the difference of two monotone absolutely continuous functions.
Rel Vribles, Fll 2014 Problem set 5 Solution suggestions Exerise 1. Let f be bsolutely ontinuous on [, b] Show tht nd T b (f) P b (f) f (x) dx [f ] +. Conlude tht if f is in AC then it is the differene
More informationGlobal alignment. Genome Rearrangements Finding preserved genes. Lecture 18
Computt onl Biology Leture 18 Genome Rerrngements Finding preserved genes We hve seen before how to rerrnge genome to obtin nother one bsed on: Reversls Knowledge of preserved bloks (or genes) Now we re
More informationLesson 2: The Pythagorean Theorem and Similar Triangles. A Brief Review of the Pythagorean Theorem.
27 Lesson 2: The Pythgoren Theorem nd Similr Tringles A Brief Review of the Pythgoren Theorem. Rell tht n ngle whih mesures 90º is lled right ngle. If one of the ngles of tringle is right ngle, then we
More informationChapter 4 StateSpace Planning
Leture slides for Automted Plnning: Theory nd Prtie Chpter 4 StteSpe Plnning Dn S. Nu CMSC 722, AI Plnning University of Mrylnd, Spring 2008 1 Motivtion Nerly ll plnning proedures re serh proedures Different
More information(a) A partition P of [a, b] is a finite subset of [a, b] containing a and b. If Q is another partition and P Q, then Q is a refinement of P.
Chpter 7: The Riemnn Integrl When the derivtive is introdued, it is not hrd to see tht the it of the differene quotient should be equl to the slope of the tngent line, or when the horizontl xis is time
More informationModelling the Electrolyte Flow in a Fullscale Copper Electrorefining Tankhouse Cell
Modelling the Eletrolyte Flow in Fullsle Copper Eletrorefining Tnkhouse Cell Andres Kemminger, Andres Ludwig Montnuniversitet Leoben Deprtment Metllurgy, Chir of Simultion nd Modelling of Metllurgil Proesses
More informationFor the flux through a surface: Ch.24 Gauss s Law In last chapter, to calculate electric filede at a give location: q For point charges: K i r 2 ˆr
Ch.24 Guss s Lw In lst hpter, to lulte eletri filed t give lotion: q For point hrges: K i e r 2 ˆr i dq For ontinuous hrge distributions: K e r 2 ˆr However, for mny situtions with symmetri hrge distribution,
More informationPhysics 505 Homework No. 11 Solutions S111
Physis 55 Homework No 11 s S111 1 This problem is from the My, 24 Prelims Hydrogen moleule Consider the neutrl hydrogen moleule, H 2 Write down the Hmiltonin keeping only the kineti energy terms nd the
More informationf (x)dx = f(b) f(a). a b f (x)dx is the limit of sums
Green s Theorem If f is funtion of one vrible x with derivtive f x) or df dx to the Fundmentl Theorem of lulus, nd [, b] is given intervl then, ording This is not trivil result, onsidering tht b b f x)dx
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06073 HIKARI Ltd, www.mhikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationThe EmissionAbsorption of Energy analyzed by QuantumRelativity. Abstract
The missionabsorption of nergy nlyzed by QuntumReltivity Alfred Bennun* & Néstor Ledesm** Abstrt The uslity horizon llows progressive quntifition, from n initil nk prtile, whih yields its energy s blk
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationMACsolutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences  Finite Elements  Finite Volumes  Boundry Elements MACsolutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be relvlues nd smooth The pproximtion of n integrl by numericl
More informationMA10207B: ANALYSIS SECOND SEMESTER OUTLINE NOTES
MA10207B: ANALYSIS SECOND SEMESTER OUTLINE NOTES CHARLIE COLLIER UNIVERSITY OF BATH These notes hve been typeset by Chrlie Collier nd re bsed on the leture notes by Adrin Hill nd Thoms Cottrell. These
More informationIntroduction to Olympiad Inequalities
Introdution to Olympid Inequlities Edutionl Studies Progrm HSSP Msshusetts Institute of Tehnology Snj Simonovikj Spring 207 Contents Wrm up nd AmGm inequlity 2. Elementry inequlities......................
More informationChem Homework 11 due Monday, Apr. 28, 2014, 2 PM
Chem 44  Homework due ondy, pr. 8, 4, P.. . Put this in eq 8.4 terms: E m = m h /m e L for L=d The degenery in the ring system nd the inresed sping per level (4x bigger) mkes the sping between the HOO
More informationQUADRATIC EQUATION. Contents
QUADRATIC EQUATION Contents Topi Pge No. Theory 004 Exerise  0509 Exerise  093 Exerise  3 45 Exerise  4 6 Answer Key 78 Syllus Qudrti equtions with rel oeffiients, reltions etween roots nd oeffiients,
More informationPAIR OF LINEAR EQUATIONS IN TWO VARIABLES
PAIR OF LINEAR EQUATIONS IN TWO VARIABLES. Two liner equtions in the sme two vriles re lled pir of liner equtions in two vriles. The most generl form of pir of liner equtions is x + y + 0 x + y + 0 where,,,,,,
More informationa) Read over steps (1) (4) below and sketch the path of the cycle on a P V plot on the graph below. Label all appropriate points.
Prole 3: Crnot Cyle of n Idel Gs In this prole, the strting pressure P nd volue of n idel gs in stte, re given he rtio R = / > of the volues of the sttes nd is given Finlly onstnt γ = 5/3 is given You
More information1.3 SCALARS AND VECTORS
Bridge Course Phy I PUC 24 1.3 SCLRS ND VECTORS Introdution: Physis is the study of nturl phenomen. The study of ny nturl phenomenon involves mesurements. For exmple, the distne etween the plnet erth nd
More informationThermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report
Therml iffusivity Pul Hughes eprtment of Physics nd Astronomy The University of nchester nchester 3 9PL Second Yer Lbortory Report Nov 4 Abstrct We investigted the therml diffusivity of cylindricl block
More informationwhere the box contains a finite number of gates from the given collection. Examples of gates that are commonly used are the following: a b
CS 2942 9/11/04 Quntum Ciruit Model, SolovyKitev Theorem, BQP Fll 2004 Leture 4 1 Quntum Ciruit Model 1.1 Clssil Ciruits  Universl Gte Sets A lssil iruit implements multioutput oolen funtion f : {0,1}
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationSIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WITH VARIOUS TYPES OF COLUMN BASES
Advned Steel Constrution Vol., No., pp. 788 () 7 SIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WIT VARIOUS TYPES OF COLUMN BASES J. ent sio Assoite Professor, Deprtment of Civil nd Environmentl
More informationLine Integrals. Partitioning the Curve. Estimating the Mass
Line Integrls Suppose we hve curve in the xy plne nd ssocite density δ(p ) = δ(x, y) t ech point on the curve. urves, of course, do not hve density or mss, but it my sometimes be convenient or useful to
More informationPhysics 9 Fall 2011 Homework 2  Solutions Friday September 2, 2011
Physics 9 Fll 0 Homework  s Fridy September, 0 Mke sure your nme is on your homework, nd plese box your finl nswer. Becuse we will be giving prtil credit, be sure to ttempt ll the problems, even if you
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More informationProving the Pythagorean Theorem
Proving the Pythgoren Theorem W. Bline Dowler June 30, 2010 Astrt Most people re fmilir with the formul 2 + 2 = 2. However, in most ses, this ws presented in lssroom s n solute with no ttempt t proof or
More informationChapter Gauss Quadrature Rule of Integration
Chpter 7. Guss Qudrture Rule o Integrtion Ater reding this hpter, you should e le to:. derive the Guss qudrture method or integrtion nd e le to use it to solve prolems, nd. use Guss qudrture method to
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationCore 2 Logarithms and exponentials. Section 1: Introduction to logarithms
Core Logrithms nd eponentils Setion : Introdution to logrithms Notes nd Emples These notes ontin subsetions on Indies nd logrithms The lws of logrithms Eponentil funtions This is n emple resoure from MEI
More informationON THE INEQUALITY OF THE DIFFERENCE OF TWO INTEGRAL MEANS AND APPLICATIONS FOR PDFs
ON THE INEQUALITY OF THE DIFFERENCE OF TWO INTEGRAL MEANS AND APPLICATIONS FOR PDFs A.I. KECHRINIOTIS AND N.D. ASSIMAKIS Deprtment of Eletronis Tehnologil Edutionl Institute of Lmi, Greee EMil: {kehrin,
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationArrow s Impossibility Theorem
Rep Voting Prdoxes Properties Arrow s Theorem Arrow s Impossiility Theorem Leture 12 Arrow s Impossiility Theorem Leture 12, Slide 1 Rep Voting Prdoxes Properties Arrow s Theorem Leture Overview 1 Rep
More informationCalculus Cheat Sheet. Integrals Definitions. where F( x ) is an antiderivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx
Clulus Chet Sheet Integrls Definitions Definite Integrl: Suppose f ( ) is ontinuous AntiDerivtive : An ntiderivtive of f ( ) on [, ]. Divide [, ] into n suintervls of is funtion, F( ), suh tht F = f.
More informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More informationForces on curved surfaces Buoyant force Stability of floating and submerged bodies
Stti Surfe ores Stti Surfe ores 8m wter hinge? 4 m ores on plne res ores on urved surfes Buont fore Stbilit of floting nd submerged bodies ores on Plne res Two tpes of problems Horizontl surfes (pressure
More information6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 CH 3. CH 3 C a. NMR spectroscopy. Different types of NMR
6.. Spetrosopy NMR spetrosopy Different types of NMR NMR spetrosopy involves intertion of mterils with the lowenergy rdiowve region of the eletromgneti spetrum NMR spetrosopy is the sme tehnology s tht
More informationMore Properties of the Riemann Integral
More Properties of the Riemnn Integrl Jmes K. Peterson Deprtment of Biologil Sienes nd Deprtment of Mthemtil Sienes Clemson University Februry 15, 2018 Outline More Riemnn Integrl Properties The Fundmentl
More informationSection 4.4. Green s Theorem
The Clulus of Funtions of Severl Vriles Setion 4.4 Green s Theorem Green s theorem is n exmple from fmily of theorems whih onnet line integrls (nd their higherdimensionl nlogues) with the definite integrls
More information6.1 Definition of the Riemann Integral
6 The Riemnn Integrl 6. Deinition o the Riemnn Integrl Deinition 6.. Given n intervl [, b] with < b, prtition P o [, b] is inite set o points {x, x,..., x n } [, b], lled grid points, suh tht x =, x n
More informationKinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)
Introduction Kinemtic Wves These re wves which result from the conservtion eqution E t + I = 0 (1) where E represents sclr density field nd I, its outer flux. The onedimensionl form of (1) is E t + I
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationJackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The twodimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero
More informationMatrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix
tries Definition of tri mtri is regulr rry of numers enlosed inside rkets SCHOOL OF ENGINEERING & UIL ENVIRONEN Emple he following re ll mtries: ), ) 9, themtis ), d) tries Definition of tri Size of tri
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationSTEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.
STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More information6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 H 3 CH3 C. NMR spectroscopy. Different types of NMR
6.. Spetrosopy NMR spetrosopy Different types of NMR NMR spetrosopy involves intertion of mterils with the lowenergy rdiowve region of the eletromgneti spetrum NMR spetrosopy is the sme tehnology s tht
More informationSECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS
SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS INSITU PROBING OF DOMAIN POLING IN Bi 4 Ti 3 O 12 THIN FILMS BY OPTICAL SECOND HARMONIC GENERATION YANIV BARAD, VENKATRAMAN GOPALAN Mterils Reserh Lortory
More informationThe RiemannStieltjes Integral
Chpter 6 The RiemnnStieltjes Integrl 6.1. Definition nd Eistene of the Integrl Definition 6.1. Let, b R nd < b. ( A prtition P of intervl [, b] is finite set of points P = { 0, 1,..., n } suh tht = 0
More informationSTURM LIOUVILLE OPERATORS WITH MEASUREVALUED COEFFICIENTS. 1. Introduction. + q(x)y(x) = zr(x)y(x)
STUM LIOUVILLE OPEATOS WITH MEASUEVALUED COEFFICIENTS JONATHAN ECKHADT AND GEALD TESCHL Abstrt. We give omprehensive tretment of Sturm Liouville opertors whose oeffiients re mesures inluding full disussion
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationEffects of Drought on the Performance of Two Hybrid Bluegrasses, Kentucky Bluegrass and Tall Fescue
TITLE: OBJECTIVE: AUTHOR: SPONSORS: Effets of Drought on the Performne of Two Hyrid Bluegrsses, Kentuky Bluegrss nd Tll Fesue Evlute the effets of drought on the visul qulity nd photosynthesis in two hyrid
More informationAlgorithm Design and Analysis
Algorithm Design nd Anlysis LECTURE 8 Mx. lteness ont d Optiml Ching Adm Smith 9/12/2008 A. Smith; sed on slides y E. Demine, C. Leiserson, S. Rskhodnikov, K. Wyne Sheduling to Minimizing Lteness Minimizing
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the righthnd side limit equls to the lefthnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationNANOSCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE
NnoSle Effets in Clindril Contts Sri et l. NANOSCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE Ö. T. Sri, G. G. Adms, S. Müftü Mehnil Engineering Deprtment Northestern
More informationThe Thermodynamics of Aqueous Electrolyte Solutions
18 The Thermodynmics of Aqueous Electrolyte Solutions As discussed in Chpter 10, when slt is dissolved in wter or in other pproprite solvent, the molecules dissocite into ions. In queous solutions, strong
More informationProblems set # 3 Physics 169 February 24, 2015
Prof. Anhordoqui Problems set # 3 Physis 169 Februry 4, 015 1. A point hrge q is loted t the enter of uniform ring hving liner hrge density λ nd rdius, s shown in Fig. 1. Determine the totl eletri flux
More informationNONDETERMINISTIC FSA
Tw o types of nondeterminism: NONDETERMINISTIC FS () Multiple strtsttes; strtsttes S Q. The lnguge L(M) ={x:x tkes M from some strtstte to some finlstte nd ll of x is proessed}. The string x = is
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationDevelopment of new equation to estimate the maximum soil depth by using the safety factor
Lndslides nd Engineered Slopes. Experiene, Theory nd Prtie Avers et l. (Eds) 2016 Assoiione Geoteni Itlin, Rome, Itly, ISBN 9781138029880 Development of ne eqution to estimte the mximum soil depth
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationIntro to Nuclear and Particle Physics (5110)
Intro to Nucler nd Prticle Physics (5110) Feb, 009 The Nucler Mss Spectrum The Liquid Drop Model //009 1 E(MeV) n n(n1)/ E/[ n(n1)/] (MeV/pir) 1 C 16 O 0 Ne 4 Mg 7.7 14.44 19.17 8.48 4 5 6 6 10 15.4.41
More informationInterpreting Integrals and the Fundamental Theorem
Interpreting Integrls nd the Fundmentl Theorem Tody, we go further in interpreting the mening of the definite integrl. Using Units to Aid Interprettion We lredy know tht if f(t) is the rte of chnge of
More informationH 4 H 8 N 2. Example 1 A compound is found to have an accurate relative formula mass of It is thought to be either CH 3.
. Spetrosopy Mss spetrosopy igh resolution mss spetrometry n e used to determine the moleulr formul of ompound from the urte mss of the moleulr ion For exmple, the following moleulr formuls ll hve rough
More informationReflection Property of a Hyperbola
Refletion Propert of Hperol Prefe The purpose of this pper is to prove nltill nd to illustrte geometrill the propert of hperol wherein r whih emntes outside the onvit of the hperol, tht is, etween the
More informationA027 Uncertainties in Local Anisotropy Estimation from Multioffset VSP Data
A07 Uncertinties in Locl Anisotropy Estimtion from Multioffset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.
More information2.4 Linear Inequalities and Problem Solving
Section.4 Liner Inequlities nd Problem Solving 77.4 Liner Inequlities nd Problem Solving S 1 Use Intervl Nottion. Solve Liner Inequlities Using the Addition Property of Inequlity. 3 Solve Liner Inequlities
More informationTerminal Velocity and Raindrop Growth
Terminl Velocity nd Rindrop Growth Terminl velocity for rindrop represents blnce in which weight mss times grvity is equl to drg force. F 3 π3 ρ L g in which is drop rdius, g is grvittionl ccelertion,
More informationBest Approximation in the 2norm
Jim Lmbers MAT 77 Fll Semester 111 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationAMATH 731: Applied Functional Analysis Fall Some basics of integral equations
AMATH 731: Applied Functionl Anlysis Fll 2009 1 Introduction Some bsics of integrl equtions An integrl eqution is n eqution in which the unknown function u(t) ppers under n integrl sign, e.g., K(t, s)u(s)
More informationDIRECT CURRENT CIRCUITS
DRECT CURRENT CUTS ELECTRC POWER Consider the circuit shown in the Figure where bttery is connected to resistor R. A positive chrge dq will gin potentil energy s it moves from point to point b through
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview IntroductionModelling Bsic concepts to understnd n ODE. Description nd properties
More informationNumbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point
GCSE C Emple 7 Work out 9 Give your nswer in its simplest form Numers n inies Reiprote mens invert or turn upsie own The reiprol of is 9 9 Mke sure you only invert the frtion you re iviing y 7 You multiply
More informationAvailable online at ScienceDirect. Procedia Engineering 120 (2015 ) EUROSENSORS 2015
Aville online t www.sienediret.om SieneDiret Proedi Engineering 10 (015 ) 887 891 EUROSENSORS 015 A Fesiility Study for SelfOsillting Loop for Three DegreeofFreedom Coupled MEMS Resontor Fore Sensor
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationPDE Notes. Paul Carnig. January ODE s vs PDE s 1
PDE Notes Pul Crnig Jnury 2014 Contents 1 ODE s vs PDE s 1 2 Section 1.2 Het diffusion Eqution 1 2.1 Fourier s w of Het Conduction............................. 2 2.2 Energy Conservtion.....................................
More information#6A&B Magnetic Field Mapping
#6A& Mgnetic Field Mpping Gol y performing this lb experiment, you will: 1. use mgnetic field mesurement technique bsed on Frdy s Lw (see the previous experiment),. study the mgnetic fields generted by
More informationPolynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230
Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given
More informationBIFURCATIONS IN ONEDIMENSIONAL DISCRETE SYSTEMS
BIFRCATIONS IN ONEDIMENSIONAL DISCRETE SYSTEMS FRANCESCA AICARDI In this lesson we will study the simplest dynmicl systems. We will see, however, tht even in this cse the scenrio of different possible
More informationOn the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou
On the Uncertinty of Sensors Bsed on Mgnetic Effects E. ristoforou, E. Kyfs, A. Kten, DM Kepptsoglou Ntionl Technicl University of Athens, Zogrfou Cmpus, Athens 1578, Greece Tel: +3177178, Fx: +3177119,
More informationIf only one fertilizer x is used, the dependence of yield z(x) on x first was given by Mitscherlich (1909) in form of the differential equation
Mitsherlih s Lw: Generliztion with severl Fertilizers Hns Shneeerger Institute of Sttistis, University of ErlngenNürnerg, Germny 00, 5 th August Astrt: It is shown, tht the ropyield z in dependene on
More informationOrthogonal Polynomials and LeastSquares Approximations to Functions
Chpter Orthogonl Polynomils nd LestSqures Approximtions to Functions **4/5/3 ET. Discrete LestSqures Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More information( ) as a fraction. Determine location of the highest
AB/ Clulus Exm Review Sheet Solutions A Prelulus Type prolems A1 A A3 A4 A5 A6 A7 This is wht you think of doing Find the zeros of f( x) Set funtion equl to Ftor or use qudrti eqution if qudrti Grph to
More informationMeasuring Electron Work Function in Metal
n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd
More information