Скачиваний:
122
Добавлен:
17.03.2015
Размер:
282.76 Кб
Скачать
(p−p0)2 2

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

Волновую функцию любого состояния можно представить в виде су-

перпозиции плоских волн де-Бройля:

 

 

ψ(~r, t) =

1

Z

d3p ϕ(p,~ t) e

i

(p~r−Et).

(3.2)

 

(2πh¯)3/2

В сущности это разложение есть не что иное, как разложения в интеграл Фурье, где коэффициенты разложения ϕ(p,~ t) определяются следующим

образом:

1

 

 

 

 

ϕ(p,~ t) =

Z

d3p ψ(~r, t) e

i

(p~r−Et).

 

(2πh¯)3/2

Функция ϕ(p,~ t) может быть интерпретирована как волновая функция частицы в импульсном представлении. Причем из условия нормировки (3.1) следует, что функция ϕ(p,~ t) также нормирована:

Z

| ϕ(p,~ t) |2 d3p = 1. (3.3)

Задание 1. Найти среднюю скорость движения и ширину волнового пакета ψ(x, t), если амплитуда ϕ(p) имеет вид гауссовой кривой:

ϕ(p) = N e,

где p0 и σ – параметры распределения. Найти область локализации частицы в момент времени t = 0.

Решение. Из условия нормировки (3.3) определим коэффициент N . Интеграл по p легко вычисляется, если ввести новую переменную интегри-

рования u = (p − p0)/σ:

 

Z

 

du = N 2 σ π = 1,

N 2

Z

e

σ2

dp = N 2σ

e−u

 

 

(p p0)2

 

2

 

 

 

 

−∞

 

 

 

−∞

 

 

 

 

q

откуда N = 1/ σ π. Используя разложение (3.2), запишем функцию

ψ(x, t):

 

 

 

 

 

 

1

 

i

 

 

 

 

 

 

 

(p p0)2

 

ψ(x, t) =

 

Z

dp e

(px−Et) e

2

.

(3.4)

2¯hσπ3/2

 

 

 

−∞

 

 

 

 

 

 

21

С учетом связи между кинетической энергией и импульсом E = p2/2m, выражение для фазы может быть преобразовано к виду:

 

i

 

 

 

 

 

 

 

 

 

p

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

2

 

 

 

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a p

 

 

 

 

 

 

 

 

 

 

 

 

(px

Et)

(p −

 

20)

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ c +

 

 

 

,

