The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics
|
|
- Fay Charles
- 6 years ago
- Views:
Transcription
1 The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics Davide Branduardi, Theoretical Molecular Biophysics Group Max Planck for Biophysics, Frankfurt am Main, Germany
2 TyOutline What is free energy and why is important Traditional QC approach Explicit sampling Enhanced sampling PLUMED plugin for free energy calculations 2
3 TyOutline What is free energy and why is important Traditional QC approach Explicit sampling Enhanced sampling PLUMED plugin for free energy calculations 3
4 Ty What is free energy and why is important F (s) = k B T ln P (s) P (s) = e V (x) k bt δ(s (x) s)dx Q conformational transitions and equilibria s= state (bound/unbound, reactant/products, phasea/phaseb) ligand binding mechanism of transporters or channels chemical reactions and catalysis phase transitions... P (out) P (near) P (in) 4
5 TyOutline What is free energy and why is important Traditional QC approach Explicit sampling Enhanced sampling PLUMED plugin for free energy calculations 5
6 TyTraditional QC approach Rigid rotor/ harmonic approximation is well known [1] F (s) =H(s) TS(s) S(s) =R + R ln (q t (s)q e (s)q r (s)q v (s)) + RT ln qt (s)q e (s)q r (s)q v (s) T trasl, elect, rot, vib partition functions where the partition functions are obtained from the calculated vibrational spectrum (Hessian), moments of inertia, electronic structure H(s) =H el (s)+h t (s)+h v (s)+h r (s) It requires electronic structure and Hessian in the local minima and transition states: 1 single conformation per state. [1] McQuarrie, Simon Physical Chemistry, a molecular approach 6
7 Ty Finding the RC (I) Static approach: look for saddle points by using chemical intuition and TS optimization. m1 m2 7910
8 Ty Finding the RC Static approach: look for saddle points by using chemical intuition and TS optimization. m1 Verify then by IRC: move forward and backward to verify to end in minima
9 Ty Finding the RC Static approach: look for saddle points by using chemical intuition and TS optimization. m1 m2 Verify then by IRC: move forward and backward to verify to end in minima
10 Ty Finding the RC Static approach: look for saddle points by using chemical intuition and TS optimization. m1 m2 Verify then by IRC: move forward and backward to verify to end in minima. G = H T S Rates may be derived from Eyring equation [2] : k(t )= k bt hc 0 e G/RT [2] Eyring, Chem. Rev vol. 17 pp
11 TyTraditional QC approach PRO: one calculation per point (3 single points calculations per rate and free energy differences) CONS: optimizing TS is not trivial (redundant coordinates), the landscape must be very simple (Free) energy ortho coor Reaction coor 11
12 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness 12
13 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness A TS B Perfectly parallel basins, Not too bad (small correction over static approaches). Typical case of symmetry in a reaction 12
14 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness 12
15 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness Perfectly parallel basins but energetically asymmetric, Can be done with static QC but need to sample all the paths 12
16 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness 12
17 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness Rough paths: vibrations are meaningless, paths from IRC are not representative of reactions: USE MD BASED METHODS 12
18 TyProblems with traditional QC approach Many problems are not smooth as before: many parallel valleys (conformers), solvent induced roughness 12
19 TyOutline What is free energy and why is important Traditional QM approach Explicit sampling Enhanced sampling PLUMED plugin for free energy calculations 13
20 TyThe ingredients of the probability picture Free energy is connected to the probability of a given state to occur in a single measurement F (s) = k B T ln P (s) P (s) = Probability e V (x) k B T δ(s (x) s)dx e V (x) k B T dx State (bond length, folded/ unfolded descr., liquid/ crystal/amorphous) 14
21 TyExplicit sampling via molecular dynamics a real system (cuvette) a a F (s) = k B T ln P (s) b a a a a b a a a a a set of independent objects (ensemble) a a a a b b MD: evolve in time with Newton s eq. of motion use ergodic hypothesis a b a a a a b a a t 2 t 3 t... all the available states can be sampled with t->infinite P (s) 15
22 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy React V 1 Prod 2 k BT RC ν = ν 0 exp ( β V ) ν 0 = s 1 16
23 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy React V 1 Prod 2 k BT RC MD (DFT) can reach hundreds of ps. Lucky if you see one single event! ν = ν 0 exp ( β V ) ν 0 = s 1 16
24 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy React V 1 Prod 2 k BT RC MD (DFT) can reach hundreds of ps. Lucky if you see one single event! ν = ν 0 exp ( β V ) property ν 0 = s 1 time 16
25 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy 100 ps React V 1 Prod 2 k BT RC MD (DFT) can reach hundreds of ps. Lucky if you see one single event! ν = ν 0 exp ( β V ) property ν 0 = s 1 time 16
26 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy 100 ps React V 1 Prod 2 k BT RC MD (DFT) can reach hundreds of ps. Lucky if you see one single event! ν = ν 0 exp ( β V ) property ν 0 = s 1 time If you want thermodynamics you must average over events 16
27 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy React 100 ps V 1 Prod 2 k BT RC MD (DFT) can reach hundreds of ps. Lucky if you see one single event! ν = ν 0 exp ( β V ) property ν 0 = s 1 time If you want thermodynamics you must average over events property time 16
28 TyLimitations Transition of 8 kcal/mol can take up to milliseconds at 300K Free Energy 100 ps V 1 MD (DFT) can reach hundreds of ps. Lucky if you see one single event! ν = ν 0 exp ( β V ) ν 0 = s 1 This 2 k BT React is not Prod your solution RC property time If you want thermodynamics you must average over events property time 16
29 TyOutline What is free energy and why is important Traditional QM approach Explicit sampling Enhanced sampling PLUMED plugin for free energy calculations 17
30 TyThe ancestor: Torrie & Valleau [3] Free energy: Free energy for a biased system: use biasing potential to overcome barriers and retrieve free energies [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
31 TyLessons from Torrie & Valleau Biasing potentials: Free Energy React Prod 1 2 k BT RC [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
32 TyLessons from Torrie & Valleau Biasing potentials: Free Energy React biasing pot Prod 1 2 k BT RC [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
33 TyLessons from Torrie & Valleau Free Energy Biasing potentials: React biasing pot 1 2 k BT resulting potential Prod 1 2 k BT RC Free Energy Metastabilities are React removed Prod RC [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
34 TyLessons from Torrie & Valleau Free Energy Biasing potentials: React biasing pot 1 2 k BT resulting potential Prod 1 2 k BT RC Free Energy Metastabilities are React removed Prod RC allow to measure free energies: useful allow to overcome barriers: cheap need only an additional term to standard DFT forces: easy [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
35 TyAdaptive sampling methods Adaptive Biasing Force Darve & Pohorille, J. Chem. Phys. 2001, vol. 115, pp Adaptive Umbrella Sampling Mezei, M. J Comput Phys 1987, vol. 68, pp Self Healing Umbrella Sampling Marsili et. al J. Chem. Phys. B vol. 110, pp MetadynamicsLaio & Parrinello, PNAS 2002, vol. 20, pp Conformational Flooding Grubmüller, Phys Rev E 1995, vol. 52, pp
36 Ty Adaptive sampling methods (e.g. metadynamics) without ES with ES Φ Ψ from milli secs... free energy is faster not black box to nano secs! 21
37 Ty Thermodynamic-Integration-like approaches The ancestor is mean force calculation through harmonic potential F (s) = k B T ln P (s) F (s) s s0 = k B T 1 e V (x) k B T δ(s (x) s)dx s k B T δ(s (x) s)dx e V (x) s0 Dirac delta function! Gaussian function (I like this trick!) F (s) s s0 k B T e V (x) k B T e 1 k 2k B T (s (x) s) 2 dx It is an average of a biased simulation! F (s) s s0 1 e V (x)+k/2(s (x) s 0 ) 2 k B T dx s e V (x) k B T e k 2k B T (s (x) s) 2 dx k(s (x) s 0 )e V (x)+k/2(s (x) s0 ) 2 k B T dx s0 = ks (x) s 0 bias,s0 (opposite of the mean force ) 22
38 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
39 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
40 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
41 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
42 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
43 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
44 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
45 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
46 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
47 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
48 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
49 TyThermodynamic Integration " 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability +WHAM etc... 23
50 TyThermodynamic-Integration-like approaches Thermodynamic integration WHAM Roux, Comp Phys Comm 1995, vol. 91, pp Free energy perturbation Beveridge, DiCapua, Annu Rev Biophys Biophys Chem 1989, vol 18, pp Beveridge, DiCapua, Annu Rev Biophys Biophys Chem 1989, vol 18, pp Jarzynski-equation based approaches (steered-md) Jarzynski Phys Rev Lett 1997, vol. 78, pp Crooks-equation based approaches (two directions steered-md) Crooks J Stat Phys 1998, vol. 90, pp SN2 reaction: if you like this movie come to the tutorial to see how to make it! 24
51 TyAdaptive sampling vs TI-like methods Adaptive sampling: The problem is multidimensional but probably the accessible phase space is limited: they explore only what you need The free energy landscape has competitive, parallel reactive paths (solvent degrees of freedom, rotations) that can be overcome at some point TI-like: In 1-d (max 2-d, via WHAM or, better in DFT, rbf fitting schemes) very effective Sometimes trivially parallel (many umbrellas at same time) 25
52 TyThe problem of the s F (s) = k B T ln P (s) A multidimensional landscape: turn into monodimensional 26
53 TyThe problem of the s F (s) = k B T ln P (s) A multidimensional landscape: turn into monodimensional CV2 b a CV1 26
54 TyThe problem of the s F (s) = k B T ln P (s) A multidimensional landscape: turn into monodimensional CV2 b a CV1 26
55 TyThe problem of the s F (s) = k B T ln P (s) A multidimensional landscape: turn into monodimensional CV2 b a CV1 Lodola, Branduardi, De Vivo, Capoferri, Mor, Piomelli, Cavalli, PlosOne, 2012, vol 7, pp. e
56 TyThe problem of the s Use machine learning approaches Tribello et. al. PNAS 2010, vol. 107, pp
57 TyOutline What is free energy and why is important Traditional QM approach Explicit sampling Enhanced sampling PLUMED plugin for free energy calculations 28
58 TyThe PLUMED plugin: Bonomi, Branduardi et al., Comp. Phys. Comm vol. 180 pp.1961 small connection with MD code slow fast Already in many classical MD codes (LAMMPS, DLPOLY, SANDER, GROMACS, NAMD) now in FHI-aims! 29
59 TyPLUMED: what it contains Methods Constrained relaxation Umbrella sampling Metadynamics (welltempered) d-afed Steered MD Harmonic walls Ratcheting Reweighting schemes String method Descriptors (with der.) Distances, angles, dihedrals Coordination numbers Contact maps Path collective variables Energy Function of variables SPRINT variables... + MetaGUI analysis interface and PLUMED preparation GUI in VMD+analysis tools! 30
60 TyIn FHI-aims compile separately: Makefile.meta in control.in plumed.dat this enables PLUMED set a distance as set anotherdistance as set repulsive boundaries set metadynamics 31
61 TyAcknowledgements Luca Ghiringhelli (FHI, Berlin) The PLUMED developers (ETHZ, USILU, Cambridge, UCSF, SISSA, Mount Sinai, Curtin, EPFL, MPI-BP) Giovanni Bussi (SISSA, Triest) Michele Parrinello (ETHZ/USILU, Lugano) José Faraldo-Gómez (MPI-BP, Frankfurt/Main) 32
62 Ty Thank you! 33
Computing free energies with PLUMED 2.0
Computing free energies with PLUMED 2.0 Davide Branduardi Formerly at MPI Biophysics, Frankfurt a.m. TyOutline of the talk Relevance of free energy computation Free-energy from histograms: issues Tackling
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationa micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics Density functional theory and beyond: Computational materials science for real materials
More informationPerforming Metadynamics Simulations Using NAMD
Performing Metadynamics Simulations Using NAMD Author: Zhaleh Ghaemi Contents 1 Introduction 2 1.1 Goals of this tutorial.......................... 2 2 Standard metadynamics simulations 3 2.1 Free energy
More informationa micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics FHI-aims Developers and Users Meeting Organized by Fritz Haber Institute, CECAM-mm1p,
More informationFree energy calculations
Free energy calculations Jochen Hub & David van der Spoel Overview Free energies and Probabilities Thermodynamic cycles (Free energy perturbation (FEP)) Thermodynamic integration (TI) (Jarzynski equality
More informationBiomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)
Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationa micro PLUMED2.0 Tutorial for ESPResSo users PLUMED2.0: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED2.0 Tutorial for ESPResSo users PLUMED2.0: portable plugin for free-energy calculations with molecular dynamics ESPResSo summer school 2013 Organized by CECAM, SimTech, ICP and University
More informationSupporting Information
1 Supporting Information 2 3 4 Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds 5 6 7 Fabio Pietrucci 1, José C. Aponte 2,3,*, Richard Starr 2,4, Andrea Pérez-Villa
More informationSampling the free energy surfaces of collective variables
Sampling the free energy surfaces of collective variables Jérôme Hénin Enhanced Sampling and Free-Energy Calculations Urbana, 12 September 2018 Please interrupt! struct bioinform phys chem theoretical
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationStatistical Mechanics for Proteins
The Partition Function From Q all relevant thermodynamic properties can be obtained by differentiation of the free energy F: = kt q p E q pd d h T V Q ), ( exp 1! 1 ),, ( 3 3 3 ),, ( ln ),, ( T V Q kt
More informationMETAGUI A VMD EXTENSION TO ANALYZE AND VISUALIZE METADYNAMICS SIMULATIONS
METAGUI A VMD EXTENSION TO ANALYZE AND VISUALIZE METADYNAMICS SIMULATIONS Alessandro Laio SISSA & DEMOCRITOS, Trieste Coworkers: Xevi Biarnes Fabio Pietrucci Fabrizio Marinelli Metadynamics (Laio A. and
More informationenhanced sampling methods recovering <a> Accelerated MD Introduction to Accelerated Molecular Dynamics an enhanced sampling method
enhanced sampling methods Introduction to Accelerated Molecular Dynamics an enhanced sampling method Yi Wang McCammon group Aug 2011, NBCR Summer Institute Umbrella sampling B. ROUX, "The calculation of
More informationCrossing the barriers - simulations of activated processes
Crossing the barriers - simulations of activated processes Mgr. Ján Hreha for 6 th Student Colloquium and School on Mathematical Physics Faculty of Mathematics, Physics and Informatics Comenius University
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationSub -T g Relaxation in Thin Glass
Sub -T g Relaxation in Thin Glass Prabhat Gupta The Ohio State University ( Go Bucks! ) Kyoto (January 7, 2008) 2008/01/07 PK Gupta(Kyoto) 1 Outline 1. Phenomenology (Review). A. Liquid to Glass Transition
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationControlling fluctuations
Controlling fluctuations Michele Parrinello Department of Chemistry and Applied Biosciences ETH Zurich and ICS, Università della Svizzera Italiana, Lugano, Switzerland Today s menu Introduction Fluctuations
More informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
More informationSupporting Information
Supporting Information The effect of strong acid functional groups on electrode rise potential in capacitive mixing by double layer expansion Marta C. Hatzell, a Muralikrishna Raju, b Valerie J. Watson,
More informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More informationModelling of Reaction Mechanisms KJE-3102
Modelling of Reaction Mechanisms KJE-3102 Kathrin H. Hopmann Kathrin.hopmann@uit.no Outline Potential energy surfaces Transition state optimization Reaction coordinates Imaginary frequencies Verification
More informationFree Energy Simulation Methods
Free Energy Simulation Methods Free energy simulation methods Many methods have been developed to compute (relative) free energies on the basis of statistical mechanics Free energy perturbation Thermodynamic
More informationMaking the best of a bad situation: a multiscale approach to free energy calculation
Making the best of a bad situation: arxiv:1901.04455v2 [physics.comp-ph] 15 Jan 2019 a multiscale approach to free energy calculation Michele Invernizzi,, and Michele Parrinello, Department of Physics,
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 informationBayesian reconstruction of free energy profiles from umbrella samples
Bayesian reconstruction of free energy profiles from umbrella samples Thomas Stecher with Gábor Csányi and Noam Bernstein 21st November 212 Motivation Motivation and Goals Many conventional methods implicitly
More informationFree energy calculations using molecular dynamics simulations. Anna Johansson
Free energy calculations using molecular dynamics simulations Anna Johansson 2007-03-13 Outline Introduction to concepts Why is free energy important? Calculating free energy using MD Thermodynamical Integration
More informationEnergy Profiles and Chemical Reactions
Energy rofiles and Chemical eactions + B C D E egensburg er-la Norrby Modeling Kinetics Molecular Modeling Stationary points Energies Barriers eaction rates DFT, ab initio, MM Ground & Transition states
More informationLandscape Approach to Glass Transition and Relaxation. Lecture # 4, April 1 (last of four lectures)
Landscape Approach to Glass Transition and Relaxation Lecture # 4, April 1 (last of four lectures) Relaxation in the glassy state Instructor: Prabhat Gupta The Ohio State University (gupta.3@osu.edu) Review
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 informationBiology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key
Biology 5357 Chemistry & Physics of Biomolecules Examination #1 Proteins Module September 29, 2017 Answer Key Question 1 (A) (5 points) Structure (b) is more common, as it contains the shorter connection
More informationQuick reference guide on PLUMED with Quantum ESPRESSO
Quick reference guide on PLUMED with Quantum ESPRESSO Changru Ma SISSA, Trieste March 30, 2011 Contents 1 Introduction 2 1.1 Overview................................... 2 1.2 Collective variables.............................
More informationAspects of nonautonomous molecular dynamics
Aspects of nonautonomous molecular dynamics IMA, University of Minnesota, Minneapolis January 28, 2007 Michel Cuendet Swiss Institute of Bioinformatics, Lausanne, Switzerland Introduction to the Jarzynski
More informationTransition Theory Abbreviated Derivation [ A - B - C] # E o. Reaction Coordinate. [ ] # æ Æ
Transition Theory Abbreviated Derivation A + BC æ Æ AB + C [ A - B - C] # E A BC D E o AB, C Reaction Coordinate A + BC æ æ Æ æ A - B - C [ ] # æ Æ æ A - B + C The rate of reaction is the frequency of
More informationA self-learning algorithm for biased molecular dynamics
A self-learning algorithm for biased molecular dynamics Tribello, G. A., Ceriotti, M., & Parrinello, M. (2010). A self-learning algorithm for biased molecular dynamics. Proceedings of the National Academy
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 informationStatistical thermodynamics for MD and MC simulations
Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems
More informationCharacterization of the free-energy landscapes of proteins by NMR-guided metadynamics
Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics B Results and Discussion BIOPHYSICS ND COMPUTTIONL BIOLOGY SI Text Inset B χ χ Inset SI Text Inset C Folding of GB3
More informationOverview & Applications. T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 04 June, 2015
Overview & Applications T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 4 June, 215 Simulations still take time Bakan et al. Bioinformatics 211. Coarse-grained Elastic
More informationMetadynamics with adaptive Gaussians
Metadynamics with adaptive Gaussians Davide Branduardi Theoretical Molecular Biophysics Group, Max Planck Institute for Biophysics, Max-von-Laue strasse 5, 60438, Frankfurt am Main, Germany Giovanni Bussi
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 informationA self-learning algorithm for biased molecular dynamics
A self-learning algorithm for biased molecular dynamics Gareth A. Tribello, Michele Ceriotti and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETHZ Zurich USI-Campus,
More informationSome Accelerated/Enhanced Sampling Techniques
Some Accelerated/Enhanced Sampling Techniques The Sampling Problem Accelerated sampling can be achieved by addressing one or more of the MD ingredients Back to free energy differences... Examples of accelerated
More informationAb initio statistical mechanics
Ab initio statistical mechanics Luca M. Ghiringhelli Fritz-Haber-Institut der MPG, Berlin Hands-on Workshop Density-Functional Theory and Beyond: Accuracy, Efficiency and Reproducibility in Computational
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationIntroduction to Theories of Chemical Reactions. Graduate Course Seminar Beate Flemmig FHI
Introduction to Theories of Chemical Reactions Graduate Course Seminar Beate Flemmig FHI I. Overview What kind of reactions? gas phase / surface unimolecular / bimolecular thermal / photochemical What
More informationEnhanced sampling of transition states
Enhanced sampling of transition states arxiv:1812.09032v1 [physics.chem-ph] 21 Dec 2018 Jayashrita Debnath,, Michele Invernizzi,, and Michele Parrinello,,, Department of Chemistry and Applied Biosciences,
More informationIntro/Review of Quantum
Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the
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 informationIntro/Review of Quantum
Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the
More informationMulti-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018
Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods
More informationA general overview of. Free energy methods. Comer et al. (2014) JPCB 119:1129. James C. (JC) Gumbart Georgia Institute of Technology, Atlanta
A general overview of Free energy methods Comer et al. (2014) JPCB 119:1129. James C. (JC) Gumbart Georgia Institute of Technology, Atlanta Computational Biophysics Workshop Georgia Tech Nov. 16 2016 Outline
More informationFree energy calculations with alchemlyb
Free energy calculations with alchemlyb Oliver Beckstein Arizona State University SPIDAL Teleconference 2019-02-01 Binding free energy measure of how strong a protein P and a ligand X stick together key
More informationFlexMD Manual. ADF Program System Release 2013
FlexMD Manual ADF Program System Release 2013 Scientific Computing & Modelling NV Vrije Universiteit, Theoretical Chemistry De Boelelaan 1083; 1081 HV Amsterdam; The Netherlands WWW: www.scm.com E-mail:
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationEnergy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding
Energy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding John K. Prentice Boulder, CO BioMed Seminar University of New Mexico Physics and Astronomy Department
More informationComputing free energy: Replica exchange
Computing free energy: Replica exchange Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls Macroscopic i a t e regime
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
More informationElectronic Supporting Information Topological design of porous organic molecules
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2017 Electronic Supporting Information Topological design of porous organic molecules Valentina Santolini,
More informationTesting Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers
computation Article Testing Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers S. Alexis Paz 1, ID and Cameron F. Abrams 3, * ID 1 Departamento de Química Teórica
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationPrinciples of Molecular Spectroscopy
Principles of Molecular Spectroscopy What variables do we need to characterize a molecule? Nuclear and electronic configurations: What is the structure of the molecule? What are the bond lengths? How strong
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 information5.62 Physical Chemistry II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 5.62 Physical Chemistry II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.62 Spring 2007 Lecture
More informationProteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability
Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Dr. Andrew Lee UNC School of Pharmacy (Div. Chemical Biology and Medicinal Chemistry) UNC Med
More informationFree energy calculations
Free energy calculations Berk Hess May 5, 2017 Why do free energy calculations? The free energy G gives the population of states: ( ) P 1 G = exp, G = G 2 G 1 P 2 k B T Since we mostly simulate in the
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationGherman Group Meeting. Thermodynamics and Kinetics and Applications. June 25, 2009
Gherman Group Meeting Thermodynamics and Kinetics and Applications June 25, 2009 Outline Calculating H f, S, G f Components which contribute to H f, S, G f Calculating ΔH, ΔS, ΔG Calculating rate constants
More informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More informationThermochemistry in Gaussian
Thermochemistry in Gaussian Joseph W. Ochterski, Ph.D. help@gaussian.com c 2000, Gaussian, Inc. June 2, 2000 Abstract The purpose of this paper is to explain how various thermochemical values are computed
More informationThermodynamics and Kinetics
Thermodynamics and Kinetics C. Paolucci University of Notre Dame Department of Chemical & Biomolecular Engineering What is the energy we calculated? You used GAMESS to calculate the internal (ground state)
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 informationEssential dynamics sampling of proteins. Tuorial 6 Neva Bešker
Essential dynamics sampling of proteins Tuorial 6 Neva Bešker Relevant time scale Why we need enhanced sampling? Interconvertion between basins is infrequent at the roomtemperature: kinetics and thermodynamics
More informationLimitations of temperature replica exchange (T-REMD) for protein folding simulations
Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research pitera@us.ibm.com Anomalies in protein folding kinetic thermodynamic 322K
More informationBiomolecular modeling III
2016, January 5 Déjà vu Enhanced sampling Biomolecular simulation Each atom x, y, z coordinates Déjà vu Enhanced sampling Expression for energy the force field = 1 2 + N i i k i (r i r 0 i ) 2 + 1 2 N
More informationFree energy, electrostatics, and the hydrophobic effect
Protein Physics 2016 Lecture 3, January 26 Free energy, electrostatics, and the hydrophobic effect Magnus Andersson magnus.andersson@scilifelab.se Theoretical & Computational Biophysics Recap Protein structure
More informationarxiv: v2 [physics.comp-ph] 4 Feb 2016
Sampling Free Energy Surfaces as Slices by Combining Umbrella Sampling and Metadynamics Shalini Awasthi, Venkat Kapil and Nisanth N. Nair 1 arxiv:1507.06764v2 [physics.comp-ph] 4 Feb 2016 Department of
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 informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
More informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationMolecular dynamics (MD) and the Monte Carlo simulation
Escaping free-energy minima Alessandro Laio and Michele Parrinello* Centro Svizzero di Calcolo Scientifico, Via Cantonale, CH-6928 Manno, Switzerland; and Department of Chemistry, Eidgenössische Technische
More informationDetermining Energy Barriers and Selectivities of a Multi-Pathway System With Infrequent Metadynamics
Determining Energy Barriers and Selectivities of a Multi-Pathway System With Infrequent Metadynamics Christopher D. Fu 1, Luiz F.L. Oliveira 1 and Jim Pfaendtner 1,a) 1 Department of Chemical Engineering,
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 informationhydrate systems Gránásy Research Institute for Solid State Physics & Optics H-1525 Budapest, POB 49, Hungary László
Phase field theory of crystal nucleation: Application to the hard-sphere and CO 2 Bjørn Kvamme Department of Physics, University of Bergen Allégaten 55, N-5007 N Bergen, Norway László Gránásy Research
More informationFigure 1: Transition State, Saddle Point, Reaction Pathway
Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces
More informationProbing the Origins of Intermolecular Vibrational and Relaxational Dynamics in Organic Solids with CP2K
Probing the Origins of Intermolecular Vibrational and Relaxational Dynamics in Organic Solids with CP2K Michael Ruggiero Department of Chemical Engineering and Biotechnology, University of Cambridge CP2K
More informationThe Potential Energy Surface (PES) Preamble to the Basic Force Field Chem 4021/8021 Video II.i
The Potential Energy Surface (PES) Preamble to the Basic Force Field Chem 4021/8021 Video II.i The Potential Energy Surface Captures the idea that each structure that is, geometry has associated with it
More informationRandom Walks A&T and F&S 3.1.2
Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will
More informationExtending Parallel Scalability of LAMMPS and Multiscale Reactive Molecular Simulations
September 18, 2012 Extending Parallel Scalability of LAMMPS and Multiscale Reactive Molecular Simulations Yuxing Peng, Chris Knight, Philip Blood, Lonnie Crosby, and Gregory A. Voth Outline Proton Solvation
More informationVan Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006) Martini Workshop 2015 Coarse Graining Basics Alex de Vries Every word or concept, clear as it may seem to be, has only a limited range of applicability
More informationMulti-scale approaches in description and design of enzymes
Multi-scale approaches in description and design of enzymes Anastassia Alexandrova and Manuel Sparta UCLA & CNSI Catalysis: it is all about the barrier The inside-out protocol: Big Aim: development of
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationA Parallel-pulling Protocol for Free Energy Evaluation INTRODUCTION B A) , p
A Parallel-pulling Protocol for Free Energy Evaluation Van Ngo Collaboratory for Advance Computing and Simulations (CACS, Department of Physics and Astronomy, University of Southern California, 3651 Watt
More informationLiquid to Glass Transition
Landscape Approach to Glass Transition and Relaxation (Lecture # 3, March 30) Liquid to Glass Transition Instructor: Prabhat Gupta The Ohio State University (gupta.3@osu.edu) PKGupta(OSU) Landscape #3
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationSupporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface
Supporting Information for Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite Edge Water Interface Zhen Jiang, Konstantin Klyukin, and Vitaly Alexandrov,, Department
More informationMCB100A/Chem130 MidTerm Exam 2 April 4, 2013
MCB1A/Chem13 MidTerm Exam 2 April 4, 213 Name Student ID True/False (2 points each). 1. The Boltzmann constant, k b T sets the energy scale for observing energy microstates 2. Atoms with favorable electronic
More informationRepresenting High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks
Energy Landscapes, Chemnitz 200 Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks Jörg Behler Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum D-44780 Bochum,
More information