STOCHASTIC CHEMICAL KINETICS
|
|
- Alan Arnold
- 5 years ago
- Views:
Transcription
1 STOCHASTIC CHEICAL KINETICS Dan Gillespie Current Support: Caltech (NIGS) Caltech (NIH) University of California at Santa Barbara (NIH) Past Support: Caltech (DARPA/AFOSR, Beckman/BNC)) University of California at Santa Barbara (DOE) olecular Sciences Institute (Sandia/DOE) Office of Naval Research ain Collaborators: Linda Petzold (UCSB) John Doyle (Caltech) A Chemically Reacting System consists of olecules of N chemical species S, S., N - Inside a volume Ω, at some temperature T. elemental reaction channels R, R., - We assume R to be a single instantaneous physical event that changes the population of at least one species. - In practice, elemental means that R must be one of two types: Unimolecular: Si something else, or Bimolecular: S + S something else. i i - All other types of reaction (trimolecular, reversible, etc.) are made up of a series of two or more elemental reactions. E.g.: S+ S2 S2 S+ S2+ S3 S4+ S5 is typically S + S S + S
2 Question: How does a spatially homogeneous (or well-stirred) chemically reacting system evolve in time? The Trad Answer: According to the reaction rate equation (RRE). A set of coupled, first-order ODEs. Derived using ad hoc, phenomenological reasoning. - Is more than the mass action equations of thermodynamics, which apply only to systems in complete equilibrium. Implies the system evolves continuously and deterministically. - Yet molecules come in integer numbers and react stochastically. Is empirically accurate for large (test tube size) systems But is sometimes not adequate for very small (cell-size) systems. * * * Doing it right : olecular Dynamics The most exact way of describing the system s evolution. The motion picture approach: Tracks the position and velocity of every molecule in the system. Simulates every collision, non-reactive as well as reactive. Shows changes in species populations and their spatial distributions. But... it s unfeasibly slow for nearly all realistic systems. 2
3 A great simplification occurs if successive reactive collisions tend to be separated in time by very many non-reactive collisions. The overall effect of the non-reactive collisions is to randomize - the velocities of the molecules (axwell-boltzmann distribution). - the positions of the molecules (spatially uniform or well-stirred), Then, instead of having to describe the system s state as the position, velocity and species of each molecule, we need only give ( ) X( t) X ( t),, X ( t), X i ( t) the number of S i molecules at time t. N But this well-stirred simplification, which... ignores the non-reactive collisions, drastically truncates the definition of the system s state,... comes at a price: X ( t) must now be viewed as a stochastic process. But in fact, the system was never deterministic to begin with. Even if molecules moved according to classical mechanics... - Unimolecular reactions always involve randomness (Q). - Bimolecular reactions usually do too. - A system of many colliding molecules is so sensitive to initial conditions that, for all practical purposes, it evolves randomly. - The system is not isolated. It s in a heat bath, which keeps it at temperature T via essentially random interactions. 3
4 For well-stirred systems, each R is completely characterized by a propensity function a ( x ) : Given the system is in state x, a ( x) dt probability that one R event will occur in the next dt. - The existence and form of a ( x ) follow from molecular physics. a state change vector ( ν,, vn ) - R induces Examples: ν i the change in ν : x x+ ν. { i } X i caused by one R event. ν the stoichiometric matrix. ( ) x ( ) c 2 : ν =,,,, ; ( ) = ( x) = S S a dt c dt x a c x ( ) c : same ν ; ( x) = 2 ( x) = 2 S S S a dt c dt x x a c x x 2 ( r2 )( v2 dt) Prob{ collision in dt} = π. Prob{ R collision } p. Ω 2 2 ( π r2 )( v2dt ) π r2 v2 p a ( x) dt p xx2 = xx2 dt= c xx2 dt Ω Ω a ( ) Prob that a typical S -S x 2 c pair reacts in next dt E R iff " collisional K. E." > E p = exp... Arrhenius! kbt Diffusional motion (well-stirred): π r v is replaced by 4 π r2 ( D + D2 ) 4
5 Two exact, rigorously derivable consequences.... The chemical master equation (CE): P( x, t x, t ) = a ( ) P(, t, t) a ( ) P(, t, t) t x x x x x x = ν ν. P( x, t x, t ) Prob { X( t) = x, given X( t ) = x } for t t. Follows from the probability statement P( x, t+ dt x, t) = P( x, t x, t) ( a ( x) dt) = = But the CE is usually too hard to solve. ( ν ) + P( x ν, t x, t ) a ( x ) dt. Averages: f ( X t ) f x P x t x t x ( ) ( ) (,, ). If we multiply the CE through by x and then sum over x, we find d X( t) = ν a ( X( t) ). dt = If there were no fluctuations, ( ( )) = ( ( ) ) = ( ( )) a X t a X t a X t, and the above would reduce to: dx( t) dt = = ( X( t) ) ν a. - This is the reaction-rate equation (RRE). - It s usually written in terms of the concentration Z( t) X( t) Ω. But as yet, we have no ustification for ignoring fluctuations. 5
6 2. The stochastic simulation algorithm (SSA): A procedure for constructing sample paths or realizations of X ( t). Idea: Generate properly distributed random numbers for - the time τ to the next reaction, - the index of that reaction. p( τ, x, t) dτ probability, given X( t ) = x, that the next reaction will occur in [ t+ τ, t+ τ+ dτ ), and will be R. = P ( τ ) a ( x ) dτ, P ( τ ) Pr( no reactions in time τ ). ( ) P ( τ + dτ ) = P ( τ ) a( x ) dτ, where a ( x) a ( ) x. dp ( τ ) a Implies ( x ) τ = a( x ) P ( τ ). Solution: P ( τ ) = e. dτ a ( x) τ a ( x) τ a ( x) p( τ, x, t) = e a ( x) = a ( x)e a( x ). Thus, - τ is an exponential random variable with mean a( x ), - is an integer random variable with probabilities a ( x) a( x ). The Direct Version of the SSA. In state x at time t, evaluate a ( x ),, a ( ) x, and a ( x) a ( x ). = 2. Draw two unit-interval uniform random numbers r and r 2, and compute τ and according to τ = ln, a( x) r = the smallest integer satisfying ak ( x) > r2 a( x ). 3. Replace t t+ τ and x x+ ν. k= 4. Record ( x, t). Return to Step, or else end the simulation. 6
7 c A Simple Example: S. a ( x ) = c x, ν =. Take X () x dx( t) RRE: = c X( t). Solution is X( t) = x e c t. dt CE: =. P( x, t x,) = c ( x + ) P( x+, t x,) xp( x, t x,) t. x! c x t c t x x = = x!( x x )! Solution: ( ) P( x, t x,) e e ( x,,, x ) which implies X ( t) x e c t ct ct =, sdev { X( t) } x e ( e ) =. SSA: Given X( t) = x, generate τ = ln, then update: c x r t t+ τ, x x. 8 S c =, X () =
8 a ( x) dt Prob that R will fire in next dt CE SSA The SSA... Is exact. It does not entail approximating dt by t. Is procedurally simple, even when the CE is intractable. Comes in a variety of implementations - Direct ethod (Gillespie, 976) - First Reaction ethod (Gillespie, 976) - Next Reaction ethod (Gibson & Bruck, 2) - First Family ethod (Lok, 23) - odified Direct ethod (Cao, Li & Petzold, 24) - Sorting Direct ethod (ccollum, et al. 26) Remains too slow for most practical problems: Simulating every reaction event one at a time ust takes too much time if any reactants are present in very large numbers. 8
9 We would be willing to sacrifice a little exactness if that would buy us a faster simulation. Tau-Leaping Approximately advances the process by a pre-selected time τ, which may encompass more than one reaction event. Key: The Poisson random variable with mean aτ can be defined: P ( aτ ) the number of events that will occur in a time τ, when the probability of an event in any dt is adt, provided a is a positive constant. With X( t ) = x, let us choose τ small enough to satisfy the Leap Condition: Each a ( x ) constant in [ t, t+ τ ]. Then the number of R firings in [ t, t τ ] = + will be P ( a ( x) τ) ( a τ) X( t τ ) x P ( x) + + ν. So = ( a τ) X( t τ ) x P ( x) + + ν - Practical Considerations for Implementing Tau-Leaping - Finding the largest τ that satisfies the Leap Condition. - Accomplished via an accuracy control parameter ε. - We estimate the largest τ for which a a ε,. Avoiding negative populations, and segueing to the SSA. - Accomplished via a second control parameter n c. - We call any reaction that is within c n firings of exhausting any reactant a critical reaction. Then we tau-leap no farther than the next firing of a critical reaction. - Becomes the SSA when all reactions are classified critical. τ 9
10 5 45 monomer X, unstable dimer X2, stable dimer X X2 X3 S--> c = S + S --> S2 c2 =.2 S2 --> S + S c3 =.5 S2 --> S3 c4 =.4 - Exact SSA Run. - Intially: X=,; X2=X3=. - 5 reactions per plotted dot - 57,67 reactions total. 5 X Time 5 45 monomer X, unstable dimer X2, stable dimer X X2 X3 S--> c = S + S --> S2 c2 =.2 S2 --> S + S c3 =.5 S2 --> S3 c4 =.4 - Explicit Tau Leaping Run - ε =.4; n c =. - leap per plotted dot leaps total. - Run time speedup over SSA > X. 5 X Time
11 a ( x) dt Prob that R will fire in next dt { a ( x ) const over τ, } CE SSA Tau-Leaping Discrete & Stochastic Speeding up Tau-Leaping: The Langevin Equation Two math facts: - If m, then P( m) N ( m, m). - 2 N ( m, σ ) = m+ σn (,). So, with X( t ) = x, suppose we can choose τ small enough to satisfy the Leap Condition, yet also large enough that a ( x) τ,. Then... X( t+ τ ) x+ P ( a ( x) τ) = + = ν ( a ( ) τ, a ( ) τ) x N x x = ν x+ a ( ) τ + a ( ) τ (,) x x N ν t+ + ν a ( ) + ν a ( ) X ( τ ) x x τ x N (,) τ. = =
12 ( ) ν ( ) X ( t+ τ ) x + ν a x τ + a x N (,) τ = = This is the Langevin leaping formula. It s faster than the ordinary tau-leaping formula, because - a ( ) τ x means lots of reaction events get leapt over in τ ; - normal random numbers can be generated faster than Poissons. It directly implies, and is entirely equivalent to, a SDE called the chemical Langevin equation (CLE): dx( t) ν a ( X( t) ) + ν a ( X( t) ) Γ ( t). dt = = N (,) - Gaussian white noise: Γ ( t) lim lim N, + +. dt dt dt dt - Satisfies Γ ( t) Γ ( t ) = δ δ ( t t ). Our discrete stochastic process X ( t) has now been approximated as a continuous stochastic process. a ( x) dt Prob that R will fire in next dt { a ( x ) const over τ, } CE SSA Tau-Leaping Discrete & Stochastic { a ( x ) τ, } CFPE * CLE Continuous & Stochastic * - J. Chem. Phys. 3:297 (2) - Am. J. Phys. 64:246 (996) - J. Phys. Chem. A 6:563 (22) 2
13 Def: All The Thermodynamic Limit X i, and Ω, with X i Ω constants. a = c x x In the thermodynamic limit, a 's grow like (system size). c a x x2 Ω = xx2 x2 all So in the thermodynamic limit, we see that in the CLE dx( t) ν a ( X( t) ) + ν a ( X( t) ) Γ ( t), dt = = - the deterministic term grows like (system size), - the stochastic term grows like (system size) /2. Rule of Thumb: Relative fluctuations die off as (system size) /2. At the thermodynamic limit the stochastic term disappears, leaving dx( t) ν a ( X( t) ) the RRE derived! dt = X ( t) has now become a continuous deterministic process. a ( x) dt Prob that R will fire in next dt { a ( x ) const over τ, } CE SSA Tau-Leaping Discrete & Stochastic { a ( x ) τ, } CFPE CLE Continuous & Stochastic X i, Ω X i Ω = const i, i RRE Continuous & Deterministic 3
14 A Spectrum of Descriptive athematical odes (for well-stirred systems) a ( x) const a ( x) τ X, Ω in [ t, t+ τ ], X i Ω const i i CE/SSA Tau- Leaping CLE / CFPE RRE exact discrete approximate approximate approximate discrete continuous stochastic stochastic stochastic continuous deterministic increasing numbers of molecules Another ulti-scale Problem Some reactions/species may be very fast, others very slow. Fast and slow are interconnected not easy to separate. Often manifests as dynamical stiffness, a known ODE problem. SSA still works, and is exact. But it s agonizingly slow. Tau-leaping remains accurate, but the Leap Condition restricts τ to the shortest (fastest) time scale of the system. So even it s too slow. One approach: Implicit Tau-Leaping - A stochastic adaptation of the implicit Euler method for ODEs. Another approach: The Slow-Scale Stochastic Simulation Algorithm - Skips over the fast reactions and simulates only the slow reactions, using specially modified propensity functions. An adaptation of the rapid equilibrium / quasi steady-state methods for RREs. 4
Stiffness 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 informationStochastic Simulation of Chemical Kinetics
Annu. Rev. Phys. Chem. 2007. 58:35 55 First published online as a Review in Advance on October 12, 2006 The Annual Review of Physical Chemistry is online at http://physchem.annualreviews.org This article
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 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 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 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 informationEfficient step size selection for the tau-leaping simulation method
THE JOURNAL OF CHEMICAL PHYSICS 24, 04409 2006 Efficient step size selection for the tau-leaping simulation method Yang Cao a Department of Computer Science, Virginia Tech, Blacksburg, Virginia 2406 Daniel
More informationStochastic Simulation.
Stochastic Simulation. (and Gillespie s algorithm) Alberto Policriti Dipartimento di Matematica e Informatica Istituto di Genomica Applicata A. Policriti Stochastic Simulation 1/20 Quote of the day D.T.
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 informationNumerical methods for stochastic simulation of biochemical systems
International Journal of Sciences and Techniques of Automatic control & computer engineering IJ-STA, Volume 4, N 2, December 2010, pp. 1268 1283. Numerical methods for stochastic simulation of biochemical
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 informationReaction time distributions in chemical kinetics: Oscillations and other weird behaviors
Introduction The algorithm Results Summary Reaction time distributions in chemical kinetics: Oscillations and other weird behaviors Ramon Xulvi-Brunet Escuela Politécnica Nacional Outline Introduction
More informationAccelerated stochastic simulation of the stiff enzyme-substrate reaction
THE JOURNAL OF CHEMICAL PHYSICS 123, 144917 2005 Accelerated stochastic simulation of the stiff enzyme-substrate reaction Yang Cao a Department of Computer Science, University of California, Santa Barbara,
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 informationApproximate accelerated stochastic simulation of chemically reacting systems
JOURNAL OF CHEICAL PHYSICS VOLUE 115, NUBER 4 22 JULY 2001 Approximate accelerated stochastic simulation of chemically reacting systems Daniel T. Gillespie a) Research Department, Code 4T4100D, Naval Air
More informationAvoiding negative populations in explicit Poisson tau-leaping
THE JOURNAL OF CHEMICAL PHYSICS 123, 054104 2005 Avoiding negative populations in explicit Poisson tau-leaping Yang Cao a Department of Computer Science, University of California, Santa Barbara, Santa
More informationAccelerated Stochastic Simulation of the Stiff Enzyme-Substrate Reaction
JCP A5.05.047 1 Accelerated Stochastic Simulation of the Stiff Enzyme-Substrate Reaction Yang Cao a) Dept. of Computer Science, Univ. of California, Santa Barbara, Santa Barbara, CA 9106 Daniel T. Gillespie
More informationAccelerated Stochastic Simulation of the Stiff Enzyme-Substrate Reaction
1 Accelerated Stochastic Simulation of the Stiff Enzyme-Substrate Reaction Yang Cao a) Dept. of Computer Science, Univ. of California, Santa Barbara, Santa Barbara, CA 9106 Daniel T. Gillespie b) Dan T
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 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 informationA New Perspective on the Linear Noise Approximation
A New Perspective on the Linear Noise Approximation Journal: IET Systems Biology anuscript ID: Draft anuscript Type: Original anuscript Date Submitted by the Author: n/a Complete List of Authors: Wallace,
More informationThe slow-scale stochastic simulation algorithm
THE JOURNAL OF CHEICAL PHYSICS 22, 046 2005 The slow-scale stochastic simulation algorithm Yang Cao a) Department of Computer Science, University of California, Santa Barbara, Santa Barbara, California
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 informationExtending the Tools of Chemical Reaction Engineering to the Molecular Scale
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
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 informationThe multinomial simulation algorithm for discrete stochastic simulation of reaction-diffusion systems
THE JOURNAL OF CHEMICAL PHYSICS 130, 0944 09 The multinomial simulation algorithm for discrete stochastic simulation of reaction-diffusion systems Sotiria Lampoudi, 1,a Dan T. Gillespie, 2 and Linda R.
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 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 informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA February 14, 2017 Introduction When we choose to model
More informationc 2005 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol. 4, No. 3, pp. 867 895 c 2005 Society for Industrial and Applied Mathematics CONSISTENCY AND STABILITY OF TAU-LEAPING SCHEMES FOR CHEMICAL REACTION SYSTEMS MURUHAN RATHINAM,
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 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 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 informationREACTION KINETICS. Catalysts substances that increase the rates of chemical reactions without being used up. e.g. enzymes.
REACTION KINETICS Study of reaction rates Why? Rates of chemical reactions are primarily controlled by 5 factors: the chemical nature of the reactants 2 the ability of the reactants to come in contact
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 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 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 informationNumerical solution of stochastic models of biochemical kinetics
Numerical solution of stochastic models of biochemical kinetics SILVANA ILIE a WAYNE H. ENRIGHT b KENNETH R. JACKSON b a Department of Mathematics, Ryerson University, Toronto, ON, M5B 2K3, Canada email:
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 informationEffect of Reactant Size on Discrete Stochastic Chemical Kinetics
Effect of Reactant Size on Discrete Stochastic Chemical Kinetics Dan T. Gillespie Dan T Gillespie Consulting, 30504 Cordoba Place, Castaic, California 9384, USA Sotiria Lampoudi a Dept. of Computer Science,
More informationStochastic Simulation Methods for Solving Systems with Multi-State Species
Stochastic Simulation Methods for Solving Systems with Multi-State Species Zhen Liu Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of
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 informationCHEMISTRY. Chapter 14 Chemical Kinetics
CHEMISTRY The Central Science 8 th Edition Chapter 14 Kozet YAPSAKLI kinetics is the study of how rapidly chemical reactions occur. rate at which a chemical process occurs. Reaction rates depends on The
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 informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
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 informationChapter 14 Chemical Kinetics
7/10/003 Chapter 14 Chemical Kinetics 14-1 Rates of Chemical Reactions 14- Reaction Rates and Concentrations 14-3 The Dependence of Concentrations on Time 14-4 Reaction Mechanisms 14-5 Reaction Mechanism
More informationTHE SLOW-SCALE SSA. Yang Cao, Dan Gillespie, Linda Petzold
THE SLOW-SCALE SSA Yang Cao, Dan Gillepie, Linda Petzold With thank to: Muruhan Rathinam, Hana El-Samad, John Doyle and E. Haeltine & J. Rawling, J. Chem. Phy 117, 6959 (2002) Chri Rao & Adam Arkin, J.
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate heory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationGeneralized Chemical Langevin Equation for Chemical Reaction Systems
ATCH Communications in athematical and in Computer Chemistry ATCH Commun. ath. Comput. Chem. 78 (207) 28-305 ISSN 0340-6253 Generalized Chemical Langevin Equation for Chemical Reaction Systems Tao Li a,
More informationElementary Reactions
Elementary Reactions Elementary reactions occur in a single encounter Unimolecular: A Rate = k[a] Bimolecular: A + B Rate = k[a][b] Termolecular: A + B + C Rate = k[a][b][c] Termolecular reactions are
More informationChapter 14 Chemical Kinetics
Chapter 14 Chemical Kinetics Thermodynamics tells us what can happen and how far towards completion a reaction will proceed. Kinetics tells us how fast the reaction will go. Study of rates of reactions
More informationOn the Interpretation of Delays in Delay Stochastic Simulation of Biological Systems
On the Interpretation of Delays in Delay Stochastic Simulation of Biological Systems Roberto Barbuti Giulio Caravagna Andrea Maggiolo-Schettini Paolo Milazzo Dipartimento di Informatica, Università di
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 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 informationChapter 14 Chemical Kinetics
4//004 Chapter 4 Chemical Kinetics 4- Rates of Chemical Reactions 4- Reaction Rates and Concentrations 4-3 The Dependence of Concentrations on Time 4-4 Reaction Mechanisms 4-5 Reaction Mechanism and Rate
More information14.1 Expressing the Reaction Rate
14.1 Expressing the Reaction Rate Differential Rate and Integrated Rate Laws Fred Omega Garces Chemistry 201 Miramar College 1 Expressing Reaction Rates Outline Kinetics: Intro Factors influencing Reaction
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 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 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 information3: Chemical Kinetics Name: HW 6: Review for Unit Test KEY Class: Date: A Products
3: Chemical Kinetics Name: HW 6: Review for Unit Test KEY Class: Date: Page 1 of 9 AP Multiple Choice Review Questions 1 16 1. The reaction rate is defined as the change in concentration of a reactant
More informationStochastic model of mrna production
Stochastic model of mrna production We assume that the number of mrna (m) of a gene can change either due to the production of a mrna by transcription of DNA (which occurs at a rate α) or due to degradation
More informationMultiscale stochastic simulation algorithm with stochastic partial equilibrium assumption for chemically reacting systems
Journal of Computational Physics 206 (2005) 395 411 www.elsevier.com/locate/jcp Multiscale stochastic simulation algorithm with stochastic partial equilibrium assumption for chemically reacting systems
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationKINETICS CHEMICAL CHEMIC. Unit. I. Multiple Choice Questions (Type-I)
Unit 4 CHEMICAL CHEMIC KINETICS I. Multiple Choice Questions (Type-I) 1. The role of a catalyst is to change. gibbs energy of reaction. enthalpy of reaction. activation energy of reaction. equilibrium
More informationBinomial leap methods for simulating stochastic chemical kinetics
Binomial leap methods for simulating stochastic chemical kinetics Tianhai Tian and Kevin Burrage Citation: The Journal of Chemical Physics 121, 10356 (2004); doi: 10.1063/1.1810475 View online: http://dx.doi.org/10.1063/1.1810475
More informationKinetics. Consider an irreversible unimolecular reaction k. -d[a]/dt = k[a] Can also describe in terms of appearance of B.
Kinetic data gives insight into reaction mechanisms kinetic analysis will describe a relationship between concentrations of all chemical species before the rate determining step in a given reaction and
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 informationThe Kramers problem and first passage times.
Chapter 8 The Kramers problem and first passage times. The Kramers problem is to find the rate at which a Brownian particle escapes from a potential well over a potential barrier. One method of attack
More informationC H E M I C N E S C I
C H E M I C A L K I N E T S C I 4. Chemical Kinetics Introduction Average and instantaneous Rate of a reaction Express the rate of a reaction in terms of change in concentration Elementary and Complex
More informationChapter 15: Phenomena
Chapter 5: Phenomena Phenomena: The reaction A(aq) + B(aq) C(aq) was studied at two different temperatures (298 K and 350 K). For each temperature the reaction was started by putting different concentrations
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 2. Chemical Kinetics. Chemical Kinetics 6/26/11. One (elementary) step reaction
Lecture Chemical Kinetics 1 One (elementary) step reaction im i i M i is the number of species i, i are the stoichiometric coefficients i i Chemical Kinetics =0ifi is not a reactant =0ifi is not a product
More information2 Reaction kinetics in gases
2 Reaction kinetics in gases October 8, 2014 In a reaction between two species, for example a fuel and an oxidizer, bonds are broken up and new are established in the collision between the species. In
More informationBrown et al, Chemistry, 2nd ed (AUS), Ch. 12:
Kinetics: Contents Brown et al, Chemistry, 2 nd ed (AUS), Ch. 12: Why kinetics? What is kinetics? Factors that Affect Reaction Rates Reaction Rates Concentration and Reaction Rate The Change of Concentration
More informationFoundations of Chemical Kinetics. Lecture 18: Unimolecular reactions in the gas phase: RRK theory
Foundations of Chemical Kinetics Lecture 18: Unimolecular reactions in the gas phase: RRK theory Marc R. Roussel Department of Chemistry and Biochemistry Frequentist interpretation of probability and chemical
More informationon-line kinetics 3!!! Chemistry 1B Fall 2013
on-line kinetics 3!!! Chemistry 1B Fall 2013 1 on-line kinetics 3!!! Chemistry 1B Fall 2013 Mechanism of a chemical reaction Elementary reactions Activation energy and reaction coordinate diagram 2 Chemistry
More informationEffect of reactant size on 2D discrete stochastic chemical kinetics
Effect of reactant size on 2D discrete stochastic chemical kinetics Sotiria Lampoudi * Dept. of Computer Science, University of California, Santa Barbara, California 93106, USA Dan T. Gillespie Dan T Gillespie
More informationRRK theory. Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge. April 3, 2009
RRK theory Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge April 3, 2009 Isomerization reactions as a test of theories of unimolecular reactions Gas-phase unimolecular
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 informationChapter 12. Chemical Kinetics
Chapter 12 Chemical Kinetics Chapter 12 Table of Contents 12.1 Reaction Rates 12.2 Rate Laws: An Introduction 12.3 Determining the Form of the Rate Law 12.4 The Integrated Rate Law 12.5 Reaction Mechanisms
More informationExtending the Stochastic Simulation Software Package StochPy with Stochastic Delays, Cell Growth and Cell Division
Masters Project Extending the Stochastic Simulation Software Package StochPy with Stochastic Delays, Cell Growth and Cell Division Author: Maxim Moinat Bioinformatics and Systems Biology VU University
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 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 informationCHAPTER 9 LECTURE NOTES
CHAPTER 9 LECTURE NOTES 9.1, 9.2: Rate of a reaction For a general reaction of the type A + 3B 2Y, the rates of consumption of A and B, and the rate of formation of Y are defined as follows: Rate of consumption
More informationChemical Kinetics. What Influences Kinetics?
Chemical Kinetics Predictions of likelihood for a reaction to occur have been based on difference in energy between products and reactants: Thermodynamics only compares reactants to products, says nothing
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 informationTHREE-BODY INTERACTIONS DRIVE THE TRANSITION TO POLAR ORDER IN A SIMPLE FLOCKING MODEL
THREE-BODY INTERACTIONS DRIVE THE TRANSITION TO POLAR ORDER IN A SIMPLE FLOCKING MODEL Purba Chatterjee and Nigel Goldenfeld Department of Physics University of Illinois at Urbana-Champaign Flocking in
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 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 informationChapter 14 Chemical Kinetics
How fast do chemical processes occur? There is an enormous range of time scales. Chapter 14 Chemical Kinetics Kinetics also sheds light on the reaction mechanism (exactly how the reaction occurs). Why
More informationarxiv: v3 [cs.cc] 28 Jan 2009
Robust Stochastic Chemical Reaction Networks and Bounded Tau-Leaping arxiv:0803.1030v3 [cs.cc] 28 Jan 2009 David Soloveichik dsolov@caltech.edu Abstract The behavior of some stochastic chemical reaction
More informationTransmission in finite populations
Transmission in finite populations Juliet Pulliam, PhD Department of Biology and Emerging Pathogens Institute University of Florida and RAPIDD Program, DIEPS Fogarty International Center US National Institutes
More informationSingle cell experiments and Gillespie s algorithm
Single cell experiments and Gillespie s algorithm Fernand Hayot Department of Neurology, Mount Sinai School of Medicine, New York, NY, 10029 email: fernand.hayot@mssm.edu August 18, 2008 1 Introduction
More informationExpress the transition state equilibrium constant in terms of the partition functions of the transition state and the
Module 7 : Theories of Reaction Rates Lecture 33 : Transition State Theory Objectives After studying this Lecture you will be able to do the following. Distinguish between collision theory and transition
More informationReaction Mechanisms Dependence of rate on temperature Activation Energy E a Activated Complex Arrhenius Equation
Kinetics Dependence of rate on Concentration (RATE LAW) Reaction Mechanisms Dependence of rate on temperature Activation Energy E a Activated Complex Arrhenius Equation Mary J. Bojan Chem 112 1 A MECHANISM
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationarxiv: v2 [q-bio.mn] 31 Aug 2007
A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays David F. Anderson Department of Mathematics, University of Wisconsin-Madison, Madison, Wi 53706
More informationv AB + C à AC + B Ø Bonds breaking
Chemistry Study Guide 3 v Kinetics reaction rates Ø Catalyst Ø Temperature Ø Concentration Ø Bonds Ø Surface area v Kinetics Ø Kinetic energy is directly proportional to the temperature Ø Gasses will react
More informationFoundations of Chemical Kinetics. Lecture 12: Transition-state theory: The thermodynamic formalism
Foundations of Chemical Kinetics Lecture 12: Transition-state theory: The thermodynamic formalism Marc R. Roussel Department of Chemistry and Biochemistry Breaking it down We can break down an elementary
More information