Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Квантовая механика / Решения задач.pdf
Скачиваний:
124
Добавлен:
21.06.2019
Размер:
720 Кб
Скачать
☆

Задача 2 (2.3.2)

Микрочастица в прямоугольной потенциальной «яме» со «стенками» бесконечной и конечной «высоты»

Вычислите волновые функции состояний un (x) и соответствующие микрочастицы массы m, находящейся энергией

стационарных связанных энергетические уровни En

в поле сил с потенциальной

∞,

− ∞ < x < 0;

 

, 0 ≤ x ≤ L;

Φ(x) = Φ0

 

>Φ0 , L < x < ∞.

Φ1

Решение

1. Как и в задаче 1, в область пространства, где потенциальная энергия бесконечна, частица попасть не может, поскольку оттуда её как бы выталкивает сила бесконечной величины. Следовательно, вероятность пребывания частицы в этой области пространства тождественно равна нулю. Поэтому при −∞ < x < 0

u(x) = 0 .

(1)

2. В области 0 ≤ x ≤ L потенциальная энергия частицы равна константе Φ0 (это значит, что здесь сила на частицу не действует,

то есть в рассматриваемой области пространства частица свободна

24

всюду, кроме границ «ямы»), и стационарное уравнение Шрёдингера, как и в задаче 1, имеет вид

−

h2

′′

(2)

 

2m

u (x) +Φ0u(x) = Eu(x) .

 

 

 

3. В области L < x < ∞, как и в 0 ≤ x ≤ L , частица свободна, однако потенциальная энергия частицы равна другой константе Φ1 ,

и стационарное уравнение Шрёдингера имеет вид, несколько отличный от (2):

−

h2

′′

(3)

 

2m

u (x) +Φ1u(x) = Eu(x) .

 

 

 

Поскольку при

E ≥Φ1 связанные состояния

частицы в

рассматриваемом поле сил невозможны (реализуются только состояния рассеяния, т.е. в терминах классической механики имеет место инфинитное движение частицы), то в соответствии с условием задачи E <Φ1. В справедливости данного утверждения мы убедимся, решив эту задачу.

В зависимости от того, в каком соотношении находятся E и Φ0 ,

возможны три варианта решения уравнения (2). Однако в условиях данной задачи, как в задаче 1, случаям E =Φ0 и E <Φ0 отвечают только тривиальные решения, т.е. состояния частицы с такими

25

энергиями не реализуются. (Студентам предлагается это проверить

самостоятельно).

Поэтому далее решение задачи выполняется для значений энергии, лежащих в диапазоне

Φ0 < E <Φ1.

Обозначив

k2 =

2m(E −Φ0 )

> 0 ,

(4)

h2

 

 

 

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

′′

2

u(x).

(5)

u (x) = −k

 

Выберем для определённости положительное значение действительной величины k:

k = 2m(E −Φ0 )

> 0 .

h2

 

Общее решение дифференциального уравнения (5):

u(x) = C1sin(kx) +C2cos(kx),

(6)

26

где C1 и C2 — произвольные постоянные.

Обозначим, далее,

κ2 =

2m(Φ1 − E)

> 0

(7)

h2

 

 

 

(κ — греческая «каппа»; κ ≠ k!), и представим уравнение (3) в виде

′′

2

u(x).

(8)

u (x) =κ

 

Выберем для определённости положительное значение действительной величины κ:

κ = 2m(Φ1 − E)

> 0.

h2

 

Общее решение дифференциального уравнения (8):

u(x) = D1exp(κx) + D2exp(−κx) ,

(9)

где D1 и D2 — произвольные постоянные.

4. Далее определим фигурирующие в (6) и (9) произвольные постоянные.

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

27

пространственной координаты, то с учётом (1) на левой границе интервала 0 ≤ x ≤ L функция (6) должна обращаться в нуль:

u(−0) = u(+0) = 0.

(10)

Подставив (6) в (10), получим соотношение

C1 0 +C2 1 = 0 ,

из которого следует, что C2 = 0. Произвольная постоянная C1

может иметь любое значение, отличное от нуля — в противном случае решение уравнения окажется тривиальным.