(3.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

4a

 

где введены следующие обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

=

 

1

 

 

+

 

it

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2mh¯

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

=

 

 

p0

 

+

 

i

 

x,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

=

 

 

 

 

p02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя результат (3.5) в выражение (3.4), получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dp e(a p−

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(x, t)

=

 

 

ec+ 4a

 

Z

 

 

 

 

2

 

)

 

 

=

 

 

 

 

 

 

(3.6)

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

2¯hσπ3/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

b2

 

 

 

=

 

 

 

 

 

 

 

ec+

4a

 

Z

 

dz e−z

 

 

=

 

 

 

 

 

 

 

ec+

4a .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3/2

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2¯hσπ

 

 

 

 

 

 

 

 

 

 

q2¯hσaπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем общее решение для ψ(x, t), выделив в нем в явном виде плоскую волну де-Бройля. Для этого добавим и вычтем в показателе экспоненты слагаемое вида (p0x −E0t), где E0 = p20/2m. С учетом явного вида коэффициентов a, b, c фаза может быть приведена к виду:

 

b2

 

i

 

σ2

c +

 

± (px − E0t) =

 

 

(p0x − E0t) −

 

(x − v0t)2,

4a

2αh¯2

где

α = 1 + itσ2 . mh¯

Подставляя этот результат в формулу (3.6), выпишем окончательное

выражение для функции ψ(x, t) :

 

 

 

 

ψ(x, t) = A e

i

(p0x−E0t),

(3.7)

где амплитуда волнового пакета A:

 

 

 

 

1

 

 

 

 

 

σ2

 

2

 

A =

 

 

 

 

e

2αh¯

2

(x−v0t) .

(3.8)

q

2αhσ¯

 

 

 

π

 

 

Здесь v0 = p0/m – скорость частицы. Из последней формулы видно, что

22

j

(x; t) j

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2 h

p

 

 

 

 

 

 

 

 

 

 

1 j

 

 

 

 

 

 

 

 

 

 

2hp j j e

 

h

p

j j

24x

 

p

j j

 

 

 

v0t

 

v0t

v0t +

h

 

x

 

 

 

 

 

 

Рис. 3. Плотность вероятности группы волн как функция x для некоторого момента времени t

центр тяжести волнового пакета (точка x = v0t) движется со скоростью частицы.

График плотности вероятности | ψ(x, t) |2 такого волнового пакета в некоторый фиксированный момент времени t изображен на рис. 3. Из этого рисунка, в частности, видно, что функция имеет максимум в точке x0 = v0t и и быстро спадает с ростом отклонения x от x0. Удвоенное расстояние от точки максимума до точек, в который квадрат модуля волновой функции убывает в e раз, можно принять за эффективную ширину волнового пакета

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

q|

α

2

 

σ2t

 

 

x

 

|

 

=

 

u1 +

 

 

.

(3.9)

 

 

 

 

σ

 

 

σ u

mh¯

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

Ширина (3.9) с течением времени растет, и, следовательно, пакет постепенно расплывается.

Чтобы найти область локализации частицы в начальный момент времени, надо найти плотность вероятности обнаружения частицы в точке x:

1

e

σ2

x

2

ρ =| ψ(x, 0) |2=

 

 

|α|22

.

2σh¯q

 

 

π | α |2

 

Максимум этой функции – точка x = 0 соответствует наиболее вероятному месторасположению частицы в момент времени t = 0.

23

Задания для самостоятельного решения

1. В момент времени t = 0 поведение частицы описывается функцией

ψ(x, 0) = N e xa22 + i p0x .

Найти нормировочный коэффициент и область локализации частицы в начальный момент времени.

2.Найти коэффициенты Фурье ϕ(p, t) для функции, приведенной в задаче 1, и определить ширину пакета в импульсном представлении.

3.Пронормировать следующие функции (здесь a – постоянная):

3.1ψ = N e−x2/2a2 ,

3.2ψ = N x e−x2/2a2 ,

3.3ψ = N (2x2 − a2) e−x2/2a2 ,

3.4ψ = N esin θ,

3.5ψ = N (3 cos2 θ − 1),

3.6ψ = N e−iϕ sin θ r e−r/2a,

3.7ψ = N cos θ r e−r/2a.

4.Используя определение волновой функции в импульсном представ-

лении

 

1

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

d3p ϕ(p,~ t) e

i

 

ψ(~r, t) =

 

 

 

 

 

(p~r−Et),

 

 

 

 

 

(2πh¯)3/2

a) доказать, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(p,~ t) =

1

Z

 

d3p ψ(~r, t) e

i

(p~r−Et),

 

 

 

(2πh¯)3/2

 

б) проверить, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

| ψ(~r, t) |2 dV

 

= Z

| ϕ(p,~ t) |2

d3p.

5. Разложить на плоские волны волновую функцию

а)

ψ(x, 0) = b,

 

 

 

 

 

 

 

 

 

 

 

 

 

−a ≤ x ≤ a,

 

 

 

 

0,

 

 

x <

a, x > a.

 

 

 

 

 

 

 

 

 

 

 

 

 

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

б) То же для функции

 

b e

i

p0x −a ≤ x ≤ a,

ψ(x, 0) =

 

 

0,

 

x <

a, x > a.

 

 

 

 

 

 

24

в) То же для функции

