<< Предыдущая стр. 7(из 18 стр.)ОГЛАВЛЕНИЕ Следующая >>
chowski, although in agreement with experiment, was clearly a highly
idealized treatment. The theory was far removed from the Newtonian
mechanics of particles. Langevin initiated a train of thought that, in
1930, culminated in a new theory of Brownian motion by L. S. Ornstein
and G. E. Uhlenbeck . For ordinary Brownian motion (e.g., carmine
particles in water) the predictions of the Ornstein-Uhlenbeck theory are
numerically indistinguishable from those of the Einstein-Smoluchowski
theory. However, the Ornstein-Uhlenbeck theory is a truly dynamical
theory and represents great progress in the understanding of Brownian
motion. Also, as we shall see later (Chapter 10), there is a Brownian
motion where the Einstein-Smoluchowski theory breaks down completely
and the Ornstein-Uhlenbeck theory is successful.
The program of reducing Brownian motion to Newtonian particle me-
chanics is still incomplete. The problem, or one formulation of it, is to
deduce each of the following theories from the one below it:

Einstein - Smoluchowski
Ornstein - Uhlenbeck
Maxwell - Boltzmann
Hamilton - Jacobi.

We shall consider the п¬Ѓrst of these reductions in detail later (Chapter 10).
Now we shall describe the Ornstein-Uhlenbeck theory for a free particle
and compare it with EinsteinвЂ™s theory.

45
46 CHAPTER 9

We let x(t) denote the position of a Brownian particle at time t and
assume that the velocity dx/dt = v exists and satisп¬Ѓes the Langevin
equation
dv(t) = в€’ОІv(t)dt + dB(t). (9.1)
Here B is a Wiener process (with variance parameter to be determined
later) and ОІ is a constant with the dimensions of frequency (inverse time).
Let m be the mass of the particle, so that we can write
d2 x dB
m 2 = в€’mОІv + m .
dt dt
This is merely formal since B is not diп¬Ђerentiable. Thus (using NewtonвЂ™s
law F = ma) we are considering the force on a free Brownian particle
as made up of two parts, a frictional force F0 = в€’mОІv with friction
coeп¬ѓcient mОІ and a п¬‚uctuating force F1 = mdB/dt which is (formally)
a Gaussian stationary process with correlation function of the form a
constant times Оґ, where the constant will be determined later.
If v(0) = v0 and x(0) = x0 , the solution of the initial value problem
is, by Theorem 8.2,
t
в€’ОІt в€’ОІt
eОІs dB(s),
v(t) = e v0 + e
0
(9.2)
t
x(t) = x0 + v(s) ds.
0

For a free particle there is no loss of generality in considering only the
case of one-dimensional motion. Let Пѓ 2 be the variance parameter of B
1
(inп¬Ѓnitesimal generator 2 Пѓ 2 d2 /dv 2 , EdB(t)2 = Пѓ 2 dt). The velocity v(t) is
Gaussian with mean
eв€’ОІt v0 ,
by (9.2). To compute the covariance, let t в‰Ґ s. Then
t s
в€’ОІt в€’ОІs
ОІt1
eОІs1 dB(s1 )
Ee e dB(t1 )e
0 0
s
= eв€’ОІ(t+s) e2ОІr Пѓ 2 dr
0
2ОІs
в€’1
в€’ОІ(t+s) 2 e
=e Пѓ .
2ОІ
47
THE ORNSTEIN-UHLENBECK THEORY OF BROWNIAN MOTION

For t = s this is
Пѓ2
(1 в€’ eв€’2ОІt ).
2ОІ
Thus, no matter what v0 is, the limiting distribution of v(t) as t в†’ в€ћ is
Gaussian with mean 0 and variance Пѓ 2 /2ОІ. Now the law of equipartition
of energy in statistical mechanics says that the mean energy of the particle
1
(in equilibrium) per degree of freedom should be 2 kT . Therefore we set

1 Пѓ2 1
m = kT.
2 2ОІ 2
That is, recalling the previous notation D = kT /mОІ, we adopt the nota-
tion
ОІkT
Пѓ2 = 2 = 2ОІ 2 D
m
for the variance parameter of B.
We summarize in the following theorem.

THEOREM 9.1 Let D and ОІ be strictly positive constants and let B be
the Wiener process on В‚ with variance parameter 2ОІ 2 D. The solution of

dv(t) = в€’ОІv(t)dt + dB(t); v(0) = v0

for t > 0 is
t
в€’ОІt
eв€’ОІ(tв€’s) dB(s).
v(t) = e v0 +
0

The random variables v(t) are Gaussian with mean

m(t) = eв€’ОІt v0

and covariance

r(t, s) = ОІD eв€’ОІ|tв€’s| в€’ eв€’ОІ(t+s) .

В‚ with in-
The v(t) are the random variables of the Markov process on
п¬Ѓnitesimal generator
d2
d 2
в€’ОІv + ОІ D 2
dv dv
48 CHAPTER 9

2
with domain including Ccom (В‚), with initial measure Оґv0 . The kernel of
the corresponding semigroup operator P t is given by

