Решение типичных задач в курсе квантовой механики Министерство образования Украины

advertisement
1
Министерство образования Украины
НТУУ “КПИ”
Решение типичных задач в курсе квантовой
механики
часть 2
Киев 2002
2
Рецензенты
Д-р ф.-м. наук Олейник В.П.
Канд. ф.-м. наук Гусева О.А.
Терентьева Ю.Г.
Решение типичных задач в курсе квантовой механики.
НТУУ «КПИ» 2002
3
Содержание
Часть 2. Типичные задачи.
Частица в одномерной прямоугольной потенциальной яме
с
бесконечными
стенками.
…..…………………….…………….…4
Частица в сферической потенциальной яме……………………...14
Туннельный эффект………………………………………………..19
Самостоятельная
расчетная
работа
«Определение
коэффициента туннелирования частицы через потенциальный
барьер»….23
4
Раздел 2. Типичные задачи
Частица в одномерной прямоугольной потенциальной яме
Постановка задачи
Пусть частица находится в одномерной прямоугольной
потенциальной яме с бесконечно высокими стенками. Ширина ямы l.
Найти:
1. Нормированные волновые функции стационарных состояний и
энергетический спектр частицы, используя уравнение Шредингера.
2. Возможные значения энергии частицы, исходя из того, что реализуются
лишь такие ее состояния, при которых в пределах данной ямы
укладывается целое число полуволн де Бройля.
3. Плотность уровней (dN/dE) (т.е. число уровней на единичный интервал
энергии), как функцию Е, при условии, что ширина ямы такова, что
уровни расположены весьма густо. Вычислить (dN/dE) для электрона
при Е = 1,0 ЭВ, и ширине ямы
l = 1,0 см.
4. С помощью соотношения неопределенностей Гейзенберга силу
давления частицы на стенки этой ямы при минимально возможной ее
энергии. Для оценок массу частицы взять равной массе электрона.
5. Вероятность пребывания частицы в области l/3<x<2l/3, если она
находится в основном энергетическом состоянии.
6. Среднее значение координаты частицы <x> и проекции импульса <px>
если частица находится в основном состоянии.
Анализ условия
- "частица массы m" означает, что по условию считается известной масса
частицы. Этой частицей может быть любой микрообъект: электрон,
протон, нейтрон или, например, -частица.
5
- "находится в потенциальной яме" - этот термин не является чисто
квантовым. Применительно к системе, описываемой классической
механикой, это означает, что частица находится в такой области
пространства, в которой действуют потенциальные силы. Причем эти
силы направлены так, чтобы удерживать частицу в данной области. В
квантовой механике понятие «силы» не существует. Квантовая
механика вообще не ищет ответа на вопрос – какая причина приводит к
изменению в состоянии частицы?, она рассматривает вероятность того,
что состояние будет именно таким. Поэтому под потенциальной ямой
подразумевается область пространства, в которой потенциальная
энергия частицы задается некоторой функцией координат –
потенциалом.
- " потенциальная яма одномерна " - следовательно, потенциальная
энергия частицы зависит только от одной-единственной координаты.
Этот термин тоже не является чисто квантовым. Если этой координатой
является одна из декартовых осей, говорят о прямоугольной
потенциальной яме (для сравнения - в классической механике потенциальная энергия силы тяготения вблизи поверхности Земли,
которая зависит только от высоты h расположения тела над землей
U(h)=mgh). Если же этой координатой является r – расстояние до
центра поля, говорят о сферической потенциальной яме (в
электростатике - потенциальная энергия поля точечного заряда
U(r)=kq/r, она зависит только от расстояния до источника поля).
- "бесконечно высокие стенки " в потенциальной яме - означает, что в
ограниченной области пространства потенциальная энергия принимает
конечные значения (в частном случае - ноль), а за пределами этой
области - потенциальная энергия бесконечно велика. Последнее
утверждение означает, что чтобы попасть в область бесконечно
большого потенциала частица должна вначале получить откуда-нибудь
дополнительную энергию, равную бесконечно большой величине. Эта
ситуация математически возможна, а физически - нет, поэтому
становится очевидно, что находиться в области бесконечно большого
потенциала частица не может вообще. Если графически изобразить
аналитическую зависимость потенциала от координаты
, x  l , x  0
U ( x)  
const , 0  x  l
получится следующая картинка:
U(x)
6
const
0
l
x
Рис.2
В точках х=0 и х=l имеется скачкообразное изменение потенциала,
который, таким образом, графически имеет форму незамкнутого
прямоугольника. Отсюда и термин "прямоугольная" потенциальная яма.
Поскольку потенциал по своей природе есть величина, определенная с
точностью до постоянной, а аналитический вид решения принципиально
не зависит от величины const, то обычно выбирают const = 0. В этом случае
расчеты упрощаются.
Решение
Задание I
Найти нормированные волновые функции стационарных состояний и
энергетический спектр частицы, используя уравнение Шредингера.
Для начала запишем стационарное уравнение Шредингера в общем
виде,

2
( x, y, z )  U ( x, y, z ) ( x, y, z )  E( x, y, z )
2m
а теперь посмотрим, какой вид оно примет с учетом условия нашей задачи.
Поскольку наша задача одномерна (потенциал зависит только от одной
координаты), то волновая функция частицы, находящейся в таком
пространстве, тоже зависит только от этой одной координаты. В таком
случае от оператора Лапласа, который в декартовой системе координат
представляет собой сумму вторых частных производных по всем
координатам, остается только вторая производная по единственной
координате х

 2 d 2  ( x)
 U ( x) ( x)  E ( x)
