Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
random books / Жакин А.И. - Лекции по оптике и квантовой механике.pdf
Скачиваний:
30
Добавлен:
07.03.2020
Размер:
3.49 Mб
Скачать

ПРЕОБРАЗОВАНИЕ ФУРЬЕ И КВАНТОВОМЕХАНИЧЕСКИЕ ОПЕРАТОРЫ

1. Вывод формулы для преобразования Фурье

Пусть дана функция ψ(x). Считая, что ее преобразование Фурье с(p) существует, представим ее в виде интеграла Фурье

 

ψ x)(= (2πh)-1/2 c(p) exp( ipx/ h) dp

(1)

−∞

Умножая это выражение на (2πh)1/2exp(-ipx/ h) и интегрируя по x, получим

(2πh)-1/2

c(p) exp[ i( p p)x/ h]dpdx =

−∞

−∞

= c(p h) [(2π)-1 exp(i(k-k )x)dx]dk

=

−∞

−∞

 

 

 

= (2π)-1/2

f(x)exp(- i p’ x)dx

(2)

−∞

где k= p/ h. Выражение, стоящее в квадратных скобках в подынтегральном выражении есть дельта-функция, поэтому равенство (2) можно записать в виде (штрих в верхнем индексе опущен)

 

с(p) = (2π)-1/2 f(x)exp(-ipx)dx

(3)

−∞

Проведенные выкладки доказывают формулу (9) в приложении 3.

Упражнение. Пользуясь изложенной методикой, вывести формулу (12) в приложении 3, считая, что задан интеграл (11).

2. Доказательство формулы

 

 

pn =

c* (p) pn c(p) dp = ψ* (x) (-i h∂/x)n ψ(x) dx,

(4)

 

−∞

−∞

 

где n —целое положительное число; (-i h∂/x)n = (-i)n hn n/xn; звездочка в верхнем индексе обозначает операцию комплексного сопряжения.

Подставляя в левую часть формулы (4) выражение для c(p) из (3) и учитывая формулу

pn exp(-ipx/ h) = (i h∂/x)n exp(-ipx/ h),

получим

212

 

 

 

1

*

n

p

n

 

 

=

 

[

ψ

(x )exp(ipx / h)dx ] [

ψ(x) (i h∂/x) exp(-ipx/ h)dx]dp (5)

 

 

 

 

 

2π −∞

−∞

 

 

 

 

−∞

 

Интегрируя второй интеграл n раз по частям, и учитывая то, что функция ψ(x) равна нулю на бесконечности, получим

∞ ∞

ψ(x) (i h∂/x)nexp(-ipx/ h )dx = exp(-ipx/ h)ψ(x) ( - i h∂/x)nψ(x)dx

−∞ −∞

Тогда выражение (5) запишется как

 

pn =

[

ψ* (x) ( - i h∂/x)nψ(x) [(2π)-1 exp(i(x-x )k)dk] dxdx

 

−∞

−∞

−∞

где k = p/ h . Выражение в скобках является дельта-функцией, поэтому из последнего равенства следует исходное утверждение (1).

Если исходная функция зависит от трех переменных x,y,z, то аналогом формул (1), (2) будут выражения

 

 

 

 

 

 

 

 

 

ψ(x,y,z) = (2πh)—3/2

c(px,py,pz) exp( i p r/ h) dpxdpydpz

(11)

 

 

 

 

−∞ −∞ −∞

 

 

 

 

 

 

 

 

 

 

 

 

c(px,py,pz)

= (2πh)—3/2

f(x,y,z) exp( -i p r/ h ) dxdydz

 

 

 

 

 

−∞ −∞ −∞

 

 

 

 

Cовершенно аналогично доказывается формула

 

 

 

 

 

 

 

 

pxn =

c*(p)pxnc(p)dpxdpydpz =

ψ* (x) ( - i h∂/x)nψ

 

 

−∞ −∞ −∞

 

 

−∞ −∞ −∞

 

(x)dxdydz

(6)

Аналогичные формулы имеют место для pyn, pzn.

Рассмотрим пример применения полученных формул. Кинетическая энергия микрочастицы выражается как

T (p) = (px2 + py2 + pz2)/2m

213

Среднее значение кинетической энергии запишется как

 

 

 

T

=

c* (p) T (p) c(p) dpxdpydpz

 

 

−∞

−∞

−∞

Разбивая интеграл на три слагаемых, и используя в каждом интеграле формулу (6) для случая n = 2, получим

 

 

 

 

 

 

T

=

ψ* (r) T

ψ(r) dxdydz ,

 

 

 

−∞

−∞

−∞

 

 

 

 

 

= 2/x2 + 2/y2 + 2/z2

(7)

T = - ( h 2/2m) ,

 

Последнее выражение является определением оператора кинетической энергии.

ПРИЛОЖЕНИЕ 5

214

СПЕЦИАЛЬНЫЕ ФУНКЦИИ

1.Классические ортогональные полиномы [7]

Вприложениях математической физики часто встречается краевая задача вида

σ(x) p′′ (x) + τ(x) p(x) + λ p(x) = 0, a x b

(1)

p(a) = 0,

p(b) = 0 ,

 

где λ - собственное число, подлежащее определению.

Для некоторого класса функциональных зависимостей σ(x), τ(x) задача

(1) имеет дискретный набор собственных значений λn и соответствующих им собственных функций pn(x), n = 1, 2, 3,... , причем λn ± λm при n m .

Введем понятие ортогональности собственных функций. Подставляя в уравнение (1) вначале λn , pn(x), а затем λm, pm(x), получим два равенства. Введем произвольную функцию r(x), называемую весовой функцией, и умножим первое равенство на r(x)pm(x), второе на r(x)pn(x) и проинтегрируем полученные соотношения по x в пределах от х=а до х=b. Беря разность выражений, полученных после интегрирования, получим

b

b

 

(λn-λm) rpnpmdx = (pmpn – pm pn)[(rσ)—rτ] dx

(2)

a

a

 

Потребуем, чтобы весовая функция r(x) удовлетворяла уравнению

 

 

(rσ)- rτ = 0

(3)

Тогда из (2) следует

 

 

b

 

 

r(x)pn(x) pm(x) dx = 0

(4)

a

Говорят, что функции pn(x), pm(x) ортогональны с весом r(x), если выполняется условие (4).

В приложениях большую роль играют уравнения, в которых функции σ(x), τ(x), r(x) имеют специальный вид, представленный следующей таблицей :

215

(a, b)

r(x)

σ

(x)

τ

p(x)

 

 

 

(x)

 

(-1, 1)

(1 - x )(1 + x)

(1 – x )

- (s+ q+2)x+q - s

Pn(s,q)(x)

 

 

 

 

 

(0, )

xs exp(-x)

x

- x + s + 1

Lns (x)

 

 

 

 

 

(- , )

exp(-x2)

1

- 2x

Hn(x)

 

 

 

 

 

 

Оказывается, что в этом случае собственные функции pn(x) являются полиномами, причем Рn(s,q)(x) называются полиномами Якоби,

Lns(x) - полиномами Лагерра, Нn(x) - полиномами Чебышева-Эрмита.

Система полиномов {pn(x)}, удовлетворяющих условиям (3),(4),

называются классическими ортогональными полиномами. Таким образом, полиномы Якоби, Лагерра, Чебышева – Эрмита являются классическими ортогональными полиномами.

Основные свойства классических ортогональных полиномов.

1. Полином pn(x) имеет

n вещественных нулей, расположенных

на

отрезке [a, b].

 

 

2.Собственные значения задачи (1) определяются соотношением

 

λn = -

n [ τ ′ (x) + 0,5 ( n - 1) σ′′ (x) ]

(5)

2.Общая формула для определения классических ортогональных полиномов имеет следующий вид

pn(x) = An

1

 

d n

(σ n (x)r(x))

(6)

 

 

 

r(x) dxn

 

Здесь коэффициенты An определяется следующими формулами:

 

 

 

 

n 1

 

An = n!an /Ann ,

Ann = (-1)n λnk

 

 

 

 

 

k =1

 

λkm = - (n-m)[(n + m - 1) σ′′ (x)/2 + τ′ (x)],

(λn0 = λn )

216

где аn —коэффициент в полиноме pn(x) при старшей степени хn. Формула (6) называется обобщенной формулой Родрига.

4. Дифференциальное уравнение для полиномов Якоби имеет вид

(1- x2) y′′ + [ q – s - (s + q + 2)x] y+ n(n + s + q + 1) y = 0,

(7)

где n = 1, 2, 3,... —целые числа. Решения этого уравнения, называемые полиномами Якоби, выражаются как

Pn (s,q)(x) =

(1)n

 

1

 

d n

((1- x)n+s (1+ x)n+q )

2n n! (1- x)s (1+ x)q dxn

 

 

При s = q = 0 полиномы Якоби называются полиномами Лежандра и выражаются как

Pn (x) =Pn ( 0,0 )(x) =

(1)n d n

((1- x2 )n )

 

 

 

2n n! dxn

 

 

5.Дифференциальное уравнения и формулы, определяющие полиномы Лаггера, имеют следующий вид

xy′′ + (1 + s —x) y+ n y = 0 ,

Ln s (x) = y (x) =

ex xs

 

d n

(ex xn+s )

 

n!

 

 

 

 

 

dxn

 

 

6. Дифференциальное уравнение

и

 

формулы, определяющие

вид

полиномов Чебышева-Эрмита, выражаются как

 

 

y′′ - 2x y+ 2 n y = 0 ,

 

 

Hn(x) = y(x) = ( - 1)n exp( x2 )

d n

(exp(-x2 ))

(11)

dxn

 

 

 

 

 

 

2.Приложение теории классических ортогональных полиномов к решению задач квантовой механики

217

2.1. Проведем исследование уравнения Шредингера, описывающего стационарные состояния водородоподобного атома:

(- h 2/2m)ψ + U(r) ψ = Eψ, U(r) = - Ze2/r, (12)

 

1

 

 

1

 

 

1

 

2

=

 

 

 

 

r2

 

 

+

 

 

1 ,

1 =

 

 

 

sinϑ

 

 

+

 

 

r

2

 

 

 

r

2

 

 

 

 

ϕ

2

 

 

 

r

r

 

 

 

 

sinϑ ϑ

ϑ

 

 

где r, ϑ ,ϕ - сферические координаты.

 

Частное решение уравнения (12) разыскиваем методом разделения

 

переменных, полагая

 

ψ = R(r) Y(ϑ ,ϕ )

(13)

Подставляя (12) в (13), получаем два уравнения

 

h2

 

1 d

2

dR

 

h2 λ

R+ U(r)R = ER

 

 

 

 

 

 

 

r

 

 

 

+

 

 

 

 

(14)

2m r

2

 

 

 

2m r

2

 

 

 

dr

 

dr

 

 

 

 

 

 

 

 

 

 

1Y = λY

 

 

 

 

 

 

(15)

Здесь λ - константа разделения.

Уравнение (15) также решаем методом разделения переменных, представляя

Y = T(ϑ ) Ф(ϕ)

(16)

Подставляя (16) в (15), получим следующие уравнения

 

 

 

Ф′′'(ϕ) = - m2 Ф(ϕ),

 

 

(17)

1 d

dT

m2sinϑ T

 

 

 

 

 

 

sinϑ

 

 

+

λT = 0

(18)

 

 

 

sinϑ dϑ

dϑ

 

 

 

 

где m —константа разделения. Используя условие периодичности, получаем решение уравнения (17) в виде

Ф(ϕ) = exp(imϕ), m = 0, ± 1, ± 2, ± 3,...

(19)

В уравнении (18) удобно перейти к переменной

x = cosϑ . Тогда

уравнение (18) примет вид

 

(1-x2)T′′ - 2 xT- [λ - m2/(1-x2)]T = 0

(20)

218

 

Соседние файлы в папке random books