Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Как понимать квантовую механику

.pdf
Скачиваний:
190
Добавлен:
28.03.2015
Размер:
4.52 Mб
Скачать

13.5. КВАЗИКЛАССИЧЕСКОЕ ПРИБЛИЖЕНИЕ

383

C+

1

 

.

(13.36)

ψ+(x)

 

sin

¯h p0

(x − a) + ϕ0

p0

Для определения фазы приравняем

логарифмические

производные

функций ψ± в точке a:

 

 

 

 

 

 

 

 

 

ψ(a)

p0

ψ+ (a)

p0

 

 

 

π

 

 

=

 

=

 

=

 

tg ϕ0

 

ϕ0 =

4 .

 

ψ(a)

¯h

ψ+(a)

¯h

13.5.4. Квазиклассическое квантование

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

Пусть частица движется в потенциальной яме, причем¨ классически разрешенная¨ область представляет собой отрезок [a, b].

Интеграл

 

b

 

1

a

p(X) dX

¯h

дает¨ приращение фазы между точками a и b. В случае бесконечновысоких стенок в точках a и b набег фазы должен быть кратен числу π (целое число полуволн). Если потенциал вблизи точек поворота близок к линейному, то, как мы определили ранее, ширина ямы с каждой стороны эффективно увеличивается на четверть полуволны и мы получаем

 

b

 

 

 

1

a

p(X) dX +

π

= n, n = 1, 2, 3, . . . .

¯h

2

Повторим те же рассуждения, более аккуратно выписывая промежуточные формулы. В классически разрешенной¨ области мы можем записать волновую функцию двумя разными способами, которые должны быть согласованы:

 

 

 

 

 

1

x

 

 

 

 

Ca

 

 

 

π

 

 

 

 

 

 

 

 

 

ψ(x) =

 

p(x)

 

sin ¯h

 

p(X) dX +

4

=

a

384

 

 

 

 

ГЛАВА 13

 

 

 

 

 

 

 

 

x

 

 

 

 

 

=

Cb

 

1

p(X) dX −

π

 

 

 

 

 

b

4

 

 

 

p(x) sin

¯h

=

 

 

 

 

 

b

 

1

x

 

 

=

Ca

 

1

p(X) dX +

 

π

 

 

 

a

 

b

4

 

p(x) sin

¯h

¯h p(X)dX +

. (13.37)

Согласованность возможна при

Ca = ±Cb,

если разность аргументов синуса составляет целое число полупериодов:

1

a

b

 

 

p(X) dX +

π

= π(n + 1), n = 0, 1, 2, . . . .

¯h

2

На фазовой плоскости (x, p) интеграл от a до b представляет собой интеграл по полупериоду — половину площади траектории, ограниченной кривой (x(t), p(t)). Пройдя из a в b частице, чтобы замкнуть период, надо еще¨ вернуться обратно, при этом импульс будет принимать те же значения с противоположным знаком −p(x). Поэтому правило квантования обычно

пишут через интеграл по периоду

8

p(X) dX = 2π¯h(n + 1 ), n = 0, 1, 2, . . .

(13.38)

2

 

Это правило называют правилом (квазиклассического) квантования Бора – Зоммерфельда. Исторически оно предшествовало созданию последовательной квантовой теории и было одним из основных положений так называемой старой квантовой механики.

Мы можем обобщить правило Бора – Зоммерфельда, записав

8

p(X) dX = ¯h(2πn + 2[π − ϕa − ϕb]), n = 0, 1, 2, . . .

(13.39)

Здесь ϕa и ϕb — фазы волновой функции вблизи точек поворота (ϕ0 в уравнении (13.31)).

13.5.5. Спектральная плотность квазиклассического спектра

Оценим интервал между соседними уровнями энергии при условии применимости правила квазиклассического квантования Бора – Зоммерфельда.

 

 

13.5. КВАЗИКЛАССИЧЕСКОЕ ПРИБЛИЖЕНИЕ

 

385

С учетом¨ параллельности dx и p вдоль траектории запишем правило

 

Бора – Зоммерфельда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2m(E − U (x)) dl = 2π¯h

 

 

J[E, x(l)] =

8

pdx =

8

|p|·|dx| =

8

 

 

n +

2 ,

 

 

 

 

 

 

 

| |

 

 

 

 

 

 

Γ

 

Γ

dl

Γ

 

p

 

 

 

 

 

здесь J[E, x(l)] — адиабатический инвариант, как функция энергии и траектории в конфигурационном пространстве.

Проварьируем это равенство:

δJ[E, x(l)] =

δJ

δE +

8Γ

δJ

