|
ние |
зависит |
только |
от величины |
интервала |
времени |
|
t'-t: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Л (/) Л (О - Ф ( / ' — 0- |
|
|
|
(81.2) |
|
|
Введенная, таким |
образом, |
корреляционная |
функция |
|
Ф |
(s) при s = 0 равна |
Л 2 , при возрастании s она |
быстро |
|
|
|
|
стремится |
к нулю. Исходя |
из до |
|
|
|
|
пущений |
о т, необходимо |
потре |
|
|
|
|
бовать, |
|
чтобы |
функция |
O(s) |
|
|
|
|
практически уже была равна ну |
|
|
|
|
лю для значений s, еще малых по |
|
|
|
|
сравнению с т. |
|
|
|
t, t' |
|
|
|
|
Поэтому |
в |
плоскости |
|
|
|
|
(рис. |
115) подынтегральная фун |
|
|
|
|
кция |
в уравнении (81.1) |
сущест |
|
Рис. |
115. К расчету двой |
венно отличается от нуля лишь в |
|
узкой |
полоске вдоль |
диагонали. |
|
ного интеграла в уравне |
|
Если |
в |
качестве |
новых |
перемен |
|
нии |
(81.1). |
|
|
|
|
|
ных |
интегрирования |
ввести |
|
|
|
|
|
|
|
|
S = |
|
f—t; |
|
|
|
|
|
|
|
U = |
t' |
+ |
t, |
|
|
|
|
|
|
то получим |
dtdt' |
dsdu |
и |
|
|
|
|
|
|
|
|
|
|
|
^du |
J O ( s ) |
ds, |
|
|
|
|
причем ввиду быстрого стремления функции Ф к нулю с ростом s интегрирование можно производить по полоске (s от —со до -j-co; и от 0 до 2т).
Таким образом, просто получаем
|
+ о |
|
- i = |
0(s)ds=A2a, |
(81.3) |
X |
J |
|
где малое по сравнению с т время а может быть обозна чено как «время когерентности» функции A (t). Оно яв ляется мерой ширины «полосы когерентности» на рис. 115. Интегралы (79.4)
о(0 = о 0 в - р ' + е " в ' J A ( S ) e M d |
о
основного уравнения v = — . A (t) можно преобразо вать, и не используя скоростей Gx , непосредственно при помощи корреляционной функции (81.2). После возведе ния в квадрат и усреднения по многим частицам из урав нения (79.4) следует:
|
|
t t |
|
|
v*(f) = vl e-W |
+er-W [\АЦ)А(V) |
em+t)dgd£'. |
|
|
о 0 |
|
|
Используя |
обозначения |
Л (g) Л (£') =Ф(£ ' — £ ) и |
£ = s , |
= |
приводим |
интеграл в пашем прибли |
жении к виду |
|
|
|
|
2/ |
+ 0 0 |
|
|
+оо |
± - j > d a \ |
G>(s)ds = ± ( e m - |
1) f 0(s)ds, |
О |
СО' |
|
ОО |
следовательно, |
|
|
|
|
При i->oo правая часть уравнения должна стремить ся к выражению kT/tn, в связи с чем для корреляцион ного интеграла получим:
Г <t>{s)ds = —$, |
|
|
(81.4) |
J |
Ш |
|
|
|
что можно было бы установить еще из уравнений |
(81.3) |
и (78.5). Корреляция скорости v(t) |
вытекает непосред |
ственно из уравнения (79.4) в результате |
умножения на |
Оо и усреднения по многим частицам: |
|
|
v(t)v(t') = |
~v~2e-m, |
|
|
или, в несколько иной записи, |
|
|
|
|
о (0 о (О = о 5 |
в _ р к " г |
' . |
|
(81.5) |
Отсюда особо простым способом определяется сред |
нее квадратичное смещение: имеем |
х= |
[w(|)rf|, |
следо- |
|
|
|
о |
|
вательпо,
li
х2 = j ' j" v(l)v(l')d£,d%.
6 6
После усреднения и использования уравнения (81.5) при £—^'—s; %-{-%'= и для больших времен ( 6 ^ 1 ) по лучаем:
~x* |
= —\du |
f e-®sXds=~t=2 |
|
— mBt=2kTBt. |
|
2 J |
J |
P |
m |
0- c o
Сучетом kTB — D приходим к хорошо известному
выражению х2 = 2Dt.
Если же условие В t^> 1 не выполнено, то в соответ ствии с (81.5) необходимо строго рассчитать величину
t t
x*^v2 f f,
о о
Это делается без затруднений и приводит к резуль тату
Это выражение совпадает с полученным выше по ме тоду Ланжевена и уже детально обсужденным выраже нием (77.6)
82. ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ ФОККЕРА — ПЛАНКА
а) Без внешнего силового поля
Еще в § 79 статистическим методом была рассчитана ве
роятность |
W(v, |
t, do) того, что скорость частицы ко вре |
мени i лежит в интервале v, |
dv, если к моменту времени |
/ = 0 она |
имела |
скорость v0. |
Для решения подобной за |
дачи изберем теперь совсем другой путь, выводя диффе ренциальное уравнение для w (v, i) из основного урав
нения v = — fiv-r-A(t). Пусть |
|
q>(0,Ti)*l |
(82.1) |
представляет вероятность того, что по истечении време ни х скорость частицы с начальной скоростью v лежит
между |
У + Т ] и v-\-r\-\-dr\. |
Тогда для распределения ско |
ростей |
к моменту |
времени |
t-\-x справедливо выражение |
|
w{v,t-\-x)= |
j w(v |
— т 1 , ^ ф ( у — r|,r])dr). |
(82.2) |
Пусть теперь время т настолько мало, что функция ф(у, г)) значительно отличается от нуля только для та ких т], значение которых мало по сравнению с гл. В та ком случае можно разложить:
|
w (v, t + т) = |
+~~ |
|
|
|
|
|
-dw |
|
1 |
.д*а>] |
|
|
|
j " w (v, f) |
ч dv |
|
— ч' |
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
dv* |
|
|
|
|
|
X Гф (v, Т)) — |
г| |
дф |
|
i |
1 2 д2 ф |
dr\. |
|
(82.3) |
|
|
|
aw |
|
т] |
— |
|
Если |
ограничиться |
|
|
|
|
2 |
|
dv2 |
|
по ц, |
то в |
квадратичными |
|
членами |
уравнение (82.3) войдут лишь следующие интегралы: |
Ф (v, ту) dt) = 1; |
т]ф (v, т)) dr) = |
т); |
|
Л ^ ^ ' ^ |
|
*1 = 5"; |
|
^ Л |
|
ф (У, Л) dr\ = T i |
|
; |
| |
^ д ф |
К |
п) |
= |
|
|
|
|
|
2 |
|
|
2 |
|
|
|
dy] |
|
|
|
|
|
|
|
|
|
|
|
|
|
да |
|
|
5а |
|
|
|
|
|
|
|
dv2 |
|
|
|
dv2 |
|
|
|
|
|
Эти средние значения можно заимствовать из преж |
них рассмотрений. |
Согласно |
обозначениям |
уравнения |
(78.2а) Г| представляет |
|
собой |
изменение скорости |
за |
время |
т, т.е. т] = U J + I — V j . |
|
|
Тем самым |
при заданном |
v = |
= Vj |
получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т) = — г/Вт и |
<Эг| |
|
|
•Вт. |
|
|
|
|
|
|
|
|
|
|
|
|
|
dv |
|
|
|
|
|
|
|
Далее, ограничиваясь членами первого порядка по т, имеем
следовательно, 6\) 2 /ди=0 .
Таким образом, уравнение (82.3) принимает вид:
w(v,t + т) — а ф , / ) = |
ауВт + а б т — |
Н |
Вт— |
|
|
|
|
да |
/Я |
dv2 |
и в пределах т->0 переходит в уравнение |
Фоккера- |
Планка: |
|
|
|
|
|
|
dw |
а |
, о dw . akT д2 ш |
• |
(82.4) |
— = |
|
Rw + By |
hp |
dt |
|
' |
dv |
m dv2 |
|
|
28-480 |
|
|
|
|
|
433 |
Если записать уравнение (82.4) в виде
|
dw |
о |
д l |
|
, кТ dw } |
, |
|
|
, а п |
|
— |
^ 6 |
— \wv-\ |
|
|
|
|
|
(82.5) |
|
|
|
|
|
|
|
|
|
|
|
|
mo" |
|
то можно заметить, что dwldt=Q, |
если w = Ce |
l k T . Рас |
пределение Максвелла, |
естественно, |
стационарно. |
|
Полное интегрирование уравнения (82.4) удается пу |
тем изящной замены переменных. Сначала |
вместо |
v в |
качестве переменной вводим |
величину p — ve$'. |
В |
этом |
случае |
имеем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
w{v,t) |
= w{pe-V,t) |
|
= |
Y{p,t), |
|
|
|
следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
dw |
_ dY |
в/. |
(Pw _ |
dW_^t |
|
|
|
|
|
do |
dp |
' |
dv2 |
|
|
dp2 |
|
|
|
|
|
|
И |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt ~~ dt |
dp\ |
dt Jv |
|
~~ |
dt |
|
P |
dp |
' |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким |
образом, из уравнения |
|
(82.4) |
получаем: |
|
|
|
|
dt |
К |
V |
|
|
.Эр2 |
|
|
|
|
|
|
где q = р — |
|
|
|
|
|
|
|
|
|
|
|
|
\ |
m |
|
|
|
|
|
|
|
|
|
|
|
|
Теперь подставим У—%е^{, откуда
Тем самым |
имеем: dt |
dt |
1 |
|
|
|
dt |
4 |
dp2 ' |
Наконец, путем |
замены |
е1^1 |
dt — dft вводим новую шка |
лу времени, |
т.е. |
полагаем: |
|
Таким образом, уравнение (82.4) приведено к при вычной форме диффузионного уравнения
д%(Р, О) |
Э 2 х ( Р . О ) |
3d |
ар2 |
Стандартное решение этого диффузионного уравне ния имеет вид:
1 |
( Р - Р . ) 3 |
Х(Р,#;Ро) = -т-Цг е |
«*> . |
4я<70' |
|
При f>->-0 это решение переходит в %(р, 0; ро)=6(р — —ро).
Если с помощью выражений
w = Xe*>; , = i L p . о |
1 |
( ^ _ , ) : |
т2р
р = ve^; р 0 = о0
снова перейти к первоначальным переменным, то полу чим:
m(v—v„ e ~ P Q 2
' O ; |
1 2zikT ( l — e ^ ) J |
Таким образом, заново получен закон (79.7) прибли жения к распределению Максвелла, на этот раз исходя из дифференциального уравнения (82.4).
б) При наличии силового поля
Пусть теперь на частицу воздействует дополнительно внешняя сила К(х) с потенциалом Ф(х), следовательно, К(х) = —дФ/дх. Аналогично предыдущим параграфам отыщем вероятность f(x, v, t)dxdv нахождения частицы в интервале х, dx и v, dv. Если в данном случае выра жение
|
|
|
<p(x,v,t,i\)dldr\ |
|
(82.6) |
представляет |
собой |
вероятность |
изменения |
х или v |
в пределах |
времени т на |, d\ и соответственно |
т), dx\, то |
применима |
аналогия с уравнением |
(82.2): |
|
f (х, v,t |
- f |
т) = |
j ' \ f (х — I, V — Т|, t) ф (х — g, V — |
|
|
-'+ |
|
|
|
|
|
-x\,l,i\)dtdi\. |
|
(82.7) |
Так же как в (82.3) мы должны разложить подынтег ральную функцию по % и г) до квадратичных членов. По-
лучающнеся при этом средние значения определим из уравнения
|
О = |
- 6 |
У + |
|
— |
|
|
+A(t), |
|
|
|
|
|
|
т |
|
|
|
|
т. е. после |
интегрирования |
|
по t |
от /т до |
( / + 1 ) т |
[см. (78.2а)]: |
|
|
|
|
|
|
|
|
|
- v. = - |
ро, т + |
|
т + G . |
(82.8) |
К тому же при т, меньшем |
времени торможения 1/|5, |
имеем: |
|
|
|
|
|
|
|
|
|
|
|
x . ^ - x ^ v . x . |
|
(82.9) |
Используя замены £ = |
X j |
+ i |
|
r |
i = U j + i — v j t |
из уравне |
ний (82.8) и |
(82.9) |
находим |
теперь средние |
значения |
|
— |
|
|
|
(усреднение производится по многим частицам с равны ми Vj = vl), ограничиваясь линейными членами по т:
Ч = ( - Р * - I - ^ L ) г,
\ т j dv
|
|
i f |
= G2 = |
|
рт, I = vx, |
|
= т. |
|
|
|
|
|
|
от |
|
|
|
di> |
|
|
|
Остальные |
средние |
значения, так же как £2 , |т) и т. п., |
при линейном приближении не вносят |
никакого |
вкла |
да |
в т. |
|
|
|
|
|
|
|
|
|
|
|
Следовательно, при разложении в ряд подынтеграль |
ного выражения (82.7) и выполнении |
интегрирования |
остается |
|
|
|
|
|
|
|
|
|
f(*,v; |
t + x) = |
f { x , v - t ) - - ^ i + ± ^ L ^ - f ^ - . |
|
|
|
|
|
|
dv |
|
2 dv2 |
|
dv |
|
|
Вводя указанные выше средние значения, произведя |
деление на т и перейдя |
к Игл т-*-0, |
|
получим |
уравнение |
Фоккера — Планка в виде |
|
|
|
|
|
|
^ |
^ - |
P |
^ + |
^ ^ |
- ^ |
- |
A |
W |
- |
f . |
(82.10) |
at |
|
Ov \ |
in |
dv j |
dx |
|
m |
dv |
|
|
|
|
Для проверки этого уравнения прежде всего устано |
вим, что известное |
распределение |
|
|
|
|
|
|
|
|
|
|
— |
0= + Ф (X) |
|
|
|
|
|
|
|
f = Ce |
kT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
действительно стационарно, т. е. обращает в нуль пра вую часть уравнения (82.10).
Если |
из уравнения |
(82.10) нужно получить |
данные |
об изменении во времени плотности р, то вначале |
напра |
шивается |
определить р с помощью выражения |
|
|
p(x,t) = |
§f(x,v,f)dv. |
(82.11) |
Если выполнить это интегрирование в уравнении (82.10), то производные по v не дадут никакого вклада, так как \(х, v, t) при |и|->оо очень быстро стремится к нулю. Получим
|
= |
— f vf(v, x)dv = |
— [pv), |
|
|
dt |
дх J v |
' |
dx K |
' |
|
т. е. хотя |
и правильный, |
но, к сожалению, крайне три |
виальный |
результат. Он выражает лишь |
закон |
сохране |
ния массы. |
|
|
|
|
|
На основании |
прежних результатов |
можно |
ожидать |
выражения для плотности в виде |
|
|
|
|
К (х) |
kT |
dp |
|
|
|
|
• р |
• |
дх |
|
|
|
|
mfi |
m(3 |
|
|
где B = l/mB — подвижность и D = £77mB—коэффициент диффузии; следовательно, надеемся получить дифферен циальное уравнение
= |
? _ Ш £ !р |
_ _ £ _ . * - ) . |
(82.12) |
dt |
дх { пф |
mf> дх j |
|
Как установил Крамере1 , уравнение (82.12) действи тельно вытекает как приближенное уравнение из выра жения (82.10):
df |
^ил |
dv |
|
L\(fV |
+ lLJL |
+ д Т |
V |
dt |
\ |
dx I \ |
dx dv |
m$ |
dx |
|
K |
(X) |
f)~ |
— |
\ m\°> |
m(3 дар |
(82.13) |
|
|
пф |
j |
dx |
|
Смысл данного преобразования выясняется при рас смотрении рис. 116: для того чтобы получить плотность р в точке х0, интегрировали функцию f(x, v) вдоль пунк тирной линии (параллельно оси v) согласно уравнению (82.11). Рассмотрим теперь наклонную прямую g, про-
1 Kramers Н. А. — «Physica», 1940, v. 7, p. 284.
ходящую через |
точку |
х0 |
и имеющую |
уравнение |
= Р(х0 —х). Вдоль этой |
прямой dv = —$dx. |
При |
движе |
нии вдоль этой прямой любая функция |
ср(х, |
v) |
изменя |
ется на |
|
|
|
|
|
|
|
|
|
|
dv |
дер |
дер |
|
|
dv |
дх |
dv |
dx |
p |
|
|
|
И Л И
|
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P |
a |
dx |
<P = |
P |
dv |
J |
вдоль |
g. |
|
|
|
|
|
|
|
|
dv |
|
|
|
|
|
|
|
|
|
|
|
Таким образом, если проинтегрировать уравнение |
|
(82.13) |
не вдоль пунктирной |
|
вертикали, |
а вдоль |
наклон |
|
|
|
|
|
|
|
ной |
прямой |
g, |
то |
первое |
|
сла |
|
|
|
|
|
|
|
гаемое в правой части урав |
|
|
|
|
|
|
|
нения |
станет |
равным |
|
нулю. |
|
|
|
|
\ 1 |
|
Выясним, |
при |
каких |
обстоя |
|
|
|
|
\ l |
|
тельствах |
подобное |
интегриро |
|
|
|
|
у |
|
вание |
по |
|
наклонной |
прямой |
|
|
|
|
хЛ |
|
может |
рассматриваться |
в |
ка |
|
|
|
|
|
l\ |
|
честве |
замены |
интегрирования |
|
|
|
|
|
|
|
при |
постоянном |
х. |
|
|
|
|
|
|
|
|
|
|
|
|
|
Функция |
f(x, |
v, |
t) |
практи |
|
|
|
|
|
|
|
чески становится равной |
нулю, |
|
Рис. 116. |
Замена |
интег |
если v превышает максималь |
|
рирования |
по |
v |
при |
по |
ное значение v0. |
На |
прямой g |
|
стоянном |
х0 |
интегриро |
в точке w0 |
координата |
х отсто |
|
ванием |
вдоль |
наклонной |
ит |
от х0 на |
у0 /Р- Но |
а0 /Р |
пред |
|
прямой |
g. |
о = |
р(х„ — |
х). |
|
ставляет |
собой |
длину |
пути |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
торможения |
частицы, |
движу |
|
щейся |
со скоростью |
v0 |
И подвергающейся только |
|
воздействию |
трения р. |
Если |
теперь К(х) |
и |
[(х, |
v) |
настолько слабо зависят от х, что на расстояниях поряд ка пути торможения самой быстрой частицы они еще мо гут рассматриваться почти постоянными, то в самом деле можно заменить интегрирование остальных членов в уравнении (82.13) вдоль g интегрированием при по
стоянном Хо. И в этом случае было бы точно |
получено |
требуемое уравнение (82.12). |
|
При |
термическом равновесии v0 имеет |
порядок |
У kT/tn. |
Следовательно, справедливости |
уравнения |
(82.12) можно ожидать лишь в том случае, когда выпол нены следующие условия:
Ж г / ^ ± |
< < К » |
J |
L |
l / |
^ i « / . |
(82.14) |
dx • V т |
& |
дх |
у |
т |
р |
|
Обсуждение |
процессов, для которых данные |
условия |
не выполняются, можно найти в работе Крамерса и Чан-
драсекара |
(см. сноску на стр. 429). |
|
8». СПЕКТРАЛЬНОЕ |
РАЗЛОЖЕНИЕ |
СТАТИСТИЧЕСКОЙ |
ФУНКЦИИ' |
В предыдущих |
параграфах мы охарактеризовали нере |
гулярное |
ускорение |
A(t) |
через среднюю |
квадратичную |
|
т |
|
|
|
|
скорость |
— ^ A { t ) d t ^ |
и |
корреляционную |
функцию |
A(t)A(t+s) |
6 |
|
|
|
|
=<D(s). |
|
|
|
Во многих случаях, в особенности если данное уско рение влияет на колебательную систему, целесообразно получить сведения о спектральном распределении. Для этой цели требуются некоторые математические вы кладки.
Под статистической функцией A(t) понимаем функ ции, подобные скорости v(t) при броуновском движении частицы или силе воздействия окружающих молекул на данную частицу. В частности, предполагаем, что среднее
значение A(t) за большой интервал времени |
равно ну |
лю, или в более строгой |
формулировке |
|
|
Ita ( |
i f Л ( • ) » ) - о |
|
|
|
о |
|
|
и |
далее, что хотя A(t) |
и изменяется |
крайне |
хаотично, |
но |
за очень большие |
промежутки |
времени |
характер |
функции остается в сущности без изменения, так что при
не слишком |
малых |
значениях т среднее |
квадратичное |
значение — |
j * A2(t)dt |
= А2 |
не зависит от t\. |
|
U |
|
|
|
Для дальнейшего |
расчета |
предположим также, что |
A(t) отличается от нуля лишь в пределах |
конечного, да- |
1 Becker R, — «Z. angew. Phys.», 1954, Bd 6, S. 23].