Processes Using Advanced Physical and

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Processes Using Advanced Physical and"

Transcription

1 156 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. ED-32, 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 pplitwo-dimensionl 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 spe-hrge pproximtionndthesolution (1 the physil models (e.g., vrying the number nd kind of ext Poisson eqution restudied by exmples of As-B 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 up-to-dte 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:nerl-purpose 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, 14-ein, 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 qusi-neutrl 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., vnies-interstitils, 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 p-n 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. ED-32, 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 field-enhnement ftor o' the diffusion oeffiient (7). h=l+&.($+l) -.5 If p-n juntions re present, the nlyti solution :' the qusi-neutrl 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 field-enhnement, 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 qusi-neutrl 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 qusi-nelltrl 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 qusi-neutrl pproximtion hs been used to reple the Poisson eqution. The boundry between the two domins wsposed t the end of the spe-h rge region. In [2] n ext solution of the Poisson equtim is presented; however, the Computtions do not involve p-n juntion. Corret implementtion of field-enhned 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 p-n 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'~ m-2 with 9 kev, LSS, Gussin; B: 4X 1'5 m-2 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' m-2 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. 3-5 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 field-ontrolled period (. 1 s) followed by n mbipolr diffusion. The strong field-ontrolled 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 p-domin 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 3-s 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 r-7 e 1. 5 I- 1 t DEPTH [pml Fig. 3. Arseni nd boron distributions fter 1-s inert diffusion t 95OoC due to the models P nd Q..1.9.ll.12 DEPTH [pl Fig. 5. Arseni nd boron distributions fter 1-s 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 field-ontrolled flux. Figs. 3-5 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 (qusi-neutrl pproximtion (4),(5),(8))leds to stronger eletril field. Therefore the field-ontrolled 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. 3-5 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 n-side of the p-n juntions lose to the juntion. Differenes of the dopnt onentrtions

5 16 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. ED-32, 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 p-n 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 p-n 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 p-n 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 n-doped side. Sine the pek vlue of the boron onentrtion is of the sme order of mgnitude s the intrinsi number (6X1Ol8 m-3 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 p-n 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 qusi-neutrl 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 p-n 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 p-n juntion no differenes n be seen. The depth of the p-n 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 qusi-neutrl pproximtion insted of the ext Poisson eqution when the durtion of the diffusion step is long enough to exeed the first, field-ontrolled 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. Time-de 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 m-2, 4 kev. Fig. 9. Trnsient behvior of the positive spe hrge density whih is defined s the mximum of nd n-p+ 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 diffusion-ontrolled 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'p-n 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 m-3 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 field-ontrolled 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 qusi-neutrl pproximtion desribes the dopnt migrtion qulittively well in the viinity of p-n 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 two-dimensionl 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. ED-32, NO. 2, FEBRUARY Fig. 12. Two-dimensionl distribution of the eletrostti potentil id'ter n inert diffusion of 1 s t 1OOO"C. Fig. 14. Two-dimensionl boron distribution fter n inert diffusion of 1 s t 1OOO"C. 3% $< t x i? 8% LI. T? Fig. 15. Two-dimensionl boron distribution fter n inert diffusion of 5 s t 1OOO"C. Fig. 13. Two-dimensionl 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. Two-dimensionl 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- U-shped 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 U-shped 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 two-dimensionl 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. Two-dimensionl 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 2-min 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 25-nm 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 qusi-neutrl 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 m-3 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. ED-32, NO. 2, FEBRUARY 1985 Fig. 2. Trnsient dependene of the lustered rseni onentrtion during 2-min nneling step t 1 C. 1 OZ2 Fig. 22. Trnsient dependene of the eletrilly tive rseni 6-min 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 x-xis shows the nneling time in logrithmi sle, the y-xis the sheet rrier onentrtion per m-2 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 m-3 in [5] nd 4X lo2 m-3 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 (s-l> 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. ) Qusi-uniformity: (only in the 1-D ode) In [ll] it is shown tht the sptil truntion error of qusi-uniform 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 qusi-uniform 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 m-3 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 p-n 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 p-n 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 two-dimensionl 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. ED-32, NO. 2, FEARIJARY 1985 [12] turns out to be fr too restritive to yield reson,[hle differenes between the qusi-neutrl pproximtion nd time steps. the ext Poisson eqution re smll in the viinity of p-n 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 qusi-neutrl 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 up-to-dte 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 blok-tridigonl (1-D) or blok-pentgonl (;!.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 (1-D) nd n itertive solver (2-D), 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 qusi-neutrl 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, Solid-stte Eletron., VO~. 23, pp , 198. P. J. Anthony, Altertion of Diffusion Profiles in semiondutors due to p-n juntions, Solid-stte 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 high-dose 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 ion-implnted silion, Springer Series in Elerrophysis, vol. 11, pp , P. H. Tsien, H. Ryssel, D. Rshenthler, nd I. Ruge, Nd:YAG lser nneling of rseni-implnted 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, G-R. 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: Springer-Verlg, [12] R. K. Bryton, F. G. Gustvson, nd G. D. Hhtel, A new effiient lgorithm for solving differentil-lgebri systems using impliit bkwrd differentition formuls, E ro. IEEE, vol. 6, pp ,1972. [13] P. Deufelhrd, A modified Newton method!.!or the solution of ill-onditioned systems of nonliner equtions with pplitions to multiple shooting, Numer. Mrh., vol. 22, pp. ;!89-315, 1974.

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1)

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 single-vrible funtions, you lerned the Fundmentl Theorem of lulus the wy most integrtion

More information

Solutions to Assignment 1

Solutions 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 information

Lecture 27: Diffusion of Ions: Part 2: coupled diffusion of cations and

Lecture 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 Nernst-Plnk Eqution Tody s topis Continue to understnd the fundmentl kinetis prmeters of diffusion of ions within n eletrilly

More information

A Mathematical Model for Unemployment-Taking an Action without Delay

A Mathematical Model for Unemployment-Taking an Action without Delay Advnes in Dynmil Systems nd Applitions. ISSN 973-53 Volume Number (7) pp. -8 Reserh Indi Publitions http://www.ripublition.om A Mthemtil Model for Unemployment-Tking n Ation without Dely Gulbnu Pthn Diretorte

More information

Lecture 1 - Introduction and Basic Facts about PDEs

Lecture 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 information

6.5 Improper integrals

6.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 information

Project 6: Minigoals Towards Simplifying and Rewriting Expressions

Project 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 information

Line Integrals and Entire Functions

Line 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 information

T 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.

T 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 information

Global alignment. Genome Rearrangements Finding preserved genes. Lecture 18

Global 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 information

Lesson 2: The Pythagorean Theorem and Similar Triangles. A Brief Review of the Pythagorean Theorem.

Lesson 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 information

Chapter 4 State-Space Planning

Chapter 4 State-Space Planning Leture slides for Automted Plnning: Theory nd Prtie Chpter 4 Stte-Spe 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.

(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 information

Modelling the Electrolyte Flow in a Full-scale Copper Electrorefining Tankhouse Cell

Modelling the Electrolyte Flow in a Full-scale Copper Electrorefining Tankhouse Cell Modelling the Eletrolyte Flow in Full-sle Copper Eletrorefining Tnkhouse Cell Andres Kemminger, Andres Ludwig Montnuniversitet Leoben Deprtment Metllurgy, Chir of Simultion nd Modelling of Metllurgil Proesses

More information

For 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

For 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 information

Physics 505 Homework No. 11 Solutions S11-1

Physics 505 Homework No. 11 Solutions S11-1 Physis 55 Homework No 11 s S11-1 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 information

f (x)dx = f(b) f(a). a b f (x)dx is the limit of sums

f (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 information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

The Emission-Absorption of Energy analyzed by Quantum-Relativity. Abstract

The Emission-Absorption of Energy analyzed by Quantum-Relativity. Abstract The mission-absorption of nergy nlyzed by Quntum-Reltivity Alfred Bennun* & Néstor Ledesm** Abstrt The uslity horizon llows progressive quntifition, from n initil nk prtile, whih yields its energy s blk

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 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 information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

Lecture 14: Quadrature

Lecture 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 rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

MA10207B: ANALYSIS SECOND SEMESTER OUTLINE NOTES

MA10207B: 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 information

Introduction to Olympiad Inequalities

Introduction to Olympiad Inequalities Introdution to Olympid Inequlities Edutionl Studies Progrm HSSP Msshusetts Institute of Tehnology Snj Simonovikj Spring 207 Contents Wrm up nd Am-Gm inequlity 2. Elementry inequlities......................

More information

Chem Homework 11 due Monday, Apr. 28, 2014, 2 PM

Chem 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 information

QUADRATIC EQUATION. Contents

QUADRATIC EQUATION. Contents QUADRATIC EQUATION Contents Topi Pge No. Theory 0-04 Exerise - 05-09 Exerise - 09-3 Exerise - 3 4-5 Exerise - 4 6 Answer Key 7-8 Syllus Qudrti equtions with rel oeffiients, reltions etween roots nd oeffiients,

More information

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

PAIR 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 information

a) 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.

a) 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 information

1.3 SCALARS AND VECTORS

1.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 information

Thermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report

Thermal 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 information

where the box contains a finite number of gates from the given collection. Examples of gates that are commonly used are the following: a b

where 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 294-2 9/11/04 Quntum Ciruit Model, Solovy-Kitev Theorem, BQP Fll 2004 Leture 4 1 Quntum Ciruit Model 1.1 Clssil Ciruits - Universl Gte Sets A lssil iruit implements multi-output oolen funtion f : {0,1}

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL 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 information

SIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WITH VARIOUS TYPES OF COLUMN BASES

SIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WITH VARIOUS TYPES OF COLUMN BASES Advned Steel Constrution Vol., No., pp. 7-88 () 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 information

Line Integrals. Partitioning the Curve. Estimating the Mass

Line 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 information

Physics 9 Fall 2011 Homework 2 - Solutions Friday September 2, 2011

Physics 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 information

Numerical integration

Numerical 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 information

Math 8 Winter 2015 Applications of Integration

Math 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 information

The Wave Equation I. MA 436 Kurt Bryan

The 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 information

Proving the Pythagorean Theorem

Proving 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 information

Chapter Gauss Quadrature Rule of Integration

Chapter 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 information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. 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 information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL 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 information

Core 2 Logarithms and exponentials. Section 1: Introduction to logarithms

Core 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 information

ON 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 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 information

Lecture 19: Continuous Least Squares Approximation

Lecture 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 information

Arrow s Impossibility Theorem

Arrow 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 information

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative 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 Anti-Derivtive : An nti-derivtive of f ( ) on [, ]. Divide [, ] into n suintervls of is funtion, F( ), suh tht F = f.

More information

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Section 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 information

Forces on curved surfaces Buoyant force Stability of floating and submerged bodies

Forces 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 information

6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 CH 3. CH 3 C a. NMR spectroscopy. Different types of NMR

6.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 information

More Properties of the Riemann Integral

More 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 information

Section 4.4. Green s Theorem

Section 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 higher-dimensionl nlogues) with the definite integrls

More information

6.1 Definition of the Riemann Integral

6.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 information

Kinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)

