Наука и техника: Электрон в слое, Курсовая работа

Министерство Образования, Молодежи и Спорта

Республики Молдова

 

Государственный университет Молдовы


Курсовая Работа

 

Тема:  Электрон в слое.


Работу выполнил

студент 3-го курса:

Радченко Андрей


Кишинёв 1997 г.
Микрочастица (электрон) в слое.

Собственно говоря, одномерная задача, которая сейчас будет рассмотрена, во многих учебных руководствах довольно подробно разобрана путём введения некоторых  упрощений.

Она состоит в следующем :

Микрочастица (электрон) движется вдоль оси x, и её движение полностью определяется следующим гамильтонианом :

 

ì -ћ2/(2m)׶2/¶x2 + U0  ,  x < -a

 Ù   ï

 H =       í -ћ2/(2m0)׶2/¶x2      ,  -a < x < a

ï

î -ћ2/(2m)׶2/¶x2 + U0  ,  x > a

Где  m   -   эффективная масса электрона в областях I , III ;

         m0 -   эффективная масса электрона в области II.

Запишем уравнение Шрёдингера для каждой области :

  ì ¶2YI/¶x2 + 2m/ћ2×(E - U0)YI = 0    ,  x £ -a

  ï

  í ¶2YII/¶x2 + 2m02×E×YI = 0         ,  -a £ x £ a

ï

  î ¶2YIII/¶x2 + 2m/ћ2×(E - U0)×YI = 0       ,  x ³ a

 

 

Область I :   

Общий вид решения уравнения Шрёдингера для 1-ой области записывается сразу :

YI(x) = A×exp(n×x) + B×exp(-n×x).

Используя свойство ограниченности волновой функции, мы придём к тому что B = 0. Значит,

YI(x) = A×exp(n×x).

Волновая функция для второй области тоже элементарно определяется :

YII(x) = C×exp(i×k×x) + D×exp(-i×k×x).

Функция состояния для третьей области выглядит так :

YIII(x) = F×exp(-n×x).

Где

               k = (2m0×E/ћ2)1/2

               n = (2m×(U0-E)/ћ2)1/2.

Стратегия наших дальнейших действий будет состоять в следующем :

¨   Напишем систему из  4  уравнений, удовлетворение которых эквивалентно удовлетворению функциями граничным условиям.

¨   В этой системе из  4  уравнений будут фигурировать неизвестные коэффициенты  A,C,D и F. Мы составим линейную однородную систему относительно них.

¨   Ясно, что существование нетривиальных решений допускается только в случае когда детерминант системы равен нулю. Как выяснится чуть позже, из этого весьма полезного факта мы извлечём уравнение, корнями которого будут возможные уровни энергии.

Приступим к осуществлению первого пункта, т.е. запишем условия сшивания волновых функций :

   YI(x=-a) = YII(x=-a)

YII(x=a) = YIII(x=a)

     YI¢(x=-a)/m = YII¢(x=-a)/m0

   YII¢(x=a)/m0 = YIII¢(x=a)/m

А в наших определениях этих функций это выглядит так :

A×exp(-n×a) = C×exp(-i×k×a) + D×exp(i×k×a)

m-1×A× n×exp(-n×a) = i×k×/m0×(C×exp(-i×k×a) - D×exp(i×k×a))

C×exp(i×k×a) + D×exp(-i×k×a) = F×exp(-n×a)

i×k×/m0×(C×exp(i×k×a) - D×exp(-i×k×a)) = - n/m×F×exp(-n×a).

Теперь составим определитель :

|exp(-n×a)     -exp(-i×k×a)                     -exp(i×k×a)                       0                         |

|m-1×n×exp(-n×a)                                   -1/m0×i×k×exp(-i×k×a)        1/m0×i×k×exp(i×k×a)   0            |

|0                  exp(i×k×a)                         exp(-i×k×a)                       -exp(-n×a)         |

|0                  1/m0×i×k×exp(i×k×a)            -1/m0×i×k×exp(-i×k×a)        1/m×n×exp(-n×a)|

