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

Вакарчук І.О. Квантова механіка

.pdf
Скачиваний:
357
Добавлен:
14.02.2016
Размер:
5.24 Mб
Скачать

При ν = 0 отримуємо вiдомий спектр енерґiї для кулонiвського потенцiалу в просторi N вимiрiв (див. §44):

me4

En,l = −2~2[nr + l + (N − 1)/2]2 .

Якщо вимiрнiсть простору N = 3, то отримаємо такий вираз2:

 

 

 

~2ν2

 

m

e2

~2ν

[l(l + 1) + 1]

2

 

 

 

 

 

En,l = −

 

 

n2

 

 

 

8m

2~2n2

2m

 

+

ν

e2 +

~2ν

[l(l + 1) + 1] ,

n = nr + l + 1.

 

2

2m

 

Приклад 1. Обчислити з умов квантування Бора–Зоммерфельда енерґетичнi рiвнi атома водню в просторi з мiнiмальною довжиною ~β, коли деформацiйна функцiя f = 1 + βPr2, Pr радiальна компонента iмпульсу.

Виходимо з “радiальної” умови квантування,

Z Z

dr dPr

= 2π~(nr + ν),

1 + βPr2

де nr = 0, 1, 2 . . ., а 0 ≤ ν < 1 i, повторюючи розрахунки, виконанi в Прикладi

3 до §30, зводимо лiву частину цього рiвняння до табличного iнтеґрала

 

 

 

 

 

π/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4pϕb Z0

 

 

 

 

cos2 x

 

 

 

 

 

dx

 

 

 

 

(1 + b sin2 x)(1 + 2βbm|E| sin2 x)

 

 

 

 

 

| |

 

 

| |

 

p

 

 

 

 

 

 

 

 

 

 

 

2πpϕ

 

me2

 

 

 

 

 

 

 

 

 

 

 

=

1 2βm E

 

pϕ

2m E

 

1 − 2βm|E| + β(me2/pϕ)2! .

 

 

 

2

2

 

2

 

p

 

 

 

 

ексцентриситет, момент

 

 

 

 

 

 

2

 

4

Тут b = (me

ǫ)

/pϕ2m|E|, ǫ =

 

1 − 2|E|pϕ/me

 

iмпульсу

 

 

~

 

 

 

 

 

2, . . . . Тепер умова квантування є такою:

 

 

 

pϕ = (l + 1/2), l = 0, 1p,

 

 

 

 

 

 

 

 

 

me2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ~(nr + ν)(1 − 2βm|E|) + pϕp1 − 2βm|E| + β(me2/pϕ)2.

 

p

2m|E|

 

Якщо ν = 1/2, то з цього рiвняння в лiнiйному наближеннi по β знаходимо

 

me4

 

me2

 

 

 

n

 

 

2l + 1

 

 

# ,

En,l = −

"1 − 2β

 

2

1

 

+

− 1

 

2~2n2

~

 

 

n2

2l + 1

4n

 

2Цей вираз уперше знайдено у статтi: C. Quesne, V. M. Tkachuk, J. Phys.

A 37, 4267 (2004).

801

де n = nr + l + 1 головне квантове число. Маємо тонку структуру спектра

деформацiя знiмає виродження за орбiтальним квантовим числом l. При

β → ∞, Enr = −e2/~β(2nr + 1).

Приклад 2. Обчислити з умов квантування Бора–Зоммерфельда енерґетичнi рiвнi частинки в кулонiвському полi U = −e2/r з масою, залежною вiд радiальної координати r, m = m(1 + a/r), m > 0, a > 0.

Запишемо вираз для енерґiї частинки, використовуючи сферичнi координати:

 

 

pr2

 

pϕ2

 

e2

E =

 

 

+

 

 

 

 

,

 

 

 

2m

2m r2

r

тут узагальнений iмпульс pϕ

= ~(l + 1/2), орбiтальне квантове число l =