(v в€’ eв€’ОІt v0 )2
в€’1
в€’2ОІt
t
p (v0 , dv) = [2ПЂОІD(1 в€’ e exp в€’
)] dv.
2
2ОІD(1 в€’ eв€’2ОІt )

The Gaussian measure Вµ with mean 0 and variance ОІD is invariant,
P tв€— Вµ = Вµ, and Вµ is the limiting distribution of v(t) as t в†’ в€ћ.

The process v is called the Ornstein-Uhlenbeck velocity process with
diп¬Ђusion coeп¬ѓcient D and relaxation time ОІ в€’1 , and the corresponding po-
sition process x given by (9.2) is called the Ornstein-Uhlenbeck process.

THEOREM 9.2 Let the v(t) be as in Theorem 9.1, and let
t
x(t) = x0 + v(s) ds.
0

Then the x(t) are Gaussian with mean

1 в€’ eв€’ОІt
m(t) = x0 +
Лњ v0
ОІ

and covariance
D
в€’2 + 2eв€’ОІt + 2eв€’ОІs в€’ eв€’ОІ|tв€’s| в€’ eв€’ОІ(t+s) .
r(t, s) = 2D min(t, s) +
Лњ
ОІ

Proof. This follows from Theorem 9.1 by integration,
t
m(t) = x0 +
Лњ m(s) ds,
0
t s
r(t, s) =
Лњ dt1 ds1 r(t1 , s1 ).
0 0

The second integration is tedious but straightforward. QED.

In particular, the variance of x(t) is

D
(в€’3 + 4eв€’ОІt в€’ eв€’2ОІt ).
2Dt +
ОІ
49
THE ORNSTEIN-UHLENBECK THEORY OF BROWNIAN MOTION

The variance in EinsteinвЂ™s theory is 2Dt. By elementary calculus, the
absolute value of the diп¬Ђerence of the two variances is less than 3DОІ в€’1 .
In the typical case of ОІ в€’1 = 10в€’8 sec ., t = 1 sec ., we make a proportional
2
в€’8
error of less than 3 Г— 10 by adopting EinsteinвЂ™s value for the variance.
The following theorem shows that the Einstein theory is a good approxi-
mation to the Ornstein-Uhlenbeck theory for a free particle.

THEOREM 9.3 Let 0 = t0 < t1 < . . . < tn , and let

в€†t = min ti в€’ tiв€’1 .
1в‰¤iв‰¤n

Let f (x1 , . . . , xn ) be the probability density function for x(t1 ), . . . , x(tn ),
where x is the Ornstein-Uhlenbeck process with x(0) = x0 , v(0) = v0 ,
diп¬Ђusion coeп¬ѓcient D and relaxation time ОІ в€’1 . Let g(x1 , . . . , xn ) be the
probability density function for w(t1 ), . . . , w(tn ), where w is the Wiener
process with w(0) = x0 and diп¬Ђusion coeп¬ѓcient D.
Let Оµ > 0. There exist N1 depending only on Оµ and n and N2 depending
only on Оµ such that if

в€†t в‰Ґ N1 ОІ в€’1 , (9.3)

2
v0
t1 в‰Ґ N2 , (9.4)
2DОІ 2
then

|f (x1 , . . . , xn ) в€’ g(x1 , . . . , xn )| dx1 . . . dxn в‰¤ Оµ. (9.5)
В‚ n

Proof. Assume, as one may without loss of generality, that x0 = 0.
Consider the non-singular linear transformation

(x1 , . . . , xn ) в†’ (Лњ1 , . . . , xn )
x Лњ

В‚n given by
on
1
xi = [2D(ti в€’ tiв€’1 )]в€’ 2 (xi в€’ xiв€’1 )
Лњ (9.6)

for i = 1, . . . , n. The random variables w(ti ) obtained when this trans-
Лњ
formation is applied to the w(ti ) are orthonormal since Ew(ti )w(tj ) =
50 CHAPTER 9

2D min(ti , tj ). Thus g , the probability density function of the w(ti ), is
Лњ Лњ
the unit Gaussian function on В‚ . Let f be the probability density func-
Лњ
n

tion of the x(ti ), where the x(ti ) are obtained by applying the linear
Лњ Лњ
transformation (9.6) to the x(ti ). The left hand side of (9.5) is unchanged
Лњ
when we replace f by f and g by g , since the total variation norm of a
Лњ
measure is unchanged under a one-to-one measurability-preserving map
such as (9.6).
We use the notation Cov for the covariance of two random variables,
Cov xy = Exy в€’ ExEy. By Theorem 9.2 and the remark following it,

Cov x(ti )x(tj ) = Cov w(ti )w(tj ) + Оµij ,

where |Оµij | в‰¤ 3DОІ в€’1 . By (9.6),

Cov x(ti )Лњ(tj ) = Оґij + Оµij ,
Лњx

 << Предыдущая стр. 7(из 18 стр.)ОГЛАВЛЕНИЕ Следующая >>