Supporting Information for The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise
|
|
- Ezra Hutchinson
- 5 years ago
- Views:
Transcription
1 TWCCC Texas Wisconsin California Control Consortium Technical report number Supporting Information for The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise Rishi Srivastava Dept. Chem. & Biol. Eng. Univ. of Wisconsin. (USA) Ethan Mastny Dept. Chem. & Biol. Eng. Univ. of Wisconsin. (USA) Eric L. Haseltine Div. Chemi. & Chem. Eng Cali. Inst. Tech., Pasadena, CA March 31, 2011 James B. Rawlings Dept. Chem. & Biol. Eng. Univ. of Wisconsin. (USA) We report in the statement of some results and the proofs omitted from the paper The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise [3]. All numberings are referred to the paper [3]. S1 Example 1: Derivation for Pap operon regulation Dividing both sides of equations (1) (4) of the paper by r 1 + r 3 and defining ɛ = 1/(r 1 + r 3 ) gives ɛ dp 1 ɛ dp 2 ɛ dp 3 ɛ dp 4 = P 1 + ɛr 2 P 2 + ɛr 4 P 3 (1) = ɛ(r 2 + r 5 )P 2 + K 1 P 1 + ɛr 6 P 4 (2) = ɛ(r 4 + r 7 )P 3 + K 3 P 1 + ɛr 8 P 4 (3) = ɛ(r 6 + r 8 )P 4 + ɛr 5 P 2 + ɛr 7 P 3 (4) (5) in which K 1 = r 1 /(r 1 + r 3 ) and K 3 = r 3 /(r 1 + r 3 ) are O(1) terms. rsrivastava2@wisc.edu. Author to whom all correspondence should be addressed. Eric L. Haseltine eamastny@wisc.edu rawlings@engr.wisc.edu 1
2 TWCCC Technical Report S1.1 Fast time scale solution In the fast time scale we rescale time t to τ = t/ɛ. Using τ we rewrite (1) (4) on the fast time scale as d ˆP 1 = ˆP 1 + ɛr 2 ˆP2 + ɛr 4 ˆP3 (6) d ˆP 2 = ɛ(r 2 + r 5 ) ˆP 2 + K 1 ˆP1 + ɛr 6 ˆP4 (7) d ˆP 3 = ɛ(r 4 + r 7 ) ˆP 3 + K 3 ˆP1 + ɛr 8 ˆP4 (8) d ˆP 4 = ɛ(r 6 + r 8 ) ˆP 4 + ɛr 5 ˆP2 + ɛr 7 ˆP3 (9) in which we denote probability on the fast time scale by ˆP i. We expand ˆP i in a power series Substituting (11) into (6) and retaining terms up to order ɛ 1 give (10) ˆP i = Ŵi0 + ɛŵi1 + ɛ 2 Ŵ i2 + O(ɛ 3 ) (11) dŵ10 + ɛ dŵ11 + = (Ŵ10 + ɛŵ11 + ) + ɛr 2 (Ŵ20 + ɛŵ21 + )+ ɛr 4 (Ŵ30 + ɛŵ31 + ) (12) Comparing O(ɛ 0 ) terms in (12) gives dŵ10 with initial condition Ŵ10(τ = 0) = 1. The solution of (13) is = Ŵ10 (13) Ŵ 10 = e τ (14) which is the equation (5) of the paper. Rewriting equation (7) using equation (11) gives dŵ20 + ɛ dŵ20 + = ɛ(r 2 + r 5 )(Ŵ20 + ɛŵ21 + ) + K 1 (Ŵ10 + ɛŵ11 + )+ ɛr 6 (Ŵ40 + ɛŵ41 + ) (15) Comparing O(ɛ 0 ) terms from (15) gives dŵ20 = K 1 Ŵ 10 (16) with initial condition Ŵ20(τ = 0) = 0. We substitute Ŵ10 from (14) into (16) and solve the differential equation to obtain Ŵ 20 = K 1 (1 e τ ) (17) which is the equation (6) of the paper. We obtain the expansions of the probabilities of states g 3 and g 4, i.e., Ŵ 30 and Ŵ40, analogously.
3 TWCCC Technical Report S1.2 Slow time scale solution In (1) we substitute expansion of P i, 1 i 4 P i = W i0 + ɛw i1 + ɛ 2 W i2 + O(ɛ 3 ) (18) giving ɛ ( dw 10 + ɛ dw 11 + ) = ( W 10 + ɛw 11 + ) + ɛr 2 (W 20 + ɛw 21 + )+ ɛr 4 (W 30 + ɛw 31 ) (19) Collecting O(ɛ 0 ) terms from (19) gives W 10 = 0 (20) which is equation (9) of the paper. Next we collect O(ɛ 1 ) terms from (19) to obtain W 11 = r 2 W 20 + r 4 W 30 (21) A better approximation for the probability of state g 1, P 1, is obtained by approximating P 1 up to O(ɛ). P 1 = W 10 + ɛw 11 (22) Using equations (20) and (21) in equation (22) gives which is equation (13) of the paper. Next we substitute (18) into (2) to obtain P 1 = ɛ(r 2 W 20 + r 4 W 30 ) (23) ɛ ( dw 20 + ɛ dw 21 + ) = ɛ(r 2 + r 5 )(W 20 + ɛw 21 + ) + K 1 (W 10 + ɛw 11 + )+ ɛr 6 (W 40 + ɛw 41 + ) (24) Comparison of O(ɛ 0 ) terms from (24) leads to equation (20). Thus comparison of O(ɛ 0 ) terms from (24) does not lead to any new information. Comparing O(ɛ 1 ) terms from (24) we obtain: dw 20 = (r 2 + r 5 )W 20 + r 6 W 40 + K 1 W 11 (25) Substituting W 11 from (21) into (25) gives dw 20 = [ r 1 + r 5 ]W 20 + r 2 W 30 + r 6 W 40 (26) which is the same as equation (10) of the paper. Evolution equations for W 30 and W 40 can be obtained similarly.
4 TWCCC Technical Report S2 Example 2: Derivations for robustness to noise in a biochemical oscillator. S2.1 sqspa reduction We consider reaction network (15) of the paper These reactions include the transcription events to form M A, translation to form A, and degradation of M A. These reactions do not lead to change in the populations of species D A and D A, hence we write P (D A, D A, M A, A) as P (M A, A).The master equation for this system is dp (M A, A) = α A D A P (M A 1, A) + α AD AP (M A 1, A) + δ MA (M A + 1)P (M A + 1, A)+ β A M A P (M A, A 1) [ α A D A + α AD A + (δ MA + β A )M A ] P (MA, A) (27) When the condition β A, δ MA α A, α A holds, the mrna species M A is a QSSA species. We further define 1 ɛ = δ MA = ɛδ MA βa = ɛβ A (β A + δ MA ) Multiplying (27) with ɛ gives us ɛ dp (M A, A) = ɛ(α A D A + α AD A)P (M A 1, A) + δ MA (M A + 1)P (M A + 1, A)+ β A M A P (M A, A 1) [ ɛ(α A D A + α AD A) + ( δ MA + β A )M A ] P (MA, A) (28) We expand P (M A, A) in a power series P (M A, A) = W 0 (M A, A) + ɛw 1 (M A, A) + ɛ 2 W 2 (M A, A) + (29) Then by grouping like powers of ɛ in (28), we obtain ɛ 0 : 0 = δ MA (M A + 1) W 0 (M A +1, A)+ β A M A W 0 (M A, A 1) ( δma + β ) A MA W 0 (M A, A) (30) Expressing this result for M A = 0, 1,... yields W 0 (M A, A) = 0 for M A > 0 for all A (31) Similarly, equating terms of O(ɛ 1 ) yields ɛ 1 dw 0 (0, A), M A = 0 : = MA W 1 (1, A) ( α A D A + α AD A) W0 (0, A) (32) ɛ 1, M A = 1 : W 1 (1, A) = β A W 1 (1, A 1) + (α A D A + α AD A)W 0 (0, A) (33) Solving the last equation for W 1 (1, A) gives A W 1 (1, A) = (α A D A + α AD A)W 0 (0, A) + β j A (α AD A + α AD A)W 0 (0, A j) Substituting this result into (32) gives j=1 dw 0 (0, A) = δ A MA β j A (α AD A + α AD A)W 0 (0, A j) β A (α A D A + α AD A)W 0 (0, A) (34) j=1 The catalytic reactions in the paper, reaction sets (16) and (17), follow from this equation.
5 TWCCC Technical Report S2.2 Deterministic quasi steady state classical (dqc) reduction Again we consider the reaction network (15) of the paper. The master equation describing the species M A can be written as: dm A = α A D A + α AD A δ MA M A β A M A (35) Under the condition β A, δ MA α A, α A, δ A, M A is a QSSA species. To apply dqc we set left hand side of equation (35) to 0 to obtain following expression for the population of M A. M A = α AD A + α A D A δ MA + β A (36) The rate of formation of A is equal to the rate of last reaction of the network da =β AM A = β A (α A D A + α AD A) (37) Since M A is already removed from the network. The only species that changes with time is species A. Equations (36) and (37) give rise to the dqc reduced reactions (18) and (19) of the paper. S3 Example 3: Derivation for Fast fluctuation We apply Ω expansion as outlined by van Kampen [4] on the master equation of reactions (20) (22) of the paper. The master equation for these reactions is dp (a, c, g) = k 1 Ω (a + 1)gP (a + 1, c 1, g) + k 2 (c + 1)(g 1)P (a 1, c + 1, g 1) Ω + k ( 3 k1 (g + 1)gP (a, c, g + 1) Ω Ω ag + k 2 Ω cg + k ) 3 g(g 1) P (a, c, g) (38) Ω There are two phases of time evolution of the species population. In the first phase, species G increases and species A and C start fluctuating. After species G has risen above a certain threshold value, we switch the simulation methodology to the SSA-Ω approach. We write g as a continuous random variable: g = Ωφ G + Ω 1/2 ξ (39) where φ G is the deterministic mean concentration of species G in the regime of SSA-Ω and Ω 1/2 ξ is the noise in G. For the transformed variable g, the probability P (a, c, g) corresponds to Π(a, c, ξ). In (38) we can write dp (a, c, g) Π(a, c, ξ) = Under this change of variables, (38) becomes Π(a, c, ξ) Ω 1/2 dφ g Π(a, c, ξ) = k 1 Ω (a+1)(ωφ G + Ω 1/2 ξ)π(a+1, c 1, ξ) Ω 1/2 dφ g Π(a, c, ξ) + k 2 Ω (c+1)[1 Ω 1/2 2 2 ](Ωφ G+Ω 1/2 ξ)π(a 1, c+1, ξ) + k 3 Ω [1 Ω 1/2 2 2 ](Ωφ G+Ω 1/2 ξ)(ωφ G +Ω 1/2 ξ 1)Π(a, c, ξ) ( k1 Ω a(ωφ G+Ω 1/2 ξ) + k 2 Ω c(ωφ G+Ω 1/2 ξ)+ k ) 3 Ω (Ωφ G+Ω 1/2 ξ)(ωφ G +Ω 1/2 ξ 1) Π(a, c, ξ) (40)
6 TWCCC Technical Report For the rate constants of the system, we have k 1, k 2 k 3. We define the following parameters K 1 = k 1 Ω = Ωk 1 2 (41) Following Mastny et al. [2], we express Π as a power series in Ω 1, Π(a, c, ξ) = W 0 (a, c, ξ) + Ω 1 W 1 (a, c, ξ) + Ω 2 W 2 (a, c, ξ) + and substitute into (40) to obtain the following expanded master equation ( Ω 1 W0 (a, c, ξ) +Ω 1 W ) 1(a, c, ξ) + Ω 1/2 dφ ( G W0 (a, c, ξ) +Ω 1 W ) 1(a, c, ξ) + = K 1 Ω 1 (a+1)(ωφ G + Ω 1/2 ξ) ( W 0 (a+1, c 1, ξ)+ω 1 W 1 (a + 1, c 1, ξ)+ ) + (Ω) 1 (c+1)[1 Ω 1/2 2 2 ](Ωφ G+Ω 1/2 ξ) ( W 0 (a 1, c+1, ξ)+ω 1 W 1 (a 1, c+1, ξ)+ ) + k 3 Ω 2 [1 Ω 1/2 2 2 ](Ωφ G+Ω 1/2 ξ)(ωφ G +Ω 1/2 ξ 1) ( W 0 (a, c, ξ) + Ω 1 W 1 (a, c, ξ) + ) ( ) (W0 K 1 Ω 1 a(ωφ G +Ω 1/2 ξ) + 1 Ω 1 c(ωφ G +Ω 1/2 ξ) (a, c, ξ) + Ω 1 W 1 (a, c, ξ) + ) k 3 Ω 2 (Ωφ G +Ω 1/2 ξ)(ωφ G +Ω 1/2 ξ 1) ( W 0 (a, c, ξ) + Ω 1 W 1 (a, c, ξ) + ) (42) Ω 0 terms. Collecting Ω 0 terms from (42) gives 0 = K 1 (a + 1)φ G W 0 (a + 1, c 1, ξ) ( K 1 aφ G + 1 cφ G ) W0 (a, c, ξ) + 1 (c + 1)φ G W 0 (a 1, c + 1, ξ) (43) Since none of the terms multiplying joint densities in (43) contain ξ, we can integrate (43) to obtain an equation for marginal density as a function of a and c only. Also we note that because of the stoichiometry of the problem a + c = N 0 and we can replace c with N 0 a. Making these substitutions in (43), we obtain 0 = K 1 (a + 1)W 0a (a + 1) K 1 aw 0a (a) 1 N 0 W 0a (a) + 1 aw 0a (a) + 1 N 0 W 0a (a 1) 1 (a 1)W 0a (a 1) (44) An analytical expression for marginal density W 0a (a) can be obtained by taking the z-transform of the equation (44). That produces equation (23) of the paper. The distribution in (23) is a binomial distribution W 0a Bin( q 1+q, N 0) in which q = 1/(K 1 ). Since a + c = N 0 equation (24) is an immediate consequence of (23). From equation (43), it is also clear that W 0 (a, ξ) = W 0a (a)w 0ξ (ξ). Ω 1/2 terms. Collecting Ω 1/2 terms from equation (42) gives dφ G W 0 (a, c, ξ) =K 1 (a+1)ξw 0 (a+1, c 1, ξ) K 1 aξw 0 (a, c, ξ)+ 1 (c+1)ξw 0 (a 1, c+1, ξ) 1 cξw 0 (a, c, ξ) + k 3 ( φ 2 G W 0(a, c, ξ) ) (45)
7 TWCCC Technical Report We recall that argument c in function W 0 is redundant because a + c = N 0 at all times. Summing equation (45) over all the values of a and c along with using the independence of a and ξ gives us dφ G = 1 c φ G + k 3 φ 2 G (46) which is (25) in the paper. Equations (23) and (46) together give an approximation of the populations of A, C and G when G has risen over a threshold value. Ω 1 terms. Collecting Ω 1 terms from equation (42) gives W 0 (a, c, ξ) =K 1 (a + 1)φ G W 1 (a + 1, c 1, ξ) K 1 aφ G W 1 (a 1, c + 1, ξ) + c + 1 (ξw 0(a 1, c + 1, ξ)) + c (φ GW 0 (a 1, c + 1, ξ) cφ G W 1 (a, c, ξ)) + 2k 3 φ G (ξw 0(a, c, ξ)) + k φ2 G 2 W 0(a, c, ξ) Summing equation (47) for all values of a and c gives 2 (47) W 0ξ (ξ) = c (ξw 0ξ(ξ))+ 1 2 c φ G 2 2 W 0ξ(ξ)+2k 3 φ G (ξw 0ξ(ξ))+ k φ2 G Rearranging equation (48) gives W 0ξ (ξ) = [( ) ] c 2k 3φ G ξw 0ξ (ξ) [( ) ] c φg k 3 φ 2 G W 0ξ (ξ) 2 W 0ξ(ξ) (48) Equation (49) is a linear Fokker-Planck equation and is the same as equation (26) of the paper. Its equivalent stochastic differential equation is given by ( ) c dξ = 2k c φ G 3φ G ξ + + k 3 φ 2 G dw (50) in which dw is a normally distributed random variable with zero mean and variance [1]. Equation (50) here is the same as equation (27) of the paper. References [1] Peter E. Kloeden and Eckhard Platen. Numerical solution of stochastic differential equations. Springer Verlag, Berlin, [2] Ethan A. Mastny, Eric L. Haseltine, and James B. Rawlings. Two classes of quasi-steady-state model reductions for stochastic kinetics. J. Chem. Phys., 127(9):094106, September [3] Rishi Srivastava, Eric L. Haseltine, Ethan A. Mastny, and James B. Rawlings. The stochastic quasi-steady-state assumption: Reducing the model but not the noise. J. Chem. Phys., Accepted for publication. [4] N. G. van Kampen. Stochastic Processes in Physics and Chemistry. Elsevier Science Publishers, Amsterdam, The Netherlands, second edition, (49)
Extending 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 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 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 informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
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 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 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 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 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 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 informationArticle (peer-reviewed)
Title Author(s) Influence of noise intensity on the spectrum of an oscillator Swain, Rabi Sankar; Gleeson, James P.; Kennedy, Michael Peter Publication date 2005-11 Original citation Type of publication
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 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 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 informationRobust stability of moving horizon estimation under bounded disturbances
WCCC eas Wisconsin California Control Consortium echnical report number 213-1 Robust stability of moving horizon estimation under bounded disturbances Luo Ji Dept. Chem. & Biol. Eng. Univ. of Wisconsin.
More informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
More informationA variational radial basis function approximation for diffusion processes
A variational radial basis function approximation for diffusion processes Michail D. Vrettas, Dan Cornford and Yuan Shen Aston University - Neural Computing Research Group Aston Triangle, Birmingham B4
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 informationSupporting Information. Methods. Equations for four regimes
Supporting Information A Methods All analytical expressions were obtained starting from quation 3, the tqssa approximation of the cycle, the derivation of which is discussed in Appendix C. The full mass
More informationNonlinear Model Predictive Control Tools (NMPC Tools)
Nonlinear Model Predictive Control Tools (NMPC Tools) Rishi Amrit, James B. Rawlings April 5, 2008 1 Formulation We consider a control system composed of three parts([2]). Estimator Target calculator Regulator
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 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 informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationthe noisy gene Biology of the Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB)
Biology of the the noisy gene Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB) day III: noisy bacteria - Regulation of noise (B. subtilis) - Intrinsic/Extrinsic
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationGENERATION OF COLORED NOISE
International Journal of Modern Physics C, Vol. 12, No. 6 (2001) 851 855 c World Scientific Publishing Company GENERATION OF COLORED NOISE LORENZ BARTOSCH Institut für Theoretische Physik, Johann Wolfgang
More informationOptimizing Economic Performance using Model Predictive Control
Optimizing Economic Performance using Model Predictive Control James B. Rawlings Department of Chemical and Biological Engineering Second Workshop on Computational Issues in Nonlinear Control Monterey,
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
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 informationAdiabatic piston - Draft for MEMPHYS
diabatic piston - Draft for MEMPHYS 4th ovember 2006 Christian GRUBER (a) and nnick LESE (b) (a) Institut de héorie des Phénomènes Physiques École Polytechnique Fédérale de Lausanne, CH-05 Lausanne, Switzerland
More informationSupporting Information: A quantitative comparison of srna-based and protein-based gene regulation
Supporting Information: A quantitative comparison of srna-based and protein-based gene regulation Pankaj Mehta a,b, Sidhartha Goyal b, and Ned S. Wingreen a a Department of Molecular Biology, Princeton
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 informationApproximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach
Approximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach Henrik T. Sykora, Walter V. Wedig, Daniel Bachrathy and Gabor Stepan Department of Applied Mechanics,
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 informationQuântica Oscilador Paramétrico
Luz e Átomos como ferramentas para Informação Quântica Oscilador Paramétrico Ótico Inst. de Física Marcelo Martinelli Lab. de Manipulação Coerente de Átomos e Luz Parametric Down Conversion Energy and
More informationAveraging II: Adiabatic Invariance for Integrable Systems (argued via the Averaging Principle)
Averaging II: Adiabatic Invariance for Integrable Systems (argued via the Averaging Principle In classical mechanics an adiabatic invariant is defined as follows[1]. Consider the Hamiltonian system with
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 informationThermodynamics of computation.
Thermodynamics of computation. Dominique Chu School of Computing University of Kent, UK D.F.Chu@kent.ac.uk C 3 Symposium December 11, 2017 Outline 1 Computation and Life itself 2 Living computers 3 Energy
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationA Generic Multivariate Distribution for Counting Data
arxiv:1103.4866v1 [stat.ap] 24 Mar 2011 A Generic Multivariate Distribution for Counting Data Marcos Capistrán and J. Andrés Christen Centro de Investigación en Matemáticas, A. C. (CIMAT) Guanajuato, MEXICO.
More information1 Introduction. Maria Luisa Guerriero. Jane Hillston
Time-Series Analysis of Circadian Clocks: Quantifying the Effect of Stochasticity on the Distribution of Phase, Period and Amplitude of Circadian Oscillations Maria Luisa Guerriero Centre for Systems Biology
More informationStochastic Processes around Central Dogma
Stochastic Processes around Central Dogma Hao Ge haoge@pu.edu.cn Beijing International Center for Mathematical Research Biodynamic Optical Imaging Center Peing University, China http://www.bicmr.org/personal/gehao/
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 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 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 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 informationMelting of Ultrathin Lubricant Film Due to Dissipative Heating of Friction Surfaces
ISSN 6-78, Technical Physics, 7, Vol. 5, No. 9, pp. 9. Pleiades Publishing, Ltd., 7. Original Russian Text.V. Khomenko, I.. Lyashenko, 7, published in Zhurnal Tekhnicheskoœ Fiziki, 7, Vol. 77, No. 9, pp.
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 informationNonlinear Stochastic Modeling and State Estimation of Weakly Observable Systems: Application to Industrial Polymerization Processes
Nonlinear Stochastic Modeling and State Estimation of Weakly Observable Systems: Application to Industrial Polymerization Processes Fernando V. Lima, James B. Rawlings and Tyler A. Soderstrom Department
More informationDETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION
.4 DETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION Peter C. Chu, Leonid M. Ivanov, and C.W. Fan Department of Oceanography Naval Postgraduate School Monterey, California.
More informationarxiv: v1 [cond-mat.stat-mech] 11 Jan 2012
arxiv:121.2361v1 [cond-mat.stat-mech] 11 Jan 212 Stochastic Description of a Bistable Frustrated Unit Ashok Garai 1, Bartlomiej Waclaw 2, Hannes Nagel 3, and Hildegard Meyer-Ortmanns 1 1 School of Engineering
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 informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
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 informationarxiv: v7 [quant-ph] 22 Aug 2017
Quantum Mechanics with a non-zero quantum correlation time Jean-Philippe Bouchaud 1 1 Capital Fund Management, rue de l Université, 75007 Paris, France. (Dated: October 8, 018) arxiv:170.00771v7 [quant-ph]
More informationPeter E. Kloeden Eckhard Platen. Numerical Solution of Stochastic Differential Equations
Peter E. Kloeden Eckhard Platen Numerical Solution of Stochastic Differential Equations Peter E. Kloeden School of Computing and Mathematics, Deakin Universit y Geelong 3217, Victoria, Australia Eckhard
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 informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part I: Linear Systems with Gaussian Noise James B. Rawlings and Fernando V. Lima Department of Chemical and Biological Engineering University of
More informationMultivariate Bayesian Linear Regression MLAI Lecture 11
Multivariate Bayesian Linear Regression MLAI Lecture 11 Neil D. Lawrence Department of Computer Science Sheffield University 21st October 2012 Outline Univariate Bayesian Linear Regression Multivariate
More informationLecture 2: Open quantum systems
Phys 769 Selected Topics in Condensed Matter Physics Summer 21 Lecture 2: Open quantum systems Lecturer: Anthony J. Leggett TA: Bill Coish 1. No (micro- or macro-) system is ever truly isolated U = S +
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 informationA scaling limit from Euler to Navier-Stokes equations with random perturbation
A scaling limit from Euler to Navier-Stokes equations with random perturbation Franco Flandoli, Scuola Normale Superiore of Pisa Newton Institute, October 208 Newton Institute, October 208 / Subject of
More informationAntibound State for Klein-Gordon Equation
International Journal of Mathematical Analysis Vol. 8, 2014, no. 59, 2945-2949 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.411374 Antibound State for Klein-Gordon Equation Ana-Magnolia
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 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 informationIntrinsic Noise in Nonlinear Gene Regulation Inference
Intrinsic Noise in Nonlinear Gene Regulation Inference Chao Du Department of Statistics, University of Virginia Joint Work with Wing H. Wong, Department of Statistics, Stanford University Transcription
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationURL: <
Citation: ngelova, Maia and en Halim, sma () Dynamic model of gene regulation for the lac operon. Journal of Physics: Conference Series, 86 (). ISSN 7-696 Published by: IOP Publishing URL: http://dx.doi.org/.88/7-696/86//7
More informationStochastic Behavior of Dissipative Hamiltonian Systems with Limit Cycles
Stochastic Behavior of Dissipative Hamiltonian Systems with Limit Cycles Wolfgang Mathis Florian Richter Richard Mathis Institut für Theoretische Elektrotechnik, Gottfried Wilhelm Leibniz Universität Hannover,
More informationDynamics and stochastic limit for phase transition problems. problems with noise. Dimitra Antonopoulou IACM-FORTH
Dynamics and stochastic limit for phase transition problems with noise IACM-FORTH International Conference on Applied Mathematics, September 16-20, 2013, ACMAC, Heraklion, Greece Table of contents 1 1.1
More informationQuantifying Intermittent Transport in Cell Cytoplasm
Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent
More informationOpen Quantum Systems and Markov Processes II
Open Quantum Systems and Markov Processes II Theory of Quantum Optics (QIC 895) Sascha Agne sascha.agne@uwaterloo.ca July 20, 2015 Outline 1 1. Introduction to open quantum systems and master equations
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationA Bayesian method for the analysis of deterministic and stochastic time series
A Bayesian method for the analysis of deterministic and stochastic time series Coryn Bailer-Jones Max Planck Institute for Astronomy, Heidelberg DPG, Berlin, March 2015 Time series modelling heteroscedastic,
More informationModelling in Biology
Modelling in Biology Dr Guy-Bart Stan Department of Bioengineering 17th October 2017 Dr Guy-Bart Stan (Dept. of Bioeng.) Modelling in Biology 17th October 2017 1 / 77 1 Introduction 2 Linear models of
More informationFE670 Algorithmic Trading Strategies. Stevens Institute of Technology
FE670 Algorithmic Trading Strategies Lecture 8. Robust Portfolio Optimization Steve Yang Stevens Institute of Technology 10/17/2013 Outline 1 Robust Mean-Variance Formulations 2 Uncertain in Expected Return
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 informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part II: Observability and Stability James B. Rawlings and Fernando V. Lima Department of Chemical and Biological Engineering University of Wisconsin
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 informationCorrespondence of regular and generalized mass action systems
Correspondence of regular and generalized mass action systems Van Vleck Visiting Assistant Professor University of Wisconsin-Madison Joint Mathematics Meetings (San Antonio, TX) Saturday, January 10, 2014
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 informationthe renormalization group (RG) idea
the renormalization group (RG) idea Block Spin Partition function Z =Tr s e H. block spin transformation (majority rule) T (s, if s i ; s,...,s 9 )= s i > 0; 0, otherwise. b Block Spin (block-)transformed
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 informationLinearizing the Hyperinflation Model
Linearizing the Hyperinflation Model Klaus Adam George Evans Seppo Honkapohja March 7, 003 Abstract This is a technical appendix to Are Stationary Hyperinflation Paths Learnable? by the same authors which
More informationIntroduction to Stochastic SIR Model
Introduction to Stochastic R Model Chiu- Yu Yang (Alex), Yi Yang R model is used to model the infection of diseases. It is short for Susceptible- Infected- Recovered. It is important to address that R
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 informationFirst Excursion Probabilities of Non-Linear Dynamical Systems by Importance Sampling. REN Limei [a],*
Progress in Applied Mathematics Vol. 5, No. 1, 2013, pp. [41 48] DOI: 10.3968/j.pam.1925252820130501.718 ISSN 1925-251X [Print] ISSN 1925-2528 [Online] www.cscanada.net www.cscanada.org First Excursion
More informationA Unified Formulation of Gaussian Versus Sparse Stochastic Processes
A Unified Formulation of Gaussian Versus Sparse Stochastic Processes Michael Unser, Pouya Tafti and Qiyu Sun EPFL, UCF Appears in IEEE Trans. on Information Theory, 2014 Presented by Liming Wang M. Unser,
More informationDynamics of a mass-spring-pendulum system with vastly different frequencies
Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate
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 informationDYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS
Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California USA DETC2005-84017
More informationAnalysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium
Analysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium Adriano W. da Silva May 30, 2013 Abstract A binary gaseous mixture with reversible reaction of type A + A = B
More information5 December 2016 MAA136 Researcher presentation. Anatoliy Malyarenko. Topics for Bachelor and Master Theses. Anatoliy Malyarenko
5 December 216 MAA136 Researcher presentation 1 schemes The main problem of financial engineering: calculate E [f(x t (x))], where {X t (x): t T } is the solution of the system of X t (x) = x + Ṽ (X s
More informationSuperparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows
Superparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows June 2013 Outline 1 Filtering Filtering: obtaining the best statistical estimation of a nature
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 informationImplicit sampling for particle filters. Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins
0/20 Implicit sampling for particle filters Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins University of California at Berkeley 2/20 Example: Try to find people in a boat in the middle of
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 informationOnline monitoring of MPC disturbance models using closed-loop data
Online monitoring of MPC disturbance models using closed-loop data Brian J. Odelson and James B. Rawlings Department of Chemical Engineering University of Wisconsin-Madison Online Optimization Based Identification
More informationAPPLICATION OF A STOCHASTIC REPRESENTATION IN NUMERICAL STUDIES OF THE RELAXATION FROM A METASTABLE STATE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 7 (1), 83-90 (2001) APPLICATION OF A STOCHASTIC REPRESENTATION IN NUMERICAL STUDIES OF THE RELAXATION FROM A METASTABLE STATE FERDINANDO DE PASQUALE 1, ANTONIO
More information