0, 1, 2, . . . . З цього виразу знаходимо радiальну компоненту iмпульсу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pr = s2m

 

2

+ E

 

 

 

 

e

pϕ2

 

 

 

,

r

 

r2

де

pϕ2 = p2ϕ − 2mae2,

e 2 = e2 + aE.

Ми звели наш приклад до стандартної квазiкласичної задачi про атом водню (див. Приклад 3 до §30). Зв’язанi стани iснують за умови, що pϕ2 > 0,

e2 + aE > 0:

m e 4

E = − 2~2 (nr + l + 1)2 ,

де радiальне квантове число nr = 0, 1, 2, . . ., а ефективне орбiтальне квантове

число

q

l = (l + 1/2)2 − 2a/aB − 1/2.

Розв’язуючи квадратне рiвняння для енерґiї E, остаточно маємо:

 

e2

~2n 2

 

 

 

− 1! ,

En,l = −

r1 +

2a

 

+

 

 

a

ma2

n 2aB

де ефективне головне квантове число

n = nr + l + 1.

Знак перед радикалом вибираємо так, щоб при a = 0 цей вираз переходив

у формулу Бора для енерґетичних рiвнiв атома водню. З умови iснування зв’язаних станiв виходить, що a < aB/8.

802

§ 103. Проблема Кеплера в теорiї Дiрака з деформацiєю

Розглянемо задачу про рух релятивiстської частинки з масою, що залежить вiд координат, у просторi з деформованою алгеброю Гайзенберґа. Крiм загального аналiзу, ми знайдемо також точний розв’язок рiвняння Дiрака для руху частинки в кулонiвському полi для певних залежностей вiд координат маси частинки та деформацiйної функцiї 3.

Почнемо з рiвняння Дiрака для частинки з потенцiальною енерґiєю U у стандартних позначеннях:

h i

ˆ ˆ 2 ˆ

(αP)c + m c β + U Ψ = EΨ,

де ˆ ˆ матрицi Дiрака, а координати та iмпульси задовольня-

α, β

ють переставнi спiввiдношення з деформованою алґеброю Гайзенберґа,

 

[x , x ] = 0,

 

 

 

 

 

 

[xj , Pˆk ] = i~δ f,

 

 

 

 

ˆ

 

ˆ

 

ˆ

ˆ

~

 

 

 

 

j

k

jk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f

 

∂f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[Pj , Pk] = −i

 

∂xj Pk ∂xk Pj , (j, k) = 1, 2, 3;

з деформацiйною функцiєю f = f(x, y, z), яка залежить лише вiд координат частинки. Припускаємо, що й маса частинки m замiнена на деяку ефективну масу m , яка також залежить вiд коор-

динат частинки:

m = mf1, f1 = f1(x, y, z).

Запровадження в рiвняння Дiрака функцiй f та f1 означає

введення додаткових сил, якi дiють на частинку, крiм тих, що представленi функцiєю U.

Уведемо новий iмпульс:

ˆ = f−1/2 ˆ f−1/2,

p P

ˆ = f1/2 ˆf1/2,

P p

3Див. також: I. O. Vakarchuk, J. Phys. A 38, 7567 (2005).

803

так, що координати та новi iмпульси є канонiчно спряженими,

[xj, xk] = 0,

[xj, pˆk] = i~δjk,

[ˆpj, pˆk] = 0.

Тепер рiвняння Дiрака стає таким: h

f

1/2

(αˆ pˆ)f

1/2

2

ˆ

 

 

c + mc f1

β

i

+ U Ψ = EΨ.

Зробимо перетворення

¯

1/2

Ψ,

Ψ = f

 

у результатi якого рiвняння Дiрака для нової функцiй ¯ набирає

Ψ

вигляду:

h i

ˆ ˆ 2 ˆ ¯ ¯ f(αp)c + mc f1β + U Ψ = EΨ.

Ми можемо дивитись на це рiвняння, як на звичайне рiвняння

Дiрака, у якому матрицi Дiрака