2m dx 2
Учитывая вид потенциала, разобьем пространство (в данном случае
ось ОХ) на три области: область I : х є (-,0), где потенциал равен
бесконечно большой величине, область II : х є (0,l), где потенциал равен
константе (мы договорились, что константа равна нулю), и область III :
7
х є (l ,) - где потенциал снова равен бесконечности. Далее следует решить
уравнение Шредингера в каждой области отдельно и, наконец, пользуясь
непрерывностью волновой функции во всем пространстве, сшить
полученные решения между собой на границах областей.
Уравнения Шредингера для каждой области выглядят следующим
образом:
 2 d 2 I ( x)
   I ( x)  EI ( x)
2m dx 2
 2 d 2 II ( x)
область II: 0  x  l

 0  II ( x)  EII ( x)
2m dx 2
 2 d 2 III ( x)
область III: x  0

   III ( x)  EIII ( x)
2m dx 2
область I:
x0

(1)
(2)
(3)
Тот факт, что результирующая волновая функция должна быть
непрерывна во всем пространстве, означает, что ее куски
I ( x), II ( x), III ( x) непрерывны каждый в своей области и, кроме того, в
пограничных точках, где стыкуются области, значения соответствующих
кусков волновой функции слева и справа должны совпадать. Поэтому
граничными условиями в точках х = 0 и х = l являют-ся следующие
соотношения
I (0)  II (0)
(4)
II (l )  III (l )
В уравнениях (1) и (3) в качестве второго слагаемого входит
бесконечно большая величина. Уравнение Шредингера, выражая
универсальный закон сохранения энергии, должно выполняться в любой
части пространства. Единственный путь удовлетворить уравнение
Шредингера в этих областях – предположить, что волновая функция здесь
тождественно равна нулю
I ( x )  0
(5)
III ( x)  0
и считать неопределенности   I ( x) и   III ( x) равными нулю.
Поскольку в дальнейшем это предположение к противоречию не приводит,
то оно является допустимым. Иными словами, в областях I и II частицы
нет.
Осталось найти волновую функцию в области II. Уравнение (2)
представляет собой хорошо известное дифференциальное уравнение
второго порядка, решение которого - гармоническая функция
II ( x)  A sin( kx)  B cos( kx)
где
k
(6)
2mE

8
А и В - константы интегрирования, найти которые нам предстоит дальше.
Кроме того, как только в процессе решения задачи будет определена
величина k, автоматически в соответствии с (6), это даст нам численное
значение полной энергии частицы Е.
Воспользовавшись первым из граничных условий - в точке х=0,
получаем
0  A sin( k  0)  B cos( k  0) ,
первое слагаемое равно нулю при любом значении А ( из-за синуса), а
второе равно нулю только если B  0 . Тогда
II ( x)  A sin( kx)
(7)
Второе граничное условие - в точке x=l дает
0  A sin( k  l )
(8)
или
k  l    n , где n-целые числа
(8а)
Следует отметить, что формально равенство (8) будет выполняться и при
А=0. Тогда волновая функция II тождественно равна нулю. А если учесть,
что и I иIII тоже равны нулю, то выходит, что во всем пространстве
частицы нет. Таким образом, вариант А=0 соответствует тривиальной
ситуации.
Возвращаясь к соотношению (8а), записываем
2mE
n
k

l
или, выражая Е,
En 
2n 2 2
2ml 2
(9)
Это есть искомое выражение для энергетического спектра частицы. Как
уже отмечалось (формула (8а)), n принимает значения целых чисел. Таким
образом, энергия частицы не непрерывна, а дискретна, то есть может
принимать ряд четко определенных значений, а квантовое число n
соответствует номеру энергетического уровня частицы.
Может ли число n быть равным нулю? Если да, тогда это означало
бы, что в области II находится частица, имеющая одновременно нулевую
потенциальную энергию (по условию) и равную нулю полную энергию
(согласно(9)), а,следовательно, нулевую кинетическую энергию. В то же
время k=0 (согласно(7)), означает, что и волновая функция частицы в
области II тоже тождественно равна нулю. Таким образом, очевидно,
значение n = 0 описывает тривиальную ситуацию - отсутствие объекта
исследования. Посему n = 0 включать в рассмотрение нет смысла, и будем
считать, что n=1,2,…
Волновая функция частицы равна
9
0, x  0

n

 ( x)   A sin(
x), n  1,2,3,...
l

0, x  0
0  x l
(10)
Остается найти нормировочный множитель А. Эта процедура является
простым математическим упражнением по вычислению интеграла.
Вспомнив, что волновая функция, будучи возведеной в квадрат, даст
плотность вероятности нахождения частицы в данной области, а после
2
интегрирования  (x) по всей области – полную вероятность нахождения
частицы во всей области, а также считая факт нахождения частицы хоть
где-нибудь внутри области достоверным, то есть таким, вероятность
которого равна единице, получаем равенства,
2

  ( x)
dx 
l
2


II
2
0


I
0

0
II
l
2
0

( x) dx  1
I
l
( x) dx  1   A 2 sin 2 (
0
n
l
2


( x) dx  0

2
2

( x) dx    II ( x) dx    III ( x) dx  1
0
l
2
l
III
( x) dx  0
l
x) dx
(11)
откуда
A
1
l
 sin
0
2
n
l

x dx
2
l
Окончательно, волновая функция, описывающая нахождение частицы в
бесконечно глубокой потенциальной яме,
 ( x) 
