Квантовая механика и интегралы по траекториям

Фейнман Ричард Филлипс

Задача 10.5. Используя метод, кратко описанный в задаче 10.2, покажите, что величины f и I(x) связаны соотношением

I(

x

)

=

12m

2h^2

1/2

^–

[V(

x

+Y)-V(

x

)]

e

^{– 6Y^2m/h^2}

dY

=

f

.

(10.65)

Обозначим

нашу приближённую функцию распределения через Z', а связанную с ней свободную энергию Гельмгольца через F', так что Z'=e^{– F'}. Тогда, применяя результаты задачи 10.5 и соотношение (10.61), получаем

Z'

=

(exp{

– [V(

x

)+I(

x

)]

})

d

x

exp

–

m

2h

_h

⁰

x^2

du

Dx(u)

.

(10.66)

Входящий в это выражение интеграл по траекториям уже вычислялся раньше; он имеет вид (10.46). Таким образом, можно записать

e

^{– F'}

=

mkT

2h^2

1/2

e

^{– U(x)}

d

x

,

(10.67)

где

U(

x

)

=

^–

V(

x

+y)

e

^{– 6Y^2m/h^2}

dy

12m

2h^2

1/2

,

(10.68)

а потенциал V(x) в явном виде не встречается.

Эти результаты означают, что свободную энергию F' можно приближённо вычислять классическим методом, т.е. с помощью выражения, подобного (10.48), и при этом получить хорошее приближение, если вместо V(x) использовать эффективный потенциал U(x), определяемый соотношением (10.68). Отметим, кстати, что эффективный потенциал зависит от температуры.

Потенциал U(x) представляет собой среднее значение потенциала V(x), полученное путём усреднения вокруг точки x с гауссовой весовой функцией, среднеквадратичное отклонение которой составляет (h^2/12m) ^1/2 . Если проанализировать различные неравенства, содержащиеся в нашем приближении, то мы найдём, что приближённое значение свободной энергии F' превышает её истинное значение F. Подробности этого обсуждаются в следующей главе [см. неравенство (11.9) и далее].

Задача 10.6. Покажите, что потенциал, определяемый соотношением (10.68), совпадает с «исправленным» потенциалом равенства (10.57) (т.е. с показателем экспоненты в этом равенстве), если в этом последнем разложить V в ряд Тейлора.

Задача 10.7. Проверьте справедливость нашего приближения на примере гармонического осциллятора, точное значение свободной энергии которого равно

F

_точное

=

kT

 ln

2sh

h

2kT

.

(10.69)

С

помощью эффективного потенциала U вычислите приближённое значение свободной энергии; покажите, что

U

=

m^2

2

x^2

+

h

12m

(10.70)

и

F

_прибл

=

kT

ln

h

kT

+

(h)^2

24kT

.

(10.71)

При различных значениях частоты определите свободную энергию или, ещё лучше, её отношение к величине kT. Предполагается, что дробь h/kT может, в частности, принимать значения 1, 2 и 4. Покажите, что, как и следовало ожидать, F' больше F и ошибка возрастает при уменьшении температуры. Обратите внимание, что даже если мы уходим очень далеко от классической области (т.е. когда отношение h/kT=2, так что вероятность пребывания системы в основном состоянии составляет 85%), приближённые результаты все ещё удивительно близки к истинным.

Сравните эти результаты с классическим приближением, где свободная энергия принимается равной kT ln(h/kT). Оно приводит к значениям 2F/h, что видно из табл. 1.

Таблица 1

h/kT

1

2

4

Точное значение

0,08266

0,8546

0,9906

Наше приближение

0,08333

0,8598

1,0264

Классический предел

0,00000

0,6931

0,6931

§ 4 Системы с несколькими переменными

Если система зависит от нескольких переменных, то (за исключением специальных задач, связанных с рассмотрением свойств симметрии) формулы, описывающие её поведение, получаются прямым обобщением уже изученных нами методов.

Жидкий гелий. В качестве примера рассмотрим задачу отыскания функции распределения в случае жидкого гелия. Предположим, что мы имеем N одинаковых атомов массы m, заключённых в некоторый объём. Предположим далее, что эти атомы взаимодействуют попарно; потенциал этого взаимодействия V(r_1,2) на больших расстояниях соответствует слабому притяжению, а на малых— очень сильному отталкиванию. Для наглядности можно представлять себе V(r) как потенциал, описывающий столкновение твёрдых шариков, т.е. положить

V(r)

=

0 при r>2,7A,

при r<2,7A.

(10.72)

Лагранжиан такой системы имеет вид

L

=

 

ⁱ

1

2

m

|R|^2

–

1

2

 

^i,j

V(r

_i,j

)

,

(10.73)

откуда следует, что функция распределения