αˆ

ˆ

множаться на деякi, зале-

, β

жнi вiд координат, масштабнi множники f, f

: αˆ

= fαˆ , βˆ

= f

β.ˆ

 

 

 

1

 

1

 

Компоненти матрицi ˆ та матриця ˆмiж собою антикомутують,

α β

квадрати компонент матрицi ˆ дорiвнюють 2, а квадрат ˆдо-

α f β

рiвнює f12.

Перш нiж починати розмову про точнi розв’язки рiвняння доцiльно перейти в рiвняннi Дiрака до нерелятивiстської межi для того, щоб з’ясувати властивостi функцiй f та f1. Щоб одержа-

ти рiвняння Шрединґера з попереднього рiвняння Дiрака при

c → ∞, уведемо нову функцiю ψ таким спiввiдношенням: h i

¯ ˆ ˆ 2 ˆ

Ψ = f(αp)c + mc f1β + E U ψ.

Тепер для функцiї ψ знаходимо таке рiвняння:

(

f(αˆ pˆ)f(αˆ pˆ)

+ m2c4f12 − (E − U)2

2m

 

2mc2

+

i~f( ˆ ˆ U)

α

2mc

)

i~cf ˆ ˆ

+ 2 β (α f1) ψ = 0.

804

Будемо вiдраховувати енерґiю вiд енерґiї спокою mc2,

 

 

 

E= E

mc2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

й пiсля простих перетворень одержуємо:

 

 

 

(

f(αˆ pˆ)f(αˆ pˆ)

+ U

(E− U)2

+

i~f(αˆ ˆ U)

 

 

 

2mc

 

2m

 

 

 

 

2mc2

 

 

mc2

 

 

i~cf

 

 

 

 

 

 

 

f12

 

 

 

βˆ (αˆ f1) )ψ = Eψ,

+

 

− 1 +

 

 

2

 

2

З останнiх двох доданкiв у фiгурних дужках цього рiвняння випливає умова на поведiнку функцiї f1 в нерелятивiстськiй межi. Справдi, для того, щоб швидкiсть свiтла c випала з нашого рiвняння при c → ∞, необхiдно, щоб f12 − 1 1/c2. Функцiя f1 може прямувати до одиницi при c → ∞ i швидше, нiж 1/c2, i тодi

вона не залишає жодного “слiду” в нерелятивiстськi межi. Якщо

f12

2

, c → ∞,

1 = mc2 U1

де U1 = U1(x, y, z) деяка функцiя координат, то з рiвняння для ψ знаходимо його нерелятивiстську межу:

 

f(

αˆ pˆ)f(αˆ pˆ)

+ U + U1 ψ = Eψ.

 

 

2m

Зробимо пiдстановку

 

 

 

 

 

 

 

ψ =

 

 

 

 

 

 

 

fϕ,

i, припускаючи, що функцiя f

pзалежить вiд довжини r радiус-

вектора r, пiсля простих перетворень, iз використанням властивостей матрицi α одержуємо таке рiвняння:

(

(f

1/2pˆf

1/2 2

+ U + U + U1) ϕ = Eϕ,

 

)

 

 

2m

 

 

 

 

 

 

U =

f df

ˆ ˆ

 

 

 

 

 

 

 

(SL),

 

 

 

 

mr dr

805

де ˆ = ~σˆ /2 оператор спiну частинки, σˆ = (ˆσx, σˆy, σˆz) матри-

S

цi Паулi, ˆ орбiтальний момент iмпульсу.

L

Отримане рiвняння можна трактувати як рiвняння Шрединґера для частинки, маса якої залежить вiд координат, m¯ = m/f2,

причому з точно визначеним розташуванням оператора iмпульсу та оберненої маси в операторi кiнетичної енерґiї:

ˆ

1

 

1

1

 

T =

2m¯ 1/4

 

1/4

.

Крiм того, якщо частинка має спiн, то в нерелятивiстськiй межi залишається величина U, назвемо її деформацiйною спiн-

