Newtonian and extended Lagrangian dynamics

Size: px
Start display at page:

Download "Newtonian and extended Lagrangian dynamics"

Transcription

1 Newtonian and extended Lagrangian dynamics Gianni Cardini and Riccardo Chelli Dipartimento di Chimica Ugo Schiff Università di Firenze, Via della Lastruccia 3, Sesto Fno, Firenze SMART January 25-29, 2016

2 The Molecular Dynamics Method Numerical experiments on model systems You need: to define a model equations of motion numerical integration A COMPUTER G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58 12

3 It is useful! Molecular dynamics methods provide essentially exact results for a model Test of theories in controlled conditions Test of a model MD can suggest new experiments Microscopic description G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58 13

4 How to? Model Ṣtucture, conditions, potential Rule to solve the equations of motions ṇumerical integration alghorithm, time step analysis {r i (t)} {v i (t)} G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58 14

5 Reminder Molecular dynamics gives exact results on a model If the model is wrong the results are wrong! If you ask a wrong question you will receive a wrong answer! G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58 15

6 Microscopic description of a chemical system Initial position and momentum of the particles Interaction law No nuclear quantum effects Classical mechanics No electronic transition Ḅorn Oppenheimer approximation G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58 16

7 Interaction law, Analitical potential V = v 1 (r i ) + v 2 (r i, r j ) + i ij ijk v 4 (r i, r j, r k, r l ) + ijkl one body:v 1 (r i ) ẹxternal field two bodies:v 2 (r i, r j ) ạtom-atom, electrostatics, stretching three bodies:v 3 (r i, r j, r k ), bending four bodies:v 4 (r i, r j, r k, r l ), torsion v 3 (r i, r j, r k )+ 17 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

8 NE ensemble Isolated system Clusters Newtonian mechanics Ẹquation of motion from II law: F = ma = r V (r) Constraints Ḥolonomic constraints [f ({r}, t) = 0] Lagrangian mechanics Ḷ({q}, { q}) = K V two scalars! G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58 18

9 Lagrangian N particles L({q}, { q}) = K V = 1 2 N j=1,n Equations of motion d L L = 0 dt q j q j M holonomic constraints: f k ({r}, t) = 0 M G j = λ k j f k ({r}, t) λ k Lagrange multipliers k=1 m i q 2 j V ({q}) d L L = G j dt q j q j 19 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

10 Hamiltonian From Lagrangian p = L(q, q) q Equations of motion H(q, p) = dq dt p = K + V dq dt dp dt = H p = H q 6N first order diff eq instead of 3N second order G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

11 Condensed phase simulations NVE ensemble Surface effects Periodic Boundary Conditions Angular Momentum no more conserved Spurious correlations ṃinimum image PBC & Electrostatics:Ewald V (ϵ = 1) = V d + V r + V s + V shape V d = 1 2 N i=1 j=1 N erfc(α rij + R ) q i q j r ij + R R=0 V r = 1 2 N i=1 j=1 V s = α N π i=1 q2 i N q i q j (4π/Λ) k 0 1 k 2 exp( k 2 /4α 2 ) cos(k r ij ) V shape = 2π 3Λ N i=1 q 2 ir i G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

12 Other Ensembles P=costant HC, J Chem Phys 72 (1980) 2384 (NPH) M Parrinello, A Rahaman, Phys Rev Lett 45 (1980) 1196 T=costant HC, J Chem Phys 72 (1980) 2384 (MC/MD) S Nosé, Mol Phys (1984) 255 () Extended Lagrangians New dynamical variables are added L = L + K (new v) V (new v) G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

13 Numerical Integration Choice of the time step, t Large as possible ṭo sample the phase space with less steps as possible but sufficiently small to sample the fastest motion and to obtain an acceptable conservation of the constants of the motion G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

14 Numerical integrators and Lioville operators Time evolution of f (x) x ({q, p}) phase space df 3N [ f dt = H f ] H q n p n p n q n n=1 Liouville operator: L 3N ıl = n=1 [ H p n H q n q n p n ] ılf (x) = df (x) dt G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

15 Numerical integrators and Lioville operators formal solution f (x t ) = e ılt f (x 0 ) classical propagator: e ılt ; qm propagator e ıht/ħ Starting point for approximate solutions ıl = ıl 1 + ıl 2 do not commute ıl 1 = ıl 2 = 3N H p n n=1 3N n=1 q n H q n [ıl 1, ıl 2 ] 0 p n 115 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