2
n
sin( x)
l
l
(12)
Задание 2.
Де Бройль предположил, что всякому микрообъекту можно
поставить в соответствие волну, длина которой обратно
пропорциональна импульсу частицы
  2 / p  2 / 2mE
(13)
(очевидно, что p2=2mE , Е – кинетическая энергия частицы). Если
рассматривать эту волну, заключенную в некоторой части
пространства, то логично сделать вывод, что стационарному (или
неизменному во времени ) состоянию должна соответствовать
10
стоячая волна. Такое состояние может реализоваться совсем не
всегда, а только когда на краях области (в нашем случае – на краях
ямы) будут узлы. И, следовательно, на ширине ямы будет
укладываться целое число полуволн де Бройля.
l  n / 2
n  1,2,3
(14)
Подставив (13) в (14) легко получаем выражение для энергии
частицы Е
En 
 2 2
2ml
2
 n2
n  1,2,3...
(15)
Последнее выражение (15) тождественно результату (9),
полученному в задании 1.
Вычислим несколько значений энергии на соседних
энергетических уровнях и оценим, как изменяется “расстояние”
между ними по мере роста номера уровня n.
Таб.1
E=|En - En-1|
En
n
 
2
2ml 2
 2 2
4
2ml 2
 2 2
9
2ml 2
 2 2
 16
2ml 2
1
2
3
4
= E/En
2
3
 2 2
2ml 2
 2 2
5
2ml 2
 2 2
7
2ml 2
3/4
5/9
7/16
Из приведенной выше таблицы видно, что абсолютное “расстояние“
между уровнями E увеличивается с ростом номера уровня, в то
время, как относительное расстояние  уменьшается, стремясь к
нулю как 1/n согласно формуле:
Lim n= Lim n(n2-(n-1)2)/n21/n0
Задание 3
Плотность уровней
dn
легко получим, взяв дифференциал от
dE
обеих частей равенства (15)
dE 
 2 2
2ml 2
 2  n  dn
11
и воспользовавшись еще раз соотношением (15) только в немного
измененной форме
 2 2
2ml
2

E
n2
Окончательный результат выглядит так
dn
n

dE 2 E
(16)
Для электрона при Е=1,0 эв и l=1,0 см
dn
l

dE 
m
10 2 0.511  10 6

 0.12  1016 (эв-1)
15
2 E 3.14  0.659  10
2
Задание 4
Соотношение неопределенностей Гейзенберга служит для
описания системы, которая находятся в пограничном состоянии
между макро- и микромиром и проявляет себя и как классический и
как квантовый объект. Поэтому понятие «сила», которое вводится в
классической механике и уже не имеющее смысла в квантовой
физике, еще может быть использовано для проведения оценочных
расчетов.
Если частицу рассматривать как классический объект, то сила
давления, как и любая сила, согласно второму закону Ньютона,
может быть вычислена как отношение приращения импульса
частицы к промежутку времени, за который это приращение
произошло,
F
p x
t
Приращение
импульса
имеет
место
вследствие
его
неопределенности, поскольку наша частица проявляет и волновые
свойства. Тогда для оценки силы давления можно использовать
приближенные (или оценочные) значения для px и t, полученные
из оотношения неопределенностей Гейзенберга, которое связывает
между собой кординату с проекцией импульса, и время с энергией
x  p  
E  t  
.
12

x


2ml 2
t 

 2
E E1
 
p x 
(17)
Здесь для простоты считаем, что частица находится в основном
энергетическом состоянии. Тогда, разделив почленно равенства (17)
одно на другое и пренебрегая несущественны множителем 2/2
получим
F
2
ml 3
(18)
Еще раз напомним, что соотношение Гейзенберга само по себе
оценочное, поэтому точность ответа в пределах одного порядка
является приемлемой.
Задание 5
Вероятность того, что частица находится где-то внутри
области 0хl, равна единице, поскольку мы исходим из того, что,
во-первых, частица существует и, во-вторых, «живет» внутри этой
области. Математическим выражением этого факта является условие
нормировки
2
 2  n 
P    dv   
sin  x  dx  1
l
 l 
v
0
l
2
Чтобы определить вероятность нахождения частицы в какойто части этой области достаточно в приведенном выше интеграле
указать координаты этой области
2
l
3
P1 / 3
2
 2  n 
1
 
sin  x  dx  1 / 3 
2
 l 
1  l
l
 4 
 2 
sin  3 n   sin  3 n 



 
3
Дальнейшие вычисления зависят от того, на каком энергетическом
уровне находится частица. По условию задания она находится в
основном состоянии, то есть n=1.
P1 / 3  1 / 3 
1
2
 4 
3
 2 
sin  3    sin  3    1 / 3  2  0.61



 
13
Таким образом, вероятность того, что частица находится в
средней трети области, равна 0,61.
Задание 6
Связь между операторами и измеряемыми физическими
величинами в квантовой механике устанавливается с помощью
формулы для среднего значения величины m в состоянии,
описываемом волновой функцией 

M    * MdV
(19)
V
Это соотношение является почти тривиальным, если принять во
внимание следующие рассуждения.
Пусть путем экспериментального измерения
некоторой
физической величины мы убеждаемся, что получаем одно и то же
численное значение (разумеется, в пределах погрешности
эксперимента), сколько бы измерений мы ни проводили. Это значит,
что данная величина имеет точное значение, которое может быть
получено как среднее по всем измерениям. С точки зрения
квантовой механики, существование точного значения физической
величины означает, что для оператора этой физической величины
волновая функция системы, находящейся в данном состоянии,
является собственной, а значит,