орбiтальною взаємодiєю.

Якщо записати рiвняння для функцiї ϕ через “старий” iмпульс з P, то маємо рiвняння Шрединґера у просторi з деформованою

алґеброю Гайзенберґа:

ˆ 2

+ U + U + U1! ϕ = Eϕ.

P

2m

Отже, якщо в нерелятивiстськiй теорiї зразу стартувати зi звичайного рiвняння Шрединґера для вивчення поведiнки частинки з деформованими переставними спiввiдношеннями для координат та iмпульсiв, то ми втрачаємо внесок вiд деформацiйної спiн-орбiтальної взаємодiї U, а також можливий “залишок” U1

вiд залежностi маси частинки вiд координати.

Розглянемо тепер рух частинки в центрально-симетричних по-

лях U, f, f1

, тобто вважаємо, що цi функцiї залежать лише вiд

 

¯

вiдстанi r. Повернемось до рiвняння Дiрака для функцiї Ψ i зве-

демо його до радiального рiвняння. Для цього вводимо оператор радiального iмпульсу

r = r−1(rpˆ − i~)

i радiальну складову матрицi αˆ ,

 

 

 

αˆr = (αˆ nˆ),

n =

r

 

.

r

Далi вводимо оператор, який свого часу запровадив ще Дiрак,

hi

~ ˆ ˆ ˆ ˆ ~ K = β (σL) + ,

806

i, обчислюючи добуток

ˆ

¯

αˆrK, перетворюємо рiвняння для Ψ в таке:

fαˆrrc +

~cf

αˆrβKˆ ˆ + mc2f1βˆ + U Ψ¯ = EΨ¯ .

 

i

 

r

 

ˆ

 

 

Оператор K є iнтеґралом руху з власними значеннями

k = ± j +

1

= ±1, ±2, . . . ,

 

2

j квантове число повного моменту iмпульсу. Тому у представ-

леннi, де оператор

ˆ

 

 

K є дiагональним, радiальне рiвняння Дiрака

має вигляд:

 

 

 

fαˆrrc +

 

~cf

αˆrβkˆ + mc2f1βˆ + U − E R¯ = 0,

i

r

причому

 

 

 

 

 

 

¯

¯

 

 

 

Ψ = Y R,

сферичний спiнор, який є власною функцiєю оператора ˆ ,

Y K

¯ радiальна функцiя. Уведемо тепер нову радiальну функцiю

R

R таким спiввiдношенням:

 

f

αˆrβkˆ + mc2f1βˆ + E − U R.

R¯

= fαˆrrc + i r

Пiдставляючи цей вираз у попереднє рiвняння, знаходимо рiвняння для R:

c2(fpˆr)2 + ~2c2kfβˆ

 

d

 

f

+ m2c4f12

 

 

 

 

 

 

dr

r

 

 

~2

2f2k2

 

 

dU

 

 

 

 

df

R = 0.

+

 

c

+ fαˆrci~

 

 

− i~mc3αˆrβfˆ

1

− (E − U)2

 

r2

dr

 

dr

Роздiлення просторових змiнних вiд змiнних, що описують внутрiшнi ступенi вiльностi, можливе тут, якщо множники бiля матриць

807

ˆ

ˆ

матимуть однаковi залежностi вiд радiальної змiнної

β, αˆr та

αˆrβ

r, тобто

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

d

 

f

 

 

dU

 

 

 

 

 

 

 

 

=

 

 

,

 

 

dr

r

dr

 

 

 

 

d

 

f

 

 

 

 

df

 

 

C2

 

 

 

=

1

,

 

 

dr

r

 

dr

де C1, C2 сталi величини.

Якщо цi умови виконано, то рiвняння для R набуває вигляду:

(c2(fpˆr)2 + ~2c2Λˆf

d

 

f

 

 

 

 

 

 

 

dr

r

 

 

 

+

~2c2f2k2

+ m2c4f12 − (E − U)2)R = 0,

 

 

 

 

 

 

 

 

 

 

