Molecular modeling for functional materials design in chemical engineering and more
|
|
- Magdalene Cox
- 5 years ago
- Views:
Transcription
1 Molecular modeling for functional materials design in chemical engineering and more Dr inż. Łukasz Radosiński Group of Bioprocess and Biomedical Engineering Wroclaw University of Science and Technology
2 Nanomaterials in technology Improved chemical resistance, Improved hardness and abrasion resistance, Increased tensile strength and flexibility, Favorable melting points, Favorable magnetic properties, Increased thermal and electrical conductivity, and Favorable surface-chemistry effects that improve the ability of a powder to be dispersed in a liquid.
3 CAD is now used in every part of industry Civil engineering, Chemical Engineering, Areonautics, Medical and pharmaceutical
4 Computer Aided Design Growing widely available computational power changed the product development process. We are able to predict properties of the system long before it was even created Golden management rule: if you want to optimize you have have numbers, if you have numbers you have to know how to calculate it, if you want to calculate you have to have tools computers + physics + chemistry + mathematics
5 Multi-scale view of the materials Buehler and Ackbarow, Materials Today, 2007
6 Molecular modelling Length [m] 10 3 DPD 10-3 Molecular Dynamics, Monte Carlo Coarse-Grained CFD 10-9 Ab Initio Time [s]
7 Ab initio methods Predicting of physical and chemical properties of materials from quantum theory Very accurate Without any empirical parameters Calculating parameters of forcefields for Molecular Dynamics and Monte Carlo methods Electronic properties of matter (charges, reactivity, dipole moments, IR spectrum and much more ) Very accurate geometrical parameters
8 Ab initio methods Size of a system is very limited (up to hundred atoms) Long computing time for even small systems We need high computing power Any movement of atoms or timedependent methods need much more computing time and power Theory is not so easy
9 Molecular dynamics Millions of atoms thousands of cores Solving of Newton equation of motion for all atoms Bigger systems (thousands of atoms) Phenomena depending on time (diffusion, cracking of materials, melting, solvation and so on ) Introduction of a temperature and pressure to the system
10 Molecular dynamics Size is still limited (up to billion of atoms) Simulated time is not enough long (bio-physical or chemical processes may take seconds) Everything depends on forcefield! Disregarding of quantum nature of matter (single atom has no energy no interaction with anything) Very large systems need high computing power and time
11 Monte Carlo method Averages over time (MD) are replaced by averages over states of a system States of a system are selected randomly
12 Monte Carlo method Cannot simulate time dependent properties Size is still limited Everything depends on forcefield! How to properly generate random state of a system, is it really random?
13 Coarse-Grained method Similar method to molecular dynamics - groups of atoms are replaced with superatoms Possibility to simulate very large systems Need to describe interactions between superatoms Disregarding of properties of separate atoms Everything depends on forcefield!
14 Macroscale Finite elements methods: CFD, FDM itp. tensions, flow of mass, temperature, etc.
15 When do we need molecular modeling? Beam deformation problem Question: displacement function d z x Governing equation: EI zz 4 d z,0 x 4 + q z = 0 Geometry: I zz = bh3 12 Boundary condition: q z = ρga Integration d z x = ρga 24EI zz x 4 d z x x How to obtain Young s modulus E?
16 How to obtain Young s modulus - experiment Laboratory measurements Molecular modeling A A L L + Δl F F = kδl E = Lk A
17 What is our goal? - To gain extensive knowledge about state-of-the-art methods of predicting and deisigning properties of material using molecular approach, - To be able to predict basic properties of the various physical systems and conntect the results with macroscopic values
18 Our way bottom up
19 Are Molecular Simulations part of Chemical Engineering?
20 Science
21 Science
22 Numerical methods numerical experiment
23 Can we simulate world?
24 Can we simulate world?
25 Can we simulate world?
26 Can we simulate world?
27 Why does it look real????
28 Why does it look real????
29 Why does it look real???? Newtonian motion!
30 Newton s Equation a F m
31 Newton s equation F dv F a v m dt m v dr dt r(t) Knowing force, initial velocities and position we can extract velocity and position in any moment of time!
32 How molecular dynamics works? V 1 1 (x1 p, y1 p,z1 p ) 2 (x2 p, y2 p,z2 p ) V 2
33 How molecular dynamics works? V 1 1 (x1 p, y1 p,z1 p ) F 12 (x1 p, y1 p,z1 p, x2 p, y2 p,z2 p ) F 21 (x1 p, y1 p,z1 p, x2 p, y2 p,z2 p ) 2 (x2 p, y2 p,z2 p ) V 2
34 How molecular dynamics works? V 1 1 (x1 p, y1 p,z1 p ) F 12 /m 1 =a 12 a = Dv Dt F 21 /m 2 =a 21 2 (x2 p, y2 p,z2 p ) V 2
35 How molecular dynamics works? ΔV 1 V 1 (t) 1 (x1 p, y1 p,z1 p ) F 12 /m 1 =a 12 a = Dv Dt V 1 (t+δt) Dv = adt V 2 (t+δt) F 21 /m 2 =a 21 2 ΔV 2 Dv = F m Dt (x2 p, y2 p,z2 p ) V 2 (t)
36 How molecular dynamics works? (x1 n, y1 n,z1 n ) 1 Dr 1 1 (x1 p, y1 p,z1 p ) Dr = v(t)dt V 1 (t+δt) (x2 p, y2 p,z2 p ) 2 Dr 2 2 (x2 n, y2 n,z2 n ) V 2 (t+δt)
37 How molecular dynamics works? (x1 n, y1 n,z1 n ) 1 V 1 (t+δt) F 21 (t + Dt) 2 V 2 (t+δt) (x2 n, y2 n,z2 n )
38 Mathematically Euler s scheme df f (t + Dt) - f (t) f (t + Dt) - f (t) dt = lim Dt 0 Dt» Dt df dt = f i+1 - f i Dt
39 Instead of solving function analytically step by step procedure f t + Δt = f t + f t Δt f Analytical solution t
40 Instead of solving function analytically step by step procedure f t + Δt = f t + f t Δt f Analytical solution f Point-by-point solution t t
41 How does it works we start from f(0) 1. Choose step Δt = 1 2. Compute next point using chain rule f t + Δt = f t + f t Δt Example f t = 2 t f t = 2
42 How does it works we start from f(0) f 0 + Δt = 1 = = 2 f f(1) 1 t
43 How does it works we start from f(0) f = = 4 f f(2) f(1) t 2 2
44 And so on point by point we reproduce solution f f(t) = 2t f(2) f(1) t 2 2
45 Problems? - Accuracy Example derivative varies (change of the function s shape) f Wrong value! t
46 Solution? More points, shorter interval Example derivative varies (change of the function s shape) f Problem? Computing efficiency! t
47 Application Solving Newton s Eq. 1. Assuming we know the Force acting on a body we write down Newton s Eq.: a = dv dt = F m
48 Application Solving Newton s Eq. 1. Assuming we know the Force acting on a body we write down Newton s Eq.: a = dv dt = F m 2. We use Euler s scheme to rewrite differentials: dv dt v t + Δt v t Δt
49 Application Solving Newton s Eq. 3. We solve our equation of motion v t + Δt v t Δt = F m v t + Δt = v t + F m Δt
50 Application Solving Newton s Eq. 3. Knowing initial condition say v(0) and setting Δt we reproduce point by point evolution of velocity using our finite difference scheme: v Δt = v 0 + F m Δt
51 Application Solving Newton s Eq. 3. Knowing initial condition say v(0) and setting Δt we reproduce point by point evolution of velocity using our finite difference scheme: v Δt = v 0 + F m Δt v 2Δt = v Δt + F m Δt
52 Application Solving Newton s Eq. 3. Knowing initial condition say v(0) and setting Δt we can reproduce point by point evolution of velocity using our finite difference scheme: v Δt = v 0 + F m Δt v 2Δt = v Δt + F m Δt v 0, v Δt, v 2Δt, v 3Δt,
53 What about position? Same scheme! 1. Having collection of v(t) we use definion of velocity as: v(t) = Δx Δt
54 What about position? Same scheme! 1. Having collection of v(t) we use definion of velocity as: v(t) = Δx Δt x t + Δt x t Δt = v(t)
55 What about position? Same scheme! 1. Having collection of v(t) we use definion of velocity as: v(t) = Δx Δt x t + Δt x t Δt = v(t) x t + Δt = x t + v t Δt
56 What about position? Same scheme! 1. Having collection of v(t) we use definion of velocity as: v(t) = Δx Δt x t + Δt x t Δt = v(t) x t + Δt = x t + v t Δt We already have this!
57 What about position? Same scheme! 1. Having collection of v(t) we use definion of velocity as: v(t) = Δx Δt x t + Δt x t Δt = v(t) x t + Δt = x t + v t Δt We already have this! x 0, x Δt, x 2Δt,
58 POWEER! Stupid pig.. Angry Birds
59 Newton`s equation
60
61 Newton`s equation algorithm
62 Newton`s equation algorithm
63 Newton`s equation algorithm example
64 Newton`s equation algorithm example
65 Harmonic motion why so important?
66 Harmonic motion why so important?
67 Harmonic motion why so important? Equilibrium length no motion Out of equilibrium - F = kx, E p = kx2 2 Hook s law Holds up to certain distance
68 Harmonic motion why so important? Equilibrium length no motion Out of equilibrium - F = kx, E p = kx2 2 Hook s law Holds up to certain distance
69 Harmonic motion why so important? Equilibrium length no motion Out of equilibrium - F = kx, E p = kx2 2 Hook s law Holds up to certain distance Equilibrium shape no motion Deformation F = EA ΔL, E L p = 0 Young s modulus Holds up to certain def. scale EA L0 ΔL2 2
70 Harmonic motion why so important? Equilibrium length no motion Out of equilibrium - F = kx, E p = kx2 2 Hook s law Holds up to certain distance Equilibrium shape no motion Deformation F = EA ΔL, E L p = 0 Young s modulus Holds up to certain def. scale EA L0 ΔL2 2
71 Harmonic motion why so important? Equilibrium length no motion Out of equilibrium - F = kx, E p = kx2 2 Hook s law Holds up to certain distance Equilibrium shape no motion Deformation F = EA ΔL, E L p = 0 Young s modulus Holds up to certain def. scale EA L0 ΔL2 2 Bond length no motion (zero temp., no quantum motion) Out of eq. F~Δl, E p ~ Δl2 2 Quantum mechanics harmonic approx.
72 Harmonic motion why so important? Why so many similarities? Physics??? Equilibrium length no motion Out of equilibrium - F = kx, E p = kx2 2 Hook s law Holds up to certain distance Equilibrium shape no motion Deformation F = EA ΔL, E L p = 0 Young s modulus Holds up to certain def. scale EA L0 ΔL2 2 Bond length no motion (zero temp., no quantum motion) Out of eq. F~Δl, E p ~ Δl2 2 Quantum mechanics harmonic approx.
73 Harmonic motion why so important? NO! Mathematics!
74 Two balls in equilibrium no motion E E x x
75 Two balls in equilibrium no motion E E No motion x x
76 Two balls in equilibrium no motion E E x No motion No forces x
77 Two balls in equilibrium no motion E E x No motion No forces Energy minimum x
78 Energy change due to movement E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0
79 Can we calculate energy E(Δx)? E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0
80 Can we calculate energy E(Δx)? E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0 de 0 E 0 + Δx = E 0 + dx Δx + 1 d 2 E 0 2 dx 2 Δx 2 +
81 Can we calculate energy E(Δx)? Minimum de 0 dx = 0! Minimum de 0 dx = 0! E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0 de 0 E 0 + Δx = E 0 + dx Δx + 1 d 2 E 0 2 dx 2 Δx 2 +
82 When Δx is small Minimum de 0 dx = 0! Minimum de 0 dx = 0! E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0 E Δx 1 2 k d 2 E 0 dx 2 Δx 2
83 When Δx is small Minimum de 0 dx = 0! Minimum de 0 dx = 0! E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0 E Δx 1 2 k d 2 E 0 dx 2 Δx 2 = kδx2 2
84 What about force? Minimum de 0 dx = 0! F = kx F = kx Minimum de 0 dx = 0! E E E 0 + Δx E 0 + Δx x E 0 = 0 x Δx Δx E 0 = 0 F = de dx = d kx2 2 dt = kx
85 Harmonic motion why so important? Chemical bonds! And every coupling with eq. distance
86 But there is sth extra! F = kx x = Asin ω 0 t ω 0 = k m
87 But there is sth extra! F = kx F ext = Asin(ωt)
88 But there is sth extra! F = kx F ext = Bsin(ωt) ma = Bsin ωt kx A = m B ω 0 ω 2 When ω ω 0 A is huge! Great response! Resonance
89 Where is chemistry? Light! F = kx F ext = qesin(ωt) Photon! When frequency = energy of photon is close to normal mode there is huge transfer of energy Adsorption!
90 And this we can calculate! O C O
91 Molecular modelling Length [m] 10 3 DPD 10-3 Molecular Dynamics, Monte Carlo Coarse-Grained CFD 10-9 Ab Initio Time [s]
92 Ab initio methods Predicting of physical and chemical properties of materials from quantum theory Very accurate Without any empirical parameters Calculating parameters of forcefields for Molecular Dynamics and Monte Carlo methods Electronic properties of matter (charges, reactivity, dipole moments, IR spectrum and much more ) Very accurate geometrical parameters
93 Ab initio methods Size of a system is very limited (up to hundred atoms) Long computing time for even small systems We need high computing power Any movement of atoms or timedependent methods need much more computing time and power Theory is not so easy
94 Molecular dynamics Millions of atoms thousands of cores Solving of Newton equation of motion for all atoms Bigger systems (thousands of atoms) Phenomena depending on time (diffusion, cracking of materials, melting, solvation and so on ) Introduction of a temperature and pressure to the system
95 Molecular dynamics Size is still limited (up to billion of atoms) Simulated time is not enough long (bio-physical or chemical processes may take seconds) Everything depends on forcefield! Disregarding of quantum nature of matter (single atom has no energy no interaction with anything) Very large systems need high computing power and time
96 Monte Carlo method Averages over time (MD) are replaced by averages over states of a system States of a system are selected randomly
97 Monte Carlo method Cannot simulate time dependent properties Size is still limited Everything depends on forcefield! How to properly generate random state of a system, is it really random?
98 Coarse-Grained method Similar method to molecular dynamics - groups of atoms are replaced with superatoms Possibility to simulate very large systems Need to describe interactions between superatoms Disregarding of properties of separate atoms Everything depends on forcefield!
99 Our way bottom up
100 What do we obtain? Trajectories Position and velocity in time
101 Molecular Dynamics Tasks Fundamental: to compute trajectories Minimazation of energy geometrical analysis of compouds To compute avarages and current physical quantities in in the investigated systems (kinetic energy, potential energy, interatomic distance, pressure, volume, atoms spatial correlation, diffusion, heat transfer, viscosity etc.)
102 Molecular Dynamics Tasks Fundamental: to compute trajectories Minimazation of energy geometrical analysis of compouds To compute avarages and current physical quantities in in the investigated systems (kinetic energy, potential energy, interatomic distance, pressure, volume, atoms spatial correlation, diffusion, heat transfer, viscosity etc.)
103 What is missing? 1. Intermolecular forces. 2. Initial positions. 3. Initial velocities.
104 What is missing? 1. Forces forcefields (quantum chemical computations), empirical potentials. 2. Initial positions crystal structure, fluids order, gas density etc. 3. Initial velocities zero, random, temperature.
105 Initial positions Initial positions of atoms should be as close to equilibrium positions (minimal energy) as possible scientific knowledge To big deviation from equilibrium positions results in unphysical behaviour large forces
106 Temperature How to introduce temperature to the system? Simplest solution is to draw velocities from Boltzmann distribution for fixed T:
107 Temperature of the system is its average kinetic energy
108 Temperature is not constant Unperfect solution of the Verlet s algorithm error cumulation
109 Problems? What is the average fluctuation in microcanonical ensamble?
110 Using Boltzmann distribution
111 Solution? Thermostat (heat bath) Typical MD simulation is in NVE regime -> transformation to NVT and NPT ensambles only in the thermodynamical limit (inifnite system size) Introduction of the heat bath separate system with fixed temperature System T(t) Heat bath T 0
112 Beredsen thermostat Weak coupling via Newton s law Coupling constant Reservouir s temperature
113 How to apply? Rescale velocities by step dependant rescalling factor
114 How to obtain lambda?
115 Beredsen s thermostat Pros: Straight forward code Robust Exponential decay to the desired temperature Cons: Results do not correspond to any ensamble, however small deviation from canonical only during equilibriation No-time reversible
116 Noose-Hoover Thermostat Introducing friction to the heat reservouir by including new variable to the equation of motion Kinetic energy of the hb Potential energy of the hb Fictious mass of the hb
117 Nose-Hoover scheme Apply Hamilton s equation of motion, Obtain coupled system of diff. eq. for r and s with coupling constant Q Solve dynamics
118 Nose-Hoover thermostat Pros: Gives canonical properties Deterministic, time reversible Conserves avarages p A, r A p, s r c Cons: In some cases extended system not ergodic entrapment in subspace slow Not implemented in NAMD!
119 Langevin s thermostat Introducing hypothetical fluid of much smaller particles each system s particle is subjected to a drag and a stochastic force Force Drag Noise
120 The drag and noise are balanced in order to obtain constant temperature Gamma is a friction coefficiant A is a random force uncorrelated in time and across particles
121 Langevin thermostat Pros: Reproduces canonical properties Ergodic Allows larger time steps as compared to non-stochastic thermostats Cons: Undeterministic Momentum trasfer not conserved not advisable to use in the system with diffusion Temperature fluctuations
122
123 Ekwilibracja Potential Energy as a Function of Time 0 Potential Energy (kcal/mol) time (ps) 123
124 Fluctuations in Total Energy for an NVE Simulation of 216 TIP3P Water Molecules (1fs time step) [the run was initiated at E = kcal/mol] -1730,2-1730,25 Total Energy (kcal/mol) -1730,3-1730, ,4-1730, , time (ps) 124
125 Thermostating & Barostating Simulation system is dynamical. In order to preserve the system specific temperature or pressure we employ in the calculation the concept of thermostat and barostat. Conceptionaly these are separate systems with defined pressure an temperature connected to our simulation balls. Mathematicaly we simply add two Newtons equations, additionaly equations defining temperature and pressure.
126 Pressure virial theorem Virial theorem: Pressure in the system:
127 Berendsen Thermostat The velocities are scaled at each step, such that the rate of change of temperature is proportional to the difference in temperature. Scaling factor:
128 Berendsen Barostat. In barostat the mechanism is similar: Scaling factor:
129 NPT V nt/p During simulation ensemble uses barostat and thermostat to equilibrate temperature and pressure. Varying, operating both: P and T results in changes of volume of the simulation box. In other words, our goal is to reach density of water about 1.0 g/cc, that is a property of our system in 298 K and 1 bar.
130 Dynamics Calculation stages to define temperature and density NPT NVT NVE equilibration of density (volume changes) equilibration of temperature final dynamics
131 NPT Box length running average
132 NPT Density
133 NVT High kinetic energy equilibrated temperature Running average Random Velocities High potential Energy Temperature Strong thermostat work
134 Real systems simulations Problems and solutions
135 Problems of simulating real systems Most if the systems contains huge number of atoms >> Typical atomic time scale is10-15 s human is living in seconds Most of interactions are nonbonded (electrostatic, dispersion etc.) infinite range of interactions
136 Periodic Boundary Conditions
137 Periodic Boundary Conditions Pros: Elimination of the boundary problem Cons: Neutral charge of the system requirement unbalanced charge creates multiple images due to PBC resulting unphysical behaviour Angular momentum in not conserved - no rotational symmetry Cannot observe oscillation of wavelength longer then PBC box gas-liquid critical point
138 Infinite system
139 Simplication to a single simulation box
140 Periodic Boundary Conditions examples Spherical Cubical
141 Goal: glucose in water no PBC Glucose PBC addition Water addition
142 Goal: glucose in water no PBC No problems with non-bonded interactions Huge number of atoms Expensive PBC Atom may interact with itself Small number of atoms Cheaper
143 Geometry optimization Each system must be optimized before your proper simulation. Every algorithm leads potential energy function to its minimum. You must be aweare that potential energy is rather multidimensional hypersurface than 3D function (seen on the pic.). Initial energy Global or local minimum
144 Geometry optimization Example of convergence criteria: 1. Convergence of the function (energy) f x k+1 f(x k ) < ε f 2. Convergence of the variable (position) x k+1 x k < ε x 3. Covergence of the gradient (force) f x k+1 f(x k ) < ε g Initial energy Global or local minimum
145 Geometry optimization Examples of most popular energy minimization algorithms: - conjugate gradients, - steepest descent. And also some damping dynamics algorithms such as: - quickmin, - fire. Initial energy Global or local minimum
146 Potential energy Geometry optimization 60º 109.5º 145º 180º 215º 215.5º 300º Angle
147 Problems of simulating real systems The results are supposed have relevance to real systems values Becouse time scale of nanoscopic systems is fs and ours time scale is s we measure AVARAGES temperature, pressure, energy, specific heat etc. Question: how long do we have to simulate the process so the avarages obtained in the simulation are thermodynamical avarages? Ergodic hypothesis.
148 Molecular dynamics protocol Initialization (positions) Energy Minimanizatiob (excess potential energy removal) Initialization of velocities (temperature) Equilibration ( smearing of kinetic energyj) Simulation Trajectory Analysis
149 Molecular dynamics protocol Initialization (positions) Energy Minimanizatiob (excess potential energy removal) Initialization of velocities (temperature) Equilibration ( smearing of kinetic energyj) Simulation Trajectory Analysis
150 Molecular dynamics protocol Initialization (positions) Energy Minimanizatiob (excess potential energy removal) Initialization of velocities (temperature) Equilibration ( smearing of kinetic energyj) Simulation Trajectory Analysis
151 Molecular dynamics protocol Initialization (positions) Energy Minimanizatiob (excess potential energy removal) Initialization of velocities (temperature) Equilibration ( smearing of kinetic energy) Simulation Trajectory Analysis
152 Molecular dynamics protocol Initialization (positions) Energy Minimanizatiob (excess potential energy removal) Initialization of velocities (temperature) Equilibration ( smearing of kinetic energy) Simulation Trajectory Analysis
153 Molecular Dynamics Forces
154 Protokół dynamiki molekularnej
155 Types of interactions Bonded interactions: Strong, Short-ranged, Typically up to the second nearest neighbour, Strongly depends on the geometry and electronic structure. Can be quantified using quantum mechanics.
156 Types of interactions nonbonded interactions: Long-ranged, Every atoms interacts with every atom in the system Weak (van der Waalsa) or strong (electrostatic)
157 Forcefield Forcefield (AMBER, CHARMM), Empirical (AIREBO, LCBOP), Quantum based (COMPASS)
158 Potential function V = V bonds +V angle +V torsion +V improper +V nonbonded
159 Funkcja potencjału empirycznego V = V bonds +V angle +V torsion +V improper +V nonbonded ( ) 2 å + k q q -q 0 V = k b b - b 0 bonds N å n=1 ( ) + K j n angles ( ) 2 å + å 1+ cos( nj -d ) + K w w - w 0 dihedrals + 4e ij s ij éë 12 æ ö ç è r - s 6 é æ ö ù å ij ê ç ij ø è r ú + i, j ë ê ij ø û ú ù û å i, j æ ç è å impropers q i q j Dr ij ö ø ( ) 2
160 Bonded interactions V = å k ( b r - r ) 2 0 bonds
161 Chemical bond Morse potential E V bonds = D e (1- exp(-a(l - l 0 ))) 2 a = w n / 2D e distance
162 Chemical bond Morse potential E V bonds = D e (1- exp(-a(l - l 0 ))) 2 a = w n / 2D e distance Pauli s exclusion principle
163 Chemical bond Morse potential E V bonds = D e (1- exp(-a(l - l 0 ))) 2 a = w n / 2D e Eq. distance distance Pauli s exclusion principle
164 Chemical bond Morse potential E V bonds = D e (1- exp(-a(l - l 0 ))) 2 a = w n / 2D e Eq. distance distance Pauli s exclusion principle Bond brakage
165 Chemical bond Morse potential E Harmonic potential V bonds = D e (1- exp(-a(l - l 0 ))) 2 a = w n / 2D e Pauli s exclusion principle Eq. distance distance Bond brakage V bonds = k 2 (l - l 0) 2 Computational simplicity
166 Bonded interaction bond angle V angle = k 2 (q -q 0) 2
167 Bonded interaction torsion angle
168 Bonded interaction torsion angle V angle = N å V n 2 n=0 (1+ cos(nw -g )) Most popular model MM2 V = ½ V tor,1 (1 cos ω ) + ½ V tor,2 (1 - cos 2 ω ) + ½ V tor,3 ( 1 - cos 3 ω ) Dipol-dipol interactions, differences in electronegativity Conjugations (alkanes) double bond Sterical effects Interactions between orbitals
169 Bonded interaction torsion angle V angle = N å V n 2 n=0 (1+ cos(nw -g ))
170 Bonded interaction improper torsion V angle = k(1- cos2w) k i l j ω Thomas W. Shattuck
171 Bonded interaction improper torsion V angle = k(1- cos2w) k i l j ω Thomas W. Shattuck Induce planar structure Example: cyclobutan
172 Bonded interaction conjugates ω ω ω Bond - bond Bond - angle V cross = k 2 (l 1 - l 10 )(l 2 - l 2 0 ) V cross = k 2 (l 1 - l 10 )(w -w 0 ) Angle - angle V cross = k 2 (w 1 -w 10 )(w 2 -w 20 )
173 Van der Waals interactions å i, j 4e ij éæ ê ç ê ë è r m r ij ö ø 12 æ - r m ç è r ij ö ø 6 ù ú ú û dipol dipol interaction type: Weak, long-ranged 1 r ij 6 dipol dipol interaction 1 r ij 12 Pauli exclusion principle two particles cannot occupy one point in space
174 Coulomb interactions
175 Coulomb interactions Each atom interacts with every atom in the system computational complexity O(n 2 )
176 Non bonded interactions ( ) 2 V = å k b b - b 0 + å k q q -q 0 + bonds N angles ( ) 2 ( n) + å K j é ë 1+ cos( nj -d ) ù û + K w w - w 0 dihedrals å n=1 12 æ s + 4e ij ö ij ç è r - s 6 é æ ö ù ij å ê ç ij ø è r ú + i, j ê ë ij ø ú û å i, j å impropers æ q i q j ö ç è ø Dr ij ( ) 2
177 Non bonded interactions ( ) 2 V = å k b b - b 0 + å k q q -q 0 + bonds N angles ( ) 2 ( n) + å K j é ë 1+ cos( nj -d ) ù û + K w w - w 0 dihedrals å n=1 12 æ s + 4e ij ö ij ç è r - s 6 é æ ö ù ij å ê ç ij ø è r ú + i, j ê ë ij ø ú û å i, j å impropers æ q i q j ö ç è ø Dr ij ( ) 2 Greatest computational complexity Every atoms interacts with every atom
178 Non bonded interactions ( ) 2 V = å k b b - b 0 + å k q q -q 0 + bonds N angles ( ) 2 ( n) + å K j é ë 1+ cos( nj -d ) ù û + K w w - w 0 dihedrals å n=1 12 æ s + 4e ij ö ij ç è r - s 6 é æ ö ù ij å ê ç ij ø è r ú + i, j ê ë ij ø ú û å i, j å impropers æ q i q j ö ç è ø Dr ij ( ) 2 Greatest computational complexity Every atoms interacts with every atom 1500 atoms
179 Non bonded interactions ( ) 2 V = å k b b - b 0 + å k q q -q 0 + bonds N angles ( ) 2 ( n) + å K j é ë 1+ cos( nj -d ) ù û + K w w - w 0 dihedrals å n=1 12 æ s + 4e ij ö ij ç è r - s 6 é æ ö ù ij å ê ç ij ø è r ú + i, j ê ë ij ø ú û å i, j å impropers æ q i q j ö ç è ø Dr ij ( ) 2 Greatest computational complexity Every atoms interacts with every atom 1500 atoms force pairs
180 Non bonded interactions Solution: Cut-off distance Ewald Summation ( ) 2 V = å k b b - b 0 + å k q q -q 0 + bonds N angles ( ) 2 ( n) + å K j é ë 1+ cos( nj -d ) ù û + K w w - w 0 dihedrals å n=1 12 æ s + 4e ij ö ij ç è r - s 6 é æ ö ù ij å ê ç ij ø è r ú + i, j ê ë ij ø ú û å i, j å impropers æ q i q j ö ç è ø Dr ij Greatest computational complexity Every atoms interacts with every atom ( ) atoms force pairs
181 Non bonded interactions cut off distance For atom A we calculate non-bonded interactions only in the sphere of radius δ Verifaction process required
182 Ewald summation method Ewald summation is an algorithm wich effciently calculates the interactions between periodic mirrors of atoms or molecules. Van der Waals interactions deacys like 1 r6 so there is no porblem usually periodic mirrors are over the cutoff. In case of Coulomb electrostatic interactions k q iq i_mirror potential r i i_mirror deacys slow. Ewald proposed 1 + which states for short range, rapidly r r r varying function and long range flat function. The f r = erfc r = 2 exp t 2 dt. π x = f r 1 f r
183 Ewald summation method Reciprocal space Real space r r Each charge is neutralised with neutralising charge Another contribution counteracts neutralisation Real space r
184 Particle-particle particle-mesh In the particle-particle particlemesh method (PPPM, P 3 M) some simplifications are made. Beyond some distance R c particles (atoms) are discetized charge density gives a contribution to neighbouring grid nodes. For close particels (atoms) interactions are calcualted by common Ewald summation. Mesh grid
185 Ewald vs. PPPM (P3M) Ewald summation Particle Particle Particle mesh Summation over every atom Summation over a discretized grid of charge denisty Uses Fourier transform Uses fast Fourier transform Faster for smaller systems Significantly faster for large systems (N > 10000) # Atoms Ewald (c. time) P3M (c. time)
186 What about large molecules? atoms
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 informationFUNCTIONALISATION, HYDRATION & DYNAMICS
FUNCTIONALISATION, HYDRATION & DYNAMICS MATERIALS STUDIO TUTORIAL WAMMBAT - DAY 4 Our main goal: MD simulation of glucose in hydrated & functionalised silica nanopore It is not a rocket science and you
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 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 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 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 informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
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 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 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 informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationT6.2 Molecular Mechanics
T6.2 Molecular Mechanics We have seen that Benson group additivities are capable of giving heats of formation of molecules with accuracies comparable to those of the best ab initio procedures. However,
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 informationBrownian Motion and Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
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 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 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 informationBiomolecules are dynamic no single structure is a perfect model
Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",
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 informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
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 informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More informationSimulations with MM Force Fields. Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi
Simulations with MM Force Fields Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi Some slides taken with permission from Howard R. Mayne Department of Chemistry University of New Hampshire Walking
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 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 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 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 informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
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 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 informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
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 informationThe micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction
ISBN 978-1-84626-081-0 Proceedings of the 2010 International Conference on Application of Mathematics and Physics Volume 1: Advances on Space Weather, Meteorology and Applied Physics Nanjing, P. R. China,
More informationPotential Energy (hyper)surface
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 informationMolecular Mechanics. Yohann Moreau. November 26, 2015
Molecular Mechanics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Mechanics, Label RFCT 2015 November 26, 2015 1 / 29 Introduction A so-called Force-Field
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
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 informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
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 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 informationInterface Resistance and Thermal Transport in Nano-Confined Liquids
1 Interface Resistance and Thermal Transport in Nano-Confined Liquids Murat Barisik and Ali Beskok CONTENTS 1.1 Introduction...1 1.2 Onset of Continuum Behavior...2 1.3 Boundary Treatment Effects on Interface
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 informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
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 informationPrinciples and Applications of Molecular Dynamics Simulations with NAMD
Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology
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 informationModeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA
HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface
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 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 informationIntermolecular Forces and Monte-Carlo Integration 열역학특수연구
Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 2003.3.28 Source of the lecture note. J.M.Prausnitz and others, Molecular Thermodynamics of Fluid Phase Equiliria Atkins, Physical Chemistry Lecture
More informationVariable Definition Notes & comments
Extended base dimension system Pi-theorem (also definition of physical quantities, ) Physical similarity Physical similarity means that all Pi-parameters are equal Galileo-number (solid mechanics) Reynolds
More informationShear Properties and Wrinkling Behaviors of Finite Sized Graphene
Shear Properties and Wrinkling Behaviors of Finite Sized Graphene Kyoungmin Min, Namjung Kim and Ravi Bhadauria May 10, 2010 Abstract In this project, we investigate the shear properties of finite sized
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationFor info and ordering all the 4 versions / languages of this book please visit: http://trl.lab.uic.edu/pon Contents Preface vii Chapter 1 Advances in Atomic and Molecular Nanotechnology Introduction 1
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 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 informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationDiffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study
Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
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 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 informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical
More informationLECTURE 6 : BASICS FOR MOLECULAR SIMULATIONS - Historical perspective - Skimming over Statistical Mechanics - General idea of Molecular Dynamics -
LECTURE 6 : BASICS FOR MOLECULAR SIMULATIONS - Historical perspective - Skimming over Statistical Mechanics - General idea of Molecular Dynamics - Force calculations, structure of MD, equations of motion
More informationAdvanced Molecular Dynamics
Advanced Molecular Dynamics Introduction May 2, 2017 Who am I? I am an associate professor at Theoretical Physics Topics I work on: Algorithms for (parallel) molecular simulations including GPU acceleration
More information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
More informationMolecular Dynamics Simulation of Methanol-Water Mixture
Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,
More informationBrief Review of Statistical Mechanics
Brief Review of Statistical Mechanics Introduction Statistical mechanics: a branch of physics which studies macroscopic systems from a microscopic or molecular point of view (McQuarrie,1976) Also see (Hill,1986;
More informationHeat Transport in Glass-Forming Liquids
Heat Transport in Glass-Forming Liquids by VINAY VAIBHAV The Institute of Mathematical Sciences CIT Campus, Taramani, Chennai 600 113, India. A thesis submitted in partial fulfillment of requirements for
More informationUniversity of Groningen. Characterization of oil/water interfaces van Buuren, Aldert Roelf
University of Groningen Characterization of oil/water interfaces van Buuren, Aldert Roelf IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it.
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
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 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 informationSupporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2016 Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers Dipak
More information1 The Nature of Fluid Flow
1 The Nature of Fluid Flow Aim. To provide the reader with basic knowledge on the treatment of fluid at scales above and below the meso scale. 1.1 INTRODUCTION The fundamentals of fluid flow on a wide
More informationTime-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 informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
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 information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.
1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product
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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationMolecular Aggregation
Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,
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 informationOrganization of NAMD Tutorial Files
Organization of NAMD Tutorial Files .1.1. RMSD for individual residues Objective: Find the average RMSD over time of each residue in the protein using VMD. Display the protein with the residues colored
More informationThe potential energy function is typically called the force field (although it is really an energy function, force is a negative gradient of energy).
7. Molecular Mechanics In molecular mechanics (MM) a molecule is described as a collection of atoms that obey the laws of classical mechanics. The potential energy of a molecule is a simple analytical
More informationCharge equilibration
Charge equilibration Taylor expansion of energy of atom A @E E A (Q) =E A0 + Q A + 1 @Q A 0 2 Q2 A @ 2 E @Q 2 A 0 +... The corresponding energy of cation/anion and neutral atom E A (+1) = E A0 + @E @Q
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationG : 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 informationChemistry 111 Syllabus
Chemistry 111 Syllabus Chapter 1: Chemistry: The Science of Change The Study of Chemistry Chemistry You May Already Know The Scientific Method Classification of Matter Pure Substances States of Matter
More informationThe GROMOS Software for (Bio)Molecular Simulation
The GROMOS Software for (Bio)Molecular Simulation Volume 2: Algorithms and Formulae for Modelling of Molecular Systems September 28, 2017 Contents Chapter 1. Introduction 2-1 Chapter 2. Basic choices
More informationJASS Modeling and visualization of molecular dynamic processes
JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St Petersburg State University, Physics faculty, Department of Computational Physics Supervisor PhD Stepanova Margarita
More informationPHY 481/581. Some classical/quantum physics for the nanometer length scale.
PHY 481/581 Some classical/quantum physics for the nanometer length scale http://creativecommons.org/licenses/by-nc-sa/3.0/ 1 What is nano-science? the science of materials whose properties scale with
More information