M  M  
(20)
где М – собственное число. Напомним, что собственное число по
своей сути как раз и является точным значением измеряемой
величины. Теперь умножим выражение (20) слева на ψ* :

* M  * M  
Поскольку в правой части равенства содержится число М –
множитель, то его запросто можно поставить на первое место, после
чего проинтегрировать равенство по всему фазовому объему

*
*

M
 dV  M   dV  M  M
Таким образом, точное значение физической величины, если таковое
существует, с одной стороны, является собственным значением
соответствующего оператора, а с другой – средним значением
14
измеряемой величины, разумеется, при достаточно большом
количестве измерений.
Вернемся к частице в потенциальной яме и определим среднее
значение ее координаты и проекции импульса на ось х.
Воспользовавшись уравнением (19), явным видом волновой
функции частицы (12), которая находится в потенциальном ящике в
основном состоянии ( n = 1 ), а также тем фактом, что оператор
координаты это умножение на координату, запишем
l

0
2
  2

2

l
sin( x) x (
sin( x)) dx   x sin 2 ( x)dx 
l
l
l
l
l
l
2
0
l
Итак, среднее значение координаты частицы в яме равно l/2.
Среднее значение импульса частицы равно нулю, что видно из
l

0
2
 
2

2
 
 2

sin( x) p x (
sin( x)) dx  
sin( x)  i (
sin( x)) dx  0
l
l
l
l
l
l 
x  l
l
0
l
Приведенные результаты оказываются почти очевидными,
если подойти к данному вопросу с классической точки зрения и
принять во внимание, что, внутри ямы координата частицы может
быть любой, а среднее значение координаты, естесственно, это
середина ямы, то есть точка l/2. Проекция же скорости (и импульса)
– это величина векторная колеблется вокруг положения равновесия,
как у маятника, изменяясь и по модулю, и по направлению. Поэтому
среднее значение проекции скорости (и импульса) равно нулю.
Следует отметить, что интерпретировать результаты
квантовомеханических расчетов с помощью классических аналогий
можно в очень редких случаях, поскольку квантовая механика, как
теория описания физических процессов вступает в свои права тогда,
когда классическая теория себя исчерпывает и не в состоянии
адекватно объяснять явления. Это происходит, например, на границе
макро и микромира. Что касается данной конкретной задачи, то ее
можно в какой-то мере рассматривать как пограничную, которая еще
как-то подчиняется классическим законам, но в то же время
квантование (дискретность) энергии не может быть объяснено
только классической теорией. Очевидно, здесь требуется
принципиально новый, квантовомеханический подход.
15
Частица в сферически-симметричной
потенциальной яме.
Потенциал, как известно, это универсальная функция координат,
характеризующая потенциальную энергию частицы в любой точке
пространства. Обычно используют прямоугольную систему
координат, тогда потенциал есть функция декартовых координат
U(x,y,z). В предыдущем разделе мы рассмотрели подробно частный
случай, когда потенциал зависит не от всех трех декартовых
координат, а лишь от одной; таким образом, «жизнь» частицы,
процессы – все происходит вдоль одного-единственного,
выделенного направления. Вообще говоря, в качестве системы
координат можно использовать любую, вариантов можно
предложить более десятка, например, хорошо известную
цилиндрическую или сферическую. Однако, бывают такие задачи, в
которых удобно использовать совершенно конкретную систему
координат, это приводит к резкому упрощению расчетов и анализа
результатов. В первую очередь выбор системы координат диктуется
симметрией потенциала. В этом случае говорят, что симметрия
задачи совпадает с его симметрией.
Чаще всего (возможно, даже чаще, чем прямоугольная) в задачах
встречается сферическая симметрия. Продиктовано это тем, что в
квантовой механике существует класс задач, посвященных
описанию поведения элементарных частиц (к примеру, электронов),
входящих в состав атомов, а, значит, находящихся в кулоновском
поле ядра. А ведь хорошо известно, что потенциальная энергия
частицы в кулоновском поле зависит не от абсолютного положения
частицы в пространстве, а от того, насколько далеко она находится
от источника поля. Поэтому говорят, что такая задача имеет
сферическую симметрию. Другими словами, потенциальная энергия
частицы имеет одно и то же значение для всех точек,
принадлежащих сферической поверхности радиуса r : U=U(r).
Что же представляет собой бесконечно глубокая сферическая
яма? В действительности она является почти полностью
придуманным, идеализированным объектом, как, впрочем, и
одномерная потенциальная яма с бесконечными стенками. Проще
всего представить ее себе, если вообразить некую точку,
окруженную на некотором расстоянии R0 абсолютно непроницаемой
сферической поверхностью. Тогда если внутри этой поверхности
находится частица, ей очень не повезло – она никогда не сможет
вырваться за пределы сферы. Тот факт, что для того, чтобы покинуть
пространство в пределах сферы R0 она должна получить бесконечно
16
большую энергию, и означает, что она находится в бесконечно
глубокой потенциальной яме.
Потенциал внутри ямы может быть не обязательно кулоновским,
в зависимости от условия задачи он может иметь любой другой
аналитический вид. Очевидно, это повлияет на состояние частицы
внутри ямы. Теперь для того, чтобы начать изучать поведение
частицы в такой потенциальной яме рассмотрим конкретное условие
задачи.
Постановка задачи
Частица массы m находится в сферически-симметричной
потенциальной яме, в которой потенциал определяется, как
0, r  R0
U (r )  
, r  R 0
Для случая, когда волновая функция зависит только от r
найти:
1.
Нормированные волновые функции частицы
2.
Возможные значения энергии частицы
3.
Для основного состояния частицы наиболее вероятное
значение расстояния rвер , а также вероятность нахождения частицы
в области r < rвер
Решение
Задания 1 и 2
Возможные значения
операторное уравнение
энергии
частицы
находим,
решая
Hˆ   E
2
ˆ
H  U (r )
2m
Поскольку потенциал имеет сферическую симметрию, то оператор
Лапласа разумно будет записать также в сферических координатах