r2

де оператор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

ˆ

 

i

 

 

 

mc

ˆ

Λ = kβ +

 

~c

αˆrC1 − i

~

αˆrβC2.

ˆ

не залежить вiд радiальної координати i його лег-

Оператор Λ

ко звести до дiагонального вигляду. Зважаючи на властивостi ма-

ˆ

 

 

 

 

 

 

 

 

триць αˆr та β, виберемо їх так:

 

 

 

βˆ =

 

I

0

 

,

αˆr =

0 −i

.

 

 

0

−I

 

 

 

i 0

 

Тодi матриця оператора ˆ дорiвнює:

Λ

ˆ

 

 

 

k

 

 

~c

 

 

 

 

 

 

 

C1

Λ =

 

~c

 

~

 

 

 

 

C1

 

mc

 

 

 

 

 

 

 

 

 

 

 

+

 

C2

 

 

 

 

+

 

~

C2

 

 

 

mc

 

 

 

 

 

 

 

,

k

 

 

 

 

 

 

 

 

 

 

 

 

 

808

айого власнi значення

λ= ±sk2 + mcC2 2 C1 2.

~~c

Якщо працювати в представленнi, де оператор ˆ є дiагональ-

Λ

ним, то наше радiальне рiвняння для R остаточно набуває такого

вигляду:

(c2(fpˆr)2 +

~2c2λf

d

 

f

 

 

 

 

dr

r

+

 

~2c2f2k2

 

+ m2c4f12 − (E − U)2)R = 0.

 

r2

Зауважимо, що оскiльки функцiї f, f1 та U зв’язанi двома умова-

ми, то лише одна з них є незалежною, наприклад, це потенцiальна енерґiя U.

Розгляньмо тепер проблему Кеплера, тобто рух зарядженої частинки в кулонiвському полi, коли потенцiальна енерґiя

e2 U = − r ,

де e2 квадрат заряду. З рiвнянь, якi зв’язують функцiї U, f, f1,

знаходимо деформацiйну функцiю

f = 1 + νr,

де ν стала величина з розмiрнiстю, оберненою до довжини, i

функцiю

a f1 = 1 + r ,

a стала, що має розмiрнiсть довжини. Причому

C1 = −e2,

C2 = a,

 

ˆ

 

 

 

 

 

а власнi значення оператора Λ

 

 

 

 

 

λ = ±s

 

 

 

 

 

 

.

k2 + ~

~c

 

 

 

mca

2

 

e2

2

 

809

Пiсля стандартної пiдстановки

R = χr ,

де χ = χ(r), радiальне рiвняння набуває такого вигляду:

~2 d2

~2

l (l + 1) −

e 2

χ = E χ,

 

 

 

+

 

 

2m

dx2

2mr2

r

де

dx = drf ,

звiдки маємо, що

 

 

xν = ln(1 + νr),

 

 

 

 

0 ≤ x < ∞.

 

 

Зiрковi величини в радiальному рiвняннi такi:

 

 

l (l + 1) = k2

 

 

 

 

 

mca

 

2

 

 

 

 

 

 

 

 

 

e2

 

2

 

+

 

 

− λ −

 

 

,

 

 

~

 

 

 

 

~c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

2

 

2

ν

 

 

 

 

 

2

ν

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

~

k

 

 

 

~

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

=

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

λ

mc a,

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mc

 

 

 

 

 

m

 

 

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

2

 

 

 

 

 

2

4

 

 

 

~2

2

ν

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m c

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

=

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2mc

 

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ефективне орбiтальне квантове число

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− 1 ,

l =

 

1

+

 

1

|

1 =

 

 

k2 − α¯2

2

2

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

− α¯

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α¯2 = α2 mca 2 ,

~

α = e2/~c стала тонкої структури; тут верхнє значення l визначає верхнiй знак для λ, а нижнє нижнiй. Отже, радiальне

810

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]