Частка в одновимірному періодичному потенціалі

Матеріал з Вікіпедії — вільної енциклопедії.
Перейти до навігації Перейти до пошуку

У квантовій механіці, частинка в одновимірному періодичному потенціалі це ідеалізована задача, яка може бути розв'язана точно (при деяких спеціального вигляду потенціалах), без спрощень. Припускається, що потенціал нескінченний та періодичний, тобто має трансляційну симетрію, що взагалі не виконується в реальних кристалах, де завжди існує як мінімум один дефект - поверхня (це приводить до іншої задачі про поверхневі стани або таммовські рівні).

Постановка задачі

[ред. | ред. код]

Розглянемо одномірну решітку позитивних іонів, віддаль між якими . Потенціал в цьому випадку буде виглядати подібно до цього:

Рівняння Шредінгера в цьому випадку буде:

з періодичним потенціалом Загальний вигляд розв'язку рівняння Шредінгера з періодичним потенціалом, згідно з теоремою Блоха

де u(x) — деяка періодична функція:

u(x + a) = u(x).

k називають квазіімпульсом, за аналогією з хвильовою функцією eikx для частки з певним імпульсом k.

Як видно, вся хвильова функція визначається певною величиною k та будь-якою ділянкою функції довжиною a.

При наближенні до країв решітки з'являється проблема граничних умов. Зручно представити решітку кільцем деякої довжини L >> a. Якщо L — довжина решітки така що L >> a, тоді число іонів в решітці настільки велике, що хвильові функції електрона в сусідніх іонів не розрізняються.

Тут замість двох граничних умов залишається тільки одне:

Якщо N — число іонів в решітці, то aN = L. Підставляючи хвильову функцію в граничну умову, получаємо квантування для k:

Модель Кроніга-Пенні

[ред. | ред. код]

Для спрощення задачі використовують прямокутний потенціал:

Використовуючи теорему Блоха можна знайти хвильову функцію у всьому просторі, проте спершу необхідно знайти розв'язок для одного періоду, та зробити його гладким на краях, тобто "зшити" значення сусідніх функцій та їхніх похідних.

Розглянемо один період потенціалу:

У нас є дві незалежні області для яких можна знайти розв'язки:

Для знаходження u(x) в кожній області необхідно зробити такі перетворення:

Аналогічним чином отримуємо

Щоб знайти повний розв'язок необхідно впевнитися в гладкості функції на границях:

та періодичності u(x) u'(x)

Ці умови дають наступну матрицю:

Для існування нетривіального розв'язку необхідно обнулення детермінанту цієї матриці. Після деяких перетворень отримуємо:

Для подальшого спрощення можна виконати такі спрощення, зміст яких полягає в переході до дельта-подібних потенціалів (гребінець Дірака) :

Тоді кінцевий результат буде:

Програмний код

[ред. | ред. код]

Код для Maple

[ред. | ред. код]

Наступний програмний код написаний на мові Maple (9.5). Він є безпосередньо графічний розв'язок .

  restart;
  with(plots):
  with(stats[statplots]):
  eq:=cos(k*a)=cos(beta*b)*cos(alpha*(a-b)) - (alpha^2+beta^2)/(2*alpha*beta)*sin(beta*b)*sin(alpha*(a-b));
  alpha:=sqrt(8*Pi^2*m*(E)*e/h^2):
  beta:=sqrt(8*Pi^2*m*(E+V)*e/h^2):
  e:=1.6*1e-19:
  a:=0.54310*1e-9:
  m:=0.19*9.1*1e-31:
  b:=1/5*a:
  h:=6.6*1e-34:
  k(E,V):=arccos(rhs(evalf(eq)));

  #Графік
  p:=plot({subs(V=10,k(E,V)),subs(V=10,-k(E,V))},E=-5..50,labels=[ka, E],color=blue):
  xyexchange(p);

  #Анімація, залежність від глибини ями
  p:=animate( plot, [{k(E,V),-k(E,V)},E=-10..50, color=blue,labels=[ka, E]], V=0..30 ):
  xyexchange(p);

На малюнках подані графічні розв'язки рівняння ( * ).

Лінії відповідають дозволеним значенням енергії. Існують області енергії, при любих значеннях хвилевого вектора неможливе існування електрону.
Лінії відповідають дозволеним значенням енергії. Показано рух закону дисперсії в залежності від глибини потенціальної ями.

На правому малюнку видно, як при деякому значенні потенційної енергії можливе утворення одновимірногобезщільового напівпровідника.

Код для Scilab

[ред. | ред. код]
Ліеії иак само відповідають дозволеним значенням енергії. Синім позначено розв"язок для моделі Кроніга-Пенні, червоним - гребінця Дірака при тих же значеннях V0b

Код нижче є фактично перекладом попередньої програми на мову Scilab, за тим виключенням, що ілюструє також і випадок переходу до гребінця Дірака.

clear all

global Pi e a m b h

Pi = 3.1415926;

step = 0.1;

e = 1.6 * 1e-19;
a = 0.54310 * 1e-9;
m = 0.19*9.1 * 1e-31;
b = 1/5 * a;
h = 6.6 * 1e-34;

function [alpha, beta] = ab(V,E)
  alpha = sqrt(8*Pi^2*m*(E)*e/h^2);
  beta  = sqrt(4*Pi^2*m*(E+V)*e/h^2);
endfunction

function r=kronigpenney(V, E)
  [alpha, beta] = ab(V,E);
  r = 1/a * acos((cos(beta*b) .* cos(alpha*(a-b)) ) - (alpha.^2+beta.^2) / (2*alpha .* beta) .* sin(beta*b) .* sin(alpha*(a-b)));
endfunction

function r=dirac(V,E)
  [alpha, beta] = ab(V,E);
  r = 1/a * acos(cos(alpha * a) - (beta.^2 * b * a) / 2 .* sin(alpha*a) / (alpha * a));
endfunction

E = [0 : step: 50];

k = kronigpenney(10, E);
plot(k, E, 'b'); plot(-k, E, 'b');

k = dirac(10, E);
plot(k, E, 'r'); plot(-k, E, 'r');

Посилання

[ред. | ред. код]