Master equation

guillefix 4th November 2016 at 2:43pm

For discrete space Stochastic processes

Discrete time master equation

For discrete time, probability to be in state nn at time t+Δtt+\Delta t is:

P(n,t+Δt)=nW(nn)P(n,t)P(n, t+\Delta t) = \sum_{n'} W(n'\rightarrow n)P(n', t)

where the W(nn)W(n'\rightarrow n) are the transition probabilities (which can be expressed as a transition matrix).

Continuous time master equation

For continuous time, we can subtract P(n,t)P(n,t) from both sides of the discrete time equation, and divide by Δt0\Delta t \rightarrow 0. Then

dP(n,t)dt=nw(nn)P(n,t)[nw(nn)]P(n,t)\frac{d P(n, t)}{dt} = \sum_{n'} w(n' \rightarrow n) P(n', t) - [\sum_{n'} w(n \rightarrow n')]P(n,t)

=nnw(nn)P(n,t)[nnw(nn)]P(n,t)= \sum_{n' \neq n} w(n' \rightarrow n) P(n', t) - [\sum_{n' \neq n} w(n \rightarrow n')]P(n,t)

where w(nn)=limΔt0W(nn)Δtw(n' \rightarrow n) = \lim_{\Delta t \rightarrow 0} \frac{W(n'\rightarrow n)}{\Delta t}, and where for the bracketed part we used that probability is conserved (i.e. the particle has to go somewhere), nw(nn)=1\sum_{n'} w(n \rightarrow n') = 1, and in the second line we cancelled the n=nn' = n terms from both terms.

Solve using Fourier series, as if it is in a (discrete) lattice. For more general networks, Fouriers may not be appropriate.. You can then use eigenvector methods