16 Numerical integrators and Lioville operators e ılt e ıl 1t e ıl 2t ọften the action of e ıl 1t and e ıl 2t can be evaluated exactly symmetric Trotter theorem (1959) given two operators such that [A, B] 0 then ) e (A+B)t = lim (e B t 2m e A t m e B t m 2m m Applying the symmetric Trotter theorem e ılt = lim (e ıl 2 t t 2m e ıl 1 m ) t m m e ıl 2 2m G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

17 Numerical integrators and Lioville operators posing t = t m e ılt = ) lim (e ıl 2 t 2 e ıl 1 t e ıl 2 t m 2 m, t 0 Hans de Raedt e Bart de Raedt: Approximate propagation [ Phys Rev A28 ( 1983) ] Tuckerman et al ṛeversible integrators e ıl t e ıl 2 t 2 e ıl 1 t e ıl 2 t 2 + O( t 3 ) G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

18 Numerical integrators and Liouville operators Example: one dimensional Hamiltonian H = p2 2m + V (q) ıl 2 = F (q) p eıl 2 t/2 = e t 2 F (q) p ıl 1 = p m q eıl1 t = e t p m q Exponential of operators e a t = (a t) k k=0 k! G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

19 Numerical integrators and Liouville operators First operator e t ( ) 2 F(q) p q(0) = p(0) k=0 ( = momentun translation Second operator e t p m ( ) q q(0) p( t 2 ) = position translation ( 1 t k! 2 F (q) ) k ( ) q(0) p p(0) ) ( ) q(0) = p ( ) t 2 q(0) p(0) + t 2 F (q) ( q(0) + t p ) m p( t 2 ) = ( ) q( t) p( t 2 ) G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

20 Numerical integrators and Lioville operators Position is changed then recompute the force Third operator ( ) ( ) e t 2 F (q( t)) p q( t) q( t) p( t 2 ) = p( t 2 ) + t 2 ( ) F (q( t)) q( t) = p( t) momentum translation This is Warning When the system is subject to olonomic constraints (SHAKE, RATTLE ) the Lagrange multipliers have to be determined to full convergence G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

21 Complex systems different time scale Potential Energy Atomic Forces U({r i }) = U intra + U inter F i = r i U = i U intra i U inter = F intra i Internal forces fast motions Intermolecular forces slow motions + F inter i G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

22 System characterized by fast and slow motions Liouville operator ıl = ıl fast + ıl slow ıl fast = p m ıl slow = F slow p Reference Hamiltonian q + Ffast p H ref = p2 2m + U(q)fast F fast (q) = du(q)fast dq q = p m fast G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

23 RESPA, reference system propagator alghorithm Tuckerman et al JCP 97(1992) 1990 Factorization of the full propagator Defining τ = t n e ıl t t ılslow = e 2 e ıl fast t t ılslow e 2 e ılfast t = [ ] e τ 2 Ffast (q) p e τ p m q e τ 2 Ffast (q) n p G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

24 full RESPA propagator e ıl t =e t 2 F slow (q) p [ ] e τ 2 F fast (q) p e τ p m q e τ 2 F fast (q) n p e t 2 F slow (q) p G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

25 Pseudocode p = p + 05 * dt * Fslow do j=1,n p = p + 05 * dt/n * Ffast q = q + dt/n * p/m! modified coordinates call FastForce p = p + 05 * dt/n * Ffast enddo! coordinates call SlowForce p = p + 05 * dt * Fslow! momenta G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

26 Canonical Ensemble N V T System in contact with an infinite thermal bath The H of the system is not conserved The H follows a Boltzmann distribution e βh G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

27 Canonical Ensemble Nosé Lagrangian SNosé, MolPhys 52 (1984) 255 A new degree of freedom: s s describes the interaction with an external bath by a velocity scaling Extended Lagrangian N m L Nose = i s 2 ( ) 2 dri U({r}) 2 dt i + Q ( ) ds 2 (3N + 1)kT eq ln(s) 2 dt G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

28 Canonical Ensemble Nosé Hamiltonian conjugate momenta from the Lagrangian p i L ṙ i = m i s 2 ṙ i p s L ṡ = Qṡ Hamiltonian of the extended system (N particles + s) H Nose = N i p 2 i 2m i s 2 U({r}) + p2 s 2Q (3N + 1)kT eq ln(s) G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

29 Canonical Ensemble Nosé equations of motion virtual variables:p, r, t dr i dt dp i dt ds dt dp s dt = H dp i = = H dr i p i m i s 2 = U({r}) r i = H = p s dp s Q = H ds = ( N i p 2 i m i s 2 (3N + 1)kT eq ) 1 s 129 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

30 Canonical Ensemble real variables r and s correspond to the real variables, while π = p/s, p s = p s /s and τ = t/s equations of motion dr i dτ = s dr i dt d π i dτ = s dp i/s dt ds dτ = s ds dt dsp s/q dτ = s Q dp s dt = p i m i s = π i m i = dp i dt = sp s Q ( N = i 1 s p i ds dt π 2 i m i (3N + 1)kT eq ) 1 Q G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

31 Canonical Ensemble Constant of the motion in real variables H Nose = N i + (sp s) 2 2Q π 2 i 2m i + U({r}) + (3N + 1)kT eq ln s transformation to real variables is not canonical H Nose It is not an Hamiltonian implementation in real variables eqm are not easy to be implemented in MD G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

32 Canonical Ensemble equations Hoover PhysRevA31 (1985)1695: non canonical change of variables in Nosé eqm π i = p i dt s,dτ = s, 1 ds s dτ = dη dt and p s = p η and posing the number of dof to 3N ṙ i = p i m ṗ i = i U({r}) p η Q p i η = p η Q N p 2 i ṗ η = 3NkT eq m i i G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

33 Canonical Ensemble η acts as a friction term ṗ η twice the difference of kinetic energy and its canonical average It is a non Hamiltonian system Conserved Energy H(r, p, η, p η ) = H(r, p) + p2 η 2Q + 3NkT eq The real Hamiltonian is H(r, p) When N i F i = 0 also e η N i F i is constant and the distribution function of the momenta is wrong f (p) 1 2πmkT e p2 2mkT 133 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

34 Canonical Ensemble Martyna et al J Chem Phys 97 (1992) 2635 The wrong distribution function is due to the existence of two conservation laws The eqm do not have a sufficient number of variables! Chain The f (p η ) must follow a Maxwell-Boltzmann distribution p η must be coupled to a thermostat The chain of thermostats should be infinite The lenght is chosen finite (M thermostats) G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

35 Canonical Ensemble Chain eqm ṙ i = p i m ṗ i = F i p η 1 Q 1 p i η j = p η j ṗ ηm = Q j ṗ η1 = N i=1,n j = 1,, M p 2 i m i 3NkT [ p 2 ηj 1 ṗ ηj = kt Q j 1 [ p 2 ηm 1 kt Q M 1 ] p η 2 Q 2 p η1 p η j+1 p ηj j = 2,, M 1 Q j+1 ] 135 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

36 Canonical Ensemble Chain Previous equations cannot be transformed in an Hamiltonian system Optimal choice: Q 1 = 3NkT τ 2, Q j = NkT τ 2 j = 2,, M τ > 10 t Conserved Energy H = H(r, p) + M j=1 p 2 η j 2Q j + 3NkT η 1 + kt M j=2 η j G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

37 Canonical Ensemble Massive Chain add a separate chain to each atom of the system much more rapid equilibration some modes are often weakly coupled Tobias et al J Phys Chem 97 (1993) rapid thermalization of a protein in solution G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

38 Canonical Ensemble Integration of Chain eqm The Liouville operator can be written as: ıl = ıl NHC + ıl 1 + ıl 2 e ıl t = e ılnhc t/2 e ıl2 t/2 e ıl1 t e ıl2 t/2 e ılnhc t/2 +O( t 3 ) e ıl H t = e ıl 2 t/2 e ıl 1 t e ıl 2 t/2 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

39 Canonical Ensemble Thermostat forces G 1 = N i=1,n p 2 i G j = p2 η j 1 Q j 1 kt 3NkT m i ıl NHC = + M 1 j=1 N i=1 p η1 p i + Q 1 p i ( G j p ηj p ηj+1 Q j+1 M j=1 p ηj Q j η j ) p ηj + G M p ηm G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

40 Canonical Ensemble Factorization of the propagator is not sufficient! The resulting integrator is not sufficiently robust The application of RESPA alone require many cycles An improvement has been obtained coupling RESPA with an higher order fatorization scheme: Suzuki-Yoshida [Phys Lett A150(1990)262,J Math Phys 32(1991)400] G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

41 Canonical Ensemble Suzuki-Yoshida Given a primitive factorization, S(λ), of an operator S(λ) = e λa 2/2 e λa 2 e λa 1/2 = e λ(a 1+A 2 ) chosen an even order of the error in the factorization 2s this gives n SY = 5 s 1 weigths, w α, such that: n SY w α = 1 α=1 The factorization is e λ(a 1+A 2 ) n SY α=1 S(w α λ) G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

42 Canonical Ensemble RESPA Setting λ = t/2 e ıl NHC t/2 n SY α=1 S(w α t/2) the operator is to be applied n times with a time step w α t 2n e ıl NHC t/2 [ nsy α=1 ( ) ] n t S w α 2n G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

43 Canonical Ensemble Factorization of e ıl NHC t/2 The choice is not unique A new kind of operator e cx x Given a function p, the action of the operator is [ e cp c k ( p p = p ) ] k p k! p k=0 = p k=0 c k k! = pec Given a function f (p), the results is f (pe c ) Terms of this kind will act scaling the momentum of the thermal bath G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

44 Isoenthalpic-Isobaric Ensemble HC, JChem Phys 72 (1980) 2384 First attempt to perform MD simulations at costant P The volume,v, is treated as a dynamical variable Cartesian coordinates {r i } are replaced by scaled coordinates { ρ i }: ρ i = r i V 1/3 Postulated Lagrangian: L({ ρ}, { ρ}, V, V ) = 1 N 2/3 mv 2 ρ i ρ i i=1 N i<j=1 U(V 1/3 ρ ij ) M P V 2 P ext V with M P piston mass and P ext external pressure G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

45 Isoenthalpic-Isobaric Ensemble Hamiltonian H = 1 2mV 2/3 N i= M P Π 2 + P ext V real variables Conserved Energy N 1 π i 2 + i=1 N U(V 1/3 ρ ij ) j>i with p V = M P V H = H 0 (r, p) + p2 V 2M P + P ext V G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

46 Isoenthalpic-Isobaric Ensemble eqm Martyna, Tobias, Klein JChemPhys 101(1994)4177 a variable ϵ = 1 3 ln(v /V V 0) with momentum p ϵ = M P 3V ṙ i = p i m i + p ϵ M P r i ṗ i = F i p ϵ M P p i V = 3Vp ϵ M P ( 1 + ṗ ϵ = 3V (P int P ext ) + 3 N dof ) 3 (3N N constr ) N p 2 i m i i=1 with P int = 1 3 Tr[Pint αβ ] and Pint αβ = 1 [ ] N piα p iβ V i=1 m i + F iα r iβ 146 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

47 Isoenthalpic-Isobaric Ensemble Parrinello - Rahman JApplPhys 52(1981)7182 Parrinello and Rahaman same estended Lagrangian of but with variable cell shape transformation matrix from scaled,s, to cartesian,r a x b x c x H = a y b y c y a z b z c z r = Hs real distances by metric matrix G = H t H V = det H, cell sides a = ax 2 + ay 2 + az 2 9 elements of the matrix are used as dynamical variables 3 sides and 3 angles of the cell 3 Euler angles!! 147 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, / 58

48 Isoenthalpic-Isobaric Ensemble Molecular case Nosé and Klein Mol Phys 50 (1983) 1055 No rotating simulation box H 11 H 12 H 13 H = 0 H 22 H H 33 H 11 = a H 12 = b cos(γ) H 13 = c cos(β) H 22 = b sin(γ) H 23 = cos(α) cos(β) cos(γ) c sin(γ) H 33 = c 2 H13 2 H2 23 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

49 Isoenthalpic-Isobaric Ensemble Molecular case Center of mass scaling Small molecules center of mass variables: R α ; P α ; F α P internal molec = 1 3V Molec α=1 [ ] P 2 α + R α M F α α G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

50 Isobaric Ensemble Phys Rev B44 (1991) 2358 Dynamical variable Stress Tensor Very useful to optimize crystal structure G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

51 es Molecules { N n } L = s 2 m iα 2 σ iα 2 + M α 2 s2 (H ρ α ) t (H ρ α ) α=1 i=1 U({ σ iα + H ρ α }) + W P 2 Tr(Ḣt Ḣ) + W s 2 ṡ2 P ext det(h) 3NkT ext ln s Lagrangian with a simple Nosé Nosé-Hoover chain separated thermostats for atoms and barostat G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

52 factorization Tuckerman et al J Phys A39 (2006) 5629 factorization scheme for an atomic system e ıl t = e ıl NHC baro t/2 e ıl NHC part t/2 e ıl g,2 t/2 p g = W p ḢH 1 ıl 1 = ıl 2 = e ıl 2 t/2 e ıl g,1 t e ıl 1 t e ıl 2 t/2 e ıl g,2 t/2 e ıl NHC part t/2 e ıl NHC baro t/2 N i=1 N i=1 [ pi + p ] g r i m i W p [ ( F i p g W p + 1 ) ] Tr[p g ] I p i N f W P p i G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

53 ıl g,1 = p gh W p H [ ıl g,1 = ıl NHC part = M 1 j=1 det[h](p int IP ext ) + 1 N f N i=1 p η1 p i + Q 1 p i ( G j p ηj p ηj+1 Q j+1 M j=1 p ηj Q j N i=1 η j p 2 i m i I ) p ηj + G M p ηm ] p g ıl NHC baro = same as ıl NHC part with p g G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

54 pseudo masses For T of particles: For T of barostat For pressure Q 1 = 3kT τ 2 p Q 1 = kt τ 2 b Q j = kt τ 2 p Q j = kt τ 2 b Martyna 1992, 1996 W P = (3N + 3)kT τ 2 b G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

55 First Principles Molecular Dynamics: Unified Approach for Molecular Dynamics and Density-Functional Theory R Car and M Parrinello Phys Rev Lett 55 (1985) 2471 extended Lagrangian L = + ij occorb i µ i ψ i (r) 2 dr + 1 M α Ṙ 2 α E[{ψ i }, R α ] 2 α ( ) λ ij ψi (r)ψ i(r)dr δ ij An intertial factor µ i, pseudomass (mass lenght 2 ), is associated to the electronic degree of freedom G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

56 First Principles Molecular Dynamics: : equations of motion Assumption: the system is on the BO surface µ ψ i = δe occorb + λ ij ψ j i = 1,, occorb (1) δ ψ i Hellman-Feynman j=1 M I RI = RI E I = 1,, Atoms (2) RI E({R I }) = Ψ 0 RI H Ψ 0 G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

57 First Principles Molecular Dynamics: : trajectories The trajectories generated by CP equations does not corresponds to the true ones unless E[{ψ i }, {R I }] is in the minimum respect {ψ i } at each time step This is obtained by choosing the value of µ to obtain a decoupling between electronic and nuclear degree of freedom G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

58 Rare events Metadynamics Laio, Parrinello PNAS 99(2002) Iannuzzi, Laio, Parrinello L 90(2003) Estended Lagrangian collective variables s α L MTD = L + α 1 2 M αṡ 2 α α 1 2 k α[s α s α ] 2 V (t, s) k α is the force constant that couple the collective variables to the system V (t, s) is a time dependent potential arising from the accumulation of repulsive gaussian hills modified every MD steps G Cardini, R Chelli (UFirenze) MD SMART January 25-29, /

7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim

7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim 69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the

More information

MD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky

MD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order

More information

Enhanced sampling via molecular dynamics II: Unbiased approaches

Enhanced sampling via molecular dynamics II: Unbiased approaches Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY

More information

Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics

Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle

More information

Density Functional Theory: from theory to Applications

Density Functional Theory: from theory to Applications Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics

More information

Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014

Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014 Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid

More information

Energy and Forces in DFT

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

What is Classical Molecular Dynamics?

What is Classical Molecular Dynamics? What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated

More information

Ab initio Molecular Dynamics Born Oppenheimer and beyond

Ab initio Molecular Dynamics Born Oppenheimer and beyond Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator

More information

Ab initio molecular dynamics

Ab initio molecular dynamics Ab initio molecular dynamics Molecular dynamics Why? allows realistic simulation of equilibrium and transport properties in Nature ensemble averages can be used for statistical mechanics time evolution

More information

Temperature and Pressure Controls

Temperature and Pressure Controls Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically

More information

Crossing the barriers - simulations of activated processes

Crossing the barriers - simulations of activated processes Crossing the barriers - simulations of activated processes Mgr. Ján Hreha for 6 th Student Colloquium and School on Mathematical Physics Faculty of Mathematics, Physics and Informatics Comenius University

More information

Free energy calculations and the potential of mean force

Free energy calculations and the potential of mean force Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science

More information

Ab initio molecular dynamics and nuclear quantum effects

Ab initio molecular dynamics and nuclear quantum effects Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and

More information

Physical Dynamics (SPA5304) Lecture Plan 2018

Physical Dynamics (SPA5304) Lecture Plan 2018 Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle

More information

Javier Junquera. Statistical mechanics

Javier Junquera. Statistical mechanics Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic

More information

Hamiltonian flow in phase space and Liouville s theorem (Lecture 5)

Hamiltonian flow in phase space and Liouville s theorem (Lecture 5) Hamiltonian flow in phase space and Liouville s theorem (Lecture 5) January 26, 2016 90/441 Lecture outline We will discuss the Hamiltonian flow in the phase space. This flow represents a time dependent

More information

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties

More information

The Particle-Field Hamiltonian

The Particle-Field Hamiltonian The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and

More information

Ab Initio Molecular Dynamic

Ab Initio Molecular Dynamic Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master / RFCT Ab Initio Molecular Dynamic p.1/67 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD)