δx(l) dl = 2π¯h δn.

δE

δx(l)

 

 

 

 

 

 

 

 

 

 

=0 на классич. x(l)

Вариация по траектории для решений классических уравнений движения дает¨ нуль. Остается¨

δJ[E, x(l)] =

δJ

 

8

 

 

 

δE = δE

2m(E − U (x)) dl =

δE

∂E

8Γ

= δE

 

 

 

m

 

 

 

dl.

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

Γ

 

2m(m

− U (x))

 

 

 

1

 

 

 

|p|

v

 

 

 

 

=

 

 

 

 

Здесь v = |mp| = dtdl — скорость.

8

dl

8

δJ[E, x(l)] = δE

v

= δE dt = δE · T = 2π¯hδn,

Γ

 

Γ

T = 2ωπ — период классического движения по траектории Γ.

Пусть δn = 1, что соответствует изменению номера уровня на один, тогда δE — расстояние между уровнями:

δE · T = δE

2π

= 2π¯h

δE = ¯hω.

ω

386

ГЛАВА 13

Спектральная плотность — число уровней на единичный интервал энергии — величина, обратная к δE:

ρ(E) = δE1 = ¯1 .

δE соответствует также энергии фотона, который должна излучить частица, чтобы перейти на уровень ниже, а ω — частота этого фотона, которая оказывается равна частоте обращения частицы. Это равенство частоты обращения частицы и частоты излучаемой электромагнитной волны естественно в классической электродинамике, но в квантовой механике частота фотона связана исключительно с его энергией. В квазиклассическом пределе эти частоты совпали, т. е. предсказания квантовой механики переходят в предсказания классической теории, как и должно быть согласно принципу соответствия.

13.5.6. Квазистационарные состояния в квазиклассике

Применяя правило квантования Бора – Зоммерфельда (13.38) или (13.39), мы можем получить «лишние» состояния дискретного спектра, которых с точки зрения квантовой механики быть не должно. Эти состояния соответствуют классическому периодическому движению с энергией, для которой возможно также убегание частицы на бесконечность (см. рис. 13.5).

U(x), E

x

Рис. 13.5. Стационарные (сплошные) и квазистационарные (пунктирные) уровни.

Эти «лишные» уровни — квазистационарные состояния. В соответствии с классической теорией помещенная¨ в квазистационарное состояние система может на протяжении длительного времени оставаться в этом состоянии, однако на больших временах проявляются квантовые свойства, и система может протуннелировать через потенциальный барьер и уйти на бесконечность.

13.5. КВАЗИКЛАССИЧЕСКОЕ ПРИБЛИЖЕНИЕ

387

Время жизни квазистационарного состояния мы можем оценить, зная вероятность туннелирования (D, мы оцениваем ее¨ в разделе 13.5.7 «Квазиклассическая вероятность туннелирования») и классический период колебаний системы (T ). Если частица может убежать через обе стенки с вероятностями D1 1 и D2 1, то за период T вероятность убегания составляет D = D1 + D2. Таким образом, вероятность убегания в единицу времени (обратная времени жизни состояния τ0)

1

=

D

 

τ0 =

T

T.

τ0

T

D

Благодаря соотношению неопределенности¨ (7.9), квазистационарный уровень имеет ширину (7.13)

¯h δE0 = 2τ0 .

Зависимость от времени квазистационарного состояния включает, помимо обычной комплексной экспоненты, еще¨ и вещественную экспоненту, обес-

печивающую экспоненциальное затухание (распад) уровня с характерным временем τ0:3

 

i

 

t

 

i

(E0−i

¯h

¯h

E0t

2τ0

= ψ0 e

¯h

2τ0

)t

ψ(t) = ψ0 e

 

 

e

 

 

 

 

.

Мы видим, что для временной эволюции квазистационарного состояния энергия получает мнимую добавку:

E = E0 − i

¯h

= E0 − i δE0.

2τ0

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

3При вычислении вероятности амплитуда возводится в квадрат, так что показатели экспоненты для амплитуды и для вероятности отличаются в два раза.

388

ГЛАВА 13

состояний продукты распада обычно вылетают с энергиями, недостаточными для преодоления потенциального барьера, т. е. они вылетают благодаря туннельному эффекту.

Правило Бора – Зоммерфельда также требует поправок, если потенциальная яма разделена барьером, через который частица может туннелировать туда-сюда. Ниже такая ситуация упоминается в разделе 13.5.8 «Несколько слов об инстантонах**».

13.5.7. Квазиклассическая вероятность туннелирования