Kinematic 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 one-dimensionl form of (1) is E t + I

More information

CS667 Lecture 6: Monte Carlo Integration 02/10/05

CS667 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 information

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 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 two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero

More information

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

Matrices 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 information

Best Approximation. Chapter The General Case

Best 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 information

STEP 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. 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 information

8 Laplace s Method and Local Limit Theorems

8 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 information

Math 270A: Numerical Linear Algebra

Math 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 information

6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 H 3 CH3 C. NMR spectroscopy. Different types of NMR

6.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 information

SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS

SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS IN-SITU 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 information

The Riemann-Stieltjes Integral

The Riemann-Stieltjes Integral Chpter 6 The Riemnn-Stieltjes 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 information

STURM LIOUVILLE OPERATORS WITH MEASURE-VALUED COEFFICIENTS. 1. Introduction. + q(x)y(x) = zr(x)y(x)

STURM LIOUVILLE OPERATORS WITH MEASURE-VALUED COEFFICIENTS. 1. Introduction. + q(x)y(x) = zr(x)y(x) STUM LIOUVILLE OPEATOS WITH MEASUE-VALUED COEFFICIENTS JONATHAN ECKHADT AND GEALD TESCHL Abstrt. We give omprehensive tretment of Sturm Liouville opertors whose oeffiients re mesures inluding full disussion