2 2 
1


1
2


(sin

)


r 2 r r r 2 sin  
r 2 sin 2   2
17
Кроме того, волновая функция частицы в потенциале, который
зависит только от одной координаты r , будет также зависеть только
от этой одной координа-ты. Тогда дифференцирование волновой
функции по остальным переменным будет давать нуль, а значит, в
выражении для оператора Лапласа можно оста-вить только те
слагаемые, которые содержат производные только по r
r 
2 2 

r 2 r r
Следуя форме потенциала, разделим пространство на две области –
первая область внутри сферы радиуса R0, а вторая - за пределами
этой сферы. Для каждой из этих областей составим и решим
уравнение Шредингера, а потом «сошьем» решения на границе
областей. Итак, уравнения
 2  2 (r )

(

2m r 2
 2  2 (r )

(

2m r 2
2  (r )
)  0  (r )  E (r )
r r
2  (r )
)    (r )  E (r )
r r
в чем-то похожи на те, что мы детально исследовали в разделе 1.
Второе уравнение содержит бесконечно большое слагаемое, поэтому
его решение – волновая функция, тождественно равная нулю.
Первое уравнение отличается от того, что мы анализировали в
разделе 1, тем, что содержит в операторе Лапласа дополнительное
слагаемое. Такое уравнение решается стандартным методом
подстановки
 (r ) 
 (r )
r
,
где функция (r) является неизвестной. После подстановки волновой
функции в таком виде в уравнение Шредингера, выполнения всех
необходимых операций дифференцирования и всех возможных
сокращений приходим к уравнению для (r)

 2  2  (r )
 E (r )
2m r 2
которое является хорошо изученным и до боли знакомым
уравнением гармонических колебаний. Его решение
18
 (r )  A sin( kr )  B cos(kr )
k
2mE

Возвращаемся к подстановке ( ). Тогда волновая функция принимает вид
 (r ) 
 (r )
r

A sin( kr)  B cos( kr)
r
Для нахождения констант интегрирования А и В воспользуемся
следующими соображениями. Вследствие деления на r , в волновой
функции в точке r =0 имеется сингулярность, или, проще говоря,
возрастание к бесконечно большой величине. Волновая функция,
однако, не имеет права быть неограниченной, поэтому, чтобы
избавиться от этой неприятности нам остается только потребовать,
чтобы числитель также стремился к нулю. Тогда неопределенность
типа
0
0
легко раскрывается (хотябы по правилу Лопиталя) и
приводит к конечной величине (r) . Тогда условие
 (0)  A sin( k  0)  B cos( k  0)  0
дает В=0 и волновая функция принимает вид
A sin( kr )
r
 (r ) 
Далее воспользуемся граничным условием. Волновая функция на
границе областей должна совпадать. Так как в области r>R0 она
равна нулю, то и на границе ( то есть в точках с r=R0) тоже должна
принимать нулевое значение
A sin( kR0 )
 ( R0 ) 
0
R0
Автоматически это дает
kR0  n
En 
 2 2
2
0
R
2m
n2
19
Таким образом, энергетический спектр частицы нами уже
получен.
Последнюю неизвестную константу интегрирования А находим
из условия нормировки, которое состоит в том, что волновая
функция, после интегрирования в квадрате по объему dV=4 r2 d r
должна дать единицу
  (r ) dV 
2
V
тогда
1
A
2R0
приобретает вид
R0

0
A 2 sin 2 (kr)
4r 2 dr  1
2
r
, и искомое выражение для волновой функции
 n (r ) 
1
2R0
sin(
n
R0
r
r)
Задание 3.
Для нахождения наиболее вероятного значения rвер необходимо
исследовать на
предмет
наличия
экстремума
функцию
1 sin 2 (kr)
 (r ) 
4r 2 , которая является плотностью вероятности
2
2R0
r
нахождения частицы в сферическом слое толщины
Дифференцирование  (r ) по r дает нам условие экстремума
2k
2 sin( kr) cos( kr)  0
R0
dr.
или, с учетом соотношения ( ), а также
того, что, по условию, частица находится в основном
энергетическом состоянии (т.е. n=1 ), определяем значения rextr,
соответствующие экстремумам функции  (r ) :
1. sin(

R0
rextr )  0
rextr1 = 0,
rextr2 = R0
соответствуют минимумам функции  (r ) , в чем легко убедиться
прямой подстановкой.
2. cos(

R0
rextr )  0
20
rextr3 = R0/2 = rвер соответствует максимуму функции плотности
вероятности и как раз является искомым наиболее вероятным
значением.
Поиск вероятности нахождения частицы в области r < rвер сводится к
вычислению интеграла

sin 2 ( r )
rextr 3
R0 / 2
R0
1
2 2 
2 R0 1
2
0 2R0 r 2 4r dr  0 R0 sin ( R0 r )dr  R0 4  2

