DFT and Molecular Dynamics
|
|
- Kelly Ford
- 6 years ago
- Views:
Transcription
1 Chapter 6 DFT and Molecular Dynamcs Let us stck to our prevous assumpton that ons behave as classcal partcles. Ths turns out to be a very good approxmaton, except for the lghtest atoms and for vbratonal modes around the equlbrum postons at low temperature. Also, let us assume the valdty of the Born-Oppenhemer (BO approxmaton and that electrons are always on the BO surface, that s, n the ground state correspondng to ther nstantaneous postons. Under these assumptons the dynamcal behavor of ons can be descrbed by a classcal Lagrangan L = 1 M Ṙ 2 E tot ({R} (6.1 2 where M are the mass of ons. The correspondng equatons of moton: d L dt Ṙ L R = 0, P = L Ṙ (6.2 are nothng but Newton s equatons. It s temptng to use Eq.6.1 as the bass for a molecular dynamcs (MD study. In classcal MD, the forces are generated by an nteratomc potental (often a sum of twobody terms lke Lennard-Jones potentals and the Newton equatons are dscretzed and numercally solved. The dscrete nterval of tme s called tme step. A sequence of atomc coordnates and veloctes s generated startng from a sutable ntal set of coordnates and veloctes. The sequence can be used to calculate thermodynamcal averages. Straghtforward MD wll sample the mcrocanoncal ensemble: constant energy at fxed volume, but t s possble to buld a dynamcs at constant temperature (canoncal ensemble usng a thermostat that smulates a thermal bath, or at constant pressure, by addng a fcttous dynamcs on the volume, and even more complex cases. MD can also be used to fnd the global mnma usng the smulated annealng technque. The confguraton space s sampled at equlbrum, then the knetc energy s gradually removed from the systems that has the possblty (but s not guaranteed to do so to reach the global mnmum. Such procedure s sometmes the only practcal way to fnd the global mnmum for especally hard problems. In mathematcal terms, easy problems are exactly solvable by computer algorthms n polynomal tme, that s, n a number of steps that s a polynomal functon of the dmenson of the problem; hard problems are solved n exponental tme. A problem s NP (nondetermnstc polynomal f ts soluton (f one exsts can be guessed and verfed n polynomal tme. Ths s the knd of problems 34
2 for whch the smulated annealng has been devsed. The determnaton of the structure n clusters s beleved to be a NP-hard problem. 6.1 Classcal Molecular Dynamcs Let us consder the most basc MD : a purely mechancal system of N atoms, enclosed n a volume V (usually wth perodcal boundary condtons, PBC, for a condensed-matter system, havng mechancal energy E = T + E p, where T = 1 2 MṘ2 s the knetc energy of ons, E p = E p ({R} s the nteratomc potental energy. Ths s known as the N V E, or mcrocanoncal, ensemble Dscretzaton of the equaton of moton The numercal soluton (ntegraton of the equatons of motons s generally performed usng the Verlet algorthm. Ths s obtaned from the followng basc and very smple equatons : R (t + δt = R (t + δtv (t + δt2 2M f (t + δt3 6 b (t + O(δt 4 (6.3 R (t δt = R (t δtv (t + δt2 2M f (t δt3 6 b (t + O(δt 4 (6.4 where V = Ṙ are veloctes, f forces actng on on. By summng and subtractng Eqs. (6.3 and (6.4 we get the Verlet algorthm: R (t + δt = 2R (t R (t δt + δt2 M f (t + O(δt 4 (6.5 V (t = 1 2δt [R (t + δt R (t δt] + O(δt 3. (6.6 The veloctes are one step behnd the postons, but they are not used to update the postons. It s possble to recast the Verlet algorthm nto an equvalent form (one gvng exactly the same trajectores n whch both veloctes and postons are updated n the same step. By combnng Eq.(6.3 wth Eq.(6.4 dsplaced n tme by +δt, one fnds V (t + δt = V (t + δt 2M [f (t + f (t + δt] (6.7 R (t + δt = R (t + δtv (t + δt2 2M f (t. (6.8 Note that the update of veloctes requres the forces for the new postons. Ths algorthm s known as Velocty Verlet. Its equvalence to the Verlet algorthm may not seem evdent, but t can be proved qute smply. In spte of hs smplcty, the Verlet algorthm, n any ncarnaton, s effcent and numercally stable. In partcular, t yelds trajectores that conserve to a very good degree of accuracy the energy E. A small loss of energy conservaton, due both to numercal errors and to the dscretzaton, s unavodable, but a systematc drft of the energy s not acceptable. In ths respect Verlet s superor to apparently better (.e. hgher-order schemes. In one of the followng sectons we wll see one reason why ths happen. 35
3 6.1.2 Thermodynamcal averages In the followng we wll use the phase space canoncal varables, collectvely ndcated as R, P, nstead of coordnates and veloctes. Let us consder an observable quantty A = A(R, P. From a practcal pont of vew, the calculaton of thermodynamcal averages n classcal MD s an average over many tme steps: A T = 1 T T A(R(t, P(tdt 1 M A(t n, t n = nδt, t M = Mδt = T. (6.9 0 M n=1 For an ergodc system (that s, one whose trajectores n a suffcently long tme pass arbtrary close to any pont n the phase space, t s beleved that: lm T A T A (6.10 where s the average over the correspondng ensemble: A = ρ(r, PA(R, PdRdP (6.11 where ρ s the probablty of a mcroscopc state. In NV E MD the mcrocanoncal ensemble s sampled: ρ NV E (R, P = g(n δ(h E (6.12 Ω where H s the Hamltonan correspondng to the Lagrangan of Eq.(6.1, E s the mechancal energy (ncludng knetc energy of ons g(n = (h 3N N! 1 for N ndstngushable atoms, and Ω, related to the entropy S by the Boltzmann relaton S = k B log Ω, s the total number of mcroscopc states: Ω = g(n drdpδ(h E. (6.13 The tme step must be as bg as possble n order to sample as much phase space as possble, but at the same tme t must be small enough to allow to follow the moton the ons wth lttle loss of accuracy (whch typcally shows up under the form of non-conservaton of the energy. Typcally δt δt max, where δt max s the perod of the fastest phonon mode: δt max = 1/ω max Verlet algorthm as untary dscretzaton of the Louvllan The tme evoluton of an observable A = A(R, P, t can be wrtten as da dt = ( A A Ṙ + R Ṗ + A = LA + A P t t (6.14 where the operator L s called the Louvllan and s related to the Hamltonan va: LA = [H, A] (6.15 where [...] stands for Posson s brakets. Assumng that A = A(R, P does not depend explctly on the tme, the Louvllan entrely determnes the tme evoluton of A: formally, A(t = e Lt A(t = 0 = U(tA(t = 0. (
4 It can be shown that L s an Hermtan operator and thus U s a untary operator,.e., U = U 1 (unsurprsngly, snce U (t = U( t and tme-reversal symmetry must hold. We can wrte L as L = ( Ṙ + R Ṗ = ( Ṙ + f (6.17 P R P and fnally as a sum of two terms, one actng on coordnates and one on momenta: L = L p + L r, where L p = f, P L r = Ṙ. (6.18 R Untl now, we have just recast the classcal equatons of moton nto a more elegant but not especally useful formalsm. Let us now dscretze the tme evoluton operator, by dvdng t nto N small ntervals δt = t/n, and apply the Trotter approxmaton: e (Lp+Lrt = [e (Lp+Lrδt] N [ N = e Lpδt/2 e Lrδt e Lpδt/2 + O(δt ] 3. (6.19 Remember that L p and L r are operators: the Trotter approxmaton s not trval. Let us apply the operator between square brackets to a pont (R (t, P (t n phase space at tme t. We wll use the known result e a / x f(x = f(x + a (6.20 f a does not depend on x. Snce L p and L r are sums of terms actng on each partcle separately, we can consder ther acton on each partcle ndependently. ( e Lpδt/2 (R, P = R, P + δt 2 f (R (R, P (6.21 e Lrδt ( ( R, P = R + δt P M, P (6.22 ( = R + δt P + δt2 M 2m f (R, P + δt 2 f (R (R, P (6.23 e Lpδt/2 ( ( R, P = R, P + δt 2 f (R (6.24 = ( R + δt M P + δt2 2M f (R, P + δt 2 [f (R + f (R ] (6.25 Notng that f (R = f (t, f (R = f (t+δt, the last expresson s nothng but the velocty Verlet algorthm for (R (t + δt, P (t + δt. In concluson: the Verlet algorthm may be derved by a dscretzaton of the tme evoluton operator that conserves untarty. Such property s crucal for any well-behaved ntegraton algorthm one can thnk of Canoncal ensemble n MD We are often nterested n systems n thermal equlbrum wth a thermal bath at temperature T : the NV T or canoncal ensemble, for whch ρ NV T (R, P = g(n Z e H(R,P/k BT (
5 where Z s the partton functon: Z = g(n drdpe H(R,P/k BT. (6.27 Integraton over P gves for the partton functon of N dentcal atoms: where λ s the thermal wavelength: Z = Z r /(N!λ 3N (6.28 λ = h 2πMkb T (6.29 and Z r s the confguratonal partton functon: Z r = dr 1... dr n e Ep(R/kBT. (6.30 In the canoncal ensemble, the temperature s related to the expectaton value of the knetc energy: N P 2 NV T = 3 2M 2 Nk BT. (6.31 =1 Ths suggests a rather straghtforward algorthm (the velocty rescalng to keep the system at the desred temperature: rescale the veloctes every tme the knetc energy departs from the target value by more than some threshold. Whle effectve to brng the system to the desred temperature, the velocty rescalng does not guarantee that the averages wll be correct, although some varants of t do qute a good job even n ths respect. The correct way to sample the canoncal ensemble s to ntroduce a sutable thermostat that smulates the contact wth a thermal bath at temperature T. A well-known one s the Nosé-Hoover thermostat. Ths ntroduces an addtonal fcttous degree of freedom, producng a dynamcal frcton force havng the effect of heatng ons when the knetc energy s lower than the desred value, coolng them n the opposte case. Specfcally, the equatons of moton become R = f M ζṙ (6.32 ζ = 1 [ N ] M Ṙ 2 3Nk b T (6.33 Q =1 where Q plays the role of thermal mass. The constant of moton for ths system s H = H + Q 2 ζ 2 + 3Nk b T ζ (6.34 but H does not generate the dynamcs (the dynamcs s non-canoncal. It can be shown that such dynamcs samples the canoncal ensemble. Although all thermodynamcal propertes could n prncple be determned from the free energy F, t s not possble to calculate drectly F from a MD smulaton. The free energy (and the partton functon and the entropy as well cannot be smply expressed as a thermodynamcal average (unlke the energy. Specalzed algorthms are needed for free energy calculatons. 38
6 6.1.5 Constant-pressure MD Very often we are nterested n smulatng systems kept at a gven pressure P rather that occupyng a fxed volume V. Constant-pressure MD can be obtaned by addng the volume V or, n a more general case, the cell parameters, to the dynamcal varables. In the smple case of a lqud, one defnes a Lagrangan: L = 1 N 2 M (V 1/3 σ Ep ({V 1/3 σ} W V 2 P V (6.35 =1 where σ = R /V 1/3 are scaled varables, P s the desred external pressure, and W s a (fcttous mass for V. For a sold, we may be nterested n knowng the equlbrum unt cell volume and form under a gven stress state (typcally a constant external hydrostatc pressure rather than workng at fxed cell and calculatng the correspondng stress. In ths case one ntroduces a matrx h, formed by the unt cell vectors a : h = (a 1, a 2, a 3, and defnes scaled varables S as S = h 1 R. The extended Lagrangan becomes L = 1 N M Ṡ GṠ E p ({hs} W Trḣt ḣ P V (6.36 =1 where G = h t h s the metrc tensor. The nterest of varable-cell dynamcs for sold-state systems resde n the possblty to smulate structural phase transtons (under appled pressure but also as a functon of temperature. 6.2 Bblography A Molecular Dynamcs Prmer, by Furo Ercoless, An ntroducton to Classcal Molecular Dynamcs, wth examples and sample codes. 39
7 6.3 Car-Parrnello Molecular Dynamcs Implementatons of MD usng frst-prncple nteratomc potental calculated from DFT, as n Eq.(6.1, are wdely used. All the MD machnery developed for classcal nteratomc potentals can be used. However these mplementaton suffer from a serous drawback. MD s qute senstve to the qualty of forces. If the forces are not the dervatves of the energy wth hgh accuracy, the MD smulaton wll have problems, appearng as a drft of quanttes that should be conserved (lke e.g. the energy from ther values. The error on DFT forces s lnear n the selfconsstency error of the charge densty (whle for the DFT energy t s quadratc. As a consequence, a very good and expensve convergence to self-consstency s requred at every tme step. In 1985 Car and Parrnello (CP proposed a dfferent approach. They ntroduced a Lagrangan for both electronc and onc degrees of freedom: L = µ dr 2 ψ k (r M Ṙ 2 E tot ({R}, {ψ} + ( Λ kl ψ 2 k(rψ l (rdr δ kl k k,l (6.37 whch generates the followng set of equatons of moton: µ ψ k = Hψ k l Λ kl ψ l, M R = E tot R (6.38 where µ s a fcttous electronc mass, and the Lagrange multplers Λ kl enforce orthonormalty constrants. The electronc degrees of freedom are, n the typcal mplementaton, expanson coeffcents of KS orbtals nto PW. The forces actng on them at each tme step are determned by the KS Hamltonan calculated from the current values of ψ k and of R. The sum over orbtals for an nsulatng system of n electrons ncludes n/2 states, assumng that spn polarzaton s neglected (every orbtal s occuped by two electrons. Most CP calculatons are done for aperodc systems or for systems havng a large unt cell (or supercell, so that typcally only the Γ pont (k = 0 s used to sample the Brlloun Zone. Note that the entre Hamltonan operator s not requred: only products Hψ are. The forces actng on ons have the Hellmann-Feynman form: E tot R = k ψ k V R ψ k (6.39 where V s the electron-on nteracton (pseudo-potental. Note however that Hellmann- Feynman theorem holds only on the exact ground state. The relaton of Car-Parrnello forces to Hellmann-Feynman forces s explaned n the next secton. Orthonormalty constrants are exactly mposed to the ψ at each tme step, usng an teratve procedure that explots the fact that the loss of orthonormalty at each tme step s small. The smulaton starts by brngng the electrons to the BO surface (that s, to the ground state at fxed ons and proceeds, usng classcal MD technology, on both electronc and onc degrees of freedom. Wth approprate values of µ and δt, the electrons always reman close to the BO surface, whle the ons follow a trajectory that s close to the trajectory they would follow n the BO approxmaton. 40
8 The Car-Parrnello dynamcs has turned out to be very successful especally n the study of low-symmetry stuatons: surfaces, clusters, lquds, dsordered materals, and for the study of chemcal reactons Why Car-Parrnello works The reasons why the Car-Parrnello dynamcs works so effectvely are qute subtle. The dynamcs for the electrons s purely classcal (and fcttous: t has nothng to do wth real electron dynamcs. As a consequence the energy would tend to equpartton between electronc and onc degrees of freedom, causng an energy transfer from onc to electronc degrees of freedom. Ths does not happen (and must not happen, otherwse the electrons wll leave the BO surface even on long smulaton tmes. If we analyze the dynamcs n terms of oscllators, we fnd that the typcal frequences assocated to the fcttous electron dynamcs are gven by ω el (ɛ ɛ j /µ, f there s a gap n the electronc spectrum. For ons, the oscllator frequences are the typcal phonon frequences. It turns out that, for reasonable values of the gap and of the fcttous electron mass µ, the maxmum phonon frequency s much smaller than the mnmum electron frequency: ωmax ph << ωmn el. The energy transfer from onc to electronc degrees of freedom s as a consequence very small even on long tmes. Ths stuaton generates a fast electron dynamcs that keeps the electrons close to the BO surface and averages out the error on the forces, so that the much slower onc dynamcs turns out to be correct (that s, very close to the BO dynamcs one would obtan from hghly converged selfconsstency. A detaled explanaton s contaned n a 1991 paper by Pastore, Smargass, and Buda. If there s no gap n the electronc spectrum, or f the gap s too small, the above pcture breaks down. It may be needed to add separate thermostats to onc and electronc degrees of freedom n order to prevent the flow of energy from the former to the latter Choce of the parameters The choce of the electronc mass µ must strke a compromse between conservaton of adabatcty (favored by small values of µ, see above and maxmum admssble tme step (that s lmted by the maxmum electronc frequency, so that the heavest µ, the smaller ωmax, el the larger δt max. Typcally µ 200 amu (1 amu=1 electron mass. For large gap systems, such as SO 2 or H 2 O, n whch adabatcty problems are mnor, µ may be ncreased up to amu and even more. Such values of µ correspond to a typcal tmestep of fs. In order to ncrease the tme step, t s customary to ntroduce the so-called mass precondtonng. In a PW bass set, the tme step s lmted by hgh-frequency components wth the largest G vector. These components are domnated by the knetc energy h 2 G 2 /2µ. Snce electronc masses are fcttous, t s advantageous to ntroduce a mass that for hgh-frequency components goes lke µ(g µ(1 + G 2. The correspondng equatons of motons are only slghtly more complex. It should be notced, however, that too heavy electron masses adversely affect the qualty of smulaton va an electron drag effect. The electron moton follow the onc moton wth some delay, thus ntroducng a drag force that appears as f the ons were heaver than ther masses. Ths mass renormalzaton must be taken nto account when 41
9 extractng vbratonal frequences from MD runs. In some cases, ths effect can ntroduce a nonneglgble devaton from the true onc dynamcs. 6.4 Bblography D. Marx and J. Hutter, n Modern Methods of Quantum Chemstry, (John von Neumann Instute for Computng, FZ Jülch, 2000 p All you wanted to know about Car-Parrnello molecular dynamcs G. Pastore, E. Smargass, F. Buda, Phys. Rev. A 44, 6334 (1991. The paper that explans why the Car-Parrnello dynamcs work. 42
The 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 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 informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationLecture 4. Macrostates and Microstates (Ch. 2 )
Lecture 4. Macrostates and Mcrostates (Ch. ) The past three lectures: we have learned about thermal energy, how t s stored at the mcroscopc level, and how t can be transferred from one system to another.
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 informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
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 informationSnce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t
8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
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 informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationSTATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1
CHE 54 A. Panagotopoulos STATSTCAL MECHACAL ESEMBLES MCROSCOPC AD MACROSCOPC ARABLES The central queston n Statstcal Mechancs can be phrased as follows: f partcles (atoms, molecules, electrons, nucle,
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 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 informationSTATISTICAL MECHANICS
STATISTICAL MECHANICS Thermal Energy Recall that KE can always be separated nto 2 terms: KE system = 1 2 M 2 total v CM KE nternal Rgd-body rotaton and elastc / sound waves Use smplfyng assumptons KE of
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
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 informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationLecture 7: Boltzmann distribution & Thermodynamics of mixing
Prof. Tbbtt Lecture 7 etworks & Gels Lecture 7: Boltzmann dstrbuton & Thermodynamcs of mxng 1 Suggested readng Prof. Mark W. Tbbtt ETH Zürch 13 März 018 Molecular Drvng Forces Dll and Bromberg: Chapters
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
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 informationTHEOREMS OF QUANTUM MECHANICS
THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn
More informationThermodynamics and statistical mechanics in materials modelling II
Course MP3 Lecture 8/11/006 (JAE) Course MP3 Lecture 8/11/006 Thermodynamcs and statstcal mechancs n materals modellng II A bref résumé of the physcal concepts used n materals modellng Dr James Ellott.1
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 informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationA how to guide to second quantization method.
Phys. 67 (Graduate Quantum Mechancs Sprng 2009 Prof. Pu K. Lam. Verson 3 (4/3/2009 A how to gude to second quantzaton method. -> Second quantzaton s a mathematcal notaton desgned to handle dentcal partcle
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More information10. Canonical Transformations Michael Fowler
10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst
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 informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationLagrangian Field Theory
Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationand Statistical Mechanics Material Properties
Statstcal Mechancs and Materal Propertes By Kuno TAKAHASHI Tokyo Insttute of Technology, Tokyo 15-855, JAPA Phone/Fax +81-3-5734-3915 takahak@de.ttech.ac.jp http://www.de.ttech.ac.jp/~kt-lab/ Only for
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want
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 informationNotes on Frequency Estimation in Data Streams
Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to
More informationDensity matrix. c α (t)φ α (q)
Densty matrx Note: ths s supplementary materal. I strongly recommend that you read t for your own nterest. I beleve t wll help wth understandng the quantum ensembles, but t s not necessary to know t n
More informationPHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions
PHYS 5C: Quantum Mechancs Sprng 07 Problem Set 3 Solutons Prof. Matthew Fsher Solutons prepared by: Chatanya Murthy and James Sully June 4, 07 Please let me know f you encounter any typos n the solutons.
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationPoisson brackets and canonical transformations
rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationCHAPTER 5: Lie Differentiation and Angular Momentum
CHAPTER 5: Le Dfferentaton and Angular Momentum Jose G. Vargas 1 Le dfferentaton Kähler s theory of angular momentum s a specalzaton of hs approach to Le dfferentaton. We could deal wth the former drectly,
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More information5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR
5.0, Prncples of Inorganc Chemstry II MIT Department of Chemstry Lecture 3: Vbratonal Spectroscopy and the IR Vbratonal spectroscopy s confned to the 00-5000 cm - spectral regon. The absorpton of a photon
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
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 informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More information14 The Postulates of Quantum mechanics
14 The Postulates of Quantum mechancs Postulate 1: The state of a system s descrbed completely n terms of a state vector Ψ(r, t), whch s quadratcally ntegrable. Postulate 2: To every physcally observable
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 informationLesson 5: Kinematics and Dynamics of Particles
Lesson 5: Knematcs and Dynamcs of Partcles hs set of notes descrbes the basc methodology for formulatng the knematc and knetc equatons for multbody dynamcs. In order to concentrate on the methodology and
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 informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationON MECHANICS WITH VARIABLE NONCOMMUTATIVITY
ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationPHYS 705: Classical Mechanics. Hamilton-Jacobi Equation
1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationLecture 20: Noether s Theorem
Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external
More informationFinding Dense Subgraphs in G(n, 1/2)
Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More information= z 20 z n. (k 20) + 4 z k = 4
Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5
More information5. Response properties in ab initio schemes
5. Response propertes n ab nto schemes A number of mportant physcal observables s expressed va dervatves of total energy (or free energy) E. Examples are: E R 2 E R a R b forces on the nucle; crtcal ponts
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More information763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.
7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
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 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints
Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or
More informationCelestial Mechanics. Basic Orbits. Why circles? Tycho Brahe. PHY celestial-mechanics - J. Hedberg
PHY 454 - celestal-mechancs - J. Hedberg - 207 Celestal Mechancs. Basc Orbts. Why crcles? 2. Tycho Brahe 3. Kepler 4. 3 laws of orbtng bodes 2. Newtonan Mechancs 3. Newton's Laws. Law of Gravtaton 2. The
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationClassical Mechanics Virtual Work & d Alembert s Principle
Classcal Mechancs Vrtual Work & d Alembert s Prncple Dpan Kumar Ghosh UM-DAE Centre for Excellence n Basc Scences Kalna, Mumba 400098 August 15, 2016 1 Constrants Moton of a system of partcles s often
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 0 Canoncal Transformatons (Chapter 9) What We Dd Last Tme Hamlton s Prncple n the Hamltonan formalsm Dervaton was smple δi δ p H(, p, t) = 0 Adonal end-pont constrants δ t ( )
More informationGrand canonical Monte Carlo simulations of bulk electrolytes and calcium channels
Grand canoncal Monte Carlo smulatons of bulk electrolytes and calcum channels Thess of Ph.D. dssertaton Prepared by: Attla Malascs M.Sc. n Chemstry Supervsor: Dr. Dezső Boda Unversty of Pannona Insttute
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 7 Specal Relatvty (Chapter 7) What We Dd Last Tme Worked on relatvstc knematcs Essental tool for epermental physcs Basc technques are easy: Defne all 4 vectors Calculate c-o-m
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationLecture 10. Reading: Notes and Brennan Chapter 5
Lecture tatstcal Mechancs and Densty of tates Concepts Readng: otes and Brennan Chapter 5 Georga Tech C 645 - Dr. Alan Doolttle C 645 - Dr. Alan Doolttle Georga Tech How do electrons and holes populate
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationThree views of mechanics
Three vews of mechancs John Hubbard, n L. Gross s course February 1, 211 1 Introducton A mechancal system s manfold wth a Remannan metrc K : T M R called knetc energy and a functon V : M R called potental
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationCredit Card Pricing and Impact of Adverse Selection
Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n
More information