ψ(x, 0) = N e2xa22 + i p0x.

4.Собственные значения и собственные функции оператора

ˆ

Для того чтобы в состоянии ψ величина, изображаемая оператором L, имела единственное значение, необходимо, чтобы выполнялось условие:

ˆ

(4.1)

L ψL = L ψL.

Волновая функция, удовлетворяющую этому условию, называется соб-

ˆ

ственной функцией оператора L, числовой коэффициент L называется

ˆ

собственным значением оператора L. Спектр собственных значений может быть дискретным, непрерывным и смешанным. Каждому собственному значению L1, L2, . . . , Ln соответствует вполне определенная волновая функция ψ1, ψ2 . . . , ψn. В этом случае говорят, что собственная функция ψn принадлежит собственному значению Ln.

Собственные функции ψn и ψm самосопряженного оператора, принадлежащие разным собственным значениям Ln и Lm, ортогональны между

собой:

Z

ψn ψm dV = 0.

Примером уравнения на собственные значения и собственные функции является стационарное уравнение Шредингера. Если гамильтониан

ˆ

системы H не зависит явно от времени, то общее решение имеет вид:

ψn(~r, t) = ψn(~r) e−iEnt/h¯ .

Здесь ψn(~r) – собственные функции оператора Гамильтона:

ˆ

H ψn(~r) = En ψn(~r),

где En собственные значения оператора полной энергии. Это уравнение называют стационарным уравнением Шредингера.

Задание 1. Найти собственные значения Mz и собственные функции

ˆ

Φ(ϕ) оператора Mz .

25

Решение:

ˆ

Оператор Mz в сферической системе координат имеет вид:

ˆ

Mz = −ih¯

∂ϕ

.

Подставляя его в уравнение на собственные функции (4.1), имеем

−ih¯ ∂ϕ Φ(ϕ) = Mz Φ(ϕ).

Решением этого уравнения является функция

Φ(ϕ) = N eiMz ϕ/h¯ .

(4.2)

Из условия однозначности функции

Φ(ϕ) = Φ(ϕ + 2π)

находим спектр собственных значений Mz :

N eiMz ϕ/h¯ = N eiMz ϕ/h¯ eiMz 2π/h¯ 1 = eiMz 2π/h¯ ,

откуда

 

Mz = mh,¯

m = 0, ±1, ±2, . . . .

ˆ

Таким образом, спектр значений оператора Mz – дискретный, проекция Mz меняется с шагом h¯. Чтобы вычислить коэффициент N , воспользуемся условием нормировки:

Z

| Φ(ϕ) |2 dϕ = 1.

0

Подставляя в него функцию (4.2), находим коэффициент нормировки

 

1

 

N 2 Z

dϕ = N 2

 

2 π = 1 N =

 

,

0

 

 

 

 

 

ˆ

с учетом которого собственная функция оператора Mz имеет вид:

Φm(ϕ) = √1 eimϕ. 2π

Квантовое число m называют азимутальным квантовым числом, оно

ˆ

определяет спектр собственных значений оператора Mz , Mz = mh¯.

26

Задание 2. Найти уровни энергии и волновые функции сферически симметричного трехмерного осциллятора с потенциальной энергией вида U (r) = µω2r2/2, где µ – масса.

Решение:

Состояние системы в сферически симметричной задачи полностью определяется тремя механическими величинами: энергией, проекцией момента импульса на ось z и квадратом момента импульса. Волновая функция, описывающая такую систему, зависит от трех квантовых чисел: n – главного квантового числа, определяющего дискретные уровни энергии, l – орбитального квантового числа, определяющего значения квадрата момента импульса, m – азимутального квантового числа, которое определяет значения проекции момента импульса на ось z.

Спектр энергии и волновые функции осциллятора найдем из стационарного уравнения Шредингера

ˆ

n,l,m = Eψn,l,m.

В рассматриваемом случае сферической симметрии оно имеет вид:

 

 