More information

3.4 Numerical integration

3.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 information

Effects of Drought on the Performance of Two Hybrid Bluegrasses, Kentucky Bluegrass and Tall Fescue

Effects 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 information

Algorithm Design and Analysis

Algorithm 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 information

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Partial 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 right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles

More information

NANO-SCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE

NANO-SCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE Nno-Sle Effets in Clindril Contts Sri et l. NANO-SCALE 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 information

The Thermodynamics of Aqueous Electrolyte Solutions

The 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 information

Problems set # 3 Physics 169 February 24, 2015

Problems 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 information

NON-DETERMINISTIC FSA

NON-DETERMINISTIC FSA Tw o types of non-determinism: NON-DETERMINISTIC FS () Multiple strt-sttes; strt-sttes S Q. The lnguge L(M) ={x:x tkes M from some strt-stte to some finl-stte nd ll of x is proessed}. The string x = is

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: 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 information

Development of new equation to estimate the maximum soil depth by using the safety factor

Development 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 978-1-138-02988-0 Development of ne eqution to estimte the mximum soil depth

More information

The 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 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 information

Intro to Nuclear and Particle Physics (5110)

Intro 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(n-1)/ E/[ n(n-1)/] (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 information

Interpreting Integrals and the Fundamental Theorem

Interpreting 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 information

H 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.

H 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 information

Reflection Property of a Hyperbola

Reflection 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 information

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data A07 Uncertinties in Locl Anisotropy Estimtion from Multi-offset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.

More information

2.4 Linear Inequalities and Problem Solving

2.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 information

Terminal Velocity and Raindrop Growth

Terminal 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 information

Best Approximation in the 2-norm

Best Approximation in the 2-norm Jim Lmbers MAT 77 Fll Semester 1-11 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 information

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations

AMATH 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 information

DIRECT CURRENT CIRCUITS

DIRECT 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 information

Ordinary differential equations

Ordinary differential equations Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties

More information

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

Numbers 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 information

Available online at ScienceDirect. Procedia Engineering 120 (2015 ) EUROSENSORS 2015

Available 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 Self-Osillting Loop for Three Degreeof-Freedom Coupled MEMS Resontor Fore Sensor

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 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 information

Lecture 1. Functional series. Pointwise and uniform convergence.

Lecture 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 information

PDE Notes. Paul Carnig. January ODE s vs PDE s 1

PDE 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&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 information

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Polynomial 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 information

BIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS

BIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS BIFRCATIONS IN ONE-DIMENSIONAL 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 information

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

On 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 information

If 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

If 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 Erlngen-Nürnerg, Germny 00, 5 th August Astrt: It is shown, tht the rop-yield z in dependene on

More information

Orthogonal Polynomials and Least-Squares Approximations to Functions

Orthogonal Polynomials and Least-Squares Approximations to Functions Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres 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

( ) 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 information

Measuring Electron Work Function in Metal

Measuring Electron Work Function in Metal n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd

More information