Таким образом, вместо (6) с учётом (10) получим:

u(x) = Csin(kx) ;

0 ≤ x ≤ L

(11)

(индекс оставшейся произвольной постоянной C1 ликвидирован за ненадобностью).

Обратившись к интервалу L < x < ∞, обратим внимание, что первое слагаемое решения (9) неограниченно возрастает при приближении к правому краю этого интервала. Поскольку волновая функция должна быть конечной (а также непрерывной и однозначной) при любых значениях своих аргументов, то произвольную постоянную D1 в (9) следует положить равной нулю и таким способом «уничтожить» недопустимое частное решение.

28

Второе же слагаемое (9) в этом пределе не только конечно, но обращается в нуль, так что D2 ≠ 0.

В результате вместо (9) получим

u(x) = Dexp(−κx) ; L < x < ∞

(12)

(индекс у оставшейся произвольной постоянной D2 также ликвидирован за ненадобностью).

Хотя искомая волновая функция, состоящая из трёх «кусков» (1), (11), (12), ещё содержит две произвольные постоянные C и D, видно, что она уже удовлетворяет краевым условиям

limu(x) = 0

;

limu(x) = 0

,

(13)

x→−∞

x→+∞

которые обеспечивают связанное состояние частицы с источником поля (центром силы). В самом деле: условия (13) означают, что обнаружить частицу вдали от удерживающего её центра силы тем менее вероятно, чем она дальше от него.

5. Волновая функция должна быть всюду непрерывной. Кроме того, непрерывной, вообще говоря, должна быть и производная от волновой функции по координатам (координате). Исключением является лишь ситуация, когда в некоторой точке потенциальная энергия скачком обращается в бесконечность. В этой точке волновая функция остаётся непрерывной и, следовательно, равной

29

нулю, т.к. частица не может попасть в область пространства, где её потенциальная энергия бесконечна. Однако производная от волновой функции может в такой точке измениться скачком от конечного значения до нуля. Рассматриваемая ситуация имеет место в данной задаче на левой границе «потенциальной ямы» x = 0

—см. (10). То же самое происходит с волновой функцией частицы

взадаче 1 на обеих границах «ямы».

На правой границе x = L «ямы», рассматриваемой в данной задаче, потенциальная энергия частицы претерпевает конечный разрыв Φ1 −Φ0 . Поэтому в этой точке должны быть непрерывными как волновая функция, так и производная от неё по координате:

u(L −0) = u(L + 0); u (L −0) = u (L + 0).

(14)

′

′

 

Соотношения (14) должны использоваться с целью «сшивки» двух «кусков» волновой функции (11), (12), т.е. подбора таких значений произвольных постоянных C и D, чтобы выполнились условия непрерывности искомой волновой функции u(x) и её производной u′(x) в точке x = L (14).

Подставив (11) и (12) в (14), получим два линейных однородных алгебраических уравнения относительно неизвестных C и D:

Csin(kL) − Dexp(−κL) = 0;

Ckcos(kL) + Dκexp(−κL) = 0.

(15)

30

Для удобства введём вместо размерных величин k и κ безразмерные

ξ = kL =

L

2m(E −Φ0 ) ;

η =κL =

L

2m(Φ1 − E)

(16)

 

h

 

 

h

 

 

и перепишем систему (15), используя обозначения (16):

Csin(ξ) − Dexp(−η) = 0;

 

Cξcos(ξ) + Dηexp(−η) = 0 .

(17)

Нетривиальные решения системы уравнений (17) существуют, если её определитель равен нулю:

det =

 

sin(ξ)

−exp(−η)

 

= 0 .

(18)

 

 

 

 

ξcos(ξ)

ηexp(−η)

 

 

 

Раскрывая определитель (18), получим следующее соотношение:

ηsin(ξ)exp(−η) +ξcos(ξ)exp(−η) = 0.

(19)

Поскольку Φ0 < E < Φ1, то величина η (16) ограничена значениями 0 < η < r, где

r =

L

2m(Φ1 −Φ0 ) .

(20)

 

h

 

 

31

Поэтому экспонента в (19) всегда отлична от нуля, и её можно «сократить»:

ηsin(ξ) +ξcos(ξ) = 0.

(21)

Выражения для ξ и η (16) включают только одну величину, которой можно варьировать, чтобы «заставить» равенство (21) выполниться — энергию частицы E (остальные величины — масса частицы m, «ширина» ямы L и её «глубина» Φ1 −Φ0 — заданы условием задачи). Однако если явно подставить выражения (16) в (21), чтобы найти подходящие значения E, получится громоздкое трансцендентное уравнение относительно E, с которым крайне неудобно работать.

6. Будем рассматривать величины ξ и η как неизвестные, связанные между собой, помимо (21), соотношением, вытекающим из определений (16) и (20):

ξ2 +η2 = r2 .

(22)

Решив систему нелинейных алгебраических уравнений (21), (22),

найдём такие значения ξ и η, при которых будет выполнено условие (18) наличия у системы уравнений (17) нетривиальных решений. Может, однако, оказаться, что уравнения (21), (22) несовместны (т.е. не имеют решений). Это будет означать, что

32

условию (18) удовлетворить невозможно, и у системы уравнений

(17)нет нетривиальных решений.

7.Если пара значений (ξ,η) — решение уравнений (21), (22), то, как следует из (21), два уравнения относительно констант C и D

(17)становятся тождественными. Действительно: выразим, например, из (21)

sin(ξ) = −η1 ξcos(ξ) ,

подставим это выражение в первое уравнение (17) и умножим полученное равенство на −η . В результате получится второе уравнение (17). Это вполне естественно, т.к. условие (18) как раз и означает, что уравнения линейно зависимы.

Из первого уравнения (17) выразим D через C:

D = Csin(ξ)exp(η) .

(23)

Подставляя (23) в (12), с учётом (1), (11) получим формулу для подсчёта волновой функции частицы:

33

 

 

 

 

 

 

 

 

 

0;

−∞ < x < 0;

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

u(x) = C sin ξ

 

;

0 ≤ x ≤ L;

(24)

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

x

 

Csin(ξ)exp η 1

−

 

;

L < x < ∞.

 

 

 

 

 

 

 

L

 

Постоянная C в (24) остаётся, разумеется, произвольной. В дальнейшем, если необходимо, её можно определить из условия нормировки.

Соответствующую (24) энергию частицы подсчитаем, используя формулы (4), (16):

E =Φ0

+

h2

ξ

2

.

(25)

2mL2

 

 

 

 

 

 

 

8. Понятно, что решить систему (21), (22) аналитически невозможно. Её численное решение требует задания конкретных значений параметров задачи (массы частицы, глубины и ширины потенциальной «ямы»), что не входит в наши планы. Можно, однако, качественно проанализировать характер решений при различных сочетаниях значений параметров, используя графоаналитический подход. Этим и займёмся.

Для удобства графического построения и последующего анализа разрешим трансцендентное уравнение (21) относительно неизвестной η, для чего разделим обе его части на sin(ξ) . Получим

34

η = −ξctg(ξ).

(26)

Рассмотренная операция возможна, только если sin(ξ) ≠ 0 , т.е.

ξ ≠πn ; n = 0, 1, 2,…

(27)

Функция (26) равна –1 при ξ = 0 (проверьте!) и обращается в нуль при

ξ =π(n +1/ 2); n = 0, 1, 2,…

(28)

В точках

ξ =πn ; n = 1, 2,…

(29)

функция (26) терпит разрывы, обращаясь в +∞ при приближении к каждой точке (29) слева и в –∞ при приближении к этой точке справа.

Значения ξ и η, являющиеся решениями системы (22), (26), графически изображаются координатами точек пересечения графика функции (26) с окружностью радиуса r, уравнение которой даёт соотношение (22). Поскольку по определению величины ξ и η неотрицательны, представляют интерес только те точки пересечения, которые расположены в первом квадранте плоскости

(ξ, η).

35

Таблица функции η = −ξctgξ

 

ξ

sinξ

cosξ

ctgξ

η = –ξctgξ

 

 

 

 

–2.41421356

 

–

–0.39269908

–0.38268343

0.92387953

–0.94805945

