Московский институт электронной техники
Кафедра КФН
Лабораторная работа 3
Туннелирование волнового пакета через прямоугольный барьер.
Общие указания.
1. В начале надо познакомиться с программой Square_ Barrier, которая моделирует туннелирование волнового пакета через прямоугольный барьер. Файлы программы находятся в папке Square_ Barrier. Программа запускается из файла main.m. Параметры задачи определяются в файле Data_Packet. В этом файле студент должен изменять параметры в соответствии с задачей своего варианта. Параметры, которые можно изменять в этом файле имеют следующий физический смысл.
1)
- средняя энергия электрона (эВ).
2)
- полуширина волнового пакета (Å).
3)
- координата центра волнового пакета
(Å).
4)
-
высота потенциального барьера (эВ).
5)
- координаты, определяющие положение
потенциального барьера (Å).
6) xBegin, xEnd – координаты, определяющие область пространства, в котором двигается волновой пакет.
7) dt – временной шаг.
2. Далее надо запустить программу Square_ Barrier из файла main.m и поэкспериментировать с управляющими элементами. Элементы управления имеют следующий смысл.
1) Start – запуск программы с начального положения.
2) Stop – остановка выполнения программы.
3) Move – возобновление работы программы с места остановки.
4) Close – грамотное завершение работы программы.
В верхней части графического окна программы выводятся изменяющиеся параметры программы, которые имеют следующий физический смысл.
1) t – текущее время, измеряется в фемтосекундах (фс).
2)
- вероятность нахождения электрона
перед барьером, которая вычисляется с
помощью следующего интеграла.
. (1)
3)
- вероятность нахождения электрона
внутри барьера, которая вычисляется с
помощью следующего интеграла.
. (2)
4)
- вероятность нахождения электрона
позади барьера, которая вычисляется с
помощью следующего интеграла.
. (3)
На Рис. 01 показана начальная стадия движения волнового пакета.
Рис.01.
Начальная стадия движения волнового пакета.
Программа Square_
Barrier моделирует
движение электрона в состоянии с волновой
функцией
имеющей
вид волнового пакета (смотри лабораторную
работу 2, приложение А). Кривые на экране
описывают распределение плотности
вероятности
в последовательные моменты времени.
Для простоты будем называть волновым
пакетом распределение плотности
вероятности
.
Синяя кривая на Рис. 01 показывает
положение пакета в начальный момент
времени. В качестве начального значения
пакета берется гауссов волновой пакет
(смотри лабораторную работу 2, приложение
А, формула (А.9)). На Рис. 01 видно, что
среднее положение пакета выбрано равным
,
а полуширина пакета равна
.
Красная кривая показывает положение
пакета в момент времени
.
В этот момент времени волновой пакет
полностью находится слева от потенциального
барьера. Поэтому вероятность найти
электрон слева от потенциального барьера
равна
.
Вероятность нахождения электрона внутри
барьера и справа от барьера
практически равна нулю.
На Рис. 01 прямоугольный потенциальный
барьер изображен синими толстыми
линиями. Показана также высота
потенциального барьера в электрон-вольтах
.
Светло-сиреневая линия показывает
среднюю энергию электрона
.
Две темно-сиреневые линия показывают
границы энергетического спектра
волнового пакета. Соответствующие
энергии
определяются следующими формулами.
(4)
Здесь
есть среднеквадратичное отклонение
волнового числа от среднего значения
,
и определяется следующей формулой.
(5)
Из формул (4), (5) видно, что чем меньше
полуширина волнового пакета
,
тем больше ширина энергетического
спектра
у этого волнового пакета.
На Рис. 02 показан момент взаимодействия волнового пакета с потенциальным барьером.
Рис.02.
Момент взаимодействия пакета с потенциальным барьером.
Красная кривая показывает положение
пакета в момент времени
.
На Рис. 02 показана часть пакета, которая
отразилась от барьера, и двигается
налево, удаляясь от барьера. Величина
этой части пакета определяет вероятность
электрону быть отраженным от барьера.
В данном примере эта вероятность равна
.
На Рис. 02 показана также часть пакета,
которая прошла через барьер, туннелировала
через барьер, и двигается направо,
удаляясь от барьера. Величина этой части
пакета определяет вероятность
туннелирования электрона через барьер.
В данном примере эта вероятность равна
.
Наконец Рис. 02 показана часть пакета,
которая локализована внутри барьера.
Величина этой части пакета определяет
вероятность нахождения электрона внутри
барьера. В данном примере эта вероятность
равна
.
3. Коэффициент отражения от барьера r (Reflection), и коэффициент прохождения через барьер t (Transmission), в задаче о туннелировании волнового пакета, определяются следующими формулами.
(6)
Заметим, что формулы (6) применяются для
достаточно больших моментов времени.
Нужно чтобы произошло взаимодействие
волнового пакета с потенциальным
барьером. Части волнового пакета
разбежались от барьера, а вероятность
нахождения электрона внутри барьера
стала равной нулю
.
Строго это выполняется в пределе
.
Реально будем выбирать такие времена,
когда вероятность
будет близка к нулю, например
.
На Рис. 03 показан момент, при котором
можно определить коэффициенты
.
Рис. 03.
Момент для определения коэффициентов
.
Отметим, что на Рис. 03 увеличены размеры
пространства, в котором движется пакет
по сравнению с Рис. 02. Это сделано, для
того чтобы в рабочем пространстве
уместились отраженная и прошедшая части
волнового пакета. В этом случае можно
считать, что волновой пакет равен нулю
за пределами рабочего пространства.
Тогда вероятности
будут совпадать с коэффициентами
отражения и прохождения.
(7)
В данном примере получаем следующие коэффициенты отражения и прохождения.
Сумма этих коэффициентов отличается
от единицы
.
Это связано с тем, что вероятность
нахождения электрона в барьере равна
примерно
.
На следующем Рис. 04 выбрано еще большее время движения пакета и большее рабочее пространство.
Рис. 04.
Момент для определения коэффициентов
.
В данном примере получаем следующие коэффициенты отражения и прохождения.
Сумма этих коэффициентов отличается
от единицы
на одну сотую. Таким образом, для
вычисления коэффициентов отражения и
прохождения с заданной точностью нужно
подобрать достаточно большое рабочее
пространство. С другой стороны, слишком
большие пространственные размеры
замедляют работу программы, и потребуется
много времени для выполнения заданий
лабораторной работы.
4. Заметим, что можно ускорить движение пакета, увеличивая временной шаг dt. Однако слишком большой шаг делает движение пакета не плавным, а рывками.
В программе Square_ Barrier предусмотрен флаг flagStop, который позволяет менять режим работы программы. Если флаг равен нулю flagStop = 0, то программа может работать в течение любого времени и части волнового пакета могут уйти за пределы рабочего пространства. Если же флаг равен единице flagStop = 1, то прекращается движение пакета в тот момент, когда части волнового пакета, достигают границ рабочего пространства.
Рекомендуется выбрать значение флага flagStop = 1.
|
Правила работы с программой Square_ Barrier. Программа
Square_ Barrier
используется, в данной лабораторной
работе, для нахождения коэффициента
прохождения T , при
заданных параметрах задачи
Чтобы получить коэффициент T с заданной точность, нужно правильно настроить работу программы.
|
.
,
то коэффициент прохождения будет
равен
с точностью 0.01. Если это условие не
выполняется нужно изменить
пространственные размеры задачи.
нужно выбрать так, чтобы волновой
пакет в начальный момент времени
располагался посередине между левой
границей и барьером. Причем пакет не
должен налезать на барьер и на границу.
На Рис. 05 (a, b)
показано неправильное начальное
расположение волнового пакета. На
Рис. 05 (c) показано
правильное расположение пакета
относительно потенциального барьера
и границ рабочего пространства.
должны
быть достаточно большими. В рабочем
пространстве должны располагаться
отраженная от барьера и прошедшая
через барьер части волнового пакета.
Причем вероятность нахождения
электрона в области барьера при этом
должна быть мала
.
На Рис. 06 (a) показан
неправильный выбор размеров рабочего
пространства для вычисления коэффициента
отражения с заданной точностью. На
Рис. 06 (b) показан
правильный выбор размеров рабочего
пространства для вычисления коэффициента
отражения с заданной точностью.
a b c
Рис. 05.
Расположение пакета в начальный момент времени, a, b – неправильное, c – правильное.
a b
Рис. 06.
Выбор пространственных параметров для вычисления коэффициента прохождения с заданной точностью, a – неправильный выбор, b – правильный выбор.
Задание 1.
Нахождение зависимости коэффициента прохождения T волнового пакета от средней энергии электрона E0.
Для своего варианта с помощью программы
Square_ Barrier
заполнить таблицу для нахождения
зависимости
.
Высоту потенциального барьера положить
равной
.
Ширину барьера a и
полуширину волнового пакета d0
взять в задании для своего варианта.
|
E0 эВ |
15.0 |
17.5 |
20.0 |
22.5 |
25.0 |
27.5 |
30.0 |
32.5 |
35.0 |
37.5 |
|
T |
|
|
|
|
|
|
|
|
|
|
|
E0 эВ |
40.0 |
42.5 |
45.0 |
47.5 |
50.0 |
52.5 |
55.0 |
57.5 |
60.0 |
|
T |
|
|
|
|
|
|
|
|
|
Для
примера с параметрами
эта таблица имеет вид.
|
E0 эВ |
15.0 |
17.5 |
20.0 |
22.5 |
25.0 |
27.5 |
30.0 |
32.5 |
35.0 |
37.5 |
|
T |
0.012 |
0.041 |
0.152 |
0.463 |
0.825 |
0.911 |
0.848 |
0.825 |
0.857 |
0.910 |
|
E0 эВ |
40.0 |
42.5 |
45.0 |
47.5 |
50.0 |
52.5 |
55.0 |
57.5 |
60.0 |
|
T |
0.959 |
0.986 |
0.992 |
0.984 |
0.971 |
0.961 |
0.957 |
0.958 |
0.964 |
График
зависимости
строится в программе Task01.m.
В этой же программе строится зависимость
для монохроматического электрона по
формулам (А.8), (А.9). Для приведенного
примера результат показан на Рис. 07.
Рис. 07.
Зависимость коэффициента прохождения от энергии электрона. Синяя линия – монохроматический электрон. Зеленая линия – волновой пакет.
После
заполнения таблицы с помощью программы
Square_ Barrier
, эти результаты использовать в программе
Task01 для построения
графика зависимости
.
Напомним, что в программе Square_ Barrier и программе Task01 нужно подставить параметры своего варианта.
Задание 2.
Нахождение зависимости коэффициента прохождения T волнового пакета от полуширины пакета d0.
Для своего варианта с помощью программы
Square_ Barrier
заполнить таблицу для нахождения
зависимости
.
Ширину потенциального барьера положить
равной
.
Высоту барьера U0
и среднюю энергию волнового пакета E0
взять в задании для своего варианта.
|
|
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
T |
|
|
|
|
|
|
|
|
|
Для
примера с параметрами
эта таблица имеет вид.
|
|
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
T |
0.235 |
0.177 |
0.136 |
0.105 |
0.080 |
0.060 |
0.043 |
0.031 |
0.022 |
График
зависимости
строится в программе Task02.m.
Для приведенного примера результат
показан на Рис. 08.
Рис. 08.
Зависимость коэффициента прохождения от полуширины волнового пакета.
Результаты,
взятые из таблицы, показаны на Рис. 08
красными кружками. На Рис.08 также показана
синяя линия зависимости
,
построенная по формулам (А.18), (А.19). Как
видно, результаты прекрасно согласуются.
Такое же согласование должно быть в
каждом варианте, если правильно настроить
программу Square_
Barrier, и с помощью
нее заполнить указанную выше таблицу.
Варианты к заданию 1.
-
Номер варианта
1
2.0
5.0
2
2.1
5.0
3
2.2
5.0
4
2.3
5.0
5
2.4
5.0
6
2.5
5.0
7
2.0
5.5
8
2.1
5.5
9
2.2
5.5
10
2.3
5.5
11
2.4
5.5
12
2.5
5.5
13
2.0
6.0
14
2.1
6.0
15
2.2
6.0
16
2.3
6.0
17
2.4
6.0
18
2.5
6.0
19
2.0
6.5
20
2.1
6.5
21
2.2
6.5
22
2.3
6.5
23
2.4
6.5
24
2.5
6.5
25
2.0
7.0
26
2.1
7.0
Варианты к заданию 2.
a =5
-
Номер варианта
E0, эВ
,
эВ1
20
20
2
19.9
20
3
19.8
20
4
19.7
20
5
19.6
20
6
30
30
7
29.9
30
8
29.8
30
9
29.7
30
10
29.6
30
11
40
40
12
39.9
40
13
39.8
40
14
39.7
40
15
39.6
40
16
50
50
17
49.9
50
18
49.8
50
19
49.7
50
20
49.6
50
21
60
60
22
59.9
60
23
59.8
60
24
59.7
60
25
59.6
60
26
59.5
60