Рассмотрим в квазиклассическом приближении одномерную задачу рассеяния. Прежде всего отметим, что в классически разрешенной¨ области квазиклассическая волновая функция (S(x) с точностью до второго члена

по ¯h) (13.23)

 

 

 

 

 

 

 

 

 

 

C

 

 

 

i

 

ψ(x)

 

 

 

 

exp

±¯h p(x) dx

 

p(x)

 

описывает частицу, которая по всей оси движется в одну сторону с постоянной плотностью потока вероятности.

Таким образом, надбарьерное отражение (E > U (x)) квазиклассическим приближением (S(x) с точность до второго члена по ¯h) не описывается.

Если высота потенциального барьера больше E, то мы можем воспользоваться квазиклассическим приближением.

Мы рассмотрим случай широкого потенциального барьера, с точками входа и выхода a и b (E = U (a) = U (b)). При этом естественный масштаб ширины барьера — длина затухания волновой функции внутри него:

¯h l(x) = |p(x)|.

Поскольку масштаб l(x) внутри барьера, как правило, переменный, критерий ширины записывается через набегающую внутри барьера фазу (мнимую) волновой функции:

L = a

b

 

1

a

b

 

dX

=

|p(X)| dX 1.

 

l(X)

 

¯h

 

 

 

 

 

 

 

мера dX

 

 

 

 

L — интервал от a до b, измеренный линейкой переменной длины l(x) (в длинах затухания).

13.5. КВАЗИКЛАССИЧЕСКОЕ ПРИБЛИЖЕНИЕ

389

Для широкого барьера мы имеем коэффициент отражения, близкий к 1, т. е. суперпозиция падающей и отраженной¨ волн приблизительно задается¨ через sin, как у границы потенциальной ямы (13.33). При этом внутри барьера, как и ранее при рассмотрении потенциальной ямы, преобладает затухающая экспонента.

Величина экспоненты внутри барьера снижается в e−L раз. Поскольку эта величина связана с амплитудой вероятности, то соответствующий вклад в коэффициент прохождения составляет

D0 = e2L.

Однако мы пока не учли вклад в коэффициент прохождения предэкспо-

1

 

ненциального множителя

 

и условий сшивки в точках входа и выхода.

p(x)

Как мы уже обсуждали ранее (13.5.1 «Как угадать и запомнить ква-

зиклассическую волновую функцию»), предэкспоненциальный множитель учитывает переменную скорость частицы, летящей в переменном потенциале, тогда как экспонента задает¨ поток частиц. Таким образом, изменение предэкспоненциального множителя не дает¨ вклада в поток и коэффициент прохождения.

Условия сшивки волновой функции в точках входа и выхода могут дать дополнительные множители порядка 1:

D = D0 · Da · Db, Da, Db 1.

Если точки входа и выхода «устроены одинаково», и в окрестностях обоих потенциал может быть приближен линейной функцией, то, в силу симметрии входа в барьер и выхода из него,

 

 

Db =

 

1

.

 

 

 

 

 

 

 

 

 

 

 

Da

 

 

 

В этом случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

b

 

D = D0 = e

2L

 

 

 

(13.40)

 

 

a

 

= exp

¯h

|p(X)| dX .

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

390

ГЛАВА 13

13.5.8. Несколько слов об инстантонах**

Внимательно рассмотрим показатель экспоненты в формуле для квазиклассического коэффициента прохождения через барьер:

 

b

 

 

a

b

 

 

 

 

a

b

 

 

 

a

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

2

 

 

 

 

¯h

|p(x)| dx = −i

¯h

 

 

2m(E − U (x)) dx =

¯h

 

2m(−E + U (x)) dx.

Последнее выражение может быть переписано как

1

умножить на действие

 

 

 

 

 

 

 

 

 

 

 

¯h

 

 

 

по периоду для колебания между точками a и b с зависимостью импульса

от координаты p(x) = ± 2m(−E + U (x)):

8b

1

p

(x) dx = 2

1

p

(x) dx.

¯h

 

 

¯h

 

 

 

 

 

 

 

a

 

 

Такая зависимость p(x) может быть получена из обычной изменением знака энергии:

E → −E, U (x) → −U (x).

Мнимое действие (интеграл от мнимого импульса) можно также описать как движение с мнимой скоростью. А поскольку перемещение между точками a и b вещественно, такая скорость соответствует мнимому изменению времени.

Если вспомнить, что амплитуда вероятности, соответствующая движению с действием S, задается¨ как (3.17)

i

e ¯h S ,

то мы получаем возможность рассматривать туннелирование через барьер, не суммируя обычные (классически запрещенные)¨ траектории (вклад которых практически компенсируется, в результате чего формула сходится очень медленно), а беря одну классически разрешенную¨ траекторию с мнимым временем движения:

 

 