Если теперь раскрыть этот определитель по обычным правилам и приравнять его к нулю, то мы получим следующее уравнение для уровней энергии:

((n/m)2 - (k/m0)2)×Sin(2×k×a) + 2×k×n/(m×m0)×Cos(2×k×a) = 0.

Это уравнение решается численным методом, а именно,  методом Ньютона.

Найдём неизвестные коэффициенты A, C, D, F для более полного описания волновой функции. Для этого воспользуемся некоторыми соотношениями, которые непосредственно вытекают из условий сшивания и условия нормировки.

C = F×exp(-n×a)×{exp(i×k×a) + exp(-3×i×k×a) ×( i×k/m0 - n/m)/(n/m + i×k/m0)}

D = C×exp(-2×i×k×a)×( i×k/m0 - n/m)/(n/m + i×k/m0)

A = exp(n×a)×(C×exp(-i×k×a) + D×exp(i×k×a)) .

Поскольку A, C и D линейно зависят от F, то целесообразно ввести обозначения :

A = RA×F

C = RC×F

D = RD×F.

RA, RC, RD - известные постоянные.

Таким образом, если мы каким-то образом узнаем константу F, то мы определим остальные константы A, C, D.  А сделаем мы это с помощью условия нормировки.

Действительно :

YI(x) = F×RA×exp(n×x)

YII(x) = F×( RC×exp(i×k×x) + RD×exp(-i×k×x)).

YIII(x) = F×exp(-n×x).

I1 + I2 + I3 = 1

Где

I1 = |F|2×|RA|2×òQexp(2×n×x)×dx = |F|2×|RA|2×(2×n)-1×exp(2×n×x) =

= |F|2×|RA|2×(2×n)-1×exp(-2×n×a)