More information

Biomolecular modeling I

Biomolecular modeling I 2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,

More information

Under evolution for a small time δt the area A(t) = q p evolves into an area

Under evolution for a small time δt the area A(t) = q p evolves into an area Physics 106a, Caltech 6 November, 2018 Lecture 11: Hamiltonian Mechanics II Towards statistical mechanics Phase space volumes are conserved by Hamiltonian dynamics We can use many nearby initial conditions

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble

More information

Physics 5153 Classical Mechanics. Canonical Transformations-1

Physics 5153 Classical Mechanics. Canonical Transformations-1 1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant

More information

Molecular Dynamics. Park City June 2005 Tully

Molecular Dynamics. Park City June 2005 Tully Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:

More information

Density Functional Theory: from theory to Applications

Density Functional Theory: from theory to Applications Density Functional Theory: from theory to Applications Uni Mainz May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials

More information

Ab initio Molecular Dynamics

Ab initio Molecular Dynamics Ab initio Molecular Dynamics Jürg Hutter Department of Chemistry, University of Zurich Overview Equations of motion and Integrators General Lagrangian for AIMD Car Parrinello MD and Born Oppenheimer MD

More information

Temperature and Pressure Controls

Temperature and Pressure Controls Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are

More information

Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism

Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger

More information