2.0
1.6
1.2
0.8
0.4
r
0.0
0.0
0.2
0.4
0.6
0.8
1.0
Рис.3
На рис.3 представлена зависимость  (r ) для R0=1, откуда очевидно, что
максимум приходится на точку rвер = 1 / 2.
Туннельный эффект
Явление туннелирования относится к группе самых ярких
явлений, экспериментально подтверждающих теоретические
положения квантовой механики. На бытовом уровне невозможно
понять, как частица умудряется пробраться через некоторое
пространство, если ее энергии для этого недостаточно? Яркий
пример тому – естесственная радиоактивность ядер : взаимодействие
между нуклонами в ядре настолько велико, что непонятно, по какой
причине время от времени из ядер, скажем, стронция вылетают
альфа-частицы. Следует сразу подчеркнуть, что вооружившись
всеми теми знаниями, которые дает квантовая механика, мы все
равно этого не узнаем, ибо вопрос «почему?» - это вопрос не по
21
адресу. Квантовая механика в состоянии ответить только на вопрос
«какова вероятность того, что данное событие при
определенных условиях будет иметь место», не вникая в причины
этого явления.
Таким образом, вопрос, а какова вероятность прохождения
частицы через данную область пространства, для которой известен
аналитический вид потенциала, является правомочным и ответ на
него мы получим, вычислив коэффициент туннелирования.
Постановка задачи
Найти для электрона, имеющего энергию Е, вероятность D
прохождения потенциального барьера шириной l и высотой U0, если
барьер имеет форму как на графиках, представленных на рис.4 а), б),
в). Уточнение: на графике в) представлена квадратичная парабола.
Для этого
1. Написать аналитическое выражение для потенциалов U(x) в
области 0xl .
2. Для каждого варианта определить координаты x1, x2,
соответствующие координатам точек входа и выхода из барьера
соответственно.
3. Вычислить коэффициент прохождения, воспользовавшись
формулой
 2 x2

D  exp   2m(U ( x)  E )dx 
  x1

U0
Е
а)
б)
0
l
Е
U0
Е
U0
в)
0
l
-l/2
Рис.4
Решение
Задания будем выполнять параллельно для всех трех вариантов .
Задание 1
а) Очевидно, U ( x)  U 0  const
0
l/2
22
б) Как видно из графика, потенциал линейно зависит от координаты
U ( x) 
U0
x
l
в) Согласно условию, потенциал квадратично зависит от
координаты. Поскольку «рожки» параболы направлены вниз,
коэффициент при старшей степени x отрицаиелен, а всилу
симметричного расположения параболы относительно оси ординат,
коэффициент при линейном члене равен нулю. Воспользовавшись
тем фактом, что парабола проходит через точки (0, U0) и (l/2, 0),
находим аналитический вид потенциала
U ( x)  U 0  U 0
4 2
x
l2
Задание 2
Точки x1 и x2 это точки пересечения графика U(x) и прямой U=Е
а) Этот случай тривиальный, поэтому сразу запишем
x1  0
x2  l
б) Координату входа в барьер находим из уравнения U 0
x1
 E, а
l
координату выхода – по графику
x1  l
E
U0
x2  l
в) Симметричные относительно оси ординат точки входа-выхода из
барьера находим из уравнения E  U 0  U 0
x1  
их значения
Задание 3
x2 
l
E
1
2
U0
l
E
1
2
U0
4 2
x1, 2
l2
23
а)
 2 x2

 2l

D  exp   2m(U ( x)  E )dx   exp   2m(U 0  E )dx  
  x1

 0

 2

exp  l 2m(U 0  E ) 
 

б)


 2 x2

 2 l

x
D  exp   2m(U ( x)  E ) dx   exp   2m(U 0  E ) dx  
l
  x1

 lE

U
0


 8 l
exp 
 3 U 0

2m(U 0  E ) 3 

в)
 2
D  exp 
 
x2

x1
l
E


1
2
U0

 2

4 2
2m(U ( x)  E ) dx   exp 
2
m
(
U

U
x

E
)
dx

0
0 2


l

l
E


 1


2
U0
  l 2m

exp 
(U 0  E )
  2 U0

