Stochastic thermodynamics of information

Size: px
Start display at page:

Download "Stochastic thermodynamics of information"

Transcription

1 Stochastic thermodynamics of information

2

3 Introduction Textbook, Review Sekimoto, K. (2010). Stochastic energetics (Vol. 799). Springer. Seifert, U. (2012). Reports on Progress in Physics, 75(12),

4 Introduction Textbook, Review Sekimoto, K. (2010). Stochastic energetics (Vol. 799). Springer. Seifert, U. (2012). Reports on Progress in Physics, 75(12),

5 Introduction Review Parrondo, J. M., Horowitz, J. M., & Sagawa, T. (2015). Thermodynamics of information. Nature physics, 11(2), 131. My Ph, D. thesis Ito, S. (2016). Information thermodynamics on causal networks and its application to biochemical signal transduction. Springer. Toyabe, S., Sagawa, T., Ueda, M., Muneyuki, E., & Sano, M. (2010). Nature physics, 6(12), 988.

6 Introduction Ito, S., & Sagawa, T. (2013). Physical review letters, 111(18), Shiraishi, N., Ito, S., Kawaguchi, K., & Sagawa, T. (2015). New Journal of Physics, 17(4), Ito, S. (2016). Scientific reports, 6. Ito, S., & Sagawa, T. (2015). Nature communications, 6. Yamamoto, S., Ito, S., Shiraishi, N., & Sagawa, T. (2016). Physical Review E, 94(5), Ito, S. (2017). arxiv preprint arxiv:

7 Introduction

8

9

10 dq = du dw dw U d dq du dw λ. dq ds T ds dq T

11 dq ds T ds bath = dq T ds + ds bath 0

12

13 m X ẍ(t) = X x U X (x(t), X (t)) + p 2 X T X X (t) ξx x(t) t mx γx UX λx TX mx γx γx=1 ẋ(t) x U X (x(t), X (t)) + p 2T X X (t)

14 dux. du X (x(t), X (t)) = x U X (x(t), X (t))dt + X X U X (x(t), X (t))dt dt dwx dqx dw X (x(t), X (t)) = X X U X (x(t), X (t))dt dq X (x(t), X (t)) = x U X (x(t), X (t))dt dq X = du X dw X

15 dq X (x(t), X (t)) = x U X (x(t), X (t))dt ẋ(t) x U X (x(t), X (t)) + p 2T X X (t) dq X =ẋ(t) ( p 2T X X (t) ẋ(t))dt jx=dqx/dt j X (t) =ẋ(t) ( p 2T X X (t) ẋ(t))

16

17 X. x. px(x) X Z H(X) = dxp X (x)lnp X (x) =h ln p X (x)i px1,,xn (x1,, xn) H(X 1,...,X n )= Z dx 1 dx n p X1,...,X n (x 1,...,x n )lnp X1,...,X n (x 1,...,x n ) = h ln p X1,...,X n (x 1,...,x n )i

18 Y X H(X Y )=H(X, Y ) H(Y ) X, Y I(X; Y )=H(X) H(X Y ) Z X Y I(X; Y Z) =H(X Z) H(X Y,Z)

19 I(X; Y )=I(Y; X) I(X; Y ) 0 I(X; Y )=0 p X,Y (x, y) =p X (x)p Y (y) x, y X Y 0 I(X; Y ) H(X Y ) H(Y X) H(X) H(Y ) I(X; Y )=H(X)+H(Y) H(X, Y ) = H(X) H(X Y ) = H(Y ) H(Y X)

20 px(x), qx(x) D(p X q X )= Z dxp X (x)ln p X(x) q X (x) I(X; Y )=D(p X,Y p X p Y )

21 D(p X q X ) 0 D(p X q X )=0 p X (x) =q X (x) x F, x G(x), px(x) Z Z F (G(x))p X (x)dx apple F G(x)p X (x)dx hf (G(x))i applef (hg(x)i) F=ln G(x)= qx(x)/px(x), D(p X q X )=hln[q X (x)/p X (x)]i appleln[hq X (x)/p X (x)i] =ln1=0

22

23 ẋ(t) x U X (x(t), X (t)) + p 2T X X (t) X (t)dt = db X t x t = x(t) X t = X (t) x t+dt x t x U X (x t, X t)dt + p 2T X db X t dbxt p(db Xt )= 1 p 2 dt exp apple (dbxt ) 2 X t+dt = 1 p 2TX xt xt+dt T (x t+dt ; x t )=p Xt+dt X t = 1 p 4 TX dt exp apple (xt+dt x t x U X (x t, X t )dt) 2 4T X dt

24 xt+dt xt. p B X t X t+dt (x t x t+dt )=T (x t ; x t+dt ) apple p Xt+dt X t (x t+dt x t ) p B X t X t+dt (x t x t+dt ) =exp jx (t) T X dt T (x t+dt ; x t )=p Xt+dt X t = apple 1 p 4 TX dt exp (xt+dt x t x U X (x t, X t )dt) 2 4T X dt j X (t) =ẋ(t) ( p 2T X X (t) ẋ(t))

25 D(p Xt+dt X t p Xt p B X t X t+dt p Xt+dt ) 0 D(p Xt+dt X t p Xt p B X t X t+dt p Xt+dt )= hj X (t)idt T X + H(X t+dt ) H(X t ) X ds X (t) =H(X t+dt ) H(X t ) ds bath (t) = hj X (t)i T X dt ds X (t)+ds bath (t) 0

26

27 λx m. ẋ(t) x U X (x(t), X (m, t)) + p 2T X X (t) m T m (x t ; x t+dt )=p Xt+dt X t,m (x t+dt x t,m)

28 m p B X t X t+dt,m (x t x t+dt,m)=t m (x t ; x t+dt ) p Xt+dt X t,m (x t+dt x t,m) p B X t X t+dt,m (x t x t+dt,m) =exp apple jx (t) T X dt j X (t) =ẋ(t) ( p 2T X X (t) ẋ(t))

29 D(p Xt+dt X t,m p Xt,M p B X t X t+dt,m p Xt+dt,M ) 0 D(p Xt+dt X t,m p Xt,M p B X t X t+dt,m p Xt+dt,M )= X hj X (t)i T X dt + H(X t+dt,m) H(X t,m) ds X (t) =H(X t+dt ) H(X t ) ds bath (t) = ds X (t)+ds bath (t) I(X t+dt ; M) I(X t ; M) hj X (t)i T X dt

30

31 Stochastic thermodynamics for 2D Langevin equations ẋ(t) x U X (x(t),y(t)) + p 2T X X (t) ẏ(t) y U Y (x(t),y(t)) + p 2T Y Y (t) ξx ξy p Xt+dt X t,y t (x t+dt x t,y t )=T X y t (x t+dt ; x t ) Ty X t (x t+dt ; x t ) Ty X t (x t ; x t+dt ) =exp apple p Yt+dt X t,y t (y t+dt x t,y t )=T Y x t (y t+dt ; y t ) jx (t)dt T X Tx Y apple t (y t+dt ; y t ) Tx Y t (y t ; y t+dt ) =exp jy (t)dt T Y j X (t) =ẋ(t) ( p 2T X X (t) ẋ(t)) j Y (t) =ẏ(t) ( p 2T Y Y (t) ẏ(t))

32 Stochastic thermodynamics for 2D Langevin equations D(p Xt+dt X t,y t p Yt+dt X t,y t p Xt,Y t p B X t X t+dt,y t+dt p B Y t X t+dt,y t+dt p Xt+dt,Y t+dt ) 0 D(p Xt+dt X t,y t p Yt+dt X t,y t p Xt,Y t p B X t X t+dt,y t+dt p B Y t X t+dt,y t+dt p Xt+dt,Y t+dt ) = hj X (t)idt T X X Y : hj Y (t)idt T Y + H(X t+dt,y t+dt ) H(X t,y t ) ds X,Y + ds bath,x + ds bath,y 0 ds X,Y = H(X t+dt,y t+dt ) H(X t,y t ) ds bath,x = ds bath,y = hj X (t)idt T X hj Y(t)idt T Y

33 Stochastic thermodynamics for 2D Langevin equations D(p Xt+dt X t,y t p Xt,Y t,y t+dt p B X t X t+dt,y t+dt p Yt,X t+dt,y) 0 D(p Xt+dt X t,y t p Xt,Y t,y t+dt p B X t X t+dt,y t+dt p Yt,X t+dt,y t+dt ) = hj X (t)idt T X + H(X t+dt ) H(X t )+I(X t ; {Y t,y t+dt }) I(X t+dt ; {Y t,y t+dt }) X: ds X (t) =H(X t+dt ) H(X t ) ds bath,x = hj X (t)idt T X ds X + ds bath,x I(X t+dt ; {Y t,y t+dt }) I(X t ; {Y t,y t+dt })=di

34 Stochastic thermodynamics for 2D Langevin equations X ds X,Y Y ds X,Y + ds bath,x + ds bath,y 0 hj X idt hj Y idt ds bath,x T X ds bath,y T Y ds X X di Y ds X + ds bath,x di hj X idt hj Y idt ds bath,x T X ds bath,y T Y

35 Summary

Maxwell's Demon in Biochemical Signal Transduction

Maxwell's Demon in Biochemical Signal Transduction Maxwell's Demon in Biochemical Signal Transduction Takahiro Sagawa Department of Applied Physics, University of Tokyo New Frontiers in Non-equilibrium Physics 2015 28 July 2015, YITP, Kyoto Collaborators

More information

Information Thermodynamics on Causal Networks

Information Thermodynamics on Causal Networks 1/39 Information Thermodynamics on Causal Networks FSPIP 2013, July 12 2013. Sosuke Ito Dept. of Phys., the Univ. of Tokyo (In collaboration with T. Sagawa) ariv:1306.2756 The second law of thermodynamics

More information

Emergent Fluctuation Theorem for Pure Quantum States

Emergent Fluctuation Theorem for Pure Quantum States Emergent Fluctuation Theorem for Pure Quantum States Takahiro Sagawa Department of Applied Physics, The University of Tokyo 16 June 2016, YITP, Kyoto YKIS2016: Quantum Matter, Spacetime and Information

More information

Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.

Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation. Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =

More information

arxiv: v2 [cond-mat.stat-mech] 16 Mar 2012

arxiv: v2 [cond-mat.stat-mech] 16 Mar 2012 arxiv:119.658v2 cond-mat.stat-mech] 16 Mar 212 Fluctuation theorems in presence of information gain and feedback Sourabh Lahiri 1, Shubhashis Rana 2 and A. M. Jayannavar 3 Institute of Physics, Bhubaneswar

More information

Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs.

Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs. Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs. D. Bernard in collaboration with M. Bauer and (in part) A. Tilloy. IPAM-UCLA, Jan 2017. 1 Strong Noise Limit of (some) SDEs Stochastic differential

More information

Math 201 Solutions to Assignment 1. 2ydy = x 2 dx. y = C 1 3 x3

Math 201 Solutions to Assignment 1. 2ydy = x 2 dx. y = C 1 3 x3 Math 201 Solutions to Assignment 1 1. Solve the initial value problem: x 2 dx + 2y = 0, y(0) = 2. x 2 dx + 2y = 0, y(0) = 2 2y = x 2 dx y 2 = 1 3 x3 + C y = C 1 3 x3 Notice that y is not defined for some

More information

arxiv: v2 [cond-mat.stat-mech] 3 Jun 2018

arxiv: v2 [cond-mat.stat-mech] 3 Jun 2018 Marginal and Conditional Second Laws of Thermodynamics Gavin E. Crooks 1 and Susanne Still 2 1 Theoretical Institute for Theoretical Science 2 University of Hawai i at Mānoa Department of Information and

More information

LANGEVIN EQUATION AND THERMODYNAMICS

LANGEVIN EQUATION AND THERMODYNAMICS LANGEVIN EQUATION AND THERMODYNAMICS RELATING STOCHASTIC DYNAMICS WITH THERMODYNAMIC LAWS November 10, 2017 1 / 20 MOTIVATION There are at least three levels of description of classical dynamics: thermodynamic,

More information

MATH 317 Fall 2016 Assignment 5

MATH 317 Fall 2016 Assignment 5 MATH 37 Fall 26 Assignment 5 6.3, 6.4. ( 6.3) etermine whether F(x, y) e x sin y îı + e x cos y ĵj is a conservative vector field. If it is, find a function f such that F f. enote F (P, Q). We have Q x

More information

Gaussian processes for inference in stochastic differential equations

Gaussian processes for inference in stochastic differential equations Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017

More information

A = (a + 1) 2 = a 2 + 2a + 1

A = (a + 1) 2 = a 2 + 2a + 1 A = (a + 1) 2 = a 2 + 2a + 1 1 A = ( (a + b) + 1 ) 2 = (a + b) 2 + 2(a + b) + 1 = a 2 + 2ab + b 2 + 2a + 2b + 1 A = ( (a + b) + 1 ) 2 = (a + b) 2 + 2(a + b) + 1 = a 2 + 2ab + b 2 + 2a + 2b + 1 3 A = (

More information

Name: SOLUTIONS Date: 11/9/2017. M20550 Calculus III Tutorial Worksheet 8

Name: SOLUTIONS Date: 11/9/2017. M20550 Calculus III Tutorial Worksheet 8 Name: SOLUTIONS Date: /9/7 M55 alculus III Tutorial Worksheet 8. ompute R da where R is the region bounded by x + xy + y 8 using the change of variables given by x u + v and y v. Solution: We know R is

More information

MATH4210 Financial Mathematics ( ) Tutorial 7

MATH4210 Financial Mathematics ( ) Tutorial 7 MATH40 Financial Mathematics (05-06) Tutorial 7 Review of some basic Probability: The triple (Ω, F, P) is called a probability space, where Ω denotes the sample space and F is the set of event (σ algebra

More information

Lecture 5: Integrals and Applications

Lecture 5: Integrals and Applications Lecture 5: Integrals and Applications Lejla Batina Institute for Computing and Information Sciences Digital Security Version: spring 2012 Lejla Batina Version: spring 2012 Wiskunde 1 1 / 21 Outline The

More information

Lecture 4: Integrals and applications

Lecture 4: Integrals and applications Lecture 4: Integrals and applications Lejla Batina Institute for Computing and Information Sciences Digital Security Version: autumn 2013 Lejla Batina Version: autumn 2013 Calculus en Kansrekenen 1 / 18

More information

Convergence & Continuity

Convergence & Continuity September 5 8, 217 Given a sequence {x n } n Ξ of elements x n X, where X is a metric space, Ξ is an indexing set. Let Ξ = N Definition (convergent sequence) A sequence {x n } is called convergent sequence

More information

1 Review of di erential calculus

1 Review of di erential calculus Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts

More information

16.2 Line Integrals. Lukas Geyer. M273, Fall Montana State University. Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall / 21

16.2 Line Integrals. Lukas Geyer. M273, Fall Montana State University. Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall / 21 16.2 Line Integrals Lukas Geyer Montana State University M273, Fall 211 Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall 211 1 / 21 Scalar Line Integrals Definition f (x) ds = lim { s i } N f (P i ) s

More information

Chapter 6 - Random Processes

Chapter 6 - Random Processes EE385 Class Notes //04 John Stensby Chapter 6 - Random Processes Recall that a random variable X is a mapping between the sample space S and the extended real line R +. That is, X : S R +. A random process

More information

Integration by Parts. MAT 126, Week 2, Thursday class. Xuntao Hu

Integration by Parts. MAT 126, Week 2, Thursday class. Xuntao Hu MAT 126, Week 2, Thursday class Xuntao Hu Recall that the substitution rule is a combination of the FTC and the chain rule. We can also combine the FTC and the product rule: d dx [f (x)g(x)] = f (x)g (x)

More information

Quasi-invariant Measures on Path Space. Denis Bell University of North Florida

Quasi-invariant Measures on Path Space. Denis Bell University of North Florida Quasi-invariant Measures on Path Space Denis Bell University of North Florida Transformation of measure under the flow of a vector field Let E be a vector space (or a manifold), equipped with a finite

More information

Math 240 Calculus III

Math 240 Calculus III Calculus III Summer 2015, Session II Monday, August 3, 2015 Agenda 1. 2. Introduction The reduction of technique, which applies to second- linear differential equations, allows us to go beyond equations

More information

Linear DifferentiaL Equation

Linear DifferentiaL Equation Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace

More information

arxiv: v4 [cond-mat.stat-mech] 3 Mar 2017

arxiv: v4 [cond-mat.stat-mech] 3 Mar 2017 Memory Erasure using Time Multiplexed Potentials Saurav Talukdar, Shreyas Bhaban and Murti V. Salapaka University of Minnesota, Minneapolis, USA. (Dated: March, 7) arxiv:69.87v [cond-mat.stat-mech] Mar

More information

macroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics

macroscopic 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 information

Gaussian Basics Random Processes Filtering of Random Processes Signal Space Concepts

Gaussian Basics Random Processes Filtering of Random Processes Signal Space Concepts White Gaussian Noise I Definition: A (real-valued) random process X t is called white Gaussian Noise if I X t is Gaussian for each time instance t I Mean: m X (t) =0 for all t I Autocorrelation function:

More information

Math 222 Spring 2013 Exam 3 Review Problem Answers

Math 222 Spring 2013 Exam 3 Review Problem Answers . (a) By the Chain ule, Math Spring 3 Exam 3 eview Problem Answers w s w x x s + w y y s (y xy)() + (xy x )( ) (( s + 4t) (s 3t)( s + 4t)) ((s 3t)( s + 4t) (s 3t) ) 8s 94st + 3t (b) By the Chain ule, w

More information

Math 266, Midterm Exam 1

Math 266, Midterm Exam 1 Math 266, Midterm Exam 1 February 19th 2016 Name: Ground Rules: 1. Calculator is NOT allowed. 2. Show your work for every problem unless otherwise stated (partial credits are available). 3. You may use

More information

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009 A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance

More information

Theoretical Tutorial Session 2

Theoretical Tutorial Session 2 1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations

More information

2.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS. differential equations with the initial values y(x 0. ; l.

2.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS. differential equations with the initial values y(x 0. ; l. Numerical Methods II UNIT.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS.1.1 Runge-Kutta Method of Fourth Order 1. Let = f x,y,z, = gx,y,z be the simultaneous first order

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

Solution to Assignment 3

Solution to Assignment 3 The Chinese University of Hong Kong ENGG3D: Probability and Statistics for Engineers 5-6 Term Solution to Assignment 3 Hongyang Li, Francis Due: 3:pm, March Release Date: March 8, 6 Dear students, The

More information

Solutions of Math 53 Midterm Exam I

Solutions of Math 53 Midterm Exam I Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior

More information

Robotics. Islam S. M. Khalil. November 15, German University in Cairo

Robotics. Islam S. M. Khalil. November 15, German University in Cairo Robotics German University in Cairo November 15, 2016 Fundamental concepts In optimal control problems the objective is to determine a function that minimizes a specified functional, i.e., the performance

More information

) k ( 1 λ ) n k. ) n n e. k!(n k)! n

) k ( 1 λ ) n k. ) n n e. k!(n k)! n k ) λ ) k λ ) λk k! e λ ) π/!. e α + α) /α e k ) λ ) k λ )! λ k! k)! ) λ k λ k! λk e λ k! λk e λ. k! ) k λ ) k k + k k k ) k ) k e k λ e k ) k EX EX V arx) X Nα, σ ) Bp) Eα) Πλ) U, θ) X Nα, σ ) E ) X α

More information

CHAPTER 3. P (B j A i ) P (B j ) =log 2. j=1

CHAPTER 3. P (B j A i ) P (B j ) =log 2. j=1 CHAPTER 3 Problem 3. : Also : Hence : I(B j ; A i ) = log P (B j A i ) P (B j ) 4 P (B j )= P (B j,a i )= i= 3 P (A i )= P (B j,a i )= j= =log P (B j,a i ) P (B j )P (A i ).3, j=.7, j=.4, j=3.3, i=.7,

More information

(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C.

(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C. Midterm 1 33B-1 015 October 1 Find the exact solution of the initial value problem. Indicate the interval of existence. y = x, y( 1) = 0. 1 + y Solution. We observe that the equation is separable, and

More information

Weak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations

Weak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations Weak Convergence of Numerical Methods for Dynamical Systems and, and a relation with Large Deviations for Stochastic Equations Mattias Sandberg KTH CSC 2010-10-21 Outline The error representation for weak

More information

First order differential equations

First order differential equations First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First

More information

The multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications

The multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications

More information

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose

More information

M343 Homework 3 Enrique Areyan May 17, 2013

M343 Homework 3 Enrique Areyan May 17, 2013 M343 Homework 3 Enrique Areyan May 17, 013 Section.6 3. Consider the equation: (3x xy + )dx + (6y x + 3)dy = 0. Let M(x, y) = 3x xy + and N(x, y) = 6y x + 3. Since: y = x = N We can conclude that this

More information

Path Integral methods for solving stochastic problems. Carson C. Chow, NIH

Path Integral methods for solving stochastic problems. Carson C. Chow, NIH Path Integral methods for solving stochastic problems Carson C. Chow, NIH Why? Often in neuroscience we run into stochastic ODEs of the form dx dt = f(x)+g(x)η(t) η(t) =0 η(t)η(t ) = δ(t t ) Examples Integrate-and-fire

More information

Math 162: Calculus IIA

Math 162: Calculus IIA Math 62: Calculus IIA Final Exam ANSWERS December 9, 26 Part A. (5 points) Evaluate the integral x 4 x 2 dx Substitute x 2 cos θ: x 8 cos dx θ ( 2 sin θ) dθ 4 x 2 2 sin θ 8 cos θ dθ 8 cos 2 θ cos θ dθ

More information

LB 220 Homework 4 Solutions

LB 220 Homework 4 Solutions LB 220 Homework 4 Solutions Section 11.4, # 40: This problem was solved in class on Feb. 03. Section 11.4, # 42: This problem was also solved in class on Feb. 03. Section 11.4, # 43: Also solved in class

More information

Stochastic Differential Equations

Stochastic Differential Equations CHAPTER 1 Stochastic Differential Equations Consider a stochastic process X t satisfying dx t = bt, X t,w t dt + σt, X t,w t dw t. 1.1 Question. 1 Can we obtain the existence and uniqueness theorem for

More information

Chapter 2 - Survival Models

Chapter 2 - Survival Models 2-1 Chapter 2 - Survival Models Section 2.2 - Future Lifetime Random Variable and the Survival Function Let T x = ( Future lifelength beyond age x of an individual who has survived to age x [measured in

More information

Math 212-Lecture 8. The chain rule with one independent variable

Math 212-Lecture 8. The chain rule with one independent variable Math 212-Lecture 8 137: The multivariable chain rule The chain rule with one independent variable w = f(x, y) If the particle is moving along a curve x = x(t), y = y(t), then the values that the particle

More information

Fluctuation relations and nonequilibrium thermodynamics II

Fluctuation relations and nonequilibrium thermodynamics II Fluctuations II p. 1/31 Fluctuation relations and nonequilibrium thermodynamics II Alberto Imparato and Luca Peliti Dipartimento di Fisica, Unità CNISM and Sezione INFN Politecnico di Torino, Torino (Italy)

More information

Heat Equation on Unbounded Intervals

Heat Equation on Unbounded Intervals Heat Equation on Unbounded Intervals MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 28 Objectives In this lesson we will learn about: the fundamental solution

More information

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12 Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

More information

Random Variables. Cumulative Distribution Function (CDF) Amappingthattransformstheeventstotherealline.

Random Variables. Cumulative Distribution Function (CDF) Amappingthattransformstheeventstotherealline. Random Variables Amappingthattransformstheeventstotherealline. Example 1. Toss a fair coin. Define a random variable X where X is 1 if head appears and X is if tail appears. P (X =)=1/2 P (X =1)=1/2 Example

More information

Separation of variables in two dimensions. Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 )

Separation of variables in two dimensions. Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 ) Separation of variables in two dimensions Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 ) Separation of variables in two dimensions Overview of method: Consider linear, homogeneous

More information

Intrinsic Noise in Nonlinear Gene Regulation Inference

Intrinsic 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 information

Starting from Heat Equation

Starting from Heat Equation Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most

More information

Find the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x

Find the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x Assignment 5 Name Find the indicated derivative. ) Find y(4) if y = sin x. ) A) y(4) = cos x B) y(4) = sin x y(4) = - cos x y(4) = - sin x ) y = (csc x + cot x)(csc x - cot x) ) A) y = 0 B) y = y = - csc

More information

= 2e t e 2t + ( e 2t )e 3t = 2e t e t = e t. Math 20D Final Review

= 2e t e 2t + ( e 2t )e 3t = 2e t e t = e t. Math 20D Final Review Math D Final Review. Solve the differential equation in two ways, first using variation of parameters and then using undetermined coefficients: Corresponding homogenous equation: with characteristic equation

More information

Integration by Substitution

Integration by Substitution November 22, 2013 Introduction 7x 2 cos(3x 3 )dx =? 2xe x2 +5 dx =? Chain rule The chain rule: d dx (f (g(x))) = f (g(x)) g (x). Use the chain rule to find f (x) and then write the corresponding anti-differentiation

More information

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior

More information

Exact Simulation of Diffusions and Jump Diffusions

Exact 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 information

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral

More information

Rough Burgers-like equations with multiplicative noise

Rough Burgers-like equations with multiplicative noise Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du

More information

WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid

WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1 I. OVERVIEW 2 The importance of splitting integrators

More information

HW2 Solutions. MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22]

HW2 Solutions. MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22] HW2 Solutions MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, 2013 Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22] Section 3.1: 1, 2, 3, 9, 16, 18, 20, 23 Section 3.2: 1, 2,

More information

for any C, including C = 0, because y = 0 is also a solution: dy

for any C, including C = 0, because y = 0 is also a solution: dy Math 3200-001 Fall 2014 Practice exam 1 solutions 2/16/2014 Each problem is worth 0 to 4 points: 4=correct, 3=small error, 2=good progress, 1=some progress 0=nothing relevant. If the result is correct,

More information

March Algebra 2 Question 1. March Algebra 2 Question 1

March Algebra 2 Question 1. March Algebra 2 Question 1 March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question

More information

Software Reliability & Testing

Software Reliability & Testing Repairable systems Repairable system A reparable system is obtained by glueing individual non-repairable systems each around a single failure To describe this gluing process we need to review the concept

More information

Calculus of Variations

Calculus of Variations ECE 68 Midterm Exam Solution April 1, 8 1 Calculus of Variations This exam is open book and open notes You may consult additional references You may even discuss the problems (with anyone), but you must

More information

Chemistry. Lecture 10 Maxwell Relations. NC State University

Chemistry. Lecture 10 Maxwell Relations. NC State University Chemistry Lecture 10 Maxwell Relations NC State University Thermodynamic state functions expressed in differential form We have seen that the internal energy is conserved and depends on mechanical (dw)

More information

Stability of Stochastic Differential Equations

Stability of Stochastic Differential Equations Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010

More information

Stochastic Processes for Physicists

Stochastic Processes for Physicists Stochastic Processes for Physicists Understanding Noisy Systems Chapter 1: A review of probability theory Paul Kirk Division of Molecular Biosciences, Imperial College London 19/03/2013 1.1 Random variables

More information

Chapter 2: Entropy and Mutual Information. University of Illinois at Chicago ECE 534, Natasha Devroye

Chapter 2: Entropy and Mutual Information. University of Illinois at Chicago ECE 534, Natasha Devroye Chapter 2: Entropy and Mutual Information Chapter 2 outline Definitions Entropy Joint entropy, conditional entropy Relative entropy, mutual information Chain rules Jensen s inequality Log-sum inequality

More information

and verify that it satisfies the differential equation:

and verify that it satisfies the differential equation: MOTIVATION: Chapter One: Basic and Review Why study differential equations? Suppose we know how a certain quantity changes with time (for example, the temperature of coffee in a cup, the number of people

More information

Linear Ordinary Differential Equations

Linear Ordinary Differential Equations MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R

More information

16. Working with the Langevin and Fokker-Planck equations

16. 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 information

C2 Differential Equations : Computational Modeling and Simulation Instructor: Linwei Wang

C2 Differential Equations : Computational Modeling and Simulation Instructor: Linwei Wang C2 Differential Equations 4040-849-03: Computational Modeling and Simulation Instructor: Linwei Wang Part II Variational Principle Calculus Revisited Partial Derivatives Function of one variable df dx

More information

Differential Equations Class Notes

Differential Equations Class Notes Differential Equations Class Notes Dan Wysocki Spring 213 Contents 1 Introduction 2 2 Classification of Differential Equations 6 2.1 Linear vs. Non-Linear.................................. 7 2.2 Seperable

More information

ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.

ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D. 4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R

More information

MATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:

MATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section: MATH 5 Final Examination May 3, 07 FORM A Name: Student Number: Section: This exam has 6 questions for a total of 50 points. In order to obtain full credit for partial credit problems, all work must be

More information

SDE Coefficients. March 4, 2008

SDE 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 information

BMIR Lecture Series on Probability and Statistics Fall 2015 Discrete RVs

BMIR Lecture Series on Probability and Statistics Fall 2015 Discrete RVs Lecture #7 BMIR Lecture Series on Probability and Statistics Fall 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University 7.1 Function of Single Variable Theorem

More information

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS Math 5B Final Exam Review Session Spring 5 SOLUTIONS Problem Solve x x + y + 54te 3t and y x + 4y + 9e 3t λ SOLUTION: We have det(a λi) if and only if if and 4 λ only if λ 3λ This means that the eigenvalues

More information

MATH20411 PDEs and Vector Calculus B

MATH20411 PDEs and Vector Calculus B MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411

More information

21-256: Partial differentiation

21-256: Partial differentiation 21-256: Partial differentiation Clive Newstead, Thursday 5th June 2014 This is a summary of the important results about partial derivatives and the chain rule that you should know. Partial derivatives

More information

An Introduction to Bessel Functions

An Introduction to Bessel Functions An Introduction to R. C. Trinity University Partial Differential Equations March 25, 2014 Bessel s equation Given p 0, the ordinary differential equation x 2 y +xy +(x 2 p 2 )y = 0, x > 0 (1) is known

More information

Distinguish between. and non-thermal energy sources.

Distinguish between. and non-thermal energy sources. Distinguish between System & Surroundings We also distinguish between thermal We also distinguish between thermal and non-thermal energy sources. P Work The gas in the cylinder is the system How much work

More information

The Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative

The Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative 1 Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f

More information

RESEARCH ARTICLE. Spectral Method for Solving The General Form Linear Fredholm Volterra Integro Differential Equations Based on Chebyshev Polynomials

RESEARCH ARTICLE. Spectral Method for Solving The General Form Linear Fredholm Volterra Integro Differential Equations Based on Chebyshev Polynomials Journal of Modern Methods in Numerical Mathematics ISSN: 29-477 online Vol, No, 2, c 2 Modern Science Publishers wwwm-sciencescom RESEARCH ARTICLE Spectral Method for Solving The General Form Linear Fredholm

More information

Comparative Statics. Autumn 2018

Comparative Statics. Autumn 2018 Comparative Statics Autumn 2018 What is comparative statics? Contents 1 What is comparative statics? 2 One variable functions Multiple variable functions Vector valued functions Differential and total

More information

1 + x 2 d dx (sec 1 x) =

1 + x 2 d dx (sec 1 x) = Page This exam has: 8 multiple choice questions worth 4 points each. hand graded questions worth 4 points each. Important: No graphing calculators! Any non-graphing, non-differentiating, non-integrating

More information

1.1. BASIC ANTI-DIFFERENTIATION 21 + C.

1.1. BASIC ANTI-DIFFERENTIATION 21 + C. .. BASIC ANTI-DIFFERENTIATION and so e x cos xdx = ex sin x + e x cos x + C. We end this section with a possibly surprising complication that exists for anti-di erentiation; a type of complication which

More information

Differentiation Shortcuts

Differentiation Shortcuts Differentiation Shortcuts Sections 10-5, 11-2, 11-3, and 11-4 Prof. Nathan Wodarz Math 109 - Fall 2008 Contents 1 Basic Properties 2 1.1 Notation............................... 2 1.2 Constant Functions.........................

More information

A Quantum Particle Undergoing Continuous Observation

A Quantum Particle Undergoing Continuous Observation A Quantum Particle Undergoing Continuous Observation V.P. Belavkin and P. Staszewski y December 1988 Published in: Physics Letters A, 140 (1989) No 7,8, pp 359 {362 Abstract A stochastic model for the

More information

or we could divide the total time T into N steps, with δ = T/N. Then and then we could insert the identity everywhere along the path.

or we could divide the total time T into N steps, with δ = T/N. Then and then we could insert the identity everywhere along the path. D. L. Rubin September, 011 These notes are based on Sakurai,.4, Gottfried and Yan,.7, Shankar 8 & 1, and Richard MacKenzie s Vietnam School of Physics lecture notes arxiv:quanth/0004090v1 1 Path Integral

More information

Sec. 1.1: Basics of Vectors

Sec. 1.1: Basics of Vectors Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x

More information

dt 2 The Order of a differential equation is the order of the highest derivative that occurs in the equation. Example The differential equation

dt 2 The Order of a differential equation is the order of the highest derivative that occurs in the equation. Example The differential equation Lecture 18 : Direction Fields and Euler s Method A Differential Equation is an equation relating an unknown function and one or more of its derivatives. Examples Population growth : dp dp = kp, or = kp

More information

Solving differential equations (Sect. 7.4) Review: Overview of differential equations.

Solving differential equations (Sect. 7.4) Review: Overview of differential equations. Solving differential equations (Sect. 7.4 Previous class: Overview of differential equations. Exponential growth. Separable differential equations. Review: Overview of differential equations. Definition

More information

AMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC

AMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC Lecture 13 Copyright y Hongyun Wang, UCSC Recap: Fitting to exact data *) Data: ( x j, y j ), j = 1,,, N y j = f x j *) Polynomial fitting Gis phenomenon *) Cuic spline Convergence of cuic spline *) Application

More information