Задача про электрон в идеальном бесконечном кристалле

в 16:23, , рубрики: квантовая механика, уравнение Шрёдингера, физика полупроводников, физика твердого тела

Здравствуйте, друзья!

Согласно законам квантовой механики, электрон в изолированном атоме может находиться только на определенных дискретных энергетических уровнях. Однако, когда атомы объединяются в периодическую кристаллическую решетку, ситуация меняется: вместо отдельных уровней возникают разрешенные и запрещенные энергетические зоны.

Как образуются зоны? Разобраться в этом поможет квантовая механика!

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

Рассмотрим предельный случай модели Кронига-Пенни. Пусть потенциал, создаваемый одномерной решёткой имеет вид:

 Uleft ( x right )=alphasum_{n=-infty }^{+infty}deltaleft ( x-na right )

График потенциала

График потенциала

Нужно найти энергетический спектр и отвечающий ему волновые функции. Для этого нужно решить следующее стационарное уравнение Шрёдингера:

-frac{hbar^2}{2m}psi^{''}+alpha sum_{n=-infty }^{infty}deltaleft ( x-na right )psi=Epsi

Но в каком виде искать решение?

Трансляционная симметрия. Теорема Блоха

Нетрудно заметить, что если данную систему сместить на величину кратную а, то система перейдёт сама в себя. Это так называемая трансляционная симметрия. Данному преобразованию отвечает оператор трансляции, имеющий вид:

hat{T_{a}}=expleft ( ifrac{hat{p}a}{hbar} right )

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

psileft ( x right )=chi left ( x right )uleft ( x right )

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

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

chi left ( x+a right )=lambda chi left ( x right )

Получили функциональное уравнение Коши, решение которого известно:

chileft ( x right )=exp left ( ikx right )

Данная функция отвечает собственному значению:

lambda=exp left ( ika right )

Возвращаясь к исходной задаче, решение уравнения Шрёдингера будем искать в виде:

psi left ( x right )=exp left ( ikx right )uleft ( x right )

Мы получили функцию Блоха. В 1928 году физик Феликс Блох сформулировал теорему, в которой установил вид волновой функции в периодическом потенциале.

Решение уравнение Шрёдингера

Перейдём к непосредственному отысканию волновых функций, подставим установленный вид решения в уравнение Шрёдингера и получим новое уравнение на функцию u(x), рассмотрим его в области -a<x<0:

begin{matrix} u^{''}+2iku^{'}+left ( beta^2-k^2 right )u=0 \ beta^2=frac{2mE}{hbar^2}end{matrix}

А его решение:

u_1left ( x right )=A_1expleft [ -i(k+beta)x right ]+B_1expleft [ -i(k-beta)x right ]

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

psi_1left ( x right )=A_1expleft [ -ibeta x right ]+B_1expleft [ ibeta x right ]

Очевидно, что в области 0<x<a решение будет иметь такой же вид:

psi_{2}left ( x right )=A_{2}expleft [ -ibeta x right ]+B_{2}expleft [ ibeta x right ]

Но как связать коэффициенты A и B из этих областей? Подействуем оператором трансляции на решение из 1-й области, но тогда оно должно совпадать с решением из 2-й области:

hat{T}_apsi_1left ( x right )=psi_{2}left ( x right )

Из данного условия и получаем связь коэффициентов:

begin{matrix} A_{2}=expleft [ ileft ( k+beta  right )a right ]A_1 \ B_{2}=expleft [ ileft ( k-beta  right )a right ]B_1end{matrix}

Тогда волновая функция во 2-й области принимает вид:

psi_2left ( x right )=expleft{ ikaright}left{ A_1expleft [ -ibetaleft ( x-a right ) right ]+B_1expleft [ ibetaleft ( x-a right ) right ]right}

Теперь нужно "сшить" полученные решения в точке x=0:

left{begin{matrix} psi_1left ( 0 right )=psi_2left ( 0 right )\ left [ psi^{'} right ]_{x=0}=frac{2malpha}{hbar^2}psi_1left ( 0 right )end{matrix}right.

Подробно останавливаться на получении условий сшивки не буду.

Получаем однородную систему линейных уравнений относительно А и B, нетривиальное решение которой существует при равенстве определителя нулю:

begin{matrix} cos left ( beta aright )+frac{gamma }{beta}sin left ( beta a right )=cos left ( ka right ) \ gamma=frac{malpha}{hbar^2}end{matrix}

Получили трансцендентное уравнение на энергетический спектр.

Квазиволновой вектор и энергетические зоны

Находя собственные функции оператора трансляции и решая уравнения Шрёдингера, мы пока не задумывались над физическом смыслом введённой постоянной k. Взглянув на вид Блоховских функций, можно сделать вывод, что постоянная k описывает распространение плоской волны с волновым вектором k. В физике твёрдого тела данный вектор носит название квазиволнового.

Чем же он отличается от обычного волнового вектора? Ответ кроется в полученном выше уравнении на энергетический спектр. В отличие от обычного волнового вектора, квазиволновой определён неоднозначно. Действительно, если к произведению ka прибавить 2π, то полученное уравнение не изменится, и совсем другому значению квазиволнового вектора отвечает тот же спектр энергий, а значит состояния с волновыми векторами k и k+2π/a физически эквивалентны. Для удобства значение k полагают: -π/a<k<π/a. Такое определение исчерпывает физически неэквивалентные состояния.

Перейдём к изучению полученного уравнения на энергетический спектр. Для начала продемонстрирую приведённую зонную схему - зависимость энергии от квазиволнового вектора, которая получается из данного уравнения.

Приведённая зонная схема

Приведённая зонная схема

Видно, что существуют чередующиеся зоны разрешённых и запрещённых энергий.

Рассмотрим следующие интересные случаи полученного уравнения:

begin{matrix} cos left ( beta aright )+frac{gamma }{beta}sin left ( beta a right )=cos left ( ka right ) \ gamma=frac{malpha}{hbar^2}end{matrix}

1) α = 0

В этом случае решение уравнения Шрёдингера в точности совпадает с решением для свободной частицы, запрещённых зон нет, разрешённые зоны заполняют все значения энергии от 0 до ∞.

2) βa >> 1

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

3) βa<<1 и ka<<1

Данный случай описывает частицу у дна разрешённой зоны, а дисперсионное уравнение принимает вид:

1-frac{left ( beta a right )^2}{2}+frac{gamma}{beta}left ( beta a-frac{left ( beta a right )^3}{6} right )=1-frac{left ( ka right )^2}{2}

Откуда находим энергию:

begin{matrix} E=E_0+frac{left ( hbar k right )^2}{2m^*}=E_0+frac{p^2}{2m^*} \ m*=mleft ( 1+frac{agamma}{3} right ) \ E_0=frac{hbar^2 gamma}{maleft ( 1+frac{agamma}{3} right )}end{matrix}

Полученный закон дисперсий схож с законом для свободной частицы, что указывает на то, что электрон с такой энергий движется "свободно" по кристаллу с импульсом p, но обладает массой m* отличной от массы свободного электрона.

В этом случае величину p = ℏk называют квазиимпульсом, а m* - эффективной массой.

Заключение

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

Стоит отметить, что полученные волновые функции ненормируемые, а значит состояния отвечающие им физически нереализуемые. Действительно, ведь в природе не существует бесконечных идеальных кристаллов) Однако данная проблема разрешается введением периодических граничных условий, то есть если кристалл имеет длину L, то волновая функция на левой и правой границах должна совпадать. В это случае квазиволновой вектор принимает уже дискретный набор значений.

Статья подготовлена сообществом.

Автор: Syr04ek

Источник

* - обязательные к заполнению поля


https://ajax.googleapis.com/ajax/libs/jquery/3.4.1/jquery.min.js