i

 

 

2

a

b

 

 

 

 

 

 

 

 

 

 

S

 

 

2

 

 

 

 

D =

e ¯h

 

 

,

S = ¯h

 

 

2m(E − U (x)) dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Движение через потенциальный барьер с мнимым временем называют

инстантонным движением.

13.6. СОХРАНЕНИЕ ВЕРОЯТНОСТИ И УРАВНЕНИЕ НЕПРЕРЫВНОСТИ

391

Если мы имеем потенциальную яму, разделенную¨ барьером на две половины, то туннелирование через барьер приводит к тому, что система, помещенная¨ в одну половину ямы, начинает колебаться, поочередно¨ туннелируя туда-сюда. Такое состояние называют инстантоном.

Инстантоны возникают в различных задачах теории конденсированного состояния и квантовой теории поля, они, в частности, могут возникать при спонтанном нарушении симметрии как колебания вакуумного поля.

13.6. Сохранение вероятности и уравнение непрерывности

Как мы уже писали ранее (5.1.1 «Унитарная эволюция и сохранение вероятности»), сохранение полной вероятности является одним из фундаментальных принципов квантовой теории. При этом сохранение полной вероятности (вместе с линейностью и обратимостью) приводит к унитарности эволюции замкнутой системы.

Однако полная вероятность может быть записана как интеграл от плотности вероятности в конфигурационном пространстве. В силу непрерывности уравнений квантовой механики представляется интересным переписать условие сохранения вероятности в дифференциальной форме, как уравнение непрерывности для плотности вероятности.

Таким образом, мы имеем вектор состояния, заданный как волновая функция ψ(Q) на конфигурационном пространстве. Здесь Q — совокупность обобщенных¨ координат Qn (координат в конфигурационном пространстве).

Мы знаем, что

 

(Q) = (Q)|2,

 

 

(Q) dQ = 1 = const

 

— плотность вероятности в конфигурационном пространстве.

 

Уравнение непрерывности должно иметь вид

 

 

 

 

 

+ div j = 0,

(13.41)

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

∂jn

 

 

j

 

 

 

где div j =

 

n

— дивергенция в конфигурационном пространстве от

 

n ∂Q

 

, которое задает¨ плотность потока веро-

вещественного векторного поля

 

ятности.

Стоящая перед нами задача — выразить j через ψ и показать, что найденное выражение удовлетворяет уравнению непрерывности (13.41).

392

ГЛАВА 13

13.6.1. Как угадать и запомнить плотность потока вероятности

Прежде чем приступать к строгим выкладкам, угадаем ответ. Для классического распределения частиц

j(Q) = (Q) v(Q),

где v(Q) — скорость частиц в данной точке. Для волны де Бройля

ˆ = = mv ψ.

Эта же формула приближенно¨ справедлива для квазиклассической волновой функции, но теперь v уже является функцией от координат:

ˆ (Q) ≈ mv(Q) ψ(Q).

Умножая полученную формулу на ψ (Q), получаем

ψ (Q) ˆ (Q) ≈ mv ψ(Q) ψ (Q) = mv ρ(Q).

Таким образом, для волн де Бройля мы можем написать

j(Q) = ρ(Q) v(Q) = m1 ψ (Q) ˆ (Q).

Эта же формула должна быть по крайней мере приближенно¨ справедлива для квазиклассических волновых функций, но выражение ψ (Q)ˆ (Q) в общем случае является комплексным. Поэтому возьмем¨ от получившегося выражения вещественную часть. Новая гипотеза такова:

j(Q) =

 

1

 

Re ψ (Q)ˆ (Q) =

 

1

 

(ψ (Q) ˆ (Q) + ψ(Q)(ˆ (Q)) ).

 

 

 

2m

 

 

m

 

 

 

 

Как мы убедимся далее, это и есть искомая формула для гамильтониана

ˆ

pˆ2

 

 

 

 

 

 

вида H =

 

+ U (Q). В учебниках по квантовой механике, с учетом¨ pˆ =

2m

= −i¯h , ее¨ обычно записывают в следующем виде:

 

 

 

 

 

j =

i¯h

(

 

ψ

ψ + ψ

ψ ).

(13.42)

 

 

 

 

 

 

 

 

 

 

 

2m

 

 

 

Если параметризовать волновую функцию через плотность вероятности ρ(x) = (x)2| и фазу ϕ(x) = arg ψ(x), то

ψ(x) = ρ(x) e(x) j = ρ(x) ¯h ϕ(x), (13.43)

m

скорость