Примечание: взятие последнего интеграла производится при помощи
стандартной тригонометрической замены
l
E
l
E
x
1
sin 
dx 
1
cos d
2
U0
2
U0
С целью проверки полученных знаний и навыков студентам
предлагается приведенная ниже самостоятельная работа
24
Самостоятельная расчетная работа «Определение
коэффициента туннелирования электрона через
потенциальный барьер ».
ЗАДАНИЕ
Электрон, имея кинетическую энергию Е, налетает на одномерный
потенциальный барьер, вид которого приведен в таблице вариантов (таб.2).
1-й уровень сложности – потенциальный барьер состоит из 2-х
частей, каждая из которых представляет собой линейную функцию
(графики 1-6).
2-й уровень сложности – потенциальный барьер состоит из 2-х
частей - комбинация линейной и квадратичной функции (графики 718).
3-й уровень сложности – потенциальный барьер состоит либо
из 3-х частей – комбинации линейных и квадратичных функций
(графики 19-30), либо из комбинации 2-х разных квадратичных
функций (графики 31-34).
1.
2.
3.
4.
5.
Вычислить коэффициент туннелирования D. Для этого
Выбрать вариант в соответствии с уровнем сложности, определенным
преподавателем.
Записать аналитическое выражение для потенциала U(x).
Определить координаты точек вхождения и выхода из барьера x1 и x2.
Получить явное выражение для коэффициента туннелирования D.
Указание: для облегчения процедуры взятия интеграла ввести новую
безразмерную переменную t=x/l и безразмерный параметр =E/U0.
Вычислить вероятность туннелирования α-частицы (массу α-частицы
взять в справочнике) для случая, когда ее кинетическая энергия
составляет Е = 0.6 U0, а ширина барьера – порядка радиуса ядра атома
урана.
Таблица функций U(x), значений x1 и x2 , значений интегралов, а
также коэффициентов D для самопроверки приведена в таб.3.
25
Таб.2 Варианты потенциальных барьеров
1
Uo/E
2
Uo/E
3
Uo/E
4
Uo/E
5
Uo/E
6
Uo/E
1.0
1.0
1.0
1.0
1.0
1.0
0.8
0.8
0.8
0.8
0.8
0.8
0.8
0.6
0.6
0.6
0.6
0.6
0.6
0.4
0.4
0.4
0.4
0.4
0.2
0.2
0.2
0.2
0.2
0.6
0.4
0.4
0.2
0.2
0.0
0.0
0.0
0.2
0.4
0.6
0.8
1.0
Uo/E
t=x/l
8
1.0
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.0
0.0
9
Uo/E
0.2
0.4
0.6
0.8
1.0
0.8
0.8
t=x/l
0.0
10
Uo/E
1.0
1.0
0.0
0.4
0.6
0.8
0.6
0.6
1.0
11
Uo/E
0.8
0.6
0.2
t=x/l
0.0
0.0
0.2
0.4
0.6
0.8
1.0
1.0
0.8
0.8
0.6
0.6
0.2
0.4
0.6
0.8
1.0
t=x/l
0.0
0.0
12
Uo/E
1.0
0.0
t=x/l
13
Uo/E
0.2
0.4
0.2
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.4
0.2
0.2
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.0
15
Uo/E
16
Uo/E
1.0
0.8
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0.0
0.0
0.2
0.4
0.6
0.8
1.0
22
0.8
1.0
t=x/l
17
0.2
0.4
0.6
0.8
1.0
23
Uo/E
1.0
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0.2
0.4
0.6
0.8
1.0
t=x/l
24
Uo/E
0.4
0.4
0.4
0.2
0.2
0.2
0.0
0.0
0.0
29
Uo/E
1.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
30
Uo/E
1.0
0.4
0.6
0.8
1.0
0.4
0.6
0.8
1.0
t=x/l
19
Uo/E
0.0
0.0
0.2
0.2
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
31
Uo/E
1.0
0.0
0.2
0.4
0.6
0.8
1.0
0.0
0.0
26
Uo/E
1.0
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0.0
0.2
0.4
0.6
0.8
Uo/E
1.0
t=x/l
0.0
32
1.0
1.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
33
Uo/E
0.6
0.4
0.2
0.2
0.2
0.2
0.2
0.2
0.0
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.8
1.0
t=x/l
27
21
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
28
Uo/E
1.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.2
0.2
0.4
0.6
0.8
1.0
t=x/l
34
Uo/E
0.4
t=x/l
0.6
1.0
0.6
1.0
t=x/l
0.4
0.0
0.4
0.8
0.4
0.0
0.6
0.6
1.0
0.6
0.4
0.4
0.8
0.8
0.6
0.2
0.2
Uo/E
1.0
0.4
0.0
0.6
Uo/E
0.0
0.6
0.0
0.4
0.2
t=x/l
0.4
t=x/l
0.2
0.4
0.6
1.0
20
14
0.0
0.0
t=x/l
0.6
0.4
0.8
0.8
t=x/l
1.0
0.8
0.8
0.6
1.0
0.8
0.6
0.8
0.4
0.8
0.6
0.8
0.8
0.2
0.6
0.4
1.0
0.8
0.0
0.4
Uo/E
0.8
0.0
0.2
1.0
0.4
t=x/l
25
Uo/E
0.2
t=x/l
0.2
1.0
0.4
1.0
0.2
0.8
0.0
0.6
0.8
0.0
1.0
0.0
0.6
0.6
18
0.0
0.0
0.6
0.4
t=x/l
0.2
0.0
0.8
0.2
1.0
0.4
0.6
0.0
0.8
0.6
0.8
0.0
0.6
Uo/E
0.8
0.8
0.4
0.6
1.0
1.0
0.2
0.8
0.2
t=x/l
0.0
1.0
0.0
Uo/E
0.6
0.4
0.0
1.0
0.4
0.6
0.0
t=x/l
0.2
Uo/E
1.0
0.2
0.0
0.0
1.0
0.4
0.2
Uo/E
1.0
0.4
0.4
7
Uo/E
1.0
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
0.0
0.0
0.2
0.4
0.6
0.8
1.0
t=x/l
26
Таб.3 Таблица ответов для самопроверки
Аналитическое
выражение для
потенциала
Примерный вид потенциала
(E / Uo =0.6)
Аналитическое выражение для
коэффиуиента туннелирования D
(t =x/l
(t =x/l )
Пределы
интегрирования
Коэффициент туннелирования D
 = E / Uo)
Uo/E
1.0
0.8
I1
 4t (t  1)
0.6
t2
0.4
0.2
x/l
0.0
0.0 0.2 0.4 0.6 0.8 1.0
e
2
 2 mU l  4t ( t 1)  dt

t1
Uo/E
1.0
t1 
1
(1  1   )
2
e
t2 =1 / 2

2
(1 ) 
2 mU 0 l

4 2

2
(1 ) 
2 mU 0 l

4 2
t1 = 1 / 2
0.8
I2
 4t (t  1)
0.6
0.4
0.2
x/l
0.0
0.0 0.2 0.4 0.6 0.8 1.0
0.0
0.0
0.5
e
t2
2

2 mU l   4t ( t 1)  dtt

t1
1.0
x/l
4t 2
0.4
0.6
e
0.8
1.0
Uo/E
t2
2
 2 mU l  4t 2  dt

