Lecture 14: Advanced Conformational Sampling
|
|
- Hollie Porter
- 6 years ago
- Views:
Transcription
1 Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy
2 Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions between basins are infrequent at room temperature. Barriers are poorly sampled.
3 Biased Sampling Biasing potentials w(x) Thermodynamic properties can be unbiased : A w x =0 = U x dxa x e U x = dxa x e w x e U x e w x = dxe dxe w x e U x e w x Q U w A x e w x U w Q U w e w x U w = A x e w x U w e w x U w
4 Biased Sampling Methods Umbrella sampling Targeted/Steered MD Local elevation, conformational flooding, metadynamics, essential dynamics...
5 Generalized Ensembles Microcanonical Ensemble E =densityof states= d [H E] Canonical ensemble S E =k B ln E p c E Q c = de E e E e E = canonical weight Multicanonical ensemble Q mu = de E W E Set so that: p mu E = W E 1 E =e S E /k B p c E = E e E Q c generalized weight function E W E =constant Q mu Often expressed as: E p mu E M E, W E =e
6 Multicanonical ensemble (cont.) Q mu = de E e M E, = ded [H E]e M E, M [H, ] = d e MD/MC with modified Hamiltonian M(Γ) = M[H(Γ)] The M function (spline) is adjusted by trial and error until the distribution of energies is constant within a given range [Emin, Emax] The multicanonical ensemble samples both low energy and high energy conformations barrier crossing. Unbiasing: A canonical = A e [H M ] mu e [H M ] mu
7 Generalized Ensembles Canonical ensemble Q c = d e H p = e H Q c Extended Ensemble: a parameter λ becomes a dynamical variable Q= d d e f e H = d e f Q c p, = e f H e Q p = d p, = e f d e H = ef Q c Q Q When p(λ) = constant f = lnq c const. Dimensionless free energy at λ
8 Extended Ensemble example: Simulated Tempering λ = temperature Q= dtd e f T e H /k BT = dte f T Q c T Want to sample temperature uniformly within a range: p T = d p,t = ef T Q c T =constant Q Therefore we seek: f T = lnq c T const. In actual applications sample a discrete set of temperatures T m Q= m d e f m T m e mh Generalized Hamiltonian is: H m =H f m T m / m
9 Conformational sampling in Simulated Tempering MC Velocities are not considered (Q Z above). Generalized energy function: Sampling distribution: U m x =U x f m T m / m x,t m e mu m x Two kinds of MC moves: 1. Change of coordinates at constant temperature: x,t m x',t m ' =e m[u x' U x ] 2. Change of temperature at fixed coordinates: x,t m x,t m' ' =e m' m U x f m' f m
10 Conformational sampling in Simulated Tempering MD p,x,t m e m[k p U x ] f m Sampling distribution: 1. Constant temperature MD for n steps at T m with potential function U x think of it as a move: p,x,t m p',x',t m 2. Attempt to temperature move at constant positions + rescaled velocities: p,x,t m T m' /T m 1/2 p,x,t m' m' K'= m' p' 2 2M = m' m m' p 2 2M = m p 2 2M = mk So: ' =e [ m'k' m K m' m U x ] f m' f m =e [ m' m U x ] f m' f m Same as in MC
11 Simulated Tempering (cont.) The weight factors f m (T m ) a.k.a. dimensionless free energies - are adjusted by trial and error until all of the temperatures are visited approximately equally. This can be a time consuming and tedious process. Temperatures can not be spaced too far apart to keep MC acceptance probabilities at a reasonable level.
12 Simulated Tempering (cont.) When the system is visiting high temperatures, barrier crossings are more likely. Then new conformations may cool down and reach the temperature of interest. The samples at the temperature of interest can be used directly to compute thermodynamic averages; each T-ensemble is canonical A T = 1 n T samplesatt A x k (can also unbias from other temperatures WHAM/MBAR, later)
13 Temperature Replica Exchange (a.k.a. Parallel Tempering) In Simulated Tempering equal visitation of temperatures is ensured by the free energy weights f m In Parallel Tempering the same is ensured by having each temperature correspond to an individual replica of the system. T 1 T 2 T 2 T 3 T 4... T n We consider the generalized canonical ensemble of the collection of replicas. Because the replicas are not interacting, the partition function of the RE ensemble is the product of the individual partition functions. Q RE =Q 1 Q 2 Q 3 Q n Q i = d i e ih i Q RE = d 1 d n e i i H i
14 T-RE The state of the RE ensemble is specified by an ordered sequence of momenta/coordinate pairs: Two kinds of moves: 1, 2, n 1. Change of coordinates in one replica (MC or MD): 1,, i,, n 1,, ' i,, n 2. Exchanges of state between a pair of replicas: 1,, i,, j,, n 1,, j,, i,, n 2a. MC: exchange coordinates (no velocities) 2b. MD: rescale velocities at the new temperature
15 Recall Metropolis MC algorithm: probability of accepting move=min ( 1, ρ j ρ i ) =min (1,e βδu ij ) T-RE/MC RE before exp 1 U 1 i U i j U j n U n RE after exp 1 U 1 i U j j U i n U n RE after RE before exp[ j i U i U j ] Accept or reject exchange attempt based on this quantity
16 T-RE/MD R.E. before exp[ i K i U i j K j U j ] R.E. after exp[ i j i K j U j j i j K i U i ] R.E. after R.E. before exp[ j i U i U j ] Same as in MC
17 T-RE R.E. after R.E. before exp[ j i U i U j ] Exchange will be accepted with 100% probability if lower temperature gets the lower energy. Otherwise the exchange has some probability to succeed if either the temperature difference is small or if the energy difference is small, or both. On average larger systems require smaller spacing of temperatures: U i U j U T T j T i =C v T j T i and C v N T 1 T 2 T 1 T 2 p(u) small N p(u) large N U U
18 Replica exchange molecular dynamics Y. Sugita, Y. Okamoto Chem. Phys. Let., 314, 261 (1999)
19 Replica exchange molecular dynamics Y. Sugita, Y. Okamoto Chem. Phys. Let., 314, 261 (1999)
20 Replica exchange molecular dynamics Y. Sugita, Y. Okamoto Chem. Phys. Let., 314, 261 (1999)
21 Replica exchange molecular dynamics Y. Sugita, Y. Okamoto Chem. Phys. Let., 314, 261 (1999)
22 Replica exchange molecular dynamics replica Y. Sugita, Y. Okamoto (1999) Chem. Phys. Let., 314:261
23 Replica exchange molecular dynamics 650 Temperature trajectory of a walker Step number (fsec)
24 Hamiltonian Replica Exchange (HREM) A RE method in which different replicas correspond to (slightly) different potential functions rather then temperature. U(x,λ)=U 0 (x)+v(x,λ) Base energy Perturbation energy HRE exp{ [U x 1, 1 U x n, n ]} The probability of exchange between two replicas in HREM U before =U x 1, 1 U x 2, 2 U after =U x 2, 1 U x 1, 2 where U 1 =U x 1, 2 U x 1, 1 after before =exp[ U 1 U 2 ] U 2 =U x 2, 1 U x 2, 2 (proposed change in energy of conformation x 1 ) (proposed change in energy of conformation x 2 )
25 Two examples of HREM applications REUS: Replica Exchange Umbrella Sampling [Sugita, Kitao, Okamoto, 2000)] ( x, ) = w( x ) V λ ; d λ w x =biasing potential Originally proposed to compute the end-to-end distance PMF of a peptide, in which case the biasing potentials are harmonic restraining potentials of the end-to-end distance d. BEDAM: Binding Energy Distribution Analysis Method [Gallicchio, Lapelosa, Levy 2010] V x, = V RL x V RL x =ligand-receptor interaction energy Replicas are distributed from λ=0 (unbound state) to λ=1 (bound state). Replicas at small λ provide good sampling of ligand conformations whereas replicas at larger λ s provide good statistics for binding free energy estimation. In either case the Weighted Histogram Analysis Method (WHAM) is used to merge the data from multiple replicas [Gallicchio, Andrec, Felts, Levy, 2005]
26 Transition Path Sampling: Simple Picture The trajectory of a Brownian particle moving in a double well First passage (waiting) time and the transition event Large conformational change is a rare and fast event
27 Probability of a Path Classical Mechanics: one dominant path and deterministic Path integral: famous in Quantum Mechanics Stochastic Process: random, statistical result from path ensemble
28 Probability of a Path: Example The probability of a single path: For one dimensional overdamped Langevin dynamics, the probability of a single step is :
29 Transition Path Sampling (Chandler and Bolhuis) Monte Carlo simulation in path space Get one successful path from state A to B From the middle of the original path, shoot to A or B shoot to the other state
30 Weighted Ensemble Method (Zuckerman et al.) Schematic illustration of the Weighted Ensemble (WE) method, using N = 3 bins and M = 2 simulations per bin.
Advanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationMulticanonical parallel tempering
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 13 1 APRIL 2002 Multicanonical parallel tempering Roland Faller, Qiliang Yan, and Juan J. de Pablo Department of Chemical Engineering, University of Wisconsin,
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 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 informationRecovering Kinetics from a Simplified Protein Folding Model Using Replica Exchange Simulations: A Kinetic Network and Effective Stochastic Dynamics
11702 J. Phys. Chem. B 2009, 113, 11702 11709 Recovering Kinetics from a Simplified Protein Folding Model Using Replica Exchange Simulations: A Kinetic Network and Effective Stochastic Dynamics Weihua
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationReplica exchange methodology. John Karanicolas June 2003
Replica exchange methodology John Karanicolas MMTSB@PSC, June 2003 Outline o Motivation o Theory o Practical considerations o MMTSB Tool Set Why? o Sampling on rugged potential energy surfaces is difficult
More informationLecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods
Statistical Thermodynamics Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Dr. Ronald M. Levy ronlevy@temple.edu Free energy calculations Free energy
More informationFree-Energy Calculations in Protein Folding by Generalized-Ensemble Algorithms
Free-Energy Calculations in Protein Folding by Generalized-Ensemble Algorithms Yuji Sugita 1 and Yuko Okamoto 2 Department of Theoretical Studies Institute for Molecular Science Okazaki, Aichi 444-8585,
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
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 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 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 informationLecture V: Multicanonical Simulations.
Lecture V: Multicanonical Simulations. 1. Multicanonical Ensemble 2. How to get the Weights? 3. Example Runs (2d Ising and Potts models) 4. Re-Weighting to the Canonical Ensemble 5. Energy and Specific
More informationMonte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects
Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT
More informationGeneralized Ensembles: Multicanonical Simulations
Generalized Ensembles: Multicanonical Simulations 1. Multicanonical Ensemble 2. How to get the Weights? 3. Example Runs and Re-Weighting to the Canonical Ensemble 4. Energy and Specific Heat Calculation
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 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 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 informationChE 503 A. Z. Panagiotopoulos 1
ChE 503 A. Z. Panagiotopoulos 1 STATISTICAL MECHANICAL ENSEMLES 1 MICROSCOPIC AND MACROSCOPIC ARIALES The central question in Statistical Mechanics can be phrased as follows: If particles (atoms, molecules,
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 informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
More informationComputing 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 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 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 informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
More informationLecture 8: Computer Simulations of Generalized Ensembles
Lecture 8: Computer Simulations of Generalized Ensembles Bernd A. Berg Florida State University November 6, 2008 Bernd A. Berg (FSU) Generalized Ensembles November 6, 2008 1 / 33 Overview 1. Reweighting
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 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 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 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 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 informationChoosing weights for simulated tempering
PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural
More informationMarkov Processes. Stochastic process. Markov process
Markov Processes Stochastic process movement through a series of well-defined states in a way that involves some element of randomness for our purposes, states are microstates in the governing ensemble
More informationClassical Monte-Carlo simulations
Classical Monte-Carlo simulations Graduate Summer Institute Complex Plasmas at the Stevens Insitute of Technology Henning Baumgartner, A. Filinov, H. Kählert, P. Ludwig and M. Bonitz Christian-Albrechts-University
More informationGeneralized-Ensemble Algorithms for Molecular Simulations of Biopolymers
Generalized-Ensemble Algorithms for Molecular Simulations of Biopolymers Ayori Mitsutake, a,b Yuji Sugita, a,b and Yuko Okamoto a,b,1 a Department of Theoretical Studies Institute for Molecular Science
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 new combination of replica exchange Monte Carlo and histogram analysis for protein folding and thermodynamics
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 3 15 JULY 2001 A new combination of replica exchange Monte Carlo and histogram analysis for protein folding and thermodynamics Dominik Gront Department of
More informationModeling & Simulation of Glass Structure
Modeling & Simulation of Glass Structure VCG Lecture 19 John Kieffer Department of Materials Science and Engineering University of Michigan 1 Overview Historical perspective Simulation methodologies Theoretical
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 informationTSTC Lectures: Theoretical & Computational Chemistry
TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez May 2009, Lecture #2 Statistical Mechanics: Structure S = k B log W! Boltzmann's headstone at the Vienna Zentralfriedhof.!! Photo
More informationMultidimensional replica-exchange method for free-energy calculations
JOURNAL OF CHEICAL PHYSICS VOLUE 113, NUBER 15 15 OCTOBER 2000 ultidimensional replica-exchange method for free-energy calculations Yuji Sugita a) Department of Theoretical Studies, Institute for olecular
More informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
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 informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More informationLarge-Scale Asynchronous and Distributed Multidimensional Replica Exchange Molecular Simulations and Efficiency Analysis
WWW.C-CHEM.ORG FULL PAPER Large-Scale Asynchronous and Distributed Multidimensional Replica Exchange Molecular Simulations and Efficiency Analysis Junchao Xia, [a] William F. Flynn, [a,b] Emilio Gallicchio,
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More 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 informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationMarkov Chain Monte Carlo Simulations and Their Statistical Analysis An Overview
Markov Chain Monte Carlo Simulations and Their Statistical Analysis An Overview Bernd Berg FSU, August 30, 2005 Content 1. Statistics as needed 2. Markov Chain Monte Carlo (MC) 3. Statistical Analysis
More informationParallel Tempering Algorithm for. Conformational Studies of Biological Molecules
Parallel Tempering Algorithm for Conformational Studies of Biological Molecules Ulrich H.E. Hansmann 1 Department of Theoretical Studies, Institute for Molecular Science Okazaki, Aichi 444, Japan ABSTRACT
More informationMarkov Chain Monte Carlo Method
Markov Chain Monte Carlo Method Macoto Kikuchi Cybermedia Center, Osaka University 6th July 2017 Thermal Simulations 1 Why temperature 2 Statistical mechanics in a nutshell 3 Temperature in computers 4
More informationMolecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods
Molecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods Charles L. Brooks III MMTSB/CTBP 2006 Summer Workshop CHARMM Simulations The flow of data and information
More informationMonte Carlo simulation of proteins through a random walk in energy space
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 16 22 APRIL 2002 Monte Carlo simulation of proteins through a random walk in energy space Nitin Rathore and Juan J. de Pablo a) Department of Chemical Engineering,
More informationImportance Sampling in Monte Carlo Simulation of Rare Transition Events
Importance Sampling in Monte Carlo Simulation of Rare Transition Events Wei Cai Lecture 1. August 1, 25 1 Motivation: time scale limit and rare events Atomistic simulations such as Molecular Dynamics (MD)
More informationLecture 5: Temperature, Adiabatic Processes
Lecture 5: Temperature, Adiabatic Processes Chapter II. Thermodynamic Quantities A.G. Petukhov, PHYS 743 September 20, 2017 Chapter II. Thermodynamic Quantities Lecture 5: Temperature, Adiabatic Processes
More informationGenerating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques
Generating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques Peter Bolhuis van t Hoff institute for Molecular Sciences University of Amsterdam, The Netherlands
More informationComputational Physics (6810): Session 13
Computational Physics (6810): Session 13 Dick Furnstahl Nuclear Theory Group OSU Physics Department April 14, 2017 6810 Endgame Various recaps and followups Random stuff (like RNGs :) Session 13 stuff
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More informationLocal Hyperdynamics. Arthur F. Voter Theoretical Division Los Alamos National Laboratory Los Alamos, NM, USA
Local Hyperdynamics Arthur F. Voter Theoretical Division National Laboratory, NM, USA Summer School in Monte Carlo Methods for Rare Events Division of Applied Mathematics Brown University Providence, RI
More informationUse of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations
26 J. Chem. Theory Comput. 2007, 3, 26-41 Use of the Weighted Histogram Analysis ethod for the Analysis of Simulated and Parallel Tempering Simulations John D. Chodera,*, William C. Swope, Jed W. Pitera,
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 informationOn the efficient sampling of pathways in the transition path ensemble
On the efficient sampling of pathways in the transition path ensemble Thijs J. H. Vlugt{ and Berend Smit* Department of Chemical Engineering, University of Amsterdam, Nieuwe Achtergracht 166, 1018 WV Amsterdam,
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 informationEquilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach
THE JOURNAL OF CHEMICAL PHYSICS 123, 124912 2005 Equilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach Chakravarthy Ayyagari, a Dmitry Bedrov, and Grant D.
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 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 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 informationThe PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics
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
More informationLecture 20: Binding Free Energy Calculations
Statistical Thermodynamics ecture 20: Binding Free Energy Calculations Dr. onald M. evy ronlevy@temple.edu Outline focus on the binding problem of protein-ligand systems. Statistical theory of non-covalent
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 informationConvergence of replica exchange molecular dynamics
THE JOURNAL OF CHEMICAL PHYSICS 123, 154105 2005 Convergence of replica exchange molecular dynamics Wei Zhang and Chun Wu Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware
More information4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania
1. Introduction Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania Bibliography Computer experiments Ensemble averages and time averages Molecular
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 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 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 informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More informationSampling the conformation space of complex systems, such as
Replica exchange with solute tempering: A method for sampling biological systems in explicit water Pu Liu*, Byungchan Kim*, Richard A. Friesner, and B. J. Berne Department of Chemistry and Center for Biomolecular
More informationNumerical methods in molecular dynamics and multiscale problems
Numerical methods in molecular dynamics and multiscale problems Two examples T. Lelièvre CERMICS - Ecole des Ponts ParisTech & MicMac project-team - INRIA Horizon Maths December 2012 Introduction The aim
More informationMethods for searching and sampling configurational space
International Spring School Statistical Thermodynamics, Santiago de Chile Tuesday, ovember 28, 2017 Lecture 18 Methods for searching and sampling configurational space Prof. Dr. Wilfred F. van Gunsteren
More informationEquilibrium sampling is the engine that drives computational
Replica exchange with nonequilibrium switches Andrew J. Ballard a and Christopher Jarzynski a,b,1 a Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742; and b
More informationComputer Simulation of Peptide Adsorption
Computer Simulation of Peptide Adsorption M P Allen Department of Physics University of Warwick Leipzig, 28 November 2013 1 Peptides Leipzig, 28 November 2013 Outline Lattice Peptide Monte Carlo 1 Lattice
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationBasic Ingredients of Free Energy Calculations: A Review
Feature Article Basic Ingredients of Free Energy Calculations: A Review CLARA D. CHRIST, 1 ALAN E. MARK, 2 WILFRED F. van GUNSTEREN 1 1 Laboratory of Physical Chemistry, Swiss Federal Institute of Technology,
More informationGeneral Construction of Irreversible Kernel in Markov Chain Monte Carlo
General Construction of Irreversible Kernel in Markov Chain Monte Carlo Metropolis heat bath Suwa Todo Department of Applied Physics, The University of Tokyo Department of Physics, Boston University (from
More informationRare Events. Transition state theory Bennett-Chandler Approach 16.2 Transition path sampling16.4
Rare Events Transition state theory 16.1-16.2 Bennett-Chandler Approach 16.2 Transition path sampling16.4 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory
More informationWang-Landau Monte Carlo simulation. Aleš Vítek IT4I, VP3
Wang-Landau Monte Carlo simulation Aleš Vítek IT4I, VP3 PART 1 Classical Monte Carlo Methods Outline 1. Statistical thermodynamics, ensembles 2. Numerical evaluation of integrals, crude Monte Carlo (MC)
More informationExchange frequency in replica exchange molecular dynamics
THE JOURNAL OF CHEMICAL PHYSICS 128, 024103 2008 Exchange frequency in replica exchange molecular dynamics Daniel Sindhikara Department of Physics, University of Florida, Gainesville, Florida 32611-8435,
More informationStatistical methods in atomistic computer simulations
Statistical methods in atomistic computer simulations Prof. Michele Ceriotti, michele.ceriotti@epfl.ch This course gives an overview of simulation techniques that are useful for the computational modeling
More informationSIMULATED TEMPERING: A NEW MONTECARLO SCHEME
arxiv:hep-lat/9205018v1 22 May 1992 SIMULATED TEMPERING: A NEW MONTECARLO SCHEME Enzo MARINARI (a),(b) and Giorgio PARISI (c) Dipartimento di Fisica, Università di Roma Tor Vergata, Via della Ricerca Scientifica,
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 0 Simple Biasing Methods David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo simulation Markov chain to generate
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 informationEnhanced sampling via molecular dynamics II: Unbiased approaches
Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY
More informationProtein folding is a fundamental problem in modern structural
Protein folding pathways from replica exchange simulations and a kinetic network model Michael Andrec, Anthony K. Felts, Emilio Gallicchio, and Ronald M. Levy* Department of Chemistry and Chemical Biology
More informationTRANSITION PATH SAMPLING: Throwing Ropes Over Rough Mountain Passes, in the Dark
Annu. Rev. Phys. Chem. 2002. 53:291 318 DOI: 10.1146/annurev.physchem.53.082301.113146 Copyright c 2002 by Annual Reviews. All rights reserved TRANSITION PATH SAMPLING: Throwing Ropes Over Rough Mountain
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 9 I. OVERVIEW Our treatment of the classical ensembles makes clear that the free energy is a quantity of particular importance in statistical mechanics.
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches
More information