We start with some important background material in classical and quantum mechanics.

We start with some important background material in classical and quantum mechanics. Chapter Basics We start with some important background material in classical and quantum mechanics.. Classical mechanics Lagrangian mechanics Compared to Newtonian mechanics, Lagrangian mechanics has the

More information

Conical Intersections. Spiridoula Matsika

Conical Intersections. Spiridoula Matsika Conical Intersections Spiridoula Matsika The Born-Oppenheimer approximation Energy TS Nuclear coordinate R ν The study of chemical systems is based on the separation of nuclear and electronic motion The

More information

arxiv:gr-qc/ v2 6 Apr 1999

arxiv:gr-qc/ v2 6 Apr 1999 1 Notations I am using the same notations as in [3] and [2]. 2 Temporal gauge - various approaches arxiv:gr-qc/9801081v2 6 Apr 1999 Obviously the temporal gauge q i = a i = const or in QED: A 0 = a R (1)

More information

HAMILTON S PRINCIPLE

HAMILTON S PRINCIPLE HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates

More information

From quantum to classical statistical mechanics. Polyatomic ideal gas.

From quantum to classical statistical mechanics. Polyatomic ideal gas. From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,

More information

Ab Initio Molecular Dynamic

Ab Initio Molecular Dynamic Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master Ab Initio Molecular Dynamic p.1/67 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD) 1 - Which

