Statistical Thermodynamics Lecture 17 Kinetics and Markov State
- Slides: 18
Statistical Thermodynamics Lecture 17: Kinetics and Markov State Models Dr. Ronald M. Levy ronlevy@temple. edu
Computational approaches to kinetics Advanced conformational sampling methods from past several lectures have primarily focused on thermodynamics (ensembles, averages, PMFs) Now we turn our interest to kinetics by differentiating microstates and macrostates There is a vast theoretical literatures on the nonequilibrium statistical mechanical aspects of kinetics which is beyond the scope of this lecture. These two references can provide you with some starting points: R. Zwanzig. Nonequilibrium Statistical Mechanics. 2001. Oxford University Press. Hänggi, Talkner & Borkovec, Rev. Mod. Phys. 62: 251 -341 (1990)
A B D Regions of space free energy C A B reaction coordinate Discrete states
A B Regions of space D C A D B C Discrete states
Kinetics between macrostates as a stochastic process with discrete states stochastic process – a random function of time and past history Markov process – a random function of time and the current (macro)state A D D B C C B A time
Any given realization of a path among the macrostates is unpredictable, but we can still write down equations that describe the time-evolution of probabilities, e. g. P(state=D, time=t | state=A, time=0) In general, a master equation describes the timeevolution of probabilities as follows, Zwanzig, J. Stat. Phys. 30: 255 (1983)
Any given realization of a path among the macrostates is unpredictable, but we can still write down equations that describe the time-evolution of probabilities, e. g. P(state=D, time=t | state=A, time=0) In matrix form,
Two-state Markov State Model A k 1 k 2 B
columns of U are eigenvectors of K eigenvalues of K
• The eigenvalues of K give the characteristic rates of the system • One eigenvalue is always 0. This represents the system in equilibrium, and the eigenvector corresponding to the 0 eigenvalue is proportional to the probabilities of the macrostates at equilibrium. • In general, the decay to equilibrium from any non-equilibrium starting point will consist of a superposition of (N-1) exponentials, where N is the number of macrostates. λi<0 depend only on rates ai depend on rates and starting condition and can be positive or negative
Two-state Markov State Model A k 1 k 2 B if P(A, 0) is 0 etc…
Simulating jump Markov processes How do we construct a “move set” over the kinetic network so that the statistics satisfy ? “Gillespie algorithm”: the amount of time spent in the current state should be an exponential random variable with rate parameter equal to the sum of the rates exiting the current state, and the next state should be chosen with probability proportional to the rate corresponding to that edge
A 1 µs-1 10 µs-1 D B 5 µs-1 C The amount of time t spent in B is a random variable with distribution where kt = 1 + 5 + 10 µs-1, i. e. the mean lifetime in state B is 1/16 µs-1 = 62. 5 ps The probabilities of next jumping to states A, C or D are 1/16, 5/16 and 10/16=5/8, respectively.
A 1 µs-1 1 fs-1 D B 1 ns-1 C The amount of time t spent in B is a random variable with distribution where kt = 1 µs-1 + 1 ns-1 + 1 fs-1 ≈ 1 fs-1, i. e. the mean lifetime in state B is approximately 1 fs. The probabilities of next jumping to states A, C or D are approximately 109, 10 -6, and (1 -10 -9 -10 -6), respectively.
How do we construct MSM from Molecular Dynamics simulation automatically? MSMBuilder 3 from V. S. Pande: https: //github. com/rmcgibbo/msmbuilder 3 Py. Emma from F. Noe: https: //github. com/markovmodel/pyemma Or you can write your own C/C++ codes.
Introduction to Kinetic Lab Goals: • Understand first passage time • Estimate rate constants between two macrostates, folded and unfolded states • Obtain an Arrhenius and anti-Arrhenius plot Zheng, Andrec, Gallicchio, Levy, J. Phys. Chem. B, 2009. DOI: 10. 1021/jp 900445 t.
Introduction to Kinetic Lab PMF along x at three temperatures Rate constants of the 2 -D potential
The end! Thank you!
- Thermodynamics and statistical mechanics
- Thermodynamics and statistical mechanics
- Statistical thermodynamics is a study of
- Statistical thermodynamics
- Steady state approximation examples
- Markov chains tutorial
- Markov chain absorbing state
- 01:640:244 lecture notes - lecture 15: plat, idah, farad
- Heat capacity is path function
- Kelvin planck statement
- Zeroth law of thermodynamics
- In a flow process the work transfer may be of which type
- State function in thermodynamics
- Second law of thermodynamics
- State statistical committee of the republic of azerbaijan
- Kinetics and equilibrium
- Explain law of mass action
- Zero order kinetics
- Planar kinetics of a rigid body force and acceleration