2

 

2

r +

Mˆ 2

+ U (r) ψn,l,m = Eψn,l,m.

(4.3)

 

 

2

2

 

2µr ∂r

 

2µr

 

 

 

 

 

ˆ 2

 

ˆ

ˆ

 

Так как операторы M

и Mz коммутируют с H (напомним, что операто-

ˆ 2

ˆ

 

 

 

 

 

 

 

 

ры M

и Mz зависят только от углов и не зависят от расстояния r), то

функция ψ(r, θ, ϕ) может быть также и собственной функцией операто-

ˆ 2

ˆ

ров M

и Mz :

ˆ 2

M ψn,l,m(r, θ, ϕ)

ˆ

Mz ψn,l,m(r, θ, ϕ)

=

M 2ψn,l,m(r, θ, ϕ),

M 2 = h¯2l(l + 1), (4.4)

=

Mz ψn,l,m(r, θ, ϕ),

Mz = m h¯.

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

ψn,l,m(r, θ, ϕ) = R(r)Yl,m(θ, ϕ),

(4.5)

где Yl,m(θ, ϕ) – сферические функции, общие для любой сферически симметричной задачи, Yl,m(θ, ϕ) = Pl,m(cos θ)eimϕ, Pl,m(cos θ) – полином Лежанра. Напомним, что эти функции являются собственными функциями

ˆ 2

ˆ

и удовлетворяют системе уравнений (4.4).

операторов M

и Mz

27

Подставляя общий вид решения (4.5) в (4.3) с учетом (4.4), получаем уравнение для функции R(r):

 

2

 

 

2

