Advanced Molecular Molecular Dynamics
|
|
- Sheryl Nelson
- 5 years ago
- Views:
Transcription
1 Advanced Molecular Molecular Dynamics Technical details May 12, 2014
2 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 ) The time scales as p m Some people don t realize this I have seen papers where people scale masses and accelerate MD!
3 Energy comparisons Euler not stable VV & LF very stable Euler not reversible VV+LF symplectic Largest timestep: 2 energy potential kinetic half step kinetic full step time step This is 3.14 steps per period!!! Advantage: great stability! Disadvantage: too (?) great stability
4 The kinetic energy is LF and VV With the same initial conditions LF and VV produce the same trajectory sequence x i v i = (v i-1/2 + v i+1/2 )/2 K = 1/2 m v 2 unless equal: v i2 < [(v i-1/2 + v i+1/2 )/2] 2 The LF velocity are the real steps in the trajectory LF <K> should be close to <U> VV <K> might underestimate <U>
5 Integrator conclusions Not considering v, LF and VV are identical LF gives a good K VV has x and v at the same time but with VV you can also do half steps for v and get the better K Use whatever you want, but be aware of the K issue with VV
6 Today Discuss technical aspects required for MD simulations
7 How to calculate pressure? The pressure is the derivative of the free-energy wrt the volume Pressure definition: Virial definition: P = 2 V (E kin W ) W (r) = 3 2 V du dv potential energy For isotropic scaling: The viral sum: dr i dv = r i 3V W (r) = 1 2 X (virial definitions can contain different factors) r i F i For pair interactions: W (r) = 1 2 i X r ij F ij i<j
8 How to calculate interactions U(r) = X bonds U bond (r)+ X U angle (r)+x angles dihs X X q i q j + + A ij 4 i j>i 0 r ij rij 12 U dih (r) B ij r 6 ij We often only need: F i = du dr i For bonded interactions: simply loop over items in the sum and calculate F (and U) Non-bonded for small molecules: do the double loop Cost O(N 2 ): prohibitive for large systems
9 Non-bonded cut-off 3 oxygen atom pair in water 2 V Coulomb /100 Cut-off interactions beyond a radius V (kj/mol) 1 0 V LJ Fine for LJ Not fine for Coulomb r (nm) Potential should be the integral of the force V co (r) = V (r) V (rc ) r<r c 0 r r c
10 Cut-off effect on Lennard-Jones Lennard-Jones potential decays a r -6 But one atom sees many others For constant density beyond the cut-off, the missing LJ energy is: U LJ = Z 1 4 r 2 1 N C 6 ij r r 6 dr =3 1 N C 6 ij c This missing attraction can be added: long-range or dispersion correction Virial correction is identical, but adds a factor of 6 The pressure correction can be significant r 3 c
11 Inhomogeneous dispersion Uniform correction does not work for inhomogeneous systems e.g. phase boundaries and lipid membranes
12 Cut-off & force fields The basis, including LJ parameters, of most biomolecular force fields is decades old Simulation were done with cut-off s of 0.8/0.9 nm Force-fields were parametrized to give the correct density and H vap with the cut-off used Using a larger cut-off, or dispersion correction will thus result in a too high density To correct this: re-parametrize all LJ interactions Advice: use the right cut-off for the force field!
13 Non-zero force at cut-off With a plain cut-off: F(r c )!=0 This could give integration errors A huge problem for Coulomb No real issue for LJ larger issues for small r Normal Force Shifted Force Shift Function Solution: f(r) 0.5 switch F to 0 shift V to r
14 Coulomb cut-off Coulomb interactions decay as 1/r Cut-off can t be used What to do? >1 r One option: High dielectric: r c weak electrostatics =1 r use a reaction-field
15 Reaction-field Developed for dipoles Linear dipole reaction F V rf = q iq j 4 0 apple 1 + kr r ij 2 c ij k = 1 r 3 c rf 1 2 rf +1 For charges: 3 rf Additional constant c = 1 r c 2 rf required for V(r c )=0 Implicit assumption: uniform background energy (kj/mol) Coulomb RF, c=0 RF with c charge r (nm)
16 Reaction-field cntd Issue: F(r c )!=0 Solution: use rf = 1 F rf (r c )= q iq j r 2 c apple 1 2 rf 2 2 rf +1 conducting or tin-foil boundary condition But shouldn t you match the dielectric of the solvent? Integration errors often worse than deviation in dielectric Also, mismatch goes as: 1/ r
17 Boundary conditions One option: end the system Spherical boundary But what happens at the boundary? apolar liquids OK water problematic What effects on the pressure?
18 Periodic boundary conditions Minimum image convention Much simpler: use r c <L/2
19 Periodic unit or primitive cell triclinic rectangular compact Different shapes can represent the same periodic boundary conditions What matters is not the shape but the periodic shift vectors Different shapes useful for different purposes
20 3D triclinic unit-cells
21 Calculating periodic interactions j j j j i j i i j j Calculating all interactions of one i with many j: you can find which j-image you need easier: move i to different periodic shifts More on this later...
22 Electrostatics in periodic systems X X X NX NX q i q j V = f 2 n x n y n z i j r ij,n We can write down the sum over all charge pairs in all periodic images But as Coulomb goes as 1/r, this sum is only conditionally convergent Direct sums have bad convergence
23 Ewald summation 8 6 1/r direct space reciprocal space 4 V = V dir + V rec + V V dir = f 2 NX i,j X n x X n y X n z q i q j erfc( r ij,n ) r ij,n V rec = f 2 V NX i,j q i q j X m x X m y X m z exp ( m/ ) 2 +2 im (r i r j ) m 2 V 0 = f X N p q 2 i i
24 Particle mesh methods Ewald summation is slow: O(N 2 ) Solution: do the reciprocal part on a mesh Particle-Particle Particle-Mesh (PPPM/P3M) Particle-Mesh Ewald (PME) Most popular SPME (smooth) by Darden et al. spread charges on grid 3D FFT solve in fourier space 3D FFT q gather forces from grid
25 PME PME Parameters cut-off smoothing parameter B spreading order S grid size: M x,m y,m z Accuracy determined by: real space error beyond cut-off: erfc(b r c )/r c spreading accuracy spreading order smoothing parameter grid spacing Computational cost: direct: O(r c3 ) spread: O(#charges*S 3 ) 3D FFT: O(N log(n)) Complex, for SPME no simple analytical formula But, important for performance and accuracy of simulations!
26 PME settings in practice Complex, but also costly! Use what others use with your system and/or software Typical settings: order 4: spread 4 3 =64 points cut-off 0.9 nm grid spacing 0.12 nm #grid point similar to #particles Soon: tools to set parameters based on force accuracy But what does a force accuracy of 0.1 kj/mol/nm 2 mean?
27 Long-range electrostatics methods All methods calculate the same potential and forces computational cost pre-factor communication cost Ewald summation PPPM / Particle Mesh Ewald Fast Multipole Method Multigrid Electrostatics O(N 3/2 ) small (FFT) high O(N log N) small (FFT) high O(N) larger low O(N) larger low
28 Fast multipole method
29 Multigrid electrostatics
30 PME and charged systems With net system charge: implicit assumption: uniform background charge no effect on F, effect on U Problems with nonuniform dielectric Always safer to neutralize with ions!
31 Adding (counter-)ions Simply add ions Na + or Cl - to neutralize Adding just a few ions can lead to sampling problems Better: add Na + and Cl - at physiological concentration, if possible What is physiological concentration locally? Some systems might need different ions DNA/RNA: Na + or Mg 2+?
32 Electrostatics at surfaces One surface is impossible, at least 2 To use PME with PBC: Add 2/3 of vacuum Use dipole correction Yeh&Berkowitz (JCP 111,3155) U z = 2 V M 2 z, F z,i = 4 q i V M z
33 Algorithms to do efficient MD
34 Recap. What do we need to calculate? bonded interactions: cheap non-bonded interactions: expensive maybe PME: expensive integration: cheap
35 Calculating non-bonded interactions Using a Verlet list: Make a Verlet pair list using radius: r list = r c + r buf Calculate interactions for n steps within cut-off r c When to update the list? Option: when a particle moved more than r buf /2 Becomes expensive for large systems
36 Charge groups In the early years of MD cut-offs were used: bad with cut-off electrostatics! partial remedy: use neutral charge groups e.g. group 3 atoms in water: only dipole Remnants of this still in the Gromacs package often used without a Verlet list bad for energy conservation with thermostat fine for most purposes
37 Order of particles & interactions Without particle ordering interactions are randomly distributed in memory: bad performance Sorting the particles on a grid groups interactions good for performance & parallelization
38 Parallel Molecular Dynamics Particle or force decomposition bad memory access characteristics Spatial or domain decomposition good memory access & communication half shell eighth shell midpoint r c r c 1 2 r c (a) (b) (c)
39 Domains & load balancing 7 3 r b r c d r c For inhomogeneous systems load balancing is required Gromacs has full 3D dynamic load balancing
40 Reading Read Frenkel&Smit part V.F on saving CPU time I will put an exercise on the site this afternoon Next lecture: November 10, 10:00 at FB55 No lecture on November 15
Computer 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 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 informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationParallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29
Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29 Outline A few words on MD applications and the GROMACS package The main work in an MD simulation Parallelization Stream computing
More informationOptimizing GROMACS for parallel performance
Optimizing GROMACS for parallel performance Outline 1. Why optimize? Performance status quo 2. GROMACS as a black box. (PME) 3. How does GROMACS spend its time? (MPE) 4. What you can do What I want to
More informationAn FPGA Implementation of Reciprocal Sums for SPME
An FPGA Implementation of Reciprocal Sums for SPME Sam Lee and Paul Chow Edward S. Rogers Sr. Department of Electrical and Computer Engineering University of Toronto Objectives Accelerate part of Molecular
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 informationThe Fast Multipole Method in molecular dynamics
The Fast Multipole Method in molecular dynamics Berk Hess KTH Royal Institute of Technology, Stockholm, Sweden ADAC6 workshop Zurich, 20-06-2018 Slide BioExcel Slide Molecular Dynamics of biomolecules
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 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 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 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 informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationLecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects
Statistical Thermodynamics Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects Dr. Ronald M. Levy ronlevy@temple.edu Bare Molecular Mechanics Atomistic Force Fields: torsion
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 informationN-body simulations. Phys 750 Lecture 10
N-body simulations Phys 750 Lecture 10 N-body simulations Forward integration in time of strongly coupled classical particles or (rigid-body composite objects) dissipation Brownian motion Very general
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 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 informationInteratomic Potentials. The electronic-structure problem
Interatomic Potentials Before we can start a simulation, we need the model! Interaction between atoms and molecules is determined by quantum mechanics: Schrödinger Equation + Born-Oppenheimer approximation
More informationMolecular Dynamics. What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4
Molecular Dynamics What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4 MSE485/PHY466/CSE485 1 The Molecular Dynamics (MD) method
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 informationLong-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC
Long-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC Axel Arnold Institute for Computational Physics Universität Stuttgart October 10, 2012 Long-range interactions gravity Coulomb interaction dipolar
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 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 informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationCoarse-Grained Models!
Coarse-Grained Models! Large and complex molecules (e.g. long polymers) can not be simulated on the all-atom level! Requires coarse-graining of the model! Coarse-grained models are usually also particles
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 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 informationMOLECULAR DYNAMICS SIMULATIONS OF BIOMOLECULES: Long-Range Electrostatic Effects
Annu. Rev. Biophys. Biomol. Struct. 1999. 28:155 79 Copyright c 1999 by Annual Reviews. All rights reserved MOLECULAR DYNAMICS SIMULATIONS OF BIOMOLECULES: Long-Range Electrostatic Effects Celeste Sagui
More informationComputation of non-bonded interactions: Part 1
Computation of non-bonded interactions: Part 1 Computation of non-local (non-bonded) interactions in the potential function E p non-local scales with the number of atoms as N. For the large molecular systems
More informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationMolecular Dynamics. The Molecular Dynamics (MD) method for classical systems (not H or He)
Molecular Dynamics What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4 1 The Molecular Dynamics (MD) method for classical systems
More informationMolecular Dynamics 9/6/16
Molecular Dynamics What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: Lesar Chpt 6, F&S Chpt 4 1 The Molecular Dynamics (MD) method for classical
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 informationThis is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail.
This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Hub, Jochen S.; Groot, Bert L. De; Grubmüller, Helmut;
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 informationMetropolis, 2D Ising model
Metropolis, 2D Ising model You can get visual understanding from the java applets available, like: http://physics.ucsc.edu/~peter/ising/ising.html Average value of spin is magnetization. Abs of this as
More informationEfficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers
Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Jaewoon Jung (RIKEN, RIKEN AICS) Yuji Sugita (RIKEN, RIKEN AICS, RIKEN QBiC, RIKEN ithes) Molecular Dynamics
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 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 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 informationForce fields, thermo- and barostats. Berk Hess
Force fields, thermo- and barostats Berk Hess What is a force field? A force field usually consists of three parts: a set of functional forms parameters for the functional forms that, usually, depend on
More informationTutorial on the smooth particle-mesh Ewald algorithm
Tutorial on the smooth particle-mesh Ewald algorithm These notes explain the smooth particle-mesh Ewald algorithm (SPME) in detail. The accompanying library libpme6 is a fully-featured implementation of
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationLecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together.
Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions Let's get together. Most gases are NOT ideal except at very low pressures: Z=1 for ideal gases Intermolecular interactions come
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 informationTime scales. Notes. Mixed Implicit/Explicit. Newmark Methods [ ] [ ] "!t = O 1 %
Notes Time scales! For using Pixie (the renderer) make sure you type use pixie first! Assignment 1 questions?! [work out]! For position dependence, characteristic time interval is "!t = O 1 % $ ' # K &!
More informationComputational Molecular Biophysics. Computational Biophysics, GRS Jülich SS 2013
Computational Molecular Biophysics Computational Biophysics, GRS Jülich SS 2013 Computational Considerations Number of terms for different energy contributions / E el ~N ~N ~N ~N 2 Computational cost for
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
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 informationSimulations of bulk phases. Periodic boundaries. Cubic boxes
Simulations of bulk phases ChE210D Today's lecture: considerations for setting up and running simulations of bulk, isotropic phases (e.g., liquids and gases) Periodic boundaries Cubic boxes In simulations
More informationBiomolecular modeling II
2015, December 16 System boundary and the solvent Biomolecule in solution typical MD simulations molecular system in aqueous solution task make the system as small as possible (reduce cost) straightforward
More informationThis semester. Books
Models mostly proteins from detailed to more abstract models Some simulation methods This semester Books None necessary for my group and Prof Rarey Molecular Modelling: Principles and Applications Leach,
More informationMultiple time step Monte Carlo simulations: Application to charged systems with Ewald summation
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 1 1 JULY 004 Multiple time step Monte Carlo simulations: Application to charged systems with Ewald summation Katarzyna Bernacki a) Department of Chemistry
More informationScientific Computing II
Technische Universität München ST 008 Institut für Informatik Dr. Miriam Mehl Scientific Computing II Final Exam, July, 008 Iterative Solvers (3 pts + 4 extra pts, 60 min) a) Steepest Descent and Conjugate
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 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 informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 17, 23) I. INTRODUCTION One frequently obtains detailed dynamical information
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 16 Dielectrics and Reaction Field Method David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Molecular models intramolecular
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 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 informationMolecular mechanics. classical description of molecules. Marcus Elstner and Tomáš Kubař. April 29, 2016
classical description of molecules April 29, 2016 Chemical bond Conceptual and chemical basis quantum effect solution of the SR numerically expensive (only small molecules can be treated) approximations
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
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 informationQuantifying Artifacts in Ewald Simulations of Inhomogeneous Systems with a Net Charge
pubs.acs.org/jctc Quantifying Artifacts in Ewald Simulations of Inhomogeneous Systems with a et Charge Jochen S. Hub,*, Bert L. de Groot, Helmut Grubmu ller, and Gerrit Groenhof*,, Institute for Microbiology
More informationFENZI: GPU-enabled Molecular Dynamics Simulations of Large Membrane Regions based on the CHARMM force field and PME
211 IEEE International Parallel & Distributed Processing Symposium : GPU-enabled Molecular Dynamics Simulations of Large Membrane Regions based on the force field and PME Narayan Ganesan, Michela Taufer
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationPractical calculations using first-principles QM Convergence, convergence, convergence
Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................
More information(3.5.1) V E x, E, (3.5.2)
Lecture 3.5 Capacitors Today we shall continue our discussion of electrostatics and, in particular, the concept of electrostatic potential energy and electric potential. The main example which we have
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More informationDraft TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL. Nikolai Shegunov, Ivan Hristov
TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL Nikolai Shegunov, Ivan Hristov March 9, Seminar in IMI,BAS,Soa 1/25 Mathematical model We consider the Hamiltonian :
More informationSeminar in Particles Simulations. Particle Mesh Ewald. Theodora Konstantinidou
Seminar in Particles Simulations Particle Mesh Ewald Theodora Konstantinidou This paper is part of the course work: seminar 2 Contents 1 Introduction 3 1.1 Role of Molecular Dynamics - MD.......................
More informationPH 548 Atomistic Simulation Techniques
PH 548 Atomistic Simulation Techniques Lectures: Lab: 2-0-2-6 Tuesday 12-1 PM Thursday 12-1 PM Monday 2-5 PM P. K. Padmanabhan PH548: Two Excellent Books How to do it? 1 2 F 12 F 23 3 F i = F ij j i F
More informationXXL-BIOMD. Large Scale Biomolecular Dynamics Simulations. onsdag, 2009 maj 13
XXL-BIOMD Large Scale Biomolecular Dynamics Simulations David van der Spoel, PI Aatto Laaksonen Peter Coveney Siewert-Jan Marrink Mikael Peräkylä Uppsala, Sweden Stockholm, Sweden London, UK Groningen,
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationRWTH Aachen University
IPCC @ RWTH Aachen University Optimization of multibody and long-range solvers in LAMMPS Rodrigo Canales William McDoniel Markus Höhnerbach Ahmed E. Ismail Paolo Bientinesi IPCC Showcase November 2016
More informationJunior-level Electrostatics Content Review
Junior-level Electrostatics Content Review Please fill out the following exam to the best of your ability. This will not count towards your final grade in the course. Do your best to get to all the questions
More informationEfficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 5 1 AUGUST 2001 Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems Ruhong Zhou IBM Thomas J. Watson
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 5 The virial theorem 1 The virial theorem We have seen that we can solve the equation of motion for the two-body problem analytically and thus obtain expressions describing the future
More informationNotes on Ewald summation techniques
February 3, 011 Notes on Ewald summation techniques Adapted from similar presentation in PHY 71 he total electrostatic potential energy of interaction between point charges {q i } at the positions {r i
More informationMOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES
MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES V.V. Korenkov 1,3, a, E.G. Nikonov 1, b, M. Popovičová 2, с 1 Joint Institute for Nuclear
More informationForce Field for Water Based on Neural Network
Force Field for Water Based on Neural Network Hao Wang Department of Chemistry, Duke University, Durham, NC 27708, USA Weitao Yang* Department of Chemistry, Duke University, Durham, NC 27708, USA Department
More informationA microsecond a day keeps the doctor away: Efficient GPU Molecular Dynamics with GROMACS
GTC 20130319 A microsecond a day keeps the doctor away: Efficient GPU Molecular Dynamics with GROMACS Erik Lindahl erik.lindahl@scilifelab.se Molecular Dynamics Understand biology We re comfortably on
More informationComputational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, South Carolina 29634, USA.
Supporting Information: Reproducing the Ensemble Average Polar Solvation Energy of a Protein from a Single Structure: Gaussian-Based Smooth Dielectric Function for Macromolecular Modeling Arghya Chakravorty
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 informationCO 2 molecule. Morse Potential One of the potentials used to simulate chemical bond is a Morse potential of the following form: O C O
CO 2 molecule The aim of this project is a numerical analysis of adsorption spectra of CO2 molecule simulated by a double Morse potential function. In the project you should achieve following tasks: 1.
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 informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationChapter 3. Crystal Binding
Chapter 3. Crystal Binding Energy of a crystal and crystal binding Cohesive energy of Molecular crystals Ionic crystals Metallic crystals Elasticity What causes matter to exist in three different forms?
More informationThe Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii
The Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii Fundamental Points About Which to Be Thinking It s clear the PES is useful, so how can I construct it for an arbitrary
More informationAccepted Manuscript. Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations
Accepted Manuscript Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations Robert D. Seel, David J. Hardy, James C. Phillips PII: S001-1(0)001-1
More information5 Classical Molecular Dynamics (MD) Simulations with Tinker
5 Classical Molecular Dynamics (MD) Simulations with Tinker In this set of exercises, you will be briey introduced to the concept of classical force elds and molecular dynamics via propagation of atomic
More informationLecture 3: methods and terminology, part I
Eritis sicut Deus, scientes bonum et malum. Lecture 3: methods and terminology, part I Molecular dynamics, internal coordinates, Z-matrices, force fields, periodic boundary conditions, Hartree-Fock theory,
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 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 informationRama Abbady. Zina Smadi. Diala Abu-Hassan
1 Rama Abbady Zina Smadi Diala Abu-Hassan (00:00) (10:00) Types of Molecules in the Cell 1. Water Molecules: a large portion of the cell mass is water (70% of total cell mass). 2. Organic molecules (carbon
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationClassical Electrostatics for Biomolecular Simulations
pubs.acs.org/cr Classical Electrostatics for Biomolecular Simulations G. Andreś Cisneros, Mikko Karttunen, Pengyu Ren, and Celeste Sagui*, Department of Chemistry, Wayne State University, Detroit, Michigan
More information