t1
Uo/E
1.0
I4
4(t  1) 2
0.8
0.6
e
0.4
0.2
0.0
0.0
0.5
t2
2

2 mU l  4 ( t 1) 2  dtt

t1
0.8
2t
0.6
0.4
0.2
0.0
0.0
0.5
t1=1 / 2
1
t2  1 

2
x/l
1.0
Uo/E
1.0
I5
1

2
t2 =1 / 2
x/l
1.0
e
t1 
0.2
I3
1
t 2  (1  1   )
2
e

t2
2
2 mU l  2t  dt

t1
1
t1  
2
1
t2 
2
e

2
1 
1 1 
2 mU 0 l  1  ln


4 
 

2
1 
1 1 
2 mU 0 l  1  ln


4 
 
e
27
e
Аналитическое
выражение (t
Примерный вид потенциала
(E / Uo =0.6)
Аналитическое выражение для
коэффиуиента туннелирования D
=x/l )
2(1  t )
(t =x/l
 = E / Uo)
t1 
Uo/E
1.0
0.8
0.6
e
0.4
I6
0.2
0.0
0.0
Пределы
интегрирования
0.5

t2
2
2 mU l  2 ( t 1) dt

t1

2
2 mU 0

x/l
1.0
1
3
Коэффиуиент туннелирования D
1
2
1
t2  1  
2
3
l (1 ) 2
3
e
2
1

2 mU 0 l (1 ) 2

3
Uo/E
1.0
t1  0
0.9
I7
0.8
2(1   )t  
t2
2
 2 mU l  2 ( t 1) dt

t1
1
t2 
2
t2
2
 2 mU l  2( t 1)(1 ) dt

t1
1
t1 
2
t2  1
0.7
0.6
0.5
0.0
x/l
0.4
0.8
e
e

Uo/E
1.0
2(1   )(1  t )
I8
0.9
0.8
e
0.7
0.6
0.5
0.0
0.8
0.8
I9
0.6
e
0.4
0.2
0.0
0.0
2
1
2 mU 0 l 1

3
x/l
0.4
Uo/E
1.0
3t
e

2
1
2 mU 0 l 1

3
0.5
x/l
1.0

t2
2
2 mU l  3t  dt

t1
t1 

3
1
t2 
3
3
e
2
2

2 mU 0 l (1 ) 2

9
28
Uo/E
1.0
(1  t )
0.8
0.6
I 10
e
0.4
0.2
0.0
0.0
Аналитическое
выражение (t
t2
2
 2 mU l  3(1t )  dt

t1
2
t1 
3
3
t2  1 

3
e
2
2

2 mU 0 l (1 ) 2

9
x/l
1.0
0.5
Примерный вид потенциала
(E / Uo =0.6)
Аналитическое выражение для
коэффиуиента туннелирования D
=x/l )
(t =x/l
Пределы
интегрирования
 = E / Uo)
Коэффиуиент туннелирования D
Uo/E
1.0
t1  0
0.9
3(1   )t  
I 11
0.8
0.7
0.6
e
x/l
0.5
0.000
0.667

t2
2
2 mU l  3t (1 )dt

t1
t2 
2
3
t2  1
Uo/E
1.0
I 12 3(1   )(1  t )  
0.8
e
Uo/E
0.7
1.00.6
0.80.5
e
t1 
t2
0.9
1
3
2
 2 mU l  3(1 )( t 1) dt

t1
e

2
2
2 mU 0 l 1

9

2
2
2 mU 0 l 1

9

2
(1 ) 
2 mU 0 l

6 2
x/l
0.000
0.667
0.6
0.4
0.2
 9t (t  2 / 3)  
0.0
0.0
x/l
0.1
0.2
0.3
0.4
I 13
e

t2
2
2 mU l  9t ( t 2 / 3) dt

t1
1 1
3
1
t2 
3
t1 
e
Uo/E
1.0
 9(t  1)(t  2 / 3)  
I 14
0.8
0.6
e
0.4
0.2
x/l
0.0
0.6
Uo/E
1.0
0.8
0.6
0.8
1.0

t2
2
2
2 mU l  9 ( t 1)( t 2 / 3) dt t1 

3
t1
t2 
2
(1  1   )
3
e

2
(1 ) 
2 mU 0 l

6 2
29
I 15
9t2
Аналитическое
выражение (t
e
Примерный вид потенциала
(E / Uo =0.6)

2
2 mU l  9t 2  dt

t1
Аналитическое выражение для
коэффиуиента туннелирования D
=x/l )
1.0
9(t-1)2
e
0.8
0.6

t1 

3
1
t2 
3
e

Пределы
интегрирования
t2
2
2 mU l  9 ( t 1)2  dt

t1
t1 
2
3

t2  1 
3
0.4
2
1
1 1

2 mU 0 l  1  ln

6


e

2
1
1 1

2 mU 0 l  1  ln

6


0.2
x/l
0.0
0.6
0.8
1.0
Uo/E
1.0

0.8
I 17
0.7
dt
t1
0.6
0.5
0.0
t1  0
t2
0.9
1
0.5
1

0.8
0.7
t1
0.6
0.5
x/l
1.0
e
t1  0
t2
0.9
0.5
0.0
1
2
x/l
1.0
Uo/E
1.0
I 18
t2 
dt
1
t2 
3





Коэффиуиент туннелирования D
 = E / Uo)
(t =x/l
Uo/E
I 16
t2
e

2
1
2 mU 0 l 1

2

2
1
2 mU 0 l 1

3





Download