More information

Molecular Dynamics at Constant Pressure: Allowing the System to Control Volume Fluctuations via a Shell Particle

Molecular Dynamics at Constant Pressure: Allowing the System to Control Volume Fluctuations via a Shell Particle Entropy 2013, 15, 3941-3969; doi:10.3390/e15093941 Review OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Molecular Dynamics at Constant Pressure: Allowing the System to Control Volume

More information

Aspects of nonautonomous molecular dynamics

Aspects of nonautonomous molecular dynamics Aspects of nonautonomous molecular dynamics IMA, University of Minnesota, Minneapolis January 28, 2007 Michel Cuendet Swiss Institute of Bioinformatics, Lausanne, Switzerland Introduction to the Jarzynski

More information

IV. Classical Molecular Dynamics

IV. Classical Molecular Dynamics IV. Classical Molecular Dynamics Basic Assumptions: 1. Born-Oppenheimer Approximation 2. Classical mechanical nuclear motion Unavoidable Additional Approximations: 1. Approximate potential energy surface

More information

Accelerated Quantum Molecular Dynamics

Accelerated Quantum Molecular Dynamics Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended

More information

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations, Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal

More information

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model

More information

Set the initial conditions r i. Update neighborlist. r i. Get new forces F i

Set the initial conditions r i. Update neighborlist. r i. Get new forces F i Set the initial conditions r i t 0, v i t 0 Update neighborlist Get new forces F i r i Solve the equations of motion numerically over time step t : r i t n r i t n + v i t n v i t n + Perform T, P scaling

