Ab Initio Molecular Dynamic
|
|
- Elvin Little
- 5 years ago
- Views:
Transcription
1 Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master Ab Initio Molecular Dynamic p.1/67
2 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD) 1 - Which one? 2 - Car-Parinello! IV. Application : Free Energy Calculations V. Conclusions Ab Initio Molecular Dynamic p.2/67
3 Introduction Goal : Compute statistical averages, as in Monte Carlo. Sample phase space Macroscopic properties Tool : Time evolution of the system. Ab Initio Molecular Dynamic p.3/67
4 Microscopic becomes Macroscopic Trajectories sampling phase space A macro = A = A exp( βe)d ri d p i exp( βe)d ri d p i Ergodic hypothesis verified : A macro = A = lim T A T = lim T 1 T T 0 A(t)dt True only if : System at equilibrium Simulation time T much larger than caracteristic times of the system. Ab Initio Molecular Dynamic p.4/67
5 Short History of MD Simulations using balls (hard and soft) Computer Simulations Super-computer MANIAC (Los Alamos in 1952) Monte Carlo : Metropolis, Rosenbluth, Rosenbluth, Teller, and Teller (1953) MD Hard Sphere : Alder and Wainwright (1957) MD with actual movement equations : Rahman (Argonne 1964) Ab Initio Molecular Dynamic Ehrenfest 1927 : First approach. C. Leforestier 1978 : Born-Oppenheimer Car and Parrinello (Sissa 1985) Ab Initio Molecular Dynamic p.5/67
6 Computer simulations Four steps : Initialization Propagation Equilibration Analysis Propagation : All schemes based on numerical solution m i r i = i V ({ r i }) Which V? Empirical : classical force fields Quantum : semi-empirical, ab initio and DFT How to obtain V? V pre-calculated and tabulated : limited to few degrees of freedom, allows for quantum nuclei ; V calculated on the fly : mostly classical nuclei Ab Initio Molecular Dynamic p.6/67
7 Time scale V Simulation duration (ps) Atoms Force Fields Semi-empirical DFT, HF post-hf 1 10 Ab Initio Molecular Dynamic p.7/67
8 Classical Molecular Dynamic Ab Initio Molecular Dynamic p.8/67
9 Equations of motion Many formalism : Newton m i r i = i V ({ r i }) Lagrange Lagrangian for a system with n degrees of freedom L({ r i },{ ri }) = K({ r i },{ ri }) V ({ r i }) Equation of motion : Euler-Lagrange d dt L L r i = 0 r i Impulsion p i = L, Force F i = L r i r i Newton is foudn again : p i = F i Hamilton Ab Initio Molecular Dynamic p.9/67
10 Equations of motion Many formalism : Newton m i r i = i V ({ r i }) Lagrange Hamilton Impulsion p i used instead of velocity Hamiltonian H({ r i }, { p i }) = 3n i=1 q i p i L r i = K({ r i }, { p i }) + V ({ r i }) Ab Initio Molecular Dynamic p.9/67
11 Equations of motion Many formalism : Newton m i r i = i V ({ r i }) Lagrange Hamilton Equations of motion q i = H p i p i = H q i H is constant along the trajectory, interpreted as the total energy dh dt = [ H r i + H ] [ H H p i H ] H = 0 r i p i r i p i p i q i i i Conclusion : all equivalent! Ab Initio Molecular Dynamic p.9/67
12 Working ingredients Initialization : { r i },{ p i } at t = 0 { r 0 i } : from cristallographic structure, PDB or by similarity to avoid very repulsive part of the potential. { p 0 i } : most of the time 0, random or taken from Maxwell-Bolztman distribution. How to compute V? Time is discrete, how to choose time step t? As large as possible Smaller than characteristic time of the system Ab Initio Molecular Dynamic p.10/67
13 Empirical potentials Force fields VSEPR extension : Geometries close to a reference : d CC 1, 54 Å, d CH 1, 09 Å, α Ethane : d CC = 1, 536 Å, Propane : d CC = 1, 526 Å, α = 112, 4, Butane : d (e) CC = 1, 533 Å, d(i) CC = 1, 533 Å, α = 112, 8. Energy close the reference energy! E = E 0 + corrections... Ab Initio Molecular Dynamic p.11/67
14 General form of a Force Field Main contributions : E MM = V bonding + V angle + V torsion + V nb +... V bonding Bonding energy V bonding = 1 2 i bondings k r,i (r i r 0 i ) 2 V angle Valence angle (bending) modification V angle = 1 2 i,j bondings k θ,ij (θ ij θ 0 ij) 2 Ab Initio Molecular Dynamic p.12/67
15 General form of a Force Field Main contributions : E MM = V bonding + V angle + V torsion + V nb +... V torsion Torsion energy (rocking) V torsion = 1 2 i [V 1 (1 + cosϕ i )+ V 2 (1 cos 2ϕ i ) + V 3 (1 + cos 3ϕ i )] Ab Initio Molecular Dynamic p.12/67
16 General form of a Force Field Main contributions : E MM = V bonding + V angle + V torsion + V nb +... V torsion Torsion energy Energie Eclipsee E(ψ) Energie Butane : V1=1., V2=0., V3=1. Syn Gauche Decalee120 ψ (angle HCCH) ψ (angle CCCC) Anti Ab Initio Molecular Dynamic p.12/67
17 General form of a Force Field Main contributions : E MM = V bonding + V angle + V torsion + V nb +... V nb Non Bonded interactions Van Der Waals : V V dw (d) = ǫ i,j atomes ( d 0 ij d ij ) 12 2 ( d 0 ij d ij ) 6 Ab Initio Molecular Dynamic p.12/67
18 General form of a Force Field Main contributions : E MM = V bonding + V angle + V torsion + V nb +... V nb Non Bonded interactions electrostatic type : V el = i,j atomes q i q j Dd ij H-bonds Ab Initio Molecular Dynamic p.12/67
19 Popular force fields MM2, MM3 Mainly for organic chemistry. MMFF94 Only for organic chemistre! AMBER, CHARMm Good for biological systems UFF Universal Force Field. Universaly avoided! ESFF Extensible and Systematic Force Field Small benchmark : Ab initio MMFF94 MM3 CHARMm AMBER E 0,4 0,4 0,7 0,8 1,1 d XH - 0,09 0,5 0,1 0,2 Ab Initio Molecular Dynamic p.13/67
20 Propagator A good propagator must : Be quick and allow for large t Time reversible Conserve mechanical energy Compute forces not too frequently Conserve phase-space volume Basic idea : Taylor expansion ri (t + t) = r i (t) + v i (t) t + ai (t) 2 t vi (t + t) = v i (t) + a i (t) t + vi (t) 2 t Ab Initio Molecular Dynamic p.14/67
21 Verlet General expression : ri (t + t) = r i (t) + v i (t) t + ai (t) 2 t ri (t t) = r i (t) ai (t) v i (t) t + 2 t thus r i (t + t) = 2 r i (t) fi (t) r i (t t) + t 2 + O( t 3 ) M i Store position and forces only Accuracy in t 3 Time reversible by construction Velocities indirectly computed (finite differences) : less accurate v i = ( r i (t + t) r i (t t))/(2 t) + O( t 2 ) Ab Initio Molecular Dynamic p.15/67
22 Leap-Frog variant Designed to improve numerical accuracy ri (t + t) = r i (t) + v i (t + t 2 ) t vi (t + t 2 ) = v i (t t 2 ) + a i (t) t vi (t) = ( vi (t + t 2 ) v i (t t 2 )) /2 + O( t 2 ) More stable Ab Initio Molecular Dynamic p.16/67
23 Velocity Verlet Improve velocity accuracy ri (t + t) = r i (t) + v i (t) t + ai (t) 2 t2 vi (t + t) = ai (t) + a i (t + t) v i (t) + t 2 Velocity accuracy in O( t 3 ) Can be seen as a predictor-corrector Ab Initio Molecular Dynamic p.17/67
24 Verlet : Summary Graphical form : Velocity verlet often used because : Simple and efficient, little force computation Accurate in t 3 Time reversible, Symplectic For realistic time step, Velocity Verlet vitesse as good Ab Initio Molecular Dynamic p.18/67
25 Constrained Dynamic Why using constraints : to get larger t (frozen d(ch)), prevent or force some evolution Principal : Lagrange multiplier Constraint defined by σ ({ r i }) = 0, For example : Distance : σ ({ r i }) = r i r j d 0 Difference of two distances : σ ({ r i }) = r i r j r k r l d 0 Coordination number : σ ({ r i }) = n i ({ r i }) n 0 with n i ({ r i }) = j i S( r i r j ) et S(r) = (1 + exp(κ(r r c ))) 1. Ab Initio Molecular Dynamic p.19/67
26 Constrained Dynamic Constraint defined by σ ({ r i }) = 0 Extended Lagrangian : L ({ r i }, { p i }) = L({ r i }, { p i }) α λ α σ α ({ r i }) Equations of motion : m i r i = V r i α λ α σ α r i SHAKE : λ α tq σ α ({ r i (t + t)}) = 0 RATTLE : λ α tq σ α ({ r i (t + t)}) = 0 Ab Initio Molecular Dynamic p.19/67
27 Which ensembles? By default, isolated system E conserved, i.e. NVE Initialize { p i } so that E tot = E target then propagate How can we obtain other ensembles? NVT : E fluctuates Q(N, V, T) = 1 N!h 3N NPT : E and V fluctuate 1 (N, P, T) V 0 N!h 3N Z Z Π i d r i d p i exp( βh) = exp ( βf(n, V, T)) Π i d r i d p i exp( β(h + PV )) = exp ( βg(n, P, T)) Ab Initio Molecular Dynamic p.20/67
28 Canonical ensemble We know tha k T (t) = 2 3Nk B k p k (t) 2 2m k p 2 k 2m k = 3 2 Nk BT leading to At equilibrium, Maxwell-Boltzmann distribution Initialize Velocities from MB distribution Or p 0 i = 0, then heating using "inversed annealing" : vi (t) α v i (t), α > 1. Equilibrate : homogenize T, to obtain T (t) = T Ab Initio Molecular Dynamic p.21/67
29 Canonical ensemble We know tha k T (t) = 2 3Nk B k p k (t) 2 2m k p 2 k 2m k = 3 2 Nk BT leading to At equilibrium, Maxwell-Boltzmann distribution Initialize Equilibrate : homogenize T, to obtain T (t) = T Rescaling : when T (t) T > T, then vi (t) v i (t) T T (t) Brutal! Thermostat : coupling with a thermal bath Ab Initio Molecular Dynamic p.21/67
30 Nosé-Hoover Extended System : a supplementary variable that mimics the thermal bath Equations ri = pi m i ξ = 1 Q ( k p 2 k 2m k gk B T ) pi = F i ξ p i = gk B Q (T T) with g number of degrees of freedom, Q mass of the thermostat Ab Initio Molecular Dynamic p.22/67
31 Nosé-Hoover Extended System Q regulates the speed of exchanges between the bath and the system Q too small : large not wanted oscillations, slow convergence Q too large : slow exchanges. Q, microcanonical ensemble! Optimum when resonating with the system, ie Q gk BT ω 2 Ab Initio Molecular Dynamic p.22/67
32 Nosé-Hoover Chains Nosé-Hoover not always ergodic, response time too slow Chain of M thermostats : ri = pi m i ξ 1 = 1 Q 1 ( k p 2 k 2m k gk B T ) ξ 1 ξ 2 pi = F i ξ 1pi ξ n = 1 Q n ( Qn 1 ξ 2 n 1 k B T ) ξ n ξ n+1 ξ M = 1 ( QM 1 ξn 1 2 k B T ) Q M Ab Initio Molecular Dynamic p.23/67
33 Test 1D oscillator H = p2 2m + mω2 x 2 2, f(p) = β 2πm e βp2 /2m, f(x) = βmω 2 2π e βω2 x 2 /2 1 thermostat 3 thermostats 4 thermostats Ab Initio Molecular Dynamic p.24/67
34 Andersen or Hybrid Monte Carlo Mix DM and Monte Carlo to simulate chocs with the thermal bath Andersen The velocity of a random particle is withdrawn from the MD distribution ; Ergodic by construction, converges to the NVT associated to the MB distribution ; Many particles can be chosen at once Re-attribution period not so important Re-attribution does not kick wafefunction too far from BO surface Hybrid Monte Carlo : all velocities are changed at once Ab Initio Molecular Dynamic p.25/67
35 Langevin Still stochastic, but white noise added directly in the equation of motion for all particles q = p t m ṗ = V (q) ξ m p + 2ξ β dw(t)/dt with W(t) a gaussian white noise Physically (mathematically?) : the most efficient to thermalize, the fastest to converge but very rare in chemistry, and not yet tested in CPMD Ab Initio Molecular Dynamic p.26/67
36 False thermostat : Berendsen Can be seen as a rescaling, with α(t + t) = = [ 1 + t τ [ 1 + ( T T (t + t) 1 t τt (t + t) )] 1/2 (T T (t + t)) T T (t + t) influe α et non α Very used in biochemistry but do not sample the canonical ensemble! Use only for equilibration! ] 1/2 Ab Initio Molecular Dynamic p.27/67
37 Isobaric-Isotherm ensemble This is the closest to chemist experiments Volumebecomes a dynamical variable No stochastic approaches Two families : Isotropic : V change, but not the box shape Hoover, Anderssen Anisotropic :all parameters of the box can change Parrinello-Rahman, Martyna-Tobias-Klein Ab Initio Molecular Dynamic p.28/67
38 Isotropic NPT : Andersen First approach (1980) Mimic the action of a piston on the box : Mass Q Kinetic energy K V = 1 2 Q V 2 Potential energy V V = PV Reduced variables : s = V 1/3 r, ṡ = V 1/3 v Hamiltonian : H V = K + K V + V + V V Ab Initio Molecular Dynamic p.29/67
39 Isotropic NPT : Andersen First approach (1980) Mimic the action of a piston on the box Reduced variables : s = V 1/3 r, ṡ = V 1/3 v Hamiltonian : H V = K + K V + V + V V Equations of motion : f s = mv 1/3 2 3 V = (P P)/Q ṡ V V with P instantaneous pressure P = 2 W = 2 1 i j>i r ijfij Internal Virial. 3V (K W) et Ab Initio Molecular Dynamic p.29/67
40 Isotropic NPT : Andersen First approach (1980) Mimic the action of a piston on the box Reduced variables : s = V 1/3 r, ṡ = V 1/3 v Hamiltonian : H V = K + K V + V + V V Equations of motion : f s = mv 1/3 2 3 V = (P P)/Q ṡ V V In fact, sampling the isobaric-iso-enthalpic ensemble, not so common. Ab Initio Molecular Dynamic p.29/67
41 Isotropic NPT : Hoover Extended system, similar to the Nosé-Hoover thermostat r = V 1/3 s si = pi m i V 1/3 ξ = gk B Q χ = V 3V (T T) pi = F i (ξ + χ) p i χ = P(t) P τ 2 P k BT with g number of degrees of freedom, Q mass of the thermostat, τ P relaxation time of the pressure fluctuations. Ab Initio Molecular Dynamic p.30/67
42 Anisotropic NPT Parrinello-Rahman Idea : reduced variable r i = h s i, with h = [ a, b, c ] and a, b et c cell parameters We introduce G = h t h P pressure tensor : P αβ = 1 V i m i r iα r iβ + i r ijα f ijβ j Hamiltonian H = K Q α β ḣ 2 αβ + V + PV Ab Initio Molecular Dynamic p.31/67
43 Anisotropic NPT Parrinello-Rahman Idea : reduced variable r i = h s i We introduce G = h t h Hamiltonian H = K Q α Equations of motion with Q mass of the box ḣ 2 αβ + V + PV m s = h 1 f mg 1 Ġ s QḦ = (P 1P)V ( h 1) t β Ab Initio Molecular Dynamic p.31/67
44 Equilibration Listening phase We check that variables that should be conserved are stable, no drift : T E RMSD : Root Mean square displacement g(r) : pair distribution function One simple criterion : averages should not depend on the initial time. Ab Initio Molecular Dynamic p.32/67
45 Analysis Production phase Compute the interesting properties : RMSD Structural parameters g(r) Fluctuations C v Vibrational spectra Kinetic constant Average on a given window, larger than fluctuation time (see Practical session!) Thermostat, Barostat are perturbing the averages Ab Initio Molecular Dynamic p.33/67
46 Shall we go beyond classical MD Pros of classical MD : Fast Good geometries, vibrations, transport properties Thermodynamical properties accurate for parameterized systems Cons Parameters! One need to know how atoms are linked : Chemical reactions cannot be studied by standard FF ReaxFF, ReBO, LOTF... No excited states No electronic properties : charge, density,... Ab Initio Molecular Dynamic p.34/67
47 Ab Initio Molecular Dynamics Ab Initio Molecular Dynamic p.35/67
48 Which ab initio MD? Electron explicitly taken into account Let us denote by : M k, Z k, R k the mass, the charge and the position of a nucleus R k ou v k the velocity of a nucleus pk = M k vk the impulsion of a nucleus m e, r i the mass and the position of an electron Ab Initio Molecular Dynamic p.36/67
49 Which ab initio MD? Electron explicitly taken into account We would like to solve ({ Rk } φ, { ) r i } ;t ({ Rk } i = H, { ) ({ Rk } r i } φ, { ) r i } ;t t with H = electrons i i 2 } {{ } T e élec. nuclei Z k elec. elec. 1 Z k Z l R kl + + r i ik r k i j>i ij k l>k }{{} V ne noy. k 2M k k }{{} Ab Initio Molecular T Dynamic N p.36/67
50 TDSCF approach Hypothesis : φ ({ Rk }, { ) r i } ;t Ψ ({ r i } ;t) χ ({ Rk } ) [ i ;t exp t t 0 dt Ẽ e (t ) with the convenient phase factor Ẽ e (t ) = d r i d R k i Ψ ({ r i } ;t) χ ({ Rk } k ) ;t H e Ψ ({ ({ Rk } ) r i } ;t)χ ;t Let us write d r = i d r i and d R = k d R k Ab Initio Molecular Dynamic p.37/67
51 TDSCF approach Hypothesis : φ ({ Rk }, { ) r i } ;t Ψ ({ r i } ;t) χ ({ Rk } ) [ i ;t exp t t 0 dt Ẽ e (t ) with Ẽ e (t ) = d r d R Ψ ({ ({ r i } ;t) χ Rk } ) ;t H e Ψ ({ ({ Rk } ) r i } ;t)χ ;t Ab Initio Molecular Dynamic p.37/67
52 TDSCF approach Hypothesis : ({ Rk } φ, { ) r i } ;t Ψ ({ r i } ;t) χ ({ Rk } ) [ i ;t exp ] dt Ẽ e (t ) t 0 t we get (project on < Ψ, < χ, use d < H > /dt = 0) : i Ψ { t = iψ + d ({ Rχ Rk } i i χ t = k 1 2M k 2 kχ + { ) ({ Rk } ;t V ne χ ;t) } Ψ d r Ψ ({ r i } ;t) H e Ψ ({ } r i } ;t) χ Ab Initio Molecular Dynamic p.37/67
53 TDSCF approach Let ({ us write (A,S : real) Rk } ) ({ Rk } ) [ χ ;t = A ;t exp is ({ Rk } ) ] ;t / we get, after separating real and imaginary part : i Ψ { t = iψ + d ({ Rχ Rk } ) ({ Rk } ;t V ne χ ;t) } Ψ i S t = 1 ( k S) 2 d r Ψ H e Ψ k A 2M k 2M k A k k A t = k 1 M k ( k S) ( k A) k 1 A ( 2 2M ks ) k Ab Initio Molecular Dynamic p.37/67
54 Ehrenfest approach Classical nuclei : = 0 S t + k 1 2M k ( k S) 2 + d r Ψ H e Ψ = 0 Hamilton-Jacobi equation analogy : S t + k 1 2M k (p k ) 2 + V e = 0 Using p k k S M k Rk (t) = k d r Ψ H e Ψ Ab Initio Molecular Dynamic p.38/67
55 Ehrenfest approach Classical nuclei : ({ Rk } χ That is lim =0 ) ;t 2 = k d ({ Rχ } Rk ({ Rk } Rk ;t) χ ( Rk δ ) R k (t) ) ;t = R k (t) i Ψ t = i 2 iψ + V ne Ψ = H e Ψ Ab Initio Molecular Dynamic p.39/67
56 Ehrenfest summary Wavefunction explicitly propagated, coupled to the nuclei M k Rk (t) = k d r Ψ H e Ψ i Ψ t = H eψ No electronic minimization, except for t = 0 Transitions between electronic states are explicitly described t imposed by electrons dynamics Very very small Ehrenfest seldom used Ab Initio Molecular Dynamic p.40/67
57 Born-Oppenheimer Wavefunction is not propagated Time { independent Schrödinger equation is solved for Rk } each : M k Rk (t) = k min Ψ 0 Ψ 0 H e Ψ 0 E 0 Ψ 0 = H e Ψ 0 Electronic minimization at each step No more transitions between electronic states t imposed by the nuclei relatively large Ab Initio Molecular Dynamic p.41/67
58 Born-Oppenheimer : HF, DFT Simplifications using independent electron methods : HF ou DFT-KS Molecular Orbitals denoted by {ϕ i } Energy minimization with orthogonal MO : : L e = Ψ 0 H e Ψ 0 + i,j Λ ij ( ϕ i ϕ j δ ij ) leads to : H eff ϕ i = j Λ ij ϕ j Ab Initio Molecular Dynamic p.42/67
59 Born-Oppenheimer : HF, DFT Simplifications using independent electron methods : HF ou DFT-KS Molecular Orbitals denoted by {ϕ i } Energy minimization with orthogonal MO : leads to : H eff ϕ i = Λ ij ϕ j j Diagonal form : H eff ϕ i = ǫ i ϕ i Ab Initio Molecular Dynamic p.42/67
60 Born-Oppenheimer : HF, DFT Simplifications using independent electron methods : HF ou DFT-KS Molecular Orbitals denoted by {ϕ i } Energy minimization with orthogonal MO : leads to : H eff ϕ i = Λ ij ϕ j j So that : M k Rk (t) = k min Ψ 0 Ψ 0 H e Ψ 0 0 = H eff ϕ i + j Λ ij ϕ j Ab Initio Molecular Dynamic p.42/67
61 Car-Parrinello (1/2) Goal : benefit from all advantages Ehrenfest : No electronic minimization BO : Large time step Tool : Adiabatic separation between fast electrons and slow nuclei How : MO described as classical variables Ab Initio Molecular Dynamic p.43/67
62 Car-Parrinello (1/2) Goal : benefit from all advantages Tool : Adiabatic separation between fast electrons and slow nuclei How : MO described as classical variables, but decoupled from the nuclei Using an extended Lagrangian formulation Ab Initio Molecular Dynamic p.43/67
63 Car-Parrinello (2/2) Goal : benefit from all advantages Tool : Adiabatic separation between fast electrons and slow nuclei How : MO described as classical variables, but decoupled from the nuclei L CP = 1 2 M kr 2 k µ i ϕ i ϕ i k i }{{} Kinetic energy µ "Fictitious" mass of the electrons Ψ 0 H e Ψ 0 }{{} Potential energy + constraints }{{} orthogonality, geometric... Ab Initio Molecular Dynamic p.44/67
64 Car-Parrinello (2/2) Goal : benefit from all advantages Tool : Adiabatic separation between fast electrons and slow nuclei Constant of motion H CP = 1 2 M kr 2 k + 1 { 2 µ Rk }) i ϕ i ϕ i + E el ({ϕ i }, k i Ab Initio Molecular Dynamic p.44/67
65 Car-Parrinello (2/2) Goal : benefit from all advantages Tool : Adiabatic separation between fast electrons and slow nuclei Equations of motion M k Rk (t) = k Ψ 0 H e Ψ 0 + k {constraints} µ ϕ i (t) = δ δϕ i Ψ 0 H e Ψ 0 + δ δϕ i {constraints} t approximately 5 to 10 times smaller than BO Ab Initio Molecular Dynamic p.44/67
66 Car-Parrinello : HF,DFT-KS We use HF or DFT-KS methods L CP = k 1 2 M kr 2 k + i 1 2 µ i ϕ i ϕ i Ψ 0 H eff Ψ 0 + i,j Λ ij ( ϕ i ϕ j δ ij ) Ab Initio Molecular Dynamic p.45/67
67 Car-Parrinello : HF,DFT-KS We use HF or DFT-KS methods Equations of motion become M k Rk (t) = k Ψ 0 H eff Ψ 0 µ ϕ i (t) = H eff ϕ i + j Λ ij ϕ j Very similar to BO : µ ϕ i (t) = 0 Ab Initio Molecular Dynamic p.45/67
68 Why does it work? Fast electrons, slow nuclei Electronic frequencies for Si 8 : f e (ω) = Ψ(t) Ψ(0) dt t=0 cos(ωt) i Triangle : latest nuclear frequency Ab Initio Molecular Dynamic p.46/67
69 Why does it work? Electronic Oscillations close to the BO surface E cons = k E phys = k 1 2 M kr 2 k µ i ϕ i ϕ i + Ψ H e Ψ i 1 2 M kr 2 k + Ψ H e Ψ = E cons T e Model System : Si FCC, 2 atoms/cell Ab Initio Molecular Dynamic p.47/67
70 Why does it work? Electronic Oscillations close to the BO surface Model System : Si FCC, 2 atoms/cell Small oscillations, stable in time Ab Initio Molecular Dynamic p.47/67
71 Why does it work? Electronic Oscillations close to the BO surface Model System : Si FCC, 2 atoms/cell Small oscillations, stable in time Ab Initio Molecular Dynamic p.47/67
72 Why does it work? Forces oscillations very small Model System : Si FCC, 2 atoms/cell Small oscillations, stable in time Oscillations averaged to 0 Ab Initio Molecular Dynamic p.48/67
73 Why does it work? Forces oscillations very small Model System : Si FCC, 2 atoms/cell Small oscillations, stable in time Oscillations averaged to 0 Ab Initio Molecular Dynamic p.48/67
74 Working Condition CP : adiabaticity Fast electrons, small nuclei : µ M k ϕ i 0 : wavefunction close to BO, always slightly above. Compromise for µ : t large, but we want electron/nuclei separation Electronic frequency ϕ i occ. and ϕ a virtual : ω e ia = 2(ǫ a ǫ i )/µ Smallest : ω e min E GAP /µ Highest : ω e max E cut /µ so that t max µ/e cut Usually : µ = , t = 5-10 au = fs. Ab Initio Molecular Dynamic p.49/67
75 Loosing adiabaticity If E GAP 0, then ω e min too close to ω nucl For example for an elongated bond, ex Sn 2 Ab Initio Molecular Dynamic p.50/67
76 Loosing adiabaticity If E GAP 0, then ω e min too close to ω nucl For example for an elongated bond, ex Sn 2 For metallic systems! Two solutions : Back to BO Electronic Thermostat Ab Initio Molecular Dynamic p.50/67
77 Comparison BO/CP BOMD Always on BO surface t τ nucl CPMD Always slightly above t τ nucl t BO 5 10 t CP Minimization at each step Only Orthogonalization Problem when deviating Stable with respect to deviations from the BO surface Works for E GAP = 0 Electronic thermostat Used for solid systems Used for liquids Ab Initio Molecular Dynamic p.51/67
78 Comparison BO/CP Ab Initio Molecular Dynamic p.52/67
79 Comparison Ehrenfest/CP Ehrenfest MD Real separation (quantum) t τ elec CPMD Fictitious separation (classical) t τ elec t CP 5 10 t Ehrenfest Rigorous Orthonormality Imposed by constraints Deviations from BO add up Stable Ab Initio Molecular Dynamic p.53/67
80 Conclusion and perspectives AIMD : includes electronic effects, allows for chemical reactions, catalytic processes... CPMD : Very used, very fashion! but... Number of atoms still limited DFT not always sufficient Perspectives Hybrid methods : QM/MM, QM/QM Post-HF and TD approachs : excited states, highly correlated materials... BO? New algorithms faster, more efficient, becomes competitive with CPMD. Ab Initio Molecular Dynamic p.54/67
81 Application : Free energy calculations Ab Initio Molecular Dynamic p.55/67
82 Free energy profile Reaction ({ coordinate q = Rk }) q Probability density P(z) = 1 d Rd P exp( βh)δ Q NV T Free energy profile F(z) = k B T ln P(z) Reaction rate : k = k BT e rg RT h ( q ({ Rk }) ) z Ab Initio Molecular Dynamic p.56/67
83 Rare event? System position for E a = 2 k B T Diffusive system. Ab Initio Molecular Dynamic p.57/67
84 Rare event? System position for E a = 5 k B T Nothing at the TS! Ab Initio Molecular Dynamic p.57/67
85 Rare event? System position for E a k B T Alternating between the two states. Long residence time, short transition time. k TST 1s 1 Slow reaction rare event. Ab Initio Molecular Dynamic p.57/67
86 Why is it so difficult? Standard MD is a "real time" method, t fs Chemical reaction time : fast around ns, typical around µs-ms, biology can be seconds Two incompatible time scales More, phase space dimension is 6N, impossible to fully sample Ab Initio Molecular Dynamic p.58/67
87 Why is it so difficult? Two incompatible time scales Solutions Work at higher T : faster sampling but might be not the same phenomena Force the reaction to occur : Biais potential Constrained dynamic Adiabatic dynamic Metadynamic In all cases, we calculate F, not F Ab Initio Molecular Dynamic p.58/67
88 Bias potential Many approaches : Umbrella sampling, adaptive biais potential, accelerated dynamics, flooding potential... Main Idea : we add a potential to "erase" the barrier ({ V Rk }) ({ Rk }) ( ({ = V + V q Rk })) Then get the non-biaised energy : O ({ Rk }) = O ({ Rk }) exp (β V) exp (β V) Problem : we do not know the barrier position, height and shape! In practice : many simulations with a moving model potential Ab Initio Molecular Dynamic p.59/67
89 Umbrella Sampling Usual bias : V ( q ({ Rk })) [ q ({ Rk }) ] 2 z Simulations for different z Ab Initio Molecular Dynamic p.60/67
90 Constrained dynamics Two families : q evolves continuously from z 0 to z f during the simulation Slow growth method : slow evolution, quasi-static Fast growth method, Jarzinsky : fast evolution, average taken on many trajectories exp ( βw) = exp ( βf) Many simulation with q set to different values of z. Ab Initio Molecular Dynamic p.61/67
91 Thermodynamic integration Change in free energy calculated by integrating free energy derivative : F(z = z 0 z = z f ) = zf z 0 F q dq Numerical evaluation using discrete values A constrained simulation is launched for each value of q(z). How to obtain F q (z)? Ab Initio Molecular Dynamic p.62/67
92 Free energy derivatives Two problems : q = z q = 0 : no sampling for the impulsion p q. Blue-Moon : Z 1/2 O cont O(z) = Z 1/2 cont with Z the reduce mass associated to q : Z = k q is non linear with respect to non unity 1 q q M k R k R k { Rk }, Jacobian matrix J Ab Initio Molecular Dynamic p.63/67
93 Free energy derivatives Two problems : q = z q = 0 : no sampling for the impulsion p q. { Rk } q is non linear with respect to, Jacobian matrix J non unity F q (z) = V q + k BT ln J q q=z Ab Initio Molecular Dynamic p.63/67
94 Free energy derivatives Two problems : q = z q = 0 : no sampling for the impulsion p q. { Rk } q is non linear with respect to, Jacobian matrix J non unity F q (z) = Z 1/2 [ λ + k B T G] Z 1/2 q=z q R k 2 q q R k R l R l, q=z with G = 1 1 Z k,l 2 M k M l λ Lagrange multiplier associated to the constraint σ ({ Rk }) = q z = 0 (SHAKE) Ab Initio Molecular Dynamic p.63/67
95 Adiabatic dynamics Idea : Decoupling motions along q from other motions Use different temperatures One (chain) thermostat associated to q One (chain) thermostat associated to the N-1 other degrees of freedom T q T nuc Change the mass of atoms concerned by q Ab Initio Molecular Dynamic p.64/67
96 Metadynamic Idea : Similar to the adiabatic dynamics, using extended Lagrangian! Reaction coordinates are included as collective variables (CV) s α Coarse MD for CV L = L 0 + α 1 2 M α s α 2 α 1 2 k α [S α (r) s α ] 2 M α large enough natural adiabatic separation Ab Initio Molecular Dynamic p.65/67
97 Metadynamic Idea : Similar to the adiabatic dynamics, using extended Lagrangian! Reaction coordinates are included as collective variables (CV) s α Coarse MD for CV on top of this, potential wells are filled with gaussians V (t) = i W exp ( s s i(t) 2 2δσ 2 Important choices for W, σ, time interval between adding two gaussians ) Ab Initio Molecular Dynamic p.65/67
98 Metadynamic Idea : Similar to the adiabatic dynamics, using extended Lagrangian! Reaction coordinates are included as collective variables (CV) s α Coarse MD for CV Ab Initio Molecular Dynamic p.65/67
99 How to choose a good RC? In practice, difficult to sample a coordinate space larger than 3D Chemist intuition might fail!! Ab Initio Molecular Dynamic p.66/67
100 Conclusion Already many tools but choosing RC is still a hot topic for large systems Classical nuclei : what about quantum effects? Corrections a posteriori : tunneling, isotopic effects Path Integral ZPE correction Ab Initio Molecular Dynamic p.67/67
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 informationMD 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 informationAb 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 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 informationAb 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 informationAb 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 informationWhat 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 informationAb 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 informationAb 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 informationGear 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 informationDensity 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 informationBiomolecular 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 informationA 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 informationSimulation of molecular systems by molecular dynamics
Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationAccelerated 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 informationMolecular 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 informationIV. 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 informationAb 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 informationDensity 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 informationUnderstanding 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 informationTemperature 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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationJavier 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 informationMolecular 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 informationClassical Molecular Dynamics
Classical Molecular Dynamics Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~mijp1 Overview of lecture n Motivation n Types of
More informationSet 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 informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationAb 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 informationFree 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 information7 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 informationComputational 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 informationConstant Pressure Langevin Dynamics: Theory and Application to the Study of Phase Behaviour in Core-Softened Systems.
Constant Pressure Langevin Dynamics: Theory and Application to the Study of Phase Behaviour in Core-Softened Systems. David Quigley A thesis submitted for the degree of Doctor of Philosophy University
More informationStructure 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 informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationScientific 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 informationMolecular dynamics: Car-Parrinello method
Formulations Potential energy Initialization Verlet algorithm MD: Steps MD: Thermo and barostats CP: Car-Parrinello Atomic units MD and CP textbo Molecular dynamics: Car-Parrinello method Víctor Luaña
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationPure and zero-point vibrational corrections to molecular properties
Pure and zero-point vibrational corrections to molecular properties Kenneth Ruud UiT The Arctic University of Norway I UNIVERSITETET TROMSØ June 30 2015 Outline Why consider vibrational effects? General
More informationAnalysis of MD Results Using Statistical Mechanics Methods. Molecular Modeling
Analysis of MD Results Using Statistical Mechanics Methods Ioan Kosztin eckman Institute University of Illinois at Urbana-Champaign Molecular Modeling. Model building. Molecular Dynamics Simulation 3.
More informationMethods of Computer Simulation. Molecular Dynamics and Monte Carlo
Molecular Dynamics Time is of the essence in biological processes therefore how do we understand time-dependent processes at the molecular level? How do we do this experimentally? How do we do this computationally?
More informationStatistical methods in atomistic computer simulations
Statistical methods in atomistic computer simulations Prof. Michele Ceriotti, michele.ceriotti@epfl.ch This course gives an overview of simulation techniques that are useful for the computational modeling
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationCPMD Tutorial Atosim/RFCT 2009/10
These exercices were inspired by the CPMD Tutorial of Axel Kohlmeyer http://www.theochem.ruhruni-bochum.de/ axel.kohlmeyer/cpmd-tutor/index.html and by other tutorials. Here is a summary of what we will
More informationLangevin 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 informationChE 210B: Advanced Topics in Equilibrium Statistical Mechanics
ChE 210B: Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson Lecture 1 Reading: 3.1-3.5 Chandler, Chapters 1 and 2 McQuarrie This course builds on the elementary concepts of statistical
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationTemperature 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 informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationCitation 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 informationPart 2: Molecular Dynamics. Literature History Practical Issues
Part 2: Molecular Dynamics Literature History Practical Issues 1 Literature A. K. Hartmann Big Practical Guide to Computer Simulations (World Scientific, 2015) M. P. Allen and D.J. Tildesley Computer Simulations
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationEnabling 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 informationBorn-Oppenheimer Approximation
Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,
More informationConical 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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationAdvanced in silico drug design
Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms
More informationQuantum Molecular Dynamics Basics
Quantum Molecular Dynamics Basics Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationIntroduction 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 informationLangevin Dynamics in Constant Pressure Extended Systems
Langevin Dynamics in Constant Pressure Extended Systems D. Quigley and M.I.J. Probert CMMP 2004 1 Talk Outline Molecular dynamics and ensembles. Existing methods for sampling at NPT. Langevin dynamics
More informationTime reversible Born Oppenheimer molecular dynamics
Time reversible Born Oppenheimer molecular dynamics Jianfeng Lu Mathematics Department Department of Physics Duke University jianfeng@math.duke.edu KI-Net Conference, CSCAMM, University of Maryland, May
More informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationIntroduction to Geometry Optimization. Computational Chemistry lab 2009
Introduction to Geometry Optimization Computational Chemistry lab 9 Determination of the molecule configuration H H Diatomic molecule determine the interatomic distance H O H Triatomic molecule determine
More informationEnhanced 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 informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationMixed quantum-classical dynamics. Maurizio Persico. Università di Pisa Dipartimento di Chimica e Chimica Industriale
Mixed quantum-classical dynamics. Maurizio Persico Università di Pisa Dipartimento di Chimica e Chimica Industriale Outline of this talk. The nuclear coordinates as parameters in the time-dependent Schroedinger
More informationQUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES
QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES 1 Tapio T. Rantala, Department of Physics, Tampere University of Technology http://www.tut.fi/semiphys CONTENTS MOTIVATION PATH INTEGRAL APPROACH
More informationModeling & Simulation of Glass Structure
Modeling & Simulation of Glass Structure VCG Lecture 21 John Kieffer Department of Materials Science and Engineering University of Michigan 1 Overview Historical perspective Simulation methodologies Theoretical
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More information5. Systems in contact with a thermal bath
5. Systems in contact with a thermal bath So far, isolated systems (micro-canonical methods) 5.1 Constant number of particles:kittel&kroemer Chap. 3 Boltzmann factor Partition function (canonical methods)
More informationStatistical 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 informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationTHERMOSTATS AND THERMAL TRANSPORT IN SOLIDS
Hands-On Tutorial Workshop, July 19 th 2011 THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS Christian Carbogno FRITZ-HABER-INSTITUT DER MAX-PLANCK-GESELLSCHAFT, BERLIN - GERMANY Hands-On Tutorial Workshop,
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationStatistical thermodynamics (mechanics)
Statistical thermodynamics mechanics) 1/15 Macroskopic quantities are a consequence of averaged behavior of many particles [tchem/simplyn.sh] 2/15 Pressure of ideal gas from kinetic theory I Molecule =
More informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationRandom Walks A&T and F&S 3.1.2
Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will
More informationFree Energy Simulation Methods
Free Energy Simulation Methods Free energy simulation methods Many methods have been developed to compute (relative) free energies on the basis of statistical mechanics Free energy perturbation Thermodynamic
More informationBiomolecular modeling I
2016, December 6 Biomolecular structure Structural elements of life Biomolecules proteins, nucleic acids, lipids, carbohydrates... Biomolecular structure Biomolecules biomolecular complexes aggregates...
More informationThermal and Mechanical Properties of Pt-Rh Alloys
arxiv:cond-mat/9611241v1 [cond-mat.mtrl-sci] 28 Nov 1996 Thermal and Mechanical Properties of Pt-Rh Alloys G. Dereli, T. Çağın, M. Uludoğan, M. Tomak Department of Physics, Middle East Technical University,
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationCopyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More information4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D.
4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro Karel Berka, Ph.D. UP Olomouc, 21.1.-25.1. 2019 Motto A theory is something nobody believes, except the person who made it An experiment
More informationIntroduction to model potential Molecular Dynamics A3hourcourseatICTP
Introduction to model potential Molecular Dynamics A3hourcourseatICTP Alessandro Mattoni 1 1 Istituto Officina dei Materiali CNR-IOM Unità di Cagliari SLACS ICTP School on numerical methods for energy,
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
More informationIMPLICIT MASS-MATRIX PENALIZATION OF HAMILTONIAN DYNAMICS WITH APPLICATION TO EXACT SAMPLING OF STIFF SYSTEMS
IMPLICIT MASS-MATRIX PENALIZATION OF HAMILTONIAN DYNAMICS WITH APPLICATION TO EXACT SAMPLING OF STIFF SYSTEMS PETR PLECHÁČ AND MATHIAS ROUSSET Résumé. An implicit mass-matrix penalization (IMMP) of Hamiltonian
More informationMolecular Modeling of Matter
Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More information