Applied Statistical Mechanics Lecture Note - 13 Molecular Dynamics Simulation
|
|
- Reginald Williams
- 5 years ago
- Views:
Transcription
1 Appled Statstcal Mechancs Lectue Note - 3 Molecula Dynamcs Smulaton 고려대학교화공생명공학과강정원
2 Contents I. Basc Molecula Dynamcs Smulaton Method II. Popetes Calculatons n MD III. MD n Othe Ensembles
3 I. Basc MD Smulaton -MC vs. MD MC MD Pobablstc smulaton technque Lmtatons eque the knowledge of an equlbum dstbuton goous samplng of lage numbe of possble phase-space gves only confguatonal popetes (not dynamc popetes!) Detemnstc smulaton technque Fully numecal fomalsm numecal soluton of N-body system A d N A( d N N exp N { βu ( )} N { βu ( )} )exp
4 I. Basc MD Smulaton - The Idea Follow the exactly same pocedue as eal expements Pepae sample pepae N patcles solve equaton of motons Connect sample to measung nstuments (e.g. themomete, vscomete, ) afte equlbaton tme, actual measuement begns Measue the popety of nteest fo a cetan tme nteval aveage popetes Example : measuement of tempeatue α mv k B T T ( t) mv ( t) k N B f
5 I. Basc MD Smulaton - Equaton of Moton Classcal Newton s equaton of moton Thee fomulaton Newtonan Lagangan Hamltonan Hamltonan pefeed fo many-body systems soluton of N dffeental equatons t m t F p p Soluton methods : Fnte Dffeence Method j j j z y x z y x p p p ),, ( ),, ( F F p p
6 I. Basc MD Smulaton - Velet Algothm Velet (967) : Vey smple, effcent and popula algothm 3 4 m m 3! 3 4 m m 3! ( t + δt) () t + p() t δt + F() t δt + () t δt + O( δt ) ( t δt) () t p() t δt + F() t δt () t δt + O( δt ) m 4 ( t + δt) + ( t δt) ( t) + F( t) δt + O( δt ) m 4 ( t + δt) ( t) ( t δt) + F( t) δt + O( δt ) featue : update wthout calculatng momentum (p)
7 I. Basc MD Smulaton - Leapfog Algothm Hockeny (970), Potte (97) Half-step leap-fog algothm Mathematcally equvalent to Velet algothm m ( t + δt) ( t) + v( t + δt) δt v( t + δt) v( t δt) + F( t) δt ( t + δt) ( t) + ( t t) v δ + F( t) δtδt m
8 . Popetes Calculaton n MD - Eneges Potental enegy Can be calculated dung foce calculaton Knetc enegy K mv
9 . Popetes Calculaton n MD - Pessues In an MD smulaton, calculaton of pessue usng tenso notaton s not the most effcent method. Fo homogeneous systems, thee s smple way to calculate pessue (Ivng and Kkwood, 950) Knetc deal gas tem Confguatonal called Val P ρk V B T + m V v ( t) v j ( t) + j ( t) F j j ( t) j ( t) F j ( t) Calculate when foce update Calculate when velocty update
10 . Popetes Calculaton n MD - Tanspot Popetes Appoaches fo tanspot popetes Method : NEMD (Non-equlbum Molecula Dynamcs) Contnuous addton and emoval of conseved quanttes Gves hgh sgnal-to-nose ato (good statstcs) Method : Equlbum molecula dynamcs Stat wth ansotopc confguaton of mass, momentum and enegy Obseve natual fluctuatons and dsspaton of mass, momentum and enegy Poo sgnal-to-nose aton (poo statstcs) All tanspot popetes can be measued at once
11 . Popetes Calculaton n MD - Tanspot Popetes Dffeental Balance Equaton Mass Enegy Momentum c(, t) t + j 0 T (, t) Dv(, t) c p + q 0 ρ + τ 0 t Dt Consttutve Equatons Fck s Law Foue s Law Newton s Law j D c q k T τ v ρv ) xy y ( x
12 . Popetes Calculaton n MD - Tanspot Popetes Pupose : Obtan tanspot coeffcent by molecula smulaton Not that the laws ae only appoxmaton that apply not-too-lage gadents In pncple tansfe coeffcents depends on c, T and v Geen-Kubo Relaton Relaton between tanspot popetes and ntegal ove tme-coelaton functon.
13 . Popetes Calculaton n MD - Tanspot Popetes Consde self-dffuson n a pue substance Consde how molecules ae dsspated when ntal confguatons ae gven as Dac delta functon Combne mass balance eqn. Wth Fck s Law Dmensonalty of gven system c(, t) D c(, t) 0 t B.C. c(, t) δ ( ) Soluton c(, t) (πdt) / d exp( ) Dt
14 . Popetes Calculaton n MD - Tanspot Popetes We do not need concentaton tself c(,t) - just dffuson coeffcent (D) < ( t) > c(, t) d c(, t) d < ( t) t > dd c(, t) D c(, t) t 0 < ( t) > ( ( t)) ddt N
15 . Popetes Calculaton n MD - Tanspot Popetes < ( t) > ( ( t)) ddt N Slope hee gves D t Plot of t vs. squae of taveled dstance gves dffuson coeffcent In 3D space, < > s mean squae dsplacement (MSD)
16 . Popetes Calculaton n MD - Tanspot Popetes An altenatve fomulaton usng velocty nstead of patcle poston t ( t) v( τ ) dτ 0 t () t v( τ ) dτ v( τ ) dτ t 0 0 dτ dτ v( τ ) v( τ ) 0 0 t τ dτ dτ v( τ ) v( τ ) 0 0 t τ dτ dτ v(0) v( τ τ ) 0 0 t t dτ dτ v(0) v( τ) 0 0 t ddt t dτ v(0) v( τ) 0 t t
17 . Popetes Calculaton n MD - Tanspot Popetes t ddt t dτ < v(0) v( τ ) > 0 t D d 0 dτ < v(0) v( τ ) > < v( 0) v( τ ) > < v( 0) v( τ ) >< v( t') v( t'' ) > Autocoelaton functon : popety dffeence between two adjacent tme steps Aea unde the cuve gves the value of self-dffuson coeffcent τ
18 . Popetes Calculaton n MD - Evaluaton of tme coelaton functons Tme consumng and eque a lot of stoage Altenatve method : FFT (Fast Foue Tansfom), Coase Ganng method A(t 0 ) A(t ) A(t ) A(t 3 ) A(t 4 ) A(t 5 ) A(t 6 ) A(t 7 ) A(t 8 ) A(t 9 ) A( Δt) A(0) A 0 A + A A + A A 3 + A 3 A 4 + A 4 A 5 + A 5 A 6 + A 6 A 7 + A 7 A 8 + A 8 A A(t 0 ) A(t ) A(t ) A(t 3 ) A(t 4 ) A(t 5 ) A(t 6 ) A(t 7 ) A(t 8 ) A(t 9 ) A( Δt) A(0) A 0 A + A A 3 + A A 4 + A 3 A 5 + A 4 A 6 + A 5 A 7 + A 6 A 8 + A 7 A A(t 0 ) A(t ) A(t ) A(t 3 ) A(t 4 ) A(t 5 ) A(t 6 ) A(t 7 ) A(t 8 ) A(t 9 ) A(5 Δ t) A(0) A 0 A 5 + A A 6 + A A 7 + A 3 A 8 + A 4 A
19 . Popetes Calculaton n MD - Tanspot Popetes Zeo-shea vscosty η VkT 0 dτ < σ xy (0) σ xy ( τ ) > σ xy x y m v v + x j f y ( j ) j Themal Conductvty d λt > q mv + u( j ) VkT dt j dτ < q(0) q( τ ) 0
20 . Popetes Calculaton n MD - Radal Dstbuton Functon Tme aveaged value of numbe densty Ensemble aveaged numbe densty g ( ) ρ ( ) ρ V g ( ) < δ( j ) N j > Just count the numbe of molecules wthn a ange
21 3. MD n Othe Ensembles - Constants Wth pope choce of g(), we can calculate useful themodynamc popetes Intenal enegy Pessue Chemcal Potental U c P 0 πnρ φ( ) g( ) d ρkt πρ 3 0 dφ( ) 3 g( ) d d 3 ρλ π dv ' μ kt ln N ' q nt 3 N V V 0 dφ( ) d g( ) 3 d
22 3. MD n Othe Ensembles Constants Hamltonan fomulaton Consevaton of knetc + potental enegy H K + U (N,V,E) ensemble Cannot be appled to othe ensemble constant T, constant P, fo example we can keep const T whle H s constant dstbuton of K and U Two types of constants Holonomc constants : may be ntegated out of equaton of moton Nonholonomc constants : non-ntegable (nvolves veloctes) Tempeatue, pessue, stess,
23 3. MD n Othe Ensembles Constants Foce momentum tempeatue to eman constant One (bad) appoach at each tme step scale momenta to foce K to desed value advance postons and momenta apply p new λp wth λ chosen to satsfy epeat equatons of moton ae evesble tanston pobabltes cannot satsfy detaled balance does not sample any well-defned ensemble
24 3. MD n Othe Ensembles Constants Gauss pncple of least constants Gaussan constants : petubatve foce ntoduced nto the equaton of moton mnmzes the devaton to classcal tajectoes of patcles fom the unpetubed tajectoes Consde a functon f, a functon of patcle acceleaton f ( ) m F m f0 : nomal Newtonan equaton of moton othewse, constaned non-newtonan equaton of moton Gauss pncple : physcal acceleaton f to be mnmum ( f ( ) ζ g( )) 0 ζ : Lagangn (Gauss) Multple
25 3. MD n Othe Ensembles Constants Constant Tempeatue constants 0 ),, ( m dt dg t g 0 3 ), ( NkT m t G ( ) 0 ) ( ) ( g f ζ 0 j j j j j j j m m m F m m F ζ Newtonan Constant foce
26 3. MD n Othe Ensembles Constants Modfed equaton of moton m m F p p ζ m F ζ one of good appoach, but tempeatue s not specfed!
27 3. MD n Othe Ensembles Nose Themostat Extended Lagangan Equaton of Moton N m( s ) N Q LNose U( ) + s gkt lns p p L ms s L s Qs s Ks Qs U gktln s
28 3. MD n Othe Ensembles Nose-Hoove Themostat Equatons of moton p m p F ξp s ξ s N p ξ gkt Q m Integaton schemes pedcto-coecto algothm s staghtfowad Velet algothm s feasble, but tcky to mplement v F t-δt t t+δt At ths step, update of ξ depends on p; update of p depends on ξ p F ξp N p ξ Q m gkt
Thermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering
Themodynamcs of solds 4. Statstcal themodynamcs and the 3 d law Kwangheon Pak Kyung Hee Unvesty Depatment of Nuclea Engneeng 4.1. Intoducton to statstcal themodynamcs Classcal themodynamcs Statstcal themodynamcs
More informationPHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle
1 PHYS 705: Classcal Mechancs Devaton of Lagange Equatons fom D Alembet s Pncple 2 D Alembet s Pncple Followng a smla agument fo the vtual dsplacement to be consstent wth constants,.e, (no vtual wok fo
More informationMechanics Physics 151
Mechancs Physcs 151 Lectue 18 Hamltonan Equatons of Moton (Chapte 8) What s Ahead We ae statng Hamltonan fomalsm Hamltonan equaton Today and 11/6 Canoncal tansfomaton 1/3, 1/5, 1/10 Close lnk to non-elatvstc
More informationEnergy in Closed Systems
Enegy n Closed Systems Anamta Palt palt.anamta@gmal.com Abstact The wtng ndcates a beakdown of the classcal laws. We consde consevaton of enegy wth a many body system n elaton to the nvese squae law and
More informationV. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. "Flux": = da i. "Force": = -Â g a ik k = X i. Â J i X i (7.
Themodynamcs and Knetcs of Solds 71 V. Pncples of Ievesble Themodynamcs 5. Onsage s Teatment s = S - S 0 = s( a 1, a 2,...) a n = A g - A n (7.6) Equlbum themodynamcs detemnes the paametes of an equlbum
More informationSet of square-integrable function 2 L : function space F
Set of squae-ntegable functon L : functon space F Motvaton: In ou pevous dscussons we have seen that fo fee patcles wave equatons (Helmholt o Schödnge) can be expessed n tems of egenvalue equatons. H E,
More informationScalars and Vectors Scalar
Scalas and ectos Scala A phscal quantt that s completel chaacteed b a eal numbe (o b ts numecal value) s called a scala. In othe wods a scala possesses onl a magntude. Mass denst volume tempeatue tme eneg
More informationPhysics 11b Lecture #2. Electric Field Electric Flux Gauss s Law
Physcs 11b Lectue # Electc Feld Electc Flux Gauss s Law What We Dd Last Tme Electc chage = How object esponds to electc foce Comes n postve and negatve flavos Conseved Electc foce Coulomb s Law F Same
More informationMultistage Median Ranked Set Sampling for Estimating the Population Median
Jounal of Mathematcs and Statstcs 3 (: 58-64 007 ISSN 549-3644 007 Scence Publcatons Multstage Medan Ranked Set Samplng fo Estmatng the Populaton Medan Abdul Azz Jeman Ame Al-Oma and Kamaulzaman Ibahm
More informationChapter 8. Linear Momentum, Impulse, and Collisions
Chapte 8 Lnea oentu, Ipulse, and Collsons 8. Lnea oentu and Ipulse The lnea oentu p of a patcle of ass ovng wth velocty v s defned as: p " v ote that p s a vecto that ponts n the sae decton as the velocty
More information24-2: Electric Potential Energy. 24-1: What is physics
D. Iyad SAADEDDIN Chapte 4: Electc Potental Electc potental Enegy and Electc potental Calculatng the E-potental fom E-feld fo dffeent chage dstbutons Calculatng the E-feld fom E-potental Potental of a
More informationStellar Astrophysics. dt dr. GM r. The current model for treating convection in stellar interiors is called mixing length theory:
Stella Astophyscs Ovevew of last lectue: We connected the mean molecula weght to the mass factons X, Y and Z: 1 1 1 = X + Y + μ 1 4 n 1 (1 + 1) = X μ 1 1 A n Z (1 + ) + Y + 4 1+ z A Z We ntoduced the pessue
More informationChapter 23: Electric Potential
Chapte 23: Electc Potental Electc Potental Enegy It tuns out (won t show ths) that the tostatc foce, qq 1 2 F ˆ = k, s consevatve. 2 Recall, fo any consevatve foce, t s always possble to wte the wok done
More informationPHY126 Summer Session I, 2008
PHY6 Summe Sesson I, 8 Most of nfomaton s avalable at: http://nngoup.phscs.sunsb.edu/~chak/phy6-8 ncludng the sllabus and lectue sldes. Read sllabus and watch fo mpotant announcements. Homewok assgnment
More informationLarge scale magnetic field generation by accelerated particles in galactic medium
Lage scale magnetc feld geneaton by acceleated patcles n galactc medum I.N.Toptygn Sant Petesbug State Polytechncal Unvesty, depatment of Theoetcal Physcs, Sant Petesbug, Russa 2.Reason explonatons The
More informationContact, information, consultations
ontact, nfomaton, consultatons hemsty A Bldg; oom 07 phone: 058-347-769 cellula: 664 66 97 E-mal: wojtek_c@pg.gda.pl Offce hous: Fday, 9-0 a.m. A quote of the week (o camel of the week): hee s no expedence
More informationPHYS Week 5. Reading Journals today from tables. WebAssign due Wed nite
PHYS 015 -- Week 5 Readng Jounals today fom tables WebAssgn due Wed nte Fo exclusve use n PHYS 015. Not fo e-dstbuton. Some mateals Copyght Unvesty of Coloado, Cengage,, Peason J. Maps. Fundamental Tools
More informationPart V: Velocity and Acceleration Analysis of Mechanisms
Pat V: Velocty an Acceleaton Analyss of Mechansms Ths secton wll evew the most common an cuently pactce methos fo completng the knematcs analyss of mechansms; escbng moton though velocty an acceleaton.
More informationDynamics of Rigid Bodies
Dynamcs of Rgd Bodes A gd body s one n whch the dstances between consttuent patcles s constant thoughout the moton of the body,.e. t keeps ts shape. Thee ae two knds of gd body moton: 1. Tanslatonal Rectlnea
More informationTian Zheng Department of Statistics Columbia University
Haplotype Tansmsson Assocaton (HTA) An "Impotance" Measue fo Selectng Genetc Makes Tan Zheng Depatment of Statstcs Columba Unvesty Ths s a jont wok wth Pofesso Shaw-Hwa Lo n the Depatment of Statstcs at
More informationReview. Physics 231 fall 2007
Reew Physcs 3 all 7 Man ssues Knematcs - moton wth constant acceleaton D moton, D pojectle moton, otatonal moton Dynamcs (oces) Enegy (knetc and potental) (tanslatonal o otatonal moton when detals ae not
More informationLearning the structure of Bayesian belief networks
Lectue 17 Leanng the stuctue of Bayesan belef netwoks Mlos Hauskecht mlos@cs.ptt.edu 5329 Sennott Squae Leanng of BBN Leanng. Leanng of paametes of condtonal pobabltes Leanng of the netwok stuctue Vaables:
More informationE For K > 0. s s s s Fall Physical Chemistry (II) by M. Lim. singlet. triplet
Eneges of He electonc ψ E Fo K > 0 ψ = snglet ( )( ) s s+ ss αβ E βα snglet = ε + ε + J s + Ks Etplet = ε + ε + J s Ks αα ψ tplet = ( s s ss ) ββ ( αβ + βα ) s s s s s s s s ψ G = ss( αβ βα ) E = ε + ε
More informationNotes on Analytical Dynamics
Notes on Analytcal Dynamcs Jan Peters & Mchael Mstry October 7, 004 Newtonan Mechancs Basc Asssumptons and Newtons Laws Lonely pontmasses wth postve mass Newtons st: Constant velocty v n an nertal frame
More informationPhysics 207 Lecture 16
Physcs 07 Lectue 6 Goals: Lectue 6 Chapte Extend the patcle odel to gd-bodes Undestand the equlbu of an extended object. Analyze ollng oton Undestand otaton about a fxed axs. Eploy consevaton of angula
More informationALL QUESTIONS ARE WORTH 20 POINTS. WORK OUT FIVE PROBLEMS.
GNRAL PHYSICS PH -3A (D. S. Mov) Test (/3/) key STUDNT NAM: STUDNT d #: -------------------------------------------------------------------------------------------------------------------------------------------
More informationIntegral Vector Operations and Related Theorems Applications in Mechanics and E&M
Dola Bagayoko (0) Integal Vecto Opeatons and elated Theoems Applcatons n Mechancs and E&M Ι Basc Defnton Please efe to you calculus evewed below. Ι, ΙΙ, andιιι notes and textbooks fo detals on the concepts
More informationPhysics 2A Chapter 11 - Universal Gravitation Fall 2017
Physcs A Chapte - Unvesal Gavtaton Fall 07 hese notes ae ve pages. A quck summay: he text boxes n the notes contan the esults that wll compse the toolbox o Chapte. hee ae thee sectons: the law o gavtaton,
More informationOne-dimensional kinematics
Phscs 45 Fomula Sheet Eam 3 One-dmensonal knematcs Vectos dsplacement: Δ total dstance taveled aveage speed total tme Δ aveage veloct: vav t t Δ nstantaneous veloct: v lm Δ t v aveage acceleaton: aav t
More informationReview of Vector Algebra and Vector Calculus Operations
Revew of Vecto Algeba and Vecto Calculus Opeatons Tpes of vaables n Flud Mechancs Repesentaton of vectos Dffeent coodnate sstems Base vecto elatons Scala and vecto poducts Stess Newton s law of vscost
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationLecture 14: Forces and Stresses
The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal
More informationCOMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS
ultscence - XXX. mcocd Intenatonal ultdscplnay Scentfc Confeence Unvesty of skolc Hungay - pl 06 ISBN 978-963-358-3- COPLEENTRY ENERGY ETHOD FOR CURVED COPOSITE BES Ákos József Lengyel István Ecsed ssstant
More informationMolecular Dynamics and Monte Carlo Methods
May 8, 1 Molecula Modelng and Smulaton Molecula Dynamcs and Monte Calo Methods Agcultual Bonfomatcs Reseach Unt, Gaduate School of Agcultual and Lfe Scences, The Unvesty of Tokyo Tohu Teada 1 Contents
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More information4.4 Continuum Thermomechanics
4.4 Contnuum Themomechancs The classcal themodynamcs s now extended to the themomechancs of a contnuum. The state aables ae allowed to ay thoughout a mateal and pocesses ae allowed to be eesble and moe
More informationSome Approximate Analytical Steady-State Solutions for Cylindrical Fin
Some Appoxmate Analytcal Steady-State Solutons fo Cylndcal Fn ANITA BRUVERE ANDRIS BUIIS Insttute of Mathematcs and Compute Scence Unvesty of Latva Rana ulv 9 Rga LV459 LATVIA Astact: - In ths pape we
More informationCSU ATS601 Fall Other reading: Vallis 2.1, 2.2; Marshall and Plumb Ch. 6; Holton Ch. 2; Schubert Ch r or v i = v r + r (3.
3 Eath s Rotaton 3.1 Rotatng Famewok Othe eadng: Valls 2.1, 2.2; Mashall and Plumb Ch. 6; Holton Ch. 2; Schubet Ch. 3 Consde the poston vecto (the same as C n the fgue above) otatng at angula velocty.
More informationDYNAMICS VECTOR MECHANICS FOR ENGINEERS: Kinematics of Rigid Bodies in Three Dimensions. Seventh Edition CHAPTER
Edton CAPTER 8 VECTOR MECANCS FOR ENGNEERS: DYNAMCS Fednand P. Bee E. Russell Johnston, J. Lectue Notes: J. Walt Ole Teas Tech Unvest Knematcs of Rgd Bodes n Thee Dmensons 003 The McGaw-ll Companes, nc.
More informationNVU dynamics. III. Simulating molecules at constant potential energy
THE JOURNAL OF CHEMICAL PHYSICS 137, 244101 (2012) NVU dynamcs. III. Smulatng molecules at constant potental enegy Tond S. Ingebgtsen a) and Jeppe C. Dye DNRF Cente Glass and Tme, IMFUFA, Depatment of
More informationMOLECULAR DYNAMICS ,..., What is it? 2 = i i
MOLECULAR DYNAMICS What s t? d d x t 2 m 2 = F ( x 1,..., x N ) =1,,N r ( x1 ( t),..., x ( t)) = v = ( x& 1 ( t ),..., x& ( t )) N N What are some uses of molecular smulatons and modelng? Conformatonal
More informationSlide 1. Quantum Mechanics: the Practice
Slde Quantum Mecancs: te Pactce Slde Remnde: Electons As Waves Wavelengt momentum = Planck? λ p = = 6.6 x 0-34 J s Te wave s an exctaton a vbaton: We need to know te ampltude of te exctaton at evey pont
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More information2/24/2014. The point mass. Impulse for a single collision The impulse of a force is a vector. The Center of Mass. System of particles
/4/04 Chapte 7 Lnea oentu Lnea oentu of a Sngle Patcle Lnea oentu: p υ It s a easue of the patcle s oton It s a vecto, sla to the veloct p υ p υ p υ z z p It also depends on the ass of the object, sla
More informationTest 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o?
Test 1 phy 0 1. a) What s the pupose of measuement? b) Wte all fou condtons, whch must be satsfed by a scala poduct. (Use dffeent symbols to dstngush opeatons on ectos fom opeatons on numbes.) c) What
More informationCapítulo. Three Dimensions
Capítulo Knematcs of Rgd Bodes n Thee Dmensons Mecánca Contents ntoducton Rgd Bod Angula Momentum n Thee Dmensons Pncple of mpulse and Momentum Knetc Eneg Sample Poblem 8. Sample Poblem 8. Moton of a Rgd
More informationEE 5337 Computational Electromagnetics (CEM)
7//28 Instucto D. Raymond Rumpf (95) 747 6958 cumpf@utep.edu EE 5337 Computatonal Electomagnetcs (CEM) Lectue #6 TMM Extas Lectue 6These notes may contan copyghted mateal obtaned unde fa use ules. Dstbuton
More informationPhysics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on
More informationCOLLEGE OF FOUNDATION AND GENERAL STUDIES PUTRAJAYA CAMPUS FINAL EXAMINATION TRIMESTER /2017
COLLEGE OF FOUNDATION AND GENERAL STUDIES PUTRAJAYA CAMPUS FINAL EXAMINATION TRIMESTER 1 016/017 PROGRAMME SUBJECT CODE : Foundaton n Engneeng : PHYF115 SUBJECT : Phscs 1 DATE : Septembe 016 DURATION :
More informationDensity Functional Theory I
Densty Functonal Theoy I cholas M. Hason Depatment of Chemsty Impeal College Lonon & Computatonal Mateals Scence Daesbuy Laboatoy ncholas.hason@c.ac.uk Densty Functonal Theoy I The Many Electon Schönge
More informationCFD Investigations of Spatial Arc Kinetic Influence on Fuel Burning- Out in the Tornado Combustor
CFD Investgatons of Spatal Ac Knetc Influence on Fuel Bunng- Out n the Tonado Combusto Igo Matveev, Appled Plasma Technology, U.S.A.,., Sehy Sebn and Anna Mostpaneno Natonal Unvesty of Shpbuldng, Uane
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationChapter I Matrices, Vectors, & Vector Calculus 1-1, 1-9, 1-10, 1-11, 1-17, 1-18, 1-25, 1-27, 1-36, 1-37, 1-41.
Chapte I Matces, Vectos, & Vecto Calculus -, -9, -0, -, -7, -8, -5, -7, -36, -37, -4. . Concept of a Scala Consde the aa of patcles shown n the fgue. he mass of the patcle at (,) can be epessed as. M (,
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationA Brief Guide to Recognizing and Coping With Failures of the Classical Regression Assumptions
A Bef Gude to Recognzng and Copng Wth Falues of the Classcal Regesson Assumptons Model: Y 1 k X 1 X fxed n epeated samples IID 0, I. Specfcaton Poblems A. Unnecessay explanatoy vaables 1. OLS s no longe
More informationLINEAR MOMENTUM. product of the mass m and the velocity v r of an object r r
LINEAR MOMENTUM Imagne beng on a skateboad, at est that can move wthout cton on a smooth suace You catch a heavy, slow-movng ball that has been thown to you you begn to move Altenatvely you catch a lght,
More informationEvent Shape Update. T. Doyle S. Hanlon I. Skillicorn. A. Everett A. Savin. Event Shapes, A. Everett, U. Wisconsin ZEUS Meeting, October 15,
Event Shape Update A. Eveett A. Savn T. Doyle S. Hanlon I. Skllcon Event Shapes, A. Eveett, U. Wsconsn ZEUS Meetng, Octobe 15, 2003-1 Outlne Pogess of Event Shapes n DIS Smla to publshed pape: Powe Coecton
More informationRotary motion
ectue 8 RTARY TN F THE RGD BDY Notes: ectue 8 - Rgd bod Rgd bod: j const numbe of degees of feedom 6 3 tanslatonal + 3 ota motons m j m j Constants educe numbe of degees of feedom non-fee object: 6-p
More informationElectron density: Properties of electron density (non-negative): => exchange-correlation functionals should respect these conditions.
lecton densty: ρ ( =... Ψ(,,..., ds d... d Pobablty of fndng one electon of abtay spn wthn a volume element d (othe electons may be anywhee. s Popetes of electon densty (non-negatve:.. 3. ρ ( d = ρ( =
More information1. Starting with the local version of the first law of thermodynamics q. derive the statement of the first law of thermodynamics for a control volume
EN10: Contnuum Mechancs Homewok 5: Alcaton of contnuum mechancs to fluds Due 1:00 noon Fda Febua 4th chool of Engneeng Bown Unvest 1. tatng wth the local veson of the fst law of themodnamcs q jdj q t and
More informationTHE MANY-BODY PROBLEM
3.30: Lectue 5 Feb 5 005 THE MANY-BODY PROBLEM Feb 5 005 3.30 Atomstc Modelng of Mateals -- Geband Cede and Ncola Maza When s a patcle lke a wave? Wavelength momentum = Planck λ p = h h = 6.6 x 0-34 J
More informationChapter IV Vector and Tensor Analysis IV.2 Vector and Tensor Analysis September 29,
hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe 9, 08 47 hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe 9, 08 48 I. ETOR AND TENOR ANALYI I... Tenso functon th Let A
More informationThe Unique Solution of Stochastic Differential Equations With. Independent Coefficients. Dietrich Ryter.
The Unque Soluton of Stochastc Dffeental Equatons Wth Independent Coeffcents Detch Ryte RyteDM@gawnet.ch Mdatweg 3 CH-4500 Solothun Swtzeland Phone +4132 621 13 07 SDE s must be solved n the ant-itô sense
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationSupersymmetry in Disorder and Chaos (Random matrices, physics of compound nuclei, mathematics of random processes)
Supesymmety n Dsoe an Chaos Ranom matces physcs of compoun nucle mathematcs of anom pocesses Lteatue: K.B. Efetov Supesymmety n Dsoe an Chaos Cambge Unvesty Pess 997999 Supesymmety an Tace Fomulae I.V.
More informationInstantaneous velocity field of a round jet
Fee shea flows Instantaneos velocty feld of a ond et 3 Aveage velocty feld of a ond et 4 Vtal ogn nozzle coe Developng egon elf smla egon 5 elf smlaty caled vaables: ~ Q ξ ( ξ, ) y δ ( ) Q Q (, y) ( )
More informationMachine Learning. Spectral Clustering. Lecture 23, April 14, Reading: Eric Xing 1
Machne Leanng -7/5 7/5-78, 78, Spng 8 Spectal Clusteng Ec Xng Lectue 3, pl 4, 8 Readng: Ec Xng Data Clusteng wo dffeent ctea Compactness, e.g., k-means, mxtue models Connectvty, e.g., spectal clusteng
More informationNuclear Chart. Takashi NAKATSUKASA Theoretical Nuclear Physics Laboratory RIKEN Nishina Center. Real-space, real-time approaches ) Few-body model
Takash NAKATSUKASA Theoetcal Nuclea Physcs Laboatoy RIKEN Nshna Cente 009.3.5-6 Mn-WS: Real-space, eal-tme appoaches DFT, TDDFT ((Q)RPA ) Few-body model (CDCC ) Nuclea Chat 70 Los Alamos Natonal Laboatoy's
More informationUNIT10 PLANE OF REGRESSION
UIT0 PLAE OF REGRESSIO Plane of Regesson Stuctue 0. Intoducton Ojectves 0. Yule s otaton 0. Plane of Regesson fo thee Vaales 0.4 Popetes of Resduals 0.5 Vaance of the Resduals 0.6 Summay 0.7 Solutons /
More information8 Baire Category Theorem and Uniform Boundedness
8 Bae Categoy Theoem and Unfom Boundedness Pncple 8.1 Bae s Categoy Theoem Valdty of many esults n analyss depends on the completeness popety. Ths popety addesses the nadequacy of the system of atonal
More information3. A Review of Some Existing AW (BT, CT) Algorithms
3. A Revew of Some Exstng AW (BT, CT) Algothms In ths secton, some typcal ant-wndp algothms wll be descbed. As the soltons fo bmpless and condtoned tansfe ae smla to those fo ant-wndp, the pesented algothms
More informationComplex atoms and the Periodic System of the elements
Complex atoms and the Peodc System of the elements Non-cental foces due to electon epulson Cental feld appoxmaton electonc obtals lft degeneacy of l E n l = R( hc) ( n δ ) l Aufbau pncple Lectue Notes
More informationCS649 Sensor Networks IP Track Lecture 3: Target/Source Localization in Sensor Networks
C649 enso etwoks IP Tack Lectue 3: Taget/ouce Localaton n enso etwoks I-Jeng Wang http://hng.cs.jhu.edu/wsn06/ png 006 C 649 Taget/ouce Localaton n Weless enso etwoks Basc Poblem tatement: Collaboatve
More informationCE 530 Molecular Simulation
CE 530 Molecular Smulaton Lecture 3 Molecular Dynamcs n Other Ensembles Davd A. Kofke Department of Chemcal Engneerng SUY Buffalo kofke@eng.buffalo.edu Revew Molecular dynamcs s a numercal ntegraton of
More informationFundamental principles
JU 07/HL Dnacs and contol of echancal sstes Date Da (0/08) Da (03/08) Da 3 (05/08) Da 4 (07/08) Da 5 (09/08) Da 6 (/08) Content Reve of the bascs of echancs. Kneatcs of gd bodes coodnate tansfoaton, angula
More information1. A body will remain in a state of rest, or of uniform motion in a straight line unless it
Pncples of Dnamcs: Newton's Laws of moton. : Foce Analss 1. A bod wll eman n a state of est, o of unfom moton n a staght lne unless t s acted b etenal foces to change ts state.. The ate of change of momentum
More informationPhysics 202, Lecture 2. Announcements
Physcs 202, Lectue 2 Today s Topcs Announcements Electc Felds Moe on the Electc Foce (Coulomb s Law The Electc Feld Moton of Chaged Patcles n an Electc Feld Announcements Homewok Assgnment #1: WebAssgn
More informationA. Thicknesses and Densities
10 Lab0 The Eath s Shells A. Thcknesses and Denstes Any theoy of the nteo of the Eath must be consstent wth the fact that ts aggegate densty s 5.5 g/cm (ecall we calculated ths densty last tme). In othe
More informationCorrespondence Analysis & Related Methods
Coespondence Analyss & Related Methods Ineta contbutons n weghted PCA PCA s a method of data vsualzaton whch epesents the tue postons of ponts n a map whch comes closest to all the ponts, closest n sense
More informationP 365. r r r )...(1 365
SCIENCE WORLD JOURNAL VOL (NO4) 008 www.scecncewoldounal.og ISSN 597-64 SHORT COMMUNICATION ANALYSING THE APPROXIMATION MODEL TO BIRTHDAY PROBLEM *CHOJI, D.N. & DEME, A.C. Depatment of Mathematcs Unvesty
More informationMinimization of the Free Energy under a Given Pressure by Natural Iteration Method
Mateals ansactons, Vol. 5, No. 3 (011) pp. 48 to 43 #011 he Japan Insttute of Metals Mnmzaton of the Fee Enegy unde a Gven Pessue by Natual Iteaton Method Naoya Kyokane 1; * and etsuo Moh ;3 1 Dvson of
More informationUser-friendly model of heat transfer in. preheating, cool down and casting
ANNUAL REPORT 2010 UIUC, August 12, 2010 Use-fendly model of heat tansfe n submeged enty nozzles dung peheatng, cool down and castng Vaun Kuma Sngh, B.G. Thomas Depatment of Mechancal Scence and Engneeng
More informationgravity r2,1 r2 r1 by m 2,1
Gavtaton Many of the foundatons of classcal echancs wee fst dscoveed when phlosophes (ealy scentsts and atheatcans) ted to explan the oton of planets and stas. Newton s ost faous fo unfyng the oton of
More informationPhysical & Interfacial Electrochemistry 2013
Physcal & Intefacal Electochemsty 013 Lectue 3. Ion-on nteactons n electolyte solutons. Module JS CH3304 MoleculaThemodynamcs and Knetcs Ion-Ion Inteactons The themodynamc popetes of electolyte solutons
More informationLASER ABLATION ICP-MS: DATA REDUCTION
Lee, C-T A Lase Ablaton Data educton 2006 LASE ABLATON CP-MS: DATA EDUCTON Cn-Ty A. Lee 24 Septembe 2006 Analyss and calculaton of concentatons Lase ablaton analyses ae done n tme-esolved mode. A ~30 s
More informationComplementi di Fisica Lecture 19
Complement d Fsca Lectue 9 Lvo Lance Unvestà d Teste Teste, 0--04 Couse Outlne - Remnde The phscs of semconducto devces: an ntoducton Quantum Mechancs: an ntoducton just a few moe comments on: (3-d) Hdogen
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationRotating Disk Electrode -a hydrodynamic method
Rotatng Dsk Electode -a hdodnamc method Fe Lu Ma 3, 0 ente fo Electochemcal Engneeng Reseach Depatment of hemcal and Bomolecula Engneeng Rotatng Dsk Electode A otatng dsk electode RDE s a hdodnamc wokng
More information19 The Born-Oppenheimer Approximation
9 The Bon-Oppenheme Appoxmaton The full nonelatvstc Hamltonan fo a molecule s gven by (n a.u.) Ĥ = A M A A A, Z A + A + >j j (883) Lets ewte the Hamltonan to emphasze the goal as Ĥ = + A A A, >j j M A
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationPhysics Exam II Chapters 25-29
Physcs 114 1 Exam II Chaptes 5-9 Answe 8 of the followng 9 questons o poblems. Each one s weghted equally. Clealy mak on you blue book whch numbe you do not want gaded. If you ae not sue whch one you do
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationChapter Fifiteen. Surfaces Revisited
Chapte Ffteen ufaces Revsted 15.1 Vecto Descpton of ufaces We look now at the vey specal case of functons : D R 3, whee D R s a nce subset of the plane. We suppose s a nce functon. As the pont ( s, t)
More informationMolecular Dynamic Simulations of Nickel Nanowires at Various Temperatures
Intenatonal Jounal of Scentfc and Innovatve Mathematcal Reseach (IJSIMR Volume 2, Issue 3, Mach 204, PP 30-305 ISS 2347-307X (Pnt & ISS 2347-342 (Onlne www.acounals.og Molecula Dynamc Smulatons of ckel
More informationMultipole Radiation. March 17, 2014
Multpole Radaton Mach 7, 04 Zones We wll see that the poblem of hamonc adaton dvdes nto thee appoxmate egons, dependng on the elatve magntudes of the dstance of the obsevaton pont,, and the wavelength,
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationMechanics Physics 151
Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and
More informationChapter IV Vector and Tensor Analysis IV.2 Vector and Tensor Analysis September 23,
hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe, 07 47 hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe, 07 48 I. ETOR AND TENOR ANALYI I... Tenso functon th Let A n n
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More information