Handbook of Stochastic Methods - 2ed - Gardiner
Handbook of Stochastic Methods - 2ed - Gardiner
+ {Aq — allt + dn(ty}?> . (1.4.28)
We now assume again that (t)dy(t)) =0 and expand, after taking averages
using the fact that (x)?
= (a + ey), — Ca + ey? = ervar{y},
to order & . (6.3.31)
The procedure can clearly be carried on to arbitrarily high order. Of course
in a one-variable system, the stationary distribution can be evaluated exactly and
the moments found by integration. But in many variable systems this is not always
possible, whereas the multivariate extension of this method is always able to be
carried out.
b) Stationary Autocorrelation Function
The autocorrelation function of x is simply related to that of y in a stationary
state by
with (8.2.10)
D=Pd. (8.2.11)
Thus, the diffusion equation is recovered. We will generalise this result shortly.
a) Kramers-Moyal or System Size Expansion Equations
We need a parameter in terms of which the numbers and transition probabilities
scale appropriately. There are two limits which are possible, both of which corres-
pond to increasing numbers of molecules:
i) limit of large cells: ! —~ co, at fixed concentration;
ii) limit of high concentration at fixed cell size.
The results are the same for pure diffusion. In either case,
tix) + 0 (8.2.12)
and a system size expansion is possible, To lowest order, this will be equivalent to
a Kramers-Moyal expansion. From (7.5.31,32) we find8.2. Multivariate Master Equation Description 309
A(x) = Y Dy x; (8.2.13)
i
Binl®) = 51m Dox + Dyx)) = DimX1 — DyuiXm » (8.2.14)
where
Dn = dn ~ 5S de (8.2.15)
and thus, in this limit, P(x, 1) obeys the Fokker-Planck equation
OP = —S2dA(e)P +} TI:DyBinl=)P - (8.2.16)
b) Continuum Form of Kramers-Moyal Expansion
The continuum form is introduced by associating a point r with a cell i and writing
PY fdr (8.2.17)
Dj —~ D(r',1) = Wr, — r+) (8.2.18)
15, (r — 1’). (8.2.19)
At this stage we make no particular symmetry assumptions on D,,, etc, so that
anisotropic inhomogeneous diffusion is included.
However, there are some requirements brought about by the meaning of the
concept “diffusion.”
i) Diffusion is observed only when a concentration gradient exists. This means that
the stationary state corresponds to constant concentration and from (8.2.13,15),
this means that
SD, =0, (8.2.20)
7
ie,
Dd = Dy. (8.2.21)
Note that detailed balance (8.2.6) implies these.
ii) Diffusion does not change the total amount of substance in the system, i.e.,
455-0 (8.2.22)
aa”
and this must be true for any value of x;. From the equation for the mean values,
this requires
xD, =0 (8.2.23)
which follows from (8.2.15) and (8.2.21)
iii) Jn the continuum notation, (8.2.20) implies that for any r,310 8. Spatially Distributed Systems
Ja Dr + 5,6) =0 (8.2.24)
and from (8.2.23), we also have
f.d°6 Hr, 6) =0. (8.2.25)
iv) If detailed balance is true, (8.2.24) is replaced by the equation obtained by
substituting (8.2.6) in the definition of D, i-e.,
Dy = Dy (8.2.26)
which gives in the continuum form
P(r + 6, —6) = Dr, 4). (8.2.27)
The derivation of a continuum form now follows in a similar way to that of the
Kramers-Moyal expansion.
We define the derivate moments
M(r) = sf P65 5D (r, 6) (8.2.28)
D(r) = $f a5 65D, 8) (8.2.29)
and it is assumed that derivate moments of higher order vanish in some appro-
priate limit, similar to those used#in the Kramers-Moyal expansion.
The detailed balance requirement (8.2.27) gives
M(r) = [ d°65O(r + 6, —8)
= f 05 [O(r, —8) + 6-PHr, —6) + ...] (8.2.30)
= —M() + 2p-D(n) +
Hence, detailed balance requires
Mr) =7-D(r). (8.2.31)
The weaker requirement (8.2.24) similarly requires the weaker condition
V-[M(r) — 7-D(r)] = 0. (8.2.32)
We now can make the continuum form of A,(x):
Aix) —~ [ 25 Hr, 6)pr + 6)
= M@)-Volr) + Die): PF pr) (8.2.33)
If detailed balance is true, we can rewrite, from (8.2.31)
Ai(x)>V [D(r)-Fa(r)) (8.2.34)8.2. Multivariate Master Equation Description 31
The general form, without detailed balance, can be obtained by defining
Jr) = M(r) — V-D(r) (8.2.35)
from (8.2.32)
V-I(r) =0 (8.2.36)
so that we can write
Ur) =0-E(), ' (8.2.37)
where E(r) is an antisymmetric tensor. Substituting, we find that by defining
H(r) = Dir) + EW), (8.2.38)
we have defined a nonsymmetric diffusion tensor H(r) and that
Afx) + 7 [H(r)-Fo(r,t)) « (8.2.39)
This means that, deterministically,
8.p(r, t) = P-LH(r)-P lr, 1)) (8.2.40)
where H(r) is symmetric, if detailed balance holds.
We now come to the fluctuation term, given by (8.2.14). To compute Bim(x),
we first consider the limit of
PY Bing f dr’ B(r, r')g(r’) (8.2.41)
where J» is an arbitrary function. By similar, but much more tedious computation,
we eventually find
PD Bingm > —2P [D(r)e(r)-73)] (8.2.42)
so that
Br, 29'P: [D(r)p(r)8(r (8.2.43)
The phenomenological theory of Sect. 8.1 now has a rational basis since (8.2.43)
is what arises from assumingsiz 8. Spatially Distributed Systems
ite.) = HOP oC, 1) — Er) (8.2.44)
in which
Cele, DEC, > = 28(t — FDL — Pole) (8.2.45)
and hence,
.o(r, 1) = 0-H(r)-Vp +P -er, t). (8.2.46)
This corresponds to a theory of inhomogeneous anisotropic diffusion without
detailed balance. This is usually simplified by setting
Din) = Dt (8.2.47)
and this gives a more familiar equation. Notice that according to (8.2.45), fluctua-
tions in different components of the current are in general correlated, unless D
is diagonal.
Comparison with Fluctuation-Dissipation Argument. The result (8.2.43) can almost
be obtained from a simple fluctuation-dissipation argument in the stationary state,
where we know the fluctuations are Poissonian. In that case,
Be AD = KOSy, (8.2.48)
corresponding to
&(r, 1) = alr), pl’) = Br — Fp) -
Since the theory is linear, we can apply (4.4.51) of Sect. 4.4.6. Here the matrices
A and AT become
A —V-H(n)-V
AT 7 Hr)". (8.2.49)
Thus,
Bir, r') = BB™
(8.2.50)
= Ao + oA* = (-7-H()-V — 9 H()-P'lelr, £')-
We note that in the stationary state, (p(r)) = dr — r')} (8.2.51)
= 279": [D(r) 0 governed by the same law of time development
of the mean [e.g.,(4.4.44)]. It is a consequence of the Markovian linear nature of
the problem. For
Ltosonde = 2 caso, 510)pde
= ((-Ax,dt + BdW(0)], x7)> (4.4.59)
and since t > 0, dW(z) is uncorrelated with x™(0), so
d
GIG.) = 4 GO). (4.4.59)
Thus, computation of G,(z) requires the knowledge of G,(0) =o and the time
development equation of the mean. This result is similar to those of Sect.3.7.4.
4.4.7 The General Single Variable Linear Equation
a) Homogeneous Case
We consider firstly the homogeneous case
t
dx = (b(d)dt + g(aW(n)|x (4.4.60)
and using the usual Ito rules, write
y =logx (4.4.61)
so that
dy = dx|x — 3(dxy}/x?
= [b(t)dt + g(t)dW(t)] — 4g(t)dt, (4.4.62)
and integrating and inverting (4.4.61), we get
x(t) = x(0) exp {f (B(¢’) — da(t’]de’ + j g(t’) w(r') (4.4.63)
0
= x(06(t) (4.4.64)
which serves to define (1).
We note that [using (4.4.8)]
Ex = You might also like