Туннелирование электрона через потенциальный барьер прямоугольной формы
Московский институт электронной техники Кафедра КФН
Лабораторная работа 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)
3)
4)
На Рис. 01 показана начальная стадия движения волнового пакета. Рис.01. Начальная стадия движения волнового пакета.
Программа Square_ Barrier моделирует движение электрона в состоянии с волновой функцией Синяя кривая на Рис. 01 показывает положение пакета в начальный момент времени. В качестве начального значения пакета берется гауссов волновой пакет (смотри лабораторную работу 2, приложение А, формула (А.9)). На Рис. 01 видно, что среднее положение пакета выбрано равным Красная кривая показывает положение пакета в момент времени На Рис. 01 прямоугольный потенциальный барьер изображен синими толстыми линиями. Показана также высота потенциального барьера в электрон-вольтах Светло-сиреневая линия показывает среднюю энергию электрона
Здесь
Из формул (4), (5) видно, что чем меньше полуширина волнового пакета
На Рис. 02 показан момент взаимодействия волнового пакета с потенциальным барьером. Рис.02. Момент взаимодействия пакета с потенциальным барьером.
Красная кривая показывает положение пакета в момент времени На Рис. 02 показана также часть пакета, которая прошла через барьер, туннелировала через барьер, и двигается направо, удаляясь от барьера. Величина этой части пакета определяет вероятность туннелирования электрона через барьер. В данном примере эта вероятность равна Наконец Рис. 02 показана часть пакета, которая локализована внутри барьера. Величина этой части пакета определяет вероятность нахождения электрона внутри барьера. В данном примере эта вероятность равна
3. Коэффициент отражения от барьера R (Reflectivity), и коэффициент пропускания через барьер T (Transmissivity), в задаче о туннелировании волнового пакета, определяются следующими формулами.
Заметим, что формулы (6) применяются для достаточно больших моментов времени. Нужно чтобы произошло взаимодействие волнового пакета с потенциальным барьером. Части волнового пакета разбежались от барьера, а вероятность нахождения электрона внутри барьера стала равной нулю Реально будем выбирать такие времена, когда вероятность На Рис. 03 показан момент, при котором можно определить коэффициенты Рис. 03. Момент для определения коэффициентов
Отметим, что на Рис. 03 увеличены размеры пространства, в котором движется пакет по сравнению с Рис. 02. Это сделано, для того чтобы в рабочем пространстве уместились отраженная и прошедшая части волнового пакета. В этом случае можно считать, что волновой пакет равен нулю за пределами рабочего пространства. Тогда вероятности
В данном примере получаем следующие коэффициенты отражения и пропускания.
Сумма этих коэффициентов отличается от единицы На следующем Рис. 04 выбрано еще большее время движения пакета и большее рабочее пространство. Рис. 04. Момент для определения коэффициентов
В данном примере получаем следующие коэффициенты отражения и пропускания. Сумма этих коэффициентов отличается от единицы
4. Заметим, что можно ускорить движение пакета, увеличивая временной шаг dt. Однако слишком большой шаг делает движение пакета не плавным, а рывками. В программе Square_ Barrier предусмотрен флаг flagStop, который позволяет менять режим работы программы. Если флаг равен нулю flagStop = 0, то программа может работать в течение любого времени и части волнового пакета могут уйти за пределы рабочего пространства. Если же флаг равен единице flagStop = 1, то прекращается движение пакета в тот момент, когда части волнового пакета, достигают границ рабочего пространства. Рекомендуется выбрать значение флага flagStop = 1.
a b c Рис. 05. Расположение пакета в начальный момент времени, a, b – неправильное, c – правильное.
a b Рис. 06. Выбор пространственных параметров для вычисления коэффициента прохождения с заданной точностью, a – неправильный выбор, b – правильный выбор. Задание 1. Нахождение зависимости коэффициента пропускания T волнового пакета от средней энергии электрона E 0. Для своего варианта с помощью программы Square_ Barrier заполнить таблицу для нахождения зависимости
Для примера с параметрами
График зависимости Рис. 07. Зависимость коэффициента пропускания от энергии электрона. Синяя линия – монохроматический электрон. Зеленая линия – волновой пакет.
После заполнения таблицы с помощью программы Square_ Barrier, эти результаты использовать в программе Task01 для построения графика зависимости Напомним, что в программе Square_ Barrier и программе Task01 нужно подставить параметры своего варианта.
Задание 2. Нахождение зависимости коэффициента пропускания T волнового пакета от полуширины пакета d 0.
Для своего варианта с помощью программы Square_ Barrier заполнить таблицу для нахождения зависимости
Для примера с параметрами
График зависимости Рис. 08. Зависимость коэффициента пропускания от полуширины волнового пакета.
Варианты к заданию 1.
Варианты к заданию 2. a =5
Приложение A. Туннелирование электрона через потенциальный барьер прямоугольной формы
Решаем уравнение Шредингера
с потенциальной энергией в виде прямоугольного барьера Здесь
Здесь Плотность вероятности начального распределения дается формулой
V (x) k 0 V 0 E 0
x 0 0 a x
После взаимодействия с барьером пакет
Другая часть проходит через барьер и уходит направо, и позволяется найти коэффициент пропускания по формуле
Еще она часть находится в области барьера
Так как волновая функция то коэффициенты R, T, W связаны соотношением
Во время работы пакета nsse в каждый момент времени вычисляются значения коэффициентов отражения R и пропускания T. Чтобы вывести их на экран, надо нажать клавишу F1. Можно убедиться, что вычисленные программой значения R, T удовлетворяют неравенству Это находится в соответствие с формулами (А.3) – (А.6)
Если рассматривать туннелирование монохроматических электронов, то для коэффициентов отражения и пропускания существуют простые аналитические формулы. Монохроматический электрон, это электрон у которого энергия, импульс и волновое число имеют определенные значения Волновая функция монохроматического электрона имеет вид
Плюс и минус перед k означает, что электрон движется направо или налево. Для электрона с заданной энергией E коэффициент пропускания T задается следующими формулами. Если
Если
Если ширину барьера a измерять в Å (ангстремах), энергии E, V0 в eV (электрон-вольтах), то вычисления в формулах (А.8), (А.9) надо выполнять по формулам
Электрон, движение которого описывается волновым пакетом, не имеет определенной энергии E. Это связано с тем, что начальный волновой пакет
с разными волновыми числами k и, соответственно, разными энергиями E = E(k) (А.7). Это можно увидеть, разложив
Здесь коэффициенты C(k) имеют физический смысл. Так |C(k)|2 – плотность вероятности того, что электрон в состоянии
Коэффициенты C(k) можно найти, выполнив обратное преобразование Фурье Тогда плотность вероятности для волнового числа будет определяться формулой
На рисунке приведены распределения |F0(x)|2 и |C(k)|2
|F0(x0)|2 /e |C(k0)|2 /e
x0–ap x0 x0+ap x k0–1/ap k0 k0+1/ap k
Из этих рисунков видно, что электрон в момент времени t = 0 локализован в интервале x0–ap < x < x0+ap, где ap - полуширина пакета. Тогда величину Dx = 2ap можно назвать неопределенность координаты электрона. Аналогично волновые числа электрона лежат в интервале k0–1/ap < k < k0+1/ap и величину Dk = 2/ap можно назвать неопределенность волнового числа электрона. Если вспомнить, что импульс электрона p и волновое число k связаны соотношением то неопределенность электрона можно найти по формуле Отсюда получаем соотношение неопределенности Гейзенберга Кроме того, благодаря связи (А.13) энергия электрона E также имеет разброс DE вокруг среднего значения E0
Если ширину пакета ap измерять в Å, энергию E в eV, то имеем вычислительную формулу Таким образом, локализованный электрон, волновая функция которого имеет вид волнового пакета, не имеет определенной энергии. Можно только говорить о средней энергии E0, и интервале DE, в котором лежит энергия электрона. Поэтому туннелирование такого электрона через потенциальный барьер, будет отличаться от результатов, поученным по формулам (А.8), (А.9) для монохроматического электрона. Результат будет зависеть от того, какие энергии принимают участие в формировании волнового пакета (А.12) Из формулы (А.15) видно, что чем шире волновой пакет, тем меньше разброс в энергиях, и тем ближе рассматриваемый электрон к монохроматическому электрону
Воспользуйтесь поиском по сайту: ![]() ©2015 - 2025 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...
|