I2 = |F|2×{ òL|RC|2×dx + òL|RD|2×dx + RC×RD*×òLexp(2×i×k×x)×dx +

+ RC*×RD×òLexp(-2×i×k×x)×dx } = |F|2×{ 2×a×(|RC|2 + |RD|2) +

((exp(2×i×k×a) - exp(-2×i×k×a))×RC×RD*/(2×i×k) +

+ i×((exp(-2×i×k×a) - exp(2×i×k×a))×RC*×RD/(2×k) }

I3 = |F|2×òWexp(-2×n×x)×dx = |F|2×(2×n)-1×exp(-2×n×a)

|F|2 = { |RA|2×(2×n)-1×exp(-2×n×a) + 2×a×(|RC|2 + |RD|2) +

((exp(2×i×k×a) - exp(-2×i×k×a))×RC×RD*/(2×i×k) +

+ i×((exp(-2×i×k×a) - exp(2×i×k×a))×RC*×RD/(2×k) + (2×n)-1×exp(-2×n×a) }-1.

Теперь, когда мы знаем  F, нетрудно определить коэффициенты A, C, D, а значит и волновую функцию, характеризующую состояние электрона.

 

Электрон в слоях

 

Задача, которая сейчас будет описана, характеризуется тем, что потенциал обладает пространственной периодичностью. Схематически это изображается так.

То есть, это ни что иное как одномерное движение электрона в периодическом поле. Графически это можно изобразить серией потенциальных барьеров или, как говорят, серией потенциальных ступенек.

Аналитически условие периодичности потенциала записывается весьма просто:

U(x)=U(x+2a)   (1)

Соотношение (1) записано в предположении, что ширина каждой потенциальной ямы равна ширине всякого потенциального барьера.

Ясно, что волновые функции, соответствующие областям I, III, удовлетворяют одному и тому же уравнению Шредингера:

2Y/¶x2 + 2m/ћ2×(E - U0)Y = 0

следовательно эти функции отличаются только постоянным множителем, который называется фазовым множителем.

Этот фазовый множитель мы будем обозначать следующим образом:

r = exp(i 2ak)

Тогда Y(x+2ma) = Y(x)×rm ,    где  m=0, ±1, ±2,...                 (2)

Оказывается, что достаточным для определения дискретного энергетического спектра (рассматривается только случай когда E<U0) и волновой функции является рассмотрение областей I, II, III. Действительно, пользуясь соотношением (2), мы определим волновую функцию на всей действительной оси.

 

Рассмотрим область I:

Уравнение Шредингера для нее записывается в виде:

2YI/¶x2 + 2m22×(E - U0)YI = 0    ,  0 > x > -a

его решение выглядит просто:

YI(x) = A×exp(n×x) + B×exp(-n×x).

Где n = (2m2 (U0-E) /ћ2)1/2

 

Рассмотрим область II:

Уравнение Шредингера для нее записывается в виде:

2YII/¶x2 + 2m12×E YII = 0    ,  a ³ x ³ 0

его решение выглядит просто:

YII(x) = C×exp(i×p×x) + D×exp(-i×p×x).

Где  p = (2m1E/ћ2)1/2

Рассмотрим область III:

2YIII/¶x2 + 2m22×(E - U0)YIII = 0    ,  2a > x > a

его решение выглядит просто:

YIII(x) = r (A×exp(n×x) + B×exp(-n×x)).

Запишем граничные условия:

YI(x=0) = YII(x=0)

YII(x=a) = YIII(x=a)

     YI¢(x=0)/m = YII¢(x=0)/m0

   YII¢(x=a)/m0 = YIII¢(x=a)/m

Подставляя волновые функции в эту систему уравнений, мы получим некоторые связи между коэффициентами A, B, C, D:

A+B=C+D

C exp(i p a)+D exp(-i p a) = exp(i 2 a k) (A exp(n a)+B exp(-n a))

(A-B) n/m2 = (C-D) i p / m1

(C exp(i p a)-D exp(-i p a)) i p / m1 = exp(i 2 a k) n/m2 (A exp(n a)-B exp(-n a))

Следуя приведённым выше соображениям, мы составим определитель :

|1                           1                                    -1                            -1                       |

|exp(i×k×2a+n×a)     exp(i×k×2a-n×a)              -exp(i×p×a)               -exp(-i×p×a)        |

|n/m2                      -n/m2                             -i×p/m1                     i×p/m1                        |

|n/m2exp(i×k×2a+n×a)                                    -n/m2×exp(i×k×2a-n×a)                         - i×p/m1×exp(i×p×a)                                        i×p/m1×exp(-i×p×a)           |

и приравняем его к нулю.

Результатом раскрытия определителя будет весьма громоздкое уравнение содержащее в качестве неизвестного энергию электрона.

Рассчитанные уровни энергии для различных эффективных масс приведены ниже.

a=10;      U=10;              m1=4;               m2=1

 

0.1135703312666857 0.6186359585387896 0.2019199605676639
0.3155348518478819 0.05047267055441365 1.263391478912778
0.4544326758658974 2.137353840637548 0.808172718170137
2.479933076698526 0.4544326758658974 6.168062551132728
5.611693924351967 1.820461802850339 1.529165865668653
1.023077302091622

                                                                                    

a=10       U=10               m1=2                m2=1

 

0.1032788024178655 0.2324238959628721 0.41331603936642
0.6460490460448886 0.930750939555283 1.26759057783714
1.656787195799296 2.098624192369327
2.593469359607937 3.141805331837109
3.744277072860902 5.887485640841992

a=10       U=10               m1=1                m2=1

0.05408120469105441 0.2163802958297131 0.4870681554965061
0.86644533469418 1.354969224117534 1.953300729714778
2.662383817919513 4.418966218448088 7.961581805911094

  

a=10       U=10               m1=0.5 m2=1

0.118992095909544 4.249561710930034 1.068004282376146
0.4754473139332004 5.78216724725356 2.955345679469631
1.895012565781256

a=10       U=10               m1=.25 m2=1

0.2898665804439349 4.30026851446248
2.479039415645616 1.132264393019809

Еще из раздела Наука и техника:


 Это интересно
 Реклама
 Поиск рефератов
 
 Афоризм
Мужчина, добившийся успеха, это тот, кто зарабатывает денег больше, чем может потратить его жена. Женщина, добившаяся успеха, это та, которая сумела найти такого мужчину.
 Гороскоп
Гороскопы
 Счётчики
bigmir)net TOP 100