Extending the Tools of Chemical Reaction Engineering to the Molecular Scale
|
|
- Virgil Burns
- 5 years ago
- Views:
Transcription
1 Extending the Tools of Chemical Reaction Engineering to the Molecular Scale Multiple-time-scale order reduction for stochastic kinetics James B. Rawlings Department of Chemical and Biological Engineering March 31, 29 Model Reduction in Reacting Flows University of Notre Dame Rawlings Molecular reaction engineering 1 / 35
2 Outline 1 Introduction to stochastic kinetics 2 Model reduction fast reactions and reactive intermediates 3 Catalyst example with fast diffusion 4 Virus example with fast fluctuation 5 Further reading Rawlings Molecular reaction engineering 2 / 35
3 Introduction to stochastic kinetics Stochastic kinetics Small species populations Species numbers are integers, reactions cause integer jumps Large fluctuations in species numbers and reaction rates Biological networks and catalyst particles Model reduction Develop reduced models from stochastic chemical reactions. These models must meet the following requirements: Simpler than the full model (fewer reactions, fewer parameters, or faster simulation times) Converge to the full model as a specified parameter goes to zero Rawlings Molecular reaction engineering 3 / 35
4 Stochastic simulation method kinetic Monte Carlo n A r 2 r 1 B A k 1 k2 k 1 = 2 k 2 = 1 n A = 6 n B = time (sec) Random number r 1 r 1 +r 2 = KMC Algorithm 1 Choose which reaction 2 Choose time step 3 Repeat 1 r 2 r 1 +r 2 = Which reaction: Time step: Sample from an exponential distribution where the distribution mean is the sum of reaction rates. Rawlings Molecular reaction engineering 4 / 35
5 KMC simulations and probability ni Multiple KMC simulations, A k1 B n B k2 Count Histogram at t=2sec, for 5 simulations n A time (sec) n A KMC simulations are samples of a probability distribution that evolves in time. We can write the evolution equation for the probability density (master equation). Rawlings Molecular reaction engineering 5 / 35
6 Chemical master equation dp(x) dt dp dt = = AP N rxn r j (x ν j )P(x ν j ) r j (x)p(x) }{{}}{{} j=1 rate into state x rate out of state x Master equation example P(n A).15 B A k 1 k2 n A = 1, n B = k 1 = 2, k 2 = 1 11 possible states 11 Coupled ODEs Time (sec) Rawlings Molecular reaction engineering 6 / 35 n A
7 Master equation Important points Chemical master equation dp(x) dt dp dt = = AP N rxn r j (x ν j )P(x ν j ) r j (x)p(x) }{{}}{{} j=1 rate into state x rate out of state x Often the dimensionality of the master equation makes direct solution infeasible The master equation shows what probability distribution is sampled in a KMC simulation A reduced master equation can lead to a new/faster simulation schemes Rawlings Molecular reaction engineering 7 / 35
8 Kinetics of multiple time scales A k 1 k 1 B k 2 C 1 Deterministic - One Time Scale k1 = 2 k 1 =.5 k2 =.5 1 Deterministic - Two Time Scales k1 = 1 k 1 = 1 k2 =.5 1 Deterministic - Two Time Scales k1 = 2 k 1 = 2 k2 = na nc 6 na nb nc 6 na nc nb time (sec) KMC - One Time Scale k1 = 2 k 1 =.5 k2 = time (sec) KMC - Two Time Scales k1 = 1 k 1 = 1 k2 =.5 1 nb time (sec) KMC - Two Time Scales k1 = 2 k 1 = 2 k2 = na nb nc 8 6 na nc 8 6 na nc time (sec) One time scale 2 nb time (sec) Reaction equilibrium 2 nb time (sec) Reactive intermediate Rawlings Molecular reaction engineering 8 / 35
9 Deterministic model reductions x non-qssa species, y QSSA species dx dt = f (x, y) ɛdy dt = g(x, y) Classical QSSA dx dt = f (x, y) = g(x, y) DAE reduced model Singular Perturbation QSSA x = X + ɛx 1 + ɛ 2 X 2 + O(ɛ 3 ) y = Y + ɛy 1 + ɛ 2 Y 2 + O(ɛ 3 ) Collect like powers of ɛ Equations for dx dt is the reduced model Separate models for fast and slow time scale Rawlings Molecular reaction engineering 9 / 35
10 SPA on the master equation Our objective Apply singular perturbation analysis to develop a reduced master equation. dp(a, b, c) dt A k 1 k 1 B k 2 C = k 1 (a + 1)P(a + 1, b 1, c) + k 1 (b + 1)P(a 1, b + 1, c) + k 2 (b + 1)P(a, b + 1, c 1) (k 1 a + k 1 b + k 2 b)p(a, b, c) P(a, b, c) = W (a, b, c) + ɛw 1 (a, b, c) + ɛ terms: W (a, b, c) = if b > In this limit b is always zero Rawlings Molecular reaction engineering 1 / 35
11 SPA on the master equation ɛ 1 terms: Reduced master equation dw (a,, c) dt Reduced mechanism = k(a + 1)W (a + 1,, c 1) kaw (a,, c) A C r = k 1k 2 k 1 +k 2 a Stochastic same as deterministic SPA mechanism Same mechanisms due to linearity First-order correction, b b = f (W (a,, c)) + O(ɛ 2 ) b = k 1 k 1 +k 2 a Rawlings Molecular reaction engineering 11 / 35
12 Comparison of mechanisms A 2B r 1 = k 1 a r 1 = k 1 b(b 1) B C r 2 = k 2 b 2 Stoch SPA A 2C ( r = k 1 k 2 k 1 2 +k c ) a na = 25, nb =, nc = k1 = 1 k 1 = 1 k2 = 1 Det SPA A 2C r = k 1 a Det QSSA A 2C k 2 + k k 1k 1 a r = k 2 4k 1 Full Model Stoch SPA Det QSSA Det SPA Time (sec) Rawlings Molecular reaction engineering 12 / 35
13 Catalyst Example Number of molecules (n) A k 1 B k 2 C B + D k 3 B + E k 2, k 3 k 1 a 2 b Time (sec) e Stoch SPA mechanism A D + A C E + C 2D + A 2E + C nd + A ne + C k 1 a r = 1 + K 3 d K 3 = k 3 k 2 K 3 d r 1 = r 1 + K 3 (d 1) r 2 = r 1 K 3 (d 1) 1 + K 3 (d 2) r n = r n 1 K 3 (d + 1 n) 1 + K 3 (d n) Deterministic SPA mechanism A C D + A E + A r = k 1a r 1 = k3k1 ad k 2 Rawlings Molecular reaction engineering 13 / 35
14 Conclusions Stochastic quasi-steady-state approximation QSSA species are removed from stochastic models with SPA Stochastic QSSA mechanisms different than deterministic QSSA mechanisms Application of stochastic QSSA: Reduces the number of kinetic parameters Speeds up KMC simulations (fewer events) Rawlings Molecular reaction engineering 14 / 35
15 Conclusions Stochastic quasi-steady-state approximation QSSA species (# of molecules) >> <1 sqspa sqspa dqc dqspa >>1 3 non QSSA species (# of molecules) Rawlings Molecular reaction engineering 15 / 35
16 Catalytic surface reaction modeling Assumptions for this talk Two dimensional surface with a lattice for adsorption, diffusion, reaction, and desorption. Square lattice, Z=4 All sites have identical properties Constant temperature Adsorbed CO molecules exhibit nearest neighbor repulsions CO O 2 CO 2 CO-black, O-gray, Empty-white. Rawlings Molecular reaction engineering 16 / 35
17 Model mechanism and time scales Adsorption CO(g)+ i α CO i O 2(g)+ i + j β O i + O j Desorption CO i γ CO(g) + i O i + O j ρ O 2(g) + i k Reaction CO i + O r j CO2(g)+ i d Diffusion CO i + 1 j i +CO j d O i + 2 j i + O j + j + j 1/sec α = 1.6 β =.8 γ =.8 ρ =.1 k r = 1 d d Rawlings Molecular reaction engineering 17 / 35
18 Singular perturbation on the master equation Surface reaction master equation dp(n, x) dt x - microscopic configuration n - number of each species X rxn = k j a j (n ν j, x ν x,j )P(n ν j, x ν x,j ) k j a j (n, x)p(n, x) j=1 X diff + d j a j (n, x ν x,j )P(n, x ν x,j ) d j a j (n, x)p(n, x) j=1 Singular perturbation P(n, x) = W (n, x) + ɛw 1 (n, x) + ɛ 2 W 2 (n, x) + ɛ = 1/d ɛ terms: Diffusion equilibration equations for W (x n) Rawlings Molecular reaction engineering 18 / 35
19 Slow time-scale evolution equation ɛ 1 terms : Reduced master equation dw (n) dt N rxn = k i s i (n ν i ) W (n ν i ) k i s i (n) W (n) i=1 What have we gained? Removed micro-states from the master equation Lattice Size Species Micro-states Coverage states N s = N s = N s = Tractable number of states, master equation can be solved Rawlings Molecular reaction engineering 19 / 35
20 Slow time-scale evolution equation dw (n) dt Reaction propensities N rxn = k i s i (n ν i ) W (n ν i ) k i s i (n) W (n) i=1 s i (x) number of reaction i on configuration x: n CO =45 black, n O =8 gray s CO O =26 s CO O =22 s CO O =26 s CO O =23 s CO O = 24.9 s i (n) = x s i(x)w (x n) Calculate with diffusion only KMC Rawlings Molecular reaction engineering 2 / 35
21 Reduced master equation solution (5x5 lattice) Steady-state probability distribution (5x5).25.4 Probability of n CO, 5x5 lattice Probability n CO Time (sec) n CO Rawlings Molecular reaction engineering 21 / 35
22 Verification of perturbation method P(n, x) = W (n) + ɛw 1 (n) + O(ɛ 2 ) x ɛ = 1/d As the diffusion rate increases P(n) approaches W (n) Steady-state probability distribution (5 5) Probability Reduced ME d=1 d=1 d=1 d=1 d=.1 d=.1 d= n CO Rawlings Molecular reaction engineering 22 / 35
23 Conclusions Surface reactions in the infinite diffusion limit SPA can be used to eliminate spatial configuration states in a reduced master equation. The reduced master equation has sufficiently few states to be simulated on small lattices. Reduced master equations of surface reactions can be used to motivate reduced KMC and reduced ODE models. Rawlings Molecular reaction engineering 23 / 35
24 Model for Vesicular Stomatitis Virus (VSV) infection!!! " # $%& # "$%& "$%& "$%& I is encapsidation of viral genome II is replication of encapsidated genome III is transcription of genome to messenger RNA " # $%&! '()* '()* '()* '()* '()* '()* Rawlings Molecular reaction engineering 24 / 35
25 Onset of fast fluctuations in the N protein Molecule level I III N mrna ( )RNA 1258 II N protein (+)RNA genomes (+)RNA 1258 ( )RNA genomes Time (hours) Features of simulation Presence of fast fluctuating and rapidly rising species Fast fluctuations slow the full KMC simulation Motivates the formulation of a simpler example to understand this phenomenon Rawlings Molecular reaction engineering 25 / 35
26 Fast fluctuation and rapid rise in VSV biology k ( )RNA + L 1 1 ( )RNA + L2 k ( )RNA + L 2 2 2( )RNA + L1 2( )RNA k 3 ( )RNA The viral genome is amplified by first two reactions The free viral proteins and messages are not amplified Values of parameters k 1, k 2 and k 3 may cause fast fluctuation in polymerases along with rapid amplification of viral genome Rawlings Molecular reaction engineering 26 / 35
27 Fast fluctuation and rapid rise Idealized problem A + G k 1 C + G r1 = 1 Ω k 1ag C + G 2G k 2 2G + A r2 = 1 Ω k 2cg k 3 G r3 = 1 Ω k g(g 1) 3 2 Species Initial number Rate constant ( m3 /mol s) A 3 k 1 = C k 2 = G 1 k 3 = Rawlings Molecular reaction engineering 27 / 35
28 The full SSA on the system No. of C molecules t (sec) No. of G molecules t (sec) time-step (sec) t (sec) Rawlings Molecular reaction engineering 28 / 35
29 The hybrid SSA - Ω technique At large population of G we want to switch to a continuous description for it: g = Ωφ G + Ω 1/2 ξ φ G is the deterministic evolution term and ξ is the continuous noise in the evolution of G We can obtain approximation for the evolution of system using hybrid SSA - Ω technique Approximation of pdf of C W (c) = (1 + q) N ( N c ) q (N c) Deterministic evolution of G dφ G dt = γ 1 c φ G k 3 2 φ2 G N q = k 2 k 1 Initial number of polymerases Ratio of rate constants Rawlings Molecular reaction engineering 29 / 35
30 Comparison of full SSA with hybrid SSA - Ω Full SSA Hybrid SSA - Ω 3 3 No. of C molecules 2 1 No. of C molecules 2 1 No. of G molecules t (sec) t (sec) No. of G molecules t (sec) t (sec) Rawlings Molecular reaction engineering 3 / 35
31 Comparison of full SSA with hybrid SSA - Ω Probability densities of C from SSA and from hybrid SSA - Ω Probability Hybrid SSA Ω SSA No. of C particles Rawlings Molecular reaction engineering 31 / 35
32 Conclusions QSSA and fast fluctuations Hybrid SSA Ω expansion matches closely the full SSA Computation speed increases by factor of 45 Application to kinetic virus infection models Rawlings Molecular reaction engineering 32 / 35
33 Acknowledgments Dr. Ethan A. Mastny, BP Alaska Dr. Eric L. Haseltine, Vertex Pharmaceuticals Rishi Srivastava, UW NSF #CNS Rawlings Molecular reaction engineering 33 / 35
34 Further reading Stochastic reaction equilibrium Haseltine, E. L. and J. B. Rawlings. Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics. J. Chem. Phys., 117(15): , October 22. Haseltine, E. L. and J. B. Rawlings. On the origins of approximations for stochastic chemical kinetics. J. Chem. Phys., 123:164115, October 25. Cao, Y., D. T. Gillespie, and L. R. Petzold. The slow-scale stochastic simulation algorithm. J. Chem. Phys., 122(1):14116, January 25. Samant, A. and D. G. Vlachos. Overcoming stiffness in stochastic simulation stemming from partial equilibrium: A multiscale Monte Carlo algorithm. J. Chem. Phys., 123:144114, 25. Salis, H. and Y. Kaznessis. Accurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions. J. Chem. Phys., 122(5):5413, February 25. Mastny, E. A., E. L. Haseltine, and J. B. Rawlings. Stochastic simulation of catalytic surface reactions in the fast diffusion limit. J. Chem. Phys., 125(19):194715, November 26. Rawlings Molecular reaction engineering 34 / 35
35 Further reading Stochastic QSSA Rao, C. V. and A. P. Arkin. Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm. J. Chem. Phys., 118(11): , March 23. van Kampen, N. G. Stochastic Processes in Physics and Chemistry. Elsevier Science Publishers, Amsterdam, The Netherlands, second edition, Mastny, E. A., E. L. Haseltine, and J. B. Rawlings. Two classes of quasi-steady-state model reductions for stochastic kinetics. J. Chem. Phys., 127(9):9416, September 27. Hensel, S., J. B. Rawlings, and J. Yin. Stochastic kinetic modeling of vesicular stomatitis virus intracellular growth. Submitted for publication in Bulletin of Math. Bio., May 28. Rawlings Molecular reaction engineering 35 / 35
Supporting Information for The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise
TWCCC Texas Wisconsin California Control Consortium Technical report number 2011-03 Supporting Information for The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise Rishi Srivastava
More informationFast Probability Generating Function Method for Stochastic Chemical Reaction Networks
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 71 (2014) 57-69 ISSN 0340-6253 Fast Probability Generating Function Method for Stochastic Chemical Reaction
More informationDerivations for order reduction of the chemical master equation
2 TWMCC Texas-Wisconsin Modeling and Control Consortium 1 Technical report number 2006-02 Derivations for order reduction of the chemical master equation Ethan A. Mastny, Eric L. Haseltine, and James B.
More informationEfficient Leaping Methods for Stochastic Chemical Systems
Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some
More informationParameter Estimation in Stochastic Chemical Kinetic Models. Rishi Srivastava. A dissertation submitted in partial fulfillment of
Parameter Estimation in Stochastic Chemical Kinetic Models by Rishi Srivastava A dissertation submitted in partial fulfillment of the requirement for the degree of Doctor of Philosophy (Chemical Engineering)
More information2008 Hotelling Lectures
First Prev Next Go To Go Back Full Screen Close Quit 1 28 Hotelling Lectures 1. Stochastic models for chemical reactions 2. Identifying separated time scales in stochastic models of reaction networks 3.
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationStochastic Modeling of Chemical Reactions
Stochastic Modeling of Chemical Reactions Eric Haseltine Department of Chemical Engineering University of Wisconsin-Madison TWMCC 27 September Texas-Wisconsin Modeling and Control Consortium 22TWMCC Outline
More informationSTOCHASTIC CHEMICAL KINETICS
STOCHASTIC CHEICAL KINETICS Dan Gillespie GillespieDT@mailaps.org Current Support: Caltech (NIGS) Caltech (NIH) University of California at Santa Barbara (NIH) Past Support: Caltech (DARPA/AFOSR, Beckman/BNC))
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationStochastic models of biochemical systems
Stochastic models of biochemical systems David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison University of Amsterdam November 14th, 212 Stochastic models
More informationKinetic Monte Carlo Simulation of Molecular Processes on Supported Metal Particles
INFAP CONICET UNSL Kinetic Monte Carlo Simulation of Molecular Processes on Supported Metal Particles Giorgio Zgrablich Universidad Nacional de San Luis San Luis, Argentina CONTENTS Motivation Non-equilibrium
More informationComputational methods for continuous time Markov chains with applications to biological processes
Computational methods for continuous time Markov chains with applications to biological processes David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Penn.
More informationNested Stochastic Simulation Algorithm for Chemical Kinetic Systems with Disparate Rates. Abstract
Nested Stochastic Simulation Algorithm for Chemical Kinetic Systems with Disparate Rates Weinan E Department of Mathematics and PACM, Princeton University, Princeton, NJ 08544, USA Di Liu and Eric Vanden-Eijnden
More informationNested Stochastic Simulation Algorithm for KMC with Multiple Time-Scales
Nested Stochastic Simulation Algorithm for KMC with Multiple Time-Scales Eric Vanden-Eijnden Courant Institute Join work with Weinan E and Di Liu Chemical Systems: Evolution of an isothermal, spatially
More informationStochastic Simulation of Biochemical Reactions
1 / 75 Stochastic Simulation of Biochemical Reactions Jorge Júlvez University of Zaragoza 2 / 75 Outline 1 Biochemical Kinetics 2 Reaction Rate Equation 3 Chemical Master Equation 4 Stochastic Simulation
More informationThe Adaptive Explicit-Implicit Tau-Leaping Method with Automatic Tau Selection
The Adaptive Explicit-Implicit Tau-Leaping Method with Automatic Tau Selection Yang Cao Department of Computer Science, 660 McBryde Hall, Virginia Tech, Blacksburg, VA 24061 Daniel T. Gillespie Dan T.
More informationStochastic Chemical Kinetics
Stochastic Chemical Kinetics Joseph K Scott November 10, 2011 1 Introduction to Stochastic Chemical Kinetics Consider the reaction I + I D The conventional kinetic model for the concentration of I in a
More informationThe Morse Potential. ChE Reactive Process Engineering. Let s have a look at a very simple model reaction :
The Morse Potential L7-1 Let s have a look at a very simple model reaction : We can describe the F + H-H energy of this system based on the relative position of the atoms. This energy is (in the most
More informationFEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH. Yiming Lou and Panagiotis D. Christofides
FEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH Yiming Lou and Panagiotis D. Christofides Department of Chemical Engineering University of California, Los Angeles IEEE 2003 Conference
More informationModelling in Systems Biology
Modelling in Systems Biology Maria Grazia Vigliotti thanks to my students Anton Stefanek, Ahmed Guecioueur Imperial College Formal representation of chemical reactions precise qualitative and quantitative
More informationSTOCHASTIC MODELING OF BIOCHEMICAL REACTIONS
STOCHASTIC MODELING OF BIOCHEMICAL REACTIONS Abhyudai Singh and João Pedro Hespanha* Department of Electrical and Computer Engineering University of California, Santa Barbara, CA 93101. Abstract The most
More informationWhen do diffusion-limited trajectories become memoryless?
When do diffusion-limited trajectories become memoryless? Maciej Dobrzyński CWI (Center for Mathematics and Computer Science) Kruislaan 413, 1098 SJ Amsterdam, The Netherlands Abstract Stochastic description
More informationCoarse bifurcation analysis of kinetic Monte Carlo simulations: A lattice-gas model with lateral interactions
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Coarse bifurcation analysis of kinetic Monte Carlo simulations: A lattice-gas model with lateral interactions Alexei G. Makeev a) Department
More informationStochastic models for chemical reactions
First Prev Next Go To Go Back Full Screen Close Quit 1 Stochastic models for chemical reactions Formulating Markov models Two stochastic equations Simulation schemes Reaction Networks Scaling limit Central
More informationStochastic modelling in Mathematical Biology
Stochastic modelling in Mathematical Biology Daniel Sánchez-Taltavull Centre de Recerca Matemàtica March 4th 2013 D. Sánchez-Taltavull (CRM) Stochastic modelling in Mathematical Biology March 4th 2013
More informationLecture 7: Simple genetic circuits I
Lecture 7: Simple genetic circuits I Paul C Bressloff (Fall 2018) 7.1 Transcription and translation In Fig. 20 we show the two main stages in the expression of a single gene according to the central dogma.
More informationTOPIC 6: Chemical kinetics
TOPIC 6: Chemical kinetics Reaction rates Reaction rate laws Integrated reaction rate laws Reaction mechanism Kinetic theories Arrhenius law Catalysis Enzimatic catalysis Fuente: Cedre http://loincognito.-iles.wordpress.com/202/04/titanic-
More informationStochastic Hybrid Systems: Modeling, analysis, and applications to networks and biology
research supported by NSF Stochastic Hybrid Systems: Modeling, analysis, and applications to networks and biology João P. Hespanha Center for Control Engineering and Computation University of California
More informationNumerical Solution of the Fokker- Planck Approximation of the Chemical Master Equation
IT Licentiate theses 2005-010 Numerical Solution of the Fokker- Planck Approximation of the Chemical Master Equation PAUL SJÖBERG UPPSALA UNIVERSITY Department of Information Technology Numerical Solution
More informationLognormal Moment Closures for Biochemical Reactions
Lognormal Moment Closures for Biochemical Reactions Abhyudai Singh and João Pedro Hespanha Abstract In the stochastic formulation of chemical reactions, the dynamics of the the first M -order moments of
More informationChemical Kinetics. Topic 7
Chemical Kinetics Topic 7 Corrosion of Titanic wrec Casón shipwrec 2Fe(s) + 3/2O 2 (g) + H 2 O --> Fe 2 O 3.H 2 O(s) 2Na(s) + 2H 2 O --> 2NaOH(aq) + H 2 (g) Two examples of the time needed for a chemical
More informationFinite-size poisoning in heterogeneous catalysis
J. Phys. A: Math. Gen. 23 (1990) L613-J-619. Printed in the UK LE ITER TO THE EDITOR Finite-size poisoning in heterogeneous catalysis D ben-avrahamt9, S Rednert, D B Considinet and P MeakinS t Center for
More informationKinetic Monte Carlo (KMC) Kinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationarxiv: v3 [math.na] 28 Jun 2013
A Convergent Reaction-Diffusion Master Equation Samuel A. Isaacson a) Department of Mathematics and Statistics Boston University Boston, MA 02215, USA arxiv:1211.6772v3 [math.na] 28 Jun 2013 The reaction-diffusion
More informationEquilibrium Time, Permutation, Multiscale and Modified Multiscale Entropies for Low-High Infection Level Intracellular Viral Reaction Kinetics
Equilibrium Time, Permutation, Multiscale and Modified Multiscale Entropies for Low-High Infection Level Intracellular Viral Reaction Kinetics Fariborz Taherkhani 1, Farid Taherkhani 2* 1 Department of
More informationFeedback Control of Surface Roughness Using Stochastic PDEs
R&D NOTE Feedback Control of Surface Roughness Using Stochastic PDEs Yiming Lou and Panagiotis D. Christofides Dept. of Chemical Engineering, University of California, Los Angeles, CA 90095 DOI 0.002/aic.0299
More informationAn Application of Perturbation Methods in Evolutionary Ecology
Dynamics at the Horsetooth Volume 2A, 2010. Focused Issue: Asymptotics and Perturbations An Application of Perturbation Methods in Evolutionary Ecology Department of Mathematics Colorado State University
More informationKINETIC MONTE CARLO MOLECULAR SIMULATIONS FOR FUEL CELL APPLICATIONS AND SURFACE REACTIONS
KINETIC MONTE CARLO MOLECULAR SIMULATIONS FOR FUEL CELL APPLICATIONS AND SURFACE REACTIONS 1 A Thesis Defense In Partial Fulfillment of the Requirements for the Degree of Master of Science in Mechanical
More informationLongtime behavior of stochastically modeled biochemical reaction networks
Longtime behavior of stochastically modeled biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison ASU Math Biology Seminar February 17th, 217 Overview
More informationA Multiple Time-Step Finite State Projection Algorithm for the Solution to the Chemical Master Equation
A Multiple Time-Step Finite State Projection Algorithm for the Solution to the Chemical Master Equation Brian Munsky and Mustafa Khammash Department of Mechanical and Environmental Engineering University
More informationA Moment Closure Method for Stochastic Chemical Reaction Networks with General Kinetics
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 70 2013 785-800 ISSN 0340-6253 A Moment Closure Method for Stochastic Chemical Reaction Networks with General
More informationIntroduction to stochastic multiscale modelling in tumour growth
Introduction to stochastic multiscale modelling in tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 3 CIMPA, Santiago de Cuba, June 2016
More informationTheoretical Models for Chemical Kinetics
Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech DRAFT v0.4a, January 16, 2011 c California Institute of Technology All rights reserved. This manuscript is for review purposes
More informationAdvanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS
Experimental Kinetics and Gas Phase Reactions Advanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS Professor Angelo R. Rossi http://homepages.uconn.edu/rossi Department of Chemistry, Room
More informationnatural development from this collection of knowledge: it is more reliable to predict the property
1 Chapter 1 Introduction As the basis of all life phenomena, the interaction of biomolecules has been under the scrutiny of scientists and cataloged meticulously [2]. The recent advent of systems biology
More informationSlow Scale Tau-leaping Method
Slow Scale Tau-leaping Method Yang Cao Linda Petzold Abstract For chemical systems involving both fast and slow scales, stiffness presents challenges for efficient stochastic simulation. Two different
More informationAn Introduction to Stochastic Simulation
Stephen Gilmore Laboratory for Foundations of Computer Science School of Informatics University of Edinburgh PASTA workshop, London, 29th June 2006 Background The modelling of chemical reactions using
More informationKinetic Monte Carlo. Heiko Rieger. Theoretical Physics Saarland University Saarbrücken, Germany
Kinetic Monte Carlo Heiko Rieger Theoretical Physics Saarland University Saarbrücken, Germany DPG school on Efficient Algorithms in Computational Physics, 10.-14.9.2012, Bad Honnef Intro Kinetic Monte
More informationMore Protein Synthesis and a Model for Protein Transcription Error Rates
More Protein Synthesis and a Model for Protein James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 3, 2013 Outline 1 Signal Patterns Example
More informationarxiv: v1 [cond-mat.stat-mech] 10 Jul 2018
Molecular search with conformational change: One-dimensional discrete-state stochastic model arxiv:1807.03740v1 cond-mat.stat-mech] 10 Jul 2018 Jaeoh Shin 1 1, 2, 3 and Anatoly B. Kolomeisky 1 Department
More informationNIH Public Access Author Manuscript Comput Biol Chem. Author manuscript; available in PMC 2010 June 1.
NIH Public Access Author Manuscript Published in final edited form as: Comput Biol Chem. 2009 June ; 33(3): 205 215. doi:10.1016/j.compbiolchem.2009.03.002. Hybrid stochastic simulations of intracellular
More informationAP CHEMISTRY CHAPTER 12 KINETICS
AP CHEMISTRY CHAPTER 12 KINETICS Thermodynamics tells us if a reaction can occur. Kinetics tells us how quickly the reaction occurs. Some reactions that are thermodynamically feasible are kinetically so
More informationLecture 4 The stochastic ingredient
Lecture 4 The stochastic ingredient Luca Bortolussi 1 Alberto Policriti 2 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste Via Valerio 12/a, 34100 Trieste. luca@dmi.units.it
More informationDynamical Monte-Carlo Simulation of Surface Kinetics
Dynamical Monte-Carlo Simulation of Surface Kinetics V. Guerra and J. Loureiro Centro de Física dos Plasmas, Instituto Superior Técnico, 1049-001 Lisboa, Portugal Abstract. This wor presents a theoretical
More informationStochastic approaches for modelling in vivo reactions
Stochastic approaches for modelling in vivo reactions T. E. Turner a, S. Schnell b,c, and K. Burrage d. a Oxford Centre for Industrial and Applied Mathematics, Mathematical Institute, 24-29 St. Giles,
More informationChemical Kinetics. The Concepts
Chemical Kinetics Copyright c 2015 by Nob Hill Publishing, LLC The purpose of this chapter is to provide a framework for determining the reaction rate given a detailed statement of the reaction chemistry.
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 informationCybergenetics: Control theory for living cells
Department of Biosystems Science and Engineering, ETH-Zürich Cybergenetics: Control theory for living cells Corentin Briat Joint work with Ankit Gupta and Mustafa Khammash Introduction Overview Cybergenetics:
More informationStochastic models, patterns formation and diffusion
Stochastic models, patterns formation and diffusion Duccio Fanelli Francesca Di Patti, Tommaso Biancalani Dipartimento di Energetica, Università degli Studi di Firenze CSDC Centro Interdipartimentale per
More informationAccurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions
THE JOURNAL OF CHEMICAL PHYSICS 122, 054103 2005 Accurate hybrid schastic simulation of a system of coupled chemical or biochemical reactions Howard Salis and Yiannis Kaznessis a) Department of Chemical
More information2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements
1. Why stochastic?. Mathematical descriptions (i) the master equation (ii) Langevin theory 3. Single cell measurements 4. Consequences Any chemical reaction is stochastic. k P d φ dp dt = k d P deterministic
More informationIntroduction to stochastic multiscale modelling of tumour growth
Introduction to stochastic multiscale modelling of tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 1 CIMPA, Santiago de Cuba, June 2016
More informationName AP CHEM / / Chapter 12 Outline Chemical Kinetics
Name AP CHEM / / Chapter 12 Outline Chemical Kinetics The area of chemistry that deals with the rate at which reactions occur is called chemical kinetics. One of the goals of chemical kinetics is to understand
More informationarxiv: v2 [q-bio.qm] 12 Jan 2017
Approximation and inference methods for stochastic biochemical kinetics - a tutorial review arxiv:1608.06582v2 [q-bio.qm] 12 Jan 2017 David Schnoerr 1,2,3, Guido Sanguinetti 2,3, and Ramon Grima 1,3,*
More informationModeling and Systems Analysis of Gene Regulatory Networks
Modeling and Systems Analysis of Gene Regulatory Networks Mustafa Khammash Center for Control Dynamical-Systems and Computations University of California, Santa Barbara Outline Deterministic A case study:
More informationarxiv:chem-ph/ v2 11 May 1995
A Monte Carlo study of temperature-programmed desorption spectra with attractive lateral interactions. A.P.J. Jansen Laboratory of Inorganic Chemistry and Catalysis arxiv:chem-ph/9502009v2 11 May 1995
More informationPart 3: Introduction to Master Equation and Complex Initial Conditions in Lattice Microbes
Part 3: Introduction to Master Equation Cells: and Complex Initial Conditions in Lattice re cells Microbes en Biophysics, and UC urgh, June 6-8, 2016 rson Joseph R. Peterson and Michael J. Hallock Luthey-Schulten
More informationOverview. kmcand DFT. An example to start with. An example to start with 03/09/2014
An Introduction of the Use of Monte Carlo Simulations in Catalysis. Tonek Jansen verview Kinetics and kinetic Monte Carlo The Master Equation Getting the Rate Constants Algorithms Examples Bangalore, 5
More informationCell Cycle Modeling for Budding Yeast with Stochastic Simulation Algorithms
Cell Cycle Modeling for Budding Yeast with Stochastic Simulation Algorithms Tae-Hyuk Ahn 1, Layne T Watson 1,2, Yang Cao 1, Clifford A Shaffer 1, and William T Baumann 3 1 Department of Computer Science,
More informationBasic modeling approaches for biological systems. Mahesh Bule
Basic modeling approaches for biological systems Mahesh Bule The hierarchy of life from atoms to living organisms Modeling biological processes often requires accounting for action and feedback involving
More informationMethods of Data Analysis Random numbers, Monte Carlo integration, and Stochastic Simulation Algorithm (SSA / Gillespie)
Methods of Data Analysis Random numbers, Monte Carlo integration, and Stochastic Simulation Algorithm (SSA / Gillespie) Week 1 1 Motivation Random numbers (RNs) are of course only pseudo-random when generated
More informationKinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationChapter 11: CHEMICAL KINETICS
Chapter : CHEMICAL KINETICS Study of the rate of a chemical reaction. Reaction Rate (fast or slow?) Igniting gasoline? Making of oil? Iron rusting? We know about speed (miles/hr). Speed Rate = changes
More informationwhere a + b = 2 (this is the general case) These all come from the fact that this is an overall second order reaction.
Chapter 7 Problems Page of 6 //007 7. Hydrolysis of ethyl acetate is as follows: EtAc + OH - Ac - + EtOH. At 5 ºC, the disappearance of OH - is used to determine the extent of the reaction, leading to
More information1. Introduction to Chemical Kinetics
1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism
More informationAdaptive Deployment of Model Reductions for Tau-leaping Simulation
Adaptive Deployment of Model Reductions for Tau-leaping Simulation Sheng Wu, 1, a) Jin Fu, 1, b) 1, c) and Linda R. Petzold Department of Computer Science, University of California Santa Barbara, Santa
More informationChemical Kinetics. Reaction Mechanisms
Chemical Kinetics Kinetics is a study of the rate at which a chemical reaction occurs. The study of kinetics may be done in steps: Determination of reaction mechanism Prediction of rate law Measurement
More informationNumerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011
Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationGene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back
Gene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back June 19, 2007 Motivation & Basics A Stochastic Approach to Gene Expression Application to Experimental
More informationSPA for quantitative analysis: Lecture 6 Modelling Biological Processes
1/ 223 SPA for quantitative analysis: Lecture 6 Modelling Biological Processes Jane Hillston LFCS, School of Informatics The University of Edinburgh Scotland 7th March 2013 Outline 2/ 223 1 Introduction
More informationChemical Kinetics. Introduction
2 / 152 Chemical Kinetics Copyright c 2018 by Nob Hill Publishing, LLC 1 / 152 Introduction The purpose of this chapter is to provide a framework for determining the reaction rate given a detailed statement
More informationNumerical Simulation for Biochemical Kinetics
Chapter 1 Numerical Simulation for Biochemical Kinetics Daniel T. Gillespie and Linda R. Petzold In chemical systems formed by living cells, the small numbers of molecules of a few reactant species can
More informationAnatoly B. Kolomeisky. Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS?
Anatoly B. Kolomeisky Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS? Motor Proteins Enzymes that convert the chemical energy into mechanical
More informationA Simple Protein Synthesis Model
A Simple Protein Synthesis Model James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 3, 213 Outline A Simple Protein Synthesis Model
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech Version 1.0b, September 14, 2014 c 2014 by Princeton University Press All rights reserved. This is the electronic edition
More informationPartial Differential Equations (PDEs)
10.34 Numerical Methods Applied to Chemical Engineering Fall 2015 Final Exam Review Partial Differential Equations (PDEs) 1. Classification (a) Many PDEs encountered by chemical engineers are second order
More informationStiffness in stochastic chemically reacting systems: The implicit tau-leaping method
JOURNAL OF CHEICAL PHYSICS VOLUE 119, NUBER 24 22 DECEBER 2003 Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method uruhan Rathinam a) Department of athematics and Statistics,
More informationResearch Article Adaptive Time-Stepping Using Control Theory for the Chemical Langevin Equation
Applied athematics Volume 215, Article ID 567275, 1 pages http://dx.doi.org/1.1155/215/567275 Research Article Adaptive Time-Stepping Using Control Theory for the Chemical Langevin Equation Silvana Ilie
More informationCatalysis at the Sub-Nanoscale: Complex CO Oxidation Chemistry on a Few Au Atoms
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 214 Catalysis at the Sub-Nanoscale: Complex CO Oxidation Chemistry on a Few Au
More informationStochastic Gene Expression: Modeling, Analysis, and Identification
Munsky; q-bio Stochastic Gene Expression: Modeling, Analysis, and Identification Mustafa Khammash University of California, Santa Barbara Center for Control, Dynamical Systems and Computation CC DC Outline
More informationRenormalization Group: non perturbative aspects and applications in statistical and solid state physics.
Renormalization Group: non perturbative aspects and applications in statistical and solid state physics. Bertrand Delamotte Saclay, march 3, 2009 Introduction Field theory: - infinitely many degrees of
More informationMultiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models
Multiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models Russel Caflisch 1, Mark Gyure 2, Bo Li 4, Stan Osher 1, Christian Ratsch 1,2, David Shao 1 and Dimitri Vvedensky 3 1 UCLA,
More informationFrom cell biology to Petri nets. Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE
From cell biology to Petri nets Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE Biology = Concentrations Breitling / 2 The simplest chemical reaction A B irreversible,
More informationReaction Mechanisms. Chemical Kinetics. Reaction Mechanisms. Reaction Mechanisms. Reaction Mechanisms. Reaction Mechanisms
Chemical Kinetics Kinetics is a study of the rate at which a chemical reaction occurs. The study of kinetics may be done in steps: Determination of reaction mechanism Prediction of rate law Measurement
More informationDiffusion. CS/CME/BioE/Biophys/BMI 279 Nov. 15 and 20, 2016 Ron Dror
Diffusion CS/CME/BioE/Biophys/BMI 279 Nov. 15 and 20, 2016 Ron Dror 1 Outline How do molecules move around in a cell? Diffusion as a random walk (particle-based perspective) Continuum view of diffusion
More informationA(γ A D A + γ R D R + γ C R + δ A )
Title: Mechanisms of noise-resistance in genetic oscillators Authors: José M. G. Vilar 1,2, Hao Yuan Kueh 1, Naama Barkai 3, and Stanislas Leibler 1,2 1 Howard Hughes Medical Institute, Departments of
More information