(

 

2

 

R(r) +

 

r2R(r) = ER(r).

(4.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

dr

 

 

r

 

 

2

 

 

q

Если ввести новую функцию Φ(r) = rR(r), обозначения r0 = h/µω¯ , λ = E/hω¯ и переменную ξ = r/r0, то уравнение преобразуется к виду:

 

2

 

+ 2λ

 

2

 

ξ2 Φ(ξ) = 0.

(4.7)

d2Φ(ξ)

l(l + 1)

 

 

 

ξ

 

 

 

Рассмотрим асимптотическое поведение функции Φ(ξ). При ξ → ∞ уравнение (4.7) сводится к уравнению

Φ′′ (ξ) − ξ2Φ(ξ) = 0,

решением которого является функция Φ(ξ) = c1e−ξ2/2 + c2eξ2/2. Так как волновая функция должна убывать на бесконечности, следует положить c2 = 0. Коэффициент c1 можно взять равным 1, поскольку функция еще не нормирована. Таким образом, функция Φ(ξ) в этой асимптотике имеет

вид:

Φ(ξ)ξ→∞ e−ξ2/2.

При ξ → 0 получаем уравнение

′′

 

l(l + 1)

 

 

Φ

(ξ) −

 

 

Φ(ξ) = 0.

(4.8)

ξ2

Решение этого уравнения будем искать в виде степенной функции Φ(ξ) = ξa. Подставляя ee в (4.8), получаем следующее соотношение:

a(a − 1) = l(l + 1),

которое удовлетворяется при двух значениях :

a = l + 1, a = −l.

Однако при a = −l функция Φ(ξ) неограниченно возрастает при ξ → 0. Конечной функция будет, только если показатель степени равен l + 1:

Φ(ξ)ξ→0 ξl+1.

Принимая во внимания асимптотическое поведение функции Φ(ξ), запишем ее виде:

Φ(ξ) = e−ξ2/2 ξl+1 v(ξ),

(4.9)

28

где v(ξ) – неизвестная функция конечной степени, которая при ξ → 0 выходит на константу, а при ξ → ∞ не ”перебивает” экспотенциальное затухание. После подстановки (4.9) в выражение (4.7), получаем уравнение на функцию v(ξ):

 

d2v(ξ)

+ 2

dv(ξ)

 

l + 1

 

ξ! + 2 λ

l

 

3

! v(ξ) = 0.

2

ξ

2

 

 

 

 

 

 

Подставляя в последнее уравнение функцию v(ξ) в виде степенного ряда

 

X

 

v(ξ) =

ck ξk,

(4.10)

k=0

причем c0 =6 0, что следует из проведенного выше исследования волновой функции в нуле, и, приравнивая коэффициенты при одинаковых степенях ξ, получаем рекуррентное соотношение:

ck+2 = ck

2(k + l + 23 − λ)

.

(4.11)

 

(k + 2)(k + 2l + 3)

 

Из этого соотношения следует, что при k → ∞ ряд превращается в eξ2 , а значит, функция Φ(ξ) eξ2/2, то есть возродилось расходящееся решение. Чтобы этого избежать, оборвем ряд. Для этого достаточно потребовать, чтобы при некотором k = kmax ≡ kr выполнялось условие:

λ = kr + l +

3

, kr = 0, 2, 4, . . . .

2

Вспоминая определение параметра λ = E/hω¯ , находим разрешенные уровни энергии сферически симметричного трехмерного осциллятора:

3 ! 1 !

En = hω¯ kr + l + 2 ≡ hω¯ n + 2 ,

где n = kr +l+1, n = 1, 2, 3, . . . . Этому уровню энергии отвечает функция

ψn,l,m = N e−ξ2/2 ξl+1 vn,l(ξ) Pl,m(cos θ)eimϕ,

где функция vn,l(ξ) определяется формулой (4.10) и представляет собой полином степени kr = n−l−1, коэффициенты которого определяются рекурентным соотношением (4.11). В силу четности kr, орбитальное квантовое число l принимает четные значения при нечетном n и нечетные значения при четном n. Заметим, что в отличие от одномерного гармонического осциллятора, состояние сферически симметричного трехмерного

29

осциллятора с заданной энергией En вырождено. Кратность вырождения s равна

s =

n−1

(2l + 1) =

n−1

(2l + 1) =

n(n + 1)

.

X

X

 

1,(l

 

2

 

 

=2p+1)

 

0,(l=2p)

 

 

 

 

 

 

Задания для самостоятельного решения 1. Найти собственные значения и собственные функции операторов :

a) x + iPˆx, б) sin Mˆ z , в) cos Mˆ z , г) e

i

ˆ

 

αMz .

 

2. С помощью повышающего и понижающего операторов Kˆ +, Kˆ най-

 

 

ˆ 2

и его

ти собственные значения оператора квадрата момента импульса M

собственные функции, отвечающие состояниям с l = 0, 1.

3.С помощью повышающего и понижающего операторов aˆ+, aˆнайти спектр энергии и первые пять волновых функций одномерного гармонического осциллятора.

4.Найти уровни энергии и волновые функции одномерного гармони-

ческого осциллятора, помещенного в электрическое поле, напряженности

~

E. Заряд частицы равен q.

5.Найти уровни энергии и волновые функции частицы, двигающейся

вцентральном поле вида

U (r) = −

Ze2

 

β

 

 

.

r

r2

6.Найти уровни энергии и волновые функции частицы, находящейся

впотенциальной яме с бесконечно высокими стенками:

а)

0, −l/2 ≤ x ≤ l/2, U = ∞, x < −l/2, x > l/2.

b)

0, 0 ≤ x ≤ a, 0 ≤ y ≤ b, 0 ≤ z ≤ c,

U= ∞, x < a, x > a; y < b, y > b; z < c, z > c.

7.Найти уровни энергии и волновые функции частицы, находящейся

впотенциальной яме вида:

U =

 

U0,

x ≤ l,

l,

 

0,

l x

 

 

 

− ≤

 

 

 

 

 

 

x ≥ l.

 

 

 

U0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8. Являются ли функции

30

Соседние файлы в папке Квантовая теория_1