More information

Rigid-Body Dynamics in the Isothermal-Isobaric Ensemble: A Test on the Accuracy and Computational Efficiency

Rigid-Body Dynamics in the Isothermal-Isobaric Ensemble: A Test on the Accuracy and Computational Efficiency Rigid-Body Dynamics in the Isothermal-Isobaric Ensemble: A Test on the Accuracy and Computational Efficiency WATARU SHIODA, MASUHIRO MIKAMI Research Institute for Computational Sciences (RICS) and Research

More information

Gear methods I + 1/18

Gear methods I + 1/18 Gear methods I + 1/18 Predictor-corrector type: knowledge of history is used to predict an approximate solution, which is made more accurate in the following step we do not want (otherwise good) methods

More information

ORAC: A Molecular Dynamics Program to Simulate Complex Molecular Systems at the atomistic level

ORAC: A Molecular Dynamics Program to Simulate Complex Molecular Systems at the atomistic level ORAC 6 User Manual Release 6.0 ORAC: A Molecular Dynamics Program to Simulate Complex Molecular Systems at the atomistic level Authors and copyright holders: Piero Procacci Massimo Marchi Contributors:

More information

Nosé-Hoover chain method for nonequilibrium molecular dynamics simulation

Nosé-Hoover chain method for nonequilibrium molecular dynamics simulation PHYSICAL REVIEW E VOLUME 61, NUMBER 5 MAY 000 Nosé-Hoover chain method for nonequilibrium molecular dynamics simulation A. C. Brańka Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego

More information

Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases

Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song

More information

Supporting Information

Supporting Information Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos

More information

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore

More information

Path Integral for Spin

