Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
|
|
- Mildred Lamb
- 5 years ago
- Views:
Transcription
1 Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk
2 The Three Lectures 1 Gillespie s algorithm and its relation to SDEs 2 Tau-leaping and multilevel approximations to Gillespie 3 Monte Carlo/SDEs in mathematical finance Vienna Des Higham Gillespie 2 / 3
3 Overview of this lecture Gillespie Chemical Master Equation Chemical Langevin Equation Thermodynamic limit Mean hitting times Eric Renshaw (Strathclyde): Watch out for the boundaries! Vienna Des Higham Gillespie 3 / 3
4 Chemical Kinetics We have some chemical species undergoing some chemical reactions Not discussing molecular dynamics: not willing or able to measure initial location and velocity of every molecule, and/or not able to compute interactions between every pair of molecules Question of interest: How does the level of each species change over time? Vienna Des Higham Gillespie 4 / 3
5 Random? Inherently stochastic? Vienna Des Higham Gillespie 5 / 3
6 Random? Inherently stochastic? NO Vienna Des Higham Gillespie 5 / 3
7 Random? Inherently stochastic? Effectively stochastic? NO Vienna Des Higham Gillespie 5 / 3
8 Random? Inherently stochastic? Effectively stochastic? NO YES Vienna Des Higham Gillespie 5 / 3
9 Random? Inherently stochastic? Effectively stochastic? NO YES Justification for stochastic modelling... Vienna Des Higham Gillespie 5 / 3
10 Protein Toggle Model is from Gene regulatory networks: A coarse-grained equation-free approach to multiscale computation, by Erban, Kevrekidis, Adalsteinsson and Elston, Journal of Chemical Physics, 124, 200. Mutually repressive proteins: mass action (ODE) dz 1 dt = κz 1 ( γ 1 + ωz 2 2 δz 1 ) For certain realistic parameter values this ODE has two linearly stable steady states: one has z 1 (t) z a 481 and the other z 1 (t) z b Unrealistic that the cell s fate is completely specified by the initial condition, and a cell cannot switch dynamically between states. Stochastic simulation (Gillespie) gives... Vienna Des Higham Gillespie / 3
11 Protein Toggle Vienna Des Higham Gillespie / 3
12 Pure decay: S c=1, initially 10 molecules CME CLE RRE Vienna Des Higham Gillespie 8 / 3
13 Pure decay: S c=1, initially 100 molecules CME CLE RRE Vienna Des Higham Gillespie 9 / 3
14 S c=1, initially K molecules CME: p i (t) is prob. of having i molecules at time t d dt p i(t) = (i + 1)p i+1 (t) ip i (t), 0 i K. Process X(t) has E[X(t)] = Ke t and Var[X(t)] = Ke t (1 e t ) For the CLE, process is an SDE dy(t) = Y(t)dt Y(t) dw(t). This has same mean and variance as X(t) For the RRE, process is an ODE with solution Ke t Vienna Des Higham Gillespie 10 / 3
15 Chemical Kinetics (Gillespie 19, 2000) Well-stirred, thermal equilibrium, fixed volume Chemical Master Equation (CME) Species have integer values ODE for prob. of every possible state at time t Chemical Langevin Equation (CLE) Species have real values SDE for concentration of each species at time t Reaction Rate Equations (RRE) Species have real values ODE for concentration of each species at time t Vienna Des Higham Gillespie 11 / 3
16 Chemical Kinetics (Gillespie 19) N chemical species, M types of reaction (e.g. A + B C) State vector X 1 (t) X 2 (t) X(t) =., X(0) = X 0 X N (t) Each reaction, 1 j M, is described by a stoichiometric vector ν j R N such that X(t) X(t) + ν j, a propensity function a j (X(t)) such that the prob. of this reaction taking place over time [t, t + dt) is a j (X(t)) dt Vienna Des Higham Gillespie 12 / 3
17 Chemical Master Equation Discrete state space, continuous time Markov chain. Let P(x, t) be the prob. that X(t) = x dp(x, t) dt = M (a j (x ν j )P(x ν j, t) a j (x)p(x, t)) j=1 (Forward Kolmogorov equation) Gillespie s stochastic simulation algorithm gives a way to compute realisations of (t, X(t)) Also called Bortz-Kalos-Lebowitz algorithm, discrete event simulation, dynamic Monte Carlo, kinetic Monte Carlo,... Takes account of every reaction expensive Vienna Des Higham Gillespie 13 / 3
18 Chemical Langevin Equation Poisson with large mean normal SDE in R N dy(t) = M ν j a j (Y(t))dt + j=1 M ν j a j (Y(t)) dw j (t) Euler Maruyama computes approximate realisations of (t, Y(t)). j=1 (Switching off the noise gives the Reaction Rate Equation) Vienna Des Higham Gillespie 14 / 3
19 Remarks: the good The CLE is a rare example where an SDE arises naturally from fundamental modelling arguments The CLE captures the fluctuations of Poisson processes about their mean accurate when the means are large The CLE gets the first and second moments exactly right for first order reaction networks We can develop hybrid/multi-scale models that combine discrete and continuous regimes In the thermodynamic limit, over a finite time [0, T], the CLE and CME converge to the deterministic RRE; but the CLE and CME are much closer to each other than to the RRE (Kurtz, 190s) Vienna Des Higham Gillespie 15 / 3
20 Remarks: the bad Looking at the modelling/approximation arguments, we see that the CLE does not make sense for small population sizes Consequently, in general the CLE solution only exists up to a stopping time In particular, we can t generally talk about the steady state of the CLE Vienna Des Higham Gillespie 1 / 3
21 Reversible Isometry c 1 X 1 X 1 X 2 c 2 X 2 X 1 (0) + X 2 (0) = V Gillespie, The chemical Langevin and Fokker Planck equations for the reversible isomerization reaction, J. Phys. Chem. A, 2002 For c 1 = c 2, CME gives binomial steady state for X 1 CLE for X 1 : dy = ( c 1 Y + c 2 (V Y)) dt c 1 YdW 1 + c 2 (V Y)dW 2 Gillespie solves the steady Fokker-Planck equation and argues that Y has a Gaussian steady state with the correct mean and variance Vienna Des Higham Gillespie 1 / 3
22 c 1 = c 2 steady state distributions? Vienna Des Higham Gillespie 18 / 3
23 Reflect on this? Recall that V denotes the number of molecules If we add reflecting boundary conditions to the CLE then we obtain a Gaussian shaped steady state, restricted to the interval [0, V] This reflected CLE steady state has the correct mean, and the variance matches the CME up to O(V 2 e V ) Vienna Des Higham Gillespie 19 / 3
24 Reflected CLE Let q( ) be the steady state probability for the CME Let p( ) be the steady state probability for the reflected CLE We can show that Also p(0) 2 q(0) = p(v/2) q(v/2) = 1 + O(1/V) 1 Vπ e( 2 + log e 2)V ( 1 + O For V = 50, we have log Vπ e( 2 e 2)V ( )) e V which agrees with Gillespie s figure (similarly for V = 500) Vienna Des Higham Gillespie 20 / 3
25 Modifying the CLE? Wilkie & Wong Positivity preserving chemical Langevin equations, Chem. Phys., 2008 Idea: delete the offending diffusion coefficients in dy 1 = ( c 1 Y 1 + c 2 Y 2 ) dt c 1 Y 1 dw 1 + c 2 Y 2 dw 2 dy 2 = (c 1 Y 1 c 2 Y 2 ) dt + c 1 Y 1 dw 1 c 2 Y 2 dw 2 But this removes all the approximation power... Vienna Des Higham Gillespie 21 / 3
26 Modifying the CLE? 1 15 V = 1 Master Equation Modified Langevin Mass Action ODE 1 15 V = 10 Master Equation Modified Langevin Mass Action ODE <X 1 2 > 11 <X 1 2 > t t Vienna Des Higham Gillespie 22 / 3
27 Focus now on mean hitting time b x a T(x) := inf (t : Z(t) = a or Z(t) = b, given Z(0) = x) Vienna Des Higham Gillespie 23 / 3
28 Mean hitting time E[T(x)] a x b Vienna Des Higham Gillespie 24 / 3
29 Markov jump/birth & death process, Z(t) Discrete states {0, 1, 2,...,M}, with 0 and M absorbing: P (Z(t + h) = i + 1 Z(t) = i) = B i h + o(h) P (Z(t + h) = i 1 Z(t) = i) = D i h + o(h) P (Z(t + h) = i Z(t) = i) = 1 (B i + D i )h + o(h) Here, B 0 = D 0 = B M = D M = 0 Starting at state Z(0) = j, the expected time to be absorbed into state 0 or M is given by U j, where 2 4 (B 1 + D 1 ) B 1 D 2 (B 2 + D 2 ) B BM 2 D M 1 (B M 1 + D M 1 ) U 1 U 2 U U M = Vienna Des Higham Gillespie 25 / 3
30 Numerical Analysis Viewpoint Linear system can be written, for 1 i M 1, B i + D i 2 (U i+1 2U i + U i 1 ) + (B i D i ) U i+1 U i 1 2 Standard finite differences on the 2 point BVP ODE B(x) + D(x) u (x) + (B(x) D(x)) u (x) = 1 2 with u(a) = u(b) = 0 Here b a = M and we have x = 1 Interesting regime is M = 1 Vienna Des Higham Gillespie 2 / 3
31 Diffusion Approximation SDE: dy(t) = (B(y(t)) D(y(t))) dt + B(y(t)) dw 1 (t) D(y(t)) dw 2 (t) Let w(x) := E[T(x)] be the expected first time to hit a or b, given that y(0) = x Then w(x) satisfies the same 2 point BVP ODE Want to show that this ODE converges to the finite difference scheme Focus on specific examples... Vienna Des Higham Gillespie 2 / 3
32 Production from a source k X B i = k and D i = 0 Mean hitting times: Jump process b x k Diffusion process [ ( 1 e 2x + e 2a e 2b ( e 2b e 2x) ) + b x k e 2b + e 2a 2 ) + (1 e 2x + e (a 2a x + e 2a ( e 2x e 2a))] e 2b + e 2a 2 Vienna Des Higham Gillespie 28 / 3
33 Convergence: fix a = 0 and let b With x = αb for fixed α (0, 1), we have b min{2(1 α),α} Jump Diffusion C e where C is independent of b Example: k = 5, a = 0, α = 1 2 : Absolute Difference in Exit Time b Vienna Des Higham Gillespie 29 / 3
34 Production c X X B i = i and D i = 0 Mean hitting times: Jump process Diffusion process 1 c b 1 s=x ( 1 e 2x e 2a b ( e 2b c e 2b e 2a x + [1 e 2x e 2a e 2b e 2a 1 s e 2l ) dl + ln b ln x l x ](ln a ln x + e 2a a e 2l l )) dl Vienna Des Higham Gillespie 30 / 3
35 Convergence With x = αb for fixed α (0, 1), we have Jump (with a = 0) lim aց0 Diffusion Cb 2 where C is independent of b Proof Uses the expansions n 1 s = ln n + γ + 1 2n + O(n 2 ), as n, s=1 x e t t dt = ln x + γ + o(1), as x 0, where the Euler-Mascheroni constant γ = , and x ( e t dt = ex ) t x x + O(x 2 ), as x Vienna Des Higham Gillespie 31 / 3
36 Example, c = 5, a = 10 3, α = Absolute Difference in Exit Time b Vienna Des Higham Gillespie 32 / 3
37 Degradation X c X B i = 0 and D i = i Mean hitting times: Jump process Diffusion process ( e 2x e 2a 1 b (e 2b e 2b e 2a c + [1 e2x e 2a e 2b e 2a] x ( 1 c 1 c x s=a+1 1 s )) e 2l dl ln b + ln x l x (ln x ln a e 2a a e 2l l )) dl Vienna Des Higham Gillespie 33 / 3
38 Convergence With x = αb for fixed α (0, 1), we have lim lim ln 2 (Jump Diffusion) = aց0 b c Proof Uses asymptotic expansions for E 1 (x) = e t x t dt, x > 0 Note: the actual hitting times grow like ln(b), so we have relative convergence like O(1/ ln b) Vienna Des Higham Gillespie 34 / 3
39 Example, c = 5, α = 1 2, a = 10 2, 10 4, Absolute Difference in Exit Time a=1e 2 a=1e 4 a=1e b ln 2 5 = Vienna Des Higham Gillespie 35 / 3
40 Key Points Three modeling regimes: CME, CLE, RRE CME is most accurate, but expensive Thermodynamic limit can help to justify the use of CLE and RRE, but must be used with care Finding cheaper alternatives to full Gillespie/CME simulations is a big area We will look at tau-leaping (Euler approximation) multi-level Monte Carlo Vienna Des Higham Gillespie 3 / 3
Simulation 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 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 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 informationPersistence and Stationary Distributions of Biochemical Reaction Networks
Persistence and Stationary Distributions of Biochemical Reaction Networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Discrete Models in Systems Biology SAMSI December
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 informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
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 informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
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 informationProblems 5: Continuous Markov process and the diffusion equation
Problems 5: Continuous Markov process and the diffusion equation Roman Belavkin Middlesex University Question Give a definition of Markov stochastic process. What is a continuous Markov process? Answer:
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 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 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 informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationBiomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems
Biomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems Processing information to coordinate activity regulatory & signaling networks
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 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 information1.3 Forward Kolmogorov equation
1.3 Forward Kolmogorov equation Let us again start with the Master equation, for a system where the states can be ordered along a line, such as the previous examples with population size n = 0, 1, 2,.
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
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 informationMean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationLecture on Parameter Estimation for Stochastic Differential Equations. Erik Lindström
Lecture on Parameter Estimation for Stochastic Differential Equations Erik Lindström Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds +
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 informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
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 information6.2.2 Point processes and counting processes
56 CHAPTER 6. MASTER EQUATIONS which yield a system of differential equations The solution to the system of equations is d p (t) λp, (6.3) d p k(t) λ (p k p k 1 (t)), (6.4) k 1. (6.5) p n (t) (λt)n e λt
More informationGene regulatory networks: A coarse-grained, equation-free approach to multiscale computation
THE JOURNAL OF CHEMICAL PHYSICS 124, 084106 2006 Gene regulatory networks: A coarse-grained, equation-free approach to multiscale computation Radek Erban a Mathematical Institute, University of Oxford,
More informationAdvanced computational methods X Selected Topics: SGD
Advanced computational methods X071521-Selected Topics: SGD. In this lecture, we look at the stochastic gradient descent (SGD) method 1 An illustrating example The MNIST is a simple dataset of variety
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
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 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 informationTable of Contents [ntc]
Table of Contents [ntc] 1. Introduction: Contents and Maps Table of contents [ntc] Equilibrium thermodynamics overview [nln6] Thermal equilibrium and nonequilibrium [nln1] Levels of description in statistical
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 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 informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationOn the fast convergence of random perturbations of the gradient flow.
On the fast convergence of random perturbations of the gradient flow. Wenqing Hu. 1 (Joint work with Chris Junchi Li 2.) 1. Department of Mathematics and Statistics, Missouri S&T. 2. Department of Operations
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 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 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 informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
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 informationLecture 5: Moment generating functions
Lecture 5: Moment generating functions Definition 2.3.6. The moment generating function (mgf) of a random variable X is { x e tx f M X (t) = E(e tx X (x) if X has a pmf ) = etx f X (x)dx if X has a pdf
More informationWeakly interacting particle systems on graphs: from dense to sparse
Weakly interacting particle systems on graphs: from dense to sparse Ruoyu Wu University of Michigan (Based on joint works with Daniel Lacker and Kavita Ramanan) USC Math Finance Colloquium October 29,
More informationUncertainty quantification and systemic risk
Uncertainty quantification and systemic risk Josselin Garnier (Université Paris Diderot) with George Papanicolaou and Tzu-Wei Yang (Stanford University) February 3, 2016 Modeling systemic risk We consider
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationLecture 1 Modeling in Biology: an introduction
Lecture 1 in Biology: an introduction 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 informationMath 180B Homework 9 Solutions
Problem 1 (Pinsky & Karlin, Exercise 5.1.3). Let X and Y be independent Poisson distributed random variables with parameters α and β, respectively. Determine the conditional distribution of X, given that
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 informationFirst Passage Time Calculations
First Passage Time Calculations Friday, April 24, 2015 2:01 PM Homework 4 will be posted over the weekend; due Wednesday, May 13 at 5 PM. We'll now develop some framework for calculating properties of
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More information1 Elementary probability
1 Elementary probability Problem 1.1 (*) A coin is thrown several times. Find the probability, that at the n-th experiment: (a) Head appears for the first time (b) Head and Tail have appeared equal number
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 informationA random perturbation approach to some stochastic approximation algorithms in optimization.
A random perturbation approach to some stochastic approximation algorithms in optimization. Wenqing Hu. 1 (Presentation based on joint works with Chris Junchi Li 2, Weijie Su 3, Haoyi Xiong 4.) 1. Department
More informationContinuum Limit of Forward Kolmogorov Equation Friday, March 06, :04 PM
Continuum Limit of Forward Kolmogorov Equation Friday, March 06, 2015 2:04 PM Please note that one of the equations (for ordinary Brownian motion) in Problem 1 was corrected on Wednesday night. And actually
More informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationSensitivity analysis of complex stochastic processes
Sensitivity analysis of complex stochastic processes Yannis Pantazis join work with Markos Katsoulakis University of Massachusetts, Amherst, USA Partially supported by DOE contract: DE-SC0002339 February,
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
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 informationNon-nested Adaptive Timesteps in Multilevel Monte Carlo Computations
Non-nested Adaptive Timesteps in Multilevel Monte Carlo Computations M.B. Giles, C. Lester, and J. Whittle Abstract This paper shows that it is relatively easy to incorporate adaptive timesteps into multilevel
More informationChemical reaction network theory for stochastic and deterministic models of biochemical reaction systems
Chemical reaction network theory for stochastic and deterministic models of biochemical reaction systems University of Wisconsin at Madison anderson@math.wisc.edu MBI Workshop for Young Researchers in
More informationStochastic Modelling of Biological Processes Lecture Notes
Stochastic Modelling of Biological Processes Lecture Notes Ruth Baker Hilary Term 218 Abstract These lecture notes have been written to accompany the Part B course Stochastic Modelling of Biological Processes.
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 informationM4A42 APPLIED STOCHASTIC PROCESSES
M4A42 APPLIED STOCHASTIC PROCESSES G.A. Pavliotis Department of Mathematics Imperial College London, UK LECTURE 1 12/10/2009 Lectures: Mondays 09:00-11:00, Huxley 139, Tuesdays 09:00-10:00, Huxley 144.
More informationSimulation of Stochastic Kinetic Models
Simulation of Stochastic Kinetic Models Andrew Golightly & Colin S. Gillespie School of Mathematics & Statistics, Newcastle University, UK A growing realisation of the importance of stochasticity in cell
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationTesting the transition state theory in stochastic dynamics of a. genetic switch
Testing the transition state theory in stochastic dynamics of a genetic switch Tomohiro Ushikubo 1, Wataru Inoue, Mitsumasa Yoda 1 1,, 3, and Masaki Sasai 1 Department of Computational Science and Engineering,
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 informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationStochastic Integration (Simple Version)
Stochastic Integration (Simple Version) Tuesday, March 17, 2015 2:03 PM Reading: Gardiner Secs. 4.1-4.3, 4.4.4 But there are boundary issues when (if ) so we can't apply the standard delta function integration
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
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 informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
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 informationSTOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY
STOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY Third edition N.G. VAN KAMPEN Institute for Theoretical Physics of the University at Utrecht ELSEVIER Amsterdam Boston Heidelberg London New York Oxford Paris
More informationStochastic Modelling in Climate Science
Stochastic Modelling in Climate Science David Kelly Mathematics Department UNC Chapel Hill dtbkelly@gmail.com November 16, 2013 David Kelly (UNC) Stochastic Climate November 16, 2013 1 / 36 Why use stochastic
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 informationChapter 2. Brownian dynamics. 2.1 Introduction
Chapter 2 Brownian dynamics 2.1 Introduction The equations of physics are deterministic rather than stochastic in nature. Stochastic differential equations are used to approximate reality. They are introduced
More informationMarkov chain Monte Carlo algorithms for SDE parameter estimation
Markov chain Monte Carlo algorithms for SDE parameter estimation Andrew Golightly and Darren J. Wilkinson Abstract This chapter considers stochastic differential equations for Systems Biology models derived
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 informationComputational complexity analysis for Monte Carlo approximations of classically scaled population processes
Computational complexity analysis for Monte Carlo approximations of classically scaled population processes David F. Anderson, Desmond J. Higham, Yu Sun arxiv:52.588v3 [math.a] 4 Jun 28 June 5, 28 Abstract
More informationStochastic Processes at Single-molecule and Single-cell levels
Stochastic Processes at Single-molecule and Single-cell levels Hao Ge haoge@pu.edu.cn Beijing International Center for Mathematical Research 2 Biodynamic Optical Imaging Center Peing University, China
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
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 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 informationExtending the multi-level method for the simulation of stochastic biological systems
Extending the multi-level method for the simulation of stochastic biological systems Christopher Lester Ruth E. Baker Michael B. Giles Christian A. Yates 29 February 216 Abstract The multi-level method
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 informationSystems Driven by Alpha-Stable Noises
Engineering Mechanics:A Force for the 21 st Century Proceedings of the 12 th Engineering Mechanics Conference La Jolla, California, May 17-20, 1998 H. Murakami and J. E. Luco (Editors) @ASCE, Reston, VA,
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationMalliavin Calculus in Finance
Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x
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 informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationNumerical Integration of SDEs: A Short Tutorial
Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx
More information