π/8

 

 

 

 

 

 

 

 

 

∞

 

0

 

0

1

–1*

 

 

 

 

2.41421356

 

π/8

0.39269908

0.38268343

0.92387953

–0.94805945

 

 

 

 

1.00000000

 

π/4

0.78539816

0.70710678

0.70710678

–0.78539816

 

 

 

 

0.41421356

 

3π/

1.17809725

0.92387953

0.38268343

–0.48798385

8

 

 

 

 

 

 

 

 

 

0

 

π/2

1.57079632

1

0

0

 

 

 

 

–0.41421356

 

5π/

1.96349541

0.92387953

–0.38268343

0.81330642

8

 

 

 

 

 

 

 

 

 

–1.00000000

 

3π/

2.35619449

0.70710678

–0.70710678

1.96349541

4

 

 

 

 

 

 

 

 

 

–2.41421356

 

7π/

2.74889357

0.38268343

–0.92387953

6.63641614

8

 

 

 

 

 

 

 

 

 

–∞

 

π

3.14159265

0

–1

∞

 

 

 

 

 

 

* ξ → 0: sinξ ≈ ξ; cosξ ≈ 1; ctgξ ≈ 1/ξ; η = –ξctgξ ≈ –1.

В любом случае, однако, 0 < ξ < r. Тогда из (25) с учётом (20) видно, что

36

Φ0 < E <Φ1.

Таким образом, энергия микрочастицы в условиях данной задачи а) всегда больше Φ0 : микрочастица, находящаяся в связанном

состоянии, с классической точки зрения не может неподвижно «лежать» на «дне ямы», потому что даже в состоянии с минимально возможной энергией её кинетическая энергия отлична от нуля (это

— уже знакомый нам эффект «нулевых колебаний); б) всегда меньше Φ1 , т.е. связанных состояний микрочастицы

при E ≥Φ1 не существует: в этой области значений энергии возможны только состояния рассеяния.

9. В зависимости от значения r (20) возможны следующие варианты.

1) 0 < r <π / 2 : функция (26) и окружность (22) не имеют точек пересечения, так что нетривиальных решений системы уравнений

(17) нет.

Это означает, что при тех сочетаниях значений параметров задачи, которые соответствуют указанному интервалу значений r, у задачи имеется только тривиальное решение u(x) = 0 , т.е.

связанное состояние микрочастицы с заданным полем силы невозможно: слишком лёгкая частица, либо слишком «мелкая» «яма», либо слишком «узкая».

37

2) π / 2 ≤ r < 3π / 2 : функция (26) и окружность (22) пересекаются в одной точке, абсцисса которой находится в интервале значений

π / 2 ≤ ξ(1)

<π ,

(30)

1

 

 

а ордината равна

η(1)

=

r2 − (ξ(1) )2 .

(31)

1

 

1

 

Это — единственное нетривиальное решение задачи для данного сочетания её параметров. Волновую функцию u1(1) (x) ,

описывающую данное единственное связанное состояние микрочастицы с источником поля, а также соответствующий энергетический уровень E1(1) подсчитаем по формулам (24), (25) для данных значений ξ1(1) и η1(1) .

Вследствие (30) функция u1(1) (x) в интервале значений координаты x между L/2 и L имеет один максимум при x1(,1m) = (π / 2ξ1(1) )L , а при x > L экспоненциально быстро убывает до нуля.

Пусть теперь r =π / 2 . При этом ξ1(1) = r =π / 2 и E1(1) =Φ1, т.е. энергия частицы принимает граничное значение между связанными состояниями частицы и состояниями рассеяния. Максимума

38

функция u(1)

(x)

достигает на правой границе «ямы», x(1)

= L :

 

1

 

 

1,m

 

u(1)

(L) = C . Но поскольку η(1)

= 0, то при x > L волновая функция

1

 

 

1

 

 

не убывает, а остаётся постоянной и равной C. Таким образом, рассматриваемое «граничное» состояние частицы является несвязанным.

 

n

max

n

max

 

3)

 

π ≤ r ≤

 

+1 π : функция (26) и окружность (22)

 

2

 

2

 

 

 

 