Path Integral for Spin Path Integral for Spin Altland-Simons have a good discussion in 3.3 Applications of the Feynman Path Integral to the quantization of spin, which is defined by the commutation relations [Ŝj, Ŝk = iɛ jk

More information

REVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!

REVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference

More information

Ab Ini'o Molecular Dynamics (MD) Simula?ons

Ab Ini'o Molecular Dynamics (MD) Simula?ons Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties

More information

Feynman s path integral approach to quantum physics and its relativistic generalization

Feynman s path integral approach to quantum physics and its relativistic generalization Feynman s path integral approach to quantum physics and its relativistic generalization Jürgen Struckmeier j.struckmeier@gsi.de, www.gsi.de/ struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme

More information

A Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential

More information

Computational Chemistry - MD Simulations

Computational Chemistry - MD Simulations Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated

More information

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately

More information

G : Statistical Mechanics

G : Statistical Mechanics G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the

More information

Physical Dynamics (PHY-304)

Physical Dynamics (PHY-304) Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.

More information

Finite-temperature Field Theory

Finite-temperature Field Theory Finite-temperature Field Theory Aleksi Vuorinen CERN Initial Conditions in Heavy Ion Collisions Goa, India, September 2008 Outline Further tools for equilibrium thermodynamics Gauge symmetry Faddeev-Popov

More information

Handout 10. Applications to Solids

Handout 10. Applications to Solids ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem

More information

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017 CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and

More information

Langevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany

Langevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace

More information

1 Hamiltonian formalism

1 Hamiltonian formalism 1 Hamiltonian formalism 1.1 Hamilton s principle of stationary action A dynamical system with a finite number n degrees of freedom can be described by real functions of time q i (t) (i =1, 2,..., n) which,

More information

Tensors, and differential forms - Lecture 2

Tensors, and differential forms - Lecture 2 Tensors, and differential forms - Lecture 2 1 Introduction The concept of a tensor is derived from considering the properties of a function under a transformation of the coordinate system. A description

More information

Path Integral Quantization of the Electromagnetic Field Coupled to A Spinor

Path Integral Quantization of the Electromagnetic Field Coupled to A Spinor EJTP 6, No. 22 (2009) 189 196 Electronic Journal of Theoretical Physics Path Integral Quantization of the Electromagnetic Field Coupled to A Spinor Walaa. I. Eshraim and Nasser. I. Farahat Department of

More information

15. Hamiltonian Mechanics

15. Hamiltonian Mechanics University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 15. Hamiltonian Mechanics Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License

More information

Generalized Coordinates, Lagrangians

Generalized Coordinates, Lagrangians Generalized Coordinates, Lagrangians Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 August 10, 2012 Generalized coordinates Consider again the motion of a simple pendulum. Since it is one

More information

Physics 212: Statistical mechanics II Lecture IV

Physics 212: Statistical mechanics II Lecture IV Physics 22: Statistical mechanics II Lecture IV Our program for kinetic theory in the last lecture and this lecture can be expressed in the following series of steps, from most exact and general to most

More information

Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges

Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges Q. Ma, NASA/GISS R. H. Tipping, Univ. of Alabama C. Boulet, Univ. Paris-Sud, France Statements The presentation is not a

More information

Lecture I: Constrained Hamiltonian systems

Lecture I: Constrained Hamiltonian systems Lecture I: Constrained Hamiltonian systems (Courses in canonical gravity) Yaser Tavakoli December 15, 2014 1 Introduction In canonical formulation of general relativity, geometry of space-time is given

More information

Ab initio molecular dynamics: Basic Theory and Advanced Met

Ab initio molecular dynamics: Basic Theory and Advanced Met Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site

More information

Appendix C extra - Concepts of statistical physics Cextra1-1

Appendix C extra - Concepts of statistical physics Cextra1-1 Appendix C extra - Concepts of statistical physics Cextra1-1 Appendix C extra - Concepts of statistical physics A conventional introductory course in physics concentrates on one-body motion: the change

More information

Quantum Interference and Selectivity through Biological Ion Channels

Quantum Interference and Selectivity through Biological Ion Channels 1 SUPPLEMENTARY INFORMATION Quantum Interference and Selectivity through Biological Ion Channels Vahid Salari 1,3,*, Hamidreza Naeij 2, Afshin Shafiee 2,3 1 Department of Physics, Isfahan University of

More information

Vector Fields. It is standard to define F µν = µ ϕ ν ν ϕ µ, so that the action may be written compactly as

Vector Fields. It is standard to define F µν = µ ϕ ν ν ϕ µ, so that the action may be written compactly as Vector Fields The most general Poincaré-invariant local quadratic action for a vector field with no more than first derivatives on the fields (ensuring that classical evolution is determined based on the

More information

On explicit integration of two non-holonomic problems

On explicit integration of two non-holonomic problems On explicit integration of two non-holonomic problems Alexey V. Borisov 1 1 Institute of Computer Sciences, Izhevsk, Russia Generalized Chaplygin systems Equations of motion of the generalized Chaplygin

More information

Chapter 1. Principles of Motion in Invariantive Mechanics

Chapter 1. Principles of Motion in Invariantive Mechanics Chapter 1 Principles of Motion in Invariantive Mechanics 1.1. The Euler-Lagrange and Hamilton s equations obtained by means of exterior forms Let L = L(q 1,q 2,...,q n, q 1, q 2,..., q n,t) L(q, q,t) (1.1)

More information

Sketchy Notes on Lagrangian and Hamiltonian Mechanics

Sketchy Notes on Lagrangian and Hamiltonian Mechanics Sketchy Notes on Lagrangian and Hamiltonian Mechanics Robert Jones Generalized Coordinates Suppose we have some physical system, like a free particle, a pendulum suspended from another pendulum, or a field

More information

Low-temperature isomers of the water hexamer

Low-temperature isomers of the water hexamer Low-temperature isomers of the water hexamer Volodymyr Babin and Francesco Paesani Department of Chemistry & Biochemistry, University of California, San Diego, La Jolla, CA 92093-0314 Hundreds of years

More information

BCMB/CHEM Spin Operators and QM Applications

BCMB/CHEM Spin Operators and QM Applications BCMB/CHEM 8190 Spin Operators and QM Applications Quantum Description of NMR Experiments Not all experiments can be described by the Bloch equations A quantum description is necessary for understanding

More information

Scientific Computing II

Scientific Computing II Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the

More information

etc., etc. Consequently, the Euler Lagrange equations for the Φ and Φ fields may be written in a manifestly covariant form as L Φ = m 2 Φ, (S.

etc., etc. Consequently, the Euler Lagrange equations for the Φ and Φ fields may be written in a manifestly covariant form as L Φ = m 2 Φ, (S. PHY 396 K. Solutions for problem set #3. Problem 1a: Let s start with the scalar fields Φx and Φ x. Similar to the EM covariant derivatives, the non-abelian covariant derivatives may be integrated by parts

More information

Attempts at relativistic QM

Attempts at relativistic QM Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and

More information

Nanoscale simulation lectures Statistical Mechanics

Nanoscale simulation lectures Statistical Mechanics Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -

More information

Lecture 4: Equations of motion and canonical quantization Read Sakurai Chapter 1.6 and 1.7

Lecture 4: Equations of motion and canonical quantization Read Sakurai Chapter 1.6 and 1.7 Lecture 4: Equations of motion and canonical quantization Read Sakurai Chapter 1.6 and 1.7 In Lecture 1 and 2, we have discussed how to represent the state of a quantum mechanical system based the superposition

More information

FYS 3120: Classical Mechanics and Electrodynamics

FYS 3120: Classical Mechanics and Electrodynamics FYS 3120: Classical Mechanics and Electrodynamics Formula Collection Spring semester 2014 1 Analytical Mechanics The Lagrangian L = L(q, q, t), (1) is a function of the generalized coordinates q = {q i

More information

Generalized valence bond molecular dynamics at constant temperature

Generalized valence bond molecular dynamics at constant temperature M o l e c u l a r Ph y s i c s, 1996, Vo l. 89, No. 5, 1265± 1276 Generalized valence bond molecular dynamics at constant temperature By DOUGLAS A. GIBSON and EM ILY A. CARTER D epartment of Chemistry

More information

Quantization of Scalar Field

Quantization of Scalar Field Quantization of Scalar Field Wei Wang 2017.10.12 Wei Wang(SJTU) Lectures on QFT 2017.10.12 1 / 41 Contents 1 From classical theory to quantum theory 2 Quantization of real scalar field 3 Quantization of

More information

Citation for published version (APA): Hess, B. (2002). Stochastic concepts in molecular simulation Groningen: s.n.

Citation for published version (APA): Hess, B. (2002). Stochastic concepts in molecular simulation Groningen: s.n. University of Groningen Stochastic concepts in molecular simulation Hess, Berk IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check

More information

Introduction and Vectors Lecture 1

Introduction and Vectors Lecture 1 1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum

More information

PHYS 352 Homework 2 Solutions

PHYS 352 Homework 2 Solutions PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables

More information

4 Quantum Mechanical Description of NMR

4 Quantum Mechanical Description of NMR 4 Quantum Mechanical Description of NMR Up to this point, we have used a semi-classical description of NMR (semi-classical, because we obtained M 0 using the fact that the magnetic dipole moment is quantized).

More information

(Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) p j . (5.1) !q j. " d dt = 0 (5.2) !p j . (5.

(Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) p j . (5.1) !q j.  d dt = 0 (5.2) !p j . (5. Chapter 5. Hamiltonian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 5.1 The Canonical Equations of Motion As we saw in section 4.7.4, the generalized

More information

Distributions of statistical mechanics

Distributions of statistical mechanics CHAPTER II Distributions of statistical mechanics The purpose of Statistical Mechanics is to explain the thermodynamic properties of macroscopic systems starting from underlying microscopic models possibly

More information

Vibrational Spectra (IR and Raman) update Tinoco has very little, p.576, Engel Ch. 18, House Ch. 6

Vibrational Spectra (IR and Raman) update Tinoco has very little, p.576, Engel Ch. 18, House Ch. 6 Vibrational Spectra (IR and Raman)- 2010 update Tinoco has very little, p.576, Engel Ch. 18, House Ch. 6 Born-Oppenheimer approx. separate electron-nuclear Assume elect-nuclear motion separate, full wave

More information

Introduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou

Introduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties

More information

Statistical Mechanics. Atomistic view of Materials

Statistical Mechanics. Atomistic view of Materials Statistical Mechanics Atomistic view of Materials What is statistical mechanics? Microscopic (atoms, electrons, etc.) Statistical mechanics Macroscopic (Thermodynamics) Sample with constrains Fixed thermodynamics

More information

Enabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package

Enabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package Enabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package Mario Fernández Pendás MSBMS Group Supervised by Bruno Escribano and Elena Akhmatskaya BCAM 18 October

More information

d 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5)

d 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5) LECTURE 12 Maxwell Velocity Distribution Suppose we have a dilute gas of molecules, each with mass m. If the gas is dilute enough, we can ignore the interactions between the molecules and the energy will

More information

Mixing Quantum and Classical Mechanics: A Partially Miscible Solution

Mixing Quantum and Classical Mechanics: A Partially Miscible Solution Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate

More information