пересекаются в nmax ≥1 точках, абсциссы которых находятся в интервалах значений

 

n −

1

(n

 

)

< nπ ;

 

 

π ≤ξ

n

max

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

а ординаты равны

ηn(nmax ) =

n = 1, 2,..., nmax ,

r2 − (ξn(nmax ) )2 .

Для данного сочетания параметров задачи имеется nmax

нетривиальных решений, нумеруемых числом n. Волновую функцию un(nmax ) (x) , описывающую n–е связанное состояние микрочастицы с источником поля, а также соответствующий энергетический уровень En(nmax ) подсчитаем по формулам (24), (25)

для данных значений ξn(nmax ) и ηn(nmax ) .

39

E(nmax ) >Φ

0

—

основной (наинизший) уровень энергии

1

 

 

 

 

частицы.

 

 

 

 

 

Поскольку

r >

nmax

π , то E(nmax ) <Φ1 — наивысший уровень

 

 

 

2

 

nmax

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

источником

поля. Состояния с

энергией,

большей

En(nmax ) ,

 

 

 

 

max

являются уже несвязанными (состояниями рассеяния).

 

Волновая

функция un(nmax ) (x)

при

0 < x < L

имеет n

чередующихся экстремумов (максимумов и минимумов) и n – 1 корней (нулей), не считая корня при x = 0.

Если r =

nmax

π , то

E(nmax ) =Φ1.

Такое граничное состояние

 

 

 

2

 

nmax

 

 

 

микрочастицы

является

состоянием

рассеяния, т.к.

при x ≥ L

u(nmax ) (x) = C

(см. случай n

max

=1). Поэтому в рассматриваемом

nmax

 

 

 

 

 

 

«граничном»

случае

количество

связанных

состояний

микрочастицы становится равным nmax −1.

10. В дополнение к тем квантовым эффектам поведения микрочастицы, находящейся в связанных состояниях, которые демонстрирует решение задачи 1, в данной задаче обнаруживается ещё один. Поскольку при x > L и E <Φ1 волновая функция частицы un (x) отлична от нуля, то конечна и вероятность обнаружить частицу в данной области пространства. Но с точки зрения классической физики эта область пространства является для частицы совершенно недоступной, поскольку кинетическая энергия

40

частицы K = E −Φ1 оказалась бы здесь отрицательной, а импульс и скорость — мнимыми! Этот факт, разумеется, является наблюдаемым и характерен для поведения реальных микрочастиц. Однако найти ему «разумное» объяснение, не совершая логического круга (т.е. не пытаясь объяснить один квантовый эффект через другой, столь же необъяснимый), невозможно.

11. При r >>1

ξn(nmax ) → nπ ; n = 1, 2,..., nmax ,

а nmax → ∞. Рассматриваемый предел, в частности, соответствует неограниченному увеличению глубины потенциальной «ямы» Φ1 −Φ0 при сохранении величин других параметров задачи

(ширины «ямы», массы частицы). Очевидно, при этом решение данной задачи должно перейти в решение задачи 1.

В самом деле, как видно из (24), при этом

0;

− ∞ < x < 0;

 

 

 

x

 

 

 

un (x) → C sin πn

 

; 0

≤ x ≤ L;

 

 

 

L

 

 

L < x < ∞.

 

0;

 

Одновременно из (25) следует, что

41

En →Φ0 + π 2h2 n2 . 2mL2

12. В том же пределе r → ∞ квантовые эффекты в поведении частицы становятся несущественными, и оказывается возможным его классическое описание, как показано в решении задачи 1.

Можно отметить, что рассматриваемый предел, как видно из (20), соответствует неравенству

h

<< L .

2m(Φ1 −Φ0 )

Если величину в знаменателе левой части этого неравенства считать оценкой порядка величины импульса p микрочастицы в

«яме», то сама левая часть с точностью до множителя 2π соответствует оценке длины волны де Бройля частицы λ = 2πh/ p.

При этом обсуждаемое неравенство приобретёт вид стандартного условия применимости классической механики для описания динамического поведения микрочастицы, которая находится в пространственной